完整建模屈曲分析命令流

goldenbrid... 发表于 2009-4-2 22:17:22 | 显示全部楼层 | 阅读模式
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
精彩评论倒序浏览

4142查看1评论

xtddpla 发表于 2012-4-6 22:29:05
楼主是特征值屈曲分析呢还是非线性屈曲分析?以及是对一个什么样的结构进行分析,说明一下我们才好决定要不要下载啊。
举报 回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

 
 
  • QQ:56984982
  • 点击这里给我发消息
    电话:13527553862
    站务咨询群桥头堡站务咨询桥梁专业交流群:
    中国桥梁专业领袖群
    工作时间
    8:00-18:00