CN110457755A - 一种轨道车辆压型件毛坯快速预示方法 - Google Patents
一种轨道车辆压型件毛坯快速预示方法 Download PDFInfo
- Publication number
- CN110457755A CN110457755A CN201910627696.3A CN201910627696A CN110457755A CN 110457755 A CN110457755 A CN 110457755A CN 201910627696 A CN201910627696 A CN 201910627696A CN 110457755 A CN110457755 A CN 110457755A
- Authority
- CN
- China
- Prior art keywords
- molded piece
- grid
- blank
- rail vehicle
- node
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000000605 extraction Methods 0.000 claims abstract description 13
- 230000007935 neutral effect Effects 0.000 claims abstract description 10
- 230000002093 peripheral effect Effects 0.000 claims abstract description 10
- 230000037237 body shape Effects 0.000 claims abstract description 4
- 230000009471 action Effects 0.000 claims description 4
- 230000007547 defect Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 239000002184 metal Substances 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000003796 beauty Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004836 empirical method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Abstract
本发明公开了压型件技术领域的一种轨道车辆压型件毛坯快速预示方法,旨在解决现有技术中轨道车辆压型件毛坯预示依靠工人经验、开发周期长、成本高的技术问题,对压型件的模型抽取能反应其主体形状的某一个表面作为分析对象;设置板料壁厚以及板料属性;对抽取的压型件的模型表面进行网格划分,并通过等距偏置算法获得模型的中性层网格;设置分析参数,建立有限元平衡方程,获得预示的轮廓线的外围节点的信息;依据外围节点的信息绘制样条线,得到展开后的毛坯预示边界线。本发明所述轨道车辆压型件毛坯快速预示方法,能够快速、准确地对压型件进行毛坯预示,减少了试模次数,大大缩短了开发周期,降低了生产成本。
Description
技术领域
本发明属于压型件技术领域,具体涉及一种轨道车辆压型件毛坯快速预示方法。
背景技术
随着轨道列车运营时间和速度的不断提高,对车辆运行安全性和舒适度提出了更高要求。而轨道车辆的结构件大多由压型件构成,压型件具有刚度大、表面平整、自重轻、外观美观、制造简单等优点。压型件的成形质量好坏对列车的稳定性和安全性有很大的影响。钣金零件毛坯的精确预示对控制零件的成形精度和制造成本有着重要意义,因此,在成形之前获得钣金零件毛坯的精确预示是保证成形质量的重要组成部分。
毛坯展开时的分析参数以及毛坯材料和零件模型的网格等都影响着毛坯预示结果的精确度,只有各种参数匹配适当时,才能得到精确的毛坯预示结果。在轨道车辆钣金零件毛坯展开方面,目前多采用经验法和实验修正相结合的方法获得毛坯展开的近似几何形状。使用经验法对钣金零件毛坯预示存在着完全依靠工人经验、开发周期长和成本高等缺点。
发明内容
本发明的目的在于提供一种轨道车辆压型件毛坯快速预示方法,以解决现有技术中轨道车辆压型件毛坯预示依靠工人经验、开发周期长、成本高的技术问题。
为达到上述目的,本发明所采用的技术方案是:一种轨道车辆压型件毛坯快速预示方法,包括以下步骤:
a、对压型件的模型抽取能反应其主体形状的某一个面作为分析对象;
b、设置板料壁厚以及板料属性;
c、对抽取的压型件的模型面进行网格划分,若抽取为内表面或外表面时需通过等距偏置算法获得模型的中性层网格,若抽取面为中性面,无需偏置;
d、设置分析参数,建立有限元平衡方程,获得预示的轮廓线的外围节点的信息;
e、依据外围节点的信息绘制样条线,得到展开后的毛坯预示边界线。
所述步骤d包括以下步骤:
da、采用几何反向变形有限元法,不考虑板料的毛坯网格为初始解的求解方法,给一步逆成形有限元法快速准确的提供一个较为符合要求的初始解;
db、根据初始解的毛坯网格、板料属性以及分析参数,得到每个网格节点的塑性应变和应力,计算得到作用在节点上的内力和外力;
dc、根据每个网格节点的外力以及预设的摩擦系数得到每个节点的摩擦力;
dd、建立有限元平衡方程;
de、以预设的收敛准则求解有限元平衡方程,对初始解进行修正,得到最终考虑成形过程和板料属性的毛坯网格。
所述收敛准则如下:
当某一迭代步计算出所有压型件网格节点残余力向量范数小于上一步相应值,但二者差值大于预定值,则进行下一步迭代,直到二者差值小于预定值,认为迭代收敛;
当某一迭代步计算出所有压型件节点残余力向量范数大于上一步相应值或迭代步数超过所设定的最大迭代步,认为迭代发散。
所述网格划分的方法是:
首先选定需要划分的零件,即压型件,并赋予网格划分的参数,包括单元类型、单元大小、最小单元大小、弦高;
再进行几何清理与几何简化工作,包括针对一些有不重要的孔洞或缺陷的曲面进行补孔操作;
最后进行网格编辑,生成网格模型。
所述等距偏置算法是:
首先计算每一个单元的法线与成形方向的单元夹角θ,当cosθ>0,保持其现有法线方向,当cosθ<0,将该单元的法线方向翻转180°;
再设定外表面或内表面的节点坐标为(x,y,z),偏置后的节点坐标为(x1,y1,z1),通过坐标变换,即原坐标分量沿法线方向移动一半厚度值,将外表面或内表面的节点偏置到中性层上。
所述板料壁厚为板料的原始厚度。
所述板料属性包括杨氏模量、泊松比、屈服强度、密度。
与现有技术相比,本发明所达到的有益效果:本发明所述轨道车辆压型件毛坯快速预示方法,能够快速、准确地对压型件进行毛坯预示,减少了试模次数,大大缩短了开发周期,降低了生产成本。
附图说明
图1是本发明实施例提供的一种轨道车辆压型件毛坯快速预示方法的流程示意图;
图2是本发明实施例提供的一种轨道车辆压型件毛坯快速预示方法的压型件的模型示意图;
图3是通过本发明实施例提供的一种轨道车辆压型件毛坯快速预示方法得到的压型件网格划分示意图;
图4是通过本发明实施例提供的一种轨道车辆压型件毛坯快速预示方法得到的压型件的展开线示意图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1所示,一种轨道车辆压型件毛坯快速预示方法主要包括以下步骤:
a、对压型件的模型抽取能反应其主体形状的某一个面作为分析对象;对于实体模型抽取面可以是外表面,也可以是内表面,对于曲面模型还可以是中性面;
b、设置板料壁厚以及板料属性;
c、对抽取的压型件的模型面进行网格划分,若抽取为内表面或外表面时需通过等距偏置算法获得模型的中性层网格,若抽取面为中性面,无需偏置;
d、设置分析参数,建立有限元平衡方程,获得预示的轮廓线的外围节点的信息;
e、依据外围节点的信息绘制样条线,得到展开后的毛坯预示边界线。
在本实施例中,选取的压型件结构形状如图2所示,应用三维软件建立压型件三维模型,将三维模型导入到毛坯快速预示分析软件中,此压型件模型为有一定厚度的曲面,对压型件的模型提取外表面。
设置压型件的板料参数及板料初始厚度,本实施例的压型件厚度为2mm,材料为Q345钢,弹性模量为204GPa,泊松比为0.332,密度为7830kg/m3,强化系数为794.33MPa,硬化指数为0.32,屈服强度340MPa左右。
网格划分首先选定需要划分的零件,即压型件,并赋予网格划分的参数,包括单元类型、单元大小、最小单元大小、弦高;再进行几何清理与几何简化工作,包括针对一些有不重要的孔洞或缺陷的曲面进行补孔操作;最后进行网格编辑,生成网格模型。本实施例考虑到毛坯预测的精度以及毛坯预测的时长,设置压型件的网格尺寸为3mm,压型件网格为四边形与三角形混合网格,最小单元大小、弦高采用软件默认值,压型件网格模型示意图如图3所示,压型件模型网格总数为7663。本实例选用的模型为压型件的外表面模型,通过等距偏置算法获得零件的中性层网格。首先进行压型件单元法线的调整,具体方法为计算每一个单元的法线与成形方向的单元夹角θ与cosθ,当cosθ>0,保持其现有法线方向,当cosθ<0,将该单元的法线方向翻转180°,最终所有单元的法线方向一致(向内或向外);设外表面的节点坐标为(x,y,z),偏置后的节点坐标为(x1,y1,z1),通过坐标变换,即原坐标分量沿法线方向移动一半厚度值,将外表面或内表面的节点偏置到中性层面上,具体可通过如下公式计算:
其中,y方向和z方向也可用上式进行求解,仅需将x1和x分别替换为y1和y,z1和z,上式中nx、ny、nz为节点法线方向分量,d为板料初始厚度。
设定分析参数,包括压边力、最大迭代步、摩擦系数等参数,然后将网格模型与分析参数提交到求解器,进行预示求解计算,生成预示的轮廓线的点的信息文件,本实施例设置压边力为2000N,最大迭代步为100步,摩擦系数为0.175。
一个与结果接近的初始解可以保证计算收敛,提高计算效率,在此提出一种考虑变形过程的几何反向变形有限元法,即与成形过程相反的逆过程,不考虑材料塑性变形的因素,将网格节点根据变形过程一一映射到初始坯料上,形成初始坯料网格,将此坯料网格作为初始解进行有限元求解。
根据初始解网格,板料参数以及预设的压边力,得到每个网格节点的塑性应变和应力,计算得到作用在节点上的内力和外力。根据每个网格节点的外力以及预设的摩擦系数得到每个节点的摩擦力。根据虚功原理,建立初始毛坯网格所有节点内力、外力和摩擦力所做的功之和为0的有限元平衡方程。
以预设的收敛准则求解有限元平衡方程,对初始解进行修正,得到最终考虑成形过程和板料参数的毛坯网格,获得预示的轮廓线的外围节点的信息。
本实施例选用收敛准则如下:当某一迭代步计算出所有压型件网格节点残余力向量范数小于上一步相应值,但二者差值大于预定值,则进行下一步迭代,直到二者差值小于预定值,可认为迭代收敛,若某一迭代步计算出所有压型件网格节点残余力向量范数大于上一步相应值或迭代步数超过所设定的最大迭代步,可认为迭代发散。当迭代失败时认为毛坯预示不满足要求,需要调整工艺参数。当计算收敛时,可得到最终的毛坯预示边界线。
依据外围节点的信息绘制样条线,得到展开后的毛坯预示边界线,如图4所示。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (7)
1.一种轨道车辆压型件毛坯快速预示方法,其特征是,包括以下步骤:
a、对压型件的模型抽取能反应其主体形状的某一个面作为分析对象;
b、设置板料壁厚以及板料属性;
c、对抽取的压型件的模型面进行网格划分,若抽取为内表面或外表面时需通过等距偏置算法获得模型的中性层网格,若抽取面为中性面,无需偏置;
d、设置分析参数,建立有限元平衡方程,获得预示的轮廓线的外围节点的信息;
e、依据外围节点的信息绘制样条线,得到展开后的毛坯预示边界线。
2.根据权利要求1所述的轨道车辆压型件毛坯快速预示方法,其特征是,所述步骤d包括以下步骤:
da、采用几何反向变形有限元法,不考虑板料的毛坯网格为初始解的求解方法,给一步逆成形有限元法快速准确的提供一个较为符合要求的初始解;
db、根据初始解的毛坯网格、板料属性以及分析参数,得到每个网格节点的塑性应变和应力,计算得到作用在节点上的内力和外力;
dc、根据每个网格节点的外力以及预设的摩擦系数得到每个节点的摩擦力;
dd、建立有限元平衡方程;
de、以预设的收敛准则求解有限元平衡方程,对初始解进行修正,得到最终考虑成形过程和板料属性的毛坯网格。
3.根据权利要求2所述的轨道车辆压型件毛坯快速预示方法,其特征是,所述收敛准则如下:
当某一迭代步计算出所有压型件网格节点残余力向量范数小于上一步相应值,但二者差值大于预定值,则进行下一步迭代,直到二者差值小于预定值,认为迭代收敛;
当某一迭代步计算出所有压型件节点残余力向量范数大于上一步相应值或迭代步数超过所设定的最大迭代步,认为迭代发散。
4.根据权利要求1所述的轨道车辆压型件毛坯快速预示方法,其特征是,所述网格划分的方法是:
首先选定需要划分的零件,即压型件,并赋予网格划分的参数,包括单元类型、单元大小、最小单元大小、弦高;
再进行几何清理与几何简化工作,包括针对一些有不重要的孔洞或缺陷的曲面进行补孔操作;
最后进行网格编辑,生成网格模型。
5.根据权利要求1所述的轨道车辆压型件毛坯快速预示方法,其特征是,所述等距偏置算法是:
首先计算每一个单元的法线与成形方向的单元夹角θ,当cosθ>0,保持其现有法线方向,当cosθ<0,将该单元的法线方向翻转180°;
再设定外表面或内表面的节点坐标为(x,y,z),偏置后的节点坐标为(x1,y1,z1),通过坐标变换,即原坐标分量沿法线方向移动一半厚度值,将外表面或内表面的节点偏置到中性层上。
6.根据权利要求1所述的轨道车辆压型件毛坯快速预示方法,其特征是,所述板料壁厚为板料的原始厚度。
7.根据权利要求1或2所述的轨道车辆压型件毛坯快速预示方法,其特征是,所述板料属性包括杨氏模量、泊松比、屈服强度、密度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910627696.3A CN110457755A (zh) | 2019-07-12 | 2019-07-12 | 一种轨道车辆压型件毛坯快速预示方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910627696.3A CN110457755A (zh) | 2019-07-12 | 2019-07-12 | 一种轨道车辆压型件毛坯快速预示方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110457755A true CN110457755A (zh) | 2019-11-15 |
Family
ID=68482708
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910627696.3A Pending CN110457755A (zh) | 2019-07-12 | 2019-07-12 | 一种轨道车辆压型件毛坯快速预示方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110457755A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111797459B (zh) * | 2020-07-21 | 2023-10-24 | 北京交通大学 | 一种有砟轨道-桥梁动力学耦合模型的构建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020177983A1 (en) * | 2001-04-17 | 2002-11-28 | Maker Bradley N. | Implicit-explicit switching for finite element analysis |
CN103970948A (zh) * | 2014-05-05 | 2014-08-06 | 中国人民解放军军械工程学院 | 一种金属橡胶优化设计与性能预报方法 |
CN104268349A (zh) * | 2014-09-30 | 2015-01-07 | 华南理工大学 | 一种复杂曲面下翻边修边线的精确控制方法 |
CN106874636A (zh) * | 2017-04-12 | 2017-06-20 | 南京航空航天大学 | 一种管材液压成形的快速预测方法 |
-
2019
- 2019-07-12 CN CN201910627696.3A patent/CN110457755A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020177983A1 (en) * | 2001-04-17 | 2002-11-28 | Maker Bradley N. | Implicit-explicit switching for finite element analysis |
CN103970948A (zh) * | 2014-05-05 | 2014-08-06 | 中国人民解放军军械工程学院 | 一种金属橡胶优化设计与性能预报方法 |
CN104268349A (zh) * | 2014-09-30 | 2015-01-07 | 华南理工大学 | 一种复杂曲面下翻边修边线的精确控制方法 |
CN106874636A (zh) * | 2017-04-12 | 2017-06-20 | 南京航空航天大学 | 一种管材液压成形的快速预测方法 |
Non-Patent Citations (4)
Title |
---|
吴建军等: "钣金零件毛坯展开计算方法研究进展", 《航空制造技术》 * |
张显藤: "橡皮囊成形钣金件初始毛坯快速设计系统的研究与实现", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
朱明华等: "飞机橡皮囊成形零件毛坯精确预示方法研究", 《航空学报》 * |
鲍益东等: "三维复杂零件修边线快速预示方法", 《机械工程学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111797459B (zh) * | 2020-07-21 | 2023-10-24 | 北京交通大学 | 一种有砟轨道-桥梁动力学耦合模型的构建方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100557614C (zh) | 一种生成三维有限元网格的方法 | |
CN102024081B (zh) | 一种汽车覆盖件模具设计方法 | |
US20170140081A1 (en) | Model setting method, forming simulation method, production method of forming tool, program, computer-readable recording medium having program recorded thereon, and finite element model | |
CN106994483B (zh) | 一种汽车覆盖件拉延模具精确型面加工的方法 | |
CN103617302A (zh) | 一种汽车覆盖件拉延模具精确型面构建的方法 | |
CN103433380B (zh) | 一种汽车覆盖件拉延模具的凸成形方法 | |
CN108704993B (zh) | 一种汽车覆盖件拉延模具研配压料面设计的方法 | |
CN102737147B (zh) | 确定板料多步成形的中间构型的几何参数的优化设计方法 | |
CN111177906B (zh) | 一种离散化模具型面精确补偿方法 | |
CN104765912A (zh) | 一种铝板冲压工艺的稳健性优化方法 | |
CN106599452A (zh) | 一种高强度钢管材三维弯曲成型回弹后形状的获取方法 | |
CN109127945A (zh) | 轻量化车身覆盖件冲压成形精度的调控方法 | |
Lan et al. | A method of constructing smooth tool surfaces for FE prediction of springback in sheet metal forming | |
CN104268349A (zh) | 一种复杂曲面下翻边修边线的精确控制方法 | |
JP2006263788A (ja) | 見込み金型形状設計システム | |
CN106021761A (zh) | 一种汽车覆盖件回弹评测方法 | |
CN110457755A (zh) | 一种轨道车辆压型件毛坯快速预示方法 | |
CN106980742B (zh) | 向有限元仿真模型引入冲压成型信息的载荷投影映射方法 | |
CN106326550B (zh) | 一种复杂蒙皮零件回弹补偿工艺模型建模方法 | |
CN109255141A (zh) | 一种汽车车身正向概念设计截面形状优化方法 | |
CN102722619A (zh) | 一种确定冲压形成汽车覆盖件的零件的材料利用率的方法 | |
Pilthammar et al. | Including die and press deformations in sheet metal forming simulations | |
CN102581183A (zh) | 确定二维整体加载成形用不等厚坯料的方法 | |
CN102632173A (zh) | 确定二维局部加载成形用不等厚坯料的方法 | |
CN110457754A (zh) | 一种轨道车辆压型件曲面翻边成形的预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191115 |
|
RJ01 | Rejection of invention patent application after publication |