ansys workbench非线性屈曲分析
(2013-08-26 21:26:29)
转载▼
标签：
ansys
很多旋转受压结构必须进行屈曲分析，常规结构屈曲分析软件有nastran、abaqus和ansys，nastran对线性大型模型分析效率较高；abaqus屈曲分析使用较少；ansys使用比较频繁，其快速建模，与CAD软件的良好借口及有限元模型前处理的便捷性（WB界面）很有吸引力，屈曲分析功能较为完善，可以进行线性、非线性和后屈曲分析。
ansys学习资料中介绍较多的是线性屈曲分析。线性屈曲分析在工业实际中预测的值偏高，有的甚至超过实际实验测试值的几十倍，线性分析唯一优势是其分析速度较快。但在实际中其预测值参考价值不大，仅给定结构屈曲失效的上限值。非线性屈曲分析考虑其他因素，包括结构加工缺陷（几何），材料非线性等，因此较为接近实际情况，但计算耗时较长。针对最艰难学习情况归纳总结非线性屈曲分析时技术要点及应注意事项。
对于规则旋转壳，承受外压载荷作用，进行非线性屈曲分析时，必须加上几何缺陷，关键步是添加APDL语句
/prep7
upgeom,0.1,1,1,file,rst
cdwrite,db,file,cdb
/solu
该步引入屈曲模态情况下的几何缺陷，缺陷为屈曲模态变形相对值的0.1倍，该值可以根据实际加工水平等其他条件确定，上述
语句保存在txt文档中，在workbench流程APDL模块调用。
分析详细流程为，static structure模块导入几何，施加载荷和边界条件，分析求解，将linear buckling拖入流程中，共享static structure模块数据，进行线性屈曲模块分析，Mechanial APDL
模块调用屈曲分析结果，并调入（addinput）上面内含几何缺陷命令语句命令的txt文件，更新，将Mechanical结果导入Finite Element modeler模块，更新，此时在缺陷附近的单元节点位置发生改变。在Finite Element modeler重构几何，导入static structure模块，此时可设置结构材料塑性参数，如屈服强度，剪切模量等。进入Mechanical分析中，设置分析选项，主要设置分析步数和稳定能，分析步应设的足够大，便于捕捉屈曲临界载荷值。
最后求解，屈曲载荷位移曲线中有突变点，该处即为屈曲临界载荷值。
以受外压的圆柱壳为例
流程如下
线性静力分析后进行屈曲分析，线性静力分析约束及载荷如下，圆柱壳两端固定约束，中间承受1MPa压力作用
屈曲分析结果如下，载荷乘子1.5155，因此结构线性屈曲载荷为1.5155MPa
将屈曲分析模态位移结果导入Mechanical APDL中，添加upgeom.txt更细结构单元节点位置，将Mechanical APDL中结果导入Finite Element Modeler中，重建模型导入static structure 中，添加边界条件，注意此时压力载荷施加应超过屈曲载荷1.5155MPa，这里取2MPa。分析设置如下图，设置分析时间2000s，打开几何分线性开关，设置载荷步初始200，最小100，最大1000.保证能有效捕捉屈曲点。打开stabilization，开启稳定性分析。
最终求解得到结构载荷位移曲线，曲线中应出现一突变段，如下图所示，在1250s时刻，位移突变，因此考虑几何非线性时，结构屈曲载荷为1.25MPa，小于线性屈曲分析结果。
注意以上分析时，缺陷放大因子取0.2。放大因子取0.1时，结构载荷位移曲线没有突变段，呈现线性关系。具体如何取值应根据实际情况确定。
查看:
2630|回复:
workbench
非线性屈曲分析实例与
问题
[复制链接]
caihaitao008
caihaitao008
caihaitao008
caihaitao008
3
主题
帖子
7
积分
学徒工
积分
7
发消息
电梯直达
发表于2013-9-5 19:15:58 |只看该作者||
WB的功能非常强大，但是自由度不如CLASSICAL。现在就遇到一个问题，WB很多教材里面关于稳定性的都是线性屈曲分析，做非线性屈曲分析的方法不是太清楚。在国
外网站上找到一些教程，但是有一些关键步骤不太清楚，所以拿出来和大家一起讨论讨论
如何用WB做非线性屈曲分析。
现在还搞不太清楚的就是提取结果这个过程。
也就是这个过程，将LINEAR BUCKLING的结果连接到MECHANICAL APDL里面，需要写入一
个INP文件，也就是WB和ANSYS之间的数据链接，不知道有哪位大神能解释里面的操作。
之后的操作就简单了，用FE修一下模型，再连到有大变形的静力分析上，通过迭代计算得到发散的解，就能得到非线性的屈曲结果。
但是现在就卡在了从LINEAR BUCKLING 到MECHANICAL APDL这一步，不知道有没有人以前
做过非线性屈曲分析，或是知道该怎么处理这一步操作的谢谢啦~
ANSYS 非线性分析指南(1) 基本过程 第一章结构静力分析 1. 1 结构分析概述 结构分析的定义： 结构分析是有限元分析方法最常用的一个应用领域。结构这个术语是一个广义的概念，它包括土木工程结构，如桥梁和建筑物；汽车结构，如车身、骨架；海洋结构，如船舶结构；航空结构，如飞机机身、机翼等，同时还包括机械零部件，如活塞传动轴等等。 在ANSYS 产品家族中有七种结构分析的类型，结构分析中计算得出的基 本未知量- 节点自由度，是位移；其他的一些未知量，如应变、应力和反力， 可通过节点位移导出。 七种结构分析的类型分别是： a. 静力分析- 用于求解静力载荷作用下结构的位移和应力等。静力分析 包括线性和非线性分析。而非线性分析涉及塑性、应力刚化、大变形、大应变、超弹性、接触面和蠕变，等。 b. 模态分析- 用于计算结构的固有频率和模态。 c. 谐波分析- 用于确定结构在随时间正弦变化的载荷作用下的响应。 d. 瞬态动力分析- 用于计算结构在随时间任意变化的载荷作用下的响应，并且可计及上述提到的静力分析中所有的非线性性质。 e. 谱分析- 是模态分析的应用拓广，用于计算由于响应谱或PSD 输入 随机振动引起的应力和应变。 f. 屈曲分析- 用于计算屈曲载荷和确定屈曲模态，ANSYS 可进行线性特征值和非线性屈曲分析。 g. 显式动力分析- ANSYS/LS-DYNA可用于计算高度非线性动力学和复 杂的接触问题。 除了前面提到的七种分析类型，还有如下特殊的分析应用： ? 断裂力学 ? 复合材料 ? 疲劳分析
? p-Method 结构分析所用的单元：绝大多数的ANSYS 单元类型可用于结构分析。单元类型从简单的杆单元和梁单元一直到较为复杂的层合壳单元和大应变实体单元 1.2 结构线性静力分析 静力分析的定义： 静力分析计算在固定不变的载荷作用下结构的响应。它不考虑惯性和阻尼的影响，如结构受随时间变化载荷的情况。可是静力分析可以计算那些固定不变的惯性载荷对结构的影响，如重力和离心力；以及那些可以近似为等价静力作用的随时间变化载荷，如通常在许多建筑规范中所定义的等价静力风载和地震载荷。 静力分析中的载荷： 静力分析用于计算由那些不包括惯性和阻尼效应的载荷作用于结构或部件上引起的位移、应力、应变和力。固定不变的载荷和响应是一种假定，即假定载荷和结构的响应随时间的变化非常缓慢，静力分析所施加的载荷包括： ? - 外部施加的作用力和压力 ? - 稳态的惯性力如中力和离心力 ? - 位移载荷 ? - 温度载荷 线性静力分析和非线性静力分析 静力分析既可以是线性的也可以是非线性的。非线性静力分析包括所有的非线性类型：大变形、塑性、蠕变、应力刚化、接触、间隙单元、超弹性单元等，本节主要讨论线性静力分析，非线性静力分析在下一节中介绍。 线性静力分析的求解步骤 1 建模 2 施加载荷和边界条件求解 3 结果评价和分析
a b a q u s压杆屈曲分析 Revised by Petrel at 2021
压杆屈曲分析1.问题描述 在钢结构中，受压杆件一般在其达到极限承载力前就会丧失稳定性，所以失稳是钢结构最为突出的问题。压杆整体失稳形式可以是弯曲、扭转和弯扭。钢构件在轴心压力作用下，弯曲失稳是常见的失稳形式。影响轴心受压构件整体稳定性的主要因素为纵向残余应力、初始弯曲、荷载初偏心及端部约束条件等。实际的轴心受压构件往往会存在上述的一种或多种缺陷，导致构件的稳定承载力降低。 本文利用abaqus对一定截面不同长细比下的H型钢构件进行屈曲分析，通过考虑材料非线性、几何非线性并引入初弯曲，得出构件发生弯曲失稳的极限荷载。通过比较不同长细比下的弯曲失稳的临界荷载得出构件荷载位移曲线，并与《规范》中的构件曲线相比较。钢构件的截面尺寸如图1-1所示。 构件的材料特性：,, 图1-1 2.长细比计算 通过计算截面几何特性，截面绕y轴的回转半径为,长细比取值及杆件长度见表1： 表1 50 60 80 100 120 150 180 （m） 1.92 2.30 3.07 3.84 4.60 5.76 6.90 3.模型分析
ABAQUS非线性屈曲分析的方法有riks法，generalstatics法（加阻尼），或者动力法。非线性屈曲分析采用riks算法实现，可以考虑材料非线性、几何非线性已及初始缺陷的影响。其中，初始缺陷可以通过屈曲模态、振型以及一般节点位移来描述。 利用abaqus进行屈曲分析，一般有两步，首先是特征值屈曲分析，此分析为线性屈曲分析，是在小变形的情况进行的，也即上面提到过的模态，目的是得出临界荷载（一般取一阶模态的eigenvalue乘以所设定的load）。其次，就是后屈曲分析，此步一般定义为非线性，原因在于是在大变形情况进行的，一般采用位移控制加修正的弧长法，可以定义材料非线性，以及几何非线性，加上初始缺陷，所以也称为非线性屈曲分析。此步分析，为了得到极限值，需要得出荷载位移曲线的下降段。缺陷较小的结构初始位移变形较小，在极值点突变，而初始缺陷较大的结构，载荷位移曲线较平滑。 4.建模计算过程 建模计算过程以长细比为50的构件为例，其余构件建模计算过程与之类似。 4.1buckle分析 1在buckle分析中创建part模块，创建的模型为三位可变形壳体单元，截面参数见图1-1，构件长度1.92。如图4-1示 图4-1 2定义材料特性及截面属性并将其赋予单元。材料定义为弹塑性，泊松比0.3，屈服强度，弹性模量；腹板和翼缘板为壳单元，厚度分别为0.008和0,01。材料定义见图4-2
Workbench 屈曲分析 1、基础概念 结构在载荷作用下由于材料弹性性能发生变形，若变形后结构上的载荷保持平衡，这种状态称为弹性平衡。如果结构在平衡状态时，受到扰动而偏离平衡位置，当扰动消除后仍能恢复原来平衡状态，这种平衡状态称为稳定平衡状态，反之，如果受到扰动而偏离平衡位置，即使扰动消除，结构仍不能恢复原来的平衡状态，而结构在新的状态下平衡，则原来的平衡状态就成为不稳定平衡状态。 当结构所受载荷达到某一值时，若增加一微小的增量，则结构平衡状态将发生很大的改变，这种现象叫做结构失稳或结构屈曲。 根据失稳的性质，结构稳定问题可分为以下三类： 第一类失稳是理想化情况，即达到某个载荷时，除结构原来的平衡状态存在外，出现第二个平衡状态，故又叫做平衡分叉失稳，数学上就是求解特征值问题，又叫做特征值屈曲分析。 第二类失稳是结构失稳，变形将大大发展，而不会出现新的变形形式，即平衡状态不发生质变，也叫极顶失稳，结构失稳时，相应载荷叫做极限载荷，理想结构或完善结构不存在，总是存在这样那样的缺陷，大多数问题属于第二类失稳问题。 第三类失稳是当在和达到某值时，结构平衡状态发生一明显跳跃，突然过渡到非临近的另一具有较大位移的平衡状态，称为跳跃失稳，跳跃失稳没有平衡分叉点，也没有极值点，如坦拱、扁壳、二力杆的失稳都属于此类。 结构弹性稳定分析属于第一类失稳对应workbench 的线性特征值分析（Eigenvalue Buckling ），考虑缺陷，非线性影响的第二类结构属于workbench 的非线性特征值分析（Eigenvalue Buckling ），第三类的失稳对应workbench 的Static Structural ，无论前屈曲平衡状态或后屈曲平衡状态均可一次计算求出，即全过程分析。 1.1屈曲分析基础理论 在平衡状态，考虑到轴向力或中面内力对弯曲变形的影响，根据势能驻值原理得到结构平衡方程为 [][](){}{}P U K K G E =+ 式中为结构弹性刚度矩阵，为结构几何刚度矩阵，也称为初应力刚度矩阵，为节点位移向量；为节点载荷向量，上式也为几何非线性分析平衡方程。 为得到随遇平衡状态，应是系统势能的二阶变分为零。即： [][](){}0=+U K K G E δ 因此必有： [][]()0K E =+G K 式中结构弹性刚度矩阵已知，结构外载荷也就是要求得屈曲载荷未知，结构几何刚度矩 阵未知，为了求得该屈曲载荷，假设有一组载荷[]0P ,对应的几何刚度矩阵为[]0G K ,并假定 屈曲时的载荷是[]0P 的λ倍，固有λ[] 0G K =[]G K ,上式可变为 []E K []G K {}U {}P
ansys学习－非线性静态分析实例 问题描述 一个子弹以给定的速度射向壁面。壁面假定是刚性的和无摩擦的。将研究子弹和壁面接触后达80微秒长的现象。目的是确定子弹的整个变形，速度历程，以及最大等效Von Mises应变。求解使用SI单位。 用轴对称单元模拟棒。求解最好能通过单一载荷步实现。在这个载荷步中，将同时施加初始速度和约束。将圆柱体末端的节点Y方向约束住以模拟一固壁面。打开自动时间分步来允许ANSYS 确定时间步长。定义分析结束的时间为8E-5秒，以确保有足够长的时间来扑捉整个变形过程。 问题详细说明 下列材料性质应用于这个问题： EX=117.0E09 (杨氏模量） DENS=8930.0 （密度） NUXY=0.35（泊松比） Yield Strength=400.0OE06（屈服强度） Tangent Modulus （剪切模量） 下列尺寸应用于这个问题： 长=32.4E-3m 直径=6.4E-3m 对于这个问题的初始速度是227.0。 图1铜圆柱体图解 求解步骤： 步骤一：设置分析标题 1、选择菜单路径：Utility Menn>File>ChangeTitle。
2、键入文字“Coppery Cylinder Impacting a Rigid Wall” 3、单击OK。 步骤二：定义单元类型 1、选择菜单路径Mail Menu>Preprocessor>Element Type>All/Edit/Delete。 2、单击Add。Library of Element Types(单元类型库）对话框出现。 3、在靠近左边的列表中，单击“Visio Solid”仅一次。 4、选靠近右边的列表中，单击“4node Plas 106”仅一次。 5、单击OK。Library of Element Types 对话框关闭。 6、单击Options (选项）。VISCO106 element type Options(visco106单元类型选项）对话框出现。 7、在关于element behavior（单元特性）的卷动柜中，卷动到“Axisymmetric” 且选中它。 8、单击OK。 9、单击Element Types (单元类型）对话框中的Close。 步骤三：定义材料性质 1、选择菜单路径Main Menu>Preprocessor>Material Props>-Constant-Isotropic. Isotropic Matersal Properties (各向同性材料性质）对话框出现。 2、单击OK来指定材料号为1。另一个I sotropic Material Properties对话框出现。 3、对杨氏模量（EX）键入117.0E09 4、对密度（DENS）键入8930。 5、对泊松比（NUXY）键入0.35。 6、单击OK。 步骤四：定义双线性各向同性强化数据表（BISO） 1、选择菜单路径Main Menu>Preprocessor>Matersal Props>Data Tables> Define/Activate . Define/Activate Data Table(定义数据表）对话柜出现。 2、在关于type of data table(数据表类型）的卷动框中，卷动到“Bilin isotr BISO”且选中它。 3、对material reference number(材料参考号）健入1。 4、对number of temperatures(温度数）键入1和单击OK。 5、选择菜单路径Main Menu>Preprocessor>Material Props>Data Tables>Edit Active. Data Table BISO对话框出现。 6、对YLD Strs(屈服应力）键入400.0e06。 7、对 Tang Mod(剪切模量）键入100.0e06。 8、选择File>Apply & Quit。 9、选择菜单路径Main Menu>Preprosessor>Material Porps>Data Tables>Graph. Graph Data Tables(图形表示数据表）对话框出现。 10、单击OK接受绘制BISO表的缺省。一个BISO表的标绘图出现在ANSYS图形窗口中。 11、在ANSYS TooLbar上单击SAVE_DB。 步骤五、产生矩形 在这一步中，你产生一个代表柱体半横截面积的矩形。
一.问题描述 在钢结构中，受压杆件一般在其达到极限承载力前就会丧失稳定性，所以失稳是钢结构最为突出的问题。压杆整体失稳形式可以是弯曲、扭转和弯扭。钢构件在轴心压力作用下，弯曲失稳是常见的失稳形式。而影响轴心受压构件整体稳定性的主要因素为纵向残余应力、初始弯曲、荷载初偏心及端部约束条件等。实际的轴心受压构件往往会存在上述的一种或多种缺陷，导致构件的稳定承载力降低。 本文主要针对任意轴对称的圆形钢管截面，利用ABAQUS有限元非线性分析软件，对其在轴心受压情况下进行特征值屈曲分析和静态及动态的非线性屈曲分析（考虑材料弹塑性和初始缺陷的影响）。通过考虑材料非线性、几何非线性并引入初弯曲，得出构件发生弯曲失稳的极限荷载,并且由弯曲失稳的临界荷载得出的构件荷载位移曲线。同时再进行非线性分析时，需要施加初始扰动，以帮助非线性分析时失稳，可以通过特征值屈曲分析得到的初始弯曲模态来定义初始缺陷；最后由可以将特征值屈曲分析得到的临界荷载作为非线性屈曲分析时所施加荷载的参考。 二.结构模型 用ABAQUS中的壳单元建立轴心受压模型，采用SI国际单位制（m）。 1.构件的材料特性： E= 2.0×1011N m2,μ=0.3, f y=2.35×
108N m2，ρ=7800kg m3，钢管半径：60mm，厚度：3mm，长度：2.5m。 2.钢管的截面尺寸及钢管受到的约束和荷载施加的模型图如图2-1及图2-2所示。 图2-1 图2-2 三.建模步骤（Buckle分析） （1）创建部件 在创建part模块中命名构件的名字为gang guan，创建的模型为三维可变形壳体单元，如图3-1所示。截面参数见图2-1，构件长度2.5m。 图3-1
ansys workbench非线性屈曲分析 (2013-08-26 21:26:29) 转载▼ 标签： ansys 很多旋转受压结构必须进行屈曲分析，常规结构屈曲分析软件有nastran、abaqus和ansys，nastran对线性大型模型分析效率较高；abaqus屈曲分析使用较少；ansys使用比较频繁，其快速建模，与CAD软件的良好借口及有限元模型前处理的便捷性（WB界面）很有吸引力，屈曲分析功能较为完善，可以进行线性、非线性和后屈曲分析。 ansys学习资料中介绍较多的是线性屈曲分析。线性屈曲分析在工业实际中预测的值偏高，有的甚至超过实际实验测试值的几十倍，线性分析唯一优势是其分析速度较快。但在实际中其预测值参考价值不大，仅给定结构屈曲失效的上限值。非线性屈曲分析考虑其他因素，包括结构加工缺陷（几何），材料非线性等，因此较为接近实际情况，但计算耗时较长。针对最艰难学习情况归纳总结非线性屈曲分析时技术要点及应注意事项。 对于规则旋转壳，承受外压载荷作用，进行非线性屈曲分析时，必须加上几何缺陷，关键步是添加APDL语句 /prep7 upgeom,0.1,1,1,file,rst cdwrite,db,file,cdb /solu 该步引入屈曲模态情况下的几何缺陷，缺陷为屈曲模态变形相对值的0.1倍，该值可以根据实际加工水平等其他条件确定，上述 语句保存在txt文档中，在workbench流程APDL模块调用。 分析详细流程为，static structure模块导入几何，施加载荷和边界条件，分析求解，将linear buckling拖入流程中，共享static structure模块数据，进行线性屈曲模块分析，Mechanial APDL
| w w w .3d s .c o m | ? D a s s a u l t S y s t èm e s | Buckling, Postbuckling, and Collapse Analysis with Abaqus | w w w .3d s .c o m | ? D a s s a u l t S y s t èm e s | Day 1 ?Lecture 1Basic Concepts and Overview ?Workshop 1Buckling and Postbuckling Analyses of a Crane Structure ?Lecture 2 Finite Element Formulation ?Lecture 3Finite Element Implementation in Abaqus ?Lecture 4Eigenvalue Buckling Analysis ?Workshop 2Eigenvalue Buckling of a Ring Subjected to External Pressure ?Workshop 3 Elastic Buckling of Ring-Supported Cylindrical Shell under Hydrostatic Pressure
| w w w .3d s .c o m | ? D a s s a u l t S y s t èm e s | Buckling, Postbuckling, and Collapse Analysis with Abaqus Day 2 ?Lecture 5 Regular and Damped Static Solution Procedures for Postbuckling Analyses ?Workshop 4Nonlinear Buckling of Ring-Supported Cylindrical Shell under Hydrostatic Pressure ?Workshop 5Static Buckling Analysis of a Circular Arch ?Lecture 6Modified Riks Static Solution Procedure for Postbuckling Analyses ?Workshop 5Static Buckling Analysis of a Circular Arch (continued)?Lecture 7Dynamic Analysis Solution Procedures for Postbuckling Analyses ?Workshop 5Static Buckling Analysis of a Circular Arch (continued)?Workshop 6Tube Crush Dynamic Analysis ?Lecture 8Putting It All Together… ?Workshop 7Capstone Workshop: Lee’s Frame Buckling Problem ?Workshop 8 Buckling and Postbuckling Analyses of a Stiffened Panel | w w w .3d s .c o m | ? D a s s a u l t S y s t èm e s | Legal Notices The Abaqus Software described in this documentation is available only under license from Dassault Systèmes and its subsidiary and may be used or reproduced only in accordance with the terms of such license. This documentation and the software described in this documentation are subject to change without prior notice. Dassault Systèmes and its subsidiaries shall not be responsible for the consequences of any errors or omissions that may appear in this documentation. No part of this documentation may be reproduced or distributed in any form without prior written permission of Dassault Systèmes or its subsidiary.? Dassault Systèmes, 2011. Printed in the United States of America Abaqus, the 3DS logo, SIMULIA and CATIA are trademarks or registered trademarks of Dassault Systèmes or its subsidiaries in the US and/or other countries. Other company, product, and service names may be trademarks or service marks of their respective owners. For additional information concerning trademarks, copyrights, and licenses, see the Legal Notices in the Abaqus 6.11 Release Notes and the notices at: https://www.doczj.com/doc/0c7540381.html,/products/products_legal.html.
!ANSYS命令流学习笔记10-利用APDL在WorkBench中进行非线性屈曲分析 !学习重点： !1、强化非线性屈曲知识 首先了解屈曲问题。在理想化情况下,当F < Fcr时, 结构处于稳定平衡状态，若引入一个小的侧向扰动力，然后卸载, 结构将返回到它的初始位置。当F > Fcr时, 结构处于不稳定平衡状态, 任何扰动力将引起坍塌。当F = Fcr时,结构处于中性平衡状态，把这个力定义为临界载荷。在实际结构中, 几何缺陷的存在或力的扰动将决定载荷路径的方向。在实际结构中, 很难达到临界载荷，因为扰动和非线性行为, 低于临界载荷时结构通常变得不稳定。 要理解非线性屈曲分析，首先要了解特征值屈曲。特征值屈曲分析预测一个理想线弹性结构的理论屈曲强度，缺陷和非线性行为阻止大多数实际结构达到理想的弹性屈曲强度，特征值屈曲一般产生非保守解, 使用时应谨慎。 !理论解，根据Euler公式。其中μ取决于固定方式。 !有限元方法， 已知在特征值屈曲问题： 求解，即可得到临界载荷 而非线性屈曲问题： 其中为结构初始刚度，为有缺陷的结构刚度，为位移矩阵，为载荷矩阵。 非线性屈曲分析时考虑结构平衡受扰动（初始缺陷、载荷扰动）的非线性静力分析，该分析时一直加载到结构极限承载状态的全过程分析，分析中可以综合考虑材料塑性、几何非线性、接触、大变形。非线性屈曲比特征值屈曲更精确，因此推荐用于设计或结构的评价。 !2、熟悉WB中非线性屈曲分析流程 (1) 前处理，施加单元载荷，进行预应力静力分析。 (2) 基于预应力静力分析，指定分析类型为特征值屈曲分析，完成特征值屈曲分析。 (3) 在APDL模块将一阶特征屈曲模态位移乘以适当系数，将此变形后的形状当做非线性分析的初始模型。
屈曲分析的过程说明： 屈曲分析是一种用于确定结构开始变得不稳定时的临介荷载和屈 曲结构发生屈曲响应时的模态形状的技术。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
目录 非线性结构分析的定义 (1) 非线性行为的原因 (1) 非线性分析的重要信息 (3) 非线性分析中使用的命令 (8) 非线性分析步骤综述 (8) 第一步：建模 (9) 第二步：加载且得到解 (9) 第三步：考察结果 (16) 非线性分析例题（GUI方法） (20) 第一步：设置分析标题 (21) 第二步：定义单元类型 (21) 第三步：定义材料性质 (22) 第四步：定义双线性各向同性强化数据表 (22) 第五步：产生矩形 (22) 1
第六步：设置单元尺寸 (23) 第七步：划分网格 (23) 第八步：定义分析类型和选项 (23) 第九步：定义初始速度 (24) 第十步：施加约束 (24) 第十一步：设置载荷步选项 (24) 第十二步：求解 (25) 第十三步：确定柱体的应变 (25) 第十四步：画等值线 (26) 第十五步：用Post26定义变量 (26) 第十六步：计算随时间变化的速度 (26) 非线性分析例题（命令流方法） (27) 非线性结构分析 非线性结构的定义 在日常生活中,会经常遇到结构非线性。例如，无论何时用钉书针钉书，金 2
属钉书钉将永久地弯曲成一个不同的形状。（看图1─1（a））如果你在一个木架上放置重物，随着时间的迁移它将越来越下垂。（看图1─1（b））。当在 汽车或卡车上装货时，它的轮胎和下面路面间接触将随货物重量的啬而变化。（看图1─1（c））如果将上面例子所载荷变形曲线画出来,你将发现它们都显示了非线性结构的基本特征--变化的结构刚性. 图1─1 非线性结构行为的普通例子 3
非线性行为的原因 引起结构非线性的原因很多，它可以被分成三种主要类型： 状态变化（包括接触） 许多普通结构的表现出一种与状态相关的非线性行为,例如，一根只能拉伸的电缆可能是松散的,也可能是绷紧的。轴承套可能是接触的,也可能是不接触的, 冻土可能是冻结的,也可能是融化的。这些系统的刚度由于系统状态的改变在不同的值之间突然变化。状态改变也许和载荷直接有关（如在电缆情况中），也可能由某种外部原因引起（如在冻土中的紊乱热力学条件）。ANSYS程序中单元的激活与杀死选项用来给这种状态的变化建模。 接触是一种很普遍的非线性行为，接触是状态变化非线性类型形中一个特殊而重要的子集。 几何非线性 如果结构经受大变形，它变化的几何形状可能会引起结构的非线性地响应。一个例的垂向刚性）。随着垂向载荷的增加，杆不断弯曲以致于动力臂明显地减少，导致杆端显示出在较高载荷下不断增长的刚性。 4
第二章 ABAQUS 基本使用方法 [2](pp15)快捷键：Ctrl+Alt+左键来缩放模型；Ctrl+Alt+中键来平移模型；Ctrl+Alt+右键来旋转模型。 ②(pp16)ABAQUS/CAE 不会自动保存模型数据，用户应当每隔一段时间自己保存模型以避免意外 丢失。 [3](pp17)平面应力问题的截面属性类型是Solid（实心体）而不是Shell（壳）。ABAQUS/CAE 推荐的建模方法是把整个数值模型（如材料、边界条件、载荷等）都直接定义在几 何模型上。 载荷类型Pressure 的含义是单位面积上的力，正值表示压力，负值表示拉力。 [4](pp22)对于应力集中问题，使用二次单元可以提高应力结果的精度。 [5](pp23)Dismiss 和Cancel 按钮的作用都是关闭当前对话框，其区别在于：前者出现在包含只读数 据的对话框中；后者出现在允许作出修改的对话框中，点击Cancel 按钮可关闭对话框，而不保存 所修改的内容。 [6](pp26)每个模型中只能有一个装配件，它是由一个或多个实体组成的，所谓的“实体”（instance） 是部件（part）在装配件中的一种映射，一个部件可以对应多个实体。材料和截面属性定义在部件 上，相互作用（interaction）、边界条件、载荷等定义在实体上，网格可以定义在部件上或实体上， 对求解过程和输出结果的控制参数定义在整个模型上。 [7](pp26) ABAQUS/CAE 中的部件有两种：几何部件（native part）和网格部件（orphan mesh part）。 创建几何部件有两种方法：（1）使用Part 功能模块中的拉伸、旋转、扫掠、倒角和放样等特征来直 接创建几何部件。（2）导入已有的CAD 模型文件，方法是：点击主菜单 File→Import→Part。网 格部件不包含特征，只包含节点、单元、 面、集合的信息。创建网格部件有三种方法：（1）导入 ODB 文件中的网格。（2）导入INP 文件中的网格。（3）把几何部件转化为网格部件，方法是：进 入Mesh 功能模块，点击主菜单Mesh→Create Mesh Part。 [8](pp31)初始分析步只有一个，名称是initial，它不能被编辑、重命名、替换、复制或删除。在初 始分析步之后，需要创建一个或多个后续分析步，主要有两大类：（1）通用分析步（general analysis step）可以用于线性或非线性分析。常用的通用分析步包含以下类型：—Static, General: ABAQUS/Standard 静力分析 —Dynamics, Implicit: ABAQUS/Standard 隐式动力分析 —Dynamics, Explicit: ABAQUS/ Explicit 显式动态分析
关于非线性分析的几点忠告 了解程序的运作方式和结构的表现行为 如果你以前没有使用过某一种特别的非线性特性，在将它用于大的，复杂的模型前，构造一个非常简单的 模型（也就是，仅包含少量单元），以及确保你理解了如何处理这种特性。 通过首先分析一个简化模型，以便使你对结构的特性有一个初步了解。对于非线性静态模型，一个初步的 线性静态分析可以使你知道模型的哪一个区域将首先经历非线性响应，以及在什么载荷范围这些非线性将 开始起作用。对于非线性瞬态分析，一个对梁，质量块及弹簧的初步模拟可以使你用最小的代价对结构的 动态有一个深入了解。在你着手最终的非线性瞬时动态分析前，初步非线性静态，线性瞬时动态，和/或模 态分析同样地可以有助于你理解你结构的非线性动态响应的不同的方面。 阅读和理解程序的输出信息和警告。至少，在你尝试后处理你的结果前，确保你的问题收敛。对于与路程 相关的问题，打印输出的平衡迭代记录在帮助你确定你的结果是有效还是无效方面是特别重的。 简化 尽可能简化最终模型。如果可以将3─D结构表示为2─D平面应力，平面应变或轴对称模型，那么这样做， 如果可以通过对称或反对称表面的使用缩减你的模型尺寸，那么这样做。（然而，如果你的模型非对称加 载，通常你不可以利用反对称来缩减非线性模型的大小。由于大位移，反对称变成不可用的。）如果你可 以忽略某个非线性细节而不影响你模型的关键区域的结果，那么这样做。 只要有可能就依照静态等效载荷模拟瞬时动态加载。 考虑对模型的线性部分建立子结构以降低中间载荷或时间增量及平衡迭代所需要的计算时间。 采用足够的网格密度 考虑到经受塑性变形的区域要求一个合理的积分点密度。每个低阶单元将提供和高阶单元所能提供的一样
Workbench （1）首先进行线性屈曲分析，得到屈曲的特征值和屈曲模态。实现方式如下： （2）添加Mechanical APDL模块 右键单击Analysis，输入模型缺陷文件：
/prep7 upgeom,0.1,1,1,file,rst cdwrite,db,file,cdb /solu UPGEOM, FACTOR, LSTEP, SBSTEP, Fname, Ext FACTOR: Multiplier for displacements being added to coordinates. The value 1.0 will add the full value of the displacements to the geometry of the finite element model. Defaults to 1.0. LSTEP: Load step number of data to be imported. Defaults to the last load step. SBSTEP: sub step number of data to be imported. Defaults to the last substep. Fname: File name and directory path (248 characters maximum, including the characters needed for the directory path). An unspecified directory path defaults to the working directory; in this case, you can use all 248 characters for the file name. The field must be input (no default). Ext:Filename extension (8 character maximum).The extension must be an RST extension. （3）添加Finite Element Modeler模块 （4）重新导入新的Static Structual模块以进行非线性屈曲分析，此时需重新建立模型的接触关系、边界条件、荷载。 本模块分析时需打开大变形（large deflection），
第25例非线性分析综合应用实例----钢板卷制成圆筒 本例介绍了综合利用ANSYS非线性分析功能模拟将钢板卷制成圆筒的方法和步骤。25.1问题描述 将钢板卷制成圆筒一般要使用卷板机。图25-1所示为对称式三辊卷板机， 该机器将钢板卷制成圆筒时分为三个步骤：首先，上辊下降使钢板发生挠曲，钢板挠曲线的最低点首先发生屈服；然后，下辊转动驱动钢板向前移动，使钢板各点发生同样的屈服形成圆筒；最后，圆筒卷制完成，上辊上升卸下筒体。 图25-1对称式三辊卷板机 用ANSYS模拟将钢板卷制成圆筒，相应地也分为三个步骤。由于第二个步骤需要模拟上、下辊转动，而ANSYS的SOLIDn单元不支持大转动，位移边界条件不能施加大的转动角度，所以上、下辊需要用壳单元建立有限元模型。上、下辊与钢板的作用需要用接触模拟，钢板卷制成圆筒材料发生屈服，产生大变形, 所以钢板卷制成圆筒包括状态非线性、材料非线性和结构非线性三种非线性。 用ANSYS模拟将钢板卷制成圆筒，计算结果可以得到圆筒直径与上辊下压量的关系，上、下辊受力大小，上、下辊的变形，下辊驱动力矩及卸载回弹等重
25.2 命令流 /CLEAR /FILNAM, EXAMPLE25 /CONFIG, NRES, 2000 /PREP7 /PNUM, VOLU, ON ET, 1, SHELL181 ET, 2, SOLID186 MP, EX, 1, 2E11 MP, DENS, 1, 7800 MP, NUXY, 1, 0.3 MP, EX, 2, 2E11 MP, DENS, 2, 7800 MP, NUXY, 2, 0.3 TB, BKIN, 2, 1 TBTEMP, 0 TBDATA,, 240E6, 0 SECTYPE, 1, SHELL SECDATA, 0.02 CYLIND, 0.38/2, 0, 0.2, 1.7, 0, 360 要数据。因为分析过程复杂，步骤较多，所以本例只采用命令流法执行命令。 !清除数据库，新建文件 ! 指定任务名为?EXAMPLE25 “ !设置最大子步数 !前处理 !进入前处理器 !打开体号 !选择单元类型，壳单元用于划分上、下辊 !实体单元用于划分钢板 !定义材料模型 1 的弹性模量 ! 定义材料模型 1 的密度 !定义材料模型 1 的泊松比 !定义材料模型 2 的弹性模量 ! 定义材料模型 2 的密度 !定义材料模型 2 的泊松比 ! 定义材料模型 2 的屈服极限、切向模量 !定义截面 !壳厚度
压杆屈曲分析 1.问题描述 在钢结构中，受压杆件一般在其达到极限承载力前就会丧失稳定性，所以失稳是钢结构最为突出的问题。压杆整体失稳形式可以是弯曲、扭转和弯扭。钢构件在轴心压力作用下，弯曲失稳是常见的失稳形式。影响轴心受压构件整体稳定性的主要因素为纵向残余应力、初始弯曲、荷载初偏心及端部约束条件等。实际 2 压杆截面尺寸（单位：m） 图1-1 2.长细比计算 通过计算截面几何特性，截面绕y轴的回转半径为i y=0.0384m ,长细比取
值及杆件长度见表1： 表1 3.模型分析 ABAQUS非线性屈曲分析的方法有riks法，general statics法（加阻尼），或者动力法。非线性屈曲分析采用riks算法实现，可以考虑材料非线性、几何非线性已及初始缺陷的影响。其中，初始缺陷可以通过屈曲模态、振型以及一般节点位移来描述。 利用abaqus进行屈曲分析，一般有两步，首先是特征值屈曲分析，此分析为线性屈曲分析，是在小变形的情况进行的，也即上面提到过的模态，目的是得出临界荷载（一般取一阶模态的eigenvalue乘以所设定的load）。其次，就是后屈曲分析，此步一般定义为非线性，原因在于是在大变形情况进行的，一般采用位移控制加修正的弧长法，可以定义材料非线性，以及几何非线性，加上初始缺陷，所以也称为非线性屈曲分析。此步分析，为了得到极限值，需要得出荷载位移曲线的下降段。缺陷较小的结构初始位移变形较小，在极值点突变，而初始缺陷较大的结构，载荷位移曲线较平滑。 4.建模计算过程 建模计算过程以长细比为50的构件为例，其余构件建模计算过程与之类似。 4.1 buckle分析 1 在buckle分析中创建part模块，创建的模型为三位可变形壳体单元，截面参数见图1-1，构件长度1.92。如图4-1示
a n s y s学习非线性静态分 析实例 Newly compiled on November 23, 2020
ansys学习－非线性静态分析实例 问题描述 一个子弹以给定的速度射向壁面。壁面假定是刚性的和无摩擦的。将研究子弹和壁面接触后达80微秒长的现象。目的是确定子弹的整个变形，速度历程，以及最大等效Von Mises应变。求解使用SI单位。 用轴对称单元模拟棒。求解最好能通过单一载荷步实现。在这个载荷步中，将同时施加初始速度和约束。将圆柱体末端的节点Y方向约束住以模拟一固壁面。打开自动时间分步来允许ANSYS确定时间步长。定义分析结束的时间为8E-5秒，以确保有足够长的时间来扑捉整个变形过程。 问题详细说明 下列材料性质应用于这个问题： EX= (杨氏模量） DENS= （密度） NUXY=（泊松比） Yield Strength=（屈服强度） Tangent Modulus （剪切模量） 下列尺寸应用于这个问题： 长=-3m 直径=-3m 对于这个问题的初始速度是。 图1铜圆柱体图解 求解步骤： 步骤一：设置分析标题 1、选择菜单路径：Utility Menn>File>ChangeTitle。 2、键入文字“Coppery Cylinder Impacting a Rigid Wall” 3、单击OK。 步骤二：定义单元类型
1、选择菜单路径Mail Menu>Preprocessor>Element Type>All/Edit/Delete。 2、单击Add。Library of Element Types(单元类型库）对话框出现。 3、在靠近左边的列表中，单击“Visio Solid”仅一次。 4、选靠近右边的列表中，单击“4node Plas 106”仅一次。 5、单击OK。Library of Element Types 对话框关闭。 6、单击Options (选项）。VISCO106 element type Options(visco106单元类型选项）对话框出现。 7、在关于element behavior（单元特性）的卷动柜中，卷动到“Axisymmetric” 且选中它。 8、单击OK。 9、单击Element Types (单元类型）对话框中的Close。 步骤三：定义材料性质 1、选择菜单路径Main Menu>Preprocessor>Material Props>-Constant-Isotropic. Isotropic Matersal Properties (各向同性材料性质）对话框出现。 2、单击OK来指定材料号为1。另一个I sotropic Material Properties对话框出现。 3、对杨氏模量（EX）键入 4、对密度（DENS）键入8930。 5、对泊松比（NUXY）键入。 6、单击OK。 步骤四：定义双线性各向同性强化数据表（BISO） 1、选择菜单路径Main Menu>Preprocessor>Matersal Props>Data Tables> Define/Activate . Define/Activate Data Table(定义数据表）对话柜出现。 2、在关于type of data table(数据表类型）的卷动框中，卷动到“Bilin isotr BISO”且选中它。 3、对material reference number(材料参考号）健入1。 4、对number of temperatures(温度数）键入1和单击OK。 5、选择菜单路径Main Menu>Preprocessor>Material Props>Data Tables>Edit Active. Data Table BISO对话框出现。 6、对YLD Strs(屈服应力）键入。 7、对 Tang Mod(剪切模量）键入。