ansys经验
Ansys之非线性收敛准则
ansys计算非线性时会绘出收敛图，其中横坐标是cumulative iteration number 纵坐标是absolute convergence norm。他们分别是累积迭代次数和绝对收敛范数，用来判断非线性分析是否收敛。
ansys在每荷载步的迭代中计算非线性的收敛判别准则和计算残差。其中计算残差是所有单元内力的范数，只有当残差小于准则时，非线性叠代才算收敛。ansys的位移收敛是基于力的收敛的,以力为基础的收敛提供了收敛量的绝对值，而以位移为基础的收敛仅提供表现收敛的相对量度。一般不单独使用位移收敛准则,否则会产生一定偏差,有些情况会造成假收敛.(ansys 非线性分析指南--基本过程Page.6) 。因此ansys官方建议用户尽量以力为基础（或力矩）的收敛误差，如果需要也可以增加以位移为基础的收敛检查。ANSYS缺省是用L2范数控制收敛。其它还有L1范数和L0范数，可用CNVTOL命令设置。在计算中L2值不断变化，若L2 由于ANSYS缺省的criterion计算是你全部变量的平方和开平方（SRSS）*valuse(你设置的值)，所以crition也有小小变化。如有需要，也可自己指定crition为某一常数，CNVTOL,F,10000,0.0001,0 就指定力的收敛控制值为10000*0.0001=1。 另外，非线性计算中用到的一个开关是SOLCONTROL 如关闭SOLCONTROL 选项，那么软件默认收敛准则：力或弯矩的收敛容差是0.001，而不考虑位移的收敛容差；如果打开SOLCONTROL 选项，同样的默认收敛准则：力或弯矩的收敛容差是0.005，而位移收敛容差是0.05。 非线性收敛非常麻烦，与网格精度、边界条件、荷载步等一系列因素有关，单元的特点对收敛的影响很大，单元的性态不好收敛则困难些；合理的步长可以使求解在真解周围不至于振荡，步长过小，计算量太大，步长过大，会由于过大的荷载步造成不收敛。网格密度适当有助于收敛，网格太密计算量太大，当然太稀计算结果会有较大的误差。究竟多少往往要针对问题进行多次试算。 如果不收敛，可以考虑一下方法改进 1.放松非线性收敛准则。 (CNVTOL #Sets convergence values for nonlinear analyses). 2.增加荷载步数。 (NSUBST #Specifies the number of substeps to be taken this load step) 3.增加每次计算的迭代次数(默认的25次) (NEQIT #Maximum number of equilibrium iterations allowed each substep) 4 重新划分单元试试，后续会得到不同的答案。 就我的经验,前两种方法效果好一些,后一种方法不一定奏效。但应注意到，放大收敛准则其实是在降低计算精度的条件下得到近似解，并且放大的收敛准则是否与实际相符或有实际意义应仔细考虑。 CNVTOL, Lab, VALUE, TOLER, NORM, MINREF ANSYS中,非线性收敛准则主要有力的收敛，位移的收敛，弯矩的收敛和转角的收敛。一般用力的控制加载时，可以使用残余力的2-范数控制收敛；而位移控制加载时，最好用位移的范数控制收敛。 When SOLCONTROL,ON, TOLER Defaults to 0.005 (0.5%) for force and moment, and 0.05 (5%) for displacement when rotational DOFs are not present. When SOLCONTROL,OFF, defaults to 0.001 (0.1%) for force and moment. 收敛精度一般可放宽至5%，以提高收敛速度。 加快收敛的方法有一下几种： 1可以增大荷载子步数,nsubst,nsbstp,nsbmn,carry 2修改收敛准则,cnvtol,lab,value,toler,norm,minref 3打开优化的非线性默认求解设置和某些强化的内部求解算法， solcontrol,key1,key2,key3,vtol（一般情况下，默认是打开的） 4重新划分网格,网格的单元不宜太大或太小, 一般在5～10厘米左右 5 检查模型的正确性 下面计算收敛过程图中的各个曲线的具体含义是什么？ 非线性计算是一个迭代计算的过程，曲线表示两次迭代之间的误差，图中分别表示力和位移 在迭代过程中的每次迭代之间的误差 关于ansys中收敛准则(cnvtol)理解 ansys中依据缺省的收敛准则，程序将对不平衡力SRSS与VALUE*TOLER的值进行比较；而VALUE的缺省值是在SRSS和MINREF中取较大值。现假如TOLER的缺省值是0.1的话， 这个准则是不是可以理解成后一次的SRSS是前一次的SRSS的01倍就收敛啦？ 请指点 我是这样理解的例如下面的命令流: cnvtol,f,5000,0.0005,0 cnvtol,u,10,0.001,2 如果不平衡力（独立的检查每一个自由度）小于等于5000*0.0005（也就是2.5），并且如果 位移的变化小于等于10*0.001时，认为子步是收敛的。 ANSYS中收敛准则，程序默认力与位移共同控制，并且收敛的控制系数好像是0.001。这样的收敛精度一般很难使塑性分析收敛，对于一般的塑性分析收敛问题，前几个荷载步（弹性阶段）用力与位移共同控制，进入塑性后用力控制或位移控制，也可以先用力后用位移控制（位移 控制比较容易收敛），至于控制系数取多少，自己根据需要逐步放大直至收敛！也有人建议最后 用能量来控制收敛， convergence value 是收敛值，convergence norm是收敛准则。ansys 可以用cnvtol命令,如：cnvtol,f,10000,0.00001,2,,其中f是指采用力结果，10000是收敛绝对值，0.00001是收敛系数，2是收敛2范数。 收敛准则应该是指选取那种结果进行收敛判定，通常有三种选择，分别是力（f），位移（u）、和能量。当然这三种形式可以单独使用也可以联合使用。收敛准则的另一层意思应该是选取什么范数形式（1、2、3范数）。一般结构通常都选 取2范数格式。 而收敛值只是收敛准则中的一部分，如cnvtol命令中的收敛绝对值与收敛系数的乘积就应该是你所指的收敛值（convergence value）。 ansys 使用收敛准则有L1，L2，L~~（无穷大）三个收敛准则。 在工程中，一般使用收敛容差（0.05）就可以拉。 建议使用位移收敛准则( cnvtol,u,0.05,,, )与力收敛准则( cnvtol,f，0.05,,, )。因为仅仅只使用一个收敛准则，会存在较大的误差。 假如你只能是使用一个收敛准则，建议你提高收敛容差（0.01以下）。 Ansys截面特性值计算 (1) 创建截面的几何模型。描述截面几何形状的面域可以在ANSYS中通过点一线一面的方式直接生成；也可以由外部文件导人。一般通过AUTO CAD来建立几何模型。在AUTO CAD中可将面域分别绘制在不同的图层上，赋予不同的颜色，通过图层开关和颜色等方式进行区分和编辑。有限元分析中，控制网格尺寸和密度对结果的分析有重要影响。在AUTOCAD中，先绘出截面的内外框线，可以用Pedit命令将多段线连成一条多义线(Polyline)，然后用region命令围成面域，也可以导人ANSYS后再形成面(AREA)。 (2) 将AUTOCAD中建立的面域另存为Sat文件，然后在ANSYS中用File—Import—sat方式导人。这种转换方式较方便，模型不会失真变形。 (3) 用Sections--->Beam--->Custom Sections--->write From Areas读取截面，然后在相同目录下用Read Sect Mesh对截面进行网格划分。面进行网格划分。 (4)sections--->Beam--->Plot Sections 即可输出截面特性。 ANSYS默认的单位系是与导人的模型一致的。在图形输出框中的坐标系是Y-Z坐标系。 也可以直接在ANSYS去建立模型去计算截面特性 注意对于箱型截面需要首先进行布尔运算才能计算截面特性值 ANSYS热应力分析例题 实例1圆简内部热应力分折： 有一无限长圆筒，其核截面结构如图13—1所示，简内壁温度为200℃，外壁温度为20℃，圆筒材料参数如表13．1所示，求圆筒内的温度场、应力场分布。 该问题属于轴对称问题。由于圆筒无限长，忽略圆筒端部的热损失。沿圆筒纵截面取宽度为10M的如图13—2所示的矩形截面作为几何模型。在求解过程中采用间接求解法和直接求解法两种方法进行求解。间接法是先选择热分析单元，对圆筒进行热分析，然后将热分析单元转化为相应的结构单元，对圆筒进行结构分析；直接法是采用热应力藕合单元，对圆筒进行热力藕合分析。 /filname,exercise1-jianjie /title,thermal stresses in a long /prep7 $Et,1,plane55 Keyopt,1,3,1 $Mp,kxx,1,70 Rectng,0.1,0.15,0,0.01 $Lsel,s,,,1,3,2 Lesize, all,,,20 $Lsel,s,,,2,4,2 Lesize,all,,,5 $Amesh,1 $Finish /solu $Antype,static Lsel,s,,,4 $Nsll,s,1 $d,all,temp,200 lsel,s,,,2 $nsll,s,1 $d,all,temp,20 allsel $outpr,basic,all solve $finish /post1 $Set,last /plopts,info,on Plnsol,temp $Finish /prep7 $Etchg,tts Keyopt,1,3,1 $Keyopt,1,6,1 Mp,ex,1,220e9 $Mp,alpx,,1,3e-6 $Mp,prxy,1,0.28 Lsel,s,,,4 $Nsll,s,1 $Cp,8,ux,all Lsel,s,,,2 $Nsll,s,1 $Cp,9,ux,all Allsel $Finish /solu $Antype,static D,all,uy,0 $Ldread,temp,,,,,,rth Allsel $Solve $Finish /post1 /title,radial stress contours Plnsol,s,x /title,axial stress contours Plnsol,s,y /title,circular stress contours Plnsol,s,z /title,equvialent stress contours Plnsol,s,eqv $finish ANSYS后处理 1.ANSYS后处理时如何按灰度输出云图？ 1）你可以到utilitymenu-plotctrls-style-colors-window colors试试2）直接utilitymenu-plotctrls－redirect plots 2 将云图输出为JPG 菜单->PlotCtrls->Redirect Plots->To JPEG Files 3.怎么在计算结果实体云图中切面? 命令流 /cplane /type 图形界面操作 <1.设置工作面为切面 <2.PlotCtrls-->Style-->Hidden line Options 将[/TYPE]选项选为section 将[/CPLANE]选项选为working plane 4.非线性计算过程中收敛曲线实时显示 solution>load step opts>output ctrls>grph solu track>on 5.运用命令流进行计算时,一个良好的习惯是: 使用SELECT COMMEND后.........其后再加上ALLSEL......... 6.应力图中左侧的文字中，SMX与SMN分别代表最大值和最小值 如你plnsolv,s,eqv 则 SMX与SMN分别代表最大值等效应力和最小值等效应力 如你要看的是plnsolv,u 则SMX与SMN分别代表位移最大值和位移最小值 不要被S迷惑 mx(max) mn(min) 7.在非线性分析中，如何根据ansys的跟踪显示来判断收敛？ 在ansys output windows 有 force convergence value值和 criterion 值当前者小于后者时，就完成一次收敛你自己可以查看两条线的意思分别是: FL2：不平衡力的2范数 FCRIT：不平衡力的收敛容差， 如果前者大于后者说明没有收敛，要继续计算，当然如果你以弯矩M为收敛准则那么就对应 M L2 和 M CRIT 希望你现在能明白 8.两个单元建成公共节点，就成了刚性连接，不是接触问题了。做为接触问题，两个互相接触的单元的节点必须是不同的。 9.接触单元主要分为有厚度和无厚度的，有厚度主要以desai 为代表，无厚度的则以goodman 为代表。尽管古得曼也提出了相应的本构关系，但是如今goodman 单元成了无厚度接触单元的代名词，相应的本构关系现在也作了 较大的改进。Ansys中接触单元并不是goodman 单元，类似于goodman单元 ansys 里面的接触单元是是通用的，而goodman是一种专业的单元。goodman单元假定两片长为L的接触面以无数微小的切向和法向弹簧所连接，接触面单元与 相邻接触面两边的单元只在结点处有力的联系。单元厚度为零，受力前两接触面完全吻合. 10.怎样检查接触单元的normal direction?是不是打开 plotctrls/symbols/esys on? ANSYS建模分析 报 告 书 课题名称ANSYS建模分析姓名 学号 院系 专业 指导老师 问题描述 在ANSYS中建立如图一所示的支承图，假定平面支架沿厚度方向受力均匀，支承架厚度为3mm。支承架由钢制成，钢的弹性模量为200Gpa，泊松比为0.3。支承架左侧边被固定，沿支承架顶面施加均匀载荷，载荷与支架共平面，载荷大小为2000N/m。要求：绘制变形图，节点位移，分析支架的主应力与等效应力。 图1 GUI操作步骤 1、定义工作文件名和工作标题 （1）定义工作文件名：执行Utility Menu> Jobname命令，在弹出【Change Jobname】对话框中输入“xuhao144139240174”。选择【New log and error files】复选框，单击OK按钮。 （2）定义工作标题：执行Utility Menu> Title命令，在弹出【Change Title】对话框中输入“This is analysis made by “xh144139240174”，单击OK按钮。 （3）重新显示：执行Utility Menu>Plot>Replot命令。 （4）关闭三角坐标符号：执行Utility Menu>PlotCtrls>Window Options命令，弹出【Window Options】对话框。在【Location of triad】下拉列表框中选择“Not Shown”选项，单击OK按钮。 2、定义单元类型和材料属性 （1）选择单元类型：执行Main Menu>Preprocessor>Element Type>Add/Edit/Delete命令，弹出【Element Type】对话框。单击Add...按钮，弹出【Library of Element Types】对话框。选择“Structural Solid”和“Quad 8node 82”选项，单击OK按钮，然后单击Close按钮。 （2）设置材料属性：执行Main Menu>Preprocessor>Material Props>Material Models命令，弹出【Define Material Models Behavior】窗口。双击【Material Model Available】列表框中的“Structural\Linear\Elastic\Isotropic”选项，弹出【Linear Isotropic Material Properties for Material Number1】对话框。在【EX】和【PRXY】文本框中分别输入“2e11”及“0.3”。单击OK按钮，然后执行Material>Exit命令，完成材料属性的设置。 3、创建几何模型 （1）生成两个矩形面：执行Main Menu>Preprocessor>Modeling>Create>Areas>Rectangle>By Dimensions命令，弹出【Create Rectangle by Dimensions】对话框。输入第一个矩形的坐标数值：“X1=0,X2=-0.25,Y1=0.025,Y2=-0.025”,单击Apply按钮。输入第二个矩形的坐标数值：“X1=-0.2,X2=-0.25,Y1=0.025,Y2=-0.225”，单击OK按钮关闭该对话框。 （2）生成两个半圆：执行Main Menu>Preprocessor>Modeling>Create>Areas>Circle>Solid Circle (1) 如果你模拟结构体中裂缝扩展过程的模拟，在Ansys中可以用全解耦损伤分析方法来近似模拟裂缝扩展，我曾用Ansys软件中提供的可以定义10,000个材料参数和单元ekill/alive 功能完成了层状路面体中表面裂缝和反射裂缝在变温作用下的扩展过程的模拟。我模拟的过程相对来说比较简单，模拟过程中我们首先要知道裂缝的可能扩展方向，这样在裂缝可能扩展的带内进行网格加密处理，加密到什么程度依据计算的问题来确定。 (2) 如果采用断裂力学理论计算含裂缝结构体的应力强度因子，建模时只需在裂尖通过命令kscon生成奇异单元即可。Ansys模块中存在的断裂力学模块可以计算I、II、III型应力强度因子(线弹性断裂力学)和J积分(弹塑性断裂力学)，在Ansys中verification里面有一个计算I型应力强度因子的例子vm143，参见该例子就可以了。 (3) 如果通过断裂力学模拟裂缝的扩展过程，需要采用动态网格划分，这方面我没有做，通过Ansys的宏命令流应该可以实现。技术参考可参阅文献：杨庆生、杨卫.断裂过程的有限元模拟.计算力学学报，1997,14(4). (4) 我现在做动荷载作用下路面结构体中应力强度因子的分布规律，我是通过位移插值得到不同时间点处的应力强度因子。如果想这样做，可参阅理论参考中关于应力强度因子计算说明。 1. 讨论两种Ansys求极限荷载的方法 （1）力加载 可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。 （2）位移加载 给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。 希望众高手讨论一下 （1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？ （2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中； 一学习ANSYS需要认识到的几点 相对于其他应用型软件而言，ANSYS作为大型权威性的有限元分析软件，对提高解决问题的能力是一个全面的锻炼过程，是一门相当难学的软件，因而，要学好ANSYS，对学习者就提出了很高的要求，一方面，需要学习者有比较扎实的力学理论基础，对ANSYS分析结果能有个比较准确的预测和判断，可以说，理论水平的高低在很大程度上决定了ANSYS使用水平；另一方面，需要学习者不断摸索出软件的使用经验不断总结以提高解决问题的效率。在学习ANSYS的方法上，为了让初学者有一个比较好的把握，特提出以下五点建议：（1）将ANSYS的学习紧密与工程力学专业结合起来 毫无疑问，刚开始接触ANSYS时，如果对有限元，单元，节点，形函数等《有限元单元法及程序设计》中的基本概念没有清楚的了解话，那么学ANSYS很长一段时间都会感觉还没入门，只是在僵硬的模仿，即使已经了解了，在学ANSYS之前，也非常有必要先反复看几遍书，加深对有限元单元法及其基本概念的理解。 作为工程力学专业的学生，虽然力学理论知识学了很多，但对许多基本概念的理解许多人基本上是只停留于一个符号的认识上，理论认识不够，更没有太多的感性认识，比如一开始学ANSYS时可能很多人都不知道钢材应输入一个多大的弹性模量是合适的。而在进行有限元数值计算时，需要对相关参数的数值有很清楚的了解，比如材料常数，直接关系到结果的正确性，一定要准确。实际上在学ANSYS时，以前学的很多基本概念和力学理论知识都忘得差不多了，因而遇到有一 定理论难度的问题可能很难下手，特别是对结果的分析，需要用到《材料力学》，《弹性力学》和《塑性力学》里面的知识进行理论上的判断，所以在这种情况下，复习一下《材料力学》，《弹性力学》和《塑性力学》是非常有必要的，加深对基本概念的理解，实际上，适当的复习并不要花很多时间，效果却很明显，不仅能勾起遥远的回忆，加深理解，又能使遇到的问题得到顺利的解决。 在涉及到复杂的非线性问题时（比如接触问题），一方面，不同的问题对应着不同的数值计算方法，求解器的选择直接关系到程序的计算代价和问题是否能顺利解决；另一方面，需要对非线性的求解过程有比较清楚的了解，知道程序的求解是如何实现的。只有这样，才能在程序的求解过程中，对计算的情况做出正确的判断。因此，要能对具体的问题选择什么计算方法做出正确判断以及对计算过程进行适当控制，对《计算方法》里面的知识必须要相当熟悉，将其理解运用到ANSYS的计算过程中来，彼此相互加强理解。要知道ANSYS是基于有限元单元法与现代数值计算方法的发展而逐步发展起来的。因此，在解决非线性问题时，千万别忘了复习一下《计算方法》。此外，对《计算固体力学》也要有所了解（一门非常难学的课），ANSYS对非线性问题处理的理论基础就是基于《计算固体力学》里面所讲到的复杂理论。 作为学工程力学的学生，提高建模能力是非常急需加强的一个方面。在做偏向于理论的分析时，可能对建模能力要求不是很高，但对于实际的工程问题，有限元模型的建立可以说是一个最重要的问题，而后 2009-04-28 14:26 ANSYS中查看截面结果的方法 一般情况下，对计算结果后处理时，显示得到的云图为结构的外表面信息。有时候，需要查看结构内部的某些截面云图，这就需要通过各种后处理技巧来获得截面的结果云图。另外，有时候需要获得截面的结果数据，也需要用到后处理的技巧。 下面对常用的查看截面结果的方法做一个介绍： 1. 通过工作平面切片查看截面云图工作平面实现。 这是比较常用的一种方法。 首先确保已经求解了问题，并得到了求解结果。 调整工作平面到需要观察的截面，可通过移动或者旋转工作平面实现。调整时注意保证工作平面与需要观察的截面平行。 在PlotCtrls菜单中设置观察类型为Section，切片平面为Working Plane。也可以通过等效的/type以及/cplane命令设置。 在通用后处理器中显示云图，得到需要查看的云图。 更简单地说，我们只需在显示云图命令前加上下面两条命令就可以了： /CPLANE,1 ! 指定截面为WP /TYPE,1,5 ! 结果显示方式选项 2. 通过定义截面查看截面云图 这种方法也需要用到工作平面与切片，步骤如下： 首先确保已经得到了求解结果。 调整工作平面到需要观察的截面。 在PlotCtrls菜单中设置观察类型为Working Plane，或者使用命令/cplane,1。通过sucr命令定义截面，选择（cplane）。 通过sumap命令定义需要查看的物理量。 通过supl命令显示结果。 3. 通过定义路径查看云图与保存数据 首先确保已经得到了求解结果。 通过path与ppath命令定义截面路径。 通过pdef命令映射路径。 通过plpath、prpath与plpagm命令显示及输出结果。 1. 讨论两种Ansys求极限荷载的方法 （1）力加载 可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。 （2）位移加载 给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。 希望众高手讨论一下 （1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？ （2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中； 3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好，因此，只要条件允许，应该尽量使用六面体单元； 4. 正确选择收敛标准，一般位移控制加载最好用位移的无穷范数控制收敛，而用力控制加载时可以用残余力的二范数控制收敛。在裂缝刚刚出现和接近破坏的阶段，可以适当放松收敛标准，保证计算的连续性； 3. 关于下降段的问题 1）在实际混凝土中都有下降段，但是在计算的时候要特别小心下降段的问题。 2）下降段很容易导致计算不收敛，有时为了计算的收敛要避免设置下降段，采用rush模型。 3）利用最大压应变准则来判断混凝土是否破坏。 4. Solid65单元中的破坏准则 1）采用Willam&Warnke五参数破坏准则 2）需要参数： 单轴抗拉强度，单轴，双轴抗压强度，围压压力，在围压作用下双轴，单轴抗压强度 5. 近来我对混凝土单元进行了一点思考，有一些想法，贴在下面，共同探讨： 1）分析混凝土结构，选择合理的材料特性是建立模型的关键，所以有必要弄清混凝土的材料特性。混凝土是脆性材料,并具有不同的拉伸和压缩特性。典型混凝土的抗拉强度只有抗压强度的8％-15%。 在ANSYS中，对于混凝土单元，材料特性ANSYS要求输入以下数据（为了清楚起见，我将几个系数均译为了中文）：弹性模量、泊松比、张开与闭合滑移面的剪切强度缩减系数、抗拉与抗压强度、极限双轴抗压强度、周围静水应力状态、静水应力状态下单轴与双轴压缩的 ansys命令流----前后处理和求解常用命令之求解与后处理.txt都是一个山的狐狸，你跟我讲什么聊斋，站在离你最近的地方，眺望你对别人的微笑，即使心是百般的疼痛只为把你的一举一动尽收眼底．刺眼的白色，让我明白什么是纯粹的伤害。3 /solu u /solu 进入求解器 3.1 加边界条件 u D, node, lab, value, value2, nend, ninc, lab2, lab3, ……lab6 定义节点位移约束Node : 预加位移约束的节点号，如果为all,则所有选中节点全加约束，此时忽略nend和ninc. Lab: ux,uy,uz,rotx,roty,rotz,all Value,value2: 自由度的数值（缺省为0） Nend, ninc: 节点范围为：node-nend，编号间隔为ninc Lab2-lab6: 将lab2-lab6以同样数值施加给所选节点。 注意：在节点坐标系中讨论 3.2 设置求解选项 u antype, status, ldstep, substep, action antype: static or 1 静力分析 buckle or 2 屈曲分析 modal or 3 模态分析 trans or 4 瞬态分析 status: new 重新分析（缺省），以后各项将忽略 rest 再分析，仅对static,full transion 有效 ldstep: 指定从哪个荷载步开始继续分析，缺省为最大的，runn数（指分析点的最后一步）substep: 指定从哪个子步开始继续分析。缺省为本目录中，runn文件中最高的子步数action, continue: 继续分析指定的ldstep,substep 说明：继续以前的分析（因某种原因中断）有两种类型 singleframe restart: 从停止点继续 需要文件：jobname.db 必须在初始求解后马上存盘 jobname.emat 单元矩阵 jobname.esav 或 .osav : 如果.esav坏了，将.osav改为.esav results file: 不必要，但如果有，后继分析的结果也将很好地附加到它后面 注意：如果初始分析生成了.rdb, .ldhi, 或rnnn 文件。必须删除再做后继分析 步骤：（1）进入anasys 以同样工作名 （2）进入求解器，并恢复数据库 （3）antype, rest （4）指定附加的荷载 （5）指定是否使用现有的矩阵（jobname.trl）（缺省重新生成） kuse: 1 用现有矩阵 （6）求解 multiframe restart:从以有结果的任一步继续（用不着） u pred,sskey, --,lskey….. 在非线性分析中是否打开预测器 sskey: off 不作预测（当有旋转自由度时或使用solid65时缺省为off） on 第一个子步后作预测（除非有旋转自由度时或使用solid65时缺省为on） -- ：未使用变量区 《大型结构分析软件的应用及开发》 学习报告 学院：建筑工程学院 专业班级：工程力学141 姓名：付贤凯 指导老师：姚激 学号：201411012111 1．模型介绍 如下图所示的一桁架结构，受一集中力大小为800N的作用，杆件的弹性模量为200GPa，泊松比为0.3。杆件的截面为正方形达长为1m，横截面面积为1m2。现求它的变形图与轴力图。 图1 桁架模型与受力简图（单位：mm） 2．建模与划分网格 利用大型有限元软件ANSYS，采用Link，2Dspar 1的单元进行模拟，通过网格的划分得到如图2所示的有限元模型。 图2 有限元模型 结合有限元模型中的约束条件为左侧在X与Y方向铰支固定，荷载条件为最右侧处施加向下的集中力P=800N。施加约束与荷载后的几何模型如图4所示。 图3 施加荷载与约束的几何模型 3．位移与轴力图 因在Y方向受力，所以主要做Y方向的位移图，又因为杆件在轴线方向有变形，故在X 方向仍有一定的位移。则图5为变形前后的板件形状。图6为模型沿Y方向的位移图，图7为模型沿X方向的位移图，图8为模型的总位移图。 图4 桁架变形前后形状图 图5 Y方向位移图 图6 X方向位移图 图7总位移图 分析所有的位移图可以看出从以看出左端变形最小，为零，右端变形最大。从总位移图可以看出最大的位移在左下点处，大小为0.164×10?5m。从X方向位移图可以看出，左下点处在X方向位移最大为0.36×10?6。从Y方向位移图可以看出最大位移在左下点处为0.164×10?5。都符合实际情况，图9为模型的轴力图。 图8 轴力图 1.ANSYS后处理时如何按灰度输出云图？ 1）你可以到utilitymenu-plotctrls-style-colors-window colors试试 2）直接utilitymenu-plotctrls－redirect plots 2 将云图输出为JPG 菜单->PlotCtrls->Redirect Plots->To JPEG Files 3.怎么在计算结果实体云图中切面? 命令流 /cplane /type 图形界面操作 <1.设置工作面为切面 <2.PlotCtrls-->Style-->Hidden line Options 将[/TYPE]选项选为section 将[/CPLANE]选项选为working plane 4.非线性计算过程中收敛曲线实时显示 solution>load step opts>output ctrls>grph solu track>on 5.运用命令流进行计算时,一个良好的习惯是: 使用SELECT COMMEND后.........其后再加上ALLSEL......... 6.应力图中左侧的文字中，SMX与SMN分别代表最大值和最小值 如你plnsolv,s,eqv 则SMX与SMN分别代表最大值等效应力和最小值等效应力 如你要看的是plnsolv,u 则SMX与SMN分别代表位移最大值和位移最小值 不要被S迷惑 mx(max) mn(min) 7.在非线性分析中，如何根据ansys的跟踪显示来判断收敛？ 在ansys output windows 有force convergenge valu 值和criterion 值当前者小于后者时，就完成一次收敛 有限元与CAE分析报告 专业： 班级： 学号： 姓名： 指导教师： 2016年 1 月 2 日 简支梁的静力分析 一、问题提出 长3m的工字型梁两端铰接中间1.5m位置处受到6KN的载荷作用，材料弹性模量E=200e9，泊松比0.28，密度7850kg/㎡ 二、建立模型 1.定义单元类型 依次单击Main Menu→Preprocessor→Elementtype→Add/Edit/Delete，出现对话框如图，单击“Add”，出现一个“Library of Element Type”对话框，在“Library of Element Type”左面的列表栏中选择“Structural Beam”，在右面的列表栏中选择3 node 189,单击“OK”。 2设置材料属性 依次单击Main Menu→Preprocessor→MaterialProps>Material Modes，出现“Define Material ModelBehavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural→Linear→Elastic→Isotropic”,出现对话框，输入弹性模量EX=2E+011,PRXY=0.28，单击“OK”。 依次单击Main Menu→Preprocessor→MaterialProps>Material Modes，出现“Define Material ModelBehavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural→Density”弹出对话框，输入DENS为7850 3.创建几何模型 1)设定梁的截面尺寸 a n s y s后处理结果图形 的处理 标准化管理处编码[BBX968T-XBB8968-NNJ668-MM9N] ansys后处理结果图形的处理 对体和面来说，ANSYS默认的结果输出格式是云图格式，而这种彩色云图打印为黑白图像时对比很不明显，无法表达清楚，这对于发表文章来说是非常不便的。发文章所用的结果图最好是等值线图，并且最好是黑白的等值线图。笔者原来进行这项工作时一般借用photoshop等第三方软件，很麻烦，并且效果不好。现通过摸索，发现通过灵活运用ansys本身也能实现这项功能。现将步骤写给大家，感谢simwe对我的帮助。 （1）将要输出的结果调出，这时为彩色云图； （2）将云图转换为等值线图的形式 GUI：plotCtrls—>Device Options—>[/DEVI]中的vector mode 选为on 命令：/DEVICE,VECTOR,1 这时结果为彩色等值线，若直接输出，打印为黑白图像时仍然不清晰，为此需进行以下几步将图像转换为黑白形式； （3）将背景变为白色 命令：jpgprf,500,100,1 /rep （4）对等值线中的等值线符号（图中为A,B,C等）的疏密进行调整 GUI：plotCtrls—>Style—>Contours—> Contours Labeling 在Key Vector mode contour label 中选中on every Nth elem，然后在N= 输入框中输入合适的数值，例如5，多试几次，直到疏密合适 命令：/clabel,1,5 （5）将彩色等值线变为黑色 GUI：plotCtrls—>Style—>Colors—>Contours Colors 将Items Numbered 1，Items Numbered 2等复选框中的颜色均选为黑色，图像即可变为黑白等值线图像命令：/color,cntr,whit,1 等等 （6）最后一步：出图 GUI：plotCtrls—>Capture Image 希望对大家能有所帮助。 一个使生成的图片在word里面比较好看的方法： 1、Plotctrls>Redirect Plots>To png file 2、选“Force White BG and Black FG",然后把Pixle resolution 换到1200！ !!!!!~~~~~!!!!!~~~~~!!!!!~~~~~!!!!!~~~~~!!!!!~~~~~!!!!! !!!!!~~~~~~~~~ansys数据处理的相关命令流~~~~~~~~~~~!!!!! !(1)数据输入的相关命令 !利用*TREAD命令读取数据文件并填充TABLE表格 *TREAD, Par, Fname, Ext, --, NSKIP !以下利用*TREAD命令读取1维数据表格 !tdata.txt文本文件含有如下内容 STRAIN STRESS 00 0.0025 0.0046 0.0067 *DIM,Ttxy,table,4,1,,TIME,ACEL *TREAD,Ttxy,tdata,txt,,1 !以下利用*TREAD命令读取2维数据表格 !要特别注意2维数据的行数 !tdata.txt文本文件含有如下内容 TIME X Y Z 0000 0.020.10.20.3 0.040.20.40.6 0.060.30.60.9 !希望输入地震波激励，X、Y、Z三个方向 *DIM,Ttxy,table,3,3,,TIME,ACEL *TREAD,Ttxy,tdata,txt,,1 !以下利用*TREAD命令读取3维数据表格 !tdata.txt文本文件含有如下内容 TEMP X Y Z 0000 0.020.10.20.3 0.040.20.40.6 0.060.30.60.9 5000 0.030.20.30.4 0.050.40.60.8 0.070.60.90.9 !希望读取不同温度下，不同时刻的泊松比 *DIM,Ttxy,table,3,3,2,TIME,NUXP,TEMP *TREAD,Ttxy,tdata,txt,,1 !利用*SREAD命令读取字符文件 *SREAD, StrArray, Fname, Ext, --, nChar, nSkip, nRead 页: 1 三维托架实体受力分析 ANSYS软件是融结构、流体、电磁场、声场和耦合场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国ANSYS公司开发，它能与多数CAD软件接口，实现数据的共享和交换，如PRO/E、UG、I-DEAS、CADDS及AutoCAD等，是现代产品设计中的高级CAD工具之一。 题目：1、三维托架实体受力分析：托架顶面承受50psi的均匀分布载荷。托架通过有孔的表面固定在墙上，托架是钢制的，弹性模量E=29×106psi，泊松比v=.试通过ANSYS输出其变形图及其托架的von Mises应力分布。 题目1的分析。先进行建模，此建模的难点在对V3的构建（既图中的红色部分）。要想构建V3，首先应将A15做出来，然后执行Main Menu>Preprocessor>Modeling>Operate>Booleans>Add>Volumes命令，将所有的实体合并为一个整体。建模后，就对模型进行网格的划分，实行Main Menu>Preprocessor>Meshing>MeshTool,先对网格尺寸进行编辑，选，然后点Meshing，Pick all进行网格划分，所得结果如图1。划分网格后，就可以对模型 施加约束并进行加载求解了。施加约束时要注意，由于三维托架只是通过两个孔进行固定，故施加约束应该只是针对两孔的内表面，执行Main Menu>Solution>Define Loads>Apply>Structrual>Displacement>Symmetry >On Areas 命令，然后拾取两孔的内表面，单击OK就行了。施加约束后，就可以对实体进行加载求解了，载荷是施加在三维托架的最顶上的表面的，加载后求解运算，托架的变形图如图2。 图1、托架网格图 图2输出的是原型托架和施加载荷后托架变形图的对比，虚线部分即为托架的原型，从图2可看出，由于载荷的作用，托架上面板明显变形了，变形最严重的就是红色部分，这是因为其离托板就远，没有任何物体与其分担载荷，故其较容易变形甚至折断。这是我们在应用托架的时候应当注意的。 学习ansys的一些心得 学习ansys的一些心得（送给初学者和没有盟币的兄弟） 1 做了布尔运算后要重画图形（删除实体）时：需拾取Utility Menu>Plot>Replot 2 标点的输入是在英文状态下，―,‖。 3 线段中点的建立：Modling>Creat>Keypoints>Fill between kps 4 还不会环形阵列。 5 所谓杆系结构指的是长度远远大于其他方向尺寸（10：1）的构件组成的结构，如连续梁，桁架，钢架等。 6 静力学分析的结果包括结构的位移，应变，应力和反作用力等，一般是使用POST1处理（普通后处理器）和查看这些结果。 7 干系结构的静力学分析—平面桁架的建模，用NODE（节点），ELEMENT(元素)创建。复杂体积的建模一般用KPS(关键点)，LINE（Straight line—直线），再生成面，再生成体。 8 如果输入的数据单位是国际单位制单位，则输出的数据单位也是国际制单位。 9 创建正六边形：Creat>Areas>Polygon>Hexagon.指定中心和半径。 10 由面沿线挤出体：Modling>Operate>Extrude>Areas>Along Lines. 11 Ansys中没有Undo命令.需及时保存数据库文件. Def Shape Only:只显示变形图.Def + Undeformed:显示未变形的图.Def + Udef egde:显示未变形的图形的边界. 13 用等高线显示：Plot Results>Contour Plot>Nodal Solu. 14 模态分析用于分析结构的振动特性，即确定结构的固有频率和振型，它也是谐响应分析，瞬态动力学分析以及谱分析等其他动力学分析的基础。 15 Ansys的模态分析是线型分析。任何非线型分析，例如，塑性，接触单元等，即使被定义了也将被忽略。 16 平面桁架：Beam(2D elastic 3) 厚壁圆筒：Solid(8 node 13)>Options(K3—Plane strain) 17 一般材料的弹性模量（EX）：2e11.泊松比(PRXY)：0.3.密度：7800 18 做完静力学分析后，再做模态分析时，要再次求解，同时预应力效果也应该打开（PSTRES,on）.可以在命令行中输入：pstres,on 也可以用菜单路径：Solution>Analysis Type>Analysis Options. 19 弹簧阻尼器单元：Combination-Spring damper 14. 20 接触问题属于状态非线性问题，是一种高度非线性行为，需要较多的计算资源。接触问题有两个基本类型：刚体-柔体的接触，柔体-柔体的接触（许多金属成型的接触问题）。在刚体-柔体的接触问题中，有的接触面与它接触的变形体相比，有较大的刚度而被当做刚体。而柔体-柔体的接触，是一种更普遍的类型，此时两个接触体具有近似的刚度，都为变形体。 21 1 点-点接触：过盈装配问题是用点点接触单元模拟面面接触的典型例子。 2 点-面接触：不必预先知道准确的接触位置，接触面之间也不需要保持一致的网格，并且允许有较大的变形和相对滑动。典型实例：模拟插头插入插座里。 3 面-面接触：刚性面作为目标面，柔性面作为接触面。 22 打开自动时间步长：Solution>Load Step Opts>Time Frequenc>Time And Substps. 有限元分析作业 作业名称扳手静态受力分析 姓名 学号 班级 宁波理工学院 题目：扳手静态受力分析: 扳手的材料参数为：弹性模量E=210GPa,泊松比u=0.3：此模型在左侧内六角施加固定位移约束，在右侧表面竖直方向上施加6 10 48 N的集中力。 模型如下图： 1-1 1.定义工作文件名和文件标题 （1）定义工作文件名:执行File-Chang Jobname-3090601048 （2）定义工作标题：执行File-Change Tile-3090601048 （3）更改工作文件储存路径：执行File-Chang Directory-E:\ANSYS 2.定义分析类型、单元类型及材料属性 （1）定义分析类型，执行Main Menu-Preferences,如下图所示： 2-1 （2）定义单元类型，执行Main Menu-Preprocessor-Element Type-Add弹出Element Type 对话框.如下图所示： 2-2 （3）定义材料属性 执行Main menu-Preprocessor-Material Props-Material models,在Define material model behavior对话框中，双击 Structual-Linear-Elastic-Isotropic.如下图所示： 2-3 3.导入几何模型 将模型导入到ANSYS，执行File-Import—PRAR…—浏览上述模型，如下图所示： 3-1 3-2 4. 网格划分 执行Main Menu-Preprocessor-meshing-Mesh Tool命令，考虑到零件的复杂性，采用智能网格划分，精度为1，其他选项为默认，如下图所示： 4-1 ANSYS使用心得体会 本次结构力学课程设计是学习使用ANSYS软件对框架结构内力进行计算，在未学习该软件前，对于此类问题，通常会采用力矩分配法来进行计算，计算过程繁复，计算量大。导致过程缓慢。 通过对ANSYS软件的学习和了解，知道了它的一些明显的优点。 相对于其他应用型软件而言，ANSYS作为大型权威性的有限元分析软件，对提高解决问题的能力是一个全面的锻炼过程，是一门相当难学的软件，因而，要学好ANSYS，对我们提出了很高的要求，一方面，需要我们有比较扎实的力学理论基础，对ANSYS分析结果能有个比较准确的预测和判断，可以说，理论水平的高低在很大程度上决定了ANSYS使用水平；另一方面，需要我们不断摸索出软件的使用经验不断总结以提高解决问题的效率。 刚开始接触ANSYS时，没有限元，单元，节点，形函数等的基本概念没有清楚的了解话，会感觉还没入门，只是在僵硬的模仿，即使已经了解了，必要先反复看几遍书，加深对有限元单元法及其基本概念的理解。 ANSYS在对结构力学的静力学分析非常方便，用来求解外载荷引起的位移、应力和力。静力分析很适合求解惯性和对结构的影响并不显著的问题。ANSYS 程序中的静力分析不仅可以进行线性分析，而且也可以进行非线性分析，如塑性、膨胀、大变形、大应变及接触分析。 但是学习的过程是充满烦恼和惊喜的，因为总是会碰到许多的新问题，需要较好的耐心去解决这些问题，这是在学习过程中遇到的最大的难题。然而，在解决问题之后，就会有恍然大悟的喜悦，可以说是痛苦和快乐并存的。所以对于初学者，缺乏经验是非常难的。必须保持良好的心态，对于不断出现的ERROR提示要坚定自己的信心，坚信自己可以解决这些问题。所有困难都会迎刃而解。 本次的学习让我认识到了提高建模能力是非常急需加强的一个方面。在做偏向于理论的分析时，可能对建模能力要求不是很高，但对于实际的工程问题，有限元模型的建立可以说是一个最重要的问题，而后面的工作变得相对简单。建模能力的提高，需要掌握好的建模思想和技巧。 ANSYS软件是一款在建模等方面非常实用的软件，本次的学习我其实并没有完全熟练地掌握它的应用，以后还要加强对它的学习，相信在以后的学习和工作中会带来巨大的便利。 有 限 元 分 析 作 业 作业名称 输气管道有限元建模分析 姓 名 陈腾飞 学 号 3070611062 班 级 07机制（2）班 宁波理工学院 题目描述： 输气管道的有限元建模与分析 计算分析模型如图1所示 承受内压：1.0e8 Pa R1=0.3 R2=0.5 管道材料参数：弹性模量E=200Gpa；泊松比v=0.26。 图1受均匀内压的输气管道计算分析模型（截面图） 题目分析： 由于管道沿长度方向的尺寸远远大于管道的直径，在计算过程中忽略管道的断面效应，认为在其方向上无应变产生。然后根据结构的对称性，只要分析其中1/4即可。此外，需注意分析过程中的单位统一。 操作步骤 1.定义工作文件名和工作标题 1.定义工作文件名。执行Utility Menu-File→Chang Jobname-3070611062，单击OK按钮。 2.定义工作标题。执行Utility Menu-File→Change Tile-chentengfei3070611062，单击OK按钮。 3.更改目录。执行Utility Menu-File→change the working directory –D/chen 2.定义单元类型和材料属性 1.设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK 2.选择单元类型。执行ANSYS Main Menu→Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 8node 82 →apply Add/Edit/Delete →Add →select Solid Brick 8node 185 →OK Options…→select K3: Plane strain →OK→Close如图2所示，选择OK接受单元类型并关闭对话框。 图2 3.设置材料属性。执行Main Menu→Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic，在EX框中输入2e11，在PRXY框中输入0.26，如图3所示，选择OK并关闭对话框。 图3 3.创建几何模型 1. 选择ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：input:1(0.3,0),2(0.5,0),3(0,0.5),4(0,0.3) →OKANSYS热应力分析经典例题
ANSYS中重要的后处理
ANSYS分析报告
ANSYS新手入门学习心得
ANSYS学习心得
ansys使用技巧(后处理)
ansys心得
ansys命令流----前后处理和求解常用命令之求解与后处理
ANSYS分析报告
ansys实用的后处理
ANSYS分析报告分析
ansys后处理结果图形的处理
ansys数据处理总结
Ansys受力分析
学习ansys的一些心得
有限元分析Ansys大作业
ANSYS使用心得体会
ansys有限元分析作业经典案例