/prep7
!创建几何模型
*afun,deg
k,1,0,2,10
k,2,0,-2,10
k,3,0,-2,-10
k,4,0,2,-10
k,5,73,2,10
k,6,73,-2,10
k,7,73,-2,-10
k,8,73,2,-10
v,1,2,3,4,5,6,7,8
K,9,0,-1,-7
K,10,0,-1,7
K,11,73,-1,-7
K,12,73,-1,7
CIRCLE,9,0.13,11,12
CIRCLE,10,0.13,12,11
AL,17,18,19,20
AL,13,14,15,16
L,9,11
L,10,12
VDRAG,7,,,,,,22,,,,,
VDRAG,8,,,,,,21,,,,,
VSBV,1,2
VSBV,4,3
!定义单元类型
ET,1,BEAM188
ET,2,SOLID285
!定义梁截面(BEAM188)
SECTYPE,1,BEAM,CSOLID,
SECOFFSET,CENT
SECDATA,0.13 !定义杆件半径
!定义1号材料参数
MP,EX,1,70000.0
MP,NUXY,1,0.3
TB,SMA,1,,,MEFF
TBDATA,1,500, 311, 120, 8.3, 0.1, 70000
TBDATA,7,0
MP,EX,2,1.8
MP,NUXY,2,0.3
!划分网格
LSEL,S,LINE,,21,22
LATT,1,1,1 !给所选择的线定义属性
LESIZE,ALL,5
LMESH,ALL
VSEL,S,VOLU,,1
VATT,2,,2
MSHAPE,1,3D
ESIZE,1
VMESH,1
SAVE
!对底部三个节点施加三向位移约束:
/SOLU
! Load step 1:
NROPT,UNSYM
NLGEOM,ON
D,2005,ALL
D,2006,ALL
D,2030,ALL
D,2010,ALL
D,1,ALL
D,17,ALL
F,2,FX,40
F,18,FX,40
SOLVE ! Solution for load step 1
! Load step 2
D,2005,ALL
D,2006,ALL
D,2030,ALL
D,2010,ALL
D,1,ALL
D,17,ALL
F,2,FX,0
F,18,FX,0
SOLVE |