ansys模拟起皱现象的主要步骤与分析说明.docVIP

模拟起皱现象的主要步骤与分析说明 模拟试验采用三维薄壳单元SHELL163进行ANSYS/LS-DYNA程序的交互式求解过程: 5.2.1启动ANSYS/LS-DYNA 启动ANSYS/LS-DYNA之后，选择ANSYS Mechanical，在File Management标签，然后设置工作目录为刚建的目录，并输入工作名，即可以进入ANSYS/LS-DYNA，如图5—2。 本实验在出现的ANSYS Launcher的界面中之所以要选用显式动力分析程序—LS-DYNA来模拟一个板料冲压的成型过程，主要是因为该软件是以Lagrange算法为主，兼有ALE和Euler算法；以显式求解为主，兼有隐式求解功能；以结构分析为主，兼有热分析、流体—结构耦合功能；以非线性动力分析为主，兼有动力分析功能（如动力分析前的预应力计算和薄板冲压成形后的起皱计算）；几乎所有的成形过程，如冲压、锻造、铸造等都可用LS-DYNA模拟，它计算的可靠性已经被无数次实验所证明，因此，我们在实际的板料冲压过程中要使用该软件分析包。 另外，在ANSYS/LS-DYNA程序显示动态分析中，选择ANSYS Mechanical/LS-D- YNA模块，该模块是应用于机械领域，主要研究板料成形过程中的应力、应变的变化，对拉深过程中的起皱、回弹以及拉裂现象都能很好的进行ANSYS/LS-DYNA的显式—隐式序列的求解。因此，我们用软件模拟冲压过程时要使用该软件模块。 5.2.2前处理建模 ⑴设置Preference选项 ①选择Structural ②选择LS-DYNA explicit 这样，以后显示的菜单完全被过滤成ANSYS/LS-DYNA的输入选项。再定义一种显式单元类型，即可激活LS-DYNA求解。 ⑵定义单元类型、实常数和材料模型 ①定义单元类型。在该对话框的单元类型库Library of Element Type中分别选择LS-DYNA Explicit、Thin Shell163壳单元，如图5—3。选择单元算法，在SHELL163壳单元选项对话框，选择Belytschko-Wong算法，即可完成单元的定义，如图5—4。 在本实验中，板料模型都采用显式壳单元SHELL163,且使用能忍弯曲能力Bely-tschko—Wong—Chiang算法。板料模型采用薄壳单元，主要是由于SHELL163是一个4节点显式结构薄壳单元，有弯曲和模特征，可加平面和法向载荷。单元在每个节点上有12个自由度：在节点X、Y和Z方向的平动，加速度，速度和绕X、Y和Z轴的转动。该单元支持显式动力学分析所有非线性特性。而Belytschko—Wong—Chiang这种单元算法是显式薄壳单元SHELL163所提供的12种算法之一。它在弯曲情况下无效，在预防沙漏方面很有效。 图5—2 启动ANSYS界面 ②定义实常数。在实常数号设置框中输入实常数号，分别在剪切因子、积分点数、壳后中输入5/6、5、0.001，即可完成外壳的壳单元实常数的定义，如图5—5。 在这里，每种单元的实常数可能有好几个实常数，也可能没有实常数。对于本实验，需要定义的实常数包括剪切因子、积分点和壳厚。在SHELL163单元中提供的如下实常数：剪切因数SHRF，它的推荐值为5/6，如果不指定则采用默认值1；NIP是通过单元厚度的积分点数值，最大值为100。如果NIP输入值为0或空，ANSYS会默认积分值为2；在本实验中NIP取经验值5；T1～T4是四个节点中每个节点处的壳厚度。本实验在T1节点处输入板料壳的厚度为0.001m。 图5—3 定义单元类型 ③定义材料模型。各材料模型的参数见表5—1。 材料模型1：凸模采用刚体。Rigid Properties for Material Number 1对话框入相应的材料参数以及限定运动的自由度，就可以完成材料模型1的设置，如图5—6，5—7。 在本实验里，对于凸模采用刚体，它限制了除Y方向平动外的所有自由度，因此只允许它在Y方向上运动。之所以要把凸模定义成刚体，主要因为刚体模型在显式动力学分析中有非常重要的意义，用刚性模型定义有限元模型中钢硬部分可以大大缩减显式分析的计算时间，这是由于定义了刚体后，缸体内所有节点的自由度都耦合到钢性体的质量中心上去了，因此，不论定义了多少节点，刚性体仅有六个自由度。作用在钢体上的力合力矩由每个时间步的节点力和力矩合成，然后计算钢性体的运动，在转换到节点位移。在定义该模型时，需要定义的参数包括密度（DENS）、弹性模量（EX）（钢体不要用不切实际的杨氏弹性模量或密度，刚体不能再变硬因为它是完全钢硬的）和波松比（NUXY），以及所有的转动约束参数(相对于整体笛卡尔坐标系，它可以约束任何一个方向、两个方向或三个方向完全被约束)和除Y向以外的