并行计算时半功倍油缸

并行计算，“时”半功倍

并行计算，“时”半功倍 2011： ——并行计算在大型汽车覆盖件冲压成形仿真分析中价值凸显由于有限元本身的缺陷，CAE仿真软件在应用中难以处理计算精度和速度这一对矛盾，而本文探讨的动态显式有限元冲压成形并行仿真系统的关键技术则可以很好地解决这一工程问题，具有重要的实际应用价值。

在汽车和模具工业中，CAE技术已被广泛应用，为提高汽车和模具的设计水平，缩短新产品的开发周期发挥着越来越重要的作用。在模具及冲压领域，板料冲压成形 CAE技术的应用可以大大提高模具精度，缩短模具新产品的开发周期，有效提高冲压件质量，在模具和汽车新产品开发中应用越来越多，应用水平也不断提高。

有限元无法克服的矛盾

目前国内、外的有限元仿真软件都是基于弹塑性有限变形理论开发的，按算法的不同可分为隐式和显式两大类，这两类软件各有特点。隐式算法要求解联立方程组，常常会遇到难以收敛的难题，但由于接触力的计算采用了拉格朗日乘子法，可以精确满足边界条件，因此计算结果的应力场比较精确，对于回弹的计算比较有利；显式算法无须求解联立方程组，因此无条件收敛，但受稳定性限制，有最小时间步长的要求，计算中常常要花费较多时间，为了减少计算时间，一般采取加大冲压速度和板料质量密度的方法。尽管显式软件有一些先天性缺陷，但是由于它没有收敛性问题，在工程中应用越来越广泛。

使用有限元仿真软件模拟实际冲压成形过程，要想得到满意的结果，必须处理好计算精度和速度这一对矛盾，而这恰恰是仿真软件所无法克服的，其原因是：

（1）通常板料划分的网格越密，有限元仿真计算的精度就越高，但是花费的计算时间就越长，特别是对于像整体侧围这样的大型汽车覆盖件，为了使计算时间不至于太长，必须控制单元总数，导致对冲压件形状拟合的精度降低；

（2）为了减少计算时间，人为地加大冲压速度，由此也带来了不真实的惯性效应，降低了计算结果的可靠性；

（3）拉延筋是冲压模具上非常重要的工艺结构，它的尺寸较小、曲率较大，要想精确地描述它，就需要很密的网格单元，因此有限元仿真软件大都采用等效拉延筋模型，但这种模型忽略了沿拉延筋方向材料间的不均匀挤压作用，另外对于压边圈闭合时不能模拟出拉延筋的成型效果，在压料面形状起伏较大时失真较大，因而不能非常正确地反映实际情况。

针对以上不足，各冲压成形CAE软件公司纷纷开发并行计算系统，如LS-DYNA、PAM-STAMP等，这些并行计算系统的开发大大提高了软件的工程适用性，有效解决了工程分析中的矛盾。

动态显式并行计算系统的关键技术

这里重点介绍冲压成形仿真分析中常用的动态显式并行计算系统的关键技术。

1.单元格式

采用4节点壳单元，在每个单元的中面定义随动坐标系。该坐标系与单元一起运动，所有物理量在随动坐标系中定义。经等参插值，可以求得得单元中的速度为：V = NT Vi，其中Vi为节点速度矢量，N为插值函数。

2.控制方程

板料成型系统必须满足Hu-Washizu变分原理，系统的控制方程为：Mü = F ext + F c - F int，其中M为质量矩阵，ü为节点加速度矢量，F ext表示节点所受外力矢量，F c表示节点接触力矢量，F int表示节点内力矢量。在显式算法中，由中心差分法可代入得到位移矢量U的方程式。3.本构方程

材料的塑性变形采用Hill的各向异性弹塑性准则。

4.接触摩擦力的计算

接触力和摩擦力的计算精度直接影响到板料变形的计算精度。接触力的计算有两种基本方法：罚函数法、拉格朗日乘子法。罚函数法是一种近似方法，它允许相互接触的边界产生穿透，并通过罚因子将接触力大小与边界穿透量大小联系起来。这种方法比较简单，适合于显式算法，罚因子的选取会影响计算结果的可靠性。拉格朗日乘子法不允许接触边界的相互穿透，是一种精确的接触力算法，但它与显式算法不相容，要求特殊的数值处理。

接触搜寻采用级域法，它针对多物体接触体系提出，避免了主从面法处理多物体接触的不便，并提高了接触搜寻效率。

板料与模具间的摩擦力采用非经典的非线性摩擦定律计算，弥补了库仑摩擦定律没有考虑摩擦力与相对滑移量之间关系的不足。

5.有限元并行化方法

通常动态显式有限元并行仿真系统是基于消息传递和区域分解算法开发的，将板料根据计算中要用的PE数分区，每个PE处理一个分区的单元计算，由于每个分区板料的相邻边界上是有作用力关系的，因此相邻分区对应的PE间在每一计算步都要进行数据通信。

（1）板料初始化分区。 并行计算执行时间Tp 的计算公式为 Tp = rsTs + Tpe + Tc，式中Ts为在一个PE上的执行时间，rs为用一个PE执行时不可并行化部分所耗时间占Ts的百分比， Tpe为并行部分的执行时间，Tc为通信所消耗的时间。通常Tpe要取各PE中最长的一个。

由上式可知，要想降低并行计算时间Tp，初始化分区必须考虑两个方面：各PE的计算任务尽可能均衡

毛毡厂

润眼仪工厂

SKY天瑞精工

成都环境试验中心

打桩机多少钱一台