×

电液成形过程中预测Cowper-Symonds本构方程材料参数的逆参数估计。 (英语) Zbl 07469792号

摘要:开发可靠的数值模拟对于理解电液成形(EHF)等高速成形过程至关重要。该数值模型应基于准确的材料特性创建。然而,应变率超过(1000s^{-1})时的动态材料特性无法通过实验方法轻易获得。因此,本研究基于逆参数估计对Cowper-Symonds本构方程中的两个材料参数进行了预测,从而使数值模拟预测的参数与从实验结果中获得的参数吻合良好。目标材料是1mm厚的Al 5052-H32片材。比较目标包括在6、7和8 kV三个输入电压下进行无EHF胀形试验时板材的最终变形形状。对于逆参数估计,两个参数的后验分布包括似然分布和先验分布。对于似然构造,提前开发了一个降阶代理模型来代替基于普通克里格法和主成分分析的数值模拟。此外,还得到了实验与降阶代理模型之间凸起高度的误差分布。将7kV处的先验分布定义为均匀分布,将7kV下的后验分布用作6kV、7kV和8kV处的前验分布。此外,采用马尔可夫链蒙特卡罗抽样,并采用Metropolis-Hastings算法获得符合后验分布的样本。在对样本独立性进行自相关计算后,选择具有(pm,0.02)区间自相关的滞后,并获得每个(text{lag}^{text{th}})样本。采集的样本总数为10^5,平均值根据采集的样本计算得出。结果表明,平均值的数值模拟与实验结果吻合良好。

MSC公司:

74-XX岁 可变形固体力学
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 布鲁诺,EJ,金属高速成形(1968),迪尔伯恩:美国工具与制造工程师协会,迪尔波恩
[2] Zittel G(2010)《高速金属成形的历史回顾》。参加:第四届高速成形国际会议
[3] Kim,J。;Song,W。;Kang,B.,利用Gurson塑性材料模型增强电磁成形的可成形性研究,Trans Korean Soc Autom Eng,2198-104(2013)·doi:10.7467/KSAE.2013.21.3.098
[4] Imbert J、Worswick M、Winkler S等人(2005)《电磁成形铝合金板材可成形性增加的分析》(编号:2005-01-0082)。SAE技术文件
[5] JM Imbert;SL温克勒;Worswick,MJ,《工具-板材相互作用对铝合金板材电磁成形中损伤演化的影响》,《工程材料技术杂志》,127145-153(2005)·数字对象标识代码:10.1115/1.1839212
[6] Balanethiram,VS;胡,X。;Altynova,M.,《超塑性:高速率下的增强成形性》,《材料加工技术杂志》,45,595-600(1994)·doi:10.1016/0924-0136(94)90404-9
[7] Altan,T。;Tekkaya,AE,《钣金成形:工艺和应用》(2012),材料园:ASM国际,材料园·doi:10.31399/asm.tb.smff.9781627083164
[8] Golowin,S。;卡马尔,M。;Shang,J.,均匀压力致动器在金属板电磁加工中的应用,Mater Eng Perform,16455-460(2007)·doi:10.1007/s11665-007-9085-4
[9] 卡马尔,M。;Daehn,GS,《成形平板的均匀压力电磁致动器》,《制造科学与工程杂志》,129369-379(2007)·数字对象标识代码:10.1115/12515481
[10] Neubauer,A。;施密克,H。;Schmidt,E.,《爆炸金属板成形温度影响的研究》,CIRP Ann,39,257-261(1990)·doi:10.1016/S0007-8506(07)61048-8
[11] Gerdooei M,Dariani B Mollaei,Liaghat GH(2009),爆炸成形中材料模型对板材成形性的影响。摘自:《世界工程大会论文集》,第2卷,第1741页
[12] Balanethiram,VS;Daehn,GS,《超塑性:高工件速度下增加的成形极限》,Scr Metal Mater,30,515-520(1994)·doi:10.1016/0956-716X(94)90613-0
[13] 斯默德,R。;温克勒,S。;Salisbury,C.,《汽车铝合金板的高应变率拉伸试验》,国际冲击工程杂志,32,541-560(2005)·doi:10.1016/j.ijimpeng.2005.04.013
[14] Yu,H。;Sun,L。;张,X。;王,S。;Li,C.,铝管电液成形和电磁成形实验,国际先进制造技术,89,9-12,3169-3176(2017)·doi:10.1007/s00170-016-9261-5
[15] Golovashchenko,旧金山;艾伦,JG;Alexander,VM,《电液成形中双相钢的成形性》,《材料加工技术杂志》,213,7,1191-1212(2013)·doi:10.1016/j.jmatprotec.2013.01.026
[16] Maris,C。;Hassannejadasl,A。;绿色,DE;Cheng,J。;Golovashchenko,旧金山;AJ吉拉德;Liang,Y.,DP600和AA5182板材的准静态和电液自由成形极限的比较,《材料加工技术杂志》,235206-219(2016)·doi:10.1016/j.jmatprotec.2016.04.028
[17] 艾哈迈德,M。;Kumar博士;Nabi,M.,《通过电液成形工艺提高AA5052合金板材的成形性》,《材料工程性能杂志》,26,1,439-452(2017)·doi:10.1007/s11665-016-2446-0
[18] 马穆托夫,AV;Golovashchenko,旧金山;马穆托夫,VS;Bonnen,JJ,钣金零件电液成形建模,材料加工技术杂志,21984-100(2015)·doi:10.1016/j.jmatprotec.2014.11.045
[19] Daehn GS、Fenton G、Vohnout VJ(2010),小位移电液成形过程中的压力不均匀性。参加:第四届高速成形国际会议,2010年3月9日至10日,美国俄亥俄州哥伦布
[20] Mane,T。;戈尔,V。;Kore,SD,板材电液成形的有限元建模,Procedia Mater Sci,6,105-114(2014)·doi:10.1016/j.mspro.2014.07.012
[21] ASTM国际(2011)ASTM E8/E8M-09金属材料拉伸试验的标准试验方法。ASTM,西康舍霍肯
[22] Cowper GR,Symonds PS(1957),悬臂梁冲击载荷中的应变硬化和应变率效应(编号TR-C11-28)。普罗维登斯布朗大学
[23] FJ泽里利;Armstrong,RW,材料动力学计算中基于位错力学的本构关系,应用物理学杂志,61,5,1816-1825(1987)·doi:10.1063/1.338024
[24] DJ斯坦伯格;科克伦,SG;Guinan,MW,适用于高应变率金属的本构模型,《应用物理学杂志》,51,3,1498-1504(1980)·数字对象标识代码:10.1063/1.327799
[25] Johnson GR,Cook WH(1983)金属在大应变、高应变率和高温下的本构模型和数据。载:第七届国际弹道学研讨会论文集,第21卷(1),第541-547页
[26] 北卡罗来纳州佩西尼奥。;Doellinger,C.,《耐撞应用轻合金动态材料特性的表征》,Mater Res,13,4,471-474(2010)·doi:10.1590/S1516-14392010000400008
[27] Jones,N.,《结构影响》(1989),剑桥:剑桥大学出版社,剑桥
[28] 曹,MC;坎贝尔,JD,《高应变率下金属和合金的塑性剪切特性》(1973),牛津:牛津大学工程科学系,牛津
[29] Calvetti D,Somersalo E(2007)《贝叶斯科学计算导论:关于主观计算的十堂课》。纽约州施普林格,第2页·Zbl 1137.65010号
[30] 吴,MA;Noh,HG;Song,WJ,使用al 5052-H34板材的电液成形数值模型的实验验证,国际先进制造技术杂志,93,1819-1828(2017)·doi:10.1007/s00170-017-0612-7
[31] 医学博士麦凯;RJ贝克曼;Conover,WJ,在计算机代码输出分析中选择输入变量值的三种方法的比较,Technometrics,42,55-61(2000)·网址:10.1080/00401706.2000.10485979
[32] Donea J、Huerta A、Ponthot JP、Rodriguez-Ferran A(2004)《任意拉格朗日-俄勒冈方法》。收录于:《计算力学百科全书》第1卷。奇切斯特·威利,第1-25页
[33] Belytschko,T。;弗拉纳根,DP;Kennedy,JM,流体-结构相互作用的用户控制网格有限元方法,计算方法应用机械工程,33,1-3,669-688(1982)·Zbl 0492.73089号 ·doi:10.1016/0045-7825(82)90127-X
[34] Gakwaya,A。;谢里菲,H。;吉洛,M。;Souli,M。;Erchiqui,F.,集成计算机辅助设计和工程系统中的ALE公式和模拟技术与工业金属成形应用,计算模型工程科学(CMES),73,3,209-266(2011)·Zbl 1223.74038号
[35] Hallquist,JO,LS-DYNA关键字用户手册970,299-800(2007),利弗莫尔:利弗莫雷软件技术公司,利弗莫
[36] Golovaschenko SF(2013)近净形状汽车面板的电液成形。Doe-FORD-GO18128号文件。福特汽车公司
[37] 马萨诸塞州奥利弗;韦伯斯特,R.,《克里格:地理信息系统的插值方法》,国际地理杂志。信息系统,4131-332(1990)·doi:10.1080/02693799008941549
[38] Lichtenstern A(2013)空间统计中的克里金方法。慕尼黑慕尼黑科技大学博士论文
[39] 马尔维?,T。;Bali,D.,《普通克里金方程中的线性和拉格朗日线性乘数》,Nafta,59,1,31-37(2009)
[40] 沃尔德,S。;Esbensen,K。;Geladi,P.,主成分分析,化学智能实验室系统,2,1-3,37-52(1987)·doi:10.1016/0169-7439(87)80084-9
[41] 阿卜迪,H。;Williams,LJ,主成分分析,Wiley Interdiscip Rev Comput Stat,2,4,433-459(2010)·doi:10.1002/wics.101
[42] Shlens J(2014)《主成分分析教程》。arXiv预打印arXiv:1404.1100
[43] Sirovich,L.,湍流和相干结构动力学。第一部分:连贯结构,Q Appl Math,45,3,561-571(1987)·Zbl 0676.76047号 ·doi:10.1090/qam/910462
[44] Newman AJ(1996)通过Karhunen-Loeve展开的模型简化,第1部分:阐述。技术报告T.R.96-32。马里兰大学系统研究所
[45] 安顿,JS;King,PI;Beran,PS,基于POD的变形网格降阶模型,数学计算模型,38,1-2,41-62(2003)·Zbl 1106.76415号 ·doi:10.1016/S0895-7177(03)90005-7
[46] 弗雷斯特,A。;Sobester,A。;Keane,A.,《通过代理建模进行工程设计:实用指南》(2008年),奇切斯特:威利·doi:10.1002/9780470770801
[47] 华盛顿州Link;Eaton,MJ,《关于MCMC中链的细化》,《方法生态进化》,3112-115(2012)·doi:10.1111/j.2041-210X.2011.00131.x
[48] 北卡罗来纳州大都会。;罗森布鲁斯,AW;罗森布鲁斯,明尼苏达州,用快速计算机器计算状态方程,化学物理杂志,211087-1092(1953)·Zbl 1431.65006号 ·数字对象标识代码:10.1063/1.1699114
[49] Geman,S。;Geman,D.,《随机松弛、吉布斯分布和图像的贝叶斯恢复》,IEEE Trans-Pattern Ana Mach Intell,6721-741(1984)·Zbl 0573.62030号 ·doi:10.1109/TPAMI.1984.4767596
[50] Gelman,A。;卡林,JB;斯特恩,HS;邓森,DB;Vehtari,A。;Rubin,DB,贝叶斯数据分析(2013),博卡拉顿:CRC出版社,博卡拉顿·doi:10.1201/b16018
[51] 罗伯茨,GO;Gelman,A。;Gilks,WR,随机行走Metropolis算法的弱收敛性和最优缩放,Ann Appl-Probab,7,110-120(1997)·Zbl 0876.60015号
[52] 罗伯茨,GO;Rosenthal,JS,各种Metropolis-Hastings算法的最佳缩放,《统计科学》,16,4,351-367(2001)·Zbl 1127.65305号 ·doi:10.1214/ss/1015346320
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。