CN113032908B - 薄壁构件空间包络成形翘曲变形预测方法 - Google Patents
薄壁构件空间包络成形翘曲变形预测方法 Download PDFInfo
- Publication number
- CN113032908B CN113032908B CN202110350910.2A CN202110350910A CN113032908B CN 113032908 B CN113032908 B CN 113032908B CN 202110350910 A CN202110350910 A CN 202110350910A CN 113032908 B CN113032908 B CN 113032908B
- Authority
- CN
- China
- Prior art keywords
- blank
- contact area
- equation
- stress
- thin
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- G06F2113/00—Details relating to the application field
- G06F2113/28—Fuselage, exterior or interior
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
本发明涉及一种薄壁构件空间包络成形翘曲变形预测方法,该方法首先确定空间包络成形过程中包络模和坯料之间时变接触区的几何尺寸,在特定坐标系下,计算坯料接触区部分任意位置的应力,在此基础上计算坯料非接触区部分任意位置的应力,构建坯料非接触区部分翘曲平衡微分方程,结合坯料非接触区部分边界条件,求解翘曲平衡微分方程,获得薄壁构件空间包络成形翘曲量。本发明可以精确预测高性能薄壁构件空间包络成形翘曲变形,进而提高薄壁构件空间包络成形精度,同时为薄壁构件空间包络成形工艺条件优化设计提供依据。
Description
技术领域
本发明涉及薄壁构件制造领域,更具体地说,涉及一种薄壁构件空间包络成形翘曲变形预测方法。
背景技术
薄壁构件具有质量轻、承载力大的特点,广泛应用于机翼、机身、燃料储箱箱体等航空航天领域。目前,薄壁构件的主要制造方法是机械加工方法,但是该方法存在性能差、加工效率低、材料利用率低、生产成本高等问题。
薄壁构件空间包络成形方法是一种薄壁构件成形新方法,具有成形载荷小,成形效率高,构件性能优良等优点。然而,薄壁构件厚度小、直径大,在空间包络成形包络模偏心载荷作用下极易发生翘曲变形,从而影响构件精度,甚至导致构件折叠破裂报废,精确高效预测翘曲变形是实现高性能薄壁构件空间包络成形的前提。目前,还没有关于薄壁构件空间包络成形翘曲变形预测方法的报道。
发明内容
本发明要解决的技术问题在于,提供一种薄壁构件空间包络成形翘曲变形预测方法,能够,为高性能薄壁构件空间包络成形精度控制和工艺优化提供依据。
本发明解决其技术问题所采用的技术方案是:构造一种薄壁构件空间包络成形翘曲变形预测方法,包括以下步骤:
S1、确定空间包络成形过程中坯料和包络模之间时变接触区的几何尺寸;其中,坯料初始厚度为h0,初始外径为R0,瞬时进给量ht的计算公式如方程(1)所示,坯料接触区平均厚度h1的计算公式如方程(2)所示,坯料非接触区平均厚度h2的计算公式如方程(3)所示,坯料接触区的瞬时半径Rt计算公式如方程(4)所示,坯料接触区的圆心角计算公式如方程(5)所示;
h1=h0-vt (2)
式中,γ为包络模轴线和z轴夹角;
S2、建立极坐标系,其中极坐标平面和下模的上端面平行,并且位于坯料非接触区部分厚度二分之一处,极坐标原点O是构件和包络模接触区的顶点,极轴Oρ将构件和包络模接触区圆心角平分,极坐标z轴经过原点O并且垂直平面;
式中,c是应变速率比,k是构件材料的剪切屈服强度,μ是坯料和模具之间的摩擦系数,σs是构件材料的屈服强度;
S4、计算坯料非接触区任意位置应力;以接触区的应力为边界条件,根据应力函数法获得非接触区任意位置应力的计算公式如方程(7)所示;
S5、构建翘曲平衡微分方程,确定边界条件;选取非接触区中任意单元体进行受力分析,构建翘曲平衡微分方程如方程(9)所示,边界条件如方程(10)所示;
式中,D是构件的弯曲刚度,w是构件的翘曲量。
S6、通过软件求解翘曲平衡微分方程,获得薄壁构件空间包络成形翘曲量。
上述方案中,步骤S3中所述的应变速率比c的计算公式如方程(11)所示:
式中,Rt2是成形时间(t-Δt)时坯料接触区瞬时半径。
上述方案中,步骤S5中所述的构件弯曲刚度D的计算公式如方程(12)所示:
式中,E是构件材料的弹性模量,μ是坯料和模具之间的摩擦系数。
上述方案中,步骤S5中所述的坯料非接触区径向力矩Mρ的计算公式如方程(13)所示:
上述方案中,步骤S5中所述的坯料非接触区边界ρa的计算公式如方程(15)所示:
实施本发明的薄壁构件空间包络成形翘曲变形预测方法,具有以下有益效果:
(1)本发明可以精确预测高性能薄壁构件空间包络成形翘曲变形,进而提高薄壁构件空间包络成形精度,同时为薄壁构件空间包络成形工艺条件优化设计提供依据。
(2)本发明可以高效预测高性能薄壁构件空间包络成形翘曲变形,进而提高薄壁构件空间包络成形工艺设计效率,缩短产品开发周期。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是薄壁构件空间包络成形初始状态示意图;
图2是薄壁构件空间包络成形翘曲变形示意图;
图3是薄壁构件空间包络成形初始状态坯料几何尺寸示意图;
图4是薄壁构件空间包络成形过程示坯料几何尺寸示意图;
图5是薄壁构件空间包络成形坯料接触区受力示意图;
图6是薄壁构件空间包络成形坯料非接触区受力示意图;
图7是构件在=120°位置的理论预测翘曲量和模拟结果对比图;
图8是构件在=180°位置的理论预测翘曲量和模拟结果对比图;
图9是构件在=180°位置的理论预测翘曲量和实验结果对比图;
图10是薄壁构件空间包络成形翘曲变形模式为凸型的示意图;
图11是薄壁构件空间包络成形翘曲变形模式为凹型的示意图;
图12是薄壁构件空间包络成形翘曲变形量随每转进给量变化图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
本发明一种薄壁构件空间包络成形翘曲变形预测方法包括以下步骤:
S1、确定薄壁构件空间包络成形原理。薄壁构件空间包络成形初始状态如图1所示,图1中1是包络模,2是坯料,3是下模。包络模同时绕自身轴线和薄壁构件中心轴线以转速n=4rad/s转动,下模以进给速度v=2mm/s带动放置于下模内的坯料向包络模靠近。随着下模进给量的增加,坯料和包络模接触区部分材料在包络模和下模的共同作用下发生塑性变形,由于包络模的转动,接触区位置不断改变,坯料每个位置都发生塑性变形,并向包络模和下模构成的型腔流动,最终坯料被成形成薄壁构件。在空间包络成形过程中,坯料和包络模非接触区部分材料在包络模施加的偏心载荷作用下容易发生翘曲变形(如图2所示),图2中4是坯料发生翘曲变形位置。
S2、确定薄壁构件空间包络成形翘曲量计算原理。确定空间包络成形过程中包络模和薄壁构件之间时变接触区的几何尺寸,在特定坐标系下,计算坯料接触区部分任意位置的应力,在此基础上计算坯料非接触区部分任意位置的应力,构建坯料非接触区部分翘曲平衡微分方程,结合坯料非接触区部分边界条件,求解翘曲平衡微分方程,获得薄壁构件空间包络成形翘曲量。
S3、确定空间包络成形过程中构件和包络模之间时变接触区的几何尺寸。坯料初始厚度h0=3mm,初始外径为R0=50mm(如图3所示)。坯料空间包络成形过程中的几何尺寸如图4所示,当成形时间t=1s时,根据方程(1)计算出瞬时进给量ht=0.5mm/rad,根据方程(2)计算出构件和包络模接触区平均厚度h1=2mm,根据方程(3)计算出坯料非接触区平均厚度h2=2.25mm,根据方程(4)计算出接触区的瞬时半径Rt=58.15mm,根据方程(5)接触区的圆心角
h1=h0-vt (2)
式中,包络模轴线和z轴夹角γ=2°。
S4、建立极坐标系,其中极坐标平面和下模的上端面平行,并且位于坯料非接触区部分厚度二分之一处,极坐标原点O是构件和包络模接触区的顶点,极轴Oρ将构件和包络模接触区圆心角平分,极坐标z轴经过原点O并且垂直平面,如图5所示。
S5、计算坯料接触区任意位置应力。根据主应力法,选取接触区中任意单元体进行受力分析,如图5所示。根据方程(6)计算出坯料接触区中z=0mm,ρ=25mm,位置的应力,其中,周向应力径向应力σρ1=144.53MPa,轴向应力σz1=180.51MPa。
式中,应变速率比c根据方程(11)计算,构件材料的剪切屈服强度k=20.2MPa,坯料和模具之间的摩擦系数μ=0.14,构件材料的屈服强度σs=35MPa。
S6、计算坯料非接触区应力。以坯料接触区的应力为边界条件,根据应力函数法获得坯料非接触区中z=0mm,ρ=25mm,位置的应力,坯料非接触区径向应力σρ=158.7MPa,坯料非接触区周向应力坯料非接触区切应力
式中,系数a1、a2、a3和a4根据方程(8)计算。
S7、构建翘曲平衡微分方程,确定边界条件。坯料非接触区的受力分析如图6所示。
构建翘曲平衡微分方程如方程(9)所示,边界条件如方程(10)所示。
式中,构件的弯曲刚度D根据方程()计算,w是构件的翘曲量。
上述步骤中,根据方程(11)计算出应变速率比c=0.83。
式中,成形时间(t-Δt)时坯料接触区瞬时半径Rt2=52.55mm。
上述步骤中,根据方程(12)计算出构件弯曲刚度D=67.8kN·mm。
式中,构件材料的弹性模量E=70GPa,坯料和模具之间的摩擦系数μ=0.14。
上述步骤中,坯料非接触区径向力矩Mρ根据方程(13)计算。
上述步骤中,根据方程(15)计算出坯料非接触区边界ρa=54.24mm。
通过薄壁构件空间包络成形翘曲变形预测方法,对构件不同位置的薄壁构件空间包络成形翘曲变形进行了预测分析,并进行了有限元模拟和实验。构件在位置的理论预测翘曲量和模拟结果如图7所示,构件在位置的理论预测翘曲量和模拟结果如图8所示,构件在位置的理论预测翘曲量和实验结果如图9所示。通过对不同工艺参数下薄壁构件空间包络成形翘曲变形进行预测分析,薄壁构件空间包络成形翘曲变形有两种变形模式,分别是凸型模式(如图10所示)和凹型模式(如图11所示),薄壁构件空间包络成形翘曲变形最大翘曲量随每转进给量的增加而增加(如图12所示)。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。
Claims (6)
1.一种薄壁构件空间包络成形翘曲变形预测方法,其特征在于,包括以下步骤:
S1、确定空间包络成形过程中坯料和包络模之间时变接触区的几何尺寸;其中,坯料初始厚度为h0,初始外径为R0,瞬时进给量ht的计算公式如方程(1)所示,坯料接触区平均厚度h1的计算公式如方程(2)所示,坯料非接触区平均厚度h2的计算公式如方程(3)所示,坯料接触区的瞬时半径Rt计算公式如方程(4)所示,坯料接触区的圆心角计算公式如方程(5)所示;
h1=h0-vt (2)
式中,γ为包络模轴线和z轴夹角;
S2、建立极坐标系,其中极坐标平面和下模的上端面平行,并且位于坯料非接触区部分厚度二分之一处,极坐标原点O是构件和包络模接触区的顶点,极轴Oρ将构件和包络模接触区圆心角平分,极坐标z轴经过原点O并且垂直平面;
式中,c是应变速率比,k是构件材料的剪切屈服强度,μ是坯料和模具之间的摩擦系数,σs是构件材料的屈服强度;
S4、计算坯料非接触区任意位置应力;以接触区的应力为边界条件,根据应力函数法获得非接触区任意位置应力的计算公式如方程(7)所示;
S5、构建翘曲平衡微分方程,确定边界条件;选取非接触区中任意单元体进行受力分析,构建翘曲平衡微分方程如方程(9)所示,边界条件如方程(10)所示;
式中,D是构件的弯曲刚度,w是构件的翘曲量。
S6、通过软件求解翘曲平衡微分方程,获得薄壁构件空间包络成形翘曲量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110350910.2A CN113032908B (zh) | 2021-03-31 | 2021-03-31 | 薄壁构件空间包络成形翘曲变形预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110350910.2A CN113032908B (zh) | 2021-03-31 | 2021-03-31 | 薄壁构件空间包络成形翘曲变形预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113032908A CN113032908A (zh) | 2021-06-25 |
CN113032908B true CN113032908B (zh) | 2022-05-06 |
Family
ID=76453184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110350910.2A Active CN113032908B (zh) | 2021-03-31 | 2021-03-31 | 薄壁构件空间包络成形翘曲变形预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113032908B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114226529A (zh) * | 2021-12-20 | 2022-03-25 | 武汉理工大学 | 基于空间包络形变-热处理协同调控的薄壁构件制造方法 |
CN114239177B (zh) * | 2021-12-20 | 2023-05-16 | 武汉理工大学 | 薄壁高筋构件空间包络成形塑性失稳预测与控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019154000A1 (zh) * | 2018-02-09 | 2019-08-15 | 中铝材料应用研究院有限公司 | 基于残余应力释放翘曲量的测量预测金属板材残余应力的方法及装置 |
CN110918847A (zh) * | 2019-11-26 | 2020-03-27 | 武汉理工大学 | 一种薄壁高筋壁板构件空间包络成形方法 |
CN111090910A (zh) * | 2019-11-26 | 2020-05-01 | 武汉理工大学 | 提高锥形薄壁高筋构件空间包络成形精度的包络模补偿方法 |
CN111375716A (zh) * | 2019-11-26 | 2020-07-07 | 武汉理工大学 | 一种薄壁高筋飞机窗框构件空间包络成形方法 |
CN112434383A (zh) * | 2019-08-26 | 2021-03-02 | 波音公司 | 复合结构翘曲控制系统 |
-
2021
- 2021-03-31 CN CN202110350910.2A patent/CN113032908B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019154000A1 (zh) * | 2018-02-09 | 2019-08-15 | 中铝材料应用研究院有限公司 | 基于残余应力释放翘曲量的测量预测金属板材残余应力的方法及装置 |
CN112434383A (zh) * | 2019-08-26 | 2021-03-02 | 波音公司 | 复合结构翘曲控制系统 |
CN110918847A (zh) * | 2019-11-26 | 2020-03-27 | 武汉理工大学 | 一种薄壁高筋壁板构件空间包络成形方法 |
CN111090910A (zh) * | 2019-11-26 | 2020-05-01 | 武汉理工大学 | 提高锥形薄壁高筋构件空间包络成形精度的包络模补偿方法 |
CN111375716A (zh) * | 2019-11-26 | 2020-07-07 | 武汉理工大学 | 一种薄壁高筋飞机窗框构件空间包络成形方法 |
Non-Patent Citations (1)
Title |
---|
《A new cylindrical ring rolling technology for manufacturing thin-walled cylindrical ring》;韩星会;《International Journal of Mechanical Sciences》;20140831;95-108 * |
Also Published As
Publication number | Publication date |
---|---|
CN113032908A (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113032908B (zh) | 薄壁构件空间包络成形翘曲变形预测方法 | |
US10814370B2 (en) | Special-purpose die for shaping aluminum-magnesium alloy by rotating extrusion | |
CN104077439B (zh) | 一种新型高强钢轮辐拉深冲孔复合工序的数值模拟方法 | |
Wong et al. | Effects of roller path and geometry on the flow forming of solid cylindrical components | |
CN107262544A (zh) | 一种轻合金杯形件旋转挤压成形模具 | |
CN103699727A (zh) | 一种强力旋压的可旋性分析数值方法 | |
CN106552852A (zh) | 一种小锥度筒形件反向旋压加工方法 | |
Tian et al. | Effects of backward path parameters on formability in conventional spinning of aluminum hemispherical parts | |
Podder et al. | Flow forming of thin-walled precision shells | |
CN110508653A (zh) | 一种筋肋板的辊辗成形方法 | |
CN103934299B (zh) | 一种消除内筒应力集中的扁挤压筒及制造方法 | |
CN101840447B (zh) | 预测旋锻过程锻压力的有限元建模方法 | |
CN101462227A (zh) | 热锻复合模具制造方法和热锻复合模具 | |
Zhu et al. | Die structure optimization for eliminating premature folding of sidewall upsetting with a controllable deformation zone | |
Gong et al. | Research on continuous multi-point forming method for rotary surface | |
CN105537479A (zh) | 一种齿条冷摆辗精密成形方法 | |
CN114239177B (zh) | 薄壁高筋构件空间包络成形塑性失稳预测与控制方法 | |
CN106734207A (zh) | 一种滚轧锻式锻件制造模具及方法 | |
CN201711480U (zh) | 一种用于长轴类锻件制坯的胎模 | |
CN101773974A (zh) | 一种提高长轴类锻件制坯锻透性的胎模 | |
Zhu et al. | The Spinning Speed Influence on the Counter-Roller Spinning Process | |
CN111046507B (zh) | 锥形薄壁高筋构件空间包络成形干涉极限判定方法 | |
Li et al. | Finite element simulation of ultrasonic incremental forming | |
Yong et al. | Effect of internal pressure on corner radius and thickness distribution of shear hydro-bending of 5A02 aluminum alloy tube | |
Li et al. | Study on the Influence of Clearance Ratio between Mandrel and Workpiece on the Spinning Stability of Cylindrical Workpiece with Extra Diameter‐Thick Ratio |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |