CN105183932A - 一种考虑模具变形的冲压计算机仿真计算方法 - Google Patents
一种考虑模具变形的冲压计算机仿真计算方法 Download PDFInfo
- Publication number
- CN105183932A CN105183932A CN201510413733.2A CN201510413733A CN105183932A CN 105183932 A CN105183932 A CN 105183932A CN 201510413733 A CN201510413733 A CN 201510413733A CN 105183932 A CN105183932 A CN 105183932A
- Authority
- CN
- China
- Prior art keywords
- mold
- mould
- die
- finite element
- computer simulation
- 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.)
- Pending
Links
Landscapes
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
本发明提出了一种考虑模具变形的冲压计算机仿真计算方法,包括以下步骤:步骤1:构建模具模型:步骤2:对步骤1构建的模具模型进行网格划分;步骤3:在有限元软件Dynaform中对模具的板料进行网格划分,设置冲压过程参数与模具的材料属性参数;步骤4:利用有限元前处理软件hypermesh对凸模、凹模、压料板的运动状态进行设置,并且对模具进行弹性体化设置,完成对冲压过程的有限元仿真,得到模具的应力应变结果以及板料成型结果。该方法相比于现有技术,在进行板料成形数值仿真时将模具设为弹塑性变形体,考虑了模具变形对板料成形数值仿真的影响,建立精度更高、更加适合先进高强钢板冲压数值模拟的分析方法。
Description
技术领域
本发明书属于汽车覆盖件冲压模具技术领域,特别涉及一种考虑模具变形的冲压计算机仿真计算方法。
背景技术
随着现代汽车的发展,人们对汽车安全性和燃油经济性有了新的要求。先进高强度钢板具有较高的强度和延伸率,适用于汽车工业关于轻量化、碰撞安全性的要求,在车身中应用更加广泛。传统普通板材的冲压CAE仿真计算中,由于板料的成型力对模具的变形影响很小,模具一般被视为刚性体模型也可以获得较高的计算精度。然而先进高强钢板由于屈服强度高,使得所需的成型力增大、模具受力增加,由此产生的模具变形对先进高强钢板成型的影响不容忽视。针对先进高强钢板的冲压数值模拟分析,由于数值模拟精度降低,使得数值模拟对冲压工艺的指导优化作用下降,增加了模具调试工作量,降低了生产效率。
通过对国内外公开发表的文献和专利进行检索,现有技术中相关的专利有“冲压模具结构分析数值模拟方法”(专利号:200710040000.4),该专利提供一种冲压模具结构分析数值模拟方法,侧重于模具的结构分析。该专利先用刚体化模具进行板料成形数值分析,得到板料与模具之间的接触应力,通过一定的载荷映射算法,将接触应力映射到模具体上,进行冲压模具结构分析,得到模具的应力应变结果。这种技术存在如下缺点:(1)没有考虑板料成形过程中模具变形对板料成形的影响;(2)在进行板料数值模拟分析时没有将模具视为弹塑性变形体。
发明内容
本发明为了克服上述现有技术的传统钢板冲压将模具体为刚性体的分析方法中,由于忽略了模具变形对成形性的影响,对于高强钢冲压CAE精度不高的问题,提供了一种考虑模具变形的先进高强钢板冲压计算机仿真方法,通过将模具体视为弹塑性体,在冲压过程中考虑模具变形对先进高强钢板成形性的影响,建立精度更高、更加适合先进高强钢板冲压数值模拟的分析方法。
一种考虑模具变形的冲压计算机仿真计算方法,包括以下步骤:
步骤1:构建模具模型:
基于模具受力对模具结构进行简化,并利用模具变形产生的实际拉延间隙值偏置拉延型面,用偏置后的型面裁剪上模,建立具有实际凸凹模间隙的模具模型;
步骤2:对步骤1构建的模具模型进行网格划分;
对模具模型进行有限元三维网格划分,并将划分后的模具导入有限元软件Dynaform中;
其中,网格尺寸范围为2-10mm;
步骤3:在有限元软件Dynaform中对模具的板料进行网格划分,板料单元模型采用BT壳单元,厚度方向积分点为7个,设置冲压过程参数与模具的材料属性参数;
步骤4:利用有限元前处理软件hypermesh对凸模、凹模、压料板的运动状态进行设置,并且对模具进行弹性体化设置,完成对冲压过程的有限元仿真,得到模具的应力应变结果以及板料成型结果。
所述步骤1中,利用模具变形产生的实际拉延间隙值偏置拉延型面,所述拉延间隙按照1.1倍的板料料厚取值
所述步骤2中对模具进行有限元网格划分时,采用局部网格加密技术保证细微型面网格单元的符型性,保证凸凹模圆角有三个以上的单元。
所述步骤4中凸模、凹模、压料板的运动状态的控制实现过程如下:
凸模的状态,通过限制凸模最底层单元6个自由度实现;
凹模与压料板的状态通过限制最上层单元x、y方向的平动自由度实现。
所述模具的弹性体化设置过程通过修改Ls-dyna关键字完成,具体如下:
在KEY文件中增加体单元控制关键字:*CONTROL_SOLID;增加弹塑性模型关键字*MAT_ISOTROPIC_ELASTIC_PLASTIC,替换刚体模型关键字Rigid-body;定义关键字*BOUNDARY_PRESCRIBED_MOTION_SET来控制弹塑性体工具的运动。
有益效果
本发明提出了一种考虑模具变形的冲压计算机仿真计算方法,包括以下步骤:步骤1:构建模具模型:步骤2:对步骤1构建的模具模型进行网格划分;步骤3:在有限元软件Dynaform中对模具的板料进行网格划分,板料单元模型采用BT壳单元,厚度方向积分点为7个,设置冲压过程参数与模具的材料属性参数;步骤4:利用有限元前处理软件hypermesh对凸模、凹模、压料板的运动状态进行设置,并且对模具进行弹性体化设置,完成对冲压过程的有限元仿真,得到模具的应力应变结果以及板料成型结果。
该方法相比于现有技术,在进行板料成形数值仿真时将模具设为弹塑性变形体,考虑了模具变形对板料成形数值仿真的影响,建立精度更高、更加适合先进高强钢板冲压数值模拟的分析方法,其优点具体如下:
(1)建立更接近实际冲压条件的仿真计算模型,得到更加精确的考虑模具变形的板料成型仿真结果,为先进高强钢板冲压成型的工艺设计、模具设计提供更具价值的参考依据;(2)得到更加精确的板料应力场结果,提高了零件回弹的仿真计算精度;(3)获得模具的应力应变状态,为探究模具受力提供参考。
附图说明
图1为本发明所述方法的流程示意图;
图2为模具爆炸图;
图3为弹性工具体板料成型图;
图4为刚性工具板料成型图;
图5为实际冲压零件图。
具体实施方式
下面将结合实例与附图做进一步的说明。
如图1所示,一种考虑模具变形的冲压计算机仿真计算方法,包括以下步骤:
步骤1:构建模具模型:
基于模具受力对模具结构进行简化,并利用模具变形产生的实际拉延间隙值偏置拉延型面,用偏置后的型面裁剪上模,建立具有实际凸凹模间隙的模具模型;
步骤2:对步骤1构建的模具模型进行网格划分;
对模具模型进行有限元三维网格划分,并将划分后的模具导入有限元软件Dynaform中;
其中,网格尺寸范围为2-10mm;
步骤3:在有限元软件Dynaform中对模具的板料进行网格划分,板料单元模型采用BT壳单元,厚度方向积分点为7个,设置冲压过程参数与模具的材料属性参数;
步骤4:利用有限元前处理软件hypermesh对凸模、凹模、压料板的运动状态进行设置,并且对模具进行弹性体化设置,完成对冲压过程的有限元仿真,得到模具的应力应变结果以及板料成型结果。
以某DP780先进高强钢前纵梁零件及其拉延模具为例,板料料厚为1mm。
简化对模具受力变形影响较小的部位,如部分标准件及其安装台,让空以及直径小于50mm的孔洞。按实际拉延间隙值偏置拉延型面,拉延间隙按照1.1倍的板料料厚取值,将型面偏置1.1mm。
用偏置后的型面裁剪上模,建立实际凸凹模间隙。划分模具体有限元三维网格,网格大小为10mm,最小为2mm,采用局部网格加密技术保证细微型面网格单元的符型性,保证凸凹模圆角有三个以上的单元,并导入Dynaform有限元软件。
图2为模具爆炸图。在Dynafrom软件中划分板料网格,板料单元模型采用BT壳单元,厚度方向积分点为7个,定义拉延筋系数、压边力等冲压工艺参数;定义模具体运动载荷边界条件,其中凹模与压料板的运动速度为1000mm/s,冲压过程分为三个阶段:1、压料阶段:压料板与凹模下行20mm,凹模与压边圈闭合;2、压料板闭合阶段:压料板、凹模、压边圈下行30mm,压料板到底。3、模具闭合阶段:凹模与压边圈继续下行30mm,拉延结束。
输出前处理产生的Ls-dyna关键字文件导入hypermesh有限元前处理软件,定义弹塑性体模具体,并提交计算,该关键字文件包括模型网格信息、运动参数、材料属性等所有信息。
计算得到模具的应力应变结果以及板料冲压成型结果。凸模、凹模、压料板最大等效应力分别为286Mpa、296Mpa、253Mpa,凸模、凹模、压料板Z向弹性变形分别为0.25mm、0.1mm、0.08mm。
图3与图4中的英文分别表示板料成型结果,crack表示破裂区,riskofcrack表示即将破裂区域,safe表示安全区域即成型结果良好,wrinkletendency表示具有起皱趋势的区域,wrinkle表示拉延起皱区域,severewrinkle表示严重起皱区域,insufficientstretch表示拉延不充分区域。
通过图3、图4、图5相同区域的拉延结果对比可知,将工具体设置为弹塑性之后得到的成型结果更接近实际零件的冲压结果,因此数值模拟具有更高的精度。
Claims (5)
1.一种考虑模具变形的冲压计算机仿真计算方法,其特征在于,包括以下步骤:
步骤1:构建模具模型:
基于模具受力对模具结构进行简化,并利用模具变形产生的实际拉延间隙值偏置拉延型面,用偏置后的型面裁剪上模,建立具有实际凸凹模间隙的模具模型;
步骤2:对步骤1构建的模具模型进行网格划分;
对模具模型进行有限元三维网格划分,并将划分后的模具导入有限元软件Dynaform中;
其中,网格尺寸范围为2-10mm;
步骤3:在有限元软件Dynaform中对模具的板料进行网格划分,板料单元模型采用BT壳单元,厚度方向积分点为7个,设置冲压过程参数与模具的材料属性参数;
步骤4:利用有限元前处理软件hypermesh对凸模、凹模、压料板的运动状态进行设置,并且对模具进行弹性体化设置,完成对冲压过程的有限元仿真,得到模具的应力应变结果以及板料成型结果。
2.根据权利要求1所述的一种考虑模具变形的冲压计算机仿真计算方法,其特征在于,所述步骤1中,利用模具变形产生的实际拉延间隙值偏置拉延型面,所述拉延间隙按照1.1倍的板料料厚取值
3.根据权利要求1所述的一种考虑模具变形的冲压计算机仿真计算方法,其特征在于,所述步骤2中对模具进行有限元网格划分时,采用局部网格加密技术保证细微型面网格单元的符型性,保证凸凹模圆角有三个以上的单元。
4.根据权利要求1所述的一种考虑模具变形的冲压计算机仿真计算方法,其特征在于,所述步骤4中凸模、凹模、压料板的运动状态的控制实现过程如下:
凸模的状态,通过限制凸模最底层单元6个自由度实现;
凹模与压料板的状态通过限制最上层单元x、y方向的平动自由度实现。
5.根据权利要求1-4任一项所述的一种考虑模具变形的冲压计算机仿真计算方法,其特征在于,所述模具的弹性体化设置过程通过修改Ls-dyna关键字完成,具体如下:
在KEY文件中增加体单元控制关键字:*CONTROL_SOLID;增加弹塑性模型关键字*MAT_ISOTROPIC_ELASTIC_PLASTIC,替换刚体模型关键字Rigid-body;定义关键字*BOUNDARY_PRESCRIBED_MOTION_SET来控制弹塑性体工具的运动。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510413733.2A CN105183932A (zh) | 2015-07-15 | 2015-07-15 | 一种考虑模具变形的冲压计算机仿真计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510413733.2A CN105183932A (zh) | 2015-07-15 | 2015-07-15 | 一种考虑模具变形的冲压计算机仿真计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105183932A true CN105183932A (zh) | 2015-12-23 |
Family
ID=54906012
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510413733.2A Pending CN105183932A (zh) | 2015-07-15 | 2015-07-15 | 一种考虑模具变形的冲压计算机仿真计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105183932A (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106175946A (zh) * | 2016-08-29 | 2016-12-07 | 北京缔佳医疗器械有限公司 | 一种牙齿矫治器的生产方法及装置 |
CN106547999A (zh) * | 2016-12-08 | 2017-03-29 | 湖南大学 | 基于显式动力学有限元方法的拉延模结构分析 |
CN106980742A (zh) * | 2017-05-09 | 2017-07-25 | 湖南天汽模汽车模具技术股份有限公司 | 一种向有限元仿真模型引入冲压成型信息的载荷投影映射方法 |
CN107025354A (zh) * | 2017-04-15 | 2017-08-08 | 贵州大学 | 一种基于极差分析的车窗升降板成形工艺优化方法 |
CN107092745A (zh) * | 2017-04-15 | 2017-08-25 | 贵州大学 | 一种基于方差分析的车窗升降板成形工艺优化方法 |
CN109598071A (zh) * | 2018-12-07 | 2019-04-09 | 厦门理工学院 | 一种半球体凸模胀形变形确定方法 |
CN109657304A (zh) * | 2018-12-03 | 2019-04-19 | 珠海格力电器股份有限公司 | 一种空调随机振动自动化仿真方法及系统 |
CN109773058A (zh) * | 2019-01-09 | 2019-05-21 | 柳州福臻车体实业有限公司 | 冲压模具符型方法 |
CN110263384A (zh) * | 2019-05-28 | 2019-09-20 | 南京理工大学 | 基于Laplacian微分域变形的三维网格曲面变厚度偏置造型方法 |
CN111238782A (zh) * | 2018-11-09 | 2020-06-05 | 华晨宝马汽车有限公司 | 用于检测模具的方法、设备、系统和可存储介质 |
CN111421043A (zh) * | 2020-03-31 | 2020-07-17 | 中国第一汽车股份有限公司 | 一种针对拉延模具压料圈部件弹性变形的补偿方法 |
CN112560319A (zh) * | 2020-12-22 | 2021-03-26 | 上海交通大学 | 一种疲劳断裂凹坑缺陷有限元模型构建方法 |
CN114547879A (zh) * | 2022-02-18 | 2022-05-27 | 福州大学 | 基于冲压成型反演算法的球罐壳板下料方法 |
CN118426904A (zh) * | 2024-03-19 | 2024-08-02 | 深圳未来技术软件有限公司 | 一种步长自适应控制的板料冲压成形仿真方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101050960A (zh) * | 2007-04-26 | 2007-10-10 | 上海交通大学 | 冲压模具结构分析数值模拟方法 |
-
2015
- 2015-07-15 CN CN201510413733.2A patent/CN105183932A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101050960A (zh) * | 2007-04-26 | 2007-10-10 | 上海交通大学 | 冲压模具结构分析数值模拟方法 |
Non-Patent Citations (3)
Title |
---|
于磊: ""考虑模具弹性变形的先进高强钢冲压回弹控制与结构拓扑优化"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
刘帅辰: ""先进高强钢冲压精确 CAE 分析及回弹补偿研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
蒋浩民,等: ""先进高强度钢板的冲压成形特性及其应用"", 《塑性工程学报》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106175946A (zh) * | 2016-08-29 | 2016-12-07 | 北京缔佳医疗器械有限公司 | 一种牙齿矫治器的生产方法及装置 |
CN106547999B (zh) * | 2016-12-08 | 2020-05-12 | 湖南大学 | 基于显式动力学有限元方法的拉延模结构分析方法 |
CN106547999A (zh) * | 2016-12-08 | 2017-03-29 | 湖南大学 | 基于显式动力学有限元方法的拉延模结构分析 |
CN107025354A (zh) * | 2017-04-15 | 2017-08-08 | 贵州大学 | 一种基于极差分析的车窗升降板成形工艺优化方法 |
CN107092745A (zh) * | 2017-04-15 | 2017-08-25 | 贵州大学 | 一种基于方差分析的车窗升降板成形工艺优化方法 |
CN106980742A (zh) * | 2017-05-09 | 2017-07-25 | 湖南天汽模汽车模具技术股份有限公司 | 一种向有限元仿真模型引入冲压成型信息的载荷投影映射方法 |
CN106980742B (zh) * | 2017-05-09 | 2020-04-24 | 湖南天汽模汽车模具技术股份有限公司 | 向有限元仿真模型引入冲压成型信息的载荷投影映射方法 |
CN111238782B (zh) * | 2018-11-09 | 2023-03-07 | 华晨宝马汽车有限公司 | 用于检测模具的方法、设备、系统和可存储介质 |
CN111238782A (zh) * | 2018-11-09 | 2020-06-05 | 华晨宝马汽车有限公司 | 用于检测模具的方法、设备、系统和可存储介质 |
CN109657304A (zh) * | 2018-12-03 | 2019-04-19 | 珠海格力电器股份有限公司 | 一种空调随机振动自动化仿真方法及系统 |
CN109657304B (zh) * | 2018-12-03 | 2021-06-18 | 珠海格力电器股份有限公司 | 一种空调随机振动自动化仿真方法及系统 |
CN109598071A (zh) * | 2018-12-07 | 2019-04-09 | 厦门理工学院 | 一种半球体凸模胀形变形确定方法 |
CN109598071B (zh) * | 2018-12-07 | 2022-10-11 | 厦门理工学院 | 一种半球体凸模胀形变形确定方法 |
CN109773058B (zh) * | 2019-01-09 | 2021-04-27 | 柳州福臻车体实业有限公司 | 冲压模具符型方法 |
CN109773058A (zh) * | 2019-01-09 | 2019-05-21 | 柳州福臻车体实业有限公司 | 冲压模具符型方法 |
CN110263384B (zh) * | 2019-05-28 | 2020-02-21 | 南京理工大学 | 基于Laplacian微分域变形的三维网格曲面变厚度偏置造型方法 |
CN110263384A (zh) * | 2019-05-28 | 2019-09-20 | 南京理工大学 | 基于Laplacian微分域变形的三维网格曲面变厚度偏置造型方法 |
CN111421043A (zh) * | 2020-03-31 | 2020-07-17 | 中国第一汽车股份有限公司 | 一种针对拉延模具压料圈部件弹性变形的补偿方法 |
CN111421043B (zh) * | 2020-03-31 | 2021-07-13 | 中国第一汽车股份有限公司 | 一种针对拉延模具压料圈部件弹性变形的补偿方法 |
CN112560319A (zh) * | 2020-12-22 | 2021-03-26 | 上海交通大学 | 一种疲劳断裂凹坑缺陷有限元模型构建方法 |
CN114547879A (zh) * | 2022-02-18 | 2022-05-27 | 福州大学 | 基于冲压成型反演算法的球罐壳板下料方法 |
CN118426904A (zh) * | 2024-03-19 | 2024-08-02 | 深圳未来技术软件有限公司 | 一种步长自适应控制的板料冲压成形仿真方法 |
CN118426904B (zh) * | 2024-03-19 | 2024-10-11 | 深圳未来技术软件有限公司 | 一种步长自适应控制的板料冲压成形仿真方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105183932A (zh) | 一种考虑模具变形的冲压计算机仿真计算方法 | |
CN104077439B (zh) | 一种新型高强钢轮辐拉深冲孔复合工序的数值模拟方法 | |
CN103480704B (zh) | 一种帽形件渐进折弯成形方法 | |
CN106547999B (zh) | 基于显式动力学有限元方法的拉延模结构分析方法 | |
CN103433380B (zh) | 一种汽车覆盖件拉延模具的凸成形方法 | |
CN104765912A (zh) | 一种铝板冲压工艺的稳健性优化方法 | |
CN101274341B (zh) | 一种降低塑性回弹的金属板材加工成形方法 | |
CN104200037A (zh) | 一种小双曲度钣金件的成形模具设计方法 | |
Do et al. | Incremental forming of 3D structured aluminum sheet | |
CN105893661A (zh) | 一种多工位级进模模具母体结构拓扑优化方法 | |
CN105426629A (zh) | 一种汽车覆盖件拉延模具结构强度计算方法 | |
CN103870641A (zh) | 一种变厚板(vrb)成形极限应力图的构建和使用方法 | |
CN103406449B (zh) | 一种汽车覆盖件拉延模具 | |
CN101976291B (zh) | 一种热交换器板片的制作方法 | |
Ahmetoglu et al. | Computer simulation for tool and process design in sheet forming | |
CN110765693A (zh) | 基于cae仿真技术的贵金属钽电容器外壳优化模具设计方法 | |
Younis et al. | Numerical and experimental investigation of parameters affect the forming load during rubber pad sheet metal forming | |
CN102672051A (zh) | 一种开关叠形弹片折弯模具 | |
JP5839323B2 (ja) | 金属板体の金型寸法自動探索システム | |
KR102599529B1 (ko) | 복합 드로잉을 이용한 테일게이트 익스텐션 제조방법 | |
CN202683775U (zh) | 一种开关叠形弹片折弯模具 | |
Younis et al. | Rubber Pad Sheet Metal Forming of Round Metal Blanks into Multi Shape Axisymmetric Cups by FEA and Experimental Methods | |
CN101539952B (zh) | 一种用于成形模具与成形装备的耦合刚度测试分析的工作方法 | |
CN102847849A (zh) | 一种解决汽车门外板的门把手周边变形的工艺方法 | |
Prakash et al. | Investigation and analysis for the wrinkling behaviour of deep drawn die sheet metal component by using fast form |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20151223 |
|
WD01 | Invention patent application deemed withdrawn after publication |