屈曲分析全过程
屈曲分析全过程
屈曲分析的过程说明:屈曲分析是一种用于确定结构开始变得不稳定时的临介荷载和屈曲结构发生屈曲响应时的模态形状的技术。
ANS 丫醍供两种结构屈曲荷载和屈曲模态分析方法:非线性屈曲分析和特征值屈曲分析。
非线性屈曲分析是在大变形效应开关打开的情况下的一种非线性静力学分析,该分析过程一直进行到结构的极限荷载或最大荷载。
非线性屈曲分析的方法是,逐步地施加一个恒定的荷载增量,直到解开始发散为止。
尤其重要的是,要一个足够小的荷载增量,来使荷载达到预期的临界屈曲荷载。
若荷载增量太大,贝屈曲分析所得到的屈曲荷载就可能不准确,在这种情况下打开自动时间步长功能,有助于避免这类问题,打开自动时间步长功能,ANS YSS序将自动寻找屈曲荷载。
特征值屈曲分析步骤为:1.建模2.获得静力解:与一般静力学分析过程一致,但必须激活预应力影响,通常只施加一个单位荷载就行了3.获得特征屈曲解:A.进入求解B.定义分析类型C.定义分析选项D.定义荷载步选项E.求解 4.扩展解之后就可以察看结果了 示例1:! ansys 7.0有限元分析实用教程 ! 3.命令流求解 ! ANSYS 命令流:! Eigenvalue BucklingK,1,0,0!创建梁实体模型K, 2,0,100L, 1,2 !创建直线 单元边长为1mmFINISH !这两行命令清除当前数据/CLEAR/TITLE,Eige nvalue Buckli ng An alysis/PREP7 !进入前处理器 ET,1,BEAM3 !选择单元 R,1,100,833.333,10 !定义实常数 MP,EX,1,200000 !弹性模量 MP,PRXY,1,0.3!泊松比ESIZE,10LMESH,ALL,ALL !划分网格FINISH !退出前处理!屈曲特征值部分/SOLU !进入求解ANTYPE,STATIC !在进行屈服分析之前,ANSYS需要从静态分析提取数据PSTRES,ON !屈服分析中采用预应力DK,1,ALL !定义约束FK,2,F Y,-1 !顶部施加载荷SOLVE !求解FINISH !退出求解/SOLU !重新进入求解模型进行屈服分析ANTYPE,BUCKLE !屈服分析类型BUCOPT, LANB,1 ! 1阶模态,子空间法I SOLVE !求解FINISH !退出求解/SOLU !重新进入求解展开模态EXPASS,ON !模态展开打开MXPAND,1 !定义需要展开的阶数SOLVE !求解FINISH !退出求解/POST1 !进入通用后处理SET,LIST !列出特征值求解结果SET,LAST !入感兴趣阶数模态结果PLDISP !显示变形后图形! Non Li near Buckl ing !非线性分析部分FINISH !这两行命令清除当前数据/CLEAR/TITLE, Non lin ear Buckli ng An alysis/PREP7 !进入前处理ET,1,BEAM3 !选择单元MP,EX,1,200000 !弹性模量MP,PRXY,1,0.3 !泊松比R,1,100,833.333,10 !定义实常数K,1,0,0,0 !底端节点K,2,0,100,0 !顶点L,1,2 !连成线ESIZE,1 !网格尺寸参数设定LMESH,ALL !划分网格FINISH !退出前处理/SOLU !进入求解ANTYPE,STATIC !静态分析类型(非屈服分析)NLGEOM,ON !打开非线性大变形设置OUTRES,ALL,ALL !选择输出数据NSUBST,20 ! 5个子步加载NEQIT,1000 ! 20步迭代AUTOTS,ON !自动时间步长LNSRCH,ON !激活线搜索选项/ESHAPE,1 !显示二维状态下变形图DK,1,ALL,0 !约束底部节点FK,2,F Y,-50000 !顶部载何稍微比特征值分析结果大FK,2,FX,-250 !施加水平扰动载荷SOLVE !求解FINISH !退出求解/POST26 !进入时间-历程后处理器RFORCE,2,1,F, Y ! 2#变量表示力/AXLAB,Y,DEFLECTION !修改 y 轴标签/REPLOT示例2:!悬臂梁受端部轴向压力作用的屈曲分析!先进行静力分析,在进行特征值屈曲分析,最后进行非线性分析!静力分析 fini/cle/filn ame,beam-flexure/tittle,beam-flexure /prep7 *set,f1,-1e6向压力荷载参数et,1,beam189NSOL,3,2,U, Y 3#变量表示y 方向位移XVAR,2将x 轴显示2 #变量PLVAR,3y 轴显示3 #变量数据/AXLAB,X , LOAD修改x 轴标签 重新显示图形设置轴mp,de ns,1,7.85e3 设置材料参数mp,ex,1,2.06e11mp,nu xy,1,0.2sectype,1,beam,l,,2 ! 设置截面参数secoffset,ce ntsecdata,0.15,0.15,0.25,0.015,0.015,0.015,0,0,0,0k,1,0k,2,2.5,0k,3,1.25,1lstr,1,2latt,1,,1,,3,,1lesize,1,,,10lmesh,1/view,1,1,1,1/eshape,1.0dk,1,,,,0,all,fk,2,fx,f1 !施加关键点压力finish/soluan type,0eqslv,spar !求解器设置稀疏矩阵直接法pstres,on !打开预应力开关solvefin ish!特征值屈曲分析/soluan type,1bucopt,lanb,6,0 !取前六阶模态分析mxpa nd,6,0,0,1,0.001solvefin ish/post1*do,i,1,6set,ipldisp,i*enddo*get,freq1,mode,1,freq fin ish!非线性屈曲分析/config,n res,200两百步的结果/prep7tb,biso,1,1,2料非线性tbtemp,0tbdata,,2.0e8,0upgeom,0.01,1,1,'beam-flexure','rst'!对有限元模型进行0.01倍的修改save,beam-flexure,dbfin ish只记录定义材阶模态的位移结果resu,beam-flexure,db/solu an type,0nlgeom,1 !打开大变形outres,all,allarclen,1,0 !弧长法设置arctrm,l !弧长法终止准则达到第一个峰值时终止计算nsubst,200,,,1fk,2,fx,f1*freq1!fk,2,fx,f1*freq1*1.2 !将轴向压力值放大,放大系数为第一阶模态的主频solve finish。
混凝土梁的屈曲分析方法
混凝土梁的屈曲分析方法一、引言混凝土梁是建筑结构中常见的一种构件,主要承受弯曲力和剪力。
在梁的设计和施工过程中,必须对其进行力学分析和计算,以确保其稳定性和安全性。
其中,屈曲是混凝土梁设计中的一个重要问题,如何进行屈曲分析是建筑工程师需要掌握的重要技能。
本文将介绍混凝土梁的屈曲分析方法。
二、混凝土梁的屈曲分析方法1.梁的截面分析梁的截面形状和尺寸对其承受弯曲力的能力有很大影响。
因此,在进行屈曲分析之前,必须对梁的截面进行分析,以确定其承受弯曲力的能力。
梁的截面分析包括以下几个步骤:(1)确定梁的几何尺寸和截面形状;(2)计算梁的惯性矩和截面模量;(3)计算梁的截面抗弯能力。
2.梁的荷载分析梁的荷载分析包括确定梁的荷载大小和作用位置,以及计算梁的剪力和弯矩。
在进行屈曲分析之前,必须对梁的荷载进行分析,以确定其在受载状态下的荷载大小和作用位置。
梁的荷载分析包括以下几个步骤:(1)确定梁的荷载大小和作用位置;(2)计算梁的剪力和弯矩;(3)确定梁的受力状态。
3.梁的稳定性分析在确定梁的截面形状、尺寸和荷载大小和作用位置后,可以进行梁的稳定性分析。
梁的稳定性分析主要是分析梁的屈曲状态和屈曲载荷。
梁的稳定性分析包括以下几个步骤:(1)计算梁的屈曲载荷;(2)确定梁的屈曲模态;(3)确定梁的屈曲状态;(4)判断梁的稳定性。
4.屈曲载荷的计算梁的屈曲载荷是指当梁受到一定弯矩作用时,梁开始发生屈曲时所承受的最大荷载。
计算梁的屈曲载荷需要考虑以下因素:(1)梁的截面形状和尺寸;(2)梁的材料性质;(3)梁的荷载大小和作用位置。
5.屈曲模态的确定梁的屈曲模态是指梁在屈曲时的挠曲形态。
梁的屈曲模态可以通过数值计算或实验测量来确定。
在进行屈曲计算时,必须考虑梁的屈曲模态对计算结果的影响。
6.屈曲状态的确定梁的屈曲状态是指梁在受到弯曲力作用下的变形状态。
梁的屈曲状态可以通过分析梁的截面形状和尺寸、材料性质和荷载大小和作用位置来确定。
MIDAS_FEA_屈曲分析教程
Step
10 几何 > 曲面> 面并集...
操作步骤
1、在选择曲面里中选 择窗口上所有的线,
2、点击确认
3、将所有的面缝合成 了一个面组。
Step
11 网格 > 自动网格划分> 自动网格面...
操作步骤
1、选择全面的面(共2 6个面)。
2、网格划分方法选择 栅格网格法,确认类型 是否选择为四边形+三 角形
3、P3,P4也按同样的步
骤进行修剪,最终将四
P1
根管生成一个面组。
注意:每次选择对象形 状前,都需要重复步骤 1,选择面组的方式。
P2
P3 P4
Step
07 几何 > 工作平面> 移动
操作步骤
1、将工作平面移动到 需要建立辅助面的地方。
2、在Dz中输入100
3、点击确认。
Step
08 几何 > 曲线 > 在工作平面上创建> 2D圆...
3、确认方向是否为法 向,是否勾选“均匀”
4、在P(P1)中输入20
5、点击确认。
Step
15 分析 > 分析工况…
操作步骤
1、添加分析工况 输入分析工况名称p1
2、选择分析类型为 “线
性屈曲”
3、分析控制中输入频 率数为5.
4、将定义好的边界条 件、荷载拖放到对应的 应用的菜单下。
5、点击确认。
Step
13 分析> 边界条件> 约束
操作步骤
1、选择前视图模式。
2、定义边界组名称: p1
3、选择主管底面所有 节点,共(40个)。
4、选择固定的约束方 式,点击确认即可。
Step
ansys 屈曲分析详细过程
非线性屈曲分析的流程图如下:
图 2 非线性屈曲分析流程图
2
3、非线性方程组求解方法 (1)增量法 增量法的实质是用分段线性的折线去代替非线性曲线。增量法求解时将荷载
分成许多级荷载增量,每次施加一个荷载增量。在一个荷载增量中假定刚度矩阵 保持不变,在不同的荷载增量中,刚度矩阵可以有不同的数值,并与应力应变关 系相对应。
图 4 荷载及拱的几何尺寸示意
几何尺寸: f/L=1/5, H/L=1/30,B/图L=图1图/15
L=89m,f=17.8m,H=2.97m,B=5.93m,R=64.5m
材料性能:拱圈采用 C40 混凝土,弹性模量为 32500MPa,泊松比为 0.3。
在 ANSYS 中进行建模分析,初始荷载为 q=100000kN/m 其过程如下:
0.5
K2
58.9
90.4
93.4
86.7
64.0
根据表 1 查得:
K2 90.4
故其理论弹性屈曲荷载为:
qcr
K2
EI x l3
3.25104 1 5933.33 2966.673
90.4
12 890003
5.38107
N
m
2、拱的弹性屈曲与非线性屈曲 对于一般的特征值屈曲分析,主要是在平衡状态,考虑到轴向力或者中面内
力对弯曲变形的影响,由最小势能原理,结构弹性屈曲分析归结为求解特征值问 题:
KD KG 0
通过特征值分析求得的解有特征值和特征向量,特征值就是临界荷载系数, 特征向量是临界荷载系数对应的屈曲模态。
特征值屈曲分析的流程图如下:
1
图 1 弹性屈曲分析流程图
非线性屈曲分析是考虑结构平衡受扰动(初始缺陷、荷载扰动)的非线性静 力分析,该分析是一直加载到结构极限承载状态的全过程分析,分析中可以综合 考虑材料塑性和几何非线性。结构非线性屈曲分析归结为求解矩阵方程:
屈曲分析流程
屈曲分析流程Bending analysis is a crucial process in the engineering and design of various structures and components. It is essential for determining the structural integrity, safety, and performance of materials under different loading conditions. 屈曲分析是工程和设计中的一个至关重要的过程,对于确定材料在不同载荷条件下的结构完整性、安全性和性能至关重要。
The process of bending analysis involves the evaluation of the stress and strain distribution in the material, as well as the determination of critical points where failure may occur. 屈曲分析的过程涉及材料中应力和应变分布的评估,以及确定可能发生失效的关键点。
By understanding the behavior of materials under bending, engineers are able to optimize the design of various structures, ensuring that they can withstand the required loads and perform effectively in their intended applications. 通过了解材料在弯曲下的行为,工程师能够优化各种结构的设计,确保它们能够承受所需的载荷并在其预期的应用中有效地发挥作用。
One of the primary aspects of bending analysis is the determination of the maximum bending moment and the corresponding stress distribution along the length of the material. 屈曲分析的一个主要方面是确定最大弯曲力矩以及材料长度沿线的相应应力分布。
abaqus屈曲分析报告实例
整个计算过程包括2个分析步,第1步做屈曲分析,第2步做极限强度分析。
第1步:屈曲分析载荷步定义如下:Step 1-InitialStep 2- Buckle并在Model-Edit Keywords的图中位置加入下面的文字,输出屈曲模态*nodefile, global=yesU,Create job 名称为“Buckling”点击continue,完成第1步的计算。
第2步:极限强度分析将“buckle”分析步替换为“riks”分析步在Basic选项卡中,Nlgeom:选择打开在Instrumentation选项卡中,定义如下参数,然后点击OK定义一个新计算工作,输入名称,点击continue在Parallelization选项卡,选择2个CPU,如下所示,点击OK。
在此编辑Model-edit keywords,删除“第1步”加入的文字“*nodefile, global=yesU,”,并在下图位置加入下段文字:*imperfection, file=buckling, step=11,点击OK,再保存文件。
最后提交计算。
提取计算结果进入visualization Module 点击 Create XY data选择 ODB filed output,点击continuePosition选择 Unique Nodal, CF:point loads选择 CF2,再点击elements/nodes选项卡,选择跨中载荷加载点,最后点击save。
重复上一步操作,Position选择 Unique Nodal, U:spatial displacement 选择 U3,再点击elements/nodes选项卡,选择板格中心点,最后点击save。
点击Create XY data, 选择operate on XY data,点击continue选择Combine(X,X)命令,横坐标选择保存的displacement曲线,纵坐标选择保存的Point load曲线,点击最后一行Create XY Data与Save as。
屈曲分析
弧型闸门支臂的屈曲分析一、问题和条件:钢闸门断面为工字钢,具体尺寸如下,试作屈曲分析:弹性模量E=200GPa泊松比μ=0.3集中力 F = 1N(单位力)二、进入ANSYS界面单击开始→程序→ansys8.0→configure ANSYS product然后在File Management中定义Working Directory(工作路径)如:class定义Job Name(工作文件名)。
如:zhamenququ三、定义单元及材料1 新建单元类型运行主菜单Preprocessor → Element Type → Add/Edit/Delete (新建/编辑/删除单元类型)命令,接着在对话框中单击“Add”按钮新建单元类型。
2 定义单元类型因为beam188单元适合于分析柔性梁结构,支持弹性,蠕变和塑性模型的计算。
先选择单元为beam,接着选择3D 2 node 188,然后单击“OK”按钮确定,完成单元类型的选择。
3 关闭单元类型对话框回到单元类型对话框,已经新建了beam单元,单击对话框中“Close”按钮关闭对话框。
4 定义截面形状及参数运行主菜单Preprocessor → Sections → Beam→Common sections命令,接着在对话框中sub-type中选择“工”字钢断面,在下面的尺寸中分别填入:W1=0.2, W2=0.2, W3=0.4,t1=0.01, t2=0.01, t3=0.01,点击OK结束。
5 设置材料属性运行主菜单Preprocessor → Material props→ Material Models (材料属性)命令。
选择材料属性命令后,系统会显示材料属性设置对话框。
6 设置杨氏弹性模量与泊松比在材料属性设置对话框右侧依次选择 Structural → Linear → Elastic → Isotropic。
完成选择命令后,接着在对话框中EX杨氏弹性模量输入2e11,PPXY泊松比输入0.3。
第17章 屈曲分析
第17章屈曲分析第1节基本知识屈曲分析是确定结构开始变得不稳定时的临界载荷和屈曲模态形状的技术。
在ANSYS 程序中提供了两种屈曲分析技术:特征值屈曲分析(线性屈曲分析)和非线性屈曲分析。
经典的屈曲分析是采用特征值屈曲分析法,所谓特征值失稳计算就是用结构的材料刚度矩阵减去荷载作用下结构的几何刚度乘以一个系数,当总刚度矩阵奇异时的就是失稳特征值。
它适用于对一个理想弹性结构的理想屈曲强度(歧点)进行预测,主要是使用特征值公式计算造成结构负刚度的应力刚度阵的比例因子。
结构在达到屈曲载荷之前其位移——变形曲线表现出线性关系,达到屈曲以后其位移——变形曲线表现出非线性关系,此种方法满足于经典教本理论。
然而,在实际的工程结构中会有一定的初始缺陷,而且在使用过程中会出现材料非线性以及大变形等非线性因素,使结构并不全是在其理想弹性屈曲强度处发生屈曲。
故特征值屈曲经常产生非保守结果即临界载荷值较大,不适用于工程结构屈曲分析,由此应运而生的是非线性屈曲分析法。
该方法是包括材料非线性、大变形等非线性因素的静力分析法,计算过程可以一直进行到结构的限制载荷或最大载荷。
因此,在实际结构的设计和估计中宜采用后一种方法。
一、特征值屈曲分析特征值屈曲分析一般由以下五个步骤:1)建立模型;2)获得静力解;3)获得特征值屈曲解;4)拓展结果;5)后处理,观察结果及输出。
在特征值屈曲分析中应注意以下几个问题:1)在前处理器PREP7中定义单元类型、实常数、材料性质和创建几何模型和有限元模型。
其中只允许线性行为,若定义了非线性单元,按线性对待。
材料的性质可以是线性、各向同性或各向异性、恒值或与温度相关的。
2)求静力解时必须激活预应力(PSTRES)选项,特征值屈曲分析需要计算应力刚度矩阵。
屈曲分析因为计算出的特征值表示屈曲荷载系数,所以一般施加单位载荷,但是ANSYS 程序的特征值限值是1 000 000,所以如果求解的特征值超过此限值的话,我们应当施加一个比较大的荷载,而不应该还是用单位力施加!计算结果为:临界力=施加力×求出来的屈曲系数,所以一般的用单位力施加的临界力就等于求出的屈曲系数,因为施加力为1。
第17章 屈曲分析
第17章屈曲分析第17章屈曲分析第1节基本知识屈曲分析是确定结构开始变得不稳定时的临界载荷和屈曲模态形状的技术。
在ANSYS程序中提供了两种屈曲分析技术:特征值屈曲分析(线性屈曲分析)和非线性屈曲分析。
经典的屈曲分析是采用特征值屈曲分析法,所谓特征值失稳计算就是用结构的材料刚度矩阵减去荷载作用下结构的几何刚度乘以一个系数,当总刚度矩阵奇异时的就是失稳特征值。
它适用于对一个理想弹性结构的理想屈曲强度(歧点)进行预测,主要是使用特征值公式计算造成结构负刚度的应力刚度阵的比例因子。
结构在达到屈曲载荷之前其位移――变形曲线表现出线性关系,达到屈曲以后其位移――变形曲线表现出非线性关系,此种方法满足于经典教本理论。
然而,在实际的工程结构中会有一定的初始缺陷,而且在使用过程中会出现材料非线性以及大变形等非线性因素,使结构并不全是在其理想弹性屈曲强度处发生屈曲。
故特征值屈曲经常产生非保守结果即临界载荷值较大,不适用于工程结构屈曲分析,由此应运而生的是非线性屈曲分析法。
该方法是包括材料非线性、大变形等非线性因素的静力分析法,计算过程可以一直进行到结构的限制载荷或最大载荷。
因此,在实际结构的设计和估计中宜采用后一种方法。
一、特征值屈曲分析特征值屈曲分析一般由以下五个步骤:1)建立模型;2)获得静力解;3)获得特征值屈曲解;4)拓展结果;5)后处理,观察结果及输出。
在特征值屈曲分析中应注意以下几个问题:1)在前处理器PREP7中定义单元类型、实常数、材料性质和创建几何模型和有限元模型。
其中只允许线性行为,若定义了非线性单元,按线性对待。
材料的性质可以是线性、各向同性或各向异性、恒值或与温度相关的。
2)求静力解时必须激活预应力(PSTRES)选项,特征值屈曲分析需要计算应力刚度矩阵。
屈曲分析因为计算出的特征值表示屈曲荷载系数,所以一般施加单位载荷,但是ANSYS程序的特征值限值是1 000 000,所以如果求解的特征值超过此限值的话,我们应当施加一个比较大的荷载,而不应该还是用单位力施加!计算结果为:临界力=施加力×求出来的屈曲系数,所以一般的用单位力施加的临界力就等于求出的屈曲系数,因为施加力为1。
MIDAS_FEA_屈曲分析教程_荣萌
4、在P(P1)中输入20
5、点击确认。
Step
15 分析 > 分析工况…
操作步骤
1、添加分析工况 输入分析工况名称p1
2、选择分析类型为 “线
性屈曲”
3、分析控制中输入频 率数为5.
4、将定义好的边界条 件、荷载拖放到对应的 应用的菜单下。
5、点击确认。
Step
13 分析> 边界条件> 约束
操作步骤
1、选择前视图模式。
2、定义边界组名称: p1
3、选择主管底面所有 节点,共(40个)。
4、选择固定的约束方 式,点击确认即可。
Step
14 分析> 荷载> 压力荷载
操作步骤
1、输入荷载组名称, 类型选择“线压力”。
2、选择主管上部面的 所有单元线(共40个)
Step
00
概要
▪屈曲分析 ▪模型
- 单位 : N mm - 各向同性弹性材料 - 实体单元 ▪荷载和边界条件 - 压力 - 约束 输出结果 - 变形 - 应力
管桁架屈曲分析
Step
01 分析 > 分析控制 – 控制 表单
操作步骤
1、分析类型 : [3D] 1
2、点击[确认]键
3、单位 : [N,mm]
柱。
Step
03 几何 >删除 > 删除子面…
操作步骤
1、选择圆柱上部面和 底面(共四个子面)
2、点击预览按钮确认 形状是否适当,然后点 击确认。
Step
04 几何 > 转换 > 旋转...
操作步骤
1、选择细圆管,用移 动的方式旋转90度,如 图2所示。
ANSYSWorkbench基础教程与工程分析详解第六章屈曲分析
工程中很多结构需要进行结构稳定性计算,如细长杆、压缩部件、真空容器等,这些构件在不稳定(屈曲)开始时,结构本质上没有变化的载荷作用下(超过一个很小的动荡),在x 方向上的微小位移会使得结构有一个很大的改变,这类问题除了要考虑强度之外,还要分析其屈曲稳定性的问题。
本章所要学习的内容包括: ¾ 了解线性屈曲分析基础¾ 掌握ANSYS Workbench 屈曲分析的操作流程 ¾ 了解线性屈曲分析的应用场合 ¾ 理解屈曲分析的结果6.1 线性屈曲分析基础特征值或线性屈曲分析预测的是理想线弹性结构的理论屈曲强度(分歧点);而非理想和非线性行为阻止许多真实的结构达到它们的理论上的弹性屈曲强度。
线弹性通常产生非保守的结果,但也是有优点的。
(1)它比非线性屈曲计算省时间,并且应当做第一步计算来评估临界载荷(屈曲开始的载荷)。
(2)线性屈曲分析可以用来作为决定产生什么样的屈曲模型形状的设计工具,为设计做指导。
线性屈曲的分析方程为:{}([][])0i i K S λΨ+=式中各个符号的含义如下。
S 表示应力刚度矩阵; i λ表示屈曲载荷乘子;i Ψ表示屈曲模态。
实际上,线性屈曲方程和自由振动方程很相似,两者都是利用相似的矩阵方法来求解特征值问题的。
线性屈曲的分析步骤与之前的静力学分析非常相似,过程如下。
(7)求解计算并保存。
ANSYS Workbench1 4.5屈曲模态分析步骤与其他有限元分析步骤大同小异,软件支持模态分析中存在接触对,但因为屈曲分析是线性分析,所以接触行为不同于非线性接触行为,接触设置的线性屈曲分析设置如表6-1所示。
表6-1 存在接触设置的线性屈曲分析设置Linear Buckling Analysis(线性屈曲分析)Contact Type (接触类型) Initially Touching (初始接触) Inside Pinball Region (Pinball 区域内) Outside Pinball Region (Pinball 区域外) Bonded (绑定) Bonded (绑定) Bonded (绑定) Free (自由) No Separation (不分离) No Separation (不分离) No Separation (不分离) Free (自由) Rough (粗糙) Bonded (绑定) Free (自由) Free (自由) Frictionless (光滑)No Separation (不分离)Free (自由)Free (自由)6.2 案例图解6.2.1 斜撑杆受压屈曲分析分析起落架中承受轴向压力的斜撑杆,杆为空心圆管,外径为52mm ,内径为44mm ,L =950mm 。
恒载加活载的特征值屈曲分析
工程上有时会遇到既有恒载 (固定不变的载荷) 又有活载 (按固定比例变化的载荷) 的特征值屈曲分析问题,比方高层建筑在自重和风载的共同作用下会发生侧弯,其受压严重的区域有发生屈曲的可能。
这里的自重就是恒载,风载就是活载。
如果做特征值屈曲分析 (以下简称屈曲分析),由于它假定了所有载荷都是按一样的比例变化的,而既有恒载又有活载的构造,恒载是不变的,这种问题应该如何进展屈曲分析呢?对这类问题,其屈曲分析的过程是一个迭代的过程,根本步骤如下:1. 如果单独恒载有可能造成构造的屈曲,应该先单独施加恒载并进展屈曲分析。
只有在单独恒载作用不会使构造发生屈曲的前提下,才需要进展恒载+活载的屈曲分析。
2. 首先施加给定的所有载荷,进展屈曲分析,求出屈曲载荷因子 (最小的正特征值);3. 恒载不变,只将活载局部除以所得的屈曲载荷因子,然后再一次进展屈曲分析,得到新的载荷因子。
4. 重复第 3 步,直到屈曲载荷因子根本等于 1,这时的活载就是在存在恒载的情况下,构造发生屈曲时的活载的大小。
下面是一个简单的例子 (构造可能不太合理,只为说明问题):一个井架,由 4 根纵向长梁和假设干水平的或斜的短梁组成,形状如图 1所示。
纵梁的截面为 30*30 mm2,短梁的截面为 10*10 mm2。
材料性能为:弹性模量 E = 202100 Mpa泊松比μ= 0.3密度ρ=7.8e-9 * 2 〔见下面关于恒载的说明〕使用BEAM188 梁单元划分网格,单元长度取50 mm。
建模过程可见附录中命令流的几何建模局部,这里不再详述。
约束井架 4 根长梁底部的所有位移自由度。
考虑两种载荷:(1)恒载井架以及固定设备的自重。
建模时不考虑设备,通过加大各杆的材料密度将设备重量分配到井架的各杆上。
这里为了简单,只定义了一种材料,即各根杆是按照原有重量等比例的放大。
如果想要对不同的杆按不同比例放大,可以定义多种材料,分别给予不同的密度。
本模型中,将材料密度放大了一倍。
ansys 屈曲分析详细过程
0.5
K2
58.9
90.4
93.4
86.7
64.0
根据表 1 查得:
K2 90.4
故其理论弹性屈曲荷载为:
qcr
K2
EI x l3
3.25104 1 5933.33 2966.673
90.4
12 890003
5.38107
N
m
2、拱的弹性屈曲与非线性屈曲 对于一般的特征值屈曲分析,主要是在平衡状态,考虑到轴向力或者中面内
力对弯曲变形的影响,由最小势能原理,结构弹性屈曲分析归结为求解特征值问 题:
KD KG 0
通过特征值分析求得的解有特征值和特征向量,特征值就是临界荷载系数, 特征向量是临界荷载系数对应的屈曲模态。
特征值屈曲分析的流程图如下:
1
图 1 弹性屈曲分析流程图
非线性屈曲分析是考虑结构平衡受扰动(初始缺陷、荷载扰动)的非线性静 力分析,该分析是一直加载到结构极限承载状态的全过程分析,分析中可以综合 考虑材料塑性和几何非线性。结构非线性屈曲分析归结为求解矩阵方程:
32510593333296667129045381089000crei2拱的弹性屈曲与非线性屈曲对于一般的特征值屈曲分析主要是在平衡状态考虑到轴向力或者中面内力对弯曲变形的影响由最小势能原理结构弹性屈曲分析归结为求解特征值问通过特征值分析求得的解有特征值和特征向量特征值就是临界荷载系数特征向量是临界荷载系数对应的屈曲模态
QSV QH / (1 y2 ) 和QSP QH y / (1 y2 ) 其中,y=df (x) / dx,当单元足够小时,可采用y y / y
图 3 曲梁均布荷载等效与分解
3
5、问题的描述 矩形截面圆弧拱桥,在梁上受均布荷载,荷载及截面尺寸如图 1 所示,要求
第17章 屈曲分析
第17章屈曲分析第1节基本知识屈曲分析是确定结构开始变得不稳定时的临界载荷和屈曲模态形状的技术。
在ANSYS 程序中提供了两种屈曲分析技术:特征值屈曲分析(线性屈曲分析)和非线性屈曲分析。
经典的屈曲分析是采用特征值屈曲分析法,所谓特征值失稳计算就是用结构的材料刚度矩阵减去荷载作用下结构的几何刚度乘以一个系数,当总刚度矩阵奇异时的就是失稳特征值。
它适用于对一个理想弹性结构的理想屈曲强度(歧点)进行预测,主要是使用特征值公式计算造成结构负刚度的应力刚度阵的比例因子。
结构在达到屈曲载荷之前其位移——变形曲线表现出线性关系,达到屈曲以后其位移——变形曲线表现出非线性关系,此种方法满足于经典教本理论。
然而,在实际的工程结构中会有一定的初始缺陷,而且在使用过程中会出现材料非线性以及大变形等非线性因素,使结构并不全是在其理想弹性屈曲强度处发生屈曲。
故特征值屈曲经常产生非保守结果即临界载荷值较大,不适用于工程结构屈曲分析,由此应运而生的是非线性屈曲分析法。
该方法是包括材料非线性、大变形等非线性因素的静力分析法,计算过程可以一直进行到结构的限制载荷或最大载荷。
因此,在实际结构的设计和估计中宜采用后一种方法。
一、特征值屈曲分析特征值屈曲分析一般由以下五个步骤:1)建立模型;2)获得静力解;3)获得特征值屈曲解;4)拓展结果;5)后处理,观察结果及输出。
在特征值屈曲分析中应注意以下几个问题:1)在前处理器PREP7中定义单元类型、实常数、材料性质和创建几何模型和有限元模型。
其中只允许线性行为,若定义了非线性单元,按线性对待。
材料的性质可以是线性、各向同性或各向异性、恒值或与温度相关的。
2)求静力解时必须激活预应力(PSTRES)选项,特征值屈曲分析需要计算应力刚度矩阵。
屈曲分析因为计算出的特征值表示屈曲荷载系数,所以一般施加单位载荷,但是ANSYS 程序的特征值限值是1 000 000,所以如果求解的特征值超过此限值的话,我们应当施加一个比较大的荷载,而不应该还是用单位力施加!计算结果为:临界力=施加力×求出来的屈曲系数,所以一般的用单位力施加的临界力就等于求出的屈曲系数,因为施加力为1。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
屈曲分析的过程说明:
屈曲分析是一种用于确定结构开始变得不稳定时的临介荷载和屈
曲结构发生屈曲响应时的模态形状的技术。
ANSYS提供两种结构屈曲荷载和屈曲模态分析方法:非线性屈曲分析和特征值屈曲分析。
非线性屈曲分析是在大变形效应开关打开的情况下的一种非线性
静力学分析,该分析过程一直进行到结构的极限荷载或最大荷载。
非
线性屈曲分析的方法是,逐步地施加一个恒定的荷载增量,直到解开
始发散为止。
尤其重要的是,要一个足够小的荷载增量,来使荷载达
到预期的临界屈曲荷载。
若荷载增量太大,则屈曲分析所得到的屈曲
荷载就可能不准确,在这种情况下打开自动时间步长功能,有助于避
免这类问题,打开自动时间步长功能,ANSYS程序将自动寻找屈曲荷载。
特征值屈曲分析步骤为:
1.建模
2.获得静力解:与一般静力学分析过程一致,但必须激活预应
力影响,通常只施加一个单位荷载就行了
3.获得特征屈曲解:
A.进入求解
B.定义分析类型
C.定义分析选项
D.定义荷载步选项
E.求解
4.扩展解
之后就可以察看结果了
示例1:
!ansys7.0有限元分析实用教程
!3.命令流求解
!ANSYS命令流:
!Eigenvalue Buckling
FINISH!这两行命令清除当前数据/CLEAR
/TITLE,Eigenvalue Buckling Analysis
/PREP7!进入前处理器
ET,1,BEAM3!选择单元
R,1,100,833.333,10!定义实常数
MP,EX,1,200000!弹性模量
MP,PRXY,1,0.3!泊松比
K,1,0,0!创建梁实体模型
K,2,0,100
L,1,2!创建直线
ESIZE,10!单元边长为1mm
LMESH,ALL,ALL!划分网格
FINISH!退出前处理
!屈曲特征值部分
/SOLU!进入求解
ANTYPE,STATIC!在进行屈服分析之前,ANSYS需要从静态分析提取数据
PSTRES,ON!屈服分析中采用预应力
DK,1,ALL!定义约束
FK,2,FY,-1!顶部施加载荷
SOLVE!求解
FINISH!退出求解
/SOLU!重新进入求解模型进行屈服分析
ANTYPE,BUCKLE!屈服分析类型
BUCOPT,LANB,1!1阶模态,子空间法
SOLVE!求解
FINISH!退出求解
/SOLU!重新进入求解展开模态
EXPASS,ON!模态展开打开
MXPAND,1!定义需要展开的阶数
SOLVE!求解
FINISH!退出求解
/POST1!进入通用后处理SET,LIST!列出特征值求解结果SET,LAST!入感兴趣阶数模态结果PLDISP!显示变形后图形
!NonLinear Buckling!非线性分析部分
FINISH!这两行命令清除当前数据
/CLEAR
/TITLE,Nonlinear Buckling Analysis
/PREP7!进入前处理
ET,1,BEAM3!选择单元
MP,EX,1,200000!弹性模量
MP,PRXY,1,0.3!泊松比
R,1,100,833.333,10!定义实常数
K,1,0,0,0!底端节点
K,2,0,100,0!顶点
L,1,2!连成线
ESIZE,1!网格尺寸参数设定
LMESH,ALL!划分网格
FINISH!退出前处理
/SOLU!进入求解
ANTYPE,STATIC!静态分析类型(非屈服分析)NLGEOM,ON!打开非线性大变形设置OUTRES,ALL,ALL!选择输出数据
NSUBST,20!5个子步加载
NEQIT,1000!20步迭代
AUTOTS,ON!自动时间步长
LNSRCH,ON!激活线搜索选项
/ESHAPE,1!显示二维状态下变形图
DK,1,ALL,0!约束底部节点
FK,2,FY,-50000!顶部载荷稍微比特征值分析结果大FK,2,FX,-250!施加水平扰动载荷
SOLVE!求解
FINISH!退出求解
/POST26!进入时间-历程后处理器RFORCE,2,1,F,Y!2#变量表示力
NSOL,3,2,U,Y!3#变量表示y方向位移
XVAR,2!将x轴显示2#变量
PLVAR,3!y轴显示3#变量数据
/AXLAB,Y,DEFLECTION!修改y轴标签
/AXLAB,X,LOAD!修改x轴标签
/REPLOT!重新显示图形
示例2:
!悬臂梁受端部轴向压力作用的屈曲分析
!先进行静力分析,在进行特征值屈曲分析,最后进行非线性分析
!静力分析
fini
/cle
/filname,beam-flexure
/tittle,beam-flexure
/prep7
*set,f1,-1e6!设置轴向压力荷载参数
et,1,beam189
mp,dens,1,7.85e3!设置材料参数
mp,ex,1,2.06e11
mp,nuxy,1,0.2
sectype,1,beam,I,,2!设置截面参数
secoffset,cent
secdata,0.15,0.15,0.25,0.015,0.015,0.015,0,0,0,0
k,1,0
k,2,2.5,0
k,3,1.25,1
lstr,1,2
latt,1,,1,,3,,1
lesize,1,,,10
lmesh,1
/view,1,1,1,1
/eshape,1.0
dk,1,,,,0,all,
fk,2,fx,f1!施加关键点压力
finish
/solu
antype,0
eqslv,spar!求解器设置稀疏矩阵直接法
pstres,on!打开预应力开关
solve
finish
!特征值屈曲分析
/solu
antype,1
bucopt,lanb,6,0!取前六阶模态分析
mxpand,6,0,0,1,0.001
solve
finish
/post1
*do,i,1,6
set,i
pldisp,i
*enddo
*get,freq1,mode,1,freq
finish
!非线性屈曲分析
/config,nres,200!只记录两百步的结果
/prep7
tb,biso,1,1,2!定义材料非线性
tbtemp,0
tbdata,,2.0e8,0
upgeom,0.01,1,1,'beam-flexure','rst'
!对有限元模型进行一阶模态的位移结果0.01倍的修改
save,beam-flexure,db
finish
resu,beam-flexure,db
/solu
antype,0
nlgeom,1!打开大变形
outres,all,all
arclen,1,0!弧长法设置
arctrm,l!弧长法终止准则达到第一个峰值时终止计算
nsubst,200,,,1
fk,2,fx,f1*freq1
!fk,2,fx,f1*freq1*1.2!将轴向压力值放大,放大系数为第一阶模态的主频
solve
finish。