ANSYS 非线性和线性对比分析的一个工程实例
二力杆几何非线性分析
土木工程中，钢模板由于制作不精细或搬运模板过程受到碰撞或者挤压等外力作用常常会造成模板某处凸起，在活、恒载作用下或搬运过程中，该处常常会突然从凸起变成凹进状态。这一现象被称为油罐效应，通常采用桁架的失稳模型进行几何非线性简化分析，因为也称为桁架的经典跳越问题。
分析模型如图1所示，采用LINK1单元构成二力杆，两端完全约束，中间节点加集中力。
图1
荷载与顶点位移理论关系为：
0(sin )(2sin )P EA x x x θθ=--
式中，E 为弹性模量，0A 为杆件初始截面积，x=V/0L ,V 为顶点的竖向位移，0L 为杆长，θ为杆件倾角。分析中取0A =102mm ,E=200Gpa, 0L =100mm, θ=0
6。
分析所用的非线性命令流如下。
!********************** 二力杆几何非线性分析************************
FINI
/CLEAR
/FILENAME,NONLINEAR_ER-LI-GAN
/TITLE,The Analysis of NONLINEAR_ER-LI-GAN /PREP7
/PNUM,LINE,1
/PNUM,KP,1
LO=100
CTA=6
*AFUN,DEG
L1=2*LO*COS(CTA)
H1=LO*SIN(CTA)
AA=10
EM=2E5
ET,1,LINK1
MP,EX,1,EM
R,1,AA
K,1
K,2,0.5*L1,H1
K,3,L1
L,1,2$L,2,3
LESIZE,ALL,,,1
LMESH,ALL
FINI
/SOLU
DK,1,ALL
DK,3,ALL
FK,2,FY,-1200
ANTYPE,0
NLGEOM,1
NSUBST,100
OUTRES,ALL,ALL
ARCLEN,ON
SOLVE
FINI
/POST1
/ESHAPE,1
EPLOT
SET,LAST
PLDISP,1 !绘制变形图
PRRSOL
FINI
/POST26
NSOL,2,2,U,Y,DISPLACEMENT
ABS,3,2
RFORCE,4,1,F,Y,F
/AXLAB,X,DISPLACEMENT
/AXLAB,Y,F
/GRID,1
XVAR,3
PLVAR,4 !绘制节点2位移和节点1的竖向支反力关系
FINI
通用后处理中可得到最后一个荷载步时节点1 和节点3受到的水平方向支座反力分别为-4759.0N和4759.0N，双杆均受拉。同时可得到结构变形图，如图2所示，并且由图可知此时节点2的位移最大为23.001mm，基本与实际情况一致。
图2
由时间—历程后处理分析可知随着时间的增长，荷载的增加，节点2 位移变大，节点1的竖向反力也不断增加，二者之间关系见图3。
由图3可知位移—荷载曲线与理论公式分析基本符合，但随着杆件倾角 的增大，曲线两个顶点附近的误差会越来越大，这与理论公式的假定有关，所以，该类模型只能模拟倾角比较的构件。
图3
以上问题是一个纯粹的几何非线性问题，分析可知在整个结构的变形过程中，尽管结构的位移很大，但是杆件的应力和应变之间一直保持着线性关系，即正比例关系，这也同实际相同，故该分析所得结果比较合理。
二力杆线性分析
分析所用的非线性命令流如下。
!*********************以下为线性分析命令流************************ FINI
/CLEAR
/FILENAME,ER-LI-GAN
/TITLE,The Analysis of er-li-gan
/PREP7
/PNUM,LINE,1
/PNUM,KP,1
LO=100
CTA=6
*AFUN,DEG
L1=2*LO*COS(CTA)
H1=LO*SIN(CTA)
AA=10
EM=2E5
ET,1,LINK1
MP,EX,1,EM
R,1,AA
K,1
K,2,0.5*L1,H1
K,3,L1
L,1,2
L,2,3
LESIZE,ALL,,,1
LMESH,ALL
DK,1,ALL
DK,3,ALL
FK,2,FY,-1200
FINI
/SOLU
DK,1,ALL
DK,3,ALL
FK,2,FY,-1200
ANTYPE,0
SOLVE
FINI
/POST1
/ESHAPE,1
EPLOT
PLDISP,1 !绘制变形图
PRRSOL !列出节点力
FINI
通用后处理中可得到节点 1 和节点3受到的水平方向支座反力分别为5708.6N和-5708.6N。同时可得到结构变形图，如图4所示，并且由图可知此时节点2的位移最大为2.75mm，节点2 还停留在X轴上方，远远小于非线性分析时的位移，同时杆件仍然受压，这与实际情况不符，所以对于这种经典的跳越式二力杆问题应该进行非线性分析才能得到比较合理的解。
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 结果评价和分析
Project1 梁的有限元建模与变形分析 计算分析模型如图1-1 所示, 习题文件名: beam。 NOTE:要求选择不同形状的截面分别进行计算。 梁承受均布载荷：1.0e5 Pa 图1-1梁的计算分析模型 梁截面分别采用以下三种截面（单位：m）： 矩形截面：圆截面：工字形截面： B=0.1, H=0.15 R=0.1 w1=0.1，w2=0.1，w3=0.2, t1=0.0114，t2=0.0114，t3=0.007 1.1进入ANSYS 程序→ANSYSED 6.1 →Interactive →change the working directory into yours →input Initial jobname: beam→Run 1.2设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK 1.3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete… →Add… →select Beam 2 node 188 →OK (back to Element Types window)→Close (the Element Type window) 1.4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural→Linear→Elastic→Isotropic→input EX:2.1e11, PRXY:0.3→OK 1.5定义截面 ANSYS Main Menu: Preprocessor →Sections →Beam →Common Sectns→分别定义矩形截面、圆截面和工字形截面：矩形截面:ID=1,B=0.1，H=0.15 →Apply →圆截面：ID=2,R=0.1 →Apply →工字形截面：ID=3,w1=0.1，w2=0.1，w3=0.2，t1=0.0114，t2=0.0114，t3=0.007→OK
1.ANSYS中混凝土模式预应力模拟的算例 !简支梁实体与预应力钢筋分析 /PREP7 egjx=2e5 !Ey agjx=140 !单根钢绞线面积 ehnt=4e4 !Eh xzxs=1.0e-5 !线胀系数 yjl=200000 !定义预加力 et,1,link8 !定义link8单元 et,2,solid95 !定义solid95单元 r,1,agjx !定义link8单元的面积 r,2 !定义第2种实常数 mp,ex,1,egjx !定义link8单元的弹性模量mp,prxy,1,0.3 !定义link8单元的泊松系数mp,alpx,1,1.0e-5 !定义线膨胀系数 mp,ex,2,ehnt !定义solid95单元的弹性模量mp,prxy,2,0.3 !定义solid95单元的泊松系数blc4, , ,100,200,3000 !定义梁体 /view,1,1,1,1 !定义ISO查看 /ang,1 vplot !绘制梁体 kwpave,6 !工作平面移动到关键点6
wpoff,-30 !工作平面移动-30mm(X) wprot,0,0,90 !工作平面旋转 vsbw,1 !分割梁体 wpoff,0,0,-40 !工作平面移动-40mm(Z) vsbw,2 !分割梁体 wpoff,0,40 !工作平面移动40mm(Y) wprot,0,90 !工作平面旋转 vsbw,all !分割梁体 wpstyl !关闭工作平面显示 nummrg,all,,,,low !整理 numcmp,all !压缩编号 esize,30 !定义网分时边长控制 lsel,s,,,28,38,10 !定义line28和38为新的选择集latt,1,1,1 !定义选择集的属性 lmesh,all !对线划分单元 allsel,all !新的选择集为所有的实体 gplot !绘制所有的实体 vsel,s,,,all !定义所有体为选择集 vatt,2,2,2 !定义选择集的属性 mshape,0,3d !将体划分单元的形状定位HEX mshkey,1 !采用MAPPED划分器 vmesh,all !对体进行划分单元
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。 步骤五、产生矩形 在这一步中，你产生一个代表柱体半横截面积的矩形。
第1章 习题 1.ANSYS软件程序包括几大功能模块？分别有什么作用？ 2.如何启动和退出ANSYS程序？ 3.ANSYS程序有哪几种文件类型？ 4.ANSYS结构有限元分析的基本过程是什么？ 5.两杆平面桁架尺寸及角度如习题图1.1所示，杆件材料的弹性模量为2.1×1011Pa，泊松 比为0.3，截面面积为10cm2，所受集中力载荷F=1000N。试采用二维杆单元LINK1计算集中力位置节点的位移和约束节点的约束反力。 习题图1.1 两杆平面桁架 第2章 习题 1.建立有限元模型有几种方法？ 2.ANSYS程序提供了哪几种坐标系供用户选择？ 3.ANSYS程序中如何平移和旋转工作平面？ 4.试分别采用自底向上的建模方法和自顶向下的建模方法建立如习题图2.1所示的平面图 形，其中没有尺寸标注的图形读者可自行假定，并试着采用布尔运算的拉伸操作将平面图形沿法向拉伸为立体图形。
习题图2.1 平面图形 5.试分别利用布尔运算建立如习题图2.2所示的立体图形，其中没有尺寸标注的图形读者 可自行假定。 习题图2.2 立体图形 6.试对习题图2.3所示的图形进行映射网格划分，并任意控制其网格尺寸，图形尺寸读者 可自行假定。 习题图2.3 映射网格划分
第3章 习题 1.试阐述ANSYS载荷类型及其加载方式。 2.试阐述ANSYS主要求解器类型及其适用范围。 3.如何进行多载荷步的创建，并进行求解？ 4.试建立如习题图3.1所示的矩形梁，并按照图形所示施加约束和载荷，矩形梁尺寸及载 荷位置大小读者可自行假定。 习题图3.1 矩形梁约束与载荷 5.试建立如习题图3.2所示的平面图形，并按照图形所示施加约束和载荷，平面图形的尺 寸及载荷大小读者可自行假定。 习题图3.2 平面图形约束与载荷 第4章 习题
郑重申明：本人能力有限，文中不可避免会有错误，欢迎朋友们批评指正，希望大家相互提高，呵呵，谢谢啦！ 12.8. Sample Rigid Body Dynamic Analysis 刚体动力学分析实例 This sample analysis demonstrates how to model a flexible component in ANSYS and export the flexible body information to a file for use in ADAMS. The example also provides brief instructions on how to perform the rigid body dynamic analysis in ADAMS, and details on how to transfer the loads from ADAMS to ANSYS in order to perform a stress analysis. 该实例演示了如果在ANSYS中制作柔性部件及输出可在ADAMS中使用的柔性体信息文件。同样该例子也提供了有关于在ADAMS中进行动力学分析的简单介绍，和如何将载荷信息从ADAMS转换到ANSYS中进行应力分析的详细介绍。 12.8.1. Problem Description 问题描述 In the linkage assembly shown below, Link3 is a flexible component. Link3 is modeled as a rectangular rod in ANSYS using SOLID45elements. The joints in ADAMS will be attached to interface points (nodes) at the middle of the holes at either end of Link3. These middle points are connected to the cylindrical joint surfaces by a spider web of BEAM4 elements. 联动装置装配如下图所示，连杆3是一个柔性部件，为矩形杆件在ANSYS中采用SOLID45单元构造。ADAMS中连接铰将连接在位于杆两端的孔中心接触节点上。这些节点会通过BEAM4单元构造的蜘蛛网格与圆柱铰表面连接。 Figure 12.5: Linkage Assembly联动装置转配图
!ANSYS命令流学习笔记10-利用APDL在WorkBench中进行非线性屈曲分析 !学习重点： !1、强化非线性屈曲知识 首先了解屈曲问题。在理想化情况下,当F < Fcr时, 结构处于稳定平衡状态，若引入一个小的侧向扰动力，然后卸载, 结构将返回到它的初始位置。当F > Fcr时, 结构处于不稳定平衡状态, 任何扰动力将引起坍塌。当F = Fcr时,结构处于中性平衡状态，把这个力定义为临界载荷。在实际结构中, 几何缺陷的存在或力的扰动将决定载荷路径的方向。在实际结构中, 很难达到临界载荷，因为扰动和非线性行为, 低于临界载荷时结构通常变得不稳定。 要理解非线性屈曲分析，首先要了解特征值屈曲。特征值屈曲分析预测一个理想线弹性结构的理论屈曲强度，缺陷和非线性行为阻止大多数实际结构达到理想的弹性屈曲强度，特征值屈曲一般产生非保守解, 使用时应谨慎。 !理论解，根据Euler公式。其中μ取决于固定方式。 !有限元方法， 已知在特征值屈曲问题： 求解，即可得到临界载荷 而非线性屈曲问题： 其中为结构初始刚度，为有缺陷的结构刚度，为位移矩阵，为载荷矩阵。 非线性屈曲分析时考虑结构平衡受扰动（初始缺陷、载荷扰动）的非线性静力分析，该分析时一直加载到结构极限承载状态的全过程分析，分析中可以综合考虑材料塑性、几何非线性、接触、大变形。非线性屈曲比特征值屈曲更精确，因此推荐用于设计或结构的评价。 !2、熟悉WB中非线性屈曲分析流程 (1) 前处理，施加单元载荷，进行预应力静力分析。 (2) 基于预应力静力分析，指定分析类型为特征值屈曲分析，完成特征值屈曲分析。 (3) 在APDL模块将一阶特征屈曲模态位移乘以适当系数，将此变形后的形状当做非线性分析的初始模型。
巷道开挖过程的有限元模型与力学分析 学院： 班级： 学号： 姓名： 指导老师：XXX老师 XX大学 2009－2010年度上学期
目录 巷道开挖过程的有限元模型与力学分析 (3) 有限元模型概述 (3) 1 建立有限元模型 (4) 1.1 定义工作文件和工作标题 (4) 1.2 定义单元类型、实常数和材料 (4) 1.3 建立几何模型 (5) 2 网格划分 (7) 3 加载与初始地应力模拟 (8) 3.1 设置分析类型 (8) 3.2施加边界条件 (9) 3.3 施加上部面压力： (10) 3.4 施加重力加速度： (10) 3.4 设置加载步骤： (10) 3.5 求解初始地应力： (11) 3.6 保存分析结果： (11) 4 浏览初始地应力的计算结果（后处理） (12) 4.1显示变形形状： (12) 4.2查看节点结果等值线图： (12) 5 开挖巷道求解 (16) 5.1杀死巷道对应的单元： (16) 5.2输入加载步文件： (16) 5.3查看杀死巷道单元后的受力情况： (16) 5.4求解开挖后的有限元模型： (16) 6 浏览查看开挖后的计算结果（后处理） (17) 6.1显示变形形状： (17) 6.2查看节点结果等值线图： (17) 7 应力集中分析 (22) 8 总结体会 (22)
巷道开挖过程的有限元模型与力学分析 某半圆形拱巷道断面，参数见图，其所处地质条件为IV级围岩，上覆盖层厚度为100米后，各材料的力学参数见表。 有限元模型概述 本题采用ANSYS有限元分析软件模拟巷道开挖过程。由于地下巷道属于细长结构物，即巷道的横断面相对于纵向的长度来说很小，可且假定在围着荷载作用下，在其纵向没有位移，只有横向发生位移. 所以，巷道的力学分析可以采用弹性力学理论中的平面应变模型进行，这是一个较复杂的非线性力学问题。采用ANSYS有限元分析软件对巷道开挖进行模拟时，应首先根据地质条件建立合适的地下有限元分析模型，由于巷道对整个地下空间来说是属于“小孔口问题”，巷道周围出现孔口应力集中，并且应力集中区域影响范围约大于1.5倍的孔口尺寸，因此建几何模型时要选择合适的尺寸，其次采用PLANE42单元类型来分析平面应变问题。接下来设定单元尺寸大小划分网格。然后分步求解载荷，第一步求解初始地应力，第二步杀死巷道对应的单元后求解开挖后的载荷，最后分析巷道周围岩石的、位移、应力、应变的变化。
目录 非线性结构分析的定义 (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
关于非线性分析的几点忠告 了解程序的运作方式和结构的表现行为 如果你以前没有使用过某一种特别的非线性特性，在将它用于大的，复杂的模型前，构造一个非常简单的 模型（也就是，仅包含少量单元），以及确保你理解了如何处理这种特性。 通过首先分析一个简化模型，以便使你对结构的特性有一个初步了解。对于非线性静态模型，一个初步的 线性静态分析可以使你知道模型的哪一个区域将首先经历非线性响应，以及在什么载荷范围这些非线性将 开始起作用。对于非线性瞬态分析，一个对梁，质量块及弹簧的初步模拟可以使你用最小的代价对结构的 动态有一个深入了解。在你着手最终的非线性瞬时动态分析前，初步非线性静态，线性瞬时动态，和/或模 态分析同样地可以有助于你理解你结构的非线性动态响应的不同的方面。 阅读和理解程序的输出信息和警告。至少，在你尝试后处理你的结果前，确保你的问题收敛。对于与路程 相关的问题，打印输出的平衡迭代记录在帮助你确定你的结果是有效还是无效方面是特别重的。 简化 尽可能简化最终模型。如果可以将3─D结构表示为2─D平面应力，平面应变或轴对称模型，那么这样做， 如果可以通过对称或反对称表面的使用缩减你的模型尺寸，那么这样做。（然而，如果你的模型非对称加 载，通常你不可以利用反对称来缩减非线性模型的大小。由于大位移，反对称变成不可用的。）如果你可 以忽略某个非线性细节而不影响你模型的关键区域的结果，那么这样做。 只要有可能就依照静态等效载荷模拟瞬时动态加载。 考虑对模型的线性部分建立子结构以降低中间载荷或时间增量及平衡迭代所需要的计算时间。 采用足够的网格密度 考虑到经受塑性变形的区域要求一个合理的积分点密度。每个低阶单元将提供和高阶单元所能提供的一样
2.1 三维实体元SOLID45 2.1.1 单元简介 1．单元形状 SOLID45单元是三维实体单元，该单元在每个节点上有三个自由度（X、Y、Z方向的位移）。 SOLID45的几何形状 2．单元参数输入 对于SOLID45单元而言，其必须或可选的参数输入包括： ●节点：I、J、K、L、M、N、O、P ●自由度：UX, UY, UZ ●材料参数 ?EX（弹性模量）、EY、EZ ?PRXY, PRYZ, PRXZ (or NUXY, NUYZ, NUXZ)（泊松比） ?ALPX, ALPY, ALPZ (or CTEX, CTEY, CTEZ or THSX, THSY, THSZ)（热膨胀系数） ?DENS（密度） ?其它参数 3．结果的输出
SOLID45单元的结果输出 对于三维实体单元SOLID45而言，工程上主要关注其在荷载作用下的变形（挠度）及应力各分量等。对于实体结构的变形，可以直接通过ANSYS中的位移来获得，而应力各分量，包括S：X, Y, Z, XY, YZ, XZ，主应力分量S：1, 2, 3，等效应力S：EQV，都可以在后处理中绘制等值线图，便于观察。 2.1.2 实例分析 考虑一根10m高的混凝土柱，EX＝2e10Pa，Prxy=0.167；Dens=2400；截面尺寸为B*H=0.6m×0.8m，在自由端作用100KN的水平推力，分析其变形和内力。 ANSYS的分析过程： ●进入ANSYS工作界面；改变工作目录Utility Menu>Change Directory… ●进入前处理模块：Main Menu>Preprocessor ●定义单元类型； 点击Element Type > Add/Edit/Delete，在弹出的对话框中添加；
平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end
第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 的屈服极限、切向模量 !定义截面 !壳厚度
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(剪切模量）键入。
/PREP7 /TITLE,Steady-state thermal analysis of pipe junction /UNITS,BIN ! 英制单位；Use U. S. Customary system of units (inches) ! /SHOW, ! Specify graphics driver for interactive run ET,1,90 ! Define 20-node, 3-D thermal solid element MP,DENS,1,.285 ! Density = .285 lbf/in^3 MPTEMP,,70,200,300,400,500 ! Create temperature table MPDATA,KXX,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12 ! 指定与温度相对应的数据材料属性；导热系数；Define conductivity values MPDATA,C,1,,.113,.117,.119,.122,.125 ! Define specific heat values（比热） MPDATA,HF,2,,426/144,405/144,352/144,275/144,221/144 ! Define film coefficient;除144是单位问题，上面的除12也是单元问题 ! Define parameters for model generation RI1=1.3 ! Inside radius of cylindrical tank RO1=1.5 ! Outside radius Z1=2 ! Length RI2=.4 ! Inside radius of pipe RO2=.5 ! Outside pipe radius Z2=2 ! Pipe length CYLIND,RI1,RO1,,Z1,,90 ! 90 degree cylindrical volume for tank WPROTA,0,-90 ! 旋转当前工作的平面；从Y到Z旋转-90度；；Rotate working plane to pipe axis CYLIND,RI2,RO2,,Z2,-90 ! 角度选择在了第四象限；90 degree cylindrical volume for pipe WPSTYL,DEFA ! 重新安排工作平面的设置；另外WPSTYL,STAT to list the status of the working plane；；Return working plane to default setting BOPT,NUMB,OFF ! 关掉布尔操作的数字警告信息；Turn off Boolean numbering warning VOVLAP,1,2 ! 交迭体；Overlap the two cylinders /PNUM,VOLU,1 ! 体编号打开；Turn volume numbers on /VIEW,,-3,-1,1
短圆柱体的热传导过程 问题：一短圆柱体，直径和高度均为1m，现在其上端面施加大小为100℃的均匀温度载荷，圆柱体下端面及侧面的温度均为0℃，试求圆柱体内部的温度场分布（假设圆柱体不与外界发生热交换）。圆柱体材料的热传导系数为30W/(m·℃)。 求解： 第一步：建立工作文件名和工作标题 在ANSYS软件中建立相应的文件夹，并选择Thermal复选框。 第二部：定义单元类型 在单元类型（element type）中选择thermal solid和quad 4node 55，在单元类型选择数字（element type reference number）输入框中输入1，在单元类型选择框里选择Axisymmetric，其余默认即可。 第三步：定义材料性能参数 在材料性能参数对话框中输入圆柱体的导热系数30. 第四步：创建几何模型、划分网格 创建数据点，输入点坐标。在第一个输入框中输入关键点编号1，并输入第一个关键点坐标0、0、0，重复输入第二个、第三个、第四个关键点，相应的坐标分别为2（0.5,0,0）；3（0.5,1,0）；4（0,1,0）。结果如下图1所示：
在模型中创建直线，选择编号为1、2的关键点生成一条直线，在选取2、3生成一条直线，同样选择编号为3、4和编号为4、1的关键点生成另外两条直线。 结果如下图2所示： 之后在plot numbering controls对话框，分别打开KP Keypoint numbers、LINE line numbers、AREA Area numbers，建立直线L1、L2、L3、L4线段。生成几 何模型，如下图所示：
a n s y s学习非线性静态 分析实例 Document serial number【NL89WT-NY98YT-NC8CB-NNUUT-NUT108】
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。