恒载加活载的特征值屈曲分析.docx
- 文档编号:10199588
- 上传时间:2023-02-09
- 格式:DOCX
- 页数:21
- 大小:335.92KB
恒载加活载的特征值屈曲分析.docx
《恒载加活载的特征值屈曲分析.docx》由会员分享,可在线阅读,更多相关《恒载加活载的特征值屈曲分析.docx(21页珍藏版)》请在冰豆网上搜索。
恒载加活载的特征值屈曲分析
工程上有时会遇到既有恒载 (固定不变的载荷) 又有活载 (按固定比例变化的载荷) 的特征值屈曲分析问题,比如高层建筑在自重和风载的共同作用下会发生侧弯,其受压严重的区域有发生屈曲的可能。
这里的自重就是恒载,风载就是活载。
如果做特征值屈曲分析 (以下简称屈曲分析),由于它假定了所有载荷都是按相同的比例变化的,而既有恒载又有活载的结构,恒载是不变的,这种问题应该如何进行屈曲分析呢?
对这类问题,其屈曲分析的过程是一个迭代的过程,基本步骤如下:
1. 如果单独恒载有可能造成结构的屈曲,应该先单独施加恒载并进行屈曲分析。
只有在单独恒载作用不会使结构发生屈曲的前提下,才需要进行恒载+活载的屈曲分析。
2. 首先施加给定的所有载荷,进行屈曲分析,求出屈曲载荷因子 (最小的正特征值);
3. 恒载不变,只将活载部分除以所得的屈曲载荷因子,然后再一次进行屈曲分析,得到新的载荷因子。
4. 重复第 3 步,直到屈曲载荷因子基本等于 1,这时的活载就是在存在恒载的情况下,结构发生屈曲时的活载的大小。
下面是一个简单的例子(结构可能不太合理,只为说明问题):
一个井架,由 4 根纵向长梁和若干水平的或斜的短梁组成,形状如图 1所示。
纵梁的截面为 30*30mm2,短梁的截面为 10*10mm2。
材料性能为:
弹性模量 E=201000Mpa
泊松比 μ=0.3
密度ρ=7.8e-9*2 (见下面关于恒载的说明)
使用 BEAM188 梁单元划分网格,单元长度取 50mm。
建模过程可见附录中命令流的几何建模部分,这里不再详述。
约束井架 4 根长梁底部的所有位移自由度。
考虑两种载荷:
(1) 恒载
井架以及固定设备的自重。
建模时不考虑设备,通过加大各杆的材料密度将设备重量分配到井架的各杆上。
这里为了简单,只定义了一种材料,即各根杆是按照原有重量等比例的放大。
如果想要对不同的杆按不同比例放大,可以定义多种材料,分别给予不同的密度。
本模型中,将材料密度放大了一倍。
(2) 活载
活载为 4 根纵梁顶端的纵向力载荷。
第一种情况是各纵梁承受相同的压力;第二种情况是两边的纵梁分别承受绝对值相同的纵向拉力和压力载荷,构成一个弯矩载荷。
分析过程如下:
(1) 只有恒载时的屈曲分析
首先对只有恒载的情况进行一次屈曲分析,以确认只有恒载不会使结构发生屈曲,否则后续分析就没有意义了。
特征值屈曲分析分为两步,第一步是在恒载作用下,执行考虑预紧力的静力分析;第二步是在静力分析基础上,执行特征值屈曲分析。
a. 静力分析的Uz结果如下:
转换到常用后处理,显示Uz的结果,变形放大系数取20:
b. 然后进行屈曲分析,重新进入Solution模块,选择求解类型为EnginBuckling,设置求解的特征值数为3,然后求解。
求解完成后,到常用后处理,首先点击:
ResultsSummary显示特征值结果如下:
***** INDEXOFDATASETSONRESULTSFILE *****
SET TIME/FREQ LOADSTEP SUBSTEP CUMULATIVE
1 16.150 1 1 1
2 16.275 1 2 2
3 16.275 1 3 3
三个特征值都是正值,其中最小的一个就是重力作用下的屈曲载荷因子。
由于该值远大于1,说明仅有重力作用是不会发生屈曲的。
读取第一组结果,显示Usum–屈曲形状,变形放大系数取1500:
(2) 恒载与活载联合作用的屈曲分析,之一
除自重外,在井架四根纵梁的顶端施加相同的压力。
第一次各施加1000N的压力,然后进行屈曲分析。
同样分为两步,其结果如下:
a. 每根纵梁受压1000N时的静力分析结果
每根纵梁承受压力1000N时的屈曲分析结果
求解完成后,到常用后处理,点击:
ResultsSummary 显示特征值结果如下:
***** INDEXOFDATASETSONRESULTSFILE *****
SET TIME/FREQ LOADSTEP SUBSTEP CUMULATIVE
1 7.0822 1 1 1
2 7.1513 1 2 2
3 7.1513 1 3 3
与第一个特征值7.082对应的屈曲形状如下:
但是,由于恒载的存在,对活载来说,7.082并不是实际的屈曲载荷因子。
为了计算实际的屈曲载荷,需要调整压力值,直到特征值近似等于1。
将压力值乘上第一个特征值,作为新的活载(近似取7000N),再次进行屈曲分析。
所得结果如下。
c. 每根纵梁受压7000N时的静力分析结果
d. 每根纵梁受压7000N时的屈曲分析结果
求解完成后,到常用后处理,点击:
ResultsSummary显示特征值结果如下:
***** INDEXOFDATASETSONRESULTSFILE *****
SET TIME/FREQ LOADSTEP SUBSTEP CUMULATIVE
1 1.6027 1 1 1
2 1.6256 1 2 2
3 1.6256 1 3 3
与第一个特征值1.6027,对应的屈曲形状如下:
图7 自重加上Z向压力7000作用下的屈曲形状Usum
同样,由于恒载的存在,对活载来说1.6027 并不是实际的屈曲载荷因子。
为了计算实际的屈曲载荷,需要调整压力值,直到特征值近似等于1。
将压力值7000乘上第一个特征值,作为新的活载(近似取7000*1.65=11550N),再次进行屈曲分析。
所得结果如下。
e. 每根纵梁受压11550N时的静力分析结果
f. 每根纵梁受压11550N时的屈曲分析结果
求解完成后,到常用后处理,点击:
ResultsSummary 显示特征值结果如下:
***** INDEXOFDATASETSONRESULTSFILE *****
SET TIME/FREQ LOADSTEP SUBSTEP CUMULATIVE
1 1.0085 1 1 1
2 1.0245 1 2 2
3 1.0245 1 3 3
与第一个特征值1.0085,对应的屈曲形状如下:
由于特征值1.0085已经非常接近1.0,因此可以停止计算,以11550N
(或11550*1.085=11648N) 作为发生屈曲时。
单根纵梁的压力载荷。
(2) 恒载与活载联合作用的屈曲分析,之二
除了自重以外,在井架四根纵梁的顶端施加两对绝对值相同,方向相反的力载荷,对井架形成弯曲载荷。
可以利用上面的结果,第一次仍各施加1000N的压力或拉力,然后进行屈曲分析。
同样分为两步,其结果如下:
a. 每根纵梁受拉力或压力1000N时的静力分析结果
b. 每根纵梁受拉力或压力1000N时的屈曲分析结果
求解完成后,到常用后处理,点击:
ResultsSummary 显示特征值结果如下:
***** INDEXOFDATASETSONRESULTSFILE *****
SET TIME/FREQ LOADSTEP SUBSTEP CUMULATIVE
1 6.9829 1 1 1
2 10.490 1 2 2
3 14.170 1 3 3
与第一个特征值6.9829对应的屈曲形状如下:
同样,由于恒载的存在,因此6.983并不是实际的屈曲载荷因子。
将压力值乘上第一个特征值,作为新的活载(近似取7000N),再次进行屈曲分析。
所得结果如下。
c.每根纵梁承受拉力或压力7000N时的静力分析结果
d. 每根纵梁承受拉力或压力7000N时的屈曲分析结果
求解完成后,到常用后处理,点击:
ResultsSummary显示特征值结果如下:
***** INDEXOFDATASETSONRESULTSFILE *****
SET TIME/FREQ LOADSTEP SUBSTEP CUMULATIVE
1 1.5656 1 1 1
2 1.8046 1 2 2
3 3.2115 1 3 3
与第一个特征值1.5656,对应的屈曲形状如下:
同样,将压力值7000乘上第一个特征值1.5656,作为新的活载(近似取7000*1.6=11200N),再次进行屈曲分析。
所得结果如下。
e. 每根纵梁承受拉力或压力11200N时的静力分析结果
f. 每根纵梁承受拉力或压力11200N时的屈曲分析结果
求解完成后,到常用后处理,点击:
Res ultsSummary 显示特征值结果如下:
***** INDEXOFDATASETSONRESULTSFILE *****
SET TIME/FREQ LOADSTEP SUBSTEP CUMULATIVE
1 1.0141 1 1 1
2 1.1418 1 2 2
3 2.0804 1 3 3
与第一个特征值1.0141,对应的屈曲形状如下:
由于特征值1.014已经接近1.0,因此可以停止计算,以11200(或11200*1.014=11357)N作为在重力和弯曲载荷作用下,发生屈曲时,单根纵梁的拉力或压力载荷。
fini
/clear
/filn,buckling_test
!
*
WPSTYLE,,,,,,,,0
!
*
/PMETH,OFF,0
KEYW,PR_SET,1
KEYW,PR_STRUC,1
!
* 变量 case_i用于不同工况,分别选择要施加什么载荷
!
case_i=1 !
Gravityonly
!
case_i=2 !
GravityonlyandForcesFz
case_i=3 !
GravityonlyandForcesFzformed Moment
!
*
f_scale2=1 !
计算工况1;活载-force的放大因子,将每次的第一个特征值连乘,初始值1
!
最终 f_scale2=1*7*1.65
f_scale3=1 !
计算工况2;活载-force的放大因子,将每次的第一个特征值连乘,初始值1
!
最终 f_scale3=1*7*1.6
!
* f_scale2用于case_i=2; f_scale3用于 case_i=3;
!
*
/PREP7
!
*
ET,1,BEAM188
!
*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,201000
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7.8e-9*2 !
计算自重需要密度。
考虑设备等的重量,将密度值放大2倍
!
*
SECTYPE, 1,BEAM,RECT,50*50,0
SECOFFSET,CENT
SECDATA,30,30,5,5,0,0,0,0,0,0
!
*
SECTYPE, 2,BEAM,RECT,10*10,0
SECOFFSET,CENT
SECDATA,10,10,4,4,0,0,0,0,0,0
!
*
/UI,BEAM,OFF
!
*
k,1,-1000,-1000,0
k,2,1000,-1000,
!
*
FLST,3,2,3,ORDE,2
FITEM,3,1
FITEM,3,-2
KGEN,3,P51X,,,,,2500,,0
!
*
/AUTO,1
/VIEW,1,,-1
/ANG,1
/REP,FAST
!
*
l,1,3
l,3,5
l,2,4
l,4,6
l,3,4
l,5,6
l,1,4
l,2,3
l,3,6
l,4,5
!
*
FLST,2,4,4,ORDE,2
FITEM,2,7
FITEM,2,-10
LOVLAP,P51X
!
*
LPLOT
numcmp,all !
以上生成了一片井架
!
*
csys,1 !
到圆柱坐标系中,准备对一片井架进行拷贝
!
*
FLST,3,14,4,ORDE,2
FITEM,3,1
FITEM,3,-14
LGEN,4,all,,,,90,,,0 !
在圆柱坐标系中,将一片井架沿圆周方向拷贝4份,围成一个井架
!
*
NUMMRG,KP,,,,LOW !
合并KP点(以及线、面等),生成完整的井架
NUMCMP,ALL
!
*
/VIEW,1,,,-1
/ANG,1
lplot
csys,0
!
*
FLST,5,4,4,ORDE,4
FITEM,5,3
FITEM,5,-4
FITEM,5,15
FITEM,5,-16
LSEL,,,,P51X
LATT,1,,1,, 15,,1
!
*
FLST,5,4,4,ORDE,4
FITEM,5,1
FITEM,5,-2
FITEM,5,27
FITEM,5,-28
LSEL,,,,P51X
LATT,1,,1,, 10,,1
!
*
FLST,5,10,4,ORDE,2
FITEM,5,29
FITEM,5,-38
LSEL,,,,P51X
LATT,1,,1,, 7,,2
!
*
FLST,5,10,4,ORDE,2
FITEM,5,5
FITEM,5,-14
LSEL,,,,P51X
LATT,1,,1,, 17,,2
!
*
FLST,5,10,4,ORDE,2
FITEM,5,17
FITEM,5,-26
LSEL,,,,P51X
LATT,1,,1,, 19,,2
!
*
FLST,5,10,4,ORDE,2
FITEM,5,39
FITEM,5,-48
LSEL,,,,P51X
LATT,1,,1,, 12,,2
!
*
allsel,all
LESIZE,ALL,20,,,,1,,,1,
!
*
LMESH,all
!
*
/VIEW, 1, 0.354987937241 , 0.874976631778 ,-0.329240729948
/ANG, 1, -60.9629091647
/auto,1
/REPLO
!
*
ALLSEL,ALL
!
*
/NUMBER,1
/PNUM,SECT,1
EPLOT
!
*
FLST,2,4,3,ORDE,4
FITEM,2,1
FITEM,2,-2
FITEM,2,9
FITEM,2,14
DK,P51X,,,,0,ALL,,,,,,
!
*
ACEL,0,0,9800,
/title, Gravityonly
/udoc,1,date,0
!
*
FINISH
!
*
/SOL
!
*
ANTYPE,0
NCNV,2,0,0,0,0
PSTRES,1
!
*
*if,case_i,eq,2,then !
自重加压力载荷
ffz=-1000*f_scale2
FLST,2,4,3,ORDE,4
FITEM,2,5
FITEM,2,-6
FITEM,2,11
FITEM,2,16
FK,P51X,FZ,ffz
/title, GravityaddFz=%ffz%
*endif
!
*
*if,case_i,eq,3,then !
自重加弯矩载荷
ffz=1000*f_scale3
FLST,2,2,3,ORDE,2
FITEM,2,6
FITEM,2,11
FK,P51X,FZ,ffz
!
*
FLST,2,2,3,ORDE,2
FITEM,2,5
FITEM,2,16
FK,P51X,FZ,-ffz
!
*
/title, Gravityadd|Fz|=%ffz%
*endif
!
*
!
*
SAVE,
!
/STATUS,SOLU
SOLVE
!
*
FINISH
/POST1
!
*
/auto,1
/USER, 1
/VIEW, 1, 0.534476049069 , 0.728260463362 ,-0.428919631721
/ANG, 1, -63.9329614215
/FOC, 1, -76.4296413144 , -127.691104361 , 2187.76597356
!
*
/udoc,1,date,0
/DSCALE,ALL,20
/EFACET,1
PLNSOL,U,Z,0,1.0
/wait,3 !
显示结果后,等待3秒,再继续后续的命令流
!
*
FINISH
save
!
*
/SOL
!
*
ANTYPE,1
!
*
BUCOPT,LANB,3,0,0,
!
*
!
/STATUS,SOLU
SOLVE
!
*
FINISH
!
*
/POST1
!
SET,LIST
SET,FIRST
!
*
/VIEW, 1, 0.827064809635 ,-0.475814768901 , 0.299272628817
/ANG, 1, -94.8821609302
!
*
/DSCALE,ALL,1500
/EFACET,1
PLNSOL,U,SUM,0,1.0
/wait,3
!
*
save
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 恒载加活载 特征值 屈曲 分析
![提示](https://static.bdocx.com/images/bang_tan.gif)