薄壁管滾珠旋壓力計(jì)算模型

趙春江，王琛，耿明超

（1.太原科技大學(xué)重型機(jī)械教育部工程研究中心，山西太原030024；2.中國(guó)電子科技集團(tuán)公司第12研究所，北京100044）

薄壁管滾珠旋壓力計(jì)算模型

趙春江1，王琛1，耿明超2

（1.太原科技大學(xué)重型機(jī)械教育部工程研究中心，山西太原030024；2.中國(guó)電子科技集團(tuán)公司第12研究所，北京100044）

給出了一種計(jì)算三維滾珠旋壓力的解析模型。根據(jù)滾珠軌跡對(duì)空間坐標(biāo)進(jìn)行旋轉(zhuǎn)，簡(jiǎn)化了滾珠與工件的空間關(guān)系表達(dá)式；基于空間解析幾何理論和旋壓件為半無(wú)限體的假設(shè)，得到了滾珠與管坯的接觸邊界方程，并在各坐標(biāo)平面上進(jìn)行投影得到二維曲線，通過(guò)對(duì)曲線的積分計(jì)算各坐標(biāo)平面內(nèi)的接觸區(qū)投影面積。變形區(qū)單位壓力通過(guò)鋼球壓入空間半無(wú)限體時(shí)成形區(qū)的平均壓力近似相等的假設(shè)給出，計(jì)算各旋壓力分量和總旋壓力；通過(guò)對(duì)旋壓力軸向分量的理論計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比，表明給出的計(jì)算模型的準(zhǔn)確性。

兵器科學(xué)與技術(shù)；滾珠旋壓；旋壓力；坐標(biāo)旋轉(zhuǎn)；解析模型；實(shí)驗(yàn)對(duì)比

0 引言

滾珠旋壓是一種強(qiáng)力精密旋壓，采用滾珠對(duì)管材進(jìn)行連續(xù)局部塑形成形［1-3］。Rotarescu［4］對(duì)滾珠旋壓的過(guò)程及其旋壓力的計(jì)算進(jìn)行了系統(tǒng)的理論推導(dǎo)和數(shù)值模擬，借鑒了擠壓成型變形區(qū)平均壓力值的計(jì)算方法，并對(duì)接觸面積的計(jì)算做了簡(jiǎn)化，但其推導(dǎo)過(guò)程復(fù)雜，且變形區(qū)接觸面積誤差較大；李茂盛等［5］忽略了滾珠軸向進(jìn)給所需的壓痕力，利用平面應(yīng)變狀態(tài)下圓弧形沖壓頭壓入半無(wú)限體時(shí)平均接觸壓力與滾珠旋壓接觸平均力相等的假設(shè)，得出了滾珠旋壓力的計(jì)算公式；文獻(xiàn)［6］按照管件的均勻變形給出了滾珠旋壓力簡(jiǎn)便的工程計(jì)算方法；Zhang等［7］對(duì)銅管內(nèi)螺紋滾珠旋壓折疊缺陷及其金相組織特性進(jìn)行了分析；Jiang等［8-9］對(duì)帶肋薄壁管的滾珠旋壓過(guò)程進(jìn)行了詳細(xì)的分析和數(shù)值模擬，研究并討論了縱向內(nèi)肋薄壁管材滾珠旋壓的變形準(zhǔn)則并對(duì)滾珠旋壓工藝參數(shù)之間的非線性建立了基于神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)模型，其成果較大地推進(jìn)了滾珠旋壓工藝和工模具設(shè)計(jì)理論的研究。馬振平等［10］對(duì)滾珠旋壓加工技術(shù)進(jìn)行了系統(tǒng)的總結(jié)。

本文給出了一種滾珠旋壓力的計(jì)算方法，相對(duì)于已有的計(jì)算方法，該方法基于對(duì)旋壓件為半無(wú)限體的假設(shè)，并按照滾珠旋壓軌跡螺旋升角將坐標(biāo)進(jìn)行旋轉(zhuǎn)，簡(jiǎn)化了工件與滾珠的空間關(guān)系表達(dá)方式；根據(jù)滾珠球面和旋件的空間解析幾何關(guān)系，結(jié)合滾珠當(dāng)前位置與上一圈的旋壓軌跡，得到了滾珠與管坯接觸區(qū)的邊界曲線表達(dá)式，從而得到了接觸區(qū)在各自對(duì)應(yīng)坐標(biāo)內(nèi)的投影面積，引入鋼球壓入空間半無(wú)限體時(shí)成形區(qū)的平均壓力近似相等的假設(shè)，從而計(jì)算各旋壓力分量；最后通過(guò)理論計(jì)算與實(shí)驗(yàn)數(shù)據(jù)的比較說(shuō)明了本文計(jì)算方法的準(zhǔn)確性。

1 工藝描述

圖1所示為滾珠旋壓原理示意圖，支撐管、外支撐圈、圓錐模環(huán)和滾珠組成了滾珠旋壓模具，旋壓模具套在工件外壁，模具和工件相對(duì)轉(zhuǎn)動(dòng)同時(shí)產(chǎn)生沿軸向的相對(duì)進(jìn)給運(yùn)動(dòng)，則套在心軸上的工件與鋼球接觸并受壓產(chǎn)生塑性變形，達(dá)到壁厚減薄的目的。

2 旋壓力空間解析模型

2.1 基本假設(shè)與坐標(biāo)設(shè)置

由圖1的滾珠旋壓工藝原理可知，在滾珠旋壓過(guò)程中，由于滾珠的曲率遠(yuǎn)大于管坯的曲率，且旋壓過(guò)程是宏觀的逐點(diǎn)擠壓變形，變形接觸區(qū)很小，且接觸區(qū)內(nèi)部變形應(yīng)力對(duì)管坯其他部分幾乎沒(méi)有影響，所以把管坯假設(shè)成半無(wú)限體。

為簡(jiǎn)化解析模型的推導(dǎo)過(guò)程，建立動(dòng)坐標(biāo)系與固定坐標(biāo)系。滾珠球心設(shè)為動(dòng)坐標(biāo)系的原點(diǎn)O，沿著滾珠軸向進(jìn)給螺旋升角的方向建立動(dòng)坐標(biāo)系x軸，過(guò)球心與管件軸線的垂直相交線作為動(dòng)坐標(biāo)系z(mì)軸，垂足記作O′，動(dòng)坐標(biāo)系y軸通過(guò)球心且垂直于Oxz坐標(biāo)平面。以O(shè)′點(diǎn)作為固定坐標(biāo)系的原點(diǎn)，固定坐標(biāo)系z(mì)′軸與動(dòng)坐標(biāo)線z軸方向相同。在管坯軸線上，作固定坐標(biāo)系y′軸，垂直于O′y′z′坐標(biāo)平面并過(guò)O′點(diǎn)建立固定坐標(biāo)系x′軸。各坐標(biāo)軸方向如圖2所示。

圖1 滾珠旋壓示意圖Fig.1 Schematic diagram of ball spinning

圖2 解析模型坐標(biāo)示意圖Fig.2 Schematic diagram of coordinates of analytic model

2.2 接觸面積的計(jì)算

根據(jù)圖2可得滾珠（球）的坐標(biāo)方程為

式中：Rb為滾珠的半徑。

由于滾珠的曲率遠(yuǎn)大于管坯的曲率，且旋壓過(guò)程是宏觀的逐點(diǎn)擠壓變形，變形接觸區(qū)很小，假設(shè)管坯外表面是半無(wú)限平面，滾珠與管坯上一道次旋壓所形成的軌跡在Oyz平面上的投影視為圓弧，其方程為

式中：f為軸向旋模比；n為滾珠個(gè)數(shù)；θ為滾珠旋壓軌跡螺旋升角。

滾珠與管坯接觸所形成的空間曲面，在動(dòng)坐標(biāo)系各坐標(biāo)平面上的投影面積計(jì)算如下。

2.2.1 沿x軸方向的投影面積

滾珠與管坯的接觸面積在Oyz坐標(biāo)平面（x軸方向）的投影面積如圖3所示的陰影部分，x軸方向面積由滾珠球面、上一道次旋壓軌跡和管件在Oyz坐標(biāo)平面內(nèi)的投影共同確定。其中：AB為滾珠與管坯表面交線的投影，實(shí)際為一段曲線，因假設(shè)管坯表面為半無(wú)限平面，故AB視為直線段。曲線為滾珠與上一道次旋壓軌跡的相交線在Oyz坐標(biāo)平面上的投影。曲線為滾珠圓的一部分，是滾珠壓入管坯內(nèi)所形成接觸區(qū)的邊界線在Oyz坐標(biāo)平面內(nèi)的投影。

圖3 x軸方向投影面積Fig.3 Projection area in x direction

滾珠在Oyz坐標(biāo)平面的投影坐標(biāo)方程為

上一道次旋壓軌跡在Oyx坐標(biāo)平面的投影坐標(biāo)方程為

根據(jù)管坯外表面為半無(wú)限體平面的假設(shè)，管件在Oyx坐標(biāo)平面內(nèi)的投影方程為

式中：tp為減薄量。

由（3）式、（4）式、（5）式，可得圖3中A點(diǎn)、B點(diǎn)、C點(diǎn)、E點(diǎn)的坐標(biāo)分別為，，則x軸方向的投影面積

2.2.2 沿y軸方向的投影面積

滾珠與管坯的接觸面積在Oxz坐標(biāo)平面（y軸方向）的投影面積如圖4所示的陰影部分，y軸方向面積由滾珠球面、上一道次旋壓軌跡和管件在Oxz坐標(biāo)平面上的投影共同決定。其中：AB為滾珠與管坯表面交線的投影，曲線為空間曲線在 Oyz坐標(biāo)平面上投影的一部分，曲線為滾珠球面圓的一部分，曲線為空間曲線在Oxz坐標(biāo)平面投影的一部分。

圖4 y軸方向投影面積Fig.4 Projection area in y direction

圖4中，AB的方程表達(dá)式為

AF的方程表達(dá)式為

根據(jù)（6）式～（9）式可得圖4的A點(diǎn)、B點(diǎn)、F點(diǎn)、G點(diǎn)在Oxz坐標(biāo)平面上的坐標(biāo)分別為A（0，tp-Rb）、，則y軸方向的面積

2.2.3 沿z軸方向的投影面積

滾珠與管坯的接觸面積在Oxy坐標(biāo)平面（z軸方向）的投影面積如圖5所示的陰影部分，z軸方向面積是由滾珠球面、上一道次旋壓軌跡和管件在Oyz坐標(biāo)平面內(nèi)的投影共同確定。其中：曲線為滾珠表面與上一道次軌跡交線在Oxy坐標(biāo)平面的投影，曲線為滾珠球面與管件外表面交線Oxy坐標(biāo)平面的投影，直線AC為滾珠壓入管坯內(nèi)所形成接觸區(qū)的邊界線在Oxy坐標(biāo)平面內(nèi)的投影。

圖5 z軸方向投影面積Fig.5 Projection area in z direction

AC的方程為

根據(jù)（10）式～（12）式可得出A點(diǎn)、B點(diǎn)、C點(diǎn)在Oxy坐標(biāo)平面上的坐標(biāo)分別為A（0，-Rb）、、.其中：，.則z軸方向的面積

2.3 旋壓力的計(jì)算

在滾珠旋壓工藝中，滾珠與管坯接觸面上的壓力在x軸、y軸、z軸3個(gè)方向上的分力可看作指向球心的均布力，由2.2節(jié)的Sx、Sy、Sz來(lái)計(jì)算三向旋壓分力，即

式中：p為平均接觸壓力。

本文采用經(jīng)實(shí)驗(yàn)和模擬驗(yàn)證可用的接觸平均壓力p計(jì)算公式［5］為

式中：σs為抗拉屈服強(qiáng)度。

由于上述計(jì)算模型的推導(dǎo)基于動(dòng)坐標(biāo)進(jìn)行，為了方便分析及與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，需將上述計(jì)算所得的各項(xiàng)旋壓力分量轉(zhuǎn)換到定坐標(biāo)系中，切向分力、軸向分力、徑向分力轉(zhuǎn)換方程分別為

3 計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比

針對(duì)上述推導(dǎo)結(jié)果，應(yīng)用VC計(jì)算機(jī)語(yǔ)言編寫了三向滾珠旋壓力的計(jì)算程序。為了與文獻(xiàn)［4］的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，選取文獻(xiàn)［4］中的材料與工藝參數(shù)。根據(jù)文獻(xiàn)［10-11］計(jì)算所用的金屬材料的機(jī)械性能參數(shù)，其計(jì)算結(jié)果和實(shí)驗(yàn)數(shù)據(jù)如圖6所示。由于實(shí)驗(yàn)條件的限制，只能對(duì)軸向旋壓力分量進(jìn)行測(cè)量，因此圖6給出的計(jì)算和實(shí)驗(yàn)數(shù)據(jù)只針對(duì)總旋壓力的軸向分量。為了充分顯示旋壓力隨進(jìn)給速度變化的非線性特性，如表1所示，第3組、第4組數(shù)據(jù)中取大跨度的進(jìn)給速度1～10 mm.

圖6 理論計(jì)算和實(shí)驗(yàn)數(shù)據(jù)比較Fig.6 Comparison between calculated and experimental data

表1 實(shí)驗(yàn)工藝數(shù)據(jù)表Tab.1 Experimental process data table

由圖6可見(jiàn)：針對(duì)各組工藝參數(shù)總旋壓力的軸向分量隨著旋模比的增加而增大，總旋壓力隨著旋模比的對(duì)數(shù)呈指數(shù)增長(zhǎng)趨勢(shì)；相同條件下，旋壓力隨著減壁量和管坯直徑的增加而增大；計(jì)算結(jié)果和實(shí)驗(yàn)數(shù)據(jù)的變化趨勢(shì)一致；在旋模比較小時(shí)，實(shí)驗(yàn)數(shù)據(jù)和計(jì)算結(jié)果吻合良好；而在旋模比較大時(shí)，實(shí)驗(yàn)數(shù)據(jù)明顯大于計(jì)算結(jié)果。原因在于旋模比增大會(huì)引起滾珠之間的金屬堆積和隆起，從而增大滾珠旋進(jìn)阻力；且大旋模比使變形速率增加，會(huì)提高工件的變形阻力，從而使旋壓力數(shù)據(jù)加大。另外，工件與芯棒的摩擦阻力會(huì)阻止工件的延伸，也會(huì)增大軸向旋壓力。

綜上所述，本文給出的解析模型能夠準(zhǔn)確反映旋壓力的變化趨勢(shì)，小旋模比時(shí)可得到高精度的旋壓力計(jì)算值。在大旋模比情況下，應(yīng)考慮實(shí)際情況對(duì)旋壓力進(jìn)行適當(dāng)修正，可保證滾珠旋壓力計(jì)算模型的精確性和可靠性。

4 結(jié)論

1）本文根據(jù)滾珠軌跡對(duì)空間坐標(biāo)的旋轉(zhuǎn)關(guān)系，基于空間解析幾何理論和旋壓件，為半無(wú)限體的假設(shè)給出了一種計(jì)算三維滾珠旋壓力的計(jì)算模型。給出了其變形區(qū)單位壓力通過(guò)鋼球壓入空間半無(wú)限體時(shí)，成形區(qū)的平均壓力近似相等的假設(shè)，相對(duì)于已有的計(jì)算模型，簡(jiǎn)化了滾珠與工件的空間關(guān)系表達(dá)式，且具有較高的精確性。

2）采用現(xiàn)有文獻(xiàn)［4］的實(shí)驗(yàn)材料和工藝參數(shù)，應(yīng)用本文給出的計(jì)算模型進(jìn)行總旋壓力的軸向分量的計(jì)算，通過(guò)對(duì)比計(jì)算結(jié)果和實(shí)驗(yàn)數(shù)據(jù)，發(fā)現(xiàn)二者的變化趨勢(shì)一致。本文給出的并考慮了大旋模比影響的計(jì)算模型具有較高的精度和可靠性。

（

）

［1］Tang Y，Chi I Y，Chen J C，et al.Experimental study of oil-filled high-speed spin forming micro-groove fin-inside tubes［J］.International Journal of Machine Tools and Manufacturing，2007，47（7/8）：1059-1068.

［2］Hua F A，Yang Y S，Zhang Y N，et al.Three-dimensional finite element analysis of tube spinning［J］.Journal of Materials Processing Technology，2005，168（1）：68-74.

［3］Wong C C，Dean T A，Lin J.A review of spinning，shear forming and flow forming processes［J］.International Journal of Machine Tools&Manufacture，2003，43（14）：1419-1435.

［4］Rotaresuc M I.A theoretical analysis of tube spinning using balls［J］. Journal of Materials Processing Technology，1995，54（1/2/3/4）：224-229.

［5］李茂盛，康達(dá)昌，張士宏，等.滾珠旋壓工藝中成形區(qū)接觸壓力的分析計(jì)算［J］.材料科學(xué)與工藝，2004，12（2）：125-128. LI Mao-sheng，KANG Da-chang，ZHANG Shi-hong，et al.Research on contacting pressure in ball-spinning［J］.Materials Science and Technology，2004，12（2）：125-128.（in Chinese）

［6］徐洪烈.強(qiáng)力旋壓技術(shù)［M］.北京：國(guó)防工業(yè)出版社，1984. XU Hong-lie.Power spinning technology［M］.Beijing：National Defense Industry Press，1984.（in Chinese）

［7］Zhang G L，Zhang S H，Zhang B，et al.Analysis on folding defects of inner grooved copper tubes during ball spin forming［J］. Journal of Materials Processing and Technology，2007，184（1）：393-400.

［8］Jiang S Y，Ren Z Y.Analysis of mechanics in ball spinning of thin-walled tube［J］.Chinese Journal of Mechanical Engineering，2008，21（1）：25-30.

［9］Jiang S Y，Ren Z Y，Xue K M，et al.Application of BPANN for prediction of backward ball spinning of thin-walled tubular part with longitudinal inner ribs［J］.Journal of Materials Processing Technology，2008，196（1/2/3）：190-196.

［10］馬振平，張濤.滾珠旋壓成形技術(shù)［M］.北京：冶金工業(yè)出版社，2011. MA Zhen-ping，ZHANG Tao.Ball spinning forming technology［M］. Beijing：Metallurgical Industry Press，2011.（in Chinese）

［11］機(jī)械設(shè)計(jì)手冊(cè)編委會(huì).機(jī)械設(shè)計(jì)手冊(cè)：第1卷［M］.北京：機(jī)械工業(yè)出版社，2004：3-313. Mechanical Design Manual Editorial Board.The mechanical design manual，volume 1［M］.Beijing：China Machine Press，2004：3-313.（in Chinese）

Calculation Model of Ball Spinning Force of Thin-walled Tubes

ZHAO Chun-jiang1，WANG Chen1，GENG Ming-chao2
（1.The Heavy Mechanical Engineering Research Center of Ministry of Education，Taiyuan University of Science and Technology，Taiyuan 030024，Shanxi，China；2.The 12th Research Institute，China Electronic Technology Group Corporation，Beijing 100044，China）

An analytical model is presented to calculate the 3D ball spinning force.The spatial relation expression of the ball and workpiece is simplified through the rotation of spatial coordinates according to the track of the ball.The contact boundary equations of the ball and the workpiece are obtained based on the theory of space analytic geometry and the assumption that the spinning workpiece is a semi-infinite body，and the 2D curves are obtained by projecting on the coordinate planes.The projection of spatial contact area on each coordinate plane is calculated.The unit pressure in the deformation zone is given based on the assumption that the average pressure is approximately equal to that while a rigid ball is pressed into the space of semi-infinite body，and then the total spinning force and its components could be calculated.The comparison between theoretical results and experimental data shows that the proposed model is accurate.

ordnance science and technology；ball spinning；spinning force；rotation of coordinate；analytical model；experimental comparison

TG316

A

1000-1093（2015）05-0928-05

10.3969/j.issn.1000-1093.2015.05.024

2014-12-29

國(guó)家自然科學(xué)基金項(xiàng)目（51375325）；山西省自然科學(xué)基金項(xiàng)目（2014011025-1）

趙春江（1975—），男，副教授。E-mail：zhaochj75@163.com