fini
/clear,all,all
/filname,gaodun
/prep7
/unit,si
SECTYPE,1, BEAM, HREC, www, 3
SECOFFSET, CENT
SECDATA,6.5,4,1.2,1.2,0.7,0.7
!!!!!!!!!!!
SECTYPE,2, BEAM, HREC, www, 3
SECOFFSET, CENT
SECDATA,6.5,1,1.2,1.2,0.3,0.3
!!!!!!!!!!!横系梁
/view,1,0,-1,0
et,1,beam188
mp,ex,1,3.50E10
mp,prxy,2,0.2
MP,DENS,2,2600 !pier材料属性
mp,ex,2,3.50E10
mp,prxy,2,0.2
MP,DENS,2,2600 !横系梁材料属性
k,1,-4
k,2,-4,,-43.5
k,3,-4,,-93.5
k,4,-4,,-143.5
k,5,4
k,6,4,,-43.5
k,7,4,,-93.5
k,8,4,,-143.5
/pnum,line,1
*do,i,1,3
l,i,i+1
*enddo
*do,i,5,7
l,i,i+1
*enddo
l,2,6
l,3,7 !墩系杆
!!!!!!!!!!!!!!!!!
lsel,all
lesize,all,5
type,1
secnum,1
mat,1
lsel,s,line,,1
lsel,a,line,,2
lsel,a,line,,3
lsel,a,line,,4
lsel,a,line,,5
lsel,a,line,,6
lmesh,all
type,1
secnum,2
mat,2
lsel,s,line,,7
lsel,a,line,,8
lmesh,all
allsel,all
nummrg,all,all
save
fini
!*********进行静力稳态求解*********!
err=0.01 !!容许误差!!
pmodi=1 !!活载系数!!
ss=0 !!迭代次数!!
*dowhile,err !!死循环控制设置,由下面的设置跳出!!
ss=ss+1
!!!!!!!!!!!!!!
/solu
antype,static
nsel,s,loc,z,-143.5
d,all,all !约束
acel,0,0,9.8 !重力加速度
sfbeam,all,2,pres,-3993.5 !横向风载
pstres,on
allsel,all
eplot
solve
fini
!进行屈曲分析
/solu
antype,buckle
BUCOPT,LANB,20,0,0
OUTPR,NSOL,ALL
solve
*get,freq1,mode,1,freq !!获得第一阶模态的屈曲荷载系数!!
*if,abs(freq1-1),le,err,then !!比较是否满足误差要求(le:≤)!!
*exit !!如果满足要求则退出循环!!
*else !!否则改变缩放系数,继续循环!!
pmodi=pmodi*freq1 !!活载系数!!
*endif
*enddo
finish
/SOLU !!对解进行扩展
EXPASS,ON
MXPAND,20,,,YES
OUTPR,NSOL,ALL
OUTRES
SOLVE
FINISH
!!!!!!!!!!!!后处理
/post1
set,first
*get,freq1,mode,1,freq
/EFACET,1
PLNSOL, U,SUM, 0,1.0
SET,NEXT
*get,freq2,mode,2,freq
/EFACET,1
PLNSOL, U,SUM, 0,1.0
SET,NEXT
*get,freq3,mode,3,freq
/EFACET,1
PLNSOL, U,SUM, 0,1.0
SET,NEXT
*get,freq4,mode,4,freq
/EFACET,1
PLNSOL, U,SUM, 0,1.0
fini
!!!!!!!!!!!!!非线性
/SOL
ANTYPE,0 !静力求解
NLGEOM,1 !打开大变形
SSTIF,ON !打开应力刚化
NEQIT,20 !指定平衡迭代的最大次数
eqslv,front,,0, !指定一个方程求解器,影响求解的速度和精度
outres,all,all !也可以使用outres,,all缺省方式
nsubst,500,,,1 !设置子步数
arclen,on !弧长法
sfbeam,all,2,pres,-freq1*3993.5 !至关重要一步
solve
fini
!!!!!!!!!!!!!绘制应力位移曲线
/POST26
FILE,'gaodun','rst','.'
/UI,COLL,1
NUMVAR,200
NSOL,2,15,U,X,tt !将要存贮点的数据内容
RFORCE,3,21,F,Y, !总的反力
!prod,4,3,,,,,,1,1,1 !给先前取值一个变量,如果以单位1开始,应乘以freq1
xvar,2
plvar,3
prvar,2,3
rforce,5,21,f,y
xvar,2
plvar,5
fini |