基于abaqus的工字钢柱的稳定性非线性分析.docx
- 文档编号:547826
- 上传时间:2022-10-11
- 格式:DOCX
- 页数:14
- 大小:402.81KB
基于abaqus的工字钢柱的稳定性非线性分析.docx
《基于abaqus的工字钢柱的稳定性非线性分析.docx》由会员分享,可在线阅读,更多相关《基于abaqus的工字钢柱的稳定性非线性分析.docx(14页珍藏版)》请在冰豆网上搜索。
基于abaqus的工字钢柱的稳定性非线性分析
【摘要】2..
【关键词】2..
」•计算方法和模型简介及分析过程3
1.1计算方法3.
1.2模型简介3.
1.3分析过程5.
1.3.1特征值屈曲分析5.
1.3.2非线性屈曲分析10
二计算结果分析1.3
2.1荷载-位移曲线1.3
2.2稳定系数比较1.4
【参考文献】16
【摘要】本文是基于Abaqus采用梁单元,对热轧工字钢在轴心受压情况下发生弯曲失稳的非线性屈曲分析。
通过考虑材料非线性、几何非线性并引入初弯曲这一几何初始缺陷,得出构件发生弯曲失稳的极限荷载。
通过保持柱子的截面不变,只改变长度,进而实现得到不同长细比的弯曲失稳的临界荷载,并得到相应的荷载位移曲线以及稳定系数,最后与规范给出的稳定系数做比较验证分析的正确性。
【关键词】特征值屈曲分析弯曲失稳非线性屈曲分析稳定系数
一.计算方法和模型简介及分析过程
1.1计算方法
利用abaqus进行屈曲分析,一般分为两步:
首先是特征值屈曲分析,此分析为线性屈曲分析,是在小变形的情况进行的,是用于预测一个理想弹性结构的理论屈曲荷载,
并在inp文件里做一定的修改,此修改是在下一步后屈曲分析所需要的初始缺陷的节点输出为.fil文件,作为下一步后屈曲分析所需要的初始缺陷的数据。
其次,就是非线性屈曲分析,
此步一般定义为非线性,原因在于是在大变形情况进行的,一般采用位移控制加修正的弧长
法(riks法),可以定义材料非线性,几何非线性,以及加上初始缺陷。
此步分析,为了得到极限值,需要得出荷载位移曲线的下降段,除了采用位移控制以及弧长法设定外,需在所
得到的inp文件中,嵌入第一步分析中的.fil节点数据。
1.2模型简介
本文利用abaqus有限元分析软件时,采用了梁单元,所选取的是热轧工字型截面,
截面参数见图一(单位m):
EditProMe
Name:
Profile-1
Shape:
I
图一热轧工字钢截面参数
又,其中的材料本构关系采用双线性塑性材料模型,所选材料为Q345钢材,故弹性模量为:
2.1e11Pa,屈服强度取:
3.45e8Pa,泊松比取0.3。
本文一共做了15个模型的分析,其中每一个模型的截面参数保持不变,通过改变杆件的计
算长度进而实现不同的长细比分析,各模型的截面惯性矩及长细比等详见表1.2模型计算参
数详表。
表1.2模型计算参数详表
编号
计算长度1
(m)
截面面积A
/2(m)
I4
1y(m)
・2
iy(m)
1
2.0
0.00568
0.00001334
0.04846643
41.27
2
2.5
0.00568
0.00001334
0.04846643
51.58
3
3.0
0.00568
0.00001334
0.04846643
61.90
4
3.5
0.00568
0.00001334
0.04846643
72.21
5
4.0
0.00568
0.00001334
0.04846643
82.53
6
4.5
0.00568
0.00001334
0.04846643
92.85
7
5.0
0.00568
0.00001334
0.04846643
103.16
8
5.5
0.00568
0.00001334
0.04846643
113.48
9
6.0
0.00568
0.00001334
0.04846643
123.8
10
6.5
0.00568
0.00001334
0.04846643
134.11
11
7.0
0.00568
0.00001334
0.04846643
144.43
12
7.5
0.00568
0.00001334
0.04846643
154.75
13
8.0
0.00568
0.00001334
0.04846643
165.06
14
8.5
0.00568
0.00001334
0.04846643
175.38
15
9.0
0.00568
0.00001334
0.04846643
185.70
注:
由于弯曲失稳是沿着弱轴方向失稳,故在此只列岀弱轴方向的相关计算参数,并不在另外列岀强轴方向的相关计算参数。
1.3分析过程
以长细比为51.58的为例;
1.3.1特征值屈曲分析:
第一步:
在Part模块定义梁单元
第二步:
在Part模块定义构件长度(2.5m)
分别输入(0,1.25),(0,-1.25)
第三步:
在Property模块定义截面属性
Name:
steel
Description:
Type:
Isotropic||
EZ1Usetemperature-dependentdata
NumDeroffieldvariables:
O(=^
Modulitimescale(forviscoelasticity):
Long-term□□Nocompression
IiNotension
(—Deto
YounQ*smodulus
Poisson*sRotio
丄210000000000
0.3
EditProfile
Name:
Profile-1
Shape:
I
I-
J
1
t
3
•
i
T
V—
b
Cancel
第四步:
在Property模块定义截面方向
点击菜单栏的"引叩并选择BeamSectionOrientation后输入(1.0,0.0,0.0)
第五步:
在Assembly模块组装模型
lInstanceType
■.•-dependent(meshonpart)i
OIndependent(meshoninstance)
Note;To匚hangeaDependentinstancedmeshjyoumustedititspart'smesh.
□Auto-ofPsetfromotherinjtanres
第六步:
在Step模块定义buckle分析步
第七步:
在Load模块定义模型的边界条件
定义轴向力
定义约束条件
第八步:
在Mesh模块定义网格划分
第九步:
在任意模块修改INP文件,修改语句如下*nodefile
第十步:
在Job模块提交并计算完成buckle分析
1.3.2非线性屈曲分析:
第一步:
拷贝前一步即特征值屈曲分析所用的模型
第二步:
在Step模块定义Riks分析步
CreateStep
Proceduretyipe:
Generala
Continue.-a
第三步:
引入初始缺陷
在第二个模型的任意模块,点击菜单栏中的model,再点击其下拉栏中的EditKeywords
修改INP文件,修改内容如下:
*innperfectionJfile=buckleb3’step=1
tle-2
2Je-3
其含义即为引入buckle分析步中的一阶模态的1%,二阶模态的1%。
作为非线性屈曲分析中
模型的初弯曲。
第五步:
在Visualization模块中进行后处理
画出荷载-位移曲线如下
提取数据至Excel进行进一步处理。
.计算结果分析
编号
1
2
3
4
5
6
7
8
9
10
长细比入
41.27
51.58
61.90
72.21
82.53
92.85
103.16
113.48
123.8
134.11
编号
11
12
13
14
15
16
17
18
19
20
长细比入
144.43
154.75
165.06
175.38
185.70
■-
■-
■-
■-
■-
2.1荷载-位移曲线
2.00E+C6
1.80E-^S
1.60E+06
1.40E-HI6
1.20E-+06
1.00E+O6
3.00E-+05
6.00E+C5
4.OOE-+C5
2.00E+05
0.00E-K0
O.OOEWLOOE-Ol2.00E-C13,OOE-Ol4.00E-015.00EF1G.00EP1
卄“虬裁*-A=51,53
A=61,90
A=72.21—A=02.53
Displacement(n)
O.OOE+OOLOOE-Ol2.OOE-Ol3*OOE-Ol4.00E-016.OOE-Ol6.OOE-OlLOOE-Ol8*OOE-Ol
DisplaGement(m)
T-A农.S5-A-103.16
k=113.40
A=123.80
十A=134.11
Displacement(m)
从以上三个荷载位移曲线可知,长细比越小,弯曲失稳的极限荷载越大;达到了弯曲失稳的
极限荷载以后,荷载位移曲线均呈现下降段。
2.2稳定系数比较
《钢结构设计规范》GB50017-2003表5.1.2-1知本文所用模型截面分类属于b类;
《钢结构设计规范》GB50017-2003中附录C轴心受压构件的稳定系数的计算公式:
当入二\fy<0.215时,「=1-內才;
ne
当40-215时,二右卜2*3入•兀!
-加2*3入•/2-4才(2-2)
其中,上述计算公式2-2中的系数坷、aa按表2.2查取:
表2.2公式2-2的系数內、a、a
截面类别
a
a
a
a类
0.41
0.986
0.152
b类
0.65
0.965
0.300
c类
4<1.05
0.73
0.906
0.595
4>1.05
1.216
0.302
d类
4<1.05
1.35
0.868
0.915
4>1.05
1.375
0.432
由于规范中的系数a、a、a是考虑了残余应力、初始弯曲的,另本文的模型并不考虑残
余应力的影响,结合上述的稳定系数计算公式2-2,为了具有一定的可比性,在选取系数a
时,取为1.0,其余系数依旧按照表2.2选取。
而本文中的模型的稳定系数按以下公式计算:
抗压强度标准值,取为3.45e8N/mm2。
由本文的数据算得的稳定系数与规范中的稳定系数的详细比较见表2.2.1;
表2.2.1稳定系数的比较
编号
长细比
失稳极限荷
载(N)
①Nu
Afy
规范中的稳定系数甲
误差
1
41.27
1751250
0.894
0.825
8.30%
2
51.58
1441980
0.736
0.765
3.76%
3
61.9
1319400
0.673
0.694
2.99%
4
72.21
1189110
0.607
0.617
1.62%
5
82.53
1064460
0.543
0.539
0.78%
6
92.85
923360
0.471
0.467
0.99%
7
103
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 abaqus 工字钢 稳定性 非线性 分析