多层平面钢框架连续倒塌仿真算例—尚晓江著作中的算例，运行通过，供大家学习参考

zhouqian1

免费的哦
梁柱截面为工字型钢，考虑底层中柱失效，模拟在失效点施加竖向荷载后的连续倒塌过程。

zhouqian1

以下为代码：
fini
/cle
/prep7

ET,1,BEAM161
KEYOPT,1,1,1
KEYOPT,1,2,2
KEYOPT,1,4,1,   
KEYOPT,1,5,2

ET,2,BEAM161
KEYOPT,2,1,1
KEYOPT,2,2,2
KEYOPT,2,4,2,   
KEYOPT,2,5,2

ET,3,BEAM161
KEYOPT,3,1,1
KEYOPT,3,2,2
KEYOPT,3,4,3,   
KEYOPT,3,5,2

*SET,_RC_SET,1,
R,1
RMODIF,1,1,0,0,0,0,0,   
RMODIF,1,6,0,0  
RMODIF,1,13,0,0,1   
RMODIF,1,16,0.3,0.016,0.3,0.01,0,0,
RMODIF,1,13,0   
RMODIF,1,14,0   

*SET,_RC_SET,2,
R,2
RMODIF,2,1,0,0,0,0,0,   
RMODIF,2,6,0,0  
RMODIF,2,13,0,0,1   
RMODIF,2,16,0.3,0.016,0.4,0.01,0,0,
RMODIF,2,13,0   
RMODIF,2,14,0   

*SET,_RC_SET,3,
R,3
RMODIF,3,1,0,0,0,0,0,   
RMODIF,3,6,0,0  
RMODIF,3,13,0,0,1   
RMODIF,3,16,0.3,0.016,0.5,0.01,0,0,
RMODIF,3,13,0   
RMODIF,3,14,0   

et,4,solid164

mp,ex,1,2.06e8
mp,dens,1,7.85
mp,nuxy,1,0.3
tb,plaw,1,,,1                                     !定义plastic kinematic模型
tbdata,,345e3,2.06e6,,,,0.01                      !定义屈服应力、切线模量及材料失效应变

k
k,,,,3.6
k,,,10,3.6
kgen,4,all,,,6
*dim,line,,7,3                                    !定义line数组
*dim,tr,,7                                        !定义tr数组
line(1,1)=1,4,7,10,2,5,8                          !第一列，直线起点
line(1,2)=2,5,8,11,5,8,11                         !第二列，直线终点
line(1,3)=3,6,9,12,3,3,3                          !第三列，直线方向点
tr(1)=2,3,3,2,1,1,1                               !tr数组赋值（单元类型及实常数）

*do,i,1,7
l,line(i,1),line(i,2)
latt,1,tr(i),tr(i),,line(i,3)
lesize,i,,,5
lmesh,i
*enddo

egen,6,100,all,,,,,,,,,,3.6
lclear,2                                         !去掉失效构件


block,-6,24,-6,6,-0.5
type,4
esize,3
vmesh,1
nummrg,node                                      !合并节点的目的：去掉刚才复制标准层时的重复点以及合并底面与框架边界的重合点

edpart,create                                    !创建part
edrd,d2r,4                                       !变形体转换为刚体

edcgen,assc                                      !定义自动单面接触及参数
edcontact                                        !定义附加接触参数

!施加位移边界条件
/solu
nsel,,loc,z,-0.5
d,all,all

!定义时间和荷载组
*dim,t,,2                                        !定义时间数组t，2行
*dim,q,,2                                        !定义荷载数组q，2行
*dim,f,,2                                        !定义荷载数组f，2行
*dim,t1,,3                                       !定义时间数组t1，2行
*dim,ft,,3                                       !定义荷载数组ft，3行
t(1)=0,1                                         !数组t赋值
q(1)=70,70                                       !数组q赋值
f(1)=2859.33,2859.33                             !数组f赋值
t1(1)=0,0.033,1                                  !数组t1赋值
ft(1)=0,-2859.33,-2859.33                        !数组ft赋值

!加载
esel,,type,,1                                    !选择全部梁单元
cm,eset,elem                                     !定义elem组元
nsel,,node,,13                                   !选择失效点
cm,nset,node                                     !定义nset组元
cm,nset1,node                                    !定义nset1组元
edload,,press,3,eset,t,q,1                       !施加均布荷载（动力松弛）
edload,,fz,,nset,t,f,1                           !施加集中力（动力松弛）
edload,,fz,,nset1,t1,ft                          !施加集中力（瞬态分析）

eddamp,,,0.576                                   !定义质量阻尼
eddamp,1,,0.000502                               !对part1定义刚度阻尼
eddamp,2,,0.000502                               !对part2定义刚度阻尼
eddamp,3,,0.000502                               !对part3定义刚度阻尼
time,1                                           !设置求解时间
edrst,,0.005                                     !设置结果文件输出间隔
edhtime,,0.005                                   !设置时程文件输出间隔
edopt,,,both                                     !设置输出文件类型

!求解
allsel
solve
fini

/post26
nsol,2,13,u,z,uz                                 !选择失效点的z向变形
prvar,2
plvar,2

zhouqian1

:):)个人觉得一个小问题就是梁柱不能能显渲染模型，比如ansys中用/eshape,1命令就能实现，但现在显示的仍为杆件，模型效果不大好，恳请高手指点。

蓝牙

很简单，toggle---》beam style 开启就可以了

cayman

楼主，我做这个例题的时候有个问题，就是nset，nset1加瞬态分析力和动力松弛力是一直提示不可以在同一个节点上重复加力，导致计算出错，为什么？

seabird08

egen,6,100,all,,,,,,,,,,3.63
lclear,2                                         !去掉失效构件

运行了这一段命令后，单元怎么显的那么乱啊

liuanxian

5# cayman

liuanxian

5# cayman
你需要两次分别选中相应的节点或节点组，分别定义或一次定义两个组元，然后才能进行分析。这样两次的荷载是加在两个不同的组元，应该就不会提示出错了。

wu870102

回复 4# 蓝牙


    lz能指导下toggle在哪吗？灰常感谢~~~

蓝牙

如图

jianqi169

回复 10# 蓝牙

好家伙，哈哈
受教

wu870102

回复 10# 蓝牙


    谢谢啦~~

yuyinghua

你好，这个连续倒塌的算例是尚晓江《Ansys/lsdyna动力分析方法和工程实例》中的例子吗，我怎么找不到啊？麻烦给我说一下吧

nan0419

谢谢楼主