*--------------------------------------------------------------------------------------
*
FEPROGRAM ADINA
*
*--------------------------------------------------------------------------------------
*
head 'B.123 FRACTURE ANALYSIS OF AN ELASTIC-PLASTIC CT-SPECIMEN'
*
*--------------------------------------------------------------------------------------
*
master idof=100111 reaction=yes frac=yes
iteration full yes
kinematics disp=large strain=large
printout vol=maximum
porthole form=yes
*
timestep default
@clear
5 1.0
TIMEFUNCTION 1
@clear
0. 0./ 5. 5.
dataend
*
subdivide default mode=division prog=arith
*
*--------------------------------------------------------------------------------------
*
coordinates / ent n y z
1 .0 .0
2 .022637 .0
3 .05035 .0
4 .0 .01725
5 .05035 .01725
6 .0 .03
7 .05035 .03
*
line st 1 1 2
surface vertex 1 2 4 1 2
surface vertex 2 2 5 4 2
surface vertex 3 2 3 5 2
surface vertex 4 4 5 7 6
subdivide surface 1 ndiv1=6 ndiv2=2 ratio1=.5
subdivide surface 2 ndiv1=6 ndiv2=4 ratio1=.5
subdivide surface 3 ndiv1=6 ndiv2=2 ratio1=.5
subdivide surface 4 ndiv1=4 ndiv2=1
*
*--------------------------------------------------------------------------------------
*
material plastic-multilinear 1 e=2.1E11 nu=.3
.001 21.28E7 / .005 25.93E7
.0125 31.25E7 / .076 43.88E7
.280 49.20E7
*
*--------------------------------------------------------------------------------------
*
* 單元類型選擇form=mixed
*
*--------------------------------------------------------------------------------------
*
egroup twodsolid 1 strain form=mixed
*
gsurface 4 node=9
*
*--------------------------------------------------------------------------------------
*
* 注意crack類型的設置
*
*--------------------------------------------------------------------------------------
*
gsurface 1 node=9 crack-type=point tip-point=2 tip-option=circular-arc,
radius=0.00001 q-point=mid cpoint1=8
1
2
3
dataend
*
*--------------------------------------------------------------------------------------
*
fixity ztran
'z-translation'
*
fixboundary line ztran
1
fixboundary point all
1
*
*--------------------------------------------------------------------------------------
*
load displacement 1 dz=0.0001
*
apply-load
1 displacement 1 point 5
*
*--------------------------------------------------------------------------------------
*
frame
boundepiction bcode=all
meshplot
loadplot
****ADD:SNAPSHOT
*SNAPSHOT 'b123.ps' APPEND=NO
****END:SNAPSHOT
****ADD:INTERACTIVE
READ END
****END:INTERACTIVE
*
*--------------------------------------------------------------------------------------
*--------------------------------------------------------------------------------------
*
fracture method=both
*
crack-prop point
2
dataend
*
j-virtual-shift ring 1 ring-num=1 vector=input vy=-0.00001
2
dataend
j-virtual-shift ring 2 ring-num=2 vector=input vy=-0.00001
2
dataend
j-virtual-shift ring 3 ring-num=3 vector=input vy=-0.00001
2
dataend
*
j-line ring 1 nrings=1 point=2
j-line ring 2 nrings=2 point=2
j-line ring 3 nrings=3 point=2
*
*--------------------------------------------------------------------------------------
*--------------------------------------------------------------------------------------
*
meshannotation
'axes' off
dataend
copy boundepiction factory default
frame
meshwindow tip size=direct xvmin=0.02260 xvmax=0.02268 yvmin=0 yvmax=0.00006
meshplot meshwindow=tip nodep=symbols_numbers subframe=topleft
modeldep geom=no or=yes def=yes dispo=len dispf=5
vsdepiction plotshift=yes shift=1
meshplot subframe=topright
vsdepiction plotshift=yes shift=2
meshplot subframe=bottomleft
vsdepiction plotshift=yes shift=3
meshplot subframe=bottomright
****ADD:SNAPSHOT
*SNAPSHOT 'b123.ps'
****END:SNAPSHOT
****ADD:INTERACTIVE
READ END
****END:INTERACTIVE
*
*--------------------------------------------------------------------------------------
*
frame
modeldep def=no
copy meshannotation factory default
meshplot
lcplot contour=line_contour-1
lcplot contour=line_contour-2
lcplot contour=line_contour-3
*
*--------------------------------------------------------------------------------------
*
adina file='b123.dat'
*
*--------------------------------------------------------------------------------------