廖嶸,楊帆,許為,張振麒
基于HyperWorks的反力架優(yōu)化設(shè)計(jì)
廖嶸,楊帆,許為,張振麒
(柳州泰姆預(yù)應(yīng)力機(jī)械有限公司,廣西 柳州 545000)
運(yùn)用HyperWorks的尺寸優(yōu)化功能對(duì)擠壓機(jī)受力部件的反力架進(jìn)行優(yōu)化分析。經(jīng)過靜力分析、第一次優(yōu)化、重新調(diào)整優(yōu)化參數(shù)和第二次優(yōu)化,最終滿足目標(biāo)函數(shù)和約束函數(shù)要求,得出了最優(yōu)的設(shè)計(jì)方案,最終優(yōu)化結(jié)果較初始方案材料減少38.0%,既滿足了強(qiáng)度與剛度的設(shè)計(jì)要求,又降低了成本。其間運(yùn)用Shell殼單元和Beam梁單元替代實(shí)體單元簡化模型,并使用其1/2模型進(jìn)行對(duì)稱分析,從而大大減少了網(wǎng)格數(shù)量,提高了運(yùn)算速度和效率;應(yīng)用第二類響應(yīng)功能創(chuàng)建約束函數(shù),較好地解決了直接調(diào)用第一類響應(yīng)無法處理的問題。
反力架;有限元;對(duì)稱分析;第二類響應(yīng)
隨著公路梁用工作錨板孔擠壓工藝的成熟及廣泛應(yīng)用,為了適應(yīng)可擠壓工作錨板2~22孔的需要,設(shè)計(jì)開發(fā)一套2工位且單工位最大可承受4.9×106N力的擠壓機(jī),因此在重載荷作用下整個(gè)反力架結(jié)構(gòu)的強(qiáng)度和變形問題顯得至關(guān)重要。柴振嶺等[1]運(yùn)用有限元數(shù)值分析與理論分析相結(jié)合的方法成功研制了拉壓雙作用梁式反力架,出于安全性和穩(wěn)定性的考慮對(duì)反力架的強(qiáng)度和剛度進(jìn)行了校驗(yàn),未對(duì)反力架進(jìn)行優(yōu)化分析;王凱等[2]提出了在有限元分析的基礎(chǔ)上,應(yīng)用粒子群優(yōu)化算法(Particle swarm optimization,PSO)求解反力架的優(yōu)化問題,其實(shí)質(zhì)是聯(lián)合了有限元分析軟件ABAQUS和多學(xué)科優(yōu)化設(shè)計(jì)軟件ISIGHT開展反力架的結(jié)構(gòu)優(yōu)化設(shè)計(jì),通過優(yōu)化分析實(shí)現(xiàn)了反力架輕量化和降低成本的目標(biāo),但其過程較為繁瑣。
HyperWorks使用先進(jìn)的優(yōu)化算法和優(yōu)化引擎讓用戶結(jié)合拓?fù)洹⑿蚊?、尺寸和形狀?yōu)化方法創(chuàng)建更好的設(shè)計(jì)方案[3],其中尺寸優(yōu)化被廣泛應(yīng)用。張積洪等[4]對(duì)車架結(jié)構(gòu)進(jìn)行拓?fù)?、形貌?yōu)化和尺寸優(yōu)化;張時(shí)鋒等[5]運(yùn)用尺寸優(yōu)化對(duì)擋板門結(jié)構(gòu)進(jìn)行輕量化設(shè)計(jì);徐晶等[6]運(yùn)用尺寸優(yōu)化對(duì)液壓機(jī)機(jī)身進(jìn)行結(jié)構(gòu)優(yōu)化。結(jié)合擠壓機(jī)的反力架實(shí)際情況,先用HyperWorks對(duì)反力架進(jìn)行靜態(tài)力學(xué)分析,有了一個(gè)初始的分析結(jié)果,根據(jù)結(jié)果再設(shè)定設(shè)計(jì)變量、目標(biāo)函數(shù)和約束函數(shù),最后進(jìn)行尺寸優(yōu)化設(shè)計(jì)。
在建立約束函數(shù)和目標(biāo)函數(shù)時(shí),可直接調(diào)用內(nèi)部的第一類響應(yīng)計(jì)算結(jié)果,如位移、應(yīng)力、作用力、接觸應(yīng)力、頻率、柔度、體積、質(zhì)量等,這些響應(yīng)都可以直接調(diào)用并可轉(zhuǎn)化為約束函數(shù)或目標(biāo)函數(shù)。當(dāng)遇到第一類響應(yīng)解決不了的問題時(shí)就可以應(yīng)用第二類響應(yīng)[4],即調(diào)用已有第一類響應(yīng)采用內(nèi)部自帶常用公式進(jìn)行算術(shù)運(yùn)算,如SUM求和、AVG平均值、SSQ平方和、MAX最大值、MIN最小值、SUMABS絕對(duì)值求和等,也可以進(jìn)行自定義函數(shù)。在對(duì)反力架的分析中將會(huì)用到以自定義函數(shù)建立第二類響應(yīng)并創(chuàng)建約束函數(shù)進(jìn)行優(yōu)化分析。
響應(yīng)對(duì)設(shè)計(jì)變量變化的靈敏度數(shù)值直觀反映設(shè)計(jì)變量對(duì)性能的影響程度和趨勢(shì),因此可以用于快速篩選出重要的變量,并對(duì)其值進(jìn)行嘗試性地和有方向性地調(diào)整,經(jīng)過不斷的迭代運(yùn)算,使其結(jié)果逐漸滿足約束函數(shù)和目標(biāo)函數(shù),這樣可以較快速度自動(dòng)的計(jì)算出最優(yōu)方案,以更少的材料成本滿足強(qiáng)度與剛度的設(shè)計(jì)要求。
擠壓機(jī)反力架主要由上箱梁、下箱梁、立柱、螺母壓板、壓墊板等組成:上、下箱梁采用鋼板拼焊而成,外形尺寸1750 mm×620 mm×400 mm;連接上、下箱梁的立柱采用36根32 mm精軋螺紋鋼,是具有高強(qiáng)度、高尺寸精度的螺紋鋼筋[7],該鋼筋在任意截面處都可擰上帶有內(nèi)螺紋的連接螺母進(jìn)行連接,上、下箱梁螺紋鋼連接兩端的長度為1700 mm;螺母壓板和壓墊板分別與上、下箱焊接成一體。反力架的三維數(shù)模如圖1所示。
圖1 反力架三維數(shù)模
上、下箱梁都是薄板件拼焊而成,因此采用Shell殼單元[8]替換薄板件實(shí)體單元,螺紋鋼采用Beam梁單元[9],截面半徑設(shè)置為16 mm。實(shí)體單元網(wǎng)格數(shù)量太大,會(huì)造成分析計(jì)算耗時(shí)過長,甚至?xí)驗(yàn)殡娔X主機(jī)性能較差而無法計(jì)算出結(jié)果。反力架模型呈對(duì)稱分布,為了減少計(jì)算量,取1/2模型進(jìn)行對(duì)稱分析[10]。要使用對(duì)稱分析必須滿足3個(gè)條件,即幾何形狀、邊界條件和載荷都具有對(duì)稱性。彭春雷等[11]通過對(duì)典型的筋整體壁板進(jìn)行優(yōu)化設(shè)計(jì)的實(shí)例說明了對(duì)稱分析的可行性,可大大節(jié)省分析時(shí)間,提高分析效率。
為便于后續(xù)優(yōu)化的進(jìn)行,將需優(yōu)化的各板件分別命名,如圖2所示,厚度根據(jù)經(jīng)驗(yàn)都設(shè)置為40 mm,工作錨板的標(biāo)準(zhǔn)厚度為80 mm,網(wǎng)格尺寸設(shè)定為10 mm,得到單元數(shù)量共89153個(gè)、節(jié)點(diǎn)數(shù)量共86294個(gè)。
圖2 各板分布圖
各部件材料屬性如表1所示。
表1 材料屬性
1.3.1 約束邊界
(1)對(duì)稱約束設(shè)置:在對(duì)稱面上建立一個(gè)局部坐標(biāo)系,將模型所有節(jié)點(diǎn)歸屬這個(gè)坐標(biāo)系,約束局部坐標(biāo)系面上所有節(jié)點(diǎn)在方向的自由度,如圖3所示;
圖3 對(duì)稱面上的局部坐標(biāo)系圖
(2)下壓墊板上環(huán)面(570 mm×210 mm)的所有節(jié)點(diǎn)與中心點(diǎn)耦合以剛性單元RBE2連接,約束中心節(jié)點(diǎn)除軸方向外其余五個(gè)方向的自由度,如圖4所示;
圖4 下壓墊板邊界約束圖
(3)上螺母壓板與上頂板、下螺母壓板與下頂板、上壓墊板與上底板、下壓墊板與下底板、工作錨板與上頂板分別采用Tie綁定連接;
(4)螺紋鋼Beam梁單元兩端的節(jié)點(diǎn)分別與上、下螺母壓板螺母墊圈面上的節(jié)點(diǎn)用剛性單元RBE3進(jìn)行耦合連接代替螺母緊固[12],如圖5所示。
1.3.2載荷施加
22孔工作錨板直徑226 mm,采用剛性單元RBE3將工作錨板面上所有節(jié)點(diǎn)與中心點(diǎn)耦合連接,選擇最大受力工況,在耦合中心節(jié)點(diǎn)處軸方向施加4.9×106N的力,如圖6所示。
圖5 RBE3單元耦合圖
HyperWorks中RBE2和RBE3都是剛性桿單元,RBE2是被耦合的所有從節(jié)點(diǎn)與一個(gè)主節(jié)點(diǎn)位移始終保持一致,RBE3則是被耦合的各個(gè)主節(jié)點(diǎn)都可獨(dú)立位移且與單個(gè)從節(jié)點(diǎn)構(gòu)成連接關(guān)系[13],載荷施加在從節(jié)點(diǎn)上,以加權(quán)方式分布在各主節(jié)點(diǎn)之上。
圖6 載荷施加
通過對(duì)反力架進(jìn)行常規(guī)的線性靜態(tài)分析,獲得了初始的靜力分析結(jié)果如圖7~9所示:軸方向最大位移3.241 mm,各板最大Mises應(yīng)力149.5 MPa,螺紋鋼最大應(yīng)力376.2 MPa。
圖7 位移云圖(1/2模型)
圖8 各板應(yīng)力云圖
圖9 螺紋鋼應(yīng)力云圖
從圖7可知:軸變形量較大,整個(gè)上箱梁整體向上位移較大的原因是螺紋鋼橫截面總面積較小的緣故,根據(jù)胡克定律有:
式中:為截面應(yīng)力;為彈性模量;為應(yīng)變。
又有:
式中:為螺紋鋼受到的拉力;為螺紋鋼的有效面積。
式中:Δ為螺紋鋼的伸長量;為螺紋鋼的兩端距離。
將式(2)、式(3)代入式(1)后變換得出螺紋鋼伸長量為:
從圖7(a)還可知,中間變形比兩端要大,這說明在工作時(shí)上箱梁已傾斜。22孔工作錨板孔1和2之間的最大距離為132 mm(圖10),調(diào)出工作錨板上1和2兩點(diǎn)沿軸方向投影在上壓墊板兩點(diǎn)1和2的軸方向位移,節(jié)點(diǎn)ID號(hào)分別為284039和283965的軸方向位移差為0.063 mm(圖11)。
按照J(rèn)B/T 7473-2010[14]對(duì)小于1000 mm上下工作臺(tái)Ⅱ級(jí)精度等級(jí)的平行度要求:
式中:L2為測量長度。
圖11 兩節(jié)點(diǎn)高度差云圖
將上壓墊板長度570 mm代入式(5)得到平行度為0.091 mm,再轉(zhuǎn)換到長度132 mm的平行度為0.021 mm。初始分析中1和2高度差0.063 mm的結(jié)果已超出此標(biāo)準(zhǔn)要求。后續(xù)將以0.021 mm作為約束條件進(jìn)行優(yōu)化。
3.1.1 建立設(shè)計(jì)變量
共建立18個(gè)板厚變量,優(yōu)化范圍設(shè)定為1~80 mm,并設(shè)置一個(gè)離散變量,輸入了常用的板料厚度規(guī)格,優(yōu)化過程中只從離散變量厚度中選取。
3.1.2 建立約束函數(shù)
(1)調(diào)用1和2兩個(gè)節(jié)點(diǎn)的軸方向的位移響應(yīng),創(chuàng)建第二類響應(yīng),建立一個(gè)函數(shù)方程式為(1,2)=1-2,創(chuàng)建約束函數(shù)為(1,2)≤0.021 mm。
(2)根據(jù)不同的材料安全系數(shù)設(shè)定最大應(yīng)力值,如表2所示。
表2 各材料的最大應(yīng)力設(shè)定值
(3)根據(jù)設(shè)定的拉應(yīng)力創(chuàng)建2個(gè)應(yīng)力約束函數(shù)。在建立應(yīng)力響應(yīng)時(shí),采用殼單元,因?yàn)闅卧袃擅妫?dāng)殼單元受力彎曲時(shí),一面為拉應(yīng)力、一面為壓應(yīng)力、中面應(yīng)力為0。Q355屬塑性好的材料,受壓一般不會(huì)被破壞[15],因此優(yōu)化時(shí)關(guān)注的只是拉應(yīng)力的大小,應(yīng)力類型選擇為1。
3.1.3 建立目標(biāo)函數(shù)
目標(biāo)函數(shù)為體積最小化。
第一次優(yōu)化經(jīng)過30次迭代收斂,體積為4.67E+08 mm3,如圖12所示,分析完成顯示為不可行設(shè)計(jì),至少違反了一個(gè)約束。
圖12 第一次優(yōu)化迭代圖
兩節(jié)點(diǎn)軸方向高度差0.031 mm(圖13),螺紋鋼最大應(yīng)力383.1 MPa(圖14),這兩項(xiàng)未滿足要求。
圖13 兩節(jié)點(diǎn)高度差云圖
圖14 螺母鋼應(yīng)力云圖
通過第一次優(yōu)化結(jié)果可以看出,即使增加箱梁的體積增大梁剛度仍無法滿足1、2兩節(jié)點(diǎn)高度差的要求,究其原因:4.9×106N頂力點(diǎn)是在端支點(diǎn)和中間支點(diǎn)的正中位置,中間的螺紋鋼支點(diǎn)需承受總的支反力與兩端支點(diǎn)合力相等,而中間螺紋鋼數(shù)量與兩端相同都是12根,而中間比兩端的變形量沒有增大1倍,是由于箱體的鋼度抵消了一部分力。優(yōu)化過程中,其算法是通過增加板的厚度提高箱梁剛度來抵消更多的力,使得中間螺紋鋼與兩端的受力盡量一致,從而滿足高度差要求,但板厚已優(yōu)化至厚度的上限(圖15),仍無法滿足要求。
圖15 第一次優(yōu)化后板厚云圖
要解決以上問題可從幾個(gè)方面考慮:
(1)增加中間的螺紋鋼的數(shù)量,但多增加一排6根螺紋鋼,3排的中間1排螺紋鋼在擰緊螺母會(huì)有困難,因工藝問題暫不考慮;
(2)加大中間的螺紋鋼直徑、減小靠兩端的直徑,從而減小中間和兩端螺紋鋼的拉伸變形量的差值;
(3)更改受力點(diǎn)位置,將受力點(diǎn)向端部偏移,會(huì)使端部的螺紋鋼多承載些力,讓中間和端部的螺紋鋼承載力差距減小。
(1)使用移動(dòng)工具將工作錨板、上壓墊板和下壓墊板的殼單元往軸方向移動(dòng)-20 mm,偏離中心向端部移動(dòng),上端載荷和下端的邊界約束自動(dòng)跟隨。
(2)增加設(shè)計(jì)變量:以螺紋鋼的半徑作為設(shè)計(jì)變量,因空間限制,最大可以用到40 mm螺紋鋼,優(yōu)化范圍設(shè)定為12.5~20 mm,再設(shè)置一個(gè)離散變量,優(yōu)化過程中只允許選用半徑規(guī)格為12.5 mm、16 mm、18 mm、20 mm的螺紋鋼。螺紋鋼分為9組,分組情況如圖16所示;
圖16 螺紋鋼分組圖
(3)其余設(shè)計(jì)變量、約束函數(shù)和目標(biāo)函數(shù)保持與第一次優(yōu)化設(shè)置不變。
第二次優(yōu)化經(jīng)過18次迭代收斂,體積為2.32×108mm3,如圖17所示,分析完成顯示為可行設(shè)計(jì),滿足所有約束。
圖17 第二次優(yōu)化迭代圖
兩節(jié)點(diǎn)軸方向高度差0.019 mm(圖18),各板最大拉應(yīng)力170.8 MPa(圖19),螺紋鋼最大應(yīng)力328.6 MPa(圖20),所有約束條件滿足要求。
圖18 兩節(jié)點(diǎn)高度差云圖
圖19 各板應(yīng)力云圖
圖20 螺紋鋼應(yīng)力云圖
各板厚度優(yōu)化前后對(duì)照如表3、圖21所示。上、下箱梁外形尺寸是一樣的,優(yōu)化后上箱梁的用料比下箱梁多很多,其主要原因是上箱梁受力面積比下箱梁小,工況有所不同。同時(shí)也說明了HyperWorks在優(yōu)化過程中可以做到在滿足約束函數(shù)的前提下,通過增加板厚來彌補(bǔ)強(qiáng)度和剛度的不足,而對(duì)起作用較小的板通過減小板厚達(dá)到體積最小化的目的。
表3 優(yōu)化前后各板厚度
注:優(yōu)化前各板厚度均為40 mm。
螺紋鋼優(yōu)化前后對(duì)照如表4、圖22所示。可以看出,中間的1~3組已經(jīng)優(yōu)化至最大,而端部的4~9組也相應(yīng)地減小,優(yōu)化出了最優(yōu)結(jié)果。9組4個(gè)變量的排列組合有很多種,即使具有一定的理論能力,按照傳統(tǒng)方法去嘗試也會(huì)非常地費(fèi)時(shí)費(fèi)力,還不一定得到最優(yōu)結(jié)果,從這可以看出HyperWorks的優(yōu)化速度的優(yōu)勢(shì)是很明顯的。
表4 優(yōu)化前后螺紋鋼半徑
注:優(yōu)化前各組螺紋鋼半徑均為16 mm。
圖21 各板厚度優(yōu)化前后對(duì)照?qǐng)D
圖22 螺紋鋼半徑優(yōu)化前后對(duì)照?qǐng)D
從表3優(yōu)化前、后對(duì)比結(jié)果可知,在滿足預(yù)期設(shè)計(jì)目標(biāo)的前提下,優(yōu)化后體積較優(yōu)化前減少了38.0%。
通過HyperWorks軟件對(duì)反力架進(jìn)行了優(yōu)化分析,最終實(shí)現(xiàn)了最優(yōu)的設(shè)計(jì)方案,并得到如下結(jié)論:
(1)HyperWorks尺寸優(yōu)化功能強(qiáng)大,通過第一次優(yōu)化失敗和第二次優(yōu)化得到最優(yōu)結(jié)果的過程可以看出:只要約束函數(shù)、設(shè)計(jì)變量和目標(biāo)函數(shù)設(shè)置合理,就能較快地得出最優(yōu)化設(shè)計(jì)方案;
(2)對(duì)于大型結(jié)構(gòu)件,盡可能使用Shell單元和Beam單元替代薄板件和桿件的實(shí)體單元,可提減少網(wǎng)格數(shù)量,提高計(jì)算效率,而且可以較好地利用HyperWorks的尺寸優(yōu)化功能進(jìn)行參數(shù)化優(yōu)化設(shè)計(jì);
(3)通過第二類響應(yīng)自定義函數(shù)可較好地解決直接調(diào)用第一類響應(yīng)無法滿足所需要的響應(yīng)的問題,對(duì)需要用第二類響應(yīng)創(chuàng)建更為復(fù)雜的約束函數(shù)或目標(biāo)函數(shù)具有一定的參考意義。
表3 優(yōu)化效果表
[1]柴振嶺,郭子雄,黃群賢. 拉壓雙作用梁式反力架的研制[J]. 實(shí)驗(yàn)室研究與探索,2015,34(3):23-26,62.
[2]王凱,張良輝,孫振川,等. 基于PSO算法的反力架設(shè)計(jì)與優(yōu)化[J]. 隧道建設(shè),2019,39(2):302-308.
[3]方獻(xiàn)軍,徐自立,熊春明. OptiStruct及HyperStudy優(yōu)化與工程應(yīng)用[M]. 北京:機(jī)械工業(yè)出版社,2021.
[4]張積洪,杜陽,龐玥,等. 基于HyperWorks的機(jī)場某拖車車架結(jié)構(gòu)的輕量化研究[J]. 機(jī)械設(shè)計(jì),2014,31(5):38-41.
[5]張時(shí)鋒,許小鵬,陶智森. 余熱鍋爐旁路擋板門殼體尺寸優(yōu)化研究[J]. 機(jī)電工程,2019,36(3):244-250.
[6]徐晶,劉基岡,張懿驊. 基于HyperMesh的液壓機(jī)機(jī)身結(jié)構(gòu)優(yōu)化[J]. 鍛壓技術(shù),2015,40(11):48-52.
[7]中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局中國國家標(biāo)準(zhǔn)化管理委員會(huì). 預(yù)應(yīng)力混凝土用螺紋鋼筋:GB/T 20065-2016[S].
[8]郭崇志,鄭彬. 基于Shell單元的壓力容器結(jié)構(gòu)安全與性能評(píng)價(jià)[J]. 石油機(jī)械,2009,37(1):73-77.
[9]董永朋,王佩艷,劉洋等. 采用Beam單元模擬復(fù)合材料加筋桁條的優(yōu)化分析[J]. 材料科學(xué)與工程學(xué)報(bào),2010,28(6):857-861,870.
[10]沈春根,孔維忠,關(guān)天龍. UG NX 11.0有限元分析基礎(chǔ)實(shí)戰(zhàn) [M]. 北京:機(jī)械工業(yè)出版社,2018.
[11]彭春雷,吳雪飛,馬本學(xué). 基于ANSYS中的對(duì)稱分析及其應(yīng)用[J]. 機(jī)械設(shè)計(jì)與制造,2007(12):68-69.
[12]馬慶豐,汪謨清,李戈操. 基于HyperWorks的全向側(cè)面叉車車架結(jié)構(gòu)有限元分析[J]. 機(jī)電工程技術(shù),2010,39(5):57-60.
[13]江鵬. 飛機(jī)外蒙皮上非受力結(jié)構(gòu)有限元模型設(shè)計(jì)淺析[J]. 科技創(chuàng)新導(dǎo)報(bào),2011(12):45-46.
[14]中華人民共和國工業(yè)和信息化部. 單雙動(dòng)薄板沖壓液壓機(jī):JB/T 7343-2010[S].
[15]劉鳴放,劉勝新. 金屬材料力學(xué)性能手冊(cè)[M]. 北京:機(jī)械工業(yè)出版社,2012.
Optimal Design of Reaction Frame Based on HyperWorks
LIAO Rong,YANG Fan,XU Wei,ZHANG Zhenqi
(LiuzhouTaimu Prestressing Force MachineryCo.,Ltd., Liuzhou 545000, China )
The size optimization function of HyperWorks is used to analyze the optimization of the reaction frame of the stressed parts of the extruder. Through the static analysis, the first optimization, readjustment of optimization parameters and the second optimization, the requirements of the objective function and the constraint function are finally met, and the optimal design scheme is obtained. Meanwhile, the shell element and the beam element are used to replace the solid element to simplify the model, and its 1/2 model is used for symmetrical analysis, which greatly reduces the number of meshes and improves the operation speed and efficiency. The second type of response function of HyperWorks is used to create constraint functions. The second type of response is the algebraic expression of the first type of response, which better solves the problem that the first type of response cannot be directly called. The material needed of the final optimization result is 38.0% less than the initial scheme, which not only meets the design requirements of strength and stiffness, but also reduces the cost.
reaction frame;finite element;symmetrical analysis;the second type of response
TH123.3
A
10.3969/j.issn.1006-0316.2023.07.011
1006-0316 (2023) 007-0072-09
2022-08-08
廖嶸(1974-),男,江西龍南人,工程師,主要從事預(yù)應(yīng)力設(shè)備及產(chǎn)品設(shè)計(jì)、研發(fā)工作,E-mail:2660255911@qq.com。