[力学][有限元][FEM]Basics of Buckling Analysis [曲屈分析基础] (《有限单元法-编程与软件应用》章节节选)

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。

《有限单元法:编程与软件应用》 曲屈分析 章节的内容节选:

11.1 稳定问题分类

结构失稳(屈曲)是指在外力作用下结构的平衡状态开始丧失,稍有扰动变形便迅速增大,最后使结构发生破坏。稳定问题一般分为两类,第一类是理想化的情况,即达到某种荷载时,除结构原来的平衡状态存在外,还可能出现第二个平衡状态,所以又称平衡分岔失稳分支点失稳 (Buckling, Bifuraction)(图1中OAB曲线),对应于数学中是求解特征值问题,故又称特征值屈曲,此类结构失稳时相应的荷载称为屈曲荷载。第二类是结构失稳时,变形将迅速增大,而不会出现新的变形形式,又称极值点失稳(Instability)(图1中OCD曲线),结构失稳时相应的荷载称为极限荷载。此外还有一种跳跃失稳(Snap-through)(图1中OEF曲线),当荷载达到某值时(E点),结构平衡状态发生一个明显的跳跃,突然过渡到非临近的另一个具有较大位移的平衡状态(F点),由于在跳跃时结构通常已经破坏,其后的状态不能被利用,所以可归入第二类失稳。

图 11‑1 稳定问题

11.2 最小势能原理

11.2.1 系统的平衡

若系统处于平衡状态,则一定存在一种位移状态,使得系统总势能为驻值,即

              (11.2‑1)

公式中

                (11.2‑2)

为系统总势能,为应变能,为外力功势能,数值上等于外力在系统变形中做的功取负值,三者都是系统位移的函数。

11.2.2 稳定平衡

如果要研究一个平衡是否稳定,则需要进一步考虑总势能的二阶变分条件,即

            (11.2‑3)

工程中感兴趣的是临界平衡位置所确定的临界荷载和相应的屈曲形态,因此,总势能对位移的二阶导数为正的条件是研究弹性稳定性的理论基础。

11.2.3 弹性稳定问题的有限元列式

分析稳定问题时,需要在总势能的计算中考虑位移的高阶项影响,此时应变能U可以分为两部分,分别是应变能的线性部分和非线性部分,即

         (11.2‑4)

外力功势能 ,则势能驻值方程可以写成

        (11.2‑5)

将应变能的线性部分和非线性部分分别写作

            (11.2‑6)

则势能驻值方程可进一步写成

             (11.2‑7)

上式中,K0 一般是不变的,KNL 与系统的位移及应力状态有关,因此该式是一个非线性方程。如果以未变形的初始位移为参考,则KNL 只与系统的应力有关,此时这种仅由应力状态所决定的单元刚度矩阵叫做几何刚度矩阵,或初应力矩阵,用表示。

11.2.4 弹性稳定问题的求解

弹性稳定问题可分为两类,一类是欧拉稳定问题,另一类是极值稳定问题。

欧拉稳定问题是指系统的初应力状态处于某种临界状态时,对于临界位移的任何挠动都可能使系统丧失稳定性。如果把欧拉失稳的临界状态看成是一种初应力状态,则欧拉稳定问题也可用公式 (11.2‑7) 求解。具体为:根据某一初始应力状态{σ} ,求得几何刚度矩阵[Kσ],同理易知当初始应力状态为λ{σ} 时,几何刚度矩阵为 λ[Kσ],则公式(11.2‑7) 可写作

       (11.2‑8)

假设此时达到临界状态,则存在一个挠动的位移状态,使得系统在外力不变的条件下也处于平衡状态,即

       (11.2‑9)

两式相减可得

          (11.2‑10)

此时,欧拉稳定问题,归结为求解一个广义特征值问题。

上式是齐次线性代数方程组,其有非零解的条件是系数行列式等于0,即

       (11.2‑11)

将公式(11.2‑10)变换为标准特征值问题,则有

          (11.2‑12)

式中广义特征值为 λi,对应屈曲分析中的屈曲因子,特征向量为 Δ{δi} ,对应屈曲分析中的屈曲模态。

对于标准特征值问题,可采用上一章模态分析中的方法进行求解,此处不再赘述。

极值稳定问题是在荷载逐渐增大的过程中达到某种状态导致变形无限增长而丧失稳定性,在非线性有限元分析中,必须用增量法求解极值稳定性问题,即是求下式的渐进解

        (11.2‑13)

极值稳定性问题一般需要考虑[Kσ]随单元应力状态和单元位移状态的变化,采用增量法求解,此处留作以后讨论。

  • 相关内容 ( Related Topics )

[00]. [Book][书] 有限单元法:编程与软件应用

[01]. 2D Truss by Matlab [Matlab 2D桁架有限元] (code by myself)

[02]. 3D Truss by MATLAB[MATLAB 3D 桁架静力分析] (code by myself)

[03]. [FEM][有限元][编程][Matlab][Code by myself] FEM Analysis: 2D Truss Element [有限元分析: 2D桁架单元]

[04]. [FEM][有限元][编程][Matlab][Code by myself] FEM Analysis: 2D Euler Beam Element [有限元分析: 2D欧拉梁单元]

[05]. [FEM][有限元][编程][Matlab][Code by myself] 2D剪切梁单元

[06]. [FEM][有限元][编程][Matlab][Code by myself] 2D Timoshenko梁单元

[07]. [FEM][有限元][编程][Matlab][Code by myself] 三角形常应变单元(CST)

[08]. [FEM][有限元][编程][Matlab][Code by myself] 平面4变形单元(Q4)

[09]. [FEM][有限元][编程][Matlab][Code by myself] 4节点四面体单元(TET4)

[10]. [FEM][有限元][编程][Matlab][Code by myself] 8节点六面体单元(C3D8)

[11]. [FEM][有限元][编程][Matlab][Code by myself] 平面8节点二次“完全积分”单元(CPS8)

[12]. [FEM][有限元][编程][Matlab][Code by myself] 平面6节点二次“完全积分”单元(CPS6)

[13]. [FEM][MATLAB][有限元] FEM Modal Analysis Programming with MATLAB (Truss Element) (桁架单元模态分析编程)

[14]. [FEM][MATLAB][有限元] FEM Modal Analysis Programming with MATLAB (Frame Elements) (框架单元模态分析编程)

[15]. [力学][有限元]Basics of Linear Buckling Analysis [线性曲屈分析基础]

[16]. [FEM][MATLAB][有限元][编程] FEM Buckling Analysis Programming with MATLAB (Truss Elements) (桁架单元曲屈分析编程)

[17]. [FEM][MATLAB][有限元][编程] FEM Buckling Analysis Programming with MATLAB (Frame Elements) (框架单元曲屈分析编程)


4 0

  • 注释 ( Comments )

  ( 如果您发现有错误,欢迎批评指正。邮箱:jidong_cui@163.com . 如果您喜欢这篇博文,请在上面给我 点个赞 吧! :-)   :-)   

  ( If you found any mistakes in the post, please let me know. Email : jidong_cui@163.com. If you like this posts, please give me a thumbs up rating on the above button! )

  • 微信公众号 ( Wechat Subscription)

WeChat_QRCode

欢迎关注 “结构之旅” 微信公众号

4 thoughts on “[力学][有限元][FEM]Basics of Buckling Analysis [曲屈分析基础] (《有限单元法-编程与软件应用》章节节选)

  1. Pingback: [FEM][MATLAB][有限元][编程] FEM Buckling Analysis Programming with MATLAB (Frame Elements) (框架单元曲屈分析编程) | 崔济东的博客 - JiDong Cui's Blog

  2. Pingback: [FEM][MATLAB][有限元][编程] 压杆稳定问题MATLAB有限元编程 (《有限单元法-编程与软件应用》章节节选) | 崔济东的博客 - JiDong Cui's Blog

  3. Pingback: [FEM][MATLAB][有限元][编程] FEM Buckling Analysis Programming with MATLAB (Truss Elements) (桁架单元曲屈分析编程) | 崔济东的博客 - JiDong Cui's Blog

  4. Pingback: [动力学] 模态分析基础 (Basics of modal analysis) | 崔济东的博客 - JiDong Cui's Blog

Leave a Reply

Your email address will not be published. Required fields are marked *

:wink: :-| :-x :twisted: :) 8-O :( :roll: :-P :oops: :-o :mrgreen: :lol: :idea: :-D :evil: :cry: 8) :arrow: :-? :?: :!:

This site uses Akismet to reduce spam. Learn how your comment data is processed.