数值模拟分析实例.docx
- 文档编号:16190089
- 上传时间:2023-07-11
- 格式:DOCX
- 页数:50
- 大小:6.10MB
数值模拟分析实例.docx
《数值模拟分析实例.docx》由会员分享,可在线阅读,更多相关《数值模拟分析实例.docx(50页珍藏版)》请在冰点文库上搜索。
数值模拟分析实例
华中科技大学体育馆数值模拟分析
6.1分析模型的建立
采用有限元软件ANSYS建立该网壳结构有限元分析模型。
整体屋盖结构共计1481个节点,4430个单元,16种截面类型。
建模时,网壳结构主体结构部分(包括主桁架、次桁架、水平支撑和檩条)采用ANSYS的LINK8杆单元建模,两侧翼的主梁、次梁和支承钢管柱均采用BEAM4梁单元,网壳结构屋面下部混凝土支承结构亦采用BEAM4梁单元。
分析时,屋面板、设备管线等荷载等效为节点荷载,施加在结构节点上。
在网壳结构有限元分析中,对于杆件采用的LINK83-DSpar单元为三维单元,假设材料为均质等直杆,且在轴向上施加载荷,可以承受单向的拉伸或者压缩,每个节点上具有三个自由度,即沿X、丫和Z坐标轴方向。
该单元具有塑性、蠕变、应力硬化和大变形等功能,能较好的模拟三维空间桁架单元。
对于两侧翼结构和下部支撑体系的柱、梁等结构采用的BEAM4单元是一个轴向拉压、扭转和弯曲单元,每个节点有三个平动自由度和三个转动自由度,具有应力刚化和大变形功能。
施工过程模拟分析时考虑时,同时考虑温度效应影响,计算时材料假定为理想弹塑性材料。
图6-1有限元分析模型
6.2分析工况选取
按照实际施工顺序,将网壳结构屋盖施工过程划分为5个工况进行施工数值模拟,计算温度取为该阶段施工完成时的环境温度。
工况1:
7榀拱形主桁架安装完毕,但临时支撑未撤除,计算温度为温度15C;
图6-2工况1中屋盖结构平面图
工况2:
两侧翼结构安装完毕,完成后拆除其临时支撑,计算温度为8C;
图6-4工况2中屋盖结构平面图
(b)长轴立面
图6-5工况2中屋盖结构立面图
工况3:
次桁架、水平支撑及楼梯安装完毕,临时支撑拆除,计算温度为29C;
(a)短轴立面
(b)长轴立面
图6-6工况2中屋盖结构平面图图6-7工况2中屋盖结构立面图
工况4:
檩条及设备管线安装完毕,计算温度为41C;
(a)短轴立面
图6-8工况2中屋盖结构平面图
(b)长轴立面
图6-9工况2中屋盖结构立面图
工况5:
屋面板及保温层等安装完毕,计算温度为16C
6.3分析结果
图6-10所示为华中科技大学网壳结构屋盖施工过程模拟分析杆件应变云图从图中可以看出,在网壳结构施工过程中,受力较大杆件都集中在支座附近。
ELEKENTSOLUTION
(a)施工工况1网壳结构杆件应变图
ELEMENTSOLUTION
5TEP=1SUE-1TIJiE=l
LEPE1(NOAVG)
DMX=・031604
SHN=-.18IE-03
SKX^158E-03
-„ieiE-Q3-・S30E-04』・178E-04・753E-05a173E-03
71DE-04_・350E-04-r301E-05・SO^E-OS
ZHJ1
(b)施工工况2网壳结构杆件应变图
ELEMEHT50LUTTDN
371P<1
SUE«1
TIME=1
IEPEL(NOAVG)
DMX■・0D7S555W
S1K
(c)施工工况3网壳结构杆件应变图
(d)
ZHJ1
(d)施工工况4网壳结构杆件应变图
(e)施工工况5网壳结构杆件应变图
图6-10施工过程中网壳结构杆件应变图
7.1空心球节点有限元分析
7.1.1分析模型
采用有限元软件ANSYS对空间节点进行受力性能进行分析。
节点建模采用ANSYS提供的4节点弹塑性壳单元,每个节点有6个自由度,包括3个线位移自由度和3个转动自由度,该单元具有大挠度小应变的特点,能比较好的模拟球节点的受力特性。
划分网格时采用映射划分和自由划分、整体划分和局部加密相结合的方法。
先对模型分区,杆与球的主体采用映射网格划分法,相贯线处由于边界复杂采用自由划分。
网格划分时通过控制单元边长,并在相贯线附近缩小控制尺寸,从而保证计算精度。
这样,整个模型共分为56681个单元,28243个结
点。
加载时,考虑节点的实际受力情况加载,节点支座底部约束所有自由度,各杆上荷载按实际受力比例轴向加载。
计算模型见图43o
图43空心球节点计算模型
分析时考虑了材料非线性和几何非线性。
材料为Q345钢,计算时材料假定为弹塑性材料,服从Von-Mises屈服准则,材料塑性按双线性等向强化考虑,弹性模量E=2.0X05N/mm2,屈服后模量取2%E,fy=370N/mm2。
分析采用增量迭代的方法,将牛顿拉斐逊方法(Newton-Raphson)、线性搜索技术(LineSearch)应用预测(Predictor)、自适应下降(AdaptiveDescent等加速收敛技术有机结合建立非线性平衡求解方法。
收敛准则为位移及不平衡节点力收敛准则。
本文不计残余应力和节点区焊缝对节点极限承载力的影响。
7.1.2分析结果
在整个加载过程中,除7、8、9三根杆外,杆上应力呈弹性状态,从受荷点至球杆交界处,应力呈比例增长,由于是轴向加载,无偏心影响,杆周应力均匀分布,基本是轴向应力,环向应力很小,其中应力最大点位于最大受力杆与球的交界处。
从计算可知,在设计荷载和检验荷载作用下,除7、8、9三根杆外,其
余各杆均未屈服。
而在1.4倍设计荷载下,7、8、9三根杆相贯处和杆底端部分区域屈服。
球面受力状态比较复杂,属于空间受力范畴。
从计算得知,球杆交界处(焊
趾附近)应力集中现象比较严重。
在1.5倍设计荷载作用下,8号、9号杆件与空心球体交界处球体出现了局部屈服。
因此,根据上述有限元分析结果,该节点极限承载力可认为是1.5倍设计荷载。
在设计荷载和检验荷载作用下,对应于各杆轴向、支座肋板及索耳板的单向应变计测点应力分析结果见表8所示,对应于杆与球面交汇区及杆相贯处的应变花测点应力结果如表9所示。
节点的等效应力云图如图44~45所示。
NODALSOLUTION
STEP=1
OCT52006
22:
25:
01
.262E-04.842E+08・168E+09・2S3E+09.337E+09
・421E+08・126E+09・211E+09・295E+09・379E+09
NODALSOLUTION
STEP=1
SUB=4
TIHE=1
SEQV|AVGJ
DWZ=,0Q3S92
SUM=.262E-04
SHX-a379E+09
fc262E-04a842E+03+16SE+09B253E+09*337E+Q9
B421E40a126E+D9.211E+09,395E+09^379E+09
ST£P=1
SUB=4
TIHE=1
SEQV|AVGJ
DMZ=,0Q3S92
SHM=.262E-04
SHX-a379E+09
fc262E-04a842E+09+16SE+09B253E+09*3ME+Q9
B421E40S126E+D9.211E+09.295E4O9^379E+09
图44设计荷载作用下节点的单元等效应力图
NODALSOLUTION
STE?
-1
SUB»S
TIME=.84375SEQV(AVG)DMX-.008209SMH».320E-04SMX=.419E+09
OCT52006
23:
05:
41
.320E-04・931E+08・186E+09・2?
9E+09・372E+09
・466E+08・140E+09.233E+09・326E+09・419E+09
NODALSOLUTION
STE?
-1
SUB»S
TIME=.84375SEQV(AVG)DMX-.008209SMH».320E-04SMX=.419E+09
OCT52006
23:
05:
41
.320E-04・931E+08・186E+09・2?
9E+09・372E+09
・466E+08・140E+09.233E+09・326E+09・419E+09
MODALSOLUTION
STE?
-1
SUB»S
TIKE=.84375
SEQ7(AVG)
DMX-.008209
SMW».320E-04
SMX=.419E+09
.320E-04.931E+08・186E+09・2?
9E+09・372E+09
・466E+08・140E+09.233E+09・326E+09・419E+09
7.2铸钢节点有限元分析
7.2.1分析模型
采用有限元软件ANSYS对铸钢空间节点进行受力性能进行分析。
节点建模采用ANSYS提供的SOLID45单元,SOLID45单元为3-D固体结构单元,由八个节点组成。
在单元每一个节点上有三个自由度,即分别沿着三个坐标轴方向。
此单元可以进行塑性、蠕变、应力硬化、大变形以及大应变分析。
能比较好的模拟实体铸钢节点的受力特性。
划分网格时采用映射划分和自由划分、整体划分和局部加密相结合的方法。
这样,整个模型共分为61530个单元,13097个结点。
分析时考虑了材料非线性和几何非线性,按双线性等向强化考虑,采用Von-Mises屈服准则。
材料参数取弹性模量E=2.10X05MPa,屈服强度fy=375MPa,屈服后弹性模量取2%E,收敛准则为位移及不平衡节点力收敛准则。
加载时,考虑节点的实际受力情况,环向索三个方向的线位移均被约束,节点可以做竖向平面内的微小转动。
A1、A2、A3方向上荷载取试验中的荷载值轴向加载,计算模型见图65。
图65节点计算模型
7.2.2分析结果
在整个加载过程中,索夹板及索耳板的应力均呈弹性状态,由应力分析图可见应力较大区域出现在索耳板上及叉耳孔周围。
从计算可知,在设计荷载和检验荷载作用下,节点未屈服。
.在设计荷载和检验荷载作用下单元的等效应力如图65和图66所示。
NODALSOLUTION
SUN=・00306
SMX=4141.599
NODALSOLUTIONSTEP“
SUB-1
TIHE=1
SEQV(AVG)
DMX«.394E-06
SHN=.00386
SMX-441.599
.0038698.136196.269294.401392.533
49.07147・202245・33S343・467441・599
1
NODALSOLUTION
STEP"SUB»1TIME=1
OCT42006
09:
21:
18
SEQV(AVG)
DMX-.394E-06
SHN=.00386
SMX-441.S99
294.401392.533
343.46*7441・599
NODALSOLUTION
STEP=1
SUE-1
OCT4Z006
09£21:
IB
TIHE=1
SEQV(AVG)
.00386
93.136196.2S9294.401392.533
14^,202245+335343.467441,599
图65设计荷载作用下节点的单元等效应力图
NODALSDLOTION
+006173134.429268+852403,275537,698
6-7,218201.641336.06447Qa-IS7604.91
HODALSOLUTION
STE?
-1
SUB-1
OCT4200609:
47:
22
TIKE=1
SEQ7(AVG)DMX-.380E-06SMM».006173
SMX=604.91
・006173134・429268・8S2403・275537・698
67・218201.641336・064470・487604・91
HODALSOLUTION
STEP-1
SUB-1
OCT42006
09:
47:
22
TIME“
SEQV(AVG)DMX-.380E-06SMM».006173
SMX=604.91
・006173134・429268・852403・27S537・698
67・218201.641336・064470・487604・91
OCT42006
09:
47:
22
BQDALSQiLUTIOIir
5TEP-1
SUE■!
T1ME=15EQV(AVG)
DMX■・3BOE-OS
SMIT二.006JL73
51K=604.91
+006173134.429268+352^03,275537.698
67,216Z01+641336.06447Qa487€04.91
图66检验设计荷载作用下节点的单元等效应力图
云阳县G25公路高架桥
附录一有限元分析结果
STEP-1
NOV262006
SUB=1
TIME-1
1.147-.836
-.991363
・525264■・214S28
.096208
.680632
.369396
.05916
工况1下桥梁挠度分析结果
1
NODALSOLUTION
AN
STEPwl
SUE-1
TIHE-1
WOV2620D6
23:
54336
EPT01(AV&)
DMX=1-14S3J®二-・55ZE-05
--5«£E-Dfi.3BSE-Ci5.133E-O-4.2£7E-0^-321E-»<
-.5L1Z-06.iaOE-04-n4E-®i4山晡EHCH
工况1下桥梁应变分析结果
HOV262006
23:
51:
25
NODALSOLUTION
远P二丄
SUB=1
-,413961--294012-.174062-+O5-4L12-065838
--353987--234037-.114087-0<05863-125812
工况2下桥梁挠度分析结果
NODALSOLUTION
AN
mp=iSUB=1TIHE^l
BOV272006
00:
03:
55
EPT01(AVG)
3fi9E" -■畫II百C-Q5 J4^E-g5: -70JE-qS .iafrE-q« .IfQE-Qf .52PE-05 aJZFE-04 -1.447-1,065-.683301-+30i291-0807L9 -1,256-,87006-.492296-.110286-271723 工况3下桥梁挠度分析结果 NODALSOLUTION SUB=1 TIKE^l EPT01(AVG]DMX■丄.449 SJffl=-.591E-05 5MX=.309E-04 AN HOV262006 22: 54: 45 -・IflEE-QS Bip«E-g«a・缶百be■©电 -53fiE-Q5・1451: -01・JZ7E-04・ NODALSOLUTION HOV272006 -2,157-1,623-1,088-.553158-.018397 -1,89-1,355-.820538-.285778-2489S3 工况4下桥梁挠度分析结果 NODALSOLUTION AN SUB=1 TIHE^l BOV27200600: 07: 02 EPT01(AVG) eJLJjEE-D4 -44&E-q^a "口—CM .338E-05 ・19 -5Z7E-04 SUB=1 BOV272006 00: 10: 34 一.65937一.62573一・39209一.15645・07519 一.14255一.50891一.27527-.04163・i9201 工况5下桥梁挠度分析结果 NODALSOLUTION AN mp=i SUB=1 BOV272006 00: ID: 51 TIHE=1 EPT01(AVG) -・iqiSE-flifi Blg? E-g«a・gfiDE-q< -54ZE-Q5・14JE-01・fZ! E-ID4・2fl? E-04 NODALSOLUTION BOV272006 00: 12: 11 --314115-.220632-.127149-+033666-0598L7 --267374-U73891-.080408-013075-106558 工况6下桥梁挠度分析结果 工况7下桥梁挠度分析结果 NODALSOLUTION AN mp=i SUB=1 BOV272006 00: 13: 51 EPT01(AVG) DHX■丄.15S SJffl=-.520E-05 -.151E-Q5 J.S5FI: -! g4・£43E-q4 -5S5E-Q5・13tE-M・fOFE-04・踽呃-阳 NODALSOLUTION SUB=1 BOV272006 00: 14: 36 -lfe70€H? 6 -1,4? 1 845708-,415746.0142L6 -1,061-.630727-.200765-229197 工况8下桥梁挠度分析结果 NODALSOLUTION BOV272006 00: 14: 57 .SS9E-04 .4.42E-0$^SSSE-O^・ .2£l£-0$・IfiTE-Ofl NODALSOLUTiaJI AN STEP-lSUBsiriKE-1UY RSYS-0DMX-1SEW*EffK=., -1.Q15^^743433471601-U99764.072073 -.079356-.607519-.335682-.063046・207991 边跨跨中偏载标准荷载下桥梁挠度分析结果 HODALSOliirriON AN STEP«1 SUB=1 TIME=-1 EPTD1(AVG) DMX=1.017 SHU=-・555E-05 SMXB.340E-D4 N0¥272006 -120E-04.2&Sr-04.4S6E-04 .7I52E-05-Jj64E-04.252E-04.340E-04 NODALSOLUTION AN -1,966-1,478 -1,722 BOV272006 00: 09: 16 989962-+502078014195 -1,234-.74602-.258136-229747 中跨跨中偏载标准荷载下桥梁挠度分析结果 NODALSOLUTION SUB=1 TIKE^l EPT01(AVG]DHX■丄.967 SJffl=-.451E-05 51K=.608E-04 AN HOV272006 00: 09: 42 ■Z75E-05 附录二部分试验照片 试验桥梁 测试工程师正在指挥加载 测试工程师正在进行应变测试 测试工程师正在进行振动测试 安装就位的拾振器 8华中科技大学体育馆有限元分析 8.1六杆节点有限元分析 8.1.1分析模型 采用有限元软件ANSYS对空间节点进行受力性能进行分析。 节点建模采用ANSYS提供的4节点弹塑性壳单元,每个节点有6个自由度,包括3个线位移自由度和3个转动自由度,该单元具有大挠度小应变的特点,能比较好的模拟薄壁钢管的受力特性。 划分网格时采用映射划分和自由划分、整体划分和局部加密相结合的方法。 先对模型分区,主管与支管的主体采用映射网格划分法,相贯线处由于边界复杂采用自由划分。 网格划分时通过控制单元边长,并在相贯线附近缩小控制尺寸,使相贯处单元尺寸与主管厚度相当,从而保证计算精度。 这样,整个模型共分为7443个单元,7439个结点。 加载时,考虑节点的实际受力情况,3号主管端三个方向的线位移均被约束,1号主管端加滑动约束并轴向加载,支管均沿轴向加载,各管上荷载按实际受力的比例加载,其中1、3、4、5、6各管 所受的荷载比例为: -2.82: -0.61: 1: 1.42: 0.39(负号表示受压力作用)。 计算 模型见图40 AN NOV2SZOO4 Zl: 03: 55 图40六杆节点计算模型 分析时考虑了材料非线性和几何非线性。 材料为Q345钢,计算时材料假定 为理想弹塑性材料,服从Von-Mises屈服准则,材料塑性按双线性等向强化考虑,材料参数取材性试验得到的数值,E=1.99X105N/mm2,fy=369N/mm2。 分析采用增量迭代的方法,将牛顿拉斐逊方法(Newton-Raphson)、线性搜索技术(Line Search)应用预测(Predictor)、自适应下降(AdaptiveDescent)等加速收敛技术有机结合建立非线性平衡求解方法。 收敛准则为位移及不平衡节点力收敛准则。 SUB=30 TIMH=・75 SEQV(NOAVG) DMX^3.0: 97 SMM=15.355 SMX-439・CQ7 15.355109.502203.649297.797391.944 62.429156,576250.723344.67439.017 图42检验荷载作用下节点的单元等效应力图 8.2四杆节点有限元分析 8.2.1分析模型 用有限元软件ANSYS对该支座节点的受力性能和变形特征进行了弹塑性分 析并得出其极限承载力。 由于该球体壁厚和球径之比在工程薄壳范围内,故采用 四边形弹塑性板壳单元,每个节点有6个自由度,包括3个线位移和3个转动自由度,这种单元具有大挠度小应变的特点,能比较好的模拟薄壁钢管及薄壁球壳的的受力特性,分析时考虑了材料非线性和几何非线性,按双线性等向强化考 虑,采用Von-
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数值 模拟 分析 实例