登录  
 加关注
查看详情
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

为梦飞翔

我希望我身边的亲人和朋友天天快乐!

 
 
 

日志

 
 

框架—剪力墙结构分析  

2010-11-14 20:32:16|  分类: Ansys |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

本例来自于《ANSYS结构有限元高级分析方法与范例应用》尚晓江 邱峰 赵海峰 李文颖等编著,已经试算过,仅有几个小错误,下为已经修改的命令流(对此模型进行了自重作用下的结构分析、风荷载分析以及模态分析,对于地震反应分析,在这里就不作讲述了):

/TITLE,frame_shear_wall
/prep7 
et,1,beam4 
et,2,shell63
r,1,0.25,0.0052,0.0052,0.5,0.5          !底层立柱实常数
r,2,0.2025,0.0034,0.0034,0.45,0.45      !其它层立柱实常数
r,3,0.06,0.0002,0.00045,0.3,0.2         !梁实常数
r,4,0.25,0.25,0.25,0.25                 !剪力墙实常数
r,5,0.1,0.1,0.1,0.1                     !楼板实常数
mp,ex,1,3.0e10 
mp,nuxy,1,0.2  
mp,dens,1,2500

!创建底层关键点 
*do,i,1,1  
k,(i-1)*20+1,0,0,(i-1)*6
k,(i-1)*20+2,6,0,(i-1)*6
k,(i-1)*20+3,12,0,(i-1)*6  
k,(i-1)*20+4,17,0,(i-1)*6  
k,(i-1)*20+5,0,5,(i-1)*6
k,(i-1)*20+6,6,5,(i-1)*6
k,(i-1)*20+7,12,5,(i-1)*6  
k,(i-1)*20+8,17,5,(i-1)*6  
k,(i-1)*20+9,0,10,(i-1)*6  
k,(i-1)*20+10,6,10,(i-1)*6 
k,(i-1)*20+11,12,10,(i-1)*6
k,(i-1)*20+12,17,10,(i-1)*6
k,(i-1)*20+14,6,1.875,(i-1)*6  
k,(i-1)*20+16,6,3.125,(i-1)*6  
k,(i-1)*20+18,6,1.875,3.9  
k,(i-1)*20+20,6,3.125,3.9  
*enddo

!创建2-11层关键点
*do,i,1,10  
k,i*20+1,0,0,(i-1)*3+6
k,i*20+2,6,0,(i-1)*3+6
k,i*20+3,12,0,(i-1)*3+6
k,i*20+4,17,0,(i-1)*3+6
k,i*20+5,0,5,(i-1)*3+6
k,i*20+6,6,5,(i-1)*3+6
k,i*20+7,12,5,(i-1)*3+6
k,i*20+8,17,5,(i-1)*3+6
k,i*20+9,0,10,(i-1)*3+6
k,i*20+10,6,10,(i-1)*3+6
k,i*20+11,12,10,(i-1)*3+6
k,i*20+12,17,10,(i-1)*3+6
k,i*20+13,0,1.875,(i-1)*3+6
k,i*20+14,6,1.875,(i-1)*3+6
k,i*20+15,0,3.125,(i-1)*3+6 
k,i*20+16,6,3.125,(i-1)*3+6
k,i*20+17,0,1.875,(i-1)*3+6+1.5 
k,i*20+18,6,1.875,(i-1)*3+6+2
k,i*20+19,0,3.125,(i-1)*3+6+1.5
k,i*20+20,6,3.125,(i-1)*3+6+2
*enddo

!创建12层关键点
*do,i,11,11
k,i*20+1,0,0,(i-1)*3+6
k,i*20+2,6,0,(i-1)*3+6
k,i*20+3,12,0,(i-1)*3+6
k,i*20+4,17,0,(i-1)*3+6
k,i*20+5,0,5,(i-1)*3+6
k,i*20+6,6,5,(i-1)*3+6
k,i*20+7,12,5,(i-1)*3+6
k,i*20+8,17,5,(i-1)*3+6
k,i*20+9,0,10,(i-1)*3+6
k,i*20+10,6,10,(i-1)*3+6
k,i*20+11,12,10,(i-1)*3+6
k,i*20+12,17,10,(i-1)*3+6
k,i*20+13,0,1.875,(i-1)*3+6
k,i*20+14,6,1.875,(i-1)*3+6
k,i*20+15,0,3.125,(i-1)*3+6 
k,i*20+16,6,3.125,(i-1)*3+6
k,i*20+17,0,1.875,(i-1)*3+6+1.8
k,i*20+18,6,1.875,(i-1)*3+6+2.5
k,i*20+19,0,3.125,(i-1)*3+6+1.8
k,i*20+20,6,3.125,(i-1)*3+6+2.5
*enddo

!创建顶层关键点
*do,i,12,12
k,i*20+1,0,0,39.5
k,i*20+2,6,0,39.5
k,i*20+3,12,0,39.5
k,i*20+4,17,0,39.5
k,i*20+5,0,5,39.5
k,i*20+6,6,5,39.5
k,i*20+7,12,5,39.5
k,i*20+8,17,5,39.5
k,i*20+9,0,10,39.5
k,i*20+10,6,10,39.5
k,i*20+11,12,10,39.5
k,i*20+12,17,10,39.5
*enddo


!生成各层梁格
*do,i,2,13
L,(i-1)*20+1,(i-1)*20+2
L,(i-1)*20+2,(i-1)*20+3
L,(i-1)*20+3,(i-1)*20+4
L,(i-1)*20+5,(i-1)*20+6
L,(i-1)*20+6,(i-1)*20+7
L,(i-1)*20+7,(i-1)*20+8
L,(i-1)*20+9,(i-1)*20+10
L,(i-1)*20+10,(i-1)*20+11
L,(i-1)*20+11,(i-1)*20+12
L,(i-1)*20+5,(i-1)*20+9
L,(i-1)*20+6,(i-1)*20+10
L,(i-1)*20+3,(i-1)*20+7
L,(i-1)*20+7,(i-1)*20+11
L,(i-1)*20+4,(i-1)*20+8
L,(i-1)*20+8,(i-1)*20+12
*enddo

!生成立柱
*do,i,1,12
L,(i-1)*20+1,i*20+1
L,(i-1)*20+2,i*20+2
L,(i-1)*20+3,i*20+3
L,(i-1)*20+4,i*20+4
L,(i-1)*20+5,i*20+5
L,(i-1)*20+6,i*20+6
L,(i-1)*20+7,i*20+7
L,(i-1)*20+8,i*20+8
L,(i-1)*20+9,i*20+9
L,(i-1)*20+10,i*20+10
L,(i-1)*20+11,i*20+11
L,(i-1)*20+12,i*20+12
*enddo

!生成底层剪力墙开孔边框
*do,i,1,1
L,(i-1)*20+14,(i-1)*20+18
L,(i-1)*20+16,(i-1)*20+20
*enddo

!生成其它层剪力墙开孔边框
*do,i,2,12
L,(i-1)*20+13,(i-1)*20+17
L,(i-1)*20+15,(i-1)*20+19
L,(i-1)*20+14,(i-1)*20+18
L,(i-1)*20+16,(i-1)*20+20
*enddo

a,1,21,25,5         !底层外侧剪力墙  
*do,i,2,12
a,(i-1)*20+1,i*20+1,i*20+5,(i-1)*20+5,(i-1)*20+15,(i-1)*20+19,(i-1)*20+17,(i-1)*20+13
*enddo
a,14,2,22,26,6,16,20,18  !内侧剪力墙
*do,i,2,12
a,(i-1)*20+2,i*20+2,i*20+6,(i-1)*20+6,(i-1)*20+16,(i-1)*20+20,(i-1)*20+18,(i-1)*20+14
*enddo

                   !生成楼板
*do,i,1,12
a,i*20+1,i*20+2,i*20+6,i*20+5
a,i*20+2,i*20+3,i*20+7,i*20+6
a,i*20+3,i*20+4,i*20+8,i*20+7
a,i*20+5,i*20+6,i*20+10,i*20+9
a,i*20+6,i*20+7,i*20+11,i*20+10
a,i*20+7,i*20+8,i*20+12,i*20+11
*enddo

!对柱与梁划分网格
Lsel,s,,,181,192              !选择底层立柱   
latt,1,1,1                   
lesize,all,,,10               !设置网格划分数ndiv=10
lmesh,all

Lsel,s,,,192,324              !选择2-12层立柱   
latt,1,2,1                 
lesize,all,,,5              !设置网格划分数ndiv=5
lmesh,all

Lsel,s,,,1,180              !选择各层梁格  
latt,1,3,1                 
lesize,all,,,8              !设置网格划分数ndiv=8
lmesh,all


!对剪力墙划分网格
Lsel,s,,,372              !选择外侧剪力墙底端  
Lsel,a,,,413              !选择内侧剪力墙底端                 
Lsel,a,,,461              !选择内侧剪力墙顶端
lesize,all,,,8            !设置网格划分数ndiv=8

Lsel,s,,,418                
Lsel,a,,,420                            
lesize,all,,,6           

                           !选择剪力墙开孔顶端
Lsel,s,,,375
Lsel,a,,,379
Lsel,a,,,383
Lsel,a,,,387
Lsel,a,,,391
Lsel,a,,,395
Lsel,a,,,399
Lsel,a,,,403
Lsel,a,,,407
Lsel,a,,,411
Lsel,a,,,415
Lsel,a,,,419
Lsel,a,,,423
Lsel,a,,,427
Lsel,a,,,431
Lsel,a,,,435
Lsel,a,,,439
Lsel,a,,,443
Lsel,a,,,447
Lsel,a,,,451
Lsel,a,,,455
Lsel,a,,,459
Lsel,a,,,463
lesize,all,,,2

asel,s,,,1,24             !选择剪力墙
aatt,1,4,2                !设置剪力墙网格划分属性
amesh,all

asel,s,,,25,96            !选择剪力墙
aatt,1,5,2                !设置剪力墙网格划分属性
amesh,all

nsel,s,loc,z,0            !底部位移约束
d,all,all
/eshape,1
finish
!***************************************
!自重作用下的结构反应
/solu
antype,static
allsel,all
acel,0,0,9.81
solve
finish

!***************************************
!风荷载作用下的结构反应
/solu
antype,static
*dim,load_1,array,12        !定义荷载数组(风荷载属表面载荷,为简化起见,将等效风荷载施加在结构外表面的梁柱汇交节点上,单位kN)
*dim,load_2,array,12
*dim,load_3,array,12
*dim,load_4,array,12
*dim,load_A,array,12
*dim,load_B,array,12
*dim,load_C,array,12
load_1(1)=3.78,2.16,2.39,2.57,2.72,2.84,2.95,3.17,3.20,3.29,3.39,3.51
load_2(1)=7.56,4.32,4.78,5.14,5.44,5.68,5.90,6.34,6.40,6.58,6.76,7.02
load_3(1)=6.93,3.96,4.38,4.71,4.99,5.21,5.41,5.81,5.87,6.03,6.20,6.44
load_4(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93
load_A(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93
load_B(1)=6.30,3.60,3.98,4.28,4.54,4.74,4.92,5.28,5.34,5.48,5.64,5.86
load_C(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93

*do,i,1,12                !轴线1、2、3、4施加荷载
fk,20*i+1,fy,load_1(i)
*enddo

*do,i,1,12
fk,20*i+2,fy,load_2(i)
*enddo

*do,i,1,12
fk,20*i+3,fy,load_3(i)
*enddo

*do,i,1,12
fk,20*i+4,fy,load_4(i)
*enddo
solve

fkdele,all,all                 !删除轴线1、2、3、4施加荷载
 
*do,i,1,12                     !轴线A、B、C施加荷载
fk,20*i+1,fx,load_A(i)
*enddo

*do,i,1,12
fk,20*i+5,fx,load_B(i)
*enddo

*do,i,1,12
fk,20*i+9,fx,load_C(i)
*enddo
solve
finish

!***************************************
!模态分析
/solu
antype,model
modopt,lanb,6                     !采用Block Lanzcos方法提取模态,提取前6个模态
modopt,lanb,6,0,0,,off            !默认设置
solve
finish

  评论这张
 
阅读(435)| 评论(0)

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018