—、引言
时程分析法是对结构动力方程直接进行逐步积分求解的一种动力分析方法。时程分析法将地震波按时段进行数值化后,输入结构体系的振动微分方程,采用直接积分法计算出结构在整个强震时域中的振动状态全过程,给出各个时刻各个杆件的内力和变形。现已成为多数国家抗震设计规范或规程的分析方法之一。
二、有限元软件ABAQUS简介
ABAQUS 是美国ABAQUS 公司(原名HKS 公司. 即Hibbitt , Karlsson & Sorensen, INC?2005年被法国达索公司收购,2007年公司更名为SIMULIA)。ABAQUS已成为国际上最先进的大型通用有限元力学分析软件之一。ABAQUS是一套功能强大的进行工程模拟的有限元软件。其解决问题的范围从相对简单的线性分析到许多复杂的非线性问题。ABAQUS拥有CAE工业领域最为广泛的材料模型,它可以模拟绝大部分工程材料的线形和非线形行为,可以进行结构的静态和动态分析,如应力、变形、振动、热传导以及对流等。也可以模拟广泛的材料性能,如金属、橡胶、塑料、弹性泡沫等,而且任何一种材料都可以和任何一种单元或复合材料的层一起用于任何合适的分析类型。
三、模型建立与求解
1、Part
Create Part: Name: Ban, 3D, Deformable, Shell , Planar,输入坐标创建一个18X9m 的壳部件,作为混凝土楼板部件;
Create Part: Name: Zhu, 3D, Deformable, W让e , Planar? 输入坐标创建一个长3m 线部件,作为柱部件;_ Create Part: Name: Liangs 3D, Deformable, Wire , Planar,输入坐标创建一个长6X3m,宽4. 5X2m的线网部件,作为梁网部件;
2、S ection
Create Material: Name: steel, General, Density 7800; Mechanical, Elasticity, Young* s Modulus 2.1ell, Poisson* Ratio 0.3;
Create Material : Name : concrete , General , Density 2500 ; Mechanical , Elasticity, Young' s Modulus 3el0, Poisson" Ratio 0.3。
Create Profile:分别创建梁剖面Liang Profile和柱剖面Zhu Pro,厚0. Im和0. 6X0. 6m,厚0. Imo
Create Section : 板截面:Name : Ban Section , Shell , Homogeneous , Shell Thickness Value: 0.1, Mat erial: concrete;梁截面:Name: Liang Section, Pro: Liang Pro name:steel;柱截面:Name: Zhu Section, Pro: Zhu Pro name: steelo
Assign Section:分别把创建好的板、梁和柱截面指派到板、梁、柱部件上,梁和柱截面指派时同时指派部件方向。
3、Assembly
Create Instance:选择板、梁、柱部件,用Dependent方式生成实体;将柱子旋转,使其与板和梁平面垂直;阵列柱子4X3二12根,将12根柱子分别移到梁网交点处,完成一层框架的构造;将构造好的一层框架阵列成3X2=6个框架层(阵列时,注意各个框架层间的距离应该稍大,这样利于平移框架),平移各层框架层,完成6层框架模型的构造。
4、Step
Create Step: Step-1, General, Static General,重力分析,结果作为时程分析数据; Step-2, General, Dynamic Implicit, Time period: 10。
5、I nteract i on
Create constraint: Tie,选择板、梁、柱接触点、面,绑定。
6、Load
Create Boundary Condition : Name : BCT , Step : Step-2 , Mechanical 、Displacement/Rotation,选择柱低端,约束住柱低端除X方向外的5个自曲度;
Create Boundary Condition : Name : BC-2 , Step : Step-2 , Mechanical , Acceleration/Angular acceleration,选择柱低端,确定;勾选Al:l (加速度单位m/s2), 在Amplitude 中创建一个Name:ac, Type:Tabular, Time span:Toal time,复制ElCentro 数据到Amplitude Data, OK, Amplitude右框下拉选择刚才建好的ac加速度,OK。
ft
7、M esh
Seed Part Instance;Assign Mesh Controls:Quard, Free, Advancing front;Assign Element Type: Standard? Linear, Reduced Integration, 壳单元S4R: Mesh Part: 最后Verify Mesh检查网格是否错误。
Create Job; Write Input (输出inp 文件),Submit;打JF Monitor 观察是否有错误、警报;点击Result进入可视化窗口Visualizationo
9、V i suaI i zat i on
Result, Step/Frame 分析步/帧,査看各时程变形情况;Tools, XY Data, Create, ODB
field output, plot顶层楼板位移时程曲线,各层楼板加速度时程曲线,各层楼板位移曲线。10、数据后处理
在Visualization, XY Data Manager中点击Edit,输出各层楼板位移数据到Excel中,利
用Excel:相邻两层的位移作差后得到层间相对位移,画出以Time为横轴、层间相对位移为纵轴画出各层层间位移时程曲线;然后层间相对位移除以层高得到层间位移角,画出以Time 为
横轴、层间位移角为纵轴画出各层层间位移角时程曲线。
四、结果展ZF
钢框架内力云图1
钢框架内力云图2
OR IMSS J>I ▼
z
z
钢框架内力云图3
smiF*W?W A
JtUfJi* 6
底端加速度时程曲线
aa
JC
aa
??
Trre
一层楼板加速度时程曲线
IJ9
7
?卩
Timo
24 3 44
Timo
s -^c b K
■%火rr ?zw ?,? A
? ?4“ ?X 、?>2t^r2 A
? ?“ — u ?;* I
C 5H <
----------m
AXJ ?t ? A>X*X WM" ? 1
>43W 5?】《C ?
V
U O =2J2Q %?:*1
?:?I Mr ?*>!-><
^i
W232」■ ?
u o e b g u
§S Q _O O O V
Tirrte
三层楼板加速度时程曲线
? ***>?>? 科
?
1 WM 5iyD2O *4J*< Time
ARM ?XSW —
四层楼板加速度时程曲线
AA1 ???EJgy W Mi P ??E*J7
?O&y
u o _@l?w
?k ?〉>u ?r ?&?o
底嵋与一层楼板层间位移
0.008 --------------------------------------------------------------------------------
0.006 --------------------------------------------------------------------------------
-0.008 -----------------------------------------------------------------------------
T&ae
底端与一层楼板层间位移
一启与二忘裕板尿间位轸
u o e b p o u p
顶层楼板加速度时程曲线
0. 002
YL 002
-0. 004
YU0&
gsy
s
.0
99【
-^.0
clclfs
lul
I
HSU 9
S 00
自.
y-so
?
gm 3SS 29
孩聲
3.0
01
01
片 9 7 9 T T Y>\ 9 J. i-
V- <71 p\CH C S G ; V
^v^VJL2 S l
riinrrriiiXiiiiiuiniiiii
9 M
M 2
0.2
0.025
-0.025 ----------------------------------------------------------------------------
TI K
三层与四层楼板层间位移
四忘弓砖樓板层何位侈
-0.015
Tine
二层与三层棲板层间位移
二层与三层楼板层间位移
三层与四层搂板层间位移
$53.
毎过?
ace aois
a 01 0.(X6 o -aixe HXOl -0.015 HX02
-0.01
LZZE .
/-
Tiwe
五层与顶层楼板层间位移
底瑞与一层楼板层间位移角
lilK
底端与一层楼板层间位移角
一层与二层楼板层间位移角
Tinr
0. 02
w u r o
c;
0. 015
0.01
0.005
0 -0.005
-0.01 -
0.015 -
0.02
0.0025 0.00Z 0.001B 0.001 0.0005 0 -0.0005 -0.001 -0.0015 -0.00Z -0.0025
0.01
0.008
0伽 O.OOQ
0.002
-0.002 -
acoa -
0.006 -
0.008
-0.01
9PS-S
i
368(
芒
篡
UJ g
891!
一层与二层楼板层间位移角
二层与三层楼板层间位移角
三层与四层棲板层间位移角
0.008 ----------------------------------------------------------------------------------
Tii??
三层与四层楼板层间位移角
四层与五层棲板层间位移角
0. 005 ---------------------------------------------------------------------------------
4 00Q ------------------------------------------------------------------------- -<1 00
5 --- --------------------------------------------------------------------
Ti?r
0.2
0.12 a co? a 002 a 001
?
-0.001 -0.002 -0.005
-4).008
Q.004 - ------------------------------------------------------------------------------- 1
TI K
五层与六层楼板层间位移角
五、总结与分析
通过分析数据:
① 该钢框架模型顶层楼板最大位移发生在9. 904615s 时,值为0.521279;根据《》:多层 柱顶位移H/5OO,对于该模型为18m/500-0. 036m,顶层楼板最大位移超出规范要求。
② 该钢框架模型最大加速度:底端发生在2. 13928s,值为3.40922, 一层楼板发生在 4. 53904s ,值为7. 02656,二层楼板发生在2.64224s,值为7. 58957 ,三层楼板发生在 3. 22835s ,值为8.61059,四层楼板发生在4.78149s,值为11.1631,五层楼板发生在 4. 78149s,值为12.4091,顶层楼板发生在4. 53904,值为13.9828。
③ 该钢框架模型层间位移:底端与一层楼板最大值为0. 006104, 一层与二层楼板最大值为 0. 025158,二层与三层楼板最大值为0.016052,三层与四层楼板最大值为0. 02146b 四层与 五层楼板最大值为0.012741,五层与六层楼板最大值为0. 016035;根据《》:多层层间位移限 值h/400,对于该模型为3m/400=0. 0. 0075m,最大层间位移超出规范要求。
④ 该钢框架模型层间位移角:底端与一层楼板最大值为0. 002035, 一层与二层楼板最大值 为0.008386,二层与三层楼板最大值为0. 005356,三层与四层楼板最大值为0.007154,四层 与五层楼板最大值为0.004247,五层与六层楼板最大值为0.003345;根据《》GB50011-2010 第5.5.1条规定:“地震作用下,多、弹性限值为1/250",最大层间位移角超出规范要求。
⑤ 由以上分析,顶层楼板位移、层间位移、层间位移角均超出规范要求,作者对此结果进 行了简单分析,造成此结果可能原因有:没有对时程地震波ElCentro 波进行系数调整,以满 足抗震规范;超过加速度时程的加速度没有按规范将其最大值调整为0. 35m/s 2:对于输入的地 震加速度时程曲线没有满足抗震规范,要满足地震动三要素:频谱特性、有效峰值和持续时间 要符合规定的要求。
⑥ 对于此钢框架结构,其抗侧刚度不足,应通过增加柱截面、梁高和梁宽等措施提高结构 抗侧刚度。
⑦ 以上分析可能不够准确,望读者批评指岀。
0. 004
0.002
-0.002
g F S
OJ
9
S 0O ,"
-0.004
I n
r r l/F
S989
&
参考文献
[1]王勖成,有限单元法?北京:清华大学出版社,2003
[2]马晓峰,ABAQUS6. 11有限元分析从入门到精通.北京:清华出版社,2013
[3]曹金凤等,ABAQUS有限元分析常见问题与解答.北京:机械工业出版社,2010
[4]徐珂.ABAQUS建筑结构分析应用.北京:中国建筑出版社出版社,2013
?0 ?0 ?0
??0
?0 ?0 ?0 ?0 ?0 ?0 ?0
?0
?
?
?0
?0 ?0 ?0 ?0 ?0 ?0 ?0
S000000UO666666UI
sooooooua .666666S&
sooooooz/moooops
SOOOOOOUO
?
一
866666
E .0
?
SOOOOOOUO
??6 ?0 ?
9f6666668
?0:G
?
946666668.0: 一
866666-
?
仃
?
9,6666668?9」
00000pr ? 9066666689
心666666Z.T
?0 ? 9Z66666689??0
? 9f66666689S666666UI
?9L6666668
0, ?666666S ?
E
? 9L66666689二 OOOOOPS
? 9f6666668.
0:一
866666
E
?
?9f6666668
?0:
?
6
?0 ?
9L6666668.0
?G
>0066666680=866666「仃 6,6666668.0
?一
00000p 「 6Z.6666668.0
.666666
』
仟 ?9L6666668.0 S666666UT
?0
?9(6666668 0??0
?9L6666668.0 S666666UI
?
9L6666668 0 ?666666」&
?9\6666668?0 二OOOOOPS
?
?0
?
?
?
0 .0 ?0
?0
?
9L6666668.0 ?一 866662
*
?
9,6666668 0
。?
6
SOOOOOOUIG
?
5S66666「仃
二
00000PB
?666666SV
?,666666;
?SOOOOOOUC 心
OOOOOOUC ?S00002S
doooooouc
心
OOOOOOLZ
Kooooooue Kooooooue ?
SOOOOOOLZ
?
?
0 心
OOOCWI.O
心
666666UI
?666666『E 二
OOOOOPS 二 866666 一
7
心
O0O0OOZ/I.6
s Tr .G
S
?
一 866666「仃
?第 二0000寻
S
?勺刁
?666666
『『 S 守
S666666U-7
r o
S666666U
一
?
666666
扶
-T00000PS
二
86666617.
VS GY 力&三二2!二龙乂 乂匚丝竺ER 舄门芯怜石E 谷WE A 気二衣滾滾E 杀处寻导等写耳寻孚壽铮孝
?0
即
6
?一
UPON 咲
UPHU 蛍
UEU ?
upd 咲
Evd
:
?
*
ONHgmuo 。.0NU 匕 o 一s 三?
ON==RPO 【
u ?
ONHOU02
?三^dald*
.u
?
9
uu v y s n b J v L*(q
P2EQUQD
*
一丄
Qpo-M Bunu 二
upow
XJZJZPWeDdusu
qof
*
*
三
50, 0., -2.70000005. 0. 51, -1.79999995, -2.70000005. 0. 52, ?3.5999999.-2.70000005, 0. 53, -5.4000001. -2.70000005, 0.
54, ■7.19999981
, -2.70000005. 0.
55, ?9?. -2.70000005. 0. 56, 9?,-4.5, 0.
57, 7」9999981. -4.5, 0. 58, 5.4000001. 45, 0. 59, 3.5999999, 45, 0. 60, 1.79999995. -4.5, 0?61, 0?,-4.5, 0.
62, ■1.79999995
, -4.5. 0.
63, -3.5999999. -4.5, 0.
64, -5.4000001. -4.5, 0.
65, ?7.19999981.-4.5, 0. 66, ?9?. -4.5, 0.
?Element, type=S4R
1, 1. 2.13.12
2, 2, 3,14.13
3, 3. 4. 15. 14
4, 4. 5,16,15
5, 5, 6.17, 16
6, 6. 7. 18. 17
7, 7, 8,19, 18
& & 9.20, 19
9. 9. 10,21,20
10, 10. 11,22,21
11,12, 13, 24. 23
12, 13. 14,25,24
13, 14. 15,26,25
14, 15. 16.27.26
15, 16. 17,28.27
16, 17. 18.29.28
17, 1& 19. 30.29
18, 19. 20.31.30
19,20,21,32,31
20,21.22.33, 32
21,23.24.35, 34
22, 24. 25. 36, 35
23, 25. 26.37, 36
24, 26. 27. 3& 37
25, 27. 28, 39. 38
26, 28. 29.40, 39
27, 29, 30.41,40
2& 30,31,42,41
29,31,32,43,42
30, 32. 33.44.43
31,34,35,46,45
32, 35. 36.47,46
33, 36. 37,48.47
34, 37. 38.49,48
35, 3& 39.50,49
36, 39.40,51,50
37,40,41,52,51
3& 41,42,53,52
39.42.43.54. 53
40.43.44.55, 54
41.45.46.57, 56
42.46. 47,5& 57
43.47.46 59, 58
44.46 49. 60, 59
45,49. 50.61.60
46. 50.51.62,61
47,51,52,63,62
4& 52. 53.64,63
49, 53. 54.65,64
50, 54. 55.66、65
*Nset, nset=Sel-k generate
1, 66, 1
*Elset, elset=Set-L generate
1, 50, 1
*Nset, nset=Sel-4. generate
1, 66, 1
*Elset, elset=Set-4, generate
1, 50, 1
*Nset, nset=Set-5. generate
1, 66, 1
?Elset, elset=Set-5, generate
1, 50, 1
** Section: Ban Section
*Shell Section, elset=Set-5. material=concrete 0.1,5
*End Part
*?
*Part. name=Liaiig
*Node
I. ?3?, 4.5. 0.
2, 3?, 4.5, 0.
3. 9?,
4.5, 0.
4, 9?,0
?,0.
5, 9?,-
4.5. 0.
6, 3?,-
4.5. 0.
7, ?3?,-
4.5. 0.
& ?9?,45. 0.
9. -9, 0
?,0.
10, -3?. 0., 0.
11, 3?,0
?,0.
12, ?9?. 4.5, 0.
13, -2.4000001. 4.5, 0. 14, -1.79999995, 4.5. 0. 15, .1.20000005, 4.5. 0.
16. -0.600000024? 4.5, 0. 17, 0?, 4.5. 0.
18, 0.600000024. 4.5, 0.
19, 1.2OOOOOO
5. 4.5, 0.
20, 1.79999995. 4.5, 0. 21, 2.400000 L 4.5, 0. 22, 3.5999999. 4.5, 0. 23, 4.19999981. 4.5, 0. 24, 4.80000019. 4.5, 0.
25, 5.4000001, 4.5, 0. 26, 6?, 4.5. 0. 27, 6.5999999, 4.5, 0. 28, 7」9999981. 4.5, 0. 29, 7.80000019. 4.5, 0. 30, &39999962. 4.5, 0. 31, 9?,3.9375. 0. 32, 9?, 3.375. 0.
33, 9?,2.8125, 0. 34, 9?, 2.25, 0. 35, 9?,1.6875, 0. 36, 9?, 1.125, 0. 37, 9?,0.5625. 0.
3& 9?,?0.5625, 0. 39, 9?,? 1.125, 0. 40, 9?,-1.6875, 0.
41, 9?,.7
*0 ? 7 . 0.
42, 9?,-2.8125, 0. 43, 9?,-3375, 0. 44, 9?,-3.9375, 0. 45, &39999962, -4.5, 0. 46, 7.80000019. -4.5, 0. 47, 7」9999981, -4.5, 0. 4& 6.5999999. -4.5. 0. 49, 6?,-4.5, 0. 50, 5.400000L -4.5. 0. 51, 4.80000019. -4.5, 0. 52, 4」9999981. -4.5, 0. 53, 3.5999999. -4.5. 0. 54, 2.400000 L -4.5. 0. 55, 1.79999995. 45. 0.
56, 1.2OOOOOO
5. -4.5, 0.
57, 0.600000024. 45, 0.
58, 0?,-4.5, 0.
59. -0.600000024. 45. 0.
60, -
1.2OOOOOO -4.5,
61, -1.79999995, -4.5,
62, -2.4000001. -4.5.
63, ?3.5999999.-4.5.
64, 419999981, -4.5,
65, -4.80000019, -4.5,
66, -5.4000001. -4.5.
67, -6., -4.5, 0. 68, ?6.5999999, -4.5. 0. 69, ?7.19999981,-4.5, 0. 70, -7.80000019, -4.5. 0. 71, ?8.39999962,-4.5, 0.
72. -9., -0.5625, 0. 73, ?9-1.125. 0. 74. ?9-1.6875, 0. 75, -9., 』?5 0. 76, -9., -2.8125, 0. 77, ?9-3.375. 0.
7& ?9-3.9375, 0.
79. -3.5999999 ?0., 0.
80. -4.19999981, 0., 0.
81, -4.80000019, 0?,0.
82, -5.4000001. 0?0.
83. -6., 0., 0.
84. -6.5999999 ?0.. 0. 85, -7.19999981, 0., 0.
86. -7.80000019, 0?0. 87, -8.39999962, 0., 0. 88, ?3-3.9375, 0. 89, ?3-3375, 0.
90. ?3-2.8125, 0. 91, ?3-? ?5 0.
92. ?3-1.6875, 0?
93. ?3-1.125. 0.
94. ?3.,-0.5625, 0?95, ?3.,0.5625. 0. 96. ?3 1.125, 0. 97, ?3., 1.6875. 0.
9& ?3 2.25, 0. 99. ?32.8125, 0. 100. -3., 3375, 0. 101. ?33.9375, 0. 102. 2.400000 L 0., 0. 103. 1.79999995, 0., 0. 104. 1.2OOOOOO5, 0?,0. 105. 0.600000024. 0.. 0. 106. 0?. 0?,0. 107. -0.600000024, 0…0.
10& -1.2OOOOOO5, 0?,0.
109. ?1.79999995,0?,0.
110, -2.4000001. 0?,0. 111. 3?,-3.9375, 0. 112, 3?,-3.375, 0. 113, 3?,-2.8125, 0. 114, 3?,225. 0. 115, 3?,-1.6875, 0. 116, 3?,?L125.0. 117, 3?,-0.5625, 0. 11& 3?,0.5625, 0. 119, 3?, 1.125, 0. 120. 3?, 1.6875, 0. 121. 3?, 2.25, 0. 122. 3?,2.8125, 0. 123. 3?, 3.375. 0. 124. 3?,3.9375, 0. 125. &39999962, 0?,0.
126. 7.80000019, 0?,0.
127. 7.19999981, 0?,0.
12& 6.5999999. 0., 0. 129. 6.. 0?,0.
130. 5.400000 L 0., 0.
131. 4.80000019, 0?0.
132. 4」9999981, 0?. 0.
133. 3.5999999. 0?,0. 134. ?9?. 3.9375. 0.
135. 1375, 0. 136. ?9?. 2.8125. 0.
137. 2.25, 0.
13& ?9?. 1.6875, 0.
139. 1.125, 0. 140. ?9?. 0.5625, 0.
141. ?8.39999962. 4.5, 0.
142. -7.80000019, 4.5, 0.
143. -7.19999981, 4.5, 0.
144. -6.5999999, 45 0.
145. -6., 4.5, 0.
146. -5.4000001, 4.5, 0.
147. -4.80000019, 4.5, 0.
14& 419999981. 4.5, 0.
149. -3.5999999, 45 0. ?Element. tvpe=B31
I. 1, 13
2.13, 14
3.1
4. 15
4, 15, 16
5,16, 17
6, 17, 18
7.1& 19
8, 19, 20
9,20,21
10,21. 2
II, 2,22
12, 22. 23
13, 23. 24
14, 24. 25
15, 25. 26
16, 26. 27
17, 27. 28
1& 28. 29
19, 29. 30
20, 30. 3
21, 3,31
22.31.32
23, 32. 33
24, 33. 34
25, 34. 35
26, 35. 36
27, 36. 37
2& 37. 4
29, 4. 38
30, 3& 39
31.39.40
32.40.41
33.41.42
34.42.43