CN101482403B - 一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法 - Google Patents
一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法 Download PDFInfo
- Publication number
- CN101482403B CN101482403B CN200910042686XA CN200910042686A CN101482403B CN 101482403 B CN101482403 B CN 101482403B CN 200910042686X A CN200910042686X A CN 200910042686XA CN 200910042686 A CN200910042686 A CN 200910042686A CN 101482403 B CN101482403 B CN 101482403B
- Authority
- CN
- China
- Prior art keywords
- resilience
- finite element
- grid
- cloud
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
本发明为一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法。该方法首先把点云数据转移到几何模型对应位置。其次,把点云数据转化为三角形或四边形网格,并输出为有限元输入格式文件,得到板料的网格;并按实际模具型面建立仿真模型,得到模具型面网格。然后,根据模具和板料网格,建立仿真模型,把回弹后的网格反向压回到回弹前的位置。最后,把有限元中的节点输出,就可得到回弹前这些测量点云的位置。采用该发明可快速准确的确定回弹补偿型面,或确定回弹后零件上的修边线。
Description
技术领域
本发明涉及一套薄板冲压成型领域,具体为一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法。
背景技术
汽车车身上主要是薄板冲压件。为了提高车身的安全性,减少车身重量,在汽车车身零件上采用越来越多高强度钢,而高强度钢的冲压回弹问题特别突出。所谓回弹,是指零件从冲压模具退出时,冲压件中弹性应力的释放,造成零件的变形。回弹后的零件与模具的几何型面不再吻合,如果按照零件几何模型加工型面,回弹导致冲压后零件的尺寸精度达不到设计尺寸要求。
为了克服回弹问题,需要准确知道各个部位的回弹量,在对应部位修改型面,这样使回弹后的零件刚好达到设计尺寸要求。问题的难点在于如何准确的得到零件各个部位的回弹量。由于有限元仿真精度的限制,对回弹不能准确计算。目前通过仿真方法,还不能准确得到冲压件的回弹量。通过对成型后零件的三维数据扫描,可以得到零件回弹后的形状。与零件设计几何模型对上后,就可以知道各个部位的回弹量,从而对回弹进行补偿。
三坐标测量设备是发现零件尺寸故障的一种有效手段,目前有各种扫描点云的方法,具有比较高的精度,能够检测加工表面的加工精度。通过各工序发生回弹后零件的测量,确定产生回弹的工序,采取措施克服回弹问题。一般是拉延序的回弹补偿不足所导致的。需要根据测量点云数据,对模型做进一步的补偿。一般在点云数据与几何模型对上后,根据点与相应部位的差值,乘一个系数,对零件型面修改,得到补偿型面。由于点数很多,该过程时间长,并且精度没有保证。
本发明是一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法。通过该方法,可根据回弹后的点云与数模的关系,得到回弹补偿面上的点。根据这些点,可快速准确的确定出回弹补偿型面。该方法提供了一种得到精确补偿型面的方法,在工业上有广泛的应用价值。
发明内容
本发明的原理如图1所示,把零件点云转化为三角形或四边形单元,输出回弹后的有限元单元和节点数据,通过模具压到其回弹前的形状,得到回弹前的单元和节点数据。这些节点数据就是回弹前的点云数据。整个过程相当于采用仿真方法把零件压到其在模具退出前的状态,也相当于模拟零件二次冲压。
本发明处理点云的流程如图2所示。首先把点云数据转移到几何模型对应位置。其次,把点云数据转化为三角形或四边形网格,并输出为有限元输入格式文件,得到板料的网格;并按实际模具型面建立仿真模型,得到模具型面网格。然后,根据模具和板料网格,建立仿真模型,把回弹后的网格反向压回到回弹前的位置。最后,把有限元中的节点输出,就可得到回弹前这些测量点云的位置。计算出的回弹前点云数据,可以保持为有限元格式,也可保持为供几何建模软件采用的其它格式。
点云数据既可以是三维测量设备直接或取的数据,也可以是经过点云预处理的数据。点云预处理如去掉噪音点、去掉重复点、光滑处理、点云转化为三角面片、处理边界等。计算出回弹前点云数据,可以保持为有限元格式,也可保持为供几何建模软件采用的其它格式。
采用该发明得到基于测量点云的回弹前后的网格,按适合有限元的回弹补偿方法,得到零件回弹补偿面上的网格,可进一步输出三角面片或数据点,根据这些数据点可对模具重新建模,得到回弹补偿后的型面。
采用该发明得到基于测量点云的回弹前后的网格。通过回弹前位置的修边线,采用网格映射方法可以得到回弹后零件修边线的位置。
以上两个都是高强度钢薄板成型中的关键问题。
附图说明
图1本发明实施原理图
图2本发明处理点云步骤
图3本发明实施实施算例中点云数据与数模的对应
图4本发明实施实施算例中点云的三角化
图5本发明实施实施算例中有限元仿真模型
图6本发明实施算例中得到回弹前后的网格
图7本发明仿真验证中回弹前后对应的网格
图8本发明仿真验证中建立的有限元模型
图9本发明仿真验证中计算的回弹前后点的位置
图10本发明输出为有限元程序能够读取的格式示例
具体实施方式
采用本发明提出的点云数据处理方法,开发了点云数据简化程序、点云数据与数模对应的程序,以及点云三角化程序等,并采用LSDYNA软件以及DYNAFORM软件来完成算例工作。
首先,将点云数据简化,每3mm左右取1个点。把简化后的点云数据与零件几何数模对上。对上的方法如下:在数模中找到三点,如零件上具有明显的几何特征点,如空圆心,凸台,棱角边缘等,这些点在点云中的位置也是可以确定的。根据其中对应的一点,就可把点云放到几何模型附近,然后,通过一条边的重合,零件位置可进一步确定,通过第三个点,就可把点云数据与几何模型数据对上。其数学处理是对点云的平移和旋转。点云与几何模型的对应关系见图3所示。
其次,把点云数据三角化,如图4所示。所谓三角化,就是由点云数据做顶点,把相关的数据点,用三角形一块一块的连接起来,这些三角形要没有重叠,并且三角形的质量要好,三角形的各边长度比较均匀。
点云三角化的目的与传统的三角化是不同的。传统的三角化是为了形成三角形面片,供几何建模软件反求建立精确的几何模型。这里是为了形成供有限元分析用的三角形网格。这些网格由单元、节点构成。
并输出为有限元程序能够读取的格式,对有限元软件LSDYNA的输入文件,需要按其输入文件格式和关键词编写数据,如*NODE,*ELEMENT_SHELL等。对节点数据,由节点号或节点坐标构成。节点号为整数,占8个字节,单元坐标为浮点数,占16个字节。对单元数据,每个单元所在的部件(part),每个面单元一般由4个节点构成,对三角形,第三个节点与第四个节点为同一节点。单元中数值都为整数,各为8个字节。
然后,把三角化后的零件转入到冲压仿真模型,调整好间隙,板料和模具不要出现穿透,把摩擦系数定义为很小,如:0.005,这样做是为了避免板料顺利进入模腔,不改变零件的局部形状。定义运动和接触边界条件,不要采用自适应网格。建立的零件二次冲压仿真模型如图5所示。首先凹模运动,与压边块接触后,向下移动,在行程模与凸模闭合。
采用有限元方法对采用点云三角化网格的零件冲压仿真模型求解,得到冲压零件前后的位置如图6所示。有限元输出结果中含有单元组成、节点坐标。这些节点就是测量点在回弹前所在的位置。
方法的验证
上面方法中,通过仿真计算的网格是否与回弹前的网格有准确对应关系是最关键的问题,如果仿真计算得到的网格与回弹前的网格有明显的差别,这说明该种方法不可用。在冲压仿真和回弹仿真后,模型中具有了回弹后的网格,这里把回弹后的网格假设为试验数据测量得到的。采用上面的方法,把仿真计算得到的网格与实际回弹前的网格相比较,就可以验证本发明提出的方法。 在零件拉深成型仿真后,继续进行回弹仿真。回弹前后的网格见图7所示(因为显示原因,后面都只取局部)。
在冲压仿真模型中把板料网格更换为回弹后的网格,调整凹凸模位置,可建立零件二次冲压仿真模型,如图8所示。在该模型中,把回弹后的网格压至回弹前的网格。
零件二次冲压仿真模型计算的回弹前网格与实际回弹前的网格的对比见图9所示。从图中可看出,网格基本重合,在法向基本没有误差,切向有少量偏移,最大偏差量在1mm以内。这证明了上节提出方法是可行的。
Claims (3)
1.一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法,其特征在于:首先把点云数据转移到几何模型对应位置;其次,把点云数据转化为三角形或四边形网格,并输出为有限元输入格式文件,得到板料的网格;并按实际模具型面建立第一仿真模型,得到模具型面网格;然后,根据模具型面网格和板料的网格,建立第二仿真模型,把回弹后的网格反向压回到回弹前的位置;最后,把有限元中的节点输出,就得到回弹前这些测量点云的位置;计算出的回弹前点云数据,保持为有限元格式,或保持为供几何建模软件采用的其它格式。
2.根据权利要求1所述的一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法,其特征在于通过建立薄板冲压第二仿真模型,对转化为有限元格式的点云数据进行冲压仿真模拟,得到这些点云数据在回弹前的位置。
3.根据权利要求1所述的一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法,其特征在于采用的点云数据是三维测量设备直接获取的数据,或是经过点云预处理的数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910042686XA CN101482403B (zh) | 2009-02-19 | 2009-02-19 | 一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910042686XA CN101482403B (zh) | 2009-02-19 | 2009-02-19 | 一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101482403A CN101482403A (zh) | 2009-07-15 |
CN101482403B true CN101482403B (zh) | 2011-06-29 |
Family
ID=40879614
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910042686XA Expired - Fee Related CN101482403B (zh) | 2009-02-19 | 2009-02-19 | 一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101482403B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105335571A (zh) * | 2015-11-24 | 2016-02-17 | 大连楼兰科技股份有限公司 | 一种碰撞仿真虚拟现实的环境搭建方法 |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101853316B (zh) * | 2010-04-20 | 2012-12-05 | 湖南大学 | 一种薄板冲压仿真中的接触状况的显示方法 |
CN101825456B (zh) * | 2010-05-11 | 2012-01-04 | 青岛再特模具有限公司 | 预变形模具中有效确定补偿量的方法 |
CN102672059B (zh) * | 2011-03-16 | 2015-09-09 | 重庆长安汽车股份有限公司 | 根据仿真冲压工件厚度确定出模具凹凸模修改型面的方法 |
CN102672060A (zh) * | 2011-03-16 | 2012-09-19 | 重庆长安汽车股份有限公司 | 一种基于板料网格的翻边修边线确定方法 |
CN102354320A (zh) * | 2011-07-19 | 2012-02-15 | 奇瑞汽车股份有限公司 | 一种基于白光扫描的梁类零件匹配方法 |
CN103128147B (zh) * | 2011-11-25 | 2015-07-08 | 北汽福田汽车股份有限公司 | 梁类模具的回弹修正加工方法和系统 |
JP6044191B2 (ja) * | 2012-08-31 | 2016-12-14 | Jfeスチール株式会社 | プレス成形品のスプリングバック抑制対策方法および解析装置 |
US9921572B2 (en) * | 2013-11-12 | 2018-03-20 | Embraer S.A. | Springback compensation in formed sheet metal parts |
CN104268349A (zh) * | 2014-09-30 | 2015-01-07 | 华南理工大学 | 一种复杂曲面下翻边修边线的精确控制方法 |
CN105335999B (zh) * | 2015-11-24 | 2019-04-19 | 大连楼兰科技股份有限公司 | 街区物体三维扫描生成物体有限元数字模型的方法 |
CN105783857B (zh) * | 2016-04-21 | 2018-06-22 | 哈尔滨工业大学 | 一种汽轮机转子挠度测量方法 |
CN105975786A (zh) * | 2016-05-13 | 2016-09-28 | 桥弘数控科技(上海)有限公司 | 一种折弯机加工部件工进深度的计算方法 |
CN106311906B (zh) * | 2016-11-01 | 2018-11-16 | 天津职业技术师范大学 | 汽车覆盖件模具及有限元计算点云装置 |
KR102326768B1 (ko) | 2018-01-30 | 2021-11-15 | 제이에프이 스틸 가부시키가이샤 | 프레스 성형 방법, 강성 향상 위치 특정 방법, 및 프레스 성형 시스템 |
CN111655391B (zh) * | 2018-01-30 | 2022-02-25 | 杰富意钢铁株式会社 | 模压成型方法、刚性提高位置确定方法、模压成型系统及模压成型品 |
CN109101736A (zh) * | 2018-08-22 | 2018-12-28 | 上海博汇模具有限公司 | 一种用于高张力板件的整体补偿方法 |
CN110814086B (zh) * | 2019-10-16 | 2021-01-01 | 唐山钢铁集团有限责任公司 | 一种板材冲压成形后回弹量的测量方法 |
CN110823689A (zh) * | 2019-10-22 | 2020-02-21 | 首钢集团有限公司 | 一种板料的摩擦系数及回弹量的测量方法 |
CN111145236B (zh) * | 2019-12-04 | 2022-10-18 | 东南大学 | 一种基于数字孪生的产品拟实物装配模型生成方法及实现框架 |
CN113221402B (zh) * | 2021-04-23 | 2023-06-02 | 湖北文理学院 | 一种冲压件回弹预测及监控方法、系统及存储介质 |
CN114429065B (zh) * | 2022-01-07 | 2024-02-23 | 西北工业大学 | 基于有限元的微铣削过程后刀面接触材料回弹标定方法 |
-
2009
- 2009-02-19 CN CN200910042686XA patent/CN101482403B/zh not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105335571A (zh) * | 2015-11-24 | 2016-02-17 | 大连楼兰科技股份有限公司 | 一种碰撞仿真虚拟现实的环境搭建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101482403A (zh) | 2009-07-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101482403B (zh) | 一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法 | |
CN107116139B (zh) | 模具型面的设计方法及覆盖件模具 | |
CN111274671B (zh) | 一种基于数字孪生的复杂产品装配过程精准修配方法及其运行系统 | |
CN104289748B (zh) | 一种大型薄壁蒙皮自适应等壁厚铣削系统及其加工方法 | |
JP4886527B2 (ja) | 金型モデルデータの修正方法 | |
CN101546349B (zh) | 一种基于曲面边缘线或截面线网格映射的回弹补偿面生成方法 | |
US7542889B2 (en) | Determination of a model of a geometry of a metal sheet forming stage | |
CN103617302A (zh) | 一种汽车覆盖件拉延模具精确型面构建的方法 | |
CN103191958B (zh) | 一种板料成形回弹的光学检测方法 | |
CN101339574A (zh) | 基于回弹补偿的混凝土搅拌叶片模具型面设计系统及方法 | |
CN103433380B (zh) | 一种汽车覆盖件拉延模具的凸成形方法 | |
CN106021761B (zh) | 一种汽车覆盖件回弹评测方法 | |
CN101949687A (zh) | 基于视觉测量的汽车车门的检测方法 | |
CN110814086B (zh) | 一种板材冲压成形后回弹量的测量方法 | |
CN110706335B (zh) | 一种凸模胀形工件回弹后三维建模方法及装置 | |
CN103727865A (zh) | 一种用于检测汽车门板的检具 | |
CN111291450B (zh) | 一种确定发动机罩回弹三维模型的方法、装置及设备 | |
CN108664689A (zh) | 用于生成模拟模型的方法 | |
CN104525668B (zh) | 一种多点成形加工曲面的冲头调形方法 | |
CN103810343B (zh) | 一种零件关键工艺表面的创成式精细化表征方法 | |
CN106294921A (zh) | 一种复杂自由曲面的抛光方法 | |
CN101968343A (zh) | 基于机器人单目视觉测量的汽车车门的检测方法 | |
CN104036085A (zh) | 用于cae仿真分析的复杂型面汽车声学部件的厚度统计方法 | |
CN109992912A (zh) | 一种基于变异系数法的最优回弹补偿系数确定方法 | |
CN109108168B (zh) | 冲压模具合模间隙的计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110629 Termination date: 20130219 |