用ANSYS建模，用程序导出ABAQUS的INP文件，再运行计算。
inp文件如下：
*HEADING
Structure analysis of Pushihe bifurcate pipe
SI UNITS(N,KG,M,S,Pa,kg/m)
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=NO,CONTACT=NO
*NODE,NSET=allnode
......
*Nset,NSET=xside
......
*Nset,NSET=Yside
......
*ELEMENT,TYPE=C3D8R,ELSET=ALLELEM
......
**===========================================================
**START MATERIAL DEFINITION
*SOLID SECTION,ELSET=ALLELEM,MATERIAL=ROCK
*MATERIAL,NAME=ROCK
*DENSITY
2700,0.0
*ELASTIC,TYPE=ISOTROPIC
200E9,0.3
**END MATERIAL DEFINITION
**===========================================
*Boundary
xside, 1, 1
xside, 2, 2
xside, 3, 3
*step
*static
*CLOAD
Yside,2,-5e+07
*restart,write,frequency=1,overlay
*endstep
另附ANSYS导入ABAQUS简单命令：
/COM,
==========================================================
/COM, Output ABAQUS/Standard input file
/COM,
====================================================
======
*CFOPEN,shench,inp !打开ABAQUS的输入INP文件
*VWRITE
('*HEADING') !项目名称
*VWRITE
('Structure analysis of Pushihe bifurcate pipe')
*vwrite
('SI UNITS(N,KG,M,S,Pa,kg/m)')
*VWRITE
('*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=NO,CONTACT=NO')
*VWRITE
('*NODE,NSET=allnode')
nsel,s,,,all !选择所有的节点
*GET,NNOD,NODE,,COUNT
NNUM = 0
*DO,I,1,NNOD,1
NNUM = NDNEXT(NNUM)
*VWRITE,NNUM,NX(NNUM),NY(NNUM),NZ(NNUM)
(F7.0,TL1,3(',',f16.9))
*ENDDO
*VWRITE
('**===========================================================')
*VWRITE
('*Nset,NSET=xside')
nsel,s,,,xside !选择边界XSIDE的节点
*GET,NNOD,NODE,,COUNT
NNUM = 0
*DO,I,1,NNOD-1,1
NNUM = NDNEXT(NNUM)
*VWRITE,NNUM
(F7.0,TL1,',')
*ENDDO
NNUM = NDNEXT(NNUM)
*VWRITE,NNUM
(F7.0,TL1,' ')
*VWRITE
('**===========================================================')
*VWRITE
('*Nset,NSET=Yside')
nsel,s,,,Yside !选择施加荷载节点YSIDE (为铅直向)
*GET,NNOD,NODE,,COUNT
NNUM = 0
*DO,I,1,NNOD-1,1
NNUM = NDNEXT(NNUM)
*VWRITE,NNUM
(F7.0,TL1,',')
*ENDDO
NNUM = NDNEXT(NNUM)
*VWRITE,NNUM
(F7.0,TL1,' ')
*VWRITE
('**===========================================================')
*VWRITE
('*ELEMENT,TYPE=C3D8R,ELSET=ALLELEM')
ESEL,S,,,ALLELEM !选择所有的单元
*GET,NELE,ELEM,,COUNT
ENUM = 0
*DO,I,1,NELE,1
ENUM = ELNEXT(ENUM)
NI = NELEM(ENUM,1)
NJ = NELEM(ENUM,2)
NK = NELEM(ENUM,3)
NL = NELEM(ENUM,4)
NM = NELEM(ENUM,5)
NN = NELEM(ENUM,6)
NP = NELEM(ENUM,7)
NQ = NELEM(ENUM,8)
*VWRITE,ENUM,NI,NJ,NK,NL,NM,NN,NP,NQ
(F7.0,TL1,8(',',F7.0,TL1),' ')
*ENDDO
*VWRITE
('**===========================================================')
!以下是ABAQUS基本命令模板
*VWRITE
('**===========================================') !两个**号在前为ABAQUS注释行
*VWRITE
('**START MATERIAL DEFINITION')
*VWRITE
('*SOLID SECTION,ELSET=ALLELEM,MATERIAL=ROCK')
*VWRITE
('*MATERIAL,NAME=ROCK')
*VWRITE
('*DENSITY')
*VWRITE
('2700,0.0')
*VWRITE
('*ELASTIC,TYPE=ISOTROPIC')
*VWRITE
('20000,0.3')
*VWRITE
('**END MATERIAL DEFINITION')
*VWRITE
('**===========================================')
*VWRITE
('*Boundary ')
*VWRITE
('*step')
*VWRITE
('*CLOAD')
*VWRITE
('*restart,write,frequency=1,overlay ') ! 为以后重启动分析作数据准备
*VWRITE
('*endstep')
*VWRITE
('*step')
*VWRITE
('*static') ! *VWRITE
('*endstep')
静力分析