CN108804783A - 一种毛坯初始残余应力的测评方法 - Google Patents

一种毛坯初始残余应力的测评方法 Download PDF

Info

Publication number
CN108804783A
CN108804783A CN201810516283.3A CN201810516283A CN108804783A CN 108804783 A CN108804783 A CN 108804783A CN 201810516283 A CN201810516283 A CN 201810516283A CN 108804783 A CN108804783 A CN 108804783A
Authority
CN
China
Prior art keywords
residual stress
pivot
matrix
deformation
stress
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.)
Granted
Application number
CN201810516283.3A
Other languages
English (en)
Other versions
CN108804783B (zh
Inventor
秦国华
林锋
王华敏
叶海潮
侯源君
韩雄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanchang Hangkong University
Original Assignee
Nanchang Hangkong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanchang Hangkong University filed Critical Nanchang Hangkong University
Priority to CN201810516283.3A priority Critical patent/CN108804783B/zh
Publication of CN108804783A publication Critical patent/CN108804783A/zh
Application granted granted Critical
Publication of CN108804783B publication Critical patent/CN108804783B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种毛坯初始残余应力的测评方法,根据零件加工变形要求,确定变形观测位置及其数量;根据精度要求设置零件划分层数,并逐层进行影响系数计算;确定调整后毛坯残余应力系数矩阵;基于旋转算法进行调整后的残余应力求解。本发明所涉及的计算模型及过程通用性好、逻辑性强,易于实现自动化,为残余应力引起加工变形的计算提供理论基础,能够基于工件变形要求对毛坯生产提出初始残余应力要求。

Description

一种毛坯初始残余应力的测评方法
技术领域
本发明涉及整体结构件制造与厚板生产领域,具体是指整体结构件制造与厚板生产领域基于单位力法与旋转算法的毛坯初始残余应力评测方法。
背景技术
航空整体结构件有着重量轻、装配时间短、生产效率高和生产成本低的优势,在航空领域被广泛使用。但由于其尺寸大且在加工过程中材料去除率高达90%以上造成刚度较低,毛坯中初始残余应力的存在容易引发加工变形。虽然通过预拉伸工艺能够大幅度地降低铝厚板中的残余应力水平,但是残留部分的残余应力仍能引发变形问题。为了研究航空整体结构件变形问题,可以通过有限元方法和解析式方法进行分析,其中有限元方法为主要的分析手段。
在有限元分析中,单元生死技术是最常用的材料去除仿真方法。Liu等人在使用层削法获得7085-T7452铝合金锻件初始残余应力的基础上,通过ANSYS有限元软件进行了挡风玻璃框架结构件变形分析,将仿真结果与实验结果比较得出有限元方法能够有效且可靠的预测初始残余应力引起的加工变形。Zhang等人在使用轮廓法获得7050-T7451铝合金厚板初始残余应力的基础上,通过有限元方法研究了T形件在不同位置处的加工变形,总结出适当的厚度方向位置可以有效地降低零件变形。通过单元生死技术进行材料去除,杀死的单元在仿真前已经设定,但是被装夹的工件在加工过程中仍存在局部变形,这将导致无法对工件的加工尺寸精度进行预测。为解决该问题,Cerutti等人以层剥法测量得到毛坯中的初始残余应力,通过在FORGE 有限元软件中加入布朗运算,同时实现了对加工整体变形和加工尺寸精度的预测,并对比实验证明了预测方法的正确性。
在解析式分析中,Nervi等人将初始残余应力拟合为空间坐标系中的函数,通过Navier-Lamè方程对7050-T7451铝合金厚板的加工变形进行了预测,并结合实验对比证明了模型的有效性。但是这种解析式分析方法同样需要结合有限元手段,难以在宏观上揭示残余应力引起加工变形的内在规律。目前的纯解析预测方法局限性较大,仅能分析材料去除方式为简单厚板层削和开单槽的弯曲变形情况。
发明内容
本发明的目的是为了进一步拓展解析式的使用范围,实现宏观上的残余应力引起变形分析,对航空整体结构件的变形进行解析式推导,而提出的结合有限元的单位力计算方法。
本发明采用以下技术方案实现上述目的。一种毛坯初始残余应力的测评方法,其步骤如下:
1)根据零件加工变形要求,确定变形观测位置;
观测位置以实际图纸要求出发,将基准位置和变形测量位置设为观测位置,或在分析中通过将基准位置设为不动点,对变形测量位置进行观测;
2)根据精度要求设置零件划分层数,并逐层进行影响系数计算;
使用有限元分析软件,将零件沿厚度方向分为n层,单元类型使用二阶单元;依次对第i 层的轧制方向和横向方向施加大小为1MPa的残余应力,计算得到影响因子DZi(x)和DZi(y);则该点的变形值其中ΣXi(z)和ΣYi(z)(1≤k≤n)分别为轧制方向和横向方向上的等效应力;
3)确定调整后毛坯残余应力系数矩阵;
根据实际的加工要求,对零件额外施加一个调整应力ΣOPTXi,记∑′Xi(z)=∑Xi(z)+∑OPTXi(z)为调整后的毛坯残余应力,控制零件变形在2ΔZ范围内,使得
由于毛坯初始残余应力具有自平衡特性,因此调整应力函数也应具有自平衡特性,即存在关系通过引入一个单位变量zn+1=1,结合夹逼原则可得AZ≤0;
式中:系数矩阵N=2n,待定的毛坯应力向量 系数dij、δZ分别为Dij=DZj(xi)、ΔZ的归一化值,应力zi=Σ′i=Σ′Xi(z),1≤i≤n;
4)基于旋转算法进行调整后的残余应力求解;
采用旋转算法求解变量Z,假定不等式解的一般形式为Z=Vν+Wω,其中ν为任意的系数向量,ω为任意的非负系数向量,V、W为基向量矩阵;记q、p分别为A=[a1,a2,...,ap]的行数和列数,则由A既可生成凸锥Aω={x|x=ω1a12a2+...+ωpap},也可生成线性空间Av={y|y=ν1a12a2+...+νpap}(ωi≥0,1≤i≤p);再记矩阵B={ai|ai∈Aπ,且-ai∈Aπ},那么求解基向量矩阵V和W的旋转算法如下:
S01:初始化基向量矩阵:
h=1;V=Iq×q,W=Φ,U=[V,W](Iq×q为单位矩阵,Φ为空集);Ij=Φ(1≤j≤q);
S02:搜索旋转枢轴;
从t=UTah中,找出tk≠0且k≤cV;其中rV、rW、rU和cV、cW、cU分别为V、W、U的行数和列数,然后执行步骤S03;否则执行步骤S04;
S03:执行旋转运算:
令枢轴号为pivot=k,将枢轴列标准化为除枢轴(即第pivot列)外,U中其它所有元素进行如下的枢轴变换,即uk=uk+tkupivot,k≠pivot;如果ah∈B,将upivot从V中剔除,即V=V-upivot如果将upivot从V中剔除,添加至W中,即将h 添加到集合Ij(j≠pivot,1≤j≤cV)中,即Ij=Ij Y{h}后,执行步骤S05;
S04:执行非旋转运算:
将h添加至所有满足tk=0的集合Ij中,即Ij=Ij Y{h}(1≤j≤cU)后,令矩阵Z={wk(i,j)=wi+wj|ti<0,tj>0,wi∈W,wj∈W},则Ik(i,j)=(IiΙIj)Y{h};从所有Ik(i,j)中找出满足的集合Ir(a,b);其中:wf是满足tf=0的列向量,则得Z的子集Z*={wi|i∈Ir(a,b)};
然后从W中剔除满足tj>0的列向量wj;若ah∈B,再从W中剔除所有满足的列向量w,并将Z*中所有向量加至W,即W={W,Z*};
S05:判断终止条件:
如果h<p,则h=h+1;计算Ij(1≤j≤rU)后转至步骤2;否则结束并返回V和W;
最后,按照旋转算法求出基向量矩阵V和W后,将zn+1=1(一维)或zN+1=1(二维)代入式Z=Vν+Wω后,经整理可得调整后的等效残余应力。
本发明所涉及的计算模型及过程通用性好、逻辑性强,易于实现自动化,为残余应力引起加工变形的计算提供理论基础,能够基于工件变形要求对毛坯生产提出初始残余应力要求。
具体实施方式
下面结合实施例对本发明进行详细描述:
本发明的工件的材料选用经过拉伸率为δ=2.4%的7075-T73铝合金板材,其弹性模量为 E=70.5Gpa,泊松比为ν=0.33。毛坯尺寸为长L=200mm,宽W=20mm,高H=6mm,而铣削成形的零件尺寸为长l'=200mm,宽w'=20mm,高h'=4.5mm,要求变形不超过0.1 mm。
步骤1:确定观测位置及数量:
由于零件具有对称性,选择长度方向55mm和100mm位置作为观测点进行分析。
步骤2:零件单元划分:
沿厚度方向Z对工件进行分层,层数n=20,每层厚度为Δz=0.3mm。若选取每层中心的残余应力作为该层的应力值,通过静力平衡条件求解出每层的应力值如表1所示。
表1毛坯内每层残余应力值
步骤3:影响系数计算:
有限元分析中单元类型采用CPS8R,网格密度设定为0.15,总共有53320个单元。选择“Static,General”分析步类型后,使用Model change进行“单元杀死”,模拟材料的去除。零件沿厚度方向Z从z=0至z=4.5mm共分为15层,考虑距离端面分别为x=55mm和x=100mm位置的A、B两点,通过逐层施加单位力,利用有限元方法计算的单位变形(即影响系数)如表2所示。
表2影响系数
步骤4:确定调整后残余应力系数矩阵:
由于薄板件的毛坯初始残余应力分布和结构均呈对称形式,55mm和100mm处两点的变形线性相关,因此需控制100mm处变形不超出ΔZ=0.05即可。结合表2的数据,根据式自平衡特性,可知应力调整模型中z16=1,且系数矩阵A如表3所示。
表3归一化的系数矩阵
步骤5:调整后残余应力求解:
初始化生成元矩阵h=1。V=I16×16,W=Φ,U=[V,W]。Ij=Φ(1≤j≤16)。
由于t=[-1,-0.8561,-0.7131,-0.5704,-0.4278,-0.2848,-0.1423,0.0003,0.1429, 0.2855,0.4281,0.5707,0.7133,0.8559,0.9985,-0.0008]T,可见t1≠0且u1在V中,故pivot=1。
执行第一次旋转变换得矩阵如表4所示。然后将h=1添加到Ij(j≠1)中,即Ij={1}(2≤j≤16)。此时将枢轴列u1从V中剔除并加至W中。
表4第一次旋转运算
由于h=1,p=4,h<p,故h=2。
又,t=[1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-0.0016]T,由于t16≠0且 u16列在V中,则枢轴povit=16。执行第二次旋转变换得矩阵如表5所示。此时将枢轴列 u16从V中剔除并加至W中。
表5第二次旋转运算
由于依然存在h<p,故h=3。
t=[0.9992,0.1440,0.2869,0.4296,0.5722,0.7152,0.8578,1.0003,1.1429,1.2855,1.4281,1.5707,1.7133,1.8559,1.9985,-0.5]T,由于t2≠0且u2在V中,则povit =2。执行第三次旋转变换得矩阵如表6所示。由于a3∈B,应将枢轴列u2从V中剔除。
表6第三次旋转运算
更新生成元矩阵,此时h=3,p=4,h<p的关系可进一步令h=4。
计算得t=[-1.1430,0,0,0,0,0,0,0,0,0,0,0,0,0,0]T,无枢轴。由于 t1<0,t16=0,故Z=Φ。因此Z*=Φ。由于a4∈B,从W中剔除u1后得W=u16。又,因此最终有W=u16,如表7所示。
表7第四次非旋转运算
至此,h=4,又p=4,故基向量矩阵最终为V=[u3,u4,...,u15]和W=u16
最后,由于z16=1,经整理后可得毛坯内的应力X应调整如下后才能保证变形在0.1mm的范围之内,即:
其中vi(1≤i≤13)为任意数,w≥0。

Claims (1)

1.一种毛坯初始残余应力的测评方法,其步骤如下:
1)根据零件加工变形要求,确定变形观测位置;
观测位置以实际图纸要求出发,将基准位置和变形测量位置设为观测位置,或在分析中通过将基准位置设为不动点,对变形测量位置进行观测;
2)根据精度要求设置零件划分层数,并逐层进行影响系数计算;
使用有限元分析软件,将零件沿厚度方向分为n层,单元类型使用二阶单元;依次对第i层的轧制方向和横向方向施加大小为1MPa的残余应力,计算得到影响因子DZi(x)和DZi(y);则该点的变形值其中ΣXi(z)和ΣYi(z)(1≤k≤n)分别为轧制方向和横向方向上的等效应力;
3)确定调整后毛坯残余应力系数矩阵;
根据实际的加工要求,对零件额外施加一个调整应力ΣOPTXi,记∑′Xi(z)=∑Xi(z)+∑OPTXi(z)为调整后的毛坯残余应力,控制零件变形在2ΔZ范围内,使得
由于毛坯初始残余应力具有自平衡特性,因此调整应力函数也应具有自平衡特性,即存在关系通过引入一个单位变量zn+1=1,结合夹逼原则可得AZ≤0;
式中:系数矩阵N=2n,待定的毛坯应力向量系数dij、δZ分别为Dij=DZj(xi)、ΔZ的归一化值,应力zi=Σ′i=Σ′Xi(z),1≤i≤n;
4)基于旋转算法进行调整后的残余应力求解;
采用旋转算法求解变量Z,假定不等式解的一般形式为Z=Vv+Wω,其中ν为任意的系数向量,ω为任意的非负系数向量,V、W为基向量矩阵;记q、p分别为A=[a1,a2,...,ap]的行数和列数,则由A既可生成凸锥Aω={x|x=ω1a12a2+...+ωpap},也可生成线性空间Aν={y|y=ν1a12a2+...+νpap}(ωi≥0,1≤i≤p);再记矩阵B={ai|ai∈Aπ,且-ai∈Aπ},那么求解基向量矩阵V和W的旋转算法如下:
S01:初始化基向量矩阵:
h=1;V=Iq×q,W=Φ,U=[V,W](Iq×q为单位矩阵,Φ为空集);Ij=Φ(1≤j≤q);
S02:搜索旋转枢轴;
从t=UTah中,找出tk≠0且k≤cV;其中rV、rW、rU和cV、cW、cU分别为V、W、U的行数和列数,然后执行步骤S03;否则执行步骤S04;
S03:执行旋转运算:
令枢轴号为pivot=k,将枢轴列标准化为除枢轴(即第pivot列)外,U中其它所有元素进行如下的枢轴变换,即uk=uk+tkupivot,k≠pivot;如果ah∈B,将upivot从V中剔除,即V=V-upivot如果将upivot从V中剔除,添加至W中,即将h添加到集合Ij(j≠pivot,1≤j≤cV)中,即Ij=IjY{h}后,执行步骤S05;
S04:执行非旋转运算:
将h添加至所有满足tk=0的集合Ij中,即Ij=IjY{h}(1≤j≤cU)后,令矩阵Z={wk(i,j)=wi+wj|ti<0,tj>0,wi∈W,wj∈W},则Ik(i,j)=(IiIIj)Y{h};从所有Ik(i,j)中找出满足的集合Ir(a,b);其中:wf是满足tf=0的列向量,则得Z的子集Z*={wi|i∈Ir(a,b)};
然后从W中剔除满足tj>0的列向量wj;若ah∈B,再从W中剔除所有满足的列向量w,并将Z*中所有向量加至W,即W={W,Z*};
S05:判断终止条件:
如果h<p,则h=h+1;计算Ij(1≤j≤rU)后转至步骤2;否则结束并返回V和W;
最后,按照旋转算法求出基向量矩阵V和W后,将zn+1=1(一维)或zN+1=1(二维)代入式Z=Vν+Wω后,经整理可得调整后的等效残余应力。
CN201810516283.3A 2018-05-25 2018-05-25 一种毛坯初始残余应力的测评方法 Active CN108804783B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810516283.3A CN108804783B (zh) 2018-05-25 2018-05-25 一种毛坯初始残余应力的测评方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810516283.3A CN108804783B (zh) 2018-05-25 2018-05-25 一种毛坯初始残余应力的测评方法

Publications (2)

Publication Number Publication Date
CN108804783A true CN108804783A (zh) 2018-11-13
CN108804783B CN108804783B (zh) 2022-06-14

Family

ID=64089062

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810516283.3A Active CN108804783B (zh) 2018-05-25 2018-05-25 一种毛坯初始残余应力的测评方法

Country Status (1)

Country Link
CN (1) CN108804783B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109492323A (zh) * 2018-11-25 2019-03-19 南昌航空大学 基于变形测试的毛坯初始残余应力高精度反求方法
CN111062095A (zh) * 2019-11-07 2020-04-24 南京航空航天大学 一种基于条件数的加工变形评价方法
CN112131683A (zh) * 2020-09-22 2020-12-25 河南科技大学 一种剖分式直齿锥齿轮的切齿加工轮坯参数确定方法
CN112378735A (zh) * 2020-10-28 2021-02-19 华南理工大学 一种考虑正交异性板残余应力效应的组合板梁单元分析法
CN113076674A (zh) * 2021-04-09 2021-07-06 华南理工大学 一种结构的建模及制造方法
CN113297758A (zh) * 2021-05-06 2021-08-24 南昌航空大学 一种大型复杂筋板件预锻成形初始毛坯优化设计方法
CN113297772A (zh) * 2021-06-07 2021-08-24 南昌航空大学 一种使用有限元软件模拟复合材料残余应力分布的方法
CN113505513A (zh) * 2021-07-21 2021-10-15 广西交科集团有限公司 一种锈蚀吊杆应力集中系数的计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006008554A1 (en) * 2004-07-23 2006-01-26 Isis Innovation Limited Method of determining residual stress within an object
CN104077442A (zh) * 2014-06-24 2014-10-01 南京晨光集团有限责任公司 基于有限元分析的大型整体薄壁件加工精度控制方法
CN105302970A (zh) * 2015-11-04 2016-02-03 沈阳黎明航空发动机(集团)有限责任公司 航空薄壁件残余应力释放-重新分布过程的预测方法
JP2017078699A (ja) * 2015-10-22 2017-04-27 株式会社神戸製鋼所 残留応力評価方法
CN107704657A (zh) * 2017-09-06 2018-02-16 北京航空航天大学 一种由残余应力释放引起的板类零件翘曲变形的计算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006008554A1 (en) * 2004-07-23 2006-01-26 Isis Innovation Limited Method of determining residual stress within an object
CN104077442A (zh) * 2014-06-24 2014-10-01 南京晨光集团有限责任公司 基于有限元分析的大型整体薄壁件加工精度控制方法
JP2017078699A (ja) * 2015-10-22 2017-04-27 株式会社神戸製鋼所 残留応力評価方法
CN105302970A (zh) * 2015-11-04 2016-02-03 沈阳黎明航空发动机(集团)有限责任公司 航空薄壁件残余应力释放-重新分布过程的预测方法
CN107704657A (zh) * 2017-09-06 2018-02-16 北京航空航天大学 一种由残余应力释放引起的板类零件翘曲变形的计算方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
万晓航等: "基于毛坯残余应力框类零件铣削加工变形研究", 《煤矿机械》 *
刘秦龙等: "开槽应力释放法优化整体结构件切削加工", 《机械设计与制造工程》 *
孙杰等: "残余应力对航空整体结构件加工变形的影响分析", 《机械工程学报》 *
王春等: "残余应力作用下的航空整体结构件加工变形仿真", 《机械研究与应用》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109492323B (zh) * 2018-11-25 2022-02-25 南昌航空大学 基于变形测试的毛坯初始残余应力高精度反求方法
CN109492323A (zh) * 2018-11-25 2019-03-19 南昌航空大学 基于变形测试的毛坯初始残余应力高精度反求方法
CN111062095A (zh) * 2019-11-07 2020-04-24 南京航空航天大学 一种基于条件数的加工变形评价方法
CN112131683A (zh) * 2020-09-22 2020-12-25 河南科技大学 一种剖分式直齿锥齿轮的切齿加工轮坯参数确定方法
CN112131683B (zh) * 2020-09-22 2024-04-19 河南科技大学 一种剖分式直齿锥齿轮的切齿加工轮坯参数确定方法
CN112378735A (zh) * 2020-10-28 2021-02-19 华南理工大学 一种考虑正交异性板残余应力效应的组合板梁单元分析法
CN112378735B (zh) * 2020-10-28 2021-11-09 华南理工大学 一种考虑正交异性板残余应力效应的组合板梁单元分析法
CN113076674A (zh) * 2021-04-09 2021-07-06 华南理工大学 一种结构的建模及制造方法
CN113297758A (zh) * 2021-05-06 2021-08-24 南昌航空大学 一种大型复杂筋板件预锻成形初始毛坯优化设计方法
CN113297772A (zh) * 2021-06-07 2021-08-24 南昌航空大学 一种使用有限元软件模拟复合材料残余应力分布的方法
CN113297772B (zh) * 2021-06-07 2022-04-01 南昌航空大学 一种使用有限元软件模拟复合材料残余应力分布的方法
CN113505513A (zh) * 2021-07-21 2021-10-15 广西交科集团有限公司 一种锈蚀吊杆应力集中系数的计算方法
CN113505513B (zh) * 2021-07-21 2022-06-07 广西交科集团有限公司 一种锈蚀吊杆应力集中系数的计算方法

Also Published As

Publication number Publication date
CN108804783B (zh) 2022-06-14

Similar Documents

Publication Publication Date Title
CN108804783A (zh) 一种毛坯初始残余应力的测评方法
Butuc et al. An experimental and theoretical analysis on the application of stress-based forming limit criterion
Habibi et al. Determination of forming limit diagram using two modified finite element models
Gorji et al. Predicting shear fracture of aluminum 6016-T4 during deep drawing: Combining Yld-2000 plasticity with Hosford–Coulomb fracture model
CN109870362B (zh) 一种高强铝合金板材的断裂成形极限图建立方法及系统
CN110147602B (zh) 一种建立折弯回弹角预测模型的方法及其应用
Cao et al. Study of fracture behavior for anisotropic 7050-T7451 high-strength aluminum alloy plate
JP2007232714A (ja) 破断限界取得方法及び装置、並びにプログラム及び記録媒体
Le Maoût et al. Classical and roll-hemming processes of pre-strained metallic sheets
Chen Experimental study on punch radius and grain size effects in V-bending process
Lawanwong et al. Double-action bending for eliminating springback in hat-shaped bending of advanced high-strength steel sheet
Haichao et al. A machining position optimization approach to workpiece deformation control for aeronautical monolithic components
Rahmaan et al. Characterization and prediction of fracture in 6000-and 7000-series aluminum alloy sheet under various stress states
Majidi et al. Effect of slide motion on springback in 2-D draw bending for AHSS
Qiang et al. Hot deformation behavior of 7A04 aluminum alloy at elevated temperature: Constitutive modeling and verification
Xue et al. Control strategy of twist springback for aluminium alloy hybrid thin-walled tube under mandrel-rotary draw bending
CN116542082A (zh) 一种膜片热压成型时变形预测方法、装置、设备和介质
Wisselink et al. Finite Element Simulation of the Stretch‐Forming of Aircraft Skins
Keshtiban et al. Die and process parameters effects on ECAP process of sheet-type samples
Chen et al. Forming limit evaluation for AA5182 aluminum alloy at warm temperatures based on M–K model
Shukla et al. Analysis of significant parameters influencing formability of titanium alloy by using over all evaluation criteria and new matrix model based on Taguchi method
Athale et al. Analytical and finite element simulation studies on earing profile of Ti-6Al-4V deep drawn cups at elevated temperatures
Xu et al. The Influence Factors of Incremental Forming Geometric Accuracy based on the Weighted Gray Correlation Analysis Method
Asgari et al. Investigation of punching parameters effect on mechanical properties of al-1100-o in incremental sheet metal hammering process
Shang et al. Anisotropic plasticity deformation during micro-deep drawing of 304 foils: An experimental and numerical investigation

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