趙雅麗，惠 蒲，崔珊珊

(西安航天發(fā)動(dòng)機(jī)有限公司，陜西 西安 710100)

0 引言

噴管類(lèi)零件是液體火箭發(fā)動(dòng)機(jī)中的常用和重要零件，其成形質(zhì)量直接影響發(fā)動(dòng)機(jī)的性能和安全。由于鈑金件生產(chǎn)成本低、 生產(chǎn)效率高，可以成形出空間結(jié)構(gòu)復(fù)雜的零件，因此這類(lèi)零件大多采用鈑金成形的工藝方法。但是由于其型面為曲母線，成形過(guò)程是一個(gè)具有幾何、材料和接觸高度非線性的復(fù)雜力學(xué)問(wèn)題，影響因素很多，生產(chǎn)實(shí)際中很難精確地控制材料的流動(dòng)，因此成形過(guò)程中會(huì)產(chǎn)生各種各樣的缺陷，影響零件的幾何精度、表面質(zhì)量和力學(xué)性能[1]。在預(yù)防成形缺陷方面，常采用理論分析和實(shí)踐驗(yàn)證相結(jié)合的方法。

目前國(guó)內(nèi)外此類(lèi)零件的工藝設(shè)計(jì)大多采用理論分析計(jì)算和實(shí)踐工程經(jīng)驗(yàn)相結(jié)合的方法，對(duì)于零件的成形缺陷預(yù)防具有明顯的局限性。有限元模擬技術(shù)的發(fā)展，為分析復(fù)雜板料的成形過(guò)程提供了比較有效的工具[2-6]。但是目前應(yīng)用有限元模擬分析計(jì)算此類(lèi)零件成形過(guò)程的文獻(xiàn)資料鮮有報(bào)道。

本文針對(duì)某型號(hào)噴管內(nèi)壁，利用PAM-STAMP鈑金數(shù)值模擬軟件對(duì)其成形過(guò)程進(jìn)行模擬仿真，直觀了解和掌握了零件成形過(guò)程中的變形規(guī)律[7-8]，找出了成形缺陷存在的原因，提出了合理的改進(jìn)措施，并生產(chǎn)出了滿足設(shè)計(jì)圖紙要求的零件。

1 零件技術(shù)要求及控制難點(diǎn)

1.1 零件尺寸及技術(shù)要求

該零件如圖1所示。其材料為1Cr18Ni9Ti，原材料壁厚2.5 mm，成形后要求最小壁厚為2.2 mm，零件型面與樣板間隙≤0.3 mm。

1.2 成形缺陷

該零件采用扇形板料對(duì)縫焊接成錐筒后使用模具液壓成形的工藝方案[9]，錐形坯料如圖2所示。成形后存在大端周圈起皺及小端喉部下方周向環(huán)狀突起的缺陷，如圖3所示。檢查發(fā)現(xiàn)與樣板間隙約0.5 mm，無(wú)法滿足設(shè)計(jì)圖紙≤0.3 mm的要求。

1.3 難點(diǎn)分析

由圖1可見(jiàn)，該零件形狀有以下幾個(gè)特點(diǎn)：①大、小端直徑比大；②型面曲率變化大。零件由錐筒液壓成形，大端主要完成收縮變形，中間部位為脹形，小端喉部為收縮變形，出口為脹形。成形過(guò)程既有脹形又有擠壓、拉伸和收縮，材料流動(dòng)和變形狀況非常復(fù)雜，難以通過(guò)理論計(jì)算掌握其成形過(guò)程中應(yīng)力應(yīng)變情況及材料流動(dòng)規(guī)律。另外，由于零件的成形過(guò)程是在封閉的模具型腔內(nèi)，無(wú)法直觀看到零件的變形過(guò)程，因此只能憑借經(jīng)驗(yàn)對(duì)工藝進(jìn)行改進(jìn)，都沒(méi)有達(dá)到理想的效果。

圖1 零件圖 圖2 坯料圖Fig.1 Part drawing Fig.2 Blank drawing

圖3 零件缺陷示意圖Fig.3 Part defect drawing

2 成形過(guò)程有限元仿真

板料沖壓成型仿真的求解方法有動(dòng)力顯式和靜力隱式兩種算法。動(dòng)力顯式算法計(jì)算時(shí)間步長(zhǎng)小，計(jì)算時(shí)間短，不涉及平衡迭代過(guò)程和解的收斂性問(wèn)題，無(wú)須構(gòu)造剛度矩陣，廣泛應(yīng)用于板料沖壓成形的分析計(jì)算[10-12]。本文采用動(dòng)力顯式算法，對(duì)噴管內(nèi)壁的成形過(guò)程進(jìn)行模擬。

2.1 有限元模型的建立

采用三維繪圖軟件UG建立陰模、陽(yáng)模及錐筒坯料模型，再將建好的模型以IGES格式導(dǎo)入PAM-STAMP中，對(duì)其進(jìn)行網(wǎng)格劃分，網(wǎng)格劃分采用自適應(yīng)網(wǎng)格優(yōu)化法。實(shí)際生產(chǎn)中，錐筒是由扇形板料對(duì)縫焊接而成，型面不是理想的錐形。為方便計(jì)算，幾何模型按照扇形展開(kāi)計(jì)算得出的理想錐體型面建模，并且忽略焊縫與基體的材料性能差異。陰、陽(yáng)模型面尺寸按照模具設(shè)計(jì)圖紙給出，陰、陽(yáng)間隙為2.6 mm。計(jì)算結(jié)果顯示，簡(jiǎn)化模型對(duì)計(jì)算結(jié)果的影響很小。計(jì)算過(guò)程中，陰、陽(yáng)模均為剛體，即不可變形體，坯料為變形體。圖4為有限元三維模型。

圖4 三維有限元模型Fig.4 Three-dimensional finite element model

2.2 成形過(guò)程分析

選取零件成形過(guò)程中幾個(gè)典型環(huán)節(jié)的剖面線圖、壁厚及應(yīng)力分布云圖進(jìn)行觀察分析。圖5是陽(yáng)模開(kāi)始接觸坯料時(shí)的狀態(tài)，可以看出，零件首先接觸的是陽(yáng)模中間凸起部位。此時(shí)零件小端還沒(méi)有與陰模接觸，大端與陰?？诓拷佑|。

圖5 初始狀態(tài)剖面圖Fig.5 Initial state profile

圖6依次顯示了零件開(kāi)始變形時(shí)的狀態(tài)、壁厚及應(yīng)力分布，此時(shí)零件被陽(yáng)模帶動(dòng)下行，由于中間接觸范圍的不斷擴(kuò)大，材料在該部位受到沿周向的拉伸力，由于此處受到的拉應(yīng)力最大，變形量最大，所以材料減薄也最嚴(yán)重。同時(shí)，由于零件大端在陰?？诓渴艿揭蚪佑|產(chǎn)生的阻力，與陽(yáng)模對(duì)零件的作用力相反，因此使得零件在兩個(gè)接觸區(qū)域之間的部位受到拉伸力作用。

圖6 零件大端進(jìn)入陰模時(shí)的狀態(tài)、壁厚及應(yīng)力分布Fig.6 The status,thickness and stress distribution when the large end of the part enters the female die

隨著陽(yáng)模下行，零件大端直徑大于陰模直徑的部分不斷進(jìn)入陰模并被陰?？诓渴湛s，隨著坯料不斷進(jìn)入陰模型腔，陰模大端口部圓角處材料在板厚方向的壓應(yīng)力不斷增加，當(dāng)壓應(yīng)力達(dá)到一定程度時(shí)，材料開(kāi)始失穩(wěn)起皺，如圖7所示。零件與陽(yáng)模接觸區(qū)域依然是壁厚最薄的區(qū)域，此處坯料發(fā)生脹形。大端起皺部位應(yīng)力最大，壁厚增加非常顯著。

圖8為零件小端開(kāi)始接觸陰模時(shí)的狀態(tài)、壁厚及應(yīng)力分布圖。此時(shí)零件大端已基本全部進(jìn)入陰模，大端起皺部位受到約束校正。此時(shí)，零件小端與陰模小口部位開(kāi)始接觸，坯料受到的摩擦阻力與接觸陽(yáng)模的中間部位受到的作用力相反，因此可以看到，兩處接觸部位的中間區(qū)域產(chǎn)生了環(huán)狀凸起，這是因?yàn)樵搮^(qū)域坯料受到兩側(cè)的壓應(yīng)力而產(chǎn)生的失穩(wěn)現(xiàn)象。零件小端與陰模圓角接觸部位料厚增加顯著。

圖9為成形結(jié)束時(shí)零件狀態(tài)、壁厚及應(yīng)力分布云圖。可以看出，大端的起皺和小端下方的環(huán)狀鼓包依然沒(méi)有消除。

2.3 成形結(jié)果分析

從圖9可以看出，成形結(jié)束時(shí)，陰陽(yáng)模間隙不均勻，說(shuō)明模具沒(méi)有完全閉合。零件的壁厚分布云圖顯示，零件小端圓角部位最大料厚為2.87 mm，已經(jīng)超過(guò)了陰陽(yáng)模間隙2.6 mm，這使得模具無(wú)法完全閉合，從而導(dǎo)致零件成形不到位。同時(shí)零件表面起皺也是導(dǎo)致型面尺寸超差的直接因素。模擬結(jié)果與實(shí)際生產(chǎn)中出現(xiàn)的狀況和缺陷完全一致。

圖7 零件大端失穩(wěn)起皺時(shí)的狀態(tài)、壁厚及應(yīng)力分布Fig.7 The status,thickness and stress distribution when the large end of the part is wrinkling

圖8 零件小端開(kāi)始接觸陰模時(shí)的狀態(tài)、壁厚及應(yīng)力分布Fig.8 The status,thickness and stress distribution when the small end of the part began to contact the female die

圖9 零件成形結(jié)束時(shí)的狀態(tài)、壁厚及應(yīng)力分布Fig.9 The status,thickness and stress distribution when the forming of the part is finished

3 原因分析及措施改進(jìn)

3.1 原因分析

起皺是板料在成形過(guò)程中受壓失穩(wěn)的一個(gè)主要表現(xiàn)形式。板料在塑性變形過(guò)程中會(huì)受到復(fù)雜的應(yīng)力狀態(tài)的作用，由于板厚度方向尺寸與其他兩個(gè)方向的尺寸相比很小，因此，板料在厚度方向上最不穩(wěn)定[13-15]。當(dāng)板面內(nèi)的壓應(yīng)力達(dá)到一定程度時(shí)，板厚方向最容易因受壓而不能維持穩(wěn)定的塑性變形，產(chǎn)生失穩(wěn)起皺現(xiàn)象。通過(guò)分析模擬過(guò)程及結(jié)果得知，零件大端起皺主要是受到周向非同軸平衡力，在力偶產(chǎn)生的切向壓應(yīng)力作用下產(chǎn)生的起皺。而零件小端的環(huán)形起皺主要是因?yàn)槭艿窖啬妇€方向的壓應(yīng)力而產(chǎn)生的。要想成形出型面尺寸合適的零件，必須有效消除起皺現(xiàn)象[16]。

3.2 措施改進(jìn)

從以上分析可以看出，由于零件成形過(guò)程中零件大端及小端喉部區(qū)域料厚增加顯著，使得模具無(wú)法完全閉合，零件成形不到位。因此，可以通過(guò)增大陰、陽(yáng)模局部間隙的方法改進(jìn)模具型面設(shè)計(jì)。模具型面改進(jìn)主要有以下兩個(gè)方面：①將陰、陽(yáng)模間隙整體增大至2.7mm；②由于零件大、小端受到與陰模接觸產(chǎn)生的摩擦阻力更大，因此通過(guò)曲線擬合，將大、小端模具間隙放大至2.9mm。

3.3 改進(jìn)后模擬結(jié)果分析

將改進(jìn)后的模具通過(guò)建立有限元模型進(jìn)行仿真，其結(jié)果如圖10所示。可以看出零件大端沒(méi)有起皺，小端凸起鼓包消除。這是因?yàn)槟＞唛g隙增加，使得零件大、小端受到的由陰模產(chǎn)生的摩擦阻力減小，從而使得零件在大端口部和小端喉部位置的料厚增加減小，反過(guò)來(lái)又促進(jìn)了零件更加順利地進(jìn)入陰模型腔。由于成形后壁厚小于陰、陽(yáng)模間隙，模具能夠完全閉合，從根本上消除了零件的起皺缺陷。

圖10 改進(jìn)后零件成形結(jié)束時(shí)的狀態(tài)、壁厚及應(yīng)力分布Fig.10 The status,thickness and stress distribution when the forming of the part is finished after improved

4 試驗(yàn)驗(yàn)證

根據(jù)仿真計(jì)算結(jié)果改進(jìn)陰、陽(yáng)模型面參數(shù)和間隙。成形后零件最小壁厚2.33 mm，符合設(shè)計(jì)圖紙中最小壁厚2.2 mm的要求；型面與樣板間隙最大為0.25 mm，滿足≤0.3 mm的設(shè)計(jì)要求，且沒(méi)有出現(xiàn)起皺及鼓包現(xiàn)象，如圖11所示。

圖11 實(shí)際零件與仿真零件外觀對(duì)比圖Fig.11 Contrast drawing of actual and simulation parts

對(duì)結(jié)果進(jìn)一步分析對(duì)比，以小端為起始，在仿真結(jié)果零件和實(shí)際零件上，沿母線方向分別均勻選取了14個(gè)位置點(diǎn)對(duì)壁厚進(jìn)行測(cè)量，測(cè)量結(jié)果如表1所示。

表1 仿真零件壁厚與實(shí)際零件壁厚對(duì)比

Tab.1 Comparison of wall thickness between simulation model and actual partmm

位置仿真零件壁厚實(shí)際零件壁厚52.562.52252.672.68452.612.62652.542.50852.482.431052.412.381252.352.331452.422.441652.462.451852.562.522052.602.632252.622.602452.652.612652.682.70

從表1可以看出：仿真零件壁厚與實(shí)際零件壁厚在相對(duì)應(yīng)位置上的最大差值僅為0.05 mm，且壁厚沿母線方向的變化趨勢(shì)完全一致。提取仿真結(jié)果零件型面曲線與理論曲線進(jìn)行比對(duì)，測(cè)量出了這14個(gè)位置點(diǎn)的零件型面間隙，與實(shí)際零件型面與樣板間隙進(jìn)行比較(表2)，可以看出：仿真零件與實(shí)際零件型面間隙最大差值僅為0.06 mm，且型面間隙沿母線方向的變化趨勢(shì)一致。可見(jiàn)仿真計(jì)算結(jié)果是準(zhǔn)確有效的。

表2 仿真零件型面間隙與實(shí)際零件型面間隙對(duì)比

Tab.2 Comparison of surface gap between simulation model and actual partmm

位置仿真零件型面間隙實(shí)際零件型面間隙50.02<0.05250.120.15450.130.10650.120.15850.150.151050.160.201250.170.201450.150.151650.190.251850.210.252050.180.22250.120.152450.100.102650.05<0.05

5 結(jié)語(yǔ)

通過(guò)有限元仿真軟件，對(duì)某型號(hào)噴管內(nèi)壁的成形過(guò)程進(jìn)行模擬，通過(guò)對(duì)成形過(guò)程及結(jié)果觀察，結(jié)合理論分析，找出了該零件成形后產(chǎn)生起皺缺陷的原因，并通過(guò)改進(jìn)模具參數(shù)設(shè)計(jì)，生產(chǎn)出了符合設(shè)計(jì)圖紙要求的零件，對(duì)今后此類(lèi)零件工藝方案的制定和實(shí)施提供了有益的參考。