【讨论】急求fracture 分析的例子！

BosFEA

Verification手册中有例题，输入文件也有。搜索一下即可。

BosFEA

手册在ADINA安装目录下的DOC子目录下。
这里面对各个例题都是详细说明，如B123就是裂纹计算的例子。
FEPROGRAM ADINA
****ADD:BATCH
*PLSYSTEM NULL
*FILEECHO OPTION=FILE FILE='b123.ilg'
*FILELOG OPTION=FILE FILE='b123.ilg'
****END:BATCH
****ADD:INTERACTIVE
FILEECHO OPTION=NONE
****END:INTERACTIVE
CONTROL PROMPT=NO
UPDATE FRAME UPPER=HEADING
*
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
*
egroup twodsolid 1 strain form=mixed
gsurface 4 node=9
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'
****ADD:SNAPSHOT
*SNAPSHOT 'b123.ps'
****END:SNAPSHOT
****ADD:BATCH
*QUIT IMMEDIATE=YES
****END:BATCH
  
别忘了读入文件时点击Batch按钮，直到完全读完才能求解。