/CLEAR

*ASK,VOF,Enter the volume fraction of fibers (0 < VOF < 0.785),0.4
FINISH
/FILNAM,Long FIberA
/Title,Long FIberA
/UIS,MSGPOP,3

PI=3.141592654
LENGTH=1
HEIGHT=0.1
EXTSIZE=1
AREA=LENGTH*LENGTH
VOF=0.40

R=SQRT(VOF*LENGTH*LENGTH/PI)

EX11=379.3e9
PRXY112=0.1

EX2=68.3e9
PRXY2=0.30


/PREP7

!---------------------Flexible function---------------------


BLC4,0,0,LENGTH/2,LENGTH/2
CYL4,0,0,R,90

AOVLAP,ALL
AGLUE,ALL
NUMCMP,AREA


! Define element type
!======================
ET,1,PLANE82   
KEYOPT,1,1,0
KEYOPT,1,6,0
KEYOPT,1,10,0  

ET,2,SOLID45

! Define material properties Orthotropic Material
!============================

! Material #1
!---------------

MP,EX,1,EX11
MP,PRXY,1,PRXY112



! Material #2
!---------------

MP,EX,2,EX2
MP,PRXY,2,PRXY2

NN=15
LESIZE,5,,,2*NN
LSEL,S,,,10,11
LESIZE,ALL,,,NN
LSEL,S,,,8,9
LESIZE,ALL,,,2*NN
LSEL,S,,,2,3
LCCAT,ALL
LESIZE,ALL,,,NN
ALLSEL,ALL

MAT,1
AMESH,1

MAT,2
AMESH,2

ARSYM,X,ALL
ARSYM,Y,ALL
NUMMRG,NODE
NUMMRG,KP
NUMCMP,NODE