命令流文件(3维桩土接触试验)

tumu123

*
BODY BLOCK NAME=4 OPTION=CENTERED POSITION=VECTOR ORIENTAT=SYSTEM,
     CX1=0.80000000000000 CX2=0.80000000000000,
     CX3=-0.500000000000000 SYSTEM=0 DX1=1.60000000000,
     DX2=1.600000000000 DX3=1.00000000000000
*
BODY BLOCK NAME=1 OPTION=CENTERED POSITION=VECTOR ORIENTAT=SYSTEM,
     CX1=0.80000000000000 CX2=0.80000000000000,
     CX3=-0.500000000000000 SYSTEM=0 DX1=0.240000000000,
     DX2=0.240000000000 DX3=1.00000000000000
*
BODY CYLINDER NAME=2 OPTION=CENTERED POSITION=VECTOR ORIENTAT=SYSTEM,
     CX1=0.80000000000000 CX2=0.80000000000000,
     CX3=-0.350000000000000 SYSTEM=0 AXIS=ZL,
     RADIUS=0.0200000000000000 LENGTH=0.700000000000000
*
BODY BLOCK NAME=3 OPTION=CENTERED POSITION=VECTOR ORIENTAT=SYSTEM,
     CX1=0.80000000000000 CX2=0.80000000000000,
     CX3=0.00500000000000000 SYSTEM=0 DX1=0.24000000000000,
     DX2=0.24000000000000 DX3=0.0100000000000000
*
BODY SUBTRACT NAME=4 KEEP-TOO=YES KEEP-IMP=YES
@CLEAR
1
@
*
BODY SUBTRACT NAME=1 KEEP-TOO=YES KEEP-IMP=YES
@CLEAR
2
@
*
BODY MERGE NAME=2 KEEP-TOO=NO
@CLEAR
3
@
*
FIXITY NAME=FX
@CLEAR
'X-TRANSLATION'
'OVALIZATION'
@
*
FIXITY NAME=FY
@CLEAR
'Y-TRANSLATION'
'OVALIZATION'
@
*
FIXITY NAME=FZ
@CLEAR
'Z-TRANSLATION'
'OVALIZATION'
@
*
FIXBOUNDARY FACES FIXITY=ALL BODY=1
@CLEAR
5  'FZ'
@
*
FIXBOUNDARY FACES FIXITY=ALL BODY=4
@CLEAR
6  'FX'
3  'FX'
4  'FY'
2  'FY'
5  'FZ'
@
*
LOAD MASS-PROPORTIONAL NAME=1 MAGNITUD=9.80000000000000,
     AX=0.00000000000000 AY=0.00000000000000 AZ=-1.00000000000000,
     INTERPRE=BODY-FORCE
*
APPLY-LOAD BODY=2
@CLEAR
1  'MASS-PROPORTIONAL' 1  'MODEL' 0 0 1 0.00000000000000 0 -1 0 0 0  ,
'NO' 0.00000000000000 0.00000000000000 1 0
@
*
TIMEFUNCTION NAME=2 IFLIB=1 FPAR1=0.00000000000000,
     FPAR2=0.00000000000000 FPAR3=0.00000000000000,
     FPAR4=0.00000000000000 FPAR5=0.00000000000000,
     FPAR6=0.00000000000000
@CLEAR
0.00000000000000 0.00000000000000
100.000000000000 1.00000000000000
1.00000000000000E+33 1.00000000000000
@
*
LOAD PRESSURE NAME=1 MAGNITUD=30000.0000000000 BETA=0.00000000000000,
     LINE=0
*
APPLY-LOAD BODY=2
@CLEAR
1  'MASS-PROPORTIONAL' 1  'MODEL' 0 0 1 0.00000000000000 0 -1 0 0 0  ,
'NO' 0.00000000000000 0.00000000000000 1 0
@
*
APPLY-LOAD BODY=2
@CLEAR
1  'MASS-PROPORTIONAL' 1  'MODEL' 0 0 1 0.00000000000000 0 -1 0 0 0  ,
'NO' 0.00000000000000 0.00000000000000 1 0
2  'PRESSURE' 1  'FACE' 4 0 2 0.00000000000000 0 -1 0 2 0  'NO',
     0.00000000000000 0.00000000000000 1 0
@
*
APPLY-LOAD BODY=2
@CLEAR
1  'MASS-PROPORTIONAL' 1  'MODEL' 0 0 1 0.00000000000000 0 -1 0 0 0  ,
'NO' 0.00000000000000 0.00000000000000 1 0
2  'PRESSURE' 1  'FACE' 4 0 2 0.00000000000000 0 -1 0 2 0  'NO',
     0.00000000000000 0.00000000000000 1 0
@
*
CGROUP CONTACT3 NAME=1 FORCES=YES TRACTION=YES NODETONO=NO,
     FRICTION=0.685300000000000 EPSN=0.00000000000000,
     EPST=0.00000000000000 DIRECTIO=NORMAL CONTINUO=YES,
     INITIAL-=ALLOWED PENETRAT=ONE DEPTH=0.00000000000000,
     OFFSET=0.00000000000000 OFFSET-T=CONSTANT CORNER-C=NO,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TIED=NO,
     TIED-OFF=0.00000000000000 HHATTMC=0.00000000000000,
     FCTMC=0.500000000000000 FTTMC=0.500000000000000 RIGID-TA=NO,
     NORMAL-S=1.00000000000000E+11 TANGENTI=0.00000000000000,
     PTOLERAN=1.00000000000000E-08 RESIDUAL=0.00100000000000000,
     LIMIT-FO=1.00000000000000 ITERATIO=2 TIME-PEN=0.00000000000000,
     CONSISTE=DEFAULT USER-FRI=NO DESCRIPT='NONE'
*
CONTACTSURFA NAME=1 PRINT=DEFAULT SAVE=DEFAULT SOLID=YES BODY=2,
     ORIENTAT=AUTOMATIC MARQUEEB=0 DESCRIPT='NONE'
@CLEAR
1 1
2 1
3 1
8 1
@
*
CONTACTSURFA NAME=2 PRINT=DEFAULT SAVE=DEFAULT SOLID=YES BODY=1,
     ORIENTAT=AUTOMATIC MARQUEEB=0 DESCRIPT='NONE'
@CLEAR
7 1
8 1
9 1
1 1
@
*
CONTACTPAIR NAME=1 TARGET=1 CONTACTO=2 FRICTION=0.685300000000000,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000,
     HHATTMC=0.00000000000000 FCTMC=0.500000000000000,
     FTTMC=0.500000000000000 NX=1 NY=1 NZ=1
*
COPY ZONE NAME1=WHOLE_MODEL NAME2=MESHPLOT00001 POSITION=1
*
CGZONE NAME=MESHPLOT00001 FEPROGRA=ADINA SUBSTRUC=0 REUSE=1,
     POSITION=1
@CLEAR
1  'ADD'
@
*
MATERIAL MOHR-COULOMB NAME=1 E=3.70000000000000E+07,
     NU=0.260000000000000 PHI=40.0000000000000 PSI=0.00000000000000,
     COH=10000.0000000000 TCUT=1.00000000000000E+10,
     DENSITY=1490.00000000000 DILATION=NO MDESCRIP='NONE'
*
MATERIAL ELASTIC NAME=2 E=2.80000000000000E+10 NU=0.200000000000000,
     DENSITY=2500.00000000000 ALPHA=0.00000000000000 MDESCRIP='NONE'
*
EGROUP THREEDSOLID NAME=1 DISPLACE=DEFAULT STRAINS=DEFAULT MATERIAL=1,
     RSINT=DEFAULT TINT=DEFAULT RESULTS=STRESSES DEGEN=NO FORMULAT=0,
     STRESSRE=GLOBAL INITIALS=NONE FRACTUR=NO CMASS=DEFAULT,
     STRAIN-F=0 UL-FORMU=DEFAULT LVUS1=0 LVUS2=0 SED=NO RUPTURE=ADINA,
     INCOMPAT=DEFAULT TIME-OFF=0.00000000000000 POROUS=NO,
     WTMC=1.00000000000000 OPTION=NONE
*
EGROUP THREEDSOLID NAME=2 DISPLACE=DEFAULT STRAINS=DEFAULT MATERIAL=2,
     RSINT=DEFAULT TINT=DEFAULT RESULTS=STRESSES DEGEN=NO FORMULAT=0,
     STRESSRE=GLOBAL INITIALS=NONE FRACTUR=NO CMASS=DEFAULT,
     STRAIN-F=0 UL-FORMU=DEFAULT LVUS1=0 LVUS2=0 SED=NO RUPTURE=ADINA,
     INCOMPAT=DEFAULT TIME-OFF=0.00000000000000 POROUS=NO,
     WTMC=1.00000000000000 OPTION=NONE
*
SUBDIVIDE BODY NAME=1 MODE=LENGTH SIZE=0.100000000000000
@CLEAR
4
@
*
SUBDIVIDE BODY NAME=2 MODE=LENGTH SIZE=0.0100000000000000
@CLEAR
@
*
SUBDIVIDE FACE NAME=7 BODY=1 MODE=LENGTH SIZE=0.0100000000000000
@CLEAR
8
9
@
*
SUBDIVIDE FACE NAME=1 BODY=2 MODE=LENGTH SIZE=0.0100000000000000
@CLEAR
2
3
@
*
GBODY NODES=4 NCOINCID=BOUNDARIES NCTOLERA=1.00000000000000E-05,
     SUBSTRUC=0 GROUP=1 PREFSHAP=AUTOMATIC COLLAPSE=NO SIZE-FUN=0,
     DELETE-S=NO ANGLE-MI=5.00000000000000 MIDNODES=CURVED,
     METHOD=DELAUNAY PATTERN=0 MESHING=FREE-FORM DEGENERA=YES,
     BOUNDARY=ADVFRONT DEG-EDGE=0 GEO-ERRO=0.00000000000000,
     SAMPLING=20 MIN-SIZE=0.00000000000000 NLAYER=1 NLTABL=0,
     AUTO-GRA=NO SIMULATE=NO PYRAMIDS=NO
@CLEAR
1 0
4 0
@
*
GBODY NODES=4 NCOINCID=BOUNDARIES NCTOLERA=1.00000000000000E-05,
     SUBSTRUC=0 GROUP=2 PREFSHAP=AUTOMATIC COLLAPSE=NO SIZE-FUN=0,
     DELETE-S=NO ANGLE-MI=5.00000000000000 MIDNODES=CURVED,
     METHOD=DELAUNAY PATTERN=0 MESHING=FREE-FORM DEGENERA=YES,
     BOUNDARY=ADVFRONT DEG-EDGE=0 GEO-ERRO=0.00000000000000,
     SAMPLING=20 MIN-SIZE=0.00000000000000 NLAYER=1 NLTABL=0,
     AUTO-GRA=NO SIMULATE=NO PYRAMIDS=NO
@CLEAR
2 0
@
*
GBODY NODES=4 NCOINCID=BOUNDARIES NCTOLERA=1.00000000000000E-05,
     SUBSTRUC=0 GROUP=2 PREFSHAP=AUTOMATIC COLLAPSE=NO SIZE-FUN=0,
     DELETE-S=NO ANGLE-MI=5.00000000000000 MIDNODES=CURVED,
     METHOD=DELAUNAY PATTERN=0 MESHING=FREE-FORM DEGENERA=YES,
     BOUNDARY=ADVFRONT DEG-EDGE=0 GEO-ERRO=0.00000000000000,
     SAMPLING=20 MIN-SIZE=0.00000000000000 NLAYER=1 NLTABL=0,
     AUTO-GRA=NO SIMULATE=NO PYRAMIDS=NO
@CLEAR
@
*
ZONE NAME=PILE NODEATTA=YES GEOMATTA=YES
@CLEAR
'ELEMENT GROUP 2'
@
*
MASTER ANALYSIS=STATIC MODEX=EXECUTE TSTART=0.00000000000000 IDOF=111,
     OVALIZAT=NONE FLUIDPOT=AUTOMATIC CYCLICPA=1 IPOSIT=STOP,
     REACTION=YES INITIALS=NO FSINTERA=NO IRINT=DEFAULT CMASS=NO,
     SHELLNDO=AUTOMATIC AUTOMATI=OFF SOLVER=SPARSE,
     CONTACT-=CONSTRAINT-FUNCTION TRELEASE=0.00000000000000,
     RESTART-=NO FRACTURE=NO LOAD-CAS=NO LOAD-PEN=NO MAXSOLME=0,
     MTOTM=2 RECL=3000 SINGULAR=YES STIFFNES=1000.00000000000,
     MAP-OUTP=NONE MAP-FORM=NO NODAL-DE='' POROUS-C=NO ADAPTIVE=0
*
TIMESTEP NAME=DEFAULT
@CLEAR
500 0.200000000000000
@
*
KINEMATICS DISPLACE=LARGE STRAINS=SMALL PRESSURE=NO INCOMPAT=NO
*
MASTER ANALYSIS=STATIC MODEX=EXECUTE TSTART=0.00000000000000 IDOF=111,
     OVALIZAT=NONE FLUIDPOT=AUTOMATIC CYCLICPA=1 IPOSIT=STOP,
     REACTION=YES INITIALS=NO FSINTERA=NO IRINT=DEFAULT CMASS=NO,
     SHELLNDO=AUTOMATIC AUTOMATI=ATS SOLVER=SPARSE,
     CONTACT-=CONSTRAINT-FUNCTION TRELEASE=0.00000000000000,
     RESTART-=NO FRACTURE=NO LOAD-CAS=NO LOAD-PEN=NO MAXSOLME=0,
     MTOTM=2 RECL=3000 SINGULAR=YES STIFFNES=1000.00000000000,
     MAP-OUTP=NONE MAP-FORM=NO NODAL-DE='' POROUS-C=NO ADAPTIVE=0
*
AUTOMATIC TIME-STEPPING MAXSUBD=100000 ACCURACY=NO,
     DISTOL=0.00100000000000000 DTMAX=3.00000000000000,
     RESTORE=AUTOMATIC RESPS=NO RESFAC=0.000100000000000000,
     DIVFAC=2.00000000000000
*
MASTER ANALYSIS=STATIC MODEX=EXECUTE TSTART=0.00000000000000 IDOF=111,
     OVALIZAT=NONE FLUIDPOT=AUTOMATIC CYCLICPA=1 IPOSIT=STOP,
     REACTION=YES INITIALS=NO FSINTERA=NO IRINT=DEFAULT CMASS=NO,
     SHELLNDO=AUTOMATIC AUTOMATI=ATS SOLVER=SPARSE,
     CONTACT-=CONSTRAINT-FUNCTION TRELEASE=0.00000000000000,
     RESTART-=NO FRACTURE=NO LOAD-CAS=NO LOAD-PEN=NO MAXSOLME=0,
     MTOTM=2 RECL=3000 SINGULAR=YES STIFFNES=1000.00000000000,
     MAP-OUTP=NONE MAP-FORM=NO NODAL-DE='' POROUS-C=NO ADAPTIVE=0

zstqq

更进一步解释说明一下吧，共同探讨

zstqq

最好吧命令流附在文件里，
直接复制的里面干扰字符太多

tumu123

这是最原始的命令流文件
不存在什么干扰
土是摩尔库仑模型
桩是弹性
复制之后就可以用

zstqq

不是呀，复制粘贴后有很都乱码呀，你试试看，
你的命令流我都复制好几天了，到现在还没不乱码三掉呢

sunnylove

乱码是网页隐含的介绍，应该是网页的问题，不是楼主的问题，楼主压缩一下，传个压缩包吧

[ 本帖最后由 sunnylove 于 2008-4-18 11:47 编辑 ]

pantherprairie

楼主发了压缩包，对adina参数设置进行说明。另外对adina模拟接触进行探讨。在adina中模拟接触是否用的是接触单元，是Goodman单元？

pantherprairie

一楼的，你的接触面法向刚度和切向刚度是怎样设置？说明你的接触算法，不要发命令流，把你参数设置说明更好

wuzhiping0118

你是怎么确定土的参数的啊

冰棱无痕

楼主，你的附件再那里啊？我怎不么看不到

luoxh523

乱码太多了
没法用