有限元非线性回弹地分析报告报告材料方法

1、基于ansys金属板料冲压成型回弹的分析方法关键字：板料冲压成型|回弹摘要：金属板料冲压成型由于板料冲压成型过程十分复杂，使得对其成型控制非常困难，随着有限元模拟技术的迅速发展，利用有限元软件模拟金属板料冲 压成型过程及分析其回弹量已成为可能，本文介绍了利用有限元软件Ansys对金属板料冲压成型过程进行模拟和回弹分析的方法。金属板料冲压成型是利用金属塑性变形的特点，通过对金属板料施加压力， 使其产生塑性变形从而获得所需要的形状1。由于板料冲压成型过程包含大位 移、大变形等十分复杂的物理现象，使得对其成型控制非常困难，以前更多的是 通过反复试验的方法制造出合乎要求的产品， 其过程需要花费大量的时

2、间和经费2。随着有限元模拟技术的迅速发展，利用有限元软件模拟金属板料冲压成型 过程及分析其回弹量已成为可能，本文介绍了利用有限元软件Ansys对金属板料冲压成型过程进行模拟和回弹分析的方法。1 Ansys软件分析回弹问题的流程图1 Ansys软件分析回弹问题的流程图有限元软件Ansys分析回弹问题的流程如图1所示3,是利用显示一一隐 式顺序求解技术来分析回弹问题，为了要得到最终的结果，需要先进行显式计算， 然后将显式计算的最终结果输入到隐式模型中进行隐式计算。2有限元分析2.1有限元模型建立有限元模型是模拟板材成形工艺过程的基础，具体的有限元模型要根据 工件、模具的实际结构形状来建立，要能够反

3、映出实际加工过程中被成形板材、 模具之间的相对位置，以及它们之间的相对运动关系。本文主要讨论的是实现冲 压成型回弹分析的模拟方法，故采用比较简单的模型来说明这个问题， 模型如图 2所示，凸模和凹模表面为非均匀的样条曲线，板料为规则的长方形薄扳。由于 是板材，建模时只需建为平面，厚度在定义单元时再输入。2.2单元和材料及网格剖分表1材料性能参数表密度弹性模量（GPa）泊松比屈服极限（MPa）强度极限（MPa）钢板7865188.780.324352598凸模、凹模78002000.3板材图3眶珞划分在实际成型加工时，相对于板材的变形，模具的变形可以不作考虑，除非研 究问题是计算模具的应力、应变，

4、在这种条件下，输入材料参数时可以把凸模、 凹模定义为刚体，这样处理可以减小很多的计算量，并且计算精度不会差多少， 但材料的弹性模量、泊松比仍需按真实值定义，不能够认为是刚体就需要定义为 很大的值，否则会在计算时出现不收敛的情况，材料的力学参数如表1所示，采 用 Plastic -Kinematic 材料模型。凸模、凹模使用LS_DYNA中的soild164三维实体单元，被加工板材使用 shell163壳单元，板料厚度为1mm。每种单元都有很多种算法可以选择，选择 适合的算法可以提高计算的速度和精度。 在定义沿板厚度方向积分点参数时，如 果仅仅是作冲压成型分析，建议设定为 3个积分点;如果要做回

5、弹分析，建议选 择设定为7个积分点，这样分析的结果更精确，但计算时间也会增加。网格划分 采用软件自适应网格划分，模型网格划分如图3所示。2.3定义接触面凸模和凹模与板材的接触都采用面一面接触方式，并选择自动接触类型，摩擦系数取0.1。需注意，板材必须定义为从面。2.4施加载荷有限元数值模拟时要尽量反映成型过程中的真实工艺情况，有了真实有效的工艺参数才能保证数值模拟结果的正确，按照实际加工时的参数来设定凸模的速 度，求解时间为0.1秒，结果输出文件.RST,输出步数为100。2.5成型计算结果伽想料刊变障耳(b)板料应变的等值线图(c) node2位移隨时间变比国(d)板料Von-Jlises效

6、宜力图图4显式计算结果实现显式求解计算有两种方式4，一种可以直接在Ansys软件中调用 Ls_Dyna显式求解器计算，优点是比较快捷方便，缺点是对前面所有步骤的设 置不能重新编辑、修改。另一种方式是在 DOS状态下，通过DOS命令来调用 Ls_dyna求解器来计算。但这种方式必须使用Ls_Dyna特定的文件格式关键字文件，关键字文件可以通过 Ansys软件导出，关键字文件的修改和编辑可 以通过操作系统中的写字板功能来完成。板料冲压成型显式计算的结果如图4所示，分别给出了板料的变形图、变形云图、板料的Vo n-Mises应力云图和板料上节点node2的位移随时间变化的曲线。通过这些图可以显示出在

7、冲压成型 的过程中，板料上应力、应变及位移的变化量。node2位移曲线显示了板料的 运动已经逐渐变小，足以进行隐式回弹分析。如果在求解结束时有一个大的振荡 运动，应该使用一定形式的阻尼重新进行显式分析，否则后面的隐式分析会出现 计算不收敛而得不到结果的现象，其原因可能是载荷的不真实或是其它的原因导 致无限大的加速度。利用Ls-dyna求解器计算的结果文件有两种，一种是以.rst和.his为结尾的Ansys的数据文件另一种是d3plot和d3thdt的Ls_dyna的结果文件。可以通 过Ansys软件的通用后处理器、时间历程处理器和Ls-Dyna模块自带的Ls-post 后处理器来分别查看结果。

8、3回弹计算调用Ansys软件隐式求解器进行回弹计算的步骤为：(1) 更改文件名，以免后面计算的结果文件将前面显式求解的结果文件覆盖。(2) 将显式单元类型转换为具有适当属性的隐式单元类型。在Ansys中，有相应的显式和隐式单元类型，在回弹计算时，所有的单元必须转化为相应的隐式 单元类型，前面的单元类型 Shell163、Solid164 转化为Shell183、Solld185 单元类型。(3) 利用Ansys软件隐式求解器进行回弹分析时，只能激活线弹性材料特性， 因此在显式分析部分板材使用的塑性材料特性必须被删除。(4) 关闭单元的形状检查，因为成型过程中板材单元可能承受极大的变形， 为了得

9、到结果必须关闭单元的形状检查。在前面显式分析时，使用到了刚性单元，在隐式计算时不需要，故需删 除。(6)重新定义边界条件，由于显式分析和隐式分析对板材约束条件不同，所 以板材的边界条件需要重新定义。(7)将软件中几何非线性开关打开，因为在隐式求解的开始，板材一般都有 较大的几何变形。完成以上的步骤后，进行隐式回弹求解。在隐式回弹求解中，其分析的起点 是显式求解的最后变形形状。图5板料回弹后应变云图4回弹计算结果分析板材回弹计算后的应变等值线云图如图 5所示，图6是板料的节点分布示 意图，为了比较板料回弹量在总变形中所占的比值，选取图 6中沿板料长度方向 直线上的33、41、49、57、65节点

10、为研究对象，分别将各节点在显式求解后变形量、隐式求解后的变形量及两者的比值在表2中列出，可以比较出，这些节点在显式求解后和隐式求解后的节点位移值及回弹应变相对于总应变的比值如 表所示:4W IS-斗:HD0ESNOH E wrrpiW -1TW50 .« «ZV 町匾厂 .5¥F - .08.9 EF &Z-BUTF»图6板料节点分布示意图表2各节点变形量及弹性变形与总变形量比值节点编号334149570.015850.02650.02960.02810.01190.005580.004030.003180.003580.00539/35.2%1

11、5.2%10.7%12.7%45.3%在选择的节点中，Node65最大的弹性变形占总变形的比值为 45.3%，其 它所选节点变形比值也在10%以上，可见回弹的存在会严重影响板料的最终形 状，因此在模具设计阶段，应充分加以考虑。5小结本文以一个简单的模型为例，介绍了利用有限元软件Ansys模拟金属板料成型过程以及回弹分析的方法，该方法不仅可以预测成型缺陷，而且可以对一些 工艺参数进行比较和优化。但在实际加工过程中，各种产品模型形状复杂度和成 型工艺参数都不相同，利用 Ansys软件提供的参数化设计语言(APDL)和用户界 面设计语言语言(UIDL )进行二次开发，开发出相应的模拟模块，可以很大的

12、提 高分析效率，是今后的发展趋势。参考文献1 王孝培冲压设计资料.重庆大学出版社.1983年.2 钟志华.薄板冲压成形过程的计算机仿真与应用.北京理工大学出版社.1998年.3 任重.ANSYS实用分析教程.北京大学出版社.2003年6月.4 北京理工大学ANSYS/LS-DYNA算法基础和使用方法北京理工大学.1999年3月.LS-DYNA使用指南第十四章发表时间：2007-7-30 作者：安世亚太 来源：e-works关键字：显式有限元LS-DYNA ANSYS第十四章 显式一隐式顺序求解某些工程处理的模拟需要既有隐式分析（ANSYS ）功能，又有显式分析（ANSYS/LS-DANA ）功

13、能。要求解这类问题，需要采用这二种解法，即一个显式求解后再进行隐式求解，或者相反。ANSYS/LS-DYNA是一个显式动力程序，用于求解实时动态问题。如果一个工程过程包括主要是静态或准静态的阶段（如在 动态阶段之前预加载，或在金属成形阶段之后线性弹性回弹），则在这个阶段最好使用ANSYS隐式程序进行分析。综合ANSYS隐式求解器和ANSYS/LS-DYNA 显式求解器的过程是一个特别有用的工具，可以用来模拟许多复杂的物理现象。本章将描述显式-隐式过程，而在第十五章将简述隐式-显式过程。14.1显式隐式顺序求解在薄板金属成形工序中，回弹变形为一个基本参数，它使成型模具设计变得大大复杂化。最好将回

14、弹变形定义为变形部分的尺寸变化，回弹变形在卸载时由弹性变形产生。在多数动态金属成型工序中，高度 非线性变形过程导致在坯料材料中产生大量的弹性应变能。在与模具动态接触时坯料中开始存储的弹性能 在成型压力消失之后释放，这个能量释放是弹性回弹的驱动力，使坯料向着原有几何构形变形。因此，在 板金属成型过程中最后部分形状不仅仅取决于模具的轮廓形状，也取决坯料部分在塑性变形时存储在该部 分的弹性能总量。因为在变形部分存储的弹性能总量是许多过程参数（如材料性质，两表面间的载荷）的 函数，在成形时期预测回弹是特别复杂的。它给设计者和分析者提岀重要问题，他们必须准确估计成型过 程将会产生的回弹量以便得到最后所需

15、的形状。此外执行一个显式-隐式连续求解会有很大好处。用显式 求解模拟动态成形过程，然后用隐式求解模拟线性回弹变形，使得薄板成形工序达到严格的设计容差。在显式隐式连续求解中，必须首先运行 ANSYS/LS-DANA 程序来模拟金属成形过程。在显式金属 成形分析中，薄板金属坯料的变形，必须由 SHELL163/SO LI D164 单元构成。然后仅将这些单元的变形形 态、应力和厚度传递给相应的隐式 ANSYS单元（SHELL181和SOLID185 ）。只要对坯料给定边界条件， 便能够执行工件的弹性回弹隐式模拟。注一一在隐式阶段只考虑弹性特性。这意味着分析从以前的状态（形状）进行纯弹性卸载。这一位

16、移 被认为是以前的形状在量上的变化。详细描述显式-隐式求解过程如下；1 采用Jobnamel运行显式分析如前所述。为了继续在ANSYS隐式分析中分析回弹影响，必须采用SHELL163对工件作网格剖分。 另外，对于SHELL163单元使用下列单元算法之一：KEYOPT （1 ）=1，2，6，7，8，9，10 ,11，12 （推荐使用10和12 ）。然后求解并终止这一分析。在开始ANSYS隐式分析之前，通常要小心地检查ANSYS/LS-DYNA 显式分析结果。特别要检查显式运行结束时是否有任何不需要的动力影响留在结构内。2 .将显式分析存储到 Jobname.DB文件内。命令：SAVEGUI: U

17、tility Menu>File>Save as注一如果不能在此时存储Jobname.DB文件，则以后将不能存储这个显式运行的数据库，只存储下面隐式运行所用的数据库文件。3 .将作业名改为Jobname 2 ，以防止覆盖显式结果文件。命令：/FILNAME ，Jobname2GUI: Utility Menu>File>Change Jobname4. 重新进入前处理命令：/PREP7GUI: Main Menu>Preprocessor5. 将显式单元转换到对应的隐式单元类型。（注意，2-D显式单元PLANE162，KEYOPT （ 3） =1的 复合材料壳单元

18、SHELL163不能在顺序求解中使用）对应的显式-隐式单元类型对是：隐式单元类型显式单元类型LINK8LINK160BEAM4BEAM161SHELL181SHELL163SOLID185SOLID164COMBIN14COMBI165MASS21MASS166LINK10LINK167虽然所有的显式单元类型都被变换，但仅仅 SHELL163单元和SOLID164单元的数据（SHELL163 , 单元应力和厚度； SOLID164，应力）传递给 SHELL181单元和SOLID185 单元。（通过 RIMPORT ;步 骤12）命令：ETCHG ，ETIGUI: Main Menu>Pr

19、eprocessor>Switch Elem Type6.重新对显式单元类型转换来的任何一个隐式单元定义关键选项、实常数、材料性质、边界条件和载 荷曲线。（对于SHELL163单元，它转换成SHELL181单元，不需要重新定义实常数，但需要重新定义其 他值。对于SOLID164单元，它转换成SOLID185单元，需要设KEYOPT（2）=1来定义统一缩减积分选 项。）显式单元的TYPE号、REAL号和MAT号保持不变，但实际的 KEYOPT选项和实常数值被重置为零 或缺省值。命令：KEYOPTRMP等等GUI: Main Menu>Preprocessor>Element T

20、ype/Real constants/Material Properties/Loads注：在ANSYS隐式阶段仅线弹性材料性质保持有效（正如 MP命令给定的）。删除从 ANSYS/LS-DYNA 运行中的非弹性性质（正如 TB命令给定的）。7 关闭形状检查，因为在显式分析时单元可能承受相当大的变形。命令：SHIPP ，OFFGUI: Main Menu>Preprocessor>Checking Ctrls>Shape Checking8. 隐式单元重新定义为变形构形命令：UPGEOMGUI: Main Menu>Preprocessor>-Moleling-U

21、pdate Geometry9. 不选择或删除任何不需要的单元（主要是在显式分析中构成任何刚性体的单元），或者将它们转换成空单元。任何未被选择、删除、转换成空单元或没有被转换成在ANSYS中保持作用的隐式单元的显式单元都将产生错误并终止分析。如果在显式分析中刚性体由SHELL163单元或SOLID164构成，这些单元在用RIMPORT 命令（见第十二步）输入应力和厚度（从 SHELL163输入到SHELL181 ）之前必须从选 择系列中排除、删除或转换成空单元，否则隐式分析将被终止。命令：ESELEDELEGUI: Utility Menu>Select>Entities 或 Ma

22、in Menu>Preprocessor>Delete>Elements10. 重新进入求解处理器命令：/SOLUGUI: Main Menu>Solution11. 通过修改或增加在显式分析中定义的边界条件在模型上设置必要的约束（例如，在金属变形分析 中，需要约束板料）。命令：D，等等GUI: Main Menu>Solution>-Loads-Apply>-Structural-Displacement>On Nodes12. 导入应力和改变壳的厚度（仅从SHELL163单元转换到SHELL181单元和从SOLID164至USOLID185）

23、。对于SOLID164 ,仅传递应力，对于SHELL163，传递应力和厚度。在转换到隐式分析之前, 将变形的积分点厚度求平均数。这个厚度值覆盖采用实常数定义的值。命令：RIMPORTGUI: Main Menu>solution>-Loads-Apply>-Structural-Other>lmport Stress13. 打开大变形效应命令：NLGEOM , ONGUI: Main Menu>Solution>Sol ' n Controls14. 求解并结束分析命令：SOLVEFINISHGUI : Main Menu>Solution&g

24、t;-Solve-Current LS Main Menu>Finish一旦求解完，可以用标准的 ANSYS后处理功能观察结果。 以下是一个实例，执行显式-隐式连续求解输入流。/batch,listresum,stamp1,db !恢复显示数据库（早先运行显式作业）/filename,stamp2 !修改作业名，使显式求解结果不被覆盖/prep7etchg,eti !将显式单元转移成隐式单元r,1,.01,.01,.01,.01 !变换实常数（用 RIMPORT得变形的壁厚）ddel,all !删除显式分析中的约束载荷tbdelete,all,all !删除材料模型（保留线性材料性质）

25、upgeom,stamp1,rst !重新为变形构形定义隐式单元d,n1,ux,0.0,uy,uz !对模型设置必要的约束（至少3个非线性节点） d,n2,ux,0.0,uyd,n3,ux,0esel,s,ename,181 !仅选择相应的SHELL181单元（非早先的刚性体） nsle,s,1 !仅选择SHELL181单元节点shpp,off !关闭形状检查finish/solutionnelgeom,on ! 接通几何非线性rimport,dyna,stamp1,rst !输入显示分析求得的应力和壁厚savesolvefinish分页14.2回弹分析疑难解答回弹有时候在数值上很难解决，引起收敛困难的典型原因和相应解决措施如下：在显式运行（LS-DYNA）中，较高的凸模速度会引起动态影响，所以在显式运行中采用较低的凸模速度， 在显式分析中可以重新定义速度来实