趙 元?jiǎng)I(yè)峰
(1.遼寧省數(shù)控機(jī)床信息物理融合與智能制造重點(diǎn)實(shí)驗(yàn)室, 遼寧 撫順 113122;2.沈陽(yáng)工學(xué)院, 遼寧 撫順 113122)
隨著現(xiàn)代工業(yè)分工與協(xié)作專業(yè)化的不斷推進(jìn),市場(chǎng)配給對(duì)傳統(tǒng)工業(yè)企業(yè)的作用越來(lái)越顯著。尤其對(duì)于大型工業(yè)產(chǎn)品來(lái)說(shuō),一家企業(yè)包辦所有零件并總裝的情況越來(lái)越少。為了集中力量發(fā)展自身核心技術(shù)專長(zhǎng),降低生產(chǎn)成本,大部分企業(yè)采取從專業(yè)化供應(yīng)商采購(gòu)成品,自己生產(chǎn)部分專用零件并總裝的模式進(jìn)行生產(chǎn),比如波音、空客等企業(yè)自身僅進(jìn)行設(shè)計(jì)、裝配和試飛等工作,將非重要零件的加工外包至全球多家供應(yīng)商,單個(gè)供應(yīng)商只生產(chǎn)獨(dú)立的零件,無(wú)法接觸總體設(shè)計(jì)思想。通用、寶馬等企業(yè)自己生產(chǎn)核心部件(發(fā)動(dòng)機(jī)、變速箱等),把總裝線設(shè)在勞動(dòng)力成本低的國(guó)家。因此這類企業(yè)生產(chǎn)過(guò)程中零部件周轉(zhuǎn)運(yùn)輸是相當(dāng)頻繁的。
由于航空產(chǎn)品一般外廓尺寸較大、結(jié)構(gòu)復(fù)雜精密,并且部分產(chǎn)品對(duì)運(yùn)輸環(huán)境要求較高,如控制環(huán)境濕度、隔絕氧氣防止腐蝕等特殊要求,因此航空產(chǎn)品從供應(yīng)商安全可靠地運(yùn)輸至總裝企業(yè)是一個(gè)急需重視的挑戰(zhàn)。
目前在航空產(chǎn)品方面,較為常用的包裝運(yùn)輸設(shè)備主要有木質(zhì)包裝箱、金屬材料包裝箱等[1],以木材料包裝箱應(yīng)用量居多,內(nèi)部輔以塑料包裝袋抽真空防腐。抽真空防腐在隔離性能上存在缺陷,輕微的漏氣即可使產(chǎn)品與氧氣直接接觸,存在腐蝕風(fēng)險(xiǎn);如采用充入高壓氮?dú)獾姆绞礁綦x,只要氮?dú)鈮毫Ω哂诖髿鈮?,即使存在輕微的漏氣,在內(nèi)壓失效前也不會(huì)讓氧氣混入,隔離效果佳。但此種內(nèi)壓包裝箱對(duì)箱體強(qiáng)度要求較高,需采用金屬材料設(shè)計(jì)箱體,并保證箱體在充壓后強(qiáng)度及穩(wěn)定性方面滿足運(yùn)輸要求。相關(guān)學(xué)者在參考?jí)毫θ萜飨嚓P(guān)標(biāo)準(zhǔn)對(duì)罐式集裝箱進(jìn)行強(qiáng)度分析[2]、鋼制壓力容器可靠性設(shè)計(jì)[3]和包裝箱結(jié)構(gòu)過(guò)渡區(qū)的強(qiáng)度優(yōu)化[4]等方面做出了有益的探索和研究,但研究對(duì)象的結(jié)構(gòu)大都較為規(guī)則,基本為圓形、橢圓形、方形等。對(duì)異型結(jié)構(gòu)的金屬包裝箱,由于其在內(nèi)壓載荷下內(nèi)部應(yīng)力分布更為復(fù)雜,對(duì)其進(jìn)行深入的研究具有重要的理論和實(shí)踐意義。
應(yīng)用國(guó)際通用的有限元分析軟件ANSYS對(duì)包裝箱進(jìn)行強(qiáng)度計(jì)算和應(yīng)力分布研究,結(jié)合包裝箱具體結(jié)構(gòu)特點(diǎn)采用三維實(shí)體有限元算法。有限元算法的求解步驟為:
(1)物理問(wèn)題及求解區(qū)域定義:依據(jù)工程實(shí)際問(wèn)題,在誤差可接受的范圍內(nèi)給出求解區(qū)域的物理特性和幾何外廓。
(2)對(duì)求解區(qū)域采取離散化處理手段:求解區(qū)域劃分成大小、形狀各異,且彼此節(jié)點(diǎn)連接的求解單元所組成的離散區(qū)域,此過(guò)程即為有限元網(wǎng)絡(luò)劃分。一般來(lái)說(shuō)單元體積越小,離散區(qū)域的仿真程度越好,計(jì)算結(jié)果就越精確,但同時(shí)其計(jì)算量和誤差可能增大,因此求解區(qū)域的離散化需考慮實(shí)際問(wèn)題的需求和可接受的誤差范圍。
(3)確定變量控制參數(shù):從工程實(shí)際問(wèn)題抽象出物理問(wèn)題,并用包含問(wèn)題變量控制參數(shù)的微分方程組表示出來(lái)。為適用于有限元計(jì)算方法,一般把微分方程組簡(jiǎn)化為與之等價(jià)的泛函數(shù)形式。
(4)單元矩陣推導(dǎo):首先創(chuàng)建一個(gè)近似的模擬解,推算出有限元列式,包括確定正確的單元坐標(biāo)系,創(chuàng)建單元函數(shù),并推導(dǎo)出各個(gè)狀態(tài)變量的離散關(guān)系,組成單元矩陣(一般稱其為剛度矩陣)。為了有限元求解過(guò)程收斂,單元矩陣需遵守多項(xiàng)要求。對(duì)工程實(shí)際問(wèn)題來(lái)說(shuō),應(yīng)關(guān)注每一單元的求解性能與邊界,單元形狀應(yīng)當(dāng)盡量簡(jiǎn)化,畸形單元不僅誤差大,而且有可能缺秩,從而使問(wèn)題求解過(guò)程無(wú)法收斂。
(5)矩陣裝配求解:?jiǎn)卧仃囇b配得到總體剛度矩陣,對(duì)總體剛度矩陣方程組求解。
包裝箱的實(shí)體模型如圖1所示,坐標(biāo)系確定為:包裝箱長(zhǎng)度方向定義為X軸,水平方向定義為Y軸,豎直方向定義為Z軸。
圖1 包裝箱實(shí)體模型圖
綜合考慮使用條件、強(qiáng)度、剛度、焊接性能、成本等相關(guān)因素,包裝箱材料選用Q460號(hào)鋼,具體參數(shù)見表1。
表1 材料數(shù)據(jù)
考慮包裝箱使用過(guò)程,選取如下兩個(gè)狀態(tài):
吊起狀態(tài):將包裝箱吊起,脫離地面,受包裝箱及產(chǎn)品重量作用。
氮封過(guò)載狀態(tài):包裝箱穩(wěn)定可靠地放置于地面,結(jié)構(gòu)受氣體壓力、包裝箱自重及產(chǎn)品重力作用,同時(shí)還需承受過(guò)載時(shí)的慣性力。
包裝箱重力載荷:包裝箱質(zhì)量為2 100 kg;
重力載荷:產(chǎn)品及其托架質(zhì)量為2 400 kg;
氣體壓力載荷:內(nèi)部氣體壓力為30 000 Pa(表壓);
過(guò)載狀態(tài)的慣性力:包裝箱需在九級(jí)海況下存放,因此需要對(duì)九級(jí)海況產(chǎn)生的慣性力進(jìn)行研究。通過(guò)調(diào)研分析,參考中國(guó)船級(jí)社(ccs)2007版《船舶與海上設(shè)施起重設(shè)備規(guī)范》中規(guī)定,物體在放置狀態(tài)下,結(jié)構(gòu)設(shè)計(jì)時(shí)應(yīng)考慮能承受下列兩種力綜合作用:
(1)垂直于甲板的加速度為±1.0 g;
前后方向平行于甲板的加速度為±0.5 g;
靜橫傾30°;
風(fēng)速55 m/s,作用于前后方向。
(2)垂直于甲板的加速度為±1.0 g;
橫向平行于甲板的加速度為±0.5 g;
靜橫傾30°;
風(fēng)速55 m/s,作用于橫向。
參考中國(guó)船級(jí)社《集裝箱檢驗(yàn)規(guī)范》2008版,集裝箱罐體在設(shè)計(jì)和制造時(shí)應(yīng)遵循壓力容器規(guī)則,并充分考慮在其規(guī)定的最大載荷下,承受運(yùn)輸過(guò)程中產(chǎn)生的如下慣性載荷:
1)在運(yùn)動(dòng)方向:總質(zhì)量的兩倍;
2)在與運(yùn)動(dòng)方向成直角的水平方向:總質(zhì)量(當(dāng)不能明確運(yùn)動(dòng)方向時(shí),可接受的最大負(fù)載應(yīng)為重量的兩倍);
3)垂直向上:總質(zhì)量;
4)垂直向下:總質(zhì)量的兩倍(包括重力)。
由于罐式集裝箱的結(jié)構(gòu)形式和工況與包裝箱更為相似,且其標(biāo)準(zhǔn)更為嚴(yán)格,因此選擇中國(guó)船級(jí)社《集裝箱檢驗(yàn)規(guī)范》2008版中關(guān)于產(chǎn)品在儲(chǔ)運(yùn)過(guò)程中所產(chǎn)生的慣性力作為考核載荷。
除包裝箱承受水平方向慣性力時(shí)載荷不具有對(duì)稱性,采用整體有限元模型進(jìn)行計(jì)算,其余狀態(tài)的結(jié)構(gòu)和載荷均具有對(duì)稱性,計(jì)算時(shí)應(yīng)用1/2有限元模型,以提升仿真效率。應(yīng)用solid186單元進(jìn)行網(wǎng)格劃分,solid186單元為高階實(shí)體帶中間節(jié)點(diǎn)的20節(jié)點(diǎn)實(shí)體單元,適用于雙自由度位移,可用于外廓形狀不規(guī)則產(chǎn)品的模擬。此單元的定義借助于20個(gè)節(jié)點(diǎn),其中每個(gè)節(jié)點(diǎn)具有沿xyz方向平移的3個(gè)自由度。SOLID186單元可模擬空間各向異性特征,支持彈塑性、超彈、蠕變、應(yīng)力強(qiáng)化、大柔度分析能力。同時(shí)其可采用混合模式,模擬不可壓縮的彈塑性模型和不可壓縮超彈材料,具有各種輸出選項(xiàng)。劃分后的1/2有限元模型,共有339 348個(gè)單元、1 089 106個(gè)節(jié)點(diǎn),如圖2所示。
約束根據(jù)不同計(jì)算狀態(tài)選取,吊起狀態(tài)約束吊耳處,氮封過(guò)載狀態(tài)施加地面約束,采用1/2有限元模型時(shí)施加對(duì)稱邊界條件。
首先對(duì)實(shí)際問(wèn)題模型進(jìn)行優(yōu)化,由于其具有結(jié)構(gòu)、負(fù)荷左右對(duì)稱的特點(diǎn),可在對(duì)稱面上施加法向?qū)ΨQ約束,水平過(guò)載狀態(tài)由于載荷無(wú)對(duì)稱型,只能采用全模型進(jìn)行仿真。然后根據(jù)吊點(diǎn)處系留方式約束吊點(diǎn)位置自由度。最后在包裝箱底部產(chǎn)品及其托架的托點(diǎn)處施加重力載荷。
考慮到有限元模型節(jié)點(diǎn)數(shù)量較多,求解規(guī)模大,為有效節(jié)省機(jī)時(shí),應(yīng)用條件共軛求解法(PCG)進(jìn)行求解。首先創(chuàng)建網(wǎng)格單元的邊界條件矩陣,自動(dòng)設(shè)置問(wèn)題初始近似解,然后通過(guò)迭代使求解結(jié)果收斂,滿足給定的容差范圍。根據(jù)結(jié)構(gòu)材料和有限元單元的差異,條件共軛求解法會(huì)自動(dòng)設(shè)置適用的初始求解條件,然后再進(jìn)行求解。條件共軛法的優(yōu)點(diǎn)在于節(jié)省機(jī)時(shí),并且求解規(guī)模越大,節(jié)省的時(shí)間越多,應(yīng)用效果良好。但需注意,病態(tài)矩陣將使求解結(jié)果不收斂,此時(shí)一方面需改善網(wǎng)格爭(zhēng)來(lái)那個(gè),重新劃分網(wǎng)格分析,另一方面考慮應(yīng)用其他求解器求解。
仿真實(shí)驗(yàn)的結(jié)果評(píng)定目前主要有四大強(qiáng)度準(zhǔn)則理論:
1)第一強(qiáng)度理論(最大拉應(yīng)力準(zhǔn)則)。
認(rèn)為造成材料脆性斷裂破壞的主導(dǎo)因素是單向最大拉應(yīng)力,與應(yīng)力狀態(tài)無(wú)關(guān)。當(dāng)結(jié)構(gòu)拉應(yīng)力達(dá)到單方向極限應(yīng)力,材料就會(huì)產(chǎn)生脆性斷裂。適用于脆性材料,例如:鑄鐵。
2)第二強(qiáng)度理論(最大線應(yīng)變準(zhǔn)則)。
認(rèn)為最大線應(yīng)變是引起斷裂的主導(dǎo)因素,與應(yīng)力狀態(tài)無(wú)關(guān)。當(dāng)最大線應(yīng)變超過(guò)單方向極限值,材料就將產(chǎn)生脆性斷裂破壞。只適用極少數(shù)脆性材料,應(yīng)用很少。
3)第三強(qiáng)度理論(最大切應(yīng)力準(zhǔn)則)。
認(rèn)為最大切應(yīng)力是產(chǎn)生屈服現(xiàn)象的主導(dǎo)因素,與應(yīng)力狀態(tài)無(wú)關(guān)。當(dāng)最大切應(yīng)力τmax大于等于單向應(yīng)力狀態(tài)下的極限切應(yīng)力τ0時(shí),材料即發(fā)生屈服破壞。主要適用于塑性材料,形式較為簡(jiǎn)單,應(yīng)用范圍廣泛。
4)第四強(qiáng)度理論(結(jié)構(gòu)改變比能準(zhǔn)則)。
認(rèn)為結(jié)構(gòu)改變比能是引起材料屈服破壞的主導(dǎo)因素,與應(yīng)力狀態(tài)無(wú)關(guān)。當(dāng)結(jié)構(gòu)內(nèi)某一點(diǎn)的改變比能超過(guò)單方向材料極限值,材料就要發(fā)生屈服破壞。適用于大多數(shù)塑性材料,計(jì)算結(jié)果比第三強(qiáng)度理論準(zhǔn)確,但求解過(guò)程不如第三強(qiáng)度理論方便。
根據(jù)實(shí)際問(wèn)題特點(diǎn),本研究依據(jù)第四強(qiáng)度理論(形狀改變比能理論)進(jìn)行結(jié)果評(píng)定,對(duì)結(jié)構(gòu)的范式等效應(yīng)力(Von Mises Stress)進(jìn)行考核,用應(yīng)力等值線來(lái)描述模型內(nèi)部的應(yīng)力分布狀態(tài),從而快速確定模型中的最危險(xiǎn)區(qū)域。
為考察包裝箱的承重能力,根據(jù)吊點(diǎn)位置分兩種情況進(jìn)行了仿真試驗(yàn)。
3.1.1 吊點(diǎn)在箱蓋中部
結(jié)構(gòu)變形集中于包裝箱底部中間位置,變形最大值為0.405 mm??紤]包裝箱整體結(jié)構(gòu)和尺寸,此變形在可接受范圍內(nèi)。結(jié)構(gòu)高等效應(yīng)力區(qū)為吊點(diǎn)周邊位置,等效應(yīng)力最大值為40.8 MPa,但遠(yuǎn)小于材料的屈服極限σ0.2(460 MPa)。此狀態(tài)結(jié)構(gòu)整體應(yīng)力水平均在較低范圍,變形值也較小,仿真試驗(yàn)結(jié)果云圖見圖3~4。
圖3 包裝箱位移分布云圖 圖4 包裝箱Mises等效應(yīng)力仿真結(jié)果
3.1.2 吊點(diǎn)在箱底上部
結(jié)構(gòu)變形在包裝箱底部中間位置較為突出,最大值為0.462 mm。考慮包裝箱整體結(jié)構(gòu)和尺寸,此變形在可接受范圍內(nèi)。包裝箱Mises等效應(yīng)力在吊點(diǎn)附近上升明顯,等效應(yīng)力最大值為83.9 MPa,但仍小于材料屈服極限σ0.2(460 MPa)。此狀態(tài)結(jié)構(gòu)整體應(yīng)力水平均在較低范圍,變形值也較小,仿真試驗(yàn)結(jié)果云圖見圖5~6。
圖5 包裝箱位移分布云圖 圖6 包裝箱Mises等效應(yīng)力仿真結(jié)果
根據(jù)《集裝箱檢驗(yàn)規(guī)范》的要求,對(duì)包裝箱在最大工作負(fù)荷下承受本文2.4節(jié)給定的各種慣性力單獨(dú)作用的情況進(jìn)行了有限元計(jì)算,根據(jù)過(guò)載方向的不同分為4種仿真狀態(tài)。
3.2.1 運(yùn)動(dòng)方向過(guò)載
結(jié)構(gòu)變形極大值出現(xiàn)在包裝箱側(cè)面中部位置,最大變形值為3.813 mm。考慮包裝箱整體結(jié)構(gòu)和尺寸,此變形在可接受范圍內(nèi),也可通過(guò)增大變形區(qū)承力梁規(guī)格或減小梁間距的方式減小變形量。
包裝箱Mises等效應(yīng)力大應(yīng)力區(qū)在吊點(diǎn)附件,最大等效應(yīng)力值為233 MPa,主體部位最大等效應(yīng)力為181 MPa,均小于材料屈服極限σ0.2(460 MPa)。此狀態(tài)結(jié)構(gòu)應(yīng)力主要由承力梁承擔(dān),應(yīng)力水平相對(duì)材料性能較低,符合設(shè)計(jì)思想,仿真試驗(yàn)結(jié)果云圖見圖7~8。
圖7 包裝箱位移分布云圖 圖8 包裝箱Mises等效應(yīng)力仿真結(jié)果
3.2.2 水平方向過(guò)載
結(jié)構(gòu)變形極大值出現(xiàn)在包裝箱側(cè)面中部位置,變形最大值為8.823 mm??紤]包裝箱整體結(jié)構(gòu)和尺寸,此變形在可接受范圍內(nèi),也可通過(guò)增大變形區(qū)承力梁規(guī)格或減小梁間距的方式減小變形量。
包裝箱Mises等效應(yīng)力大應(yīng)力區(qū)為吊點(diǎn)周邊位置,等效應(yīng)力最大值為297 MPa,主體部位最大等效應(yīng)力為231 MPa,均小于材料屈服極限σ0.2(460 MPa)。此狀態(tài)結(jié)構(gòu)應(yīng)力主要由承力梁承擔(dān),應(yīng)力水平相對(duì)材料性能較低,符合設(shè)計(jì)思想,仿真試驗(yàn)結(jié)果云圖見圖9~10。
圖9 包裝箱位移分布云圖 圖10 包裝箱Mises等效應(yīng)力仿真結(jié)果
3.2.3 垂直向上過(guò)載
結(jié)構(gòu)變形在包裝箱側(cè)面中部位置比較明顯,最大變形量為3.393 mm??紤]包裝箱整體結(jié)構(gòu)和尺寸,此變形在可接受范圍內(nèi),也可通過(guò)增大變形區(qū)承力梁規(guī)格或減小梁間距的方式減小變形量。
包裝箱Mises等效應(yīng)力在吊點(diǎn)周邊達(dá)到最大,等效應(yīng)力值最大可達(dá)225 MPa,主體部位最大等效應(yīng)力為175 MPa,均小于材料屈服極限σ0.2(460 MPa)。此狀態(tài)結(jié)構(gòu)應(yīng)力主要由承力梁承擔(dān),應(yīng)力水平相對(duì)材料性能較低,符合設(shè)計(jì)思想,仿真試驗(yàn)結(jié)果云圖見圖11~12。
圖11 包裝箱位移分布云圖 圖12 包裝箱Mises等效應(yīng)力仿真結(jié)果
3.2.4 垂直向下過(guò)載
結(jié)構(gòu)變形在包裝箱側(cè)面中部位置較為明顯,最大變形值為3.526 mm。考慮包裝箱整體結(jié)構(gòu)和尺寸,此變形在可接受范圍內(nèi),也可通過(guò)增大變形區(qū)承力梁規(guī)格或減小梁間距的方式減小變形量。
包裝箱Mises等效應(yīng)力較大區(qū)域?yàn)榈觞c(diǎn)周邊位置,最大等效應(yīng)力值為226 MPa,主體部位最大等效應(yīng)力為176 MPa,均小于材料屈服極限σ0.2(460 MPa)。此狀態(tài)結(jié)構(gòu)應(yīng)力主要由承力梁承擔(dān),應(yīng)力水平相對(duì)材料性能較低,符合設(shè)計(jì)思想,仿真試驗(yàn)結(jié)果云圖見圖13~14。
圖13 包裝箱位移分布云圖 圖14 包裝箱Mises等效應(yīng)力仿真結(jié)果
1)兩種吊點(diǎn)位置吊起狀態(tài)的最大位移均出現(xiàn)在箱底中部,吊點(diǎn)處為大應(yīng)力區(qū)。包裝箱變形及等效應(yīng)力最大值見表2。從表中數(shù)據(jù)可看出,吊起狀態(tài)包裝箱承受自身重力及重力載荷時(shí)的應(yīng)力及位移均較小,包裝箱的承重能力符合要求。
表2 吊起狀態(tài)計(jì)算結(jié)果表
2)氮封過(guò)載狀態(tài)水平方向過(guò)載時(shí)大應(yīng)力區(qū)和最大位移在箱體側(cè)面系留點(diǎn)處,其余方向過(guò)載時(shí)最大位移出現(xiàn)在箱體側(cè)面,大應(yīng)力區(qū)為箱底側(cè)面方鋼及與之接觸的部位。不同過(guò)載方向時(shí)的包裝箱最大位移、最大等效應(yīng)力和主體等效應(yīng)力值見表3,安全系數(shù)見表4。
表3 氮封過(guò)載狀態(tài)計(jì)算結(jié)果表
表4 氮封過(guò)載狀態(tài)安全系數(shù)表
《集裝箱檢驗(yàn)規(guī)范》中規(guī)定應(yīng)通過(guò)實(shí)驗(yàn)測(cè)試或有限元計(jì)算確認(rèn)罐式集裝箱在最大工作負(fù)荷下承受各種慣性力單獨(dú)作用下,其罐體、框架及其之間連接構(gòu)件的Mises等效應(yīng)力安全系數(shù)如下:
(1)對(duì)于屈服點(diǎn)明確的金屬,按已確定的屈服應(yīng)力,取安全系數(shù)1.5;
(2)對(duì)于屈服點(diǎn)尚未確定的金屬材料,一般以0.2%來(lái)規(guī)定屈服應(yīng)力,取安全系數(shù)1.5。
通過(guò)表4中數(shù)據(jù)可知,結(jié)構(gòu)承受不同方向過(guò)載力時(shí)的安全系數(shù)均大于1.5,滿足結(jié)構(gòu)靜強(qiáng)度設(shè)計(jì)要求。
1)因氣體自身張力作用,為保證箱體在內(nèi)部壓力載荷作用下可靠工作,包裝箱蒙皮應(yīng)盡量設(shè)計(jì)成圓弧狀。
2)在受尺寸約束限制,包裝箱整體不能設(shè)計(jì)成圓筒形時(shí),造成了圓弧過(guò)渡不平滑區(qū)應(yīng)力集中較明顯,對(duì)這些區(qū)域應(yīng)用較多的梁進(jìn)行加強(qiáng),對(duì)圓弧過(guò)渡較平滑區(qū)域(箱蓋上部和兩端堵頭)用較少的梁加強(qiáng)即可。
3)在梁的尺寸規(guī)格選取方面,考慮到箱體中部受載荷作用產(chǎn)生的變形相對(duì)較大,吊點(diǎn)位置在吊起狀態(tài)時(shí)承受大部分重力載荷,因此上述兩個(gè)部位采用尺寸規(guī)格較大的梁進(jìn)行加強(qiáng),其余部位采用尺寸規(guī)格較小的梁進(jìn)行加強(qiáng)。