基于弹塑性压力中性层偏移的棒材二辊矫直回弹模型-马自勇.pdf

《基于弹塑性压力中性层偏移的棒材二辊矫直回弹模型-马自勇.pdf》由会员分享，可在线阅读，更多相关《基于弹塑性压力中性层偏移的棒材二辊矫直回弹模型-马自勇.pdf（11页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第47卷第9期2叭6年9月中南大学学报(自然科学版)Journal of Central South UniVersity(Science and 1echnology)Vbl47 No9S印2016DoI：1011817巧issn1672-7207201609016基于弹塑性压力中性层偏移的棒材二辊矫直回弹模型马自勇1，马立峰1一，黄庆学1，马立东1，楚志兵1，黄志权(1太原科技大学山西省冶金设备设计理论与技术重点实验室，山西太原，030024；2吉林大学超塑性与塑性研究所，吉林长春，130000)摘要：基于二辊矫直变形工艺特点和弹塑性压力中性层偏移理论，确定应力中性层的位置；在考虑中性层偏

2、移和材料硬化的条件下，结合二辊矫直理论和矫直变形应力一应变关系，提出新的双线性矫直拟合方法：在此基础上利用微元法对棒料内部的金属进行弹塑性变形分析，给出新弯矩比公式的同时，提出中性层偏移半径比的概念；最后，结合纯弯回弹理论建立棒料矫直挠度回弹数学模型，并就中性层偏移对棒料矫直回弹影响进行实验分析和验证。实验结果表明：回弹理论计算值与实验数据吻合较好，其最大相对误差为645，明显低于忽略中性层偏移时的相对误差3165，证明了理论分析的正确性和模型的有效性。关键词：棒材矫直；中性层偏移半径比；回弹模型；实验分析中图分类号：TG3863+l 文献标志码：A 文章编号：16727207(2016109

3、302011Neutral layer offset rebound model baseiI On elastic-plastic pressureMA Ziyon91，MA Lifen91一，HUANG Qingxuel，MA Lidon91，CHU zhibin91，HUANG Zhiquanl(1Metallu唱ical Equipment Design Theory and 11echnology Key Laboratory of Shanxi ProVince，Taiyuan UniVerSity of Science and T-echn0109y，Taiyuan 030024

4、，China；2Superplastic时&Plasticity Institute，Jilin UniVersi吼Changch咖1 30000，China)Abstract：Ihe locatlon ot sness neutral layer was detemmed consldenng tlle pmcess chamctenstlcs ot似o。rolls仃aightening and the neu仃a11ayer o仃set meory of elasticplastic pressureConsidering the neutrallayer o低et andmaterial

5、 hardening conditions，the new bilinear 6tting method which was combined with tworoll straightening meo叮and s仃ess-strain was proposed0n the basis of fitting method and neu仃allayer ofrset，the elasticplastic defomation ofelement was analyzed，the formula of new moment ratio was given，and a new concept o

6、f neutrallayer ofIset mdius ratiowas also put fbnvardFinallythe calculation model of bar strai曲tening springback denection combined with the purebending sprin曲ack meory was established，and then the influence of neutrallayer o凰et rebound to accuracy of t、oroUs仃aightening was analyzed and Verified by

7、experimentsThe results show that meorctical data are in good agreement withexperimental data；the maximum relative ermr is 645，which is lower man 3 165when the neutral layer ofrset isignoredThe validity of theoretical analysis and model is proVedKey words：bar stmightening；neutrallayer oflset radius r

8、atjo；sprin曲ack model；experimental analysis收稿日期：201505一14；修回日期：2015_0919基金项目(Foundationitem)：国家重点基础研究发展规划(973计划)项口(2叭2CB722800)：中国博士后科学基金资助项目(2012M520677)：国家自然科学基金资助项目(51105264，51404160)：山两省科技攻关项I二I(20130321叭O一03)(Project(2012cB722800)supported by the National BaSicRcsearch Development Pm孕am(973 Prog

9、ram)of chilla：Project(20l 2M520677)supponedby the National Science Foundation for Post-doct啪1scientists of China；Projects(51 105264，51404160)supponed by Ihe NaIjonaI Natural science Foundan两of china；Project(20130321010-03)supponedby the Minis田ofMajor Science&Techn0109y of Shanxi Province)通信作者：马立峰，教授

10、，从事轧钢工艺与设备研究：Email：malifcngfqh163com万方数据第9期 马自勇，等：基于弹塑性压力中性层偏移的棒材二辊矫直回弹模型 30越随着石油、船舶、汽车等制造行业的飞速发展，对高强度合金钢棒材的需求量和质量要求日益提高，其中直线度要求达到l删：Tlm，甚至小于05 mrIlm，否则后续加工中材料的去除量非常大，损耗严重，成本急剧增加，因此，对棒料矫后直线度的要求十分严格。二辊矫直作为棒料精整的最后一道工序，是保证棒料矫后直线度指标的关键技术。二辊矫直是一个复杂的弹塑性变形过程，存在部分弹性变形，在卸载后弹性变形会立即恢复原态，出现回弹现象。回弹的结果使得棒料的残余曲率发生

11、严重变化，降低棒料的直线度，得到的棒料矫后直线度大于l删Tlm。由矫直曲率比公式【2】可知，无论是板料还是管棒材的矫直均是利用回弹来实现的。回弹量决定了矫直辊辊形的设计和工艺参数的制定，回弹量的预测精度也通过反弯量来影响棒料的残余挠度，最终得到高质量的矫直棒料，因此，准确有效地预测回弹是获得稳定的高精度棒料的基础。对于回弹问题，诸多学者进行了大量的研列3-14，但在预测棒料回弹量时都未考虑中性层偏移对回弹的影响，得到的回弹计算模型精度不高。官英平等【1 5】研究表明，在一定相对圆角半径下，中性层内移对板料弯曲回弹的相对误差可达70以上。相对于纯弯曲中性层偏移计算，运用本文的弹塑性压力中性层偏移

12、理论计算精度至少提高了346倍。因此，本文作者基于棒料二辊矫直变形特点和弹塑性压力中性层偏移理论，结合材料硬化规律和棒料弹塑性变形的力学关系，推导出棒料矫直过程中的挠度回弹计算模型，提出中性层偏移半径比的概念，分析了中性层偏移对同弹计算精度的影响，为后续深入研究高强度合金钢棒材矫直机理和矫直辊辊型优化设计与工艺模型的建立提供数捌基础。1棒料矫直挠度回弹模型的建立11基本假设1)变形前后棒料的横截面仍保持为平面，垂直于棒材轴线，且截面不发生畸变，截面上应变为线性分布。2)棒料是连续均质的弹塑性体，且弹性变形为线性，符合Hooke定律并认为塑性流动符合经典弹塑性理论的卸载规律。31棒料在矫商过程中

13、忽略截面的变形，并认为应变中性层与应力中性层重合。4)满足塑性变形体积不变原则，并且拉压变形部分金属的硬化规律相同，即应力一应变关系一致。矫直过程中，棒材旋转前进受凸辊压弯，并在压弯处产生塑性压缩(下阴影部分)与拉伸(上阴影部分)变形，如图1所示。在弹塑性变形区外层，切向应变为拉应变，切向应力为拉应力；变形区内层，切向应变为压应变，切向应力为压应力16。由假设3)知径向应变占，=0；同时各纤维间压应力可以忽略不计，因此，棒材矫直过程应力一应变状态如图2所示。棒十图l棒材矫直弹塑性变形区域Fig1 Elas“cplasric defonnation region of bar straighte

14、ning图2棒材矫直变形区应力一应变示意图Fig2 Diagram of bar stmightening stress and stmin indefbmation zone12基于弹塑-陛压力的中性层偏移分析121 中性层偏移几何应变分析进入矫直工序前，棒料在长度范围内基本都存在径向任意方位的弯曲，其弯曲程度各不同，但它们都会在矫直辊中央段统一为相同方位的弯曲，呈现扁担弯形态，为了便于理论分析，一般将初始弯曲简化为扁担形式的弯曲。现假设微元体爿BCD起初为拉伸状态，则矫直反弯后该位置处于压缩状态，如图2所示。由于被矫直的棒料原始弯曲程度不大，可认为棒料中性层与几何巾心轴是重合的，则棒料中性

15、层处金属纤维原始长度f0为=R岛 (1)式中：凰为初始弯曲半径；晶为微元体初始弯曲角；万方数据3022 中南大学学报(自然科学版) 第47卷微元体彳BCD处金属纤维初始长度，0D为f0D=场岛 (2)式中：托为微元体初始弯曲半径。棒材微元体被反弯后中性层长度fw为k=p氏 (3)其中：户为反弯后中性层半径；氏为矫直反弯角。此时微元体4BcD处的金属纤维长度0D为=，氏 (4)式中：r为微元体所处位置半径。因此微元体在切线方向的真应变为白地鲁地篇 OD 70c70棒材在反弯前后中性层的长度不变，即民岛=p-氏最终 勺地志式中：尺。为反弯半径。122 中性层偏移塑性变形应力一应变分析由增量理谢17

16、(6)(7)峪萼卜如)峥*一吉(叫 中等卜妒)可得到塑性变形应力关系：q=去 (9)式中：和盯，分别为微元体切线方向和径向的应力。由塑性体积不变条件得：幻=一，则塑性区等效应力万为庐击厄丽而i丽=孚1I(10)塑性区等效应变i为F=孚厄习瓦i丽=去I勺I(11)塑性变形拉伸区切向应力分布函数M，为=井陪等一 ，压缩区切向应力分布函数y为嘞=斗如悟半+， ，式中：v为中性层偏移量；吼为材料屈服极限；足为矫直弯曲形成的弹性芯半径；五为硬化系数；z为截面纵向纤维至中性层距离。半空间集中力的应力示意图如图3所示，取c(孝，即)为加载区S内表面上一点，而彳0，y，力为固体内的任意点1 91，则距离风为!

17、风=(孝一x)2+(叩一y)2+z22 (14)定义势函数日为日=肛p(善，印)臼d嗣77 (15)式中：Q=zln(po+z)一砌。贝0何=警=肌砌灿(风+z)(16)y=罢=肌础荆，去删77 (17)z方向应力函数盯：为旷去隆z别 郴，边界条件为一 fp(，77)，s域内吒2 1o，贼外当集中力P垂直作用在表面原点上时，!风=(z2+y2+z2)2，并且P=mp(孝，77)d善d，7定义的势函数可简化为lf，：掣：! (19)lf，2-=一2一 LlyJ凸辊对棒材施加的矫直压力P与塑性区径向应力以的关系为q一一芸，专 ，万方数据第9期 马自勇，等：基于弹塑性压力中性层偏移的棒材二辊矫直回弹

18、模型 3023pO I z，二f西，7、了仁，、二 ；二图3半空间集中力的应力示意图Fig3 Stress diagmm of half space concentration123基于弹塑性压力的中性层位置确定棒材矫直过程中在弹性区与塑性区的分界面上有厅=q，即q=压联立消去，得：p=3P2兀Ir， n d、|2r氏+i|n d凡+i一(21)(22)(23)俐竺竺二型 ，胪矿 Q4热爿=孚(盯丁= ，为与材料有关的系数；d为棒材直径；q为屈服极限；B为塑性硬化系数；胛为硬化指数；修正系数77与材料性能及弯凸程度有关，由大量模拟实验回归拟合确定。13内弯矩M分析棒材矫直过程中应力一应变分布如

19、图4所示。I2尺!：一二二陟一 f 厶1一一一7oj夕 久 丫， i f。肇辎7孓l _。 彬 |一。，形一纱芽髟7纠lq慷【卜乱图4棒材矫直应力一应变图Fig4 Stresstrain diagram ofbar straightening棒材断回的仕葸商度z处明断回苋发y：2“再，棒材半径为R，矫直弯曲形成的弹性芯半径为R，此截面的弯矩M为M=二”y(一zV)出+R。，十)2出+F心川2协，月J兄一，y(z+，)也 (25)将式(25)化简得M=尺3蚋1-(f+4)2_+R3q局(f+4)Fi歹了孑+R3吼届arcsin(f+4)+R3q詈所坝3q哆1一(f一4)2j+ R3即岛(一锄、F

20、而+尺3q岛a11csin(f一4)+R3q詈儿 (26)其弯矩比厨表达式为厨：兰a1l一(f+4)2l+二届(f+4)1一(f+4)2+二届arcsin(f+4)+4厂I 42乃+姜口：1一(f一4)2i+ 丁c L J!屈(f一4)扛而+昙愿arcsin(f删+耽 (27)尢 2Vll，8|吼时同巫她愿万方数据3024 中南大学学报(自然科学版) 第47卷2去兄+砉卜一甜专兄邕+静届=去(一五一詈旯卜+(一去旯专+等；吃2击肌南卜+甜去五铋一筹仍一去(一兄+詈兄卜+(一砉五专+等；4 4旯y1 2万_了；y2 2万了：面：尝；弹性极限弯矩M。：盟； M 。 4弹区比f=等；中性层偏移半径比

21、彳=云。将纯弯曲情形所得到的弯矩一曲率关系M=西(r)应用于横向载荷作用下的弯曲【201。根据基本假设2)，肖景荣等21】研究表明：材料弯曲后的卸载过程可以认为是在塑性弯曲相反的方向对被弯曲材料施加一个假想的弹性弯矩M。，并使其产生完全的弹性弯曲变形占。的过程，M。与M相等，即尥=M。由于棒材在进入矫直机前存在一定的原始挠度，根据平面假设知，棒材截面z处有初始应变岛占o=z (28)l一l蛙 一一旦式中：为棒材微段的初始曲率：z为截面纵向纤维至中性层距离。此应变为初始当量应变221。棒材微段弯曲过程示意图如图5所示。根据基本假设1)和初始当量应变定义，棒材微段截面上任意点的弹性弯曲应变占。、加

22、载后应变占及卸载后残余应变晶分别为岛=kz一确z (29)s=盯z一z (30)邻=-z一确z (31)式中：k为弹性弯矩加载时棒材微段中心轴曲率；K为加载后棒材微段中性层曲率；k为卸载后棒材微段中性层曲率，同时规定k，茁，盯。与方向一致时符号相同，方向相反时符号相反。由截面的弯矩平衡可得M。=I。E巳。z。出=日(k一) (32)式中：E为弹性模量；，为棒材断面对y轴的惯性矩。在加载和卸载过程中，占2毛+邰 (33)将式(29)(31)代入式(33)可得：茁2K一+ (34)联立式(32)和(34)及毛=M，可得棒材微段的残余曲率方程为一茜 (35)1 JK，) (ry 一_(d) (e)(

23、a)初始曲率：(b)反弯曲率；(c)残余曲率(与初始曲率反向)；(d)残余曲线(与初始曲率同向)：(e)弹性曲率图5棒材微段弯曲过程示意图Fig5 Bending process diagr锄of bar lIlicro-segment万方数据第9期 马自勇，篁!茎rf塑望丝垦查竺星堡鳖塑鳖堑三堡堑皇回翌堡型堕二二=一。因此，棒材矫直回弹挠度计算模型为M一面碉+万=万一万=露一再式中：为三点弯曲中两端点间的距离。2棒材矫直回弹挠度实验(36)21实验路线及目的在WAWl000万能试验机进行拉伸试验，得到40Cr和42CrMo合金钢2种棒材的室温拉伸应力一应变曲线，利用Origin软件对所得曲线

24、进行线性拟合，为挠度回弹计算模型构建提供参数支持。弯曲试验在10 t液压式压力试验机上进行，通过测量棒材跨中反弯和回弹后的挠度值来验证挠度回弹计算模型的准确性。22室温单向拉伸实验通过拉伸实验可以获得材料的屈服强度、伸长率和抗拉强度等参数，这些参数是影响棒材回弹的重要因素。本文以40Cr和42CrMo棒材为例进行拉伸实验，得到其应力一应变曲线如图6所示。传统双线性拟合都是将屈服极限吼到最大抗拉强度盯。之间的数据进行拟合，对于存在屈服平台的材料来说，拟合精度很差，如图6所示；而从矫直原理考虑，当弯曲条材弹塑性变形深度达到全断面高度的80后不产生裂纹时，任何金属条材都可以得到矫直【231，因此，棒

25、料矫直变形最大应变为5占。，&为弹性极限应变。利用0啦in对40Cr和42CrMo的拉伸实验数据进行函数拟合，其中拟合的数据取自屈服点至5之间，拟合发现：40Cr和42CrMo进入塑性变形后，采用新线性拟合的校正决定系数分别为0936 1和0983 7，拟合效果较好，因此，线性硬化关系万=正+五E(瓦一日)能描述40Cr和42CrMo棒材的应力一应变关系。2种材料线性拟合结果见表1。23弯曲实验采用棒材弯曲回弹测试对棒材挠度进行测试，测试方案如图7所示，百分表用于测定棒材跨中的位移。40Cr压下量分别为10，15，20和30 mm；42CrMo压下量分别为10，15和20 mm。每次压下速度为

26、10删州111in，压下动作完毕后，用百分表测量跨中位移3次并记录取平均值，测量记录完成就开始卸载，卸载速度为100舢Ilmin。应变(a)40Cr；(b)42CrMo图6应力一应变关系曲线Fig6 Relationship cues of stresss仃ain表1不同材质的力学性能Thble l Mechanical pmpert)r of diffbrent material pammeters尸棒材图7挠度测量示意图Fig7 Diagram of denection measurement对棒材进行弯曲试验，载荷一位移曲线如图8所示，在开始的弹性阶段，载荷和位移呈明显的线性变化关系。随

27、着载荷的增加，棒材跨中位移不断增加，当载荷增大到18 kN左右时，棒材开始屈服，之后棒万方数据3026 中南大学学报(自然科学版) 第47卷位移m111图8载荷一位移曲线Fig8 Loaddisplacement curves材跨中位移的增加速度迅速增大，而载荷的增加速度则减缓。当位移增加至30 mm和20 mm时，开始卸载。将实验所得数据整理计算即可得到棒材弯曲的回弹值，结果见表2。3结果分析与讨论31 中性层偏移对回弹的影响二辊矫直过程中性层偏移会改变棒材截面应力一应变的状态，从而改变棒材内弯矩的计算，最终影响回弹和棒材的矫后直线度。40Cr和42CrMo的M一4一f曲线如图9所示。表2基

28、于试验所测得的回弹值1able 2 Springback measured by test日丑跹蒉L淤堡嗵趔三(a)40Cr；(b)42CrMo1彳=0时砑一f曲线：2彳0时面一f曲线；34一f曲线。图9棒材的面一4一f曲线Fig9面一彳一fcll】ves ofbars目釜也弁滁堡啮越耳万方数据第9期 马自勇，等：基于弹塑性压力中性层偏移的棒材二辊矫直回弹模型 3027Mf曲线表现的是棒材弹塑性变形段某一截面在不同弯曲程度时的弯矩变化规律；4一f曲线表现的是棒材弹塑性变形段某一截面在不同弯曲程度时中性层的偏移情况。由图9可知：在弹区比f07时，中性层偏移值较小且对弯矩比的影响也较小，最大相对误

29、差约为74，这是因为此时棒材塑性变形深度值小，中性层偏移现象不明显。棒材弹塑性变形深度达到50即可实现矫直241，同时，由矫直原理可知，当棒材弹塑性变形深度达到80后不产生裂纹，则任何金属可以矫直，因此，棒材矫直所需要施加的弹区比f为0205，在此范围内中性层偏移对弯矩比的影响较大，40Cr与42CrMo最小相对误差分别达1265和1176，见表3；中性层偏移对棒材回弹的影响也较大，随着弯曲程度增大，回弹值的误差增大，在f=01时，40Cr与42CrMo相对误差分别为4159和1848，如图10所示。从图10可见：在弹区比f07且考虑材料硬化时，中性层偏移对棒材的回弹挠度影响较小，最大相对误差

30、均在75左右，而在矫直所需弹区比=02o5时，对棒材的回弹挠度影响较大，随着弹区比f的减小而增大。这是因为中性层偏移是由塑性变形导致，弹区比越小，棒材弯曲程度越大，棒材塑性变形层深度越大，中性层偏移也愈加明显，最终对回弹挠度的影响也愈大。32材料的力学性能对回弹的影响实验所用的是材料分别是40Cr和42CrMo，它们的力学性能参数见表1，图12所示为2种棒材弯曲回弹后的形状。图12所示为不同材料的挠度与弹区比关系曲线。从图12可知：在考虑材料硬化条件下，考虑中性层偏表3中性层偏移对弯矩的影响Table 3 Inlpact of neu仃al layer o凰et on momentE三吣型器慧

31、E(a)40Cr；(b)42CrMo1彳=0：2彳0。图10回弹挠度与弹区比的关系Fig10 Relationship between of denection and elastic region value万方数据3028 中南大学学报(自然科学版) 第47卷兰E吣蜊嚣教巨(a)40Cr；(b)42CrMo图11弯曲回弹后形状Fig11 Sh叩e ofbending and springback弹区比f140Cr回弹理论值；240Cr回弹实验值；342CrMo回弹理论值；442CrMo回弹实验值。图12不同材料的挠度与弹区比关系曲线Fig12 Relation curve 0f denec

32、tion and elastic region Valuei11 dif色rent materials移的棒材矫直挠度回弹计算理论值与实验值基本一致。它们随着弹区比f增大而减小，其中，42CrMo理论值与实验值的相对误差很小，在212363范围内；而40Cr理论值与实验值相对误差为322645，这主要是由于双线性拟合模型更接近42CrMo的应力一应变曲线，40Cr的应力一应变存在屈服平台，与42CrMo相比，其双线性硬化拟合精度较低。同时，从图12可以看出：42CrMo回弹挠度比40Cr棒材的大，这说明金属材料性能对回弹挠度影响较大。这是因为在相同的弹区比f时，屈服强度高的42CrMo弯曲所需

33、要的力矩大，最终卸载后回弹弯矩也大【251，因此42CrMo棒材回弹挠度比40Cr的大。3-3曲率K对回弹的影响挠度与曲率关系曲线如图13所示。从图13可以看出：随着曲率K的增大，棒材回弹挠度越大；这是因为曲率K增大，反弯半径就越小，棒材弯曲程度愈大，进而回弹挠度增大。同时，2种材料的回弹挠度理论值与实验值基本一致，40Cr与42Crlo最大相对误差分别为645和363，低于不考虑中性层偏移时的3165和991，表明式(30)计算精度较高。三妄吣型遥：墨L菩弯曲盐率彤l一40cr回弹理论值；240Cr回弹实验值：342CrMo回弹理论值；4-叫2CrMo回弹实验值。图13挠度与曲率关系曲线Fi

34、g13 Cu“e of denection and curvature4结论11建立了棒材矫直回弹挠度计算模型，可对高精度棒材二辊矫直过程中回弹挠度进行计算和分析。实验证明了该模型的正确性，对提高棒材矫直精度及完善其弯曲弹塑性变形理论提供了理论参考。2)基于二辊矫直变形特点和矫直理论，得到棒材矫直所需要的最大应变不超过5的结论，提出了在应变至5之间进行拟合的新双线性矫直拟合方法。与传统的拟合方法相比，40Cr与42CrMo的拟合校正决定系数由0833 1和0946 5提高到0936 l和0983 7，拟合精度分别提高了1030和372。3)棒材回弹挠度不仅与矫直过程棒材受力状态、材料力学性能和

35、曲率有关，还与棒材同矫直辊的接触位置有关。棒材回弹挠度随着中性层偏移值、材料屈万方数据第9期 马自勇，等：基于弹塑性压力中性层偏移的棒材二辊矫直回弹模型 3029服强度和曲率增大而增大。40Cr与42Crlo 2种棒材最大相对误差分别为645和363，低于不考虑中性层偏移和材料硬化时的3165和991。4)对于矫直过程弹区比f07的棒材，中性层偏移对回弹影响很小，工程分析时可忽略；但是对矫直过程弹区比f必须满足的0205范围内的棒材，中性层偏移量对回弹的影响较大，分析时不可忽略。参考文献：1】 王云棒材二辊矫直过程及关键技术研究D秦皇岛：燕山大学机械工程学院，2013：3-4wANG YuSt

36、udy on pmcess aIld key tecllllology of t、】l，omllstra噜htening【D】QiIlllllangdao：Y妇shan uniVersi哆School ofMechanic“Engineering，2013：342】 崔甫矫直技术与理论的新探索【M】北京：冶金工业出版社，2010：89cUI Fu111e New Exploration of stm噜hteIling technology alldtheoryMBe巧ing：Metallurgical Indusny Press，2010：89【3 LI Heng，SHI Kaipeng，YA

37、NG He，et a1Springback law ofthinwalled 606lT4 Alalloy tube upon bendingJ1hnsactionsof Nonfe丌Dus Metals Society of China，2012，22(suppl 2)：3573634】 岳睿带状金属板材弯曲的回弹J】塑性工程学报，2009，1 6(2)：2933YUE RuiStudy on the springback for naITow steel sheet【JJo啪al of Plasticity Engineering，2009，l 6(2)：2932【5 张冬娟，崔振山，李玉

38、强，等宽板大曲率半径纯弯曲回弹量理论分析J工程力学，2006，23(10)：7781zHANG Don卣uan，CuI zhenshan，LIYuqiang，et a1Thesprin曲ack of wide metal sheet a舭r large radius pure bending【J】Engineering Mech锄ics，2006，23(10)：77816 GU R J，YANG H，zHAN M，et a1Research on the springbackof thin walled tube NC bending based on the numericalsimulat

39、ion of the whole pmcessJComputational MaterialsScience，2008，42(4)：5375497 ZHENG Q，AOYAMA TSHIMIZU T，et a1Expe“mental andnumerical analysis of springback behaVior under eleVatedtemperatIlres in micro bending assisted by resistance heatingJ】Procedia Engineering，2014，81：1481一14868 刘婧瑶，唐承统，宁汝新薄壁管纯弯曲塑性成形

40、分析及回弹计算【J塑性工程学报，2009，16(2)：59，14LIU Qianyao，TANG Chen班ong， Nn、IG Ruxin Plasticfo肌ation analysis and springback calculating of血in-wall mbepurebendingJ】Jo啪al of Plasticity Engineering，2009，16(2)：5914【9】 BURANATHITI T，CAO JianAn e钎色ctiVe ana坼ical model for【10】12】13【14】【15】16】17【18【19】springback pIedict

41、ion in sn_aight nanging processJ】IIltemational Jo啪al of Material alld Pmduct Technolo既2004，21(123)：137153贾美慧，唐承统不锈钢管材弯曲成形回弹预测模型研究J】北京理工大学学报，2012，32(9)：910914JIA Meihui，TANG ChengtongSpringback prediction model ontube bending fo册ation of stainless steel【J】1协sactions ofBe日ing hlstitIlte ofTechllolo蹦20

42、12，32(9)：9lo-914张深，吴建军空间弯管的回弹预测【J_航空学报，201l，32(5)：953960zHANG shen，、Vu JianjllllSpringback prediction of nonplanar tIlbe bendingJ】Acta Aeronauticaet Astmnautic Sinica，20 Il，32(5)：953959ZHU Y X，LIU Y L，LI H P，et a1Springback prediction forrotarydmw bending of rectallgular H96 tube b鹊ed on isonDpic，In

43、ixed and Y0shida-Uemori锕o-su血ce hardeIling models【JMaterials and DesigIl，2013，47：200_209刘海燕，金霞板料成形的回弹预测方法研究【J机械制造与研究，2008，37(6)：40-44LIU Haiyan，JIN XiaSpringback prediction method research ofsheer metal foming【JJo啪aI of Machine DesigIl and Research，2008，37(6)：40一44潘昌平管材弯曲回弹及弯管机控制系统的研究D长春：吉林大学机械科学与工程

44、学院，2008：1725PAN Changping Research on the bending springback ofminwalled tube and the con仃ol system of tube bendingmchineDChangchun：Jilin UniVersityC01lege of MechanicalScience and Engineering，2008：1 725官英平，张庆，赵军中性层内移对弯曲回弹的影响J锻压技术，2007，32(2)：2627GUAN Yingping，ZHANG Qing，ZHA0 Jun1nfluence ofneutral l

45、ayer inside displacement on bending spring_back【J】Fo唱ing&Stamping Tbchnology，2007，32(2)：2627陈石保金属圆棒料弯曲回弹的探讨J锻压技术，199l(4)：3840CHEN ShibaoThe discussion on bending springback of metalbarJ】Fo唱ing&Stamping 1_echnology199l(4)：3840俞汉清，陈金德金属塑性成形原理M北京：机械工业出版社1999：12卜122YU Hanqing，CHEN JindeThe principle of

46、metal plasticfomingMBe日ing：china Machine Press，1999：12卜122张伟玮，韩聪，谢文才，等管材弯曲回弹对内高压成形的影响及补偿方法J】哈尔滨工业大学学报，2014，46(7)：3639ZHANG Weiwei，HAN Cong，XIE Wencai，et al_T11e efrect ofspringback of CNC bending on hydrofomed sub-仃ame andcompensation methodsJ】 Jotlnlal of Harbin Institute ofTeclIlology2014，46(7)：36

47、39JOHNs0N K L，徐秉业，罗学富接触力学【M】北京：高等万方数据3030 中南大学学报(自然科学版) 第47卷一一教育出版社，1992：5卜58 【23】崔甫矫直原理与矫直机械M北京：冶金工业出版社，JOHNSON K L， XU Bingye， LUO Xuem Contact 2007：89mechanicsM】BeUing：Hi曲er Education Press，l 992：5卜58 cul FuS仃ai曲tening pmciples锄d stm谵htening mchineM20余同希，章亮炽塑性弯曲理论及其应用【M】北京：科学出 Be日ing：Metallu唱ical Indus仃y Press，2007：89版社，1992：4344 【24闫长林m8