一种提高管材弯曲成型过程仿真效率的方法及仿真方法技术

本发明专利技术属于机械加工技术领域，公开了一种提高管材弯曲成型过程仿真效率的方法及仿真方法，该方法对用于描述弯曲成型过程的运动方程的时间域进行当量变换，将实际的弯曲时间域[0, T]变换为标准的时间域[0，1]，空间域保持不变；同时用时间域变换的变换系数对管材密度、初始速度和初始角速度进行变换，获得标准时间域下的管材密度、初始速度与初始角速度；再基于变换后获得的等价运动方程和变换后管材密度、初始速度和初始角速度，对管材弯曲成型过程进行仿真求解。本发明专利技术方法可缩短仿真计算的时间，可以适用于管材弯曲成型模具与工艺的设计中，具有很好的工程应用前景。

【技术实现步骤摘要】
一种提高管材弯曲成型过程仿真效率的方法及仿真方法
本专利技术属于机械加工
，涉及弯管成型过程仿真技术，关于一种关于管材成型仿真的计算方法，特别指一种弯管绕弯成型过程仿真分析与提高仿真效率的方法。
技术介绍
以管材为毛坯，通过塑性加工手段，制造管材零件的工艺称为管材塑性加工。管材塑性加工，由于易于实现产品的轻量化、强韧化和满足低耗高效、精确制造等方面的要求，并且能获得形状复杂的制品，具有精确稳定、无废料绿色环保等优点。管材塑性成型得到的弯管零件具有轻量化、结构稳定、高刚度等优良的力学性能，而成为塑性加工技术发展的一个重要方向。管材弯曲成型技术是管材塑性加工的重要组成部分。用管材制造的弯曲件，具有重量轻、吸震力强、介质流通量大等一系列优良性能，广泛应用于普通气体、液体的输送管路及汽车、航空航天、化工和其它高新
在管材的弯曲成型过程中，容易出现回弹、起皱、截面畸变及减薄拉裂等问题。弯曲成型过程中，管材的外侧壁受拉而减薄甚至拉裂，内侧壁受压而增厚起皱，内、外侧受力不匀导致截面畸变，卸载后又存在回弹现象。这些问题一直是管材弯曲成型加工行业未能有效解决的技术难题，也是当今国内外塑性加工学科研究的难点和热点。采用计算机辅助工程分析，即CAE技术，对管材弯曲成型过程进行仿真分析可以有效预测弯曲成型过程中出现的各种缺陷，是目前流行的设计手段与分析工具。对管材弯曲成型问题采用计算机辅助工程分析的CAE技术研究弯曲模具结构及弯曲成型工艺，涉及到管材的几何大变形、材料的弹塑性变形及管材与弯曲模具的接触等问题。这是一个高度非线性的仿真模拟过程，成型过程仿真分析的计算量非常大。为缩短仿真时间、提高分析效率，工程上采用一些特殊的方法来加速模拟过程，如在板料冲压成形的模拟中，常用速度加速法，通过人为提高冲头的运动速度，即将几十毫米/秒的实际冲压速度提高到3000mm/s左右(称为虚拟冲压速度)，来缩短弯曲过程仿真的时间。又如质量缩放法，即将某些特定单元的质量密度放大来加大时间步长，但是缩放系数的确定需要试探。这种单方面调整某一物理量的方法尽管提高了计算效率但容易产生误差，且不够严谨、缺乏理论基础的支持。
技术实现思路
本专利技术目的是提出一种管材绕弯成型过程仿真分析与提高仿真效率的计算方法。该方法可以对管材弯曲成型过程进行仿真计算，并且缩短了仿真计算的时间，因此可以有效应用于管材弯曲成型模具与工艺的设计中，具有很好的工程应用前景。基于上述目的，本专利技术首先提出一种提高管材弯曲成型过程仿真的效率的方法，该方法对用于描述弯曲成型过程的运动方程的时间域进行当量变换，将实际的弯曲时间域[0,T]变换为标准的时间域[0，1]，空间域保持不变；同时用时间域变换的变换系数对管材密度、初始速度和初始角速度进行变换，获得标准时间域下的管材密度、初始速度与初始角速度；再基于变换后获得的等价运动方程和变换后管材密度、初始速度和初始角速度，对管材弯曲成型过程进行仿真求解。具体地说，上述方法采用变换方程t1＝αt对运动方程进行变换，其中取α＝1/T。管材密度的变换公式为ρ1＝ρ/T2，初始速度变换公式为v1＝v0T，初始角速度变换公式为ω1＝ω0T＝v0T/R。作为改进，为了使求解递推公式的临界时间步长时间域变换后与变换前保持不变，在仿真求解前，还需要选择管材或弯曲模具的某一零件，如芯棒，对其质量密度进行密度缩放，密度缩放公式为ρi＝ρ1/α2＝ρ1T2。基于上述方法，本专利技术进一步提出一种基于当量变换的管材弯曲成型过程仿真方法，包括如下步骤：步骤一，将用于描述管材弯曲成型过程的运动方程的时间域进行当量变换，将实际的弯曲时间域[0,T]变换为标准的时间域[0，1]，空间域保持不变；步骤二，用步骤1时间域变换中的变换系数对管材密度、初始速度、初始角速度进行等价变换，获得标准时间域下的管材密度、初始速度和初始角速度；步骤三，选择管材或弯曲模具的某一零件，如芯棒，对其质量密度进行缩放，使求解递推公式的临界时间步长与时间域变换前保持不变；步骤四，基于步骤1-3变换后的标准时间域下的等价运动方程、管材密度、初始速度、初始角速度，对管材弯曲成型过程进行仿真求解。上述方法中，步骤1采用变换方程t1＝αt对运动方程进行变换，其中取α＝1/T。步骤2中管材密度的变换公式为ρ1＝ρ/T2，初始速度变换公式为v1＝v0T，初始角速度变换公式为ω1＝ω0T＝v0T/R。步骤3中采用密度缩放公式ρi＝ρ1/α2＝ρ1T2对所选择的管材或弯曲模具的某一零件的密度进行缩放。下面分步对本专利技术方法做出详细的解释说明。步骤1：弯曲成型基本方程变换。管材弯曲成型问题属于结构动力学问题。结构动力学问题的控制方程可以通过运动方程描述，它是借助牛顿(Newton)第二定律这一基本力学原理建立的数学表达式，描述了结构中力与位移(包括速度和加速度)的关系。结构动力学问题的运动方程一般形式：σij,j-ρui,tt-μui,t+fi＝0t∈[0,T],i＝x,y,z∈Ω(1)其中，σij是应力张量，ui是位移矢量，ρ是材料密度，μ是阻尼系数，T是弯曲成型所需的时间，ui,t是位移矢量对时间的一阶导数即速度，ui,tt是位移矢量对时间的二阶导数即加速度，μui,t是粘性力，ρui,tt是惯性力。方程(1)边界条件中的初始条件为：t＝0,u＝u0,v＝v0(2)其中，u0是初始位移，v0是初始速度。方程(1)的求解一般使用数值方法，如采用有限元方法对方程(1)进行离散求解。由于方程(1)包含空间坐标变量与时间坐标变量，因此有限元离散求解的做法是：空间域网格划分建立有限元模型，时间域采用差分格式建立递推公式。本专利技术的做法是对公式(1)所描述的弯曲成型基本方程的时间域进行变换，即时间变量的变换表达式如下：t1＝αt(3)其中，α为变换因子，又称缩放系数。这样将实际弯曲成型过程的时间域[0,T]变换到了当量仿真时间域[0，T1]＝[0，αT]。很显然，若取α＝1/T，则有T1＝1，即将实际弯曲成型时间T变换到当量仿真时间1。α也可取其它小于1的值，这样可变换到更小的时间域。称当量仿真时间域[0，1]为标准时间域。时间域变换的同时对弯曲问题(1)中的空间变量(x,y,z)不作变换，这样(1)式等价为初始条件(2)也变换为：t1＝0,u1＝u0,v1＝v0/α(5)这样通过(3)式的时间域变换将(1)式所描述的管材弯曲成型问题等价变换为(4)描述的不同时间域的管材弯曲成型问题。很显然问题(1)与问题(4)的管材及模具的尺寸与拓扑关系一致，因此有限元网格模型也不变，由于变换而发生改变的是管材及模具的材料参数与边界条件中的初始条件。这样的变换方便工程实际应用。步骤2：参数计算。经过第一步的时间变换后，用于管材弯曲仿真的主要参数，即材料密度与初始速度可以从原始数据变换计算得到(取α＝1/T)：1)材料密度：ρ1＝ρ/T2(6)2)初始速度：v1＝v0T(7)使用本方法时须注意各个物理量的量纲，即量纲要统一。如采用国际单位制SI，基本量纲为：毫米mm，千克kg，秒s，则力的量纲为N(牛顿)、材料密度量纲为kg/mm3，速度的量纲为mm/s，弹性模量E与应力σ的量纲为MPa(N/mm2)。很显然，变换因子α是无量纲的，本文档来自技高网...

【技术保护点】
一种提高管材弯曲成型过程仿真效率的方法，其特征在于，该方法对用于描述弯曲成型过程的运动方程的时间域进行当量变换，将实际的弯曲时间域[0,T]变换为标准的时间域[0，1]，空间域保持不变；同时用时间域变换的变换系数对管材密度、初始速度、初始角速度进行变换，获得标准时间域下的管材密度、初始速度和始角速度，再基于变换后获得的等价运动方程和变换后管材密度、初始速度与初始角速度，对管材弯曲成型过程进行仿真求解。

【技术特征摘要】
1.一种提高管材弯曲成型过程仿真效率的方法，其特征在于，该方法对用于描述弯曲成型过程的运动方程的时间域进行当量变换，将实际的弯曲时间域[0,T]变换为标准的时间域[0，1]，空间域保持不变；同时用时间域变换的变换系数对管材密度、初始速度和初始角速度进行变换，获得标准时间域下的管材密度、初始速度和始角速度，再基于变换后获得的等价运动方程和变换后管材密度、初始速度与初始角速度，对管材弯曲成型过程进行仿真求解。2.根据权利要求1所述的提高管材弯曲成型过程仿真效率的方法，其特征在于，选择芯棒，对其进行密度缩放。3.一种基于当量变换的管材弯曲成型过程仿真方法，其特征在于，包括如下步骤：步骤1，将用于...

【专利技术属性】
技术研发人员：刘红，蒋兰芳，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：浙江;33