*do,i,1,41
Vf%i%=0.1+2.4*(XX%i%-XX%i%**2)
El%i%=388e9*Vf%i%+3.5e9*(1-Vf%i%)
Et%i%=3.5e9*(1+2*r1*Vf%i%)/(1-r1*Vf%i%)
Glt%i%=1.3e9*(1+r2*Vf%i%)/(1-r2*Vf%i%)
Gtt%i%=1.3e9*(1+r3*Vf%i%)/(1-r3*Vf%i%)
vlt%i%=0.23*Vf%i%+0.35*(1-Vf%i%)
vtt%i%=Et%i%/(2*Gtt%i%)-1
vtl%i%=vlt%i%*Et%i%/El%i%
C11%i%=El%i%*(1-vtt%i%**2)/(1-2*vlt%i%*vtl%i%-vtt%i%**2-2*vlt%i%*vtl%i%*vtt%i%)
C22%i%=Et%i%*(1-vlt%i%*vtl%i%)/(1-2*vlt%i%*vtl%i%-vtt%i%**2-2*vlt%i%*vtl%i%*vtt%i%)
C12%i%=El%i%*(vtl%i%+vtt%i%*vtl%i%)/(1-2*vlt%i%*vtl%i%-vtt%i%**2-2*vlt%i%*vtl%i%*vtt%i%)
C23%i%=Et%i%*(vtt%i%+vlt%i%*vtl%i%)/(1-2*vlt%i%*vtl%i%-vtt%i%**2-2*vlt%i%*vtl%i%*vtt%i%)
*enddo
et,1,186
blc4,,,1,1,h
TB,ELAS,1,,,AELS
*DO,I,1,41
!0°
!TBFIELD,XCOR,XX%I%
!90°
TBFIELD,YCOR,XX%I%
tbdata,1,C11%I%,C12%I%,C12%I%
tbdata,7,C22%I%,C23%I%
tbdata,12,C22%I%
tbdata,16,Glt%I%
tbdata,19,Gtt%I%
tbdata,21,Glt%I%
*enddo
esize,0.0125
!lesize,12,,,8
!lesize,12,,,12
!lesize,12,,,16
!lesize,12,,,20
lesize,12,,,24
vmesh,all
nsel,s,loc,y,0
nsel,a,loc,y,1
d,all,uz,0,,,,ux
nsel,s,loc,x,0
nsel,a,loc,x,1
d,all,uz,0,,,,uy
ALLSEL