Workbe nch屈曲分析1、基础概念结构在载荷作用下由于材料弹性性能发生变形,若变形后结构上的载荷保持平衡,这种状态称为弹性平衡。
如果结构在平衡状态时,受到扰动而偏离平衡位置,当扰动消除后仍能恢复原来平衡状态,这种平衡状态称为稳定平衡状态,反之,如果受到扰动而偏离平衡位置,即使扰动消除,结构仍不能恢复原来的平衡状态,而结构在新的状态下平衡,则原来的平衡状态就成为不稳定平衡状态。
当结构所受载荷达到某一值时,若增加一微小的增量,则结构平衡状态将发生很大的改变,这种现象叫做结构失稳或结构屈曲。
根据失稳的性质,结构稳定问题可分为以下三类:第一类失稳是理想化情况,即达到某个载荷时,除结构原来的平衡状态存在外,出现第二个平衡状态,故又叫做平衡分叉失稳,数学上就是求解特征值问题,又叫做特征值屈曲分析。
第二类失稳是结构失稳,变形将大大发展,而不会出现新的变形形式,即平衡状态不发生质变,也叫极顶失稳,结构失稳时,相应载荷叫做极限载荷,理想结构或完善结构不存在,总是存在这样那样的缺陷,大多数问题属于第二类失稳问题。
第三类失稳是当在和达到某值时,结构平衡状态发生一明显跳跃,突然过渡到非临近的另一具有较大位移的平衡状态,称为跳跃失稳,跳跃失稳没有平衡分叉点,也没有极值点,如坦拱、扁壳、二力杆的失稳都属于此类。
结构弹性稳定分析属于第一类失稳对应workbench的线性特征值分析(Eigenvalue Buckling),考虑缺陷,非线性影响的第二类结构属于workbe nch的非线性特征值分析( Eige nvalue Buckling),第三类的失稳对应workbench的Static Structural,无论前屈曲平衡状态或后屈曲平衡状态均可一次计算求出,即全过程分析。
1.1屈曲分析基础理论在平衡状态,考虑到轴向力或中面内力对弯曲变形的影响,根据势能驻值原理得到结构平衡方程为kJ K G〕U—p:式中K E 1为结构弹性刚度矩阵,K G I为结构几何刚度矩阵,也称为初应力刚度矩阵,<U '为节点位移向量;"P*为节点载荷向量,上式也为几何非线性分析平衡方程。
为得到随遇平衡状态,应是系统势能的二阶变分为零。
即:kJ - 'K J lu因此必有:k J ■ kJ = o式中结构弹性刚度矩阵已知,结构外载荷也就是要求得屈曲载荷未知,结构几何刚度矩阵未知,为了求得该屈曲载荷,假设有一组载荷P0 \对应的几何刚度矩阵为K G [并假定屈曲时的载荷是P0的■倍,固有’k; = k G1,上式可变为K E 1- kG 丨-0写成特征值的方式为k E「%.W o式中■为第■阶的特征值,「匚为■对应的特征向量,是该阶载荷下结构的变形形状, 曲模态或失稳即屈模态。
在workbench中计算出的是•和即屈曲载荷系数和模态,而屈曲载荷为■2.1、Linear-based Eigenvalue Buckling Analysis线性屈曲分析应注意以下几点线性屈曲分析只能在静力分析模块中定义边界通过特征值屈曲分析计算的结果是在静态结构分析中应用所有载荷的屈曲载荷因子。
例如,如果在静态分析中对结构应用10 N压缩负载,如果特征值屈曲分析计算负载系数为1500,则预测的屈曲载荷为1500 X 10 = 15000N。
因此,在屈曲分析之前的静态分析中应用单位载荷是一种典型的方法。
在静态分析中所使用的所有载荷都适用于屈曲负载系数请注意,负载系数表示所有负载的比例因子。
如果某些负载是恒定的(例如,自重重力负载),而其他负载是可变的(例如,外部施加的负载),则需要采取特殊步骤确保准确的结果。
为了实现这一目,可以使用一个策略,就是是迭代特征值,调整可变载荷,直到屈曲因子变为1.0 (或接近1.0,在一些收敛公差内)③ I A= 111A= 100Load Factor = 100 F= 100 + 100 W o Load Factor = 11:F = 110+ 1 1 W oLoad Factor = 0 99: F = 110 + 0 99 W o特征值屈曲分析案例材料:结构钢模型:r=1mm L=50mm的圆柱边界:一端固定,一端施加10N集中力。
1.创建分析系统▼F-----------1-----------1▼1芳Static Structural1F Eigenvalue Budding2療Engheefing Data ■/ 上-12• En^neerino Data / j 3* Geometry /3* Geametry ? j4& Model V『-----------1K常Model g j5海Setup 'f15电Setup 習j6yy Solutiori V j'6鎰| Solution 审j 7• Results? j70 Results V jStatic Structural Eigenvalue Buckling首先创建一个结构静力分析分析系统,再创建特征值分析系统将他们数据共享。
2、静力分析边界D: Static StruiftuirailStatic Structural Tim奚L冷[A] Fixed Support■ FurceMO^N3. 求解静力分析OaO 15 JOO 帥加处mm)7.500 22500D: :anicliirilTtrtpl D@fo-nnditionTyp&: T ata 口e白匸DHUnit; mmTime1:14^12339 Max0J»3 959660J MB2257cjMEestrojsoMsas(JJT41H2SOJ«DZ741^oz节wOMin3.求解特征值E:KgmwBus BiK klingTtrtpl D@fo-nnditionTy口殳带T益1•白DigEriEi占-D門Lcdd ^ulti jlier [Liic^rj:lb.^34畑徉gKI 57 1^124 0^57760.72 31S057B52Q4汩旳仏2朋细C.1M&3 OMin在总变形中可以查看1阶变形模态和1阶特征值,可以看出一阶特征值为15.534,则屈曲『00 2I D J3D*O,CD (mrrt)10J)0 .304)0D J OO I?15^CD »^0D (mrri)750C 22300载荷为10*15.534=155.34N ,如果将静力分析中集中力改为155.34 ,计算出特征值为约等于0.9997,1,集中力155.34就是此结构的屈曲载荷。
2.2、Nonlinear-based Eigenvalue Buckling Analysis非线性屈曲分析要点至少有一个非线性属性在静力分析中被定义。
除了在静态结构分析中定义的荷载之外,还必须在屈曲分析中至少定义一个载荷来进行求解。
要启用此功能,将保持预应力加载模式”属性设置为是”(默认设置)将在特征值屈曲”分析中保留静态结构分析中的加载模式。
将属性设置为否需要您定义特征值屈曲分析的新加载模式。
这种新的加载模式可以与预应力分析完全不同在基于非线性的特征值屈曲分析中,负载乘数仅对屈曲分析中的负载进行了缩放。
在估计结构的极限屈曲载荷时,必须考虑静态结构中的载荷和特征值分析。
用于计算非线性特征值屈曲的极限屈曲载荷的方程是F BUCKLING = F RESTART + 入• P E RTRUB 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列示例那么直接,这是因为FRESTAR和FPERTRU的值基本上是分别在静态和屈曲分析中的有效载荷值。
举个例子,一个悬臂梁的理论极限弯曲强度为1000N,它受到了影响对250 n的压缩力(a)。
根据负载因素计算极限屈曲载荷(F)的过程用力学方法对线性和非线性特征值屈曲分析进行了计算,如下图所示示意图A = 25dN^=250NStatkB = 2MN Budkling777TUtiear 歸总A/i^lysis Force in $4atk=2S0M Lo»dl Factor = jIF =4X250 = 1DDON Ikiear Eto Anal阀$ r^rcebSutic = 250Nloadl Fade =也F =4X250 = 1DOONNonlinear 加电Force in 5tflUc = 250NKw Pre^tress LMd-^atternForte in Bucklimg = 250NLoad Factor = 3F = 2W* JX250=1000NNonlir^flr Base 人窗備Fwce in Static =25ONK&ep P『L$tre幻Load冲吐term =他Force m Budkling = 250Nlo聞Factcw - 3F= 250 * 3 X250 = 1DOONNofllinor Base AnalysisReep Pre-stress LoadH1沉阳n = NoForce in SUtic4A)=250NForce in Budklimg (B| = 2 SONload F*NX.A. = 3a 25D*3X2W=1DDON非线性屈曲分析案例分析模型与前文特征值分析一样,只是在求解设置中打开大变形开关1.静力分析结果0: Statk StruduralOffgrTTiazionType; T De耳口Ei由i:onLin rt: mmTirnv:1CJO079CQ6Max 。
加网湖0JQD&19150J>3D53Q71Oj&304422fiOJXI0359BC J»0D2*535fljNDH769a3451e-5Q Mln Normal La^nge Contact (rKxiIirHar)D J OO I?IS^CD MJOD (mnri)T50C 225f02.3.屈曲分析设置Details of Analysis Settings'- OptionsMax Modes to Find2. Keep Pr^-Str«5 Load-Pattern Vec=Solver Controls Salve-r Type Proq ram Control led设置中Keep Pre-Stress Load-Pattern 为YES 这时只能在静力分析中施加载荷,不允许在屈曲 分析中施加载荷。