Workbench 屈曲分析
1、基础概念
结构在载荷作用下由于材料弹性性能发生变形, 若变形后结构上的载荷保持平衡, 这种状态称为 弹性平衡 。如果结构在平衡状态时, 受到扰动而偏离平衡位置, 当扰动消除后仍能恢复原来平衡状态, 这种平衡状态称为 稳定平衡状态 ,反之,如果受到扰动而偏离平衡位置,即使扰动消除, 结构仍不能恢复原来的平衡状态, 而结构在新的状态下平衡, 则原来的平衡状态就成为 不稳定平衡状态 。
当结构所受载荷达到某一值时, 若增加一微小的增量, 则结构平衡状态将发生很大的改变,这种现象叫做 结构失稳或结构屈曲 。
根据失稳的性质,结构稳定问题可分为以下三类:
第一类失稳是理想化情况, 即达到某个载荷时, 除结构原来的平衡状态存在外, 出现第二个平衡状态, 故又叫做 平衡分叉失稳 ,数学上就是求解特征值问题, 又叫做特征值屈曲分析。
第二类失稳是结构失稳, 变形将大大发展, 而不会出现新的变形形式, 即平衡状态不发生质变,也叫 极顶失稳 ,结构失稳时, 相应载荷叫做极限载荷, 理想结构或完善结构不存在,总是存在这样那样的缺陷,大多数问题属于第二类失稳问题。
第三类失稳是当在和达到某值时, 结构平衡状态发生一明显跳跃, 突然过渡到非临近的
另一具有较大位移的平衡状态,称为 跳跃失稳 ,跳跃失稳没有平衡分叉点,也没有极值点,如坦拱、扁壳、二力杆的失稳都属于此类。
结构弹性稳定分析属于第一类失稳对应 workbench 的线性特征值分析 ( Eigenvalue Buckling),
考虑缺陷,非线性影响的第二类结构属于 workbench 的非线性特征值分析( Eigenvalue
Buckling),第三类的失稳对应 workbench 的 Static Structural ,无论前屈曲平衡状态或后屈曲
平衡状态均可一次计算求出,即全过程分析。
1.1 屈曲分析基础理论
在平衡状态, 考虑到轴向力或中面内力对弯曲变形的影响,平衡方程为
根据势能驻值原理得到结构
K E
K G
U
P
式中
K E
为结构弹性刚度矩阵,
K G
为结构几何刚度矩阵,也称为初应力刚度矩阵,
U
为节点位移向量; P 为节点载荷向量,上式也为几何非线性分析平衡方程。
为得到随遇平衡状态,应是系统势能的二阶变分为零。即:
K E
K G
U
0
因此必有:
K E
K G
0
式中结构弹性刚度矩阵已知,
结构外载荷也就是要求得屈曲载荷未知,
结构几何刚度矩
阵未知, 为了求得该屈曲载荷,
假设有一组载荷
P0
,对应的几何刚度矩阵为
K G0
,并假定屈
曲时的载荷是
P0
的
倍,固有
K G0
= K G
,上式可变为
K E
K G0
0
写成特征值的方式为
K E
K G
0
式中
为第 阶的特征值,
为 对应的特征向量,是该阶载荷下结构的变形形状,
即屈
曲模态或失稳模态。
在
workbench
中计算出的是
和
,即屈曲载荷系数和模态,而屈曲载荷为
P0
.
2.1 、 Linear-based Eigenvalue Buckling Analysis
线性屈曲分析应注意以下几点
线性屈曲分析只能在静力分析模块中定义边界
通过特征值屈曲分析计算的结果是在静态结构分析中应用所有载荷的屈曲载荷因
子。例如,如果在静态分析中对结构应用 10 N压缩负载, 如果特征值屈曲分析计算
负载系数为 1500,则预测的屈曲载荷为 1500 × 10 = 15000N。因此,在屈曲分析之
前的静态分析中应用单位载荷是一种典型的方法。
在静态分析中所使用的所有载荷都适用于屈曲负载系数
请注意,负载系数表示所有负载的比例因子。 如果某些负载是恒定的(例如,自重
重力负载),而其他负载是可变的 (例如, 外部施加的负载) ,则需要采取特殊步骤
确保准确的结果。
为了实现这一目,可以使用一个策略,就是是迭代特征值,调整可变载荷, 直到屈
曲因子变为 1.0(或接近 1.0,在一些收敛公差内)
特征值屈曲分析案例
材料:结构钢
模型: r=1mm L=50mm 的圆柱
边界:一端固定,一端施加 10N集中力。
创建分析系统
首先创建一个结构静力分析分析系统,再创建特征值分析系统将他们数据共享。
2、静力分析边界
求解静力分析
求解特征值
在总变形中可以查看 1阶变形模态和 1阶特征值,可以看出一阶特征值为 15.534,则屈曲
载荷为 10*15.534=155.34N ,如果将静力分析中集中力改为
155.34 ,计算出特征值为
0.9997,
约等于 1,集中力 155.34就是此结构的屈曲载荷。
2.2 、 Nonlinear-based Eigenvalue Buckling Analysis
非线性屈曲分析要点
至少有一个非线性属性在静力分析中被定义。
除了在静态结构分析中定义的荷载之外, 还必须在屈曲分析中至少定义一个载荷来
进行求解。要启用此功能,将 “保持预应力加载模式 ”属性设置为 “是 ”(默认设置)
将在 “特征值屈曲 ”分析中保留静态结构分析中的加载模式。 将属性设置为否需要您
定义特征值屈曲分析的新加载模式。 这种新的加载模式可以与预应力分析完全不同
在基于非线性的特征值屈曲分析中,负载乘数仅对屈曲分析中的负载进行了缩放。
在估计结构的极限屈曲载荷时,
必须考虑静态结构中的载荷和特征值分析。
用于计
算非线性特征值屈曲的极限屈曲载荷的方程是
FBUCKLING = F RESTART + λi· PERTRUBF where:
FBUCKLING = The ultimate buckling load for the structure.
FRESTART = Total loads in Static Structural analysis at the specified restart load step.
λ i = Buckling load factor for the "i'th" mode.
FPERTRUB = Perturbation loads applied in buckling analysis
例如:如果在静力分析中施加
100N集中力,在屈曲分析中加
10N力,你得到载荷因
子位 15,则结构的极限屈曲力位
100+( 15*10 ) =250
注意:可以使用一维柱的屈曲来验证上述方程的极限屈曲载荷。 然而,对于在静态
结构和特征值屈曲分析中应用的不同负载组合,计算 2D和 3D问题的极限屈曲载荷
可能不如
1D列示例那么直接,这是因为
FRESTART和 FPERTRUB的值基本上是分别在
静态和屈曲分析中的有效载荷值。
举个例子,一个悬臂梁的理论极限弯曲强度为
1000N,它受到了影响
对 250 n的压
缩力 (a)。根据负载因素计算极限屈曲载荷 (F)的过程 用力学方法对线性和非线性特
征值屈曲分析进行了计算,如下图所示 示意图
非线性屈曲分析案例
分析模型与前文特征值分析一样,只是在求解设置中打开大变形开关。
1.静力分析结果
2.屈曲分析设置
设置中 Keep Pre-Stress Load-Pattern 为 YES,这时只能在静力分析中施加载荷,不允许在屈曲分析中施加载荷。
屈曲分析结果
在总变形中可以查看 1 阶变形模态和 1 阶特征值, 可以看出一阶特征值为 14.535,则屈曲载荷为 10+10*14.534=155.34N ,如果将静力分析中集中力改为 155.34,计算出特征值为 7.877e-4 155.34+155.34*7.877e-4 155.34,集中力 155.34 就是此结构的屈曲载荷。
若屈曲分析设置中 Keep Pre-Stress Load-Pattern 为 NO,此时在静力分析中可以施加载荷,允屈曲分析中也可以施加载荷。屈曲分析中的载荷只能施加节点力。
首先在施加集中力的端面创建一个节点集合
在屈曲分析中设置集中力
屈曲分析结果
在总变形中可以查看 1 阶变形模态和 1 阶特征值, 可以看出一阶特征值为 7.2672,则屈曲载荷为 10+20*7.2672=155.34N 。
推荐访问:计划总结 屈曲 计划 解析 Workbench屈曲解析总结计划总结计划