CN105893661A - 一种多工位级进模模具母体结构拓扑优化方法 - Google Patents
一种多工位级进模模具母体结构拓扑优化方法 Download PDFInfo
- Publication number
- CN105893661A CN105893661A CN201610191125.6A CN201610191125A CN105893661A CN 105893661 A CN105893661 A CN 105893661A CN 201610191125 A CN201610191125 A CN 201610191125A CN 105893661 A CN105893661 A CN 105893661A
- Authority
- CN
- China
- Prior art keywords
- optimization
- mould
- optimize
- mold
- topology
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
本发明公开了一种多工位级进模模具母体结构拓扑优化方法,1、根据模具母体以外的凸凹模结构要素,确定母体的优化区域;2、采用有限元分析方法对分离工序进行数值模拟,获得冲压过程中模具的受力;3、基于变密度法对模具母体结构进行拓扑优化,将步骤2中的模具受力作为力的边界条件映射至拓扑优化有限元模型中,以模具结构允许变形为优化约束条件进行拓扑优化设计有限元模拟,选择不同的单元密度值进行比较,根据最大优化效益指标确定优化结果;4、根据优化结果对模具母体进行设计并检验所设计的模具是否满足变形约束条件。本发明能够有效提高模具结构刚度,使冲压件的质量更加稳定。同时,使模具达到轻量化的效果。
Description
技术领域
本发明涉及机械工程塑性加工方法及模具领域,尤其涉及一种多工位级进模模具母体结构拓扑优化方法。
背景技术
模具母体是指工作零件(凸凹模、压边圈等)之外的结构。模具母体占模具总质量的比率可达80%以上,传统的模具母体结构设计主要是凭借经验,缺乏科学理论指导。为了保证模具结构的强度刚度,模具母体结构尺寸取值往往偏大,安全系数较大,但增加模具结构体积并不一定能有效提高冲压质量的稳定性,体积增大也会造成生产成本的提高,增加冲床工作时的能源消耗。
尤其对于高精度多工位级进模,冲压过程中模具结构产生的变形量对冲压件的成形质量会产生重要影响,对于精度要求高的工位,传统模具结构不能保证其具有足够的刚度。较大的模具结构变形将会造成冲压件质量不稳定,反复修模,延长模具的生产制造周期。
发明内容
本发明的目的在于克服上述现有技术的缺点和不足,提供一种多工位级进模模具母体结构拓扑优化方法。提高模具刚度及冲压件冲压质量的稳定性,降低模具结构变形对冲压件质量的影响,减少修模次数,减轻模具重量,缩短模具生产制造周期,降低加工成本。
本发明通过下述技术方案实现:
一种多工位级进模模具母体结构拓扑优化方法,包括如下步骤:
步骤(1):根据模具母体以外的结构要素(包括凸凹模、叉车位、落料槽等),对模具母体的优化区域进行设计,确定模具母体的优化区域;
步骤(2):采用有限元分析方法对冲压工序进行有限元模拟,获得冲压过程中模具与坯料的接触力,其中坯料设置为弹塑性体,凸凹模、压边圈等冲压工具设置为刚性体;
步骤(3):建立多工位级进模模具母体结构拓扑优化有限元模型,定义凸凹模、模具母体等为弹性体,以接触力作用点的线性位移(模具变形量)为响应,以优化区域的体积最小化(最小化是指:在保证模具结构变形量在优化约束范围内的前提下,使优化区域的材料分布更加合理,无多余材料的存在,避免造成材料浪费,获得满足约束条件的最少材料分布,即体积最小化)为优化目标,以结构变形的允许范围为优化约束,以优化区域的单元密度为优化变量,基于变密度法实现模具母体结构的拓扑优化;
步骤(4):根据拓扑优化结果获得优化区域的单元密度云图,通过计算最大优化效益指标确定合适的密度值,根据优化结果设计模具母体结构,并对所设计的模具结构进行变形约束条件检验。
上述步骤(4),通过拓扑优化运算获得模具母体的单元密度云图,根据计算最大优化效益指标选择合理的密度值(通过拓扑优化运算可以获得优化区域的单元密度云图,不同的密度值对应不同的材料分布,通过对比不同密度值对应的优化效益指标选择密度值,优化效益指标越大则表明该密度值对应的材料分布越趋于合理),确定优化结果,参考此时的单元分布完成模具母体结构的设计。
上述步骤(3),建立拓扑优化有限元模型时需定义拓扑优化三要素,包括优化目标、优化变量和优化约束;优化目标是满足优化约束时,优化区域的材料合理分布、实现优化区域的体积最小化;优化变量是优化区域的单元密度;优化约束即为模具结构的允许变形量。
上述步骤(2),采用有限元分析方法(Dynaform)对拉深、弯曲等成形工序进行数值模拟,在*.dyn(后处理文件)文件中增加*DATABASE_NCFORC语句(关键字文件),获得模具与坯料的接触力;基于Deform-3D(有限元分析软件)对冲孔、冲裁等分离工序进行数值模拟,获得冲压过程中模具与坯料的接触力文件(*.KEY),其中坯料设置为弹塑性体,凸凹模、压边圈等冲压工具设置为刚性体;然后通过载荷映射将接触力作为力的边界条件加载至拓扑优化有限元模型中。
为了反映结构减重和结构刚度两种因素之间的相互影响,在选定单元密度值时引入优化效益指标公式,即:
定义结构优化体积比Vfrac=V/V0,其中V0为模具原始结构设计空间的初始体积,V为优化后结构的体积。优化体积比Vfrac则代表了结构减重情况,Vfrac越小表示结构减重效果越好。定义结构优化变形量比Sfrac=S/S0,其中S0为原结构初始变形量,S为优化后结构的变形量,Sfrac越小表示结构变形越小。优化效益指标OPT越大,表明优化所获得的结构越合理。
本发明相对于现有技术,具有如下的优点及效果:
(1)本发明是基于有限元数值模拟方法,可预先对级进模模具结构进行受力分析和拓扑优化设计,能够较为准确预测模具结构的弹性变形量。
(2)本发明采用的多工位级进模模具结构拓扑优化方法,能够有效提高模具结构刚度,减小模具结构弹性变形量,从而使冲压件的质量更加稳定。同时,可以提高材料的利用率,有效减小模具的体积和重量,达到轻量化的效果。
附图说明
图1是级进模模具母体结构拓扑优化设计流程图
图2是某多工位级进模下模可优化区域
图3是不同的单元密度值对应的优化效益值
图4是优化后下模垫脚结构
图5是某多工位级进模上模可优化区域
图6是优化后上模垫脚结构
上图中:1垫脚(可优化区域);2下模模座;3下模垫脚;4上模垫脚可优化区域;5上模模座;6上模垫脚
具体实施方式
下面结合具体实施例对本发明作进一步具体详细描述。
如图1至6所示。本发明公开了一种多工位级进模模具母体结构拓扑优化方法。
图1所示为多工位级进模模具母体结构拓扑优化设计的流程图。
主要分为以下步骤:1、根据模具母体以外的凸凹模等零件、叉车位、落料槽等结构要素,确定母体的优化区域;2、采用有限元分析方法对拉深等成形工序、冲裁等分离工序进行数值模拟,获得冲压过程中模具的受力,其中坯料设置为弹塑性体,凸凹模、压边圈等冲压工具设置为刚性体;3、基于变密度法对模具母体结构进行拓扑优化,设置凸凹模、模具母体等为弹性体,将步骤2中的模具受力作为力的边界条件映射至拓扑优化有限元模型中,以模具结构允许变形为优化约束条件进行拓扑优化设计有限元模拟,选择不同的单元密度值进行比较,根据最大优化效益指标确定优化结果;4、根据优化结果对模具母体进行设计并检验所设计的模具是否满足变形约束条件。本发明能够有效提高模具结构刚度,使冲压件的质量更加稳定。同时,使模具达到轻量化的效果。
以下通过实施例举例说明本发明多工位级进模模具母体结构拓扑优化方法。
1)下模模具母体结构拓扑优化
采用本发明提出的多工位级进模模具母体结构拓扑优化方法,对该级进模下模模具母体结构进行优化设计。首先是对母体可优化空间进行设计,如图2所示,确定母体可优化空间。然后,基于有限元法分别采用Dynaform、Deform-3D对成形工序和分离工序进行数值模拟,定义凸凹模、压边圈等冲压工具为刚性体,坯料为弹塑性体,获得各冲压工序中模具与板料的接触力。然后,定义拓扑优化三要素:优化目标、优化变量、优化约束,以优化区域的体积最小化为优化目标,以优化区域的单元密度为优化变量,以接触力作用点的线性位移(弹性变形量)为优化约束条件,即保证最大线性位移不超过某一限定值。在该实施例中限定点的最大线性位移小于0.1mm。在Hyperworks中建立模具母体结构拓扑优化有限元模型,并将接触力作为边界条件映射至拓扑优化有限元模型上,通过迭代运算获得拓扑优化结果。优化变量(单元密度值)的变化区间是0~1,分别选择单元密度值大于0.1、0.2、0.3、0.4时的单元分布进行对比,并计算其对应的优化效益值,按照单元密度大于0.10、0.20、0.30、0.40设计后结构的优化效益值分别为1.5608、1.5049、1.5322、1.4934,如图3所示。可见单元密度大于0.1时的优化效益指标最佳,因此,下模垫脚结构选取单元密度大于0.10作为参考进行结构设计最为合理。
图4所示为按单元密度大于0.10进行设计的下模母体(即垫脚)结构,具有9条垫脚。该套模具原先已按经验进行设计,需要11条垫脚,即优化前后分别为11和9条垫脚。分别对这两种结构施加所需要的接触力,得到优化前后的最大变形量分别为0.082mm和0.077mm,均满足最大变形小于0.1mm的变形约束条件,而且与优化前相比,优化后变形量减小了6.1%,重量减轻了13.2%。
2)上模模具母体结构拓扑优化
首先是对上模母体可优化空间进行设计,由于是上模,不存在落料槽及叉车位,因此只要根据凸凹模的安装需求,将不可优化区域确定并去除之后,便是模具母体可优化空间,如图5所示。然后按照上一实施例的步骤进行拓扑优化(优化约束仍为限定点的最大线性位移小于0.1mm),对比单元密度值大于0.10、0.20、0.30、0.40时的单元分布,并同样通过对优化效益指标进行计算,选取优化效益指标最佳的单元密度大于0.10的情况作为结构设计的依据。
图6所示为按照单元密度大于0.10进行设计的下模母体(即垫脚)结构,具有8条垫脚。该套模具原先已按经验进行设计,需要11条垫脚,即优化前后分别为11和8条垫脚。分别对这两种结构施加所需要的接触力,模拟计算得到优化前后的最大变形量基本相等,均满足变形约束条件,但与优化前相比,优化后的重量减轻了26.6%。
如上所述,便可较好地实现本发明。
本发明的实施方式并不受上述实施例的限制,其他任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (5)
1.一种多工位级进模模具母体结构拓扑优化方法,其特征在于包括如下步骤:
步骤(1):根据模具母体以外的结构要素,对模具母体的优化区域进行设计,确定模具母体的优化区域;
步骤(2):采用有限元分析方法对冲压工序进行有限元模拟,获得冲压过程中模具与坯料的接触力,其中坯料设置为弹塑性体,凸凹模、压边圈设置为刚性体;
步骤(3):建立多工位级进模模具母体结构拓扑优化有限元模型,定义凸凹模、模具母体为弹性体,以接触力作用点的线性位移为响应,以优化区域的体积为优化目标,以结构变形的允许范围为优化约束,以优化区域的单元密度为优化变量,基于变密度法实现模具母体结构的拓扑优化;
步骤(4):根据拓扑优化结果获得优化区域的单元密度云图,通过计算最大优化效益指标确定合适的密度值,根据优化结果设计模具母体结构,并对所设计的模具结构进行变形约束条件检验。
2.根据权利要求1所述多工位级进模模具母体结构拓扑优化方法,其特征在于:步骤(4),通过拓扑优化运算获得模具母体的单元密度云图,根据计算最大优化效益指标选择密度值,确定优化结果,参考此时的单元分布完成模具母体结构的设计。
3.根据权利要求1所述多工位级进模模具母体结构拓扑优化方法,其特征在于:步骤(3),建立拓扑优化有限元模型时需定义拓扑优化三要素,包括优化目标、优化变量和优化约束;优化目标是满足优化约束时,优化区域的材料合理分布、实现优化区域的体积最小化;优化变量是优化区域的单元密度;优化约束即为模具结构的允许变形量。
4.根据权利要求1所述多工位级进模模具母体结构拓扑优化方法,其特征在于:步骤(2),采用有限元软件Dynaform、Deform-3D分别对成形及分离工序进行数值模拟,通过在*.dyn文件中增加*DATABASE_NCFORC语句及*.KEY文件,获得模具与坯料的接触力;其中坯料设置为弹塑性体,凸凹模、压边圈设置为刚性体;然后通过载荷映射将接触力作为力的边界条件加载至拓扑优化有限元模型中。
5.根据权利要求2所述多工位级进模模具母体结构拓扑优化方法,其特征在于:
为了反映结构减重和结构刚度两种因素之间的相互影响,在选定单元密度值时引入优化效益指标公式,即:
定义结构优化体积比Vfrac=V/V0,其中V0为模具原始结构设计空间的初始体积,V为优化后结构的体积;优化体积比Vfrac则代表了结构减重情况,Vfrac越小表示结构减重效果越好;定义结构优化变形量比Sfrac=S/S0,其中S0为原结构初始变形量,S为优化后结构的变形量,Sfrac越小表示结构变形越小;优化效益指标OPT越大,表明优化所获得的结构越合理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610191125.6A CN105893661A (zh) | 2016-03-30 | 2016-03-30 | 一种多工位级进模模具母体结构拓扑优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610191125.6A CN105893661A (zh) | 2016-03-30 | 2016-03-30 | 一种多工位级进模模具母体结构拓扑优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105893661A true CN105893661A (zh) | 2016-08-24 |
Family
ID=57014062
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610191125.6A Pending CN105893661A (zh) | 2016-03-30 | 2016-03-30 | 一种多工位级进模模具母体结构拓扑优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105893661A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108009381A (zh) * | 2017-12-25 | 2018-05-08 | 北京航空航天大学 | 一种位移和全局应力混合约束下的连续体结构可靠性拓扑优化方法 |
CN108171731A (zh) * | 2017-09-28 | 2018-06-15 | 中国矿业大学(北京) | 一种顾及拓扑几何多要素约束的最小影像集自动优选方法 |
CN108256215A (zh) * | 2018-01-15 | 2018-07-06 | 广东省智能制造研究所 | 一种基于结构拓扑优化的龙门加工中心滑枕及其优化设计方法 |
CN108629086A (zh) * | 2018-04-12 | 2018-10-09 | 华中科技大学 | 一种适用于增材制造零件的应力调整方法 |
CN109583078A (zh) * | 2018-11-27 | 2019-04-05 | 上海交通大学 | 接触力自动加载方法 |
CN111054829A (zh) * | 2019-11-11 | 2020-04-24 | 泊头市金键模具有限责任公司 | 一种冲压模具的模面细化方法 |
CN111539145A (zh) * | 2020-04-22 | 2020-08-14 | 安徽工程大学 | 一种汽车衣帽架模压成型模具优化方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN204135191U (zh) * | 2014-09-28 | 2015-02-04 | 苏州工业职业技术学院 | 一种支撑件多工位级进模具 |
-
2016
- 2016-03-30 CN CN201610191125.6A patent/CN105893661A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN204135191U (zh) * | 2014-09-28 | 2015-02-04 | 苏州工业职业技术学院 | 一种支撑件多工位级进模具 |
Non-Patent Citations (5)
Title |
---|
XIA QIN-XIANG等: "Numerical Simulation and Experimental Research on Multi-position Progressive Stamping Process of Automotive Structural Part", 《APPLIED MECHANICS AND MATERIALS》 * |
夏琴香等: "某高强板汽车结构件多工位级进模结构分析", 《模具工业》 * |
汤禹成: "高强度钢板料拉深成形工艺稳健设计及模具结构拓扑优化", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
白金泽: "《LS-DYNA3D理论基础与实例分析》", 31 May 2005 * |
郭鸣骥: "多工位级进模模具母体结构优化设计关键技术研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108171731A (zh) * | 2017-09-28 | 2018-06-15 | 中国矿业大学(北京) | 一种顾及拓扑几何多要素约束的最小影像集自动优选方法 |
CN108171731B (zh) * | 2017-09-28 | 2021-12-14 | 中国矿业大学(北京) | 一种顾及拓扑几何多要素约束的最小影像集自动优选方法 |
CN108009381A (zh) * | 2017-12-25 | 2018-05-08 | 北京航空航天大学 | 一种位移和全局应力混合约束下的连续体结构可靠性拓扑优化方法 |
CN108009381B (zh) * | 2017-12-25 | 2021-05-25 | 北京航空航天大学 | 一种位移和全局应力混合约束下的连续体结构可靠性拓扑优化方法 |
CN108256215A (zh) * | 2018-01-15 | 2018-07-06 | 广东省智能制造研究所 | 一种基于结构拓扑优化的龙门加工中心滑枕及其优化设计方法 |
CN108629086A (zh) * | 2018-04-12 | 2018-10-09 | 华中科技大学 | 一种适用于增材制造零件的应力调整方法 |
CN108629086B (zh) * | 2018-04-12 | 2020-06-02 | 华中科技大学 | 一种适用于增材制造零件的应力调整方法 |
CN109583078A (zh) * | 2018-11-27 | 2019-04-05 | 上海交通大学 | 接触力自动加载方法 |
CN111054829A (zh) * | 2019-11-11 | 2020-04-24 | 泊头市金键模具有限责任公司 | 一种冲压模具的模面细化方法 |
CN111054829B (zh) * | 2019-11-11 | 2022-07-22 | 泊头市金键模具有限责任公司 | 一种冲压模具的模面细化方法 |
CN111539145A (zh) * | 2020-04-22 | 2020-08-14 | 安徽工程大学 | 一种汽车衣帽架模压成型模具优化方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105893661A (zh) | 一种多工位级进模模具母体结构拓扑优化方法 | |
CN101339574B (zh) | 基于回弹补偿的混凝土搅拌叶片模具型面设计***及方法 | |
CN103752685B (zh) | 多工位高效高精度汽车覆盖件的冲压方法 | |
CN111177906B (zh) | 一种离散化模具型面精确补偿方法 | |
CN103433380B (zh) | 一种汽车覆盖件拉延模具的凸成形方法 | |
CN104200037B (zh) | 一种小双曲度钣金件的成形模具设计方法 | |
CN102096748B (zh) | 一种轻量化大刚度精冲压力机机身的优化设计方法 | |
CN105183932A (zh) | 一种考虑模具变形的冲压计算机仿真计算方法 | |
CN104615809B (zh) | 采用逆向反求因子的回弹补偿方法 | |
CN104077439A (zh) | 一种新型高强钢轮辐拉深冲孔复合工序的数值模拟方法 | |
CN101811156A (zh) | 基于cae分析法获取板材成型回弹量的方法 | |
CN104698969A (zh) | 一种基于拟合法的回弹补偿方法 | |
CN105528484A (zh) | 工件成型力的cae仿真分析方法以及成型装置优化方法 | |
CN104809304A (zh) | 一种基于变间隙压边的铝板冲压成形工艺优化方法 | |
CN104156546A (zh) | 基于t样条的汽车覆盖件模具的形面再设计方法 | |
CN103577635A (zh) | 一种叶片型面数据的拟合方法 | |
CN106980742B (zh) | 向有限元仿真模型引入冲压成型信息的载荷投影映射方法 | |
CN105426629B (zh) | 一种汽车覆盖件拉延模具结构强度计算方法 | |
CN101944143A (zh) | 双层圆挤压筒过盈量参数化设计方法 | |
CN107025354A (zh) | 一种基于极差分析的车窗升降板成形工艺优化方法 | |
CN111014415A (zh) | 一种基于cae技术的拉延筋制造方法 | |
CN102722619B (zh) | 一种确定冲压形成汽车覆盖件的零件的材料利用率的方法 | |
CN103406449B (zh) | 一种汽车覆盖件拉延模具 | |
CN101789031B (zh) | 基于压力中心优化的级进模带料设计*** | |
CN203459528U (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160824 |
|
WD01 | Invention patent application deemed withdrawn after publication |