搜索
第11章屈曲分析11.1屈曲分析概述静力分析方法认为杆件的破坏取决于材料的强度，当杆件承受的应力小于其许用应力时，杆件便可安全工作，对于细长受压杆件这却并不一定正确。压杆在承受的应力小于其许用应力时，杆件会发生变形而失去承载能力，这类问题称为压杆屈曲问题，或者压杆失稳问题。工程中许多细长构件如发动机中的连杆、液压缸中的活塞杆和订书机中的订书针等，以及其他受压零件，如承受外压的薄壁圆筒等，在工作的过程中，都面临着压杆屈曲的问题。临界载荷是受压杆件承受压力时保持杆件形状的载荷上限。压杆承受临界载荷或更大载荷时会发生弯曲，如图11-1所示。经典材料力学使用Euler公式求取临界载荷：22lEJFcr（11-1）图11-1临界载荷下压杆发生屈曲该公式在长细比超过100有效。针对不同的压杆约束形式，参数的取值如表11-1所示。表11-1Euler公式中参数的取值约束情况一端固定，一端自由一端固定，一端绞支两端绞支两端固定两端固定，一端可横向移动群20.710.51对于压杆屈曲问题，ANSYS中一方面可以使用线性分析方法求解Euler临界载荷，另一方面可以使用非线性方法求取更为安全的临界载荷。ANSYS提供两种技术来分析屈曲问题，分别为非线性屈曲分析法和线性屈曲分析法（也称为特征值法）。因为这两种方法的结果可能截然不同（见图11-2），故需要理解它们的差异：非线性屈曲分析法通常较线性屈曲分析法更符合工程实际．使用载荷逐渐增大的非线性静力学分析，来求解破坏结构稳定的临界载荷。使用非线性屈曲分析法，甚至可以分析屈曲后的结构变化模式。线性屈曲分析法可以求解理想线性弹性理想结构的临界载荷，其结果与Euler方程求得的基本一致。图11-2不同分析方法的屈曲分析结果11.2线性屈曲分析步骤由于线性屈曲分析基于线性弹性理想结构的假设进行分析，所以该方法的结果安全性不佳，那么在设计中不宜直接采用分析结果。线性屈曲分析包含以下步骤。11.2.1前处理建立模型，包括：（1）定义单元类型，截面结构、单元常数等。在线性屈曲分析中，ANSYS对单元采取线性化处理，故即使定义了非线性的高次单元，在运行中也将被线性化处理。（2）定义材料，可以采用线性各向同性或线性正交各向异性材料，因求解刚性矩阵的需要，必须定义材料的杨氏模量。（3）建立有限元模型，包括几何建模与网格化处理。11.2.2求取静态解求取静态解，包括：（1）进入求解器，并设定求解类型为Static。（2）激活预应力效应（在求解过程中必须激活）。命令方式：PSTRES,ONGUI方式：选择MainMenuSolutionAnalysisTypeAnalysisOptions命令，找到PSTRES并选中，将其设置为打开状态。（3）施加约束和载荷：可以施加一个单位载荷，也可取一个较大的载荷（特别在求解模型的临界载荷很大时）。（4）求解并退出求解器。11.2.3求取屈曲解求取临界载荷值和屈曲模态，包括：（1）进入求解器，并设定求解类型为EigenBuckling。命令方式：ANTYPE,BUCKLEGUI方式：选择MainMenuSolutionAnalysisType-NewAnalysis命令，在弹出的对话框中，将EigenBuckling前的单选框选中。（2）设置求解选项。命令方式：BUCOPT,Method,NMODE,SHIFT,LDMULTE,RangeKey其中：Method指定临界载荷提取的方法，可为LAMB指定BlockLanczos方法，或SUBSP指定子空间迭代法。NMODE指定临界载荷提取的数目。SHIFT指定临界载荷计算起始点，默认为0.0。LDMULTE指定临界载荷计算终止点，默认为正无穷。RangeKey控制特征值提取方法的计算模式，可为CENTER或RANGE;默认为CENTER，计算范围为（SHIFTLDMULTE，SHIFT+LDMULTE），采用RANGE的计算范围为（SHIFT,LDMULTE）。GUI方式：选择MainMenuSolutionAnalysisTypeAnalysisOptions命令，在弹出的对话框中，输入命令中的各项参数。（3）设置载荷步骤、输出选项和需要扩展的模态。扩展模态的方式如下。命令方式：MXPAND,NMODE,FREQB,FREQE,Elcalc,其中：NMODE指定需要扩展的模态数目，默认为ALL，扩展求解范围内的所有模态。如果为-1,不扩展模态，而且不将模态写入结果文件中。FREQB指定特征值模态扩展的下限，如果与FREQE均默认，则扩展并写出指定求解范围内的模态。FREQE指定特征值模态扩展的上限。Elcalc网格单元计算开关，如果为NO，则不计算网格单元结果、相互作用力和能量等结果；如果为YES，计算网格单元结果、相互作用力、能量等；默认为NO。SIGNIF指定阈值，只有大于阈值的特征值模态才能被扩展。MSUPkey指定网格单元计算结果是否写入模态文件中。GUI方式：选择MainMenuSolutionLoadStepOptsExpansionPassSingleExpandExpandModes命令，在弹出的对话框中，输入命令中的各项参数。11.2.4后处理查看结果。（1）查看特征值。（2）查看屈曲变形图。11.3非线性屈曲分析步骤非线性屈曲分析属于大变形的静力学分析，在分析中将压力扩展到结构承受极限载荷。如果使用塑性材料，结构在承受载荷时可能会发生其他非线性效应，如塑性变形等。从图11-2中可以看到，使用非线性屈曲分析方法得到的临界载荷一般较线性方法小，因此在非线性分析中通常使用线性分析中的临界载荷为加载起点，分析结果出现屈曲后的变化形态。11.3.1前处理建立模型，包括：（1）定义单元类型、截面结构、单元常数等。（2）定义材料，可以采用线性各向同性或线性正交各向异性材料，因求解刚性矩阵的需要，必须定义材料的杨氏模量。（3）建立有限元模型，包括几何建模与网格化处理。11.3.2加载与求解加载并求解，包括：（1）进入求解器，并设定求解类型为static。（2）激活大变形效应。命令方式：NLGEOM,ONGUI方式：选择MainMenuSolutionAnalysisTypeSol'sControl命令，弹出SolutionControls对话框，在对话框中的AnalysisOption框下选择LargeDisplacementStatic项。（3）设置子载荷的时间步长。使用非线性屈服分析方法是逐渐增大载荷直到结果开始发散，如果载荷增量过大，得到的分析结果可能不准确。打开二分法选项和自动时间步长选项有利于避免这样的问题。打开自动时间步长选项时，程序自动求出屈服载荷。在求解时，一旦时间步长设置过大导致结果不收敛，程序将自动二分载荷步长，在小的步长下继续求解，直到能获得收敛结果。在屈曲分析中，当载荷大于等于屈曲临界载荷时，结果将不收敛。一般而言，程序将收敛到临界载荷。（4）施加约束和载荷，可从小到大依次逐步将载荷施加到模型上，不要一次施加过大的载荷，以免在求解过程中出现不收敛的现象。在施加载荷时，施加一个小的扰动，使结构屈曲发生。（5）求解并退出求解器。11.3.3后处理查看结果，包括：（1）进入通用后处理器查看变形。（2）进入时间历程后处理器查看参数随时间的变化等。11.4中间铰支增强稳定性线性分析11.4.1问题描述与分析问题描述：两端铰支的细长杆在承受压力时容易发生失稳线性（屈曲效应），工程上为了提高细长杆的稳定性，常在杆中间增加铰支提高杆的抗屈曲能力。图11-3所示为杆件在两端铰支和添加中间绞支情况下发生失稳现象的示意图。图11-3杆件受压失稳示意图求解增加中间铰支后的压杆临界载荷，验证添加中间铰支后的稳定性增强效应。有关的几何参数与和材料参数如表11-2所示。问题分析：对细长杆，可采用二维分析，使用梁单元建模，简化有限元模型。杆的约束情况为，杆长垂直方向3个铰支点位移为0，杆长方向一端固定，另一端承受压力载荷。表11-2几何参数与和材料参数几何参数材料参数杆长200杆截面正方形0.5X0.5杨氏模量30000000（注：本问题中没有给参数定义单位，但在ANSYS系统中不影响分析。）11.4.2前处理1.设定工作目录、项目名称，可使用ANSYS14MechanicalAPDLProductLauncher14.0登录，输入WorkingDirectory和JobName。可根据需要任意填写，但注意不要使用中文。2.定义单元属性。（1）选择MainMenuPreprocessorElementTypeAdd/Edit/Delete命令，在弹出的对话框中单击Add按钮，如图11-4（a）所示。（2）弹出LibraryofElementTypes，在图11-4（b）中选中Beam,2node188，单击按钮OK确认，回到ElementTypes对话框。（a）（b）图11-4定义单元（3）选中前一步定义的单元后，单击Options按钮；弹出BEAM188elementtypeoptions对话框，将第三项K3改为CubicForm，使梁单元沿长度方向为三次曲线，如图11-5所示，单击按钮OK确认，关闭对话框。（4）选择MainMenuPreprocessorSectionsBeamCommonSections命令，弹出BeamTool对话框，在对话框中设置ID为1，选择矩形截面，设置B和H为0.5，如图11-6所示，单击按钮OK确认，关闭对话框。图11-5设置梁单元关键选项图11-6设置梁截面3.定义材料特性。（1）选择MainMenuPreprocessorMaterialPropsMaterialModels命令，弹出DefineMaterialModelBehavior对话框，如图11-7（a）所示。（2）在对话框右栏中选择StructuralLinearElasticIsotropic命令，弹出对话框，在对话框中设置EX为3E+007，如图9.7(b）所示，单击OK按钮确认。关闭弹出的提示PRXY为0对话框，并关闭DefineMaterialModelBehavior对话框。（a）添加材料特性（b）设置材料特性图11-7定义材料特性4.建立有限元模型，采用直接生成网格单元的方法建立有限元模型。（1）选择MainMenuPreprocessorModelingCreateNodesInActiveCS命令，弹出CreateNodesinActiveCoordinateSystem对话框。（2）在对话框中，输入如图11-8（a）所示的数据，单击Apply按钮确认，建立节点1。（3）继续在对话框中，输入如图11-8（b）所示的数据，单击OK按钮确认，建立节点21。（4）选择MainMenuPreprocessorModelingCreateNodesFillbetweenNds命令，弹出实体选择对话框，如图11-8（c）所示，依次选择节点1、节点21，单击OK按钮确认，弹出CreateNodesBetween2Nodes对话框。（5）在弹出的对话框中设置参数如图11-8（d）所示，单击OK按钮，生成均匀分布的节点2~20，如图11-8（e）所示。（a）创建节点1（b）创建节点21（c）选取节点（d）填充节点（e）填充后的节点图11-8直接生成节点（6）选择MainMenuPreprocessorModelingCreateElementsAutoNumberedThruNodes命令，弹出实体选取对话框，如图11-9（a）所示，选择节点1和节点2，单击OK按钮确认，生成网格单元1。（7）选择MainMenuPreprocessorModelingCopyElementsAutoNumbered命令，弹出实体选取对话框，如图11-9（b）所示，单击PickAll按钮，弹出对话框CopyElements（Au