林冠屹 管殿柱 張開(kāi)拓 宋占杰

摘要：??針對(duì)某公司設(shè)計(jì)的一款龍門式折彎?rùn)C(jī)，本文在其正式投入生產(chǎn)與工作之前進(jìn)行計(jì)算機(jī)模擬分析，在驗(yàn)證設(shè)計(jì)方案可行性的同時(shí)進(jìn)行輕量化設(shè)計(jì)。首先通過(guò)Solid?Works對(duì)折彎?rùn)C(jī)進(jìn)行三維建模，模型簡(jiǎn)化后導(dǎo)入ANSYS?Workbench15.0，通過(guò)對(duì)主要受力部位施加載荷與壓強(qiáng)，劃分網(wǎng)格，在最大應(yīng)力小于材料屈服極限值，即保證材料性能的前提下，得出整體變形的分析結(jié)果。根據(jù)折彎?rùn)C(jī)變形情況，設(shè)計(jì)結(jié)構(gòu)優(yōu)化方案，將工作臺(tái)側(cè)板厚度、工作臺(tái)支撐板厚度、工作臺(tái)頂板厚度、折彎?rùn)C(jī)拉桿直徑與上壓板筋板厚度5個(gè)對(duì)分析結(jié)果影響顯著的結(jié)構(gòu)參數(shù)添加進(jìn)來(lái)，進(jìn)行正交試驗(yàn)，根據(jù)試驗(yàn)組號(hào)重新建模進(jìn)行有限元分析。以折彎?rùn)C(jī)質(zhì)量最輕為目標(biāo)函數(shù)，設(shè)置狀態(tài)變量及尺寸變量的邊界條件，使用SPSS統(tǒng)計(jì)分析得出系數(shù)，建立折彎?rùn)C(jī)輕量化數(shù)學(xué)模型，運(yùn)用Matlab優(yōu)化工具箱進(jìn)行求解，將各個(gè)參數(shù)圓整后重新建模并進(jìn)行有限元分析，最終使折彎?rùn)C(jī)整機(jī)質(zhì)量減輕了3.36%，同時(shí)保證了折彎?rùn)C(jī)的工作性能，優(yōu)化效果顯著。該研究具有一定的實(shí)際應(yīng)用價(jià)值。

關(guān)鍵詞：??正交試驗(yàn)法;?輕量化設(shè)計(jì);?龍門式折彎?rùn)C(jī);?ANSYS?Workbench;?有限元分析

中圖分類號(hào)：?TH122?文獻(xiàn)標(biāo)識(shí)碼：?A

收稿日期：?2019-09-30;?修回日期：?2019-11-15

作者簡(jiǎn)介：??林冠屹（1993-），男，碩士研究生，主要研究方向?yàn)閿?shù)字化設(shè)計(jì)及虛擬樣機(jī)技術(shù)。

通信作者：??管殿柱（1969-），男，碩士，教授，主要研究方向?yàn)閿?shù)字化設(shè)計(jì)及虛擬樣機(jī)技術(shù)。Email：?gdz_zero@126.com

折彎?rùn)C(jī)作為鈑金加工業(yè)的一種重要母機(jī)，廣泛應(yīng)用于鈑金彎曲成型加工，在汽車、金屬結(jié)構(gòu)、電子、簡(jiǎn)述裝潢等領(lǐng)域廣泛應(yīng)用[1-4]。折彎?rùn)C(jī)的工作原理是將作用力作用于壓板或滑塊上，使其連續(xù)上下運(yùn)動(dòng)，將固定工作臺(tái)上的板料或工件強(qiáng)制折彎成特定形狀。折彎?rùn)C(jī)的驅(qū)動(dòng)方式分為氣動(dòng)式、機(jī)械式和液壓式，其中應(yīng)用最廣泛的是液壓式折彎?rùn)C(jī)[5]。目前，市場(chǎng)上最常見(jiàn)的是C型折彎?rùn)C(jī)，其結(jié)構(gòu)設(shè)計(jì)簡(jiǎn)單，易于制造，但只能對(duì)某種截面形狀的材料進(jìn)行折彎，且折彎時(shí)喉口處存在應(yīng)力過(guò)大的問(wèn)題[6]。利用傳統(tǒng)設(shè)計(jì)方法自主設(shè)計(jì)的龍門式液壓折彎?rùn)C(jī)，雖然具有較大的安全系數(shù)，但其結(jié)構(gòu)設(shè)計(jì)比較笨重，增加了生產(chǎn)成本，需要進(jìn)行輕量化設(shè)計(jì)[7]。基于此，本文以某公司設(shè)計(jì)的一款龍門式折彎?rùn)C(jī)為模型，在對(duì)模型進(jìn)行靜力學(xué)分析的基礎(chǔ)上，設(shè)置目標(biāo)函數(shù)和邊界條件，采用正交試驗(yàn)方法，運(yùn)用Matlab軟件進(jìn)行求解，進(jìn)而得到折彎?rùn)C(jī)輕量化結(jié)構(gòu)，最終使折彎?rùn)C(jī)整機(jī)質(zhì)量減輕了3.36%，優(yōu)化效果明顯。該研究具有一定的實(shí)際應(yīng)用價(jià)值。

1?有限元分析

1.1?折彎?rùn)C(jī)三維模型分析

本文所用龍門式折彎?rùn)C(jī)使用單缸壓力為8?000?kN的液壓缸為動(dòng)力元件，設(shè)計(jì)的最大工作壓力為80?000?kN。折彎?rùn)C(jī)工作臺(tái)和四周的立柱及上部的油缸機(jī)架構(gòu)成了折彎?rùn)C(jī)外部結(jié)構(gòu)，起固定、支撐及連接作用[8]。超大型折彎?rùn)C(jī)的固定方式比較特殊，工作臺(tái)通過(guò)地腳螺栓與地面連接固定，工作臺(tái)底部嵌入地下，看作全約束。工作臺(tái)與油缸機(jī)架通過(guò)4個(gè)立柱連接，每個(gè)立柱內(nèi)部有4根拉桿，為折彎?rùn)C(jī)工作提供有效拉力，有效防止折彎?rùn)C(jī)機(jī)身形變位移過(guò)大。底模與工作臺(tái)頂板固定連接，上壓板與液壓缸固定連接，通過(guò)控制液壓缸的進(jìn)給量將置于上壓板與底模之間的板材折彎成不同角度。龍門式折彎?rùn)C(jī)三維模型示意圖如圖1所示。

1.2?折彎?rùn)C(jī)有限元模型分析

對(duì)折彎?rùn)C(jī)模型進(jìn)行有限元分析時(shí)，需要對(duì)折彎?rùn)C(jī)的計(jì)算機(jī)輔助設(shè)計(jì)（computer?aided?design，CAD）模型進(jìn)行適當(dāng)簡(jiǎn)化和修改，否則在轉(zhuǎn)換為計(jì)算機(jī)輔助求解（computer?aided?engineering，CAE）模型時(shí)可能會(huì)出現(xiàn)某些問(wèn)題[9]。針對(duì)龍門式折彎?rùn)C(jī)結(jié)構(gòu)的特點(diǎn)，對(duì)其模型作如下簡(jiǎn)化。刪除部分與折彎?rùn)C(jī)強(qiáng)度和剛度無(wú)關(guān)的螺紋及螺紋孔等結(jié)構(gòu)，將工作臺(tái)底部約束視為固定約束，將焊接連接強(qiáng)度視為與整體材料相同，去除對(duì)整體受力影響較小的部件。在不影響分析誤差的前提下進(jìn)行簡(jiǎn)化，確保分析的正常進(jìn)行。

對(duì)模型進(jìn)行網(wǎng)格劃分時(shí)應(yīng)遵循均勻、應(yīng)力區(qū)粗劃、應(yīng)力梯度細(xì)化的原則[10-11]，合理的網(wǎng)格劃分可提高有限元分析結(jié)果的準(zhǔn)確性，但網(wǎng)格數(shù)量的增加會(huì)加大計(jì)算機(jī)的運(yùn)算負(fù)擔(dān)，從而增加運(yùn)算時(shí)間。本文對(duì)折彎?rùn)C(jī)模型的網(wǎng)格劃分采用四面體網(wǎng)格劃分的方法，有效避免以上問(wèn)題。網(wǎng)格劃分完成后，得到模型單元數(shù)為129?693，節(jié)點(diǎn)數(shù)為550?029。

1.3?折彎?rùn)C(jī)材料及參數(shù)設(shè)置

龍門式折彎?rùn)C(jī)主體機(jī)身材料為Q235，鋼板采用焊接方式連接，折彎?rùn)C(jī)拉桿材料采用40CrMo，經(jīng)查閱手冊(cè)確定其材料參數(shù)，材料基本參數(shù)如表1所示。

折彎?rùn)C(jī)工作時(shí)，所受載荷實(shí)際上隨時(shí)間發(fā)生線性變化，在進(jìn)行靜力學(xué)分析時(shí)，僅考慮滿載時(shí)的結(jié)構(gòu)受力情況，為保證其最大受力極限工況下的安全性能，取折彎?rùn)C(jī)滿載時(shí)的工況為分析場(chǎng)景，根據(jù)初試設(shè)計(jì)要求，在折彎?rùn)C(jī)工作時(shí)，由固定在上機(jī)架上的10個(gè)公稱力為8?000?kN的液壓缸向下推動(dòng)壓塊完成板材折彎[12]，上機(jī)架承受液壓缸底部向上的反作用力F為

F=800?000?kg×9.8?N/kg=7.84×106?N

根據(jù)該龍門型折彎?rùn)C(jī)的實(shí)際技術(shù)參數(shù)，分別計(jì)算液壓缸所受壓強(qiáng)P1，凸臺(tái)所受壓強(qiáng)P2，壓板所受壓強(qiáng)P3，工作臺(tái)所受壓強(qiáng)P4。

根據(jù)壓強(qiáng)計(jì)算公式

P=F/A（1）

可得液壓缸所受壓強(qiáng)為

P1=7.84×106?N565?072.34?mm2≈13.874?MPa

凸臺(tái)壓強(qiáng)為

P2=7.84×106?N501?976.82?mm2≈15.618?MPa

壓板壓強(qiáng)為

P3=7.84×107?N25×121?500?mm2≈25.811?MPa

工作臺(tái)壓強(qiáng)為

P4=7.84×107?N9?500?000?mm2≈8.253?MPa

根據(jù)拉桿預(yù)緊力公式

F0≤0.5～0.6σsA1（2）

式中，σs為螺栓材料屈服極限，MPa;A1為螺栓危險(xiǎn)截面面積，A1≈（πd21/4）mm2。

分析驗(yàn)證結(jié)果表明，預(yù)緊力系數(shù)的值在0.05～0.5之間時(shí)，應(yīng)力隨系數(shù)的增加而增加?？紤]到折彎?rùn)C(jī)的震動(dòng)，建議以拉緊螺桿公稱力的1%～5%倍在裝配時(shí)進(jìn)行預(yù)緊[13]，拉桿直徑d=200?mm，取預(yù)緊力系數(shù)為1%，通過(guò)計(jì)算得拉桿預(yù)緊力約為292?000?N。

1.4?折彎?rùn)C(jī)靜力學(xué)分析結(jié)果

將折彎?rùn)C(jī)模型導(dǎo)入ANSYS?Workbench模塊，參數(shù)設(shè)置完成后，進(jìn)行靜力學(xué)分析，得到龍門式折彎?rùn)C(jī)機(jī)身分析云圖，折彎?rùn)C(jī)靜力學(xué)分析云圖如圖2所示。由圖2可以看出，折彎?rùn)C(jī)上壓板部分在額定工作載荷下變形量最大，最大值為5.03?mm;工作臺(tái)固定板與支撐板連接處有應(yīng)力集中現(xiàn)象，應(yīng)力最大值為209.72?MPa，符合實(shí)際工況，但機(jī)身質(zhì)量較大，需要進(jìn)一步優(yōu)化。

2?正交試驗(yàn)

2.1?正交試驗(yàn)設(shè)計(jì)

正交試驗(yàn)設(shè)計(jì)方法利用正交性挑選具有代表性的參數(shù)作為變量，建立試驗(yàn)因素水平表進(jìn)行組合試驗(yàn)，通過(guò)正交試驗(yàn)法可有效減少了試驗(yàn)次數(shù)[14-16]。由折彎?rùn)C(jī)靜力學(xué)分析可知，工作臺(tái)側(cè)板厚度X1、工作臺(tái)支撐板厚度X2、工作臺(tái)頂板厚度X3的選取直接影響折彎?rùn)C(jī)變形和應(yīng)力，同時(shí)上壓板筋板厚度X4、拉桿直徑X5取值的不同也會(huì)造成機(jī)身應(yīng)力和最大變形的變化，對(duì)其選取合理的厚度至關(guān)重要。因此，使用X1～X5這5個(gè)參數(shù)進(jìn)行正交試驗(yàn)，試驗(yàn)因素水平表如表2所示。選用L16（45）型正交表[17]，使用ANSYS?Workbench依次對(duì)16種結(jié)構(gòu)進(jìn)行靜力學(xué)分析，正交試驗(yàn)分析結(jié)果如表3所示。

2.2?回歸函數(shù)的構(gòu)造

通過(guò)正交試驗(yàn)，得到不同尺寸模型對(duì)應(yīng)的最大應(yīng)力、最大變形和機(jī)身質(zhì)量數(shù)據(jù)，設(shè)置因變量與自變量分別得到基于最小二乘法的最大應(yīng)力、最大變形、機(jī)身質(zhì)量的線性回歸系數(shù)，將表3的數(shù)據(jù)輸入SPSS18.0中[18]，得到各線性擬合方程。

機(jī)身質(zhì)量擬合方程為

y1=730?218.25+774.45x1+95.35x2+329.675x3+283.5x4+537.75x5（3）

最大應(yīng)力擬合方程為

y2=327.566-0.505x1+0.594x2-0.585x3+1.233x4-0.608x5（4）

最大變形擬合方程為

y3=4.901+0.02x1+0.03x2-0.01x3-0.05x4+0.01x5（5）

3?輕量化模型設(shè)計(jì)

3.1?設(shè)計(jì)變量與目標(biāo)函數(shù)的確定

根據(jù)本次正交試驗(yàn)選擇的設(shè)計(jì)變量，其集合為x={x1，x2，x3，x4，x5}，本次試驗(yàn)的最終目標(biāo)是求得符合應(yīng)力要求和形變要求的最小機(jī)身質(zhì)量，根據(jù)折彎?rùn)C(jī)機(jī)身質(zhì)量擬合方程（3），得到其目標(biāo)函數(shù)為

minF（x）=730?218.25+774.45x1+95.35x2+329.675x3+283.5x4+537.75x5（6）

1）?強(qiáng)度狀態(tài)變量約束。龍門式折彎?rùn)C(jī)主要材料為Q235，屈服極限為235?MPa，取安全系數(shù)為1.05，得到在進(jìn)行優(yōu)化設(shè)計(jì)時(shí)機(jī)身最大應(yīng)力不能超過(guò)223.81?MPa，因此其應(yīng)力約束為

g（1）=y2-223.81=-0.505x1+0.594x2-0.585x3+1.233x4-0.608x5+113.926（7）

2）?剛度狀態(tài)變量約束。在保證材料強(qiáng)度滿足要求的情況下，本次優(yōu)化的龍門型折彎?rùn)C(jī)要求最小變形小于5.5?mm，因此其變形約束為

g（2）=y3-5.5=0.02x1+0.03x2-0.01x3-0.05x4+0.01x5-0.599（8）

3）?幾何約束為

g（3）=x1-160，?g（4）=130-x1，?g（5）=x2-160，?g（6）=130-x2

g（7）=x3-170，?g（8）=110-x3，?g（9）=x4-55，?g（10）=40-x4（9）

g（11）=x5-210，?g（12）=180-x5

3.2?優(yōu)化函數(shù)數(shù)學(xué)模型的建立

由以上條件，得到折彎?rùn)C(jī)機(jī)身輕量化的優(yōu)化函數(shù)數(shù)學(xué)模型為

min?730?218.25+774.45x1+95.35x2+329.675x3+283.5x4+537.75x5

s.t.?g（i）≤0，?i=1，2，…，12xj>0，?j=1，2，…，5（10）

3.3?輕量化計(jì)算

為滿足輕量化要求，根據(jù)KuhnTucker條件[19]，得

ΔfX-∑12i=1ωiΔSX=0ωiSX=0，??i=1，2，…，12;?ωi≥0，?i=1，2，…，12（11）

利用Matlab軟件進(jìn)行數(shù)學(xué)模型計(jì)算，得出其計(jì)算結(jié)果，圓整后選用X1=130?mm，X2=132?mm，X3=112?mm，X4=42?mm，X5=181?mm。

4?優(yōu)化結(jié)果分析

按照優(yōu)化圓整后的尺寸修改模型，進(jìn)行有限元靜力學(xué)分析，得到優(yōu)化后的龍門式折彎?rùn)C(jī)分析結(jié)果云圖如圖3所示。

由圖3可以看出，該參數(shù)下的折彎?rùn)C(jī)模型應(yīng)力符合材料屈服強(qiáng)度的要求，變形量增加較小，機(jī)身質(zhì)量減小。將折彎?rùn)C(jī)輕量化設(shè)計(jì)前后的有限元分析結(jié)果進(jìn)行對(duì)比，計(jì)算得優(yōu)化結(jié)果對(duì)比表如表4所示。

5?結(jié)束語(yǔ)

本文通過(guò)對(duì)折彎?rùn)C(jī)進(jìn)行有限元分析，確定對(duì)折彎?rùn)C(jī)性能影響較大的結(jié)構(gòu)參數(shù)，提出正交試驗(yàn)的輕量化設(shè)計(jì)思路，設(shè)置5個(gè)變量作為正交試驗(yàn)參數(shù)，得到機(jī)身參數(shù)的回歸模型，建立優(yōu)化數(shù)學(xué)模型，對(duì)折彎?rùn)C(jī)進(jìn)行了輕量化設(shè)計(jì)。由分析設(shè)計(jì)結(jié)果可以看出，優(yōu)化之后的模型質(zhì)量減輕3.36%，應(yīng)力、變形量有一定增加，符合設(shè)計(jì)要求，達(dá)到了輕量化的目的，降低了企業(yè)的生產(chǎn)成本。該設(shè)計(jì)將多個(gè)設(shè)計(jì)變量作為優(yōu)化參數(shù)，得到輕量化結(jié)果較為顯著，該方法普遍適用于機(jī)械設(shè)備的輕量化設(shè)計(jì)。

參考文獻(xiàn)：

[1]?顧忠新，?胡智明.?基于ANSYS折彎?rùn)C(jī)機(jī)架的有限元分析[J].?鍛壓裝備與制造技術(shù)，?2014，?49（1）：?26-28.

[2]?周歡，?王麗麗，?王金榮，?等.?折彎?rùn)C(jī)液壓補(bǔ)償裝置對(duì)工作臺(tái)撓度曲線的影響分析與優(yōu)化[J].?機(jī)床與液壓，?2016，?44（20）：?40-43.

[3]?周婉，?錢怡，?平東良.?折彎?rùn)C(jī)的有限元建模及強(qiáng)度和變形分析[J].?江南大學(xué)學(xué)報(bào)：?自然科學(xué)版，?2014，?13（2）：?179-183.

[4]?田萬(wàn)英，?高建和，?潘志華.?折彎?rùn)C(jī)壓力不均勻分布研究[J].?鍛壓裝備與制造技術(shù)，?2011，?46（6）：?44-47.

[5]?范宏才.?現(xiàn)代鍛壓機(jī)械[M].?北京：?機(jī)械工業(yè)出版社，?1994.

[6]?嚴(yán)慶光，?李明哲，?蔡中義.?C型號(hào)多點(diǎn)成形壓機(jī)機(jī)架結(jié)構(gòu)有限元分析和優(yōu)化設(shè)計(jì)[J].?鍛壓機(jī)械，?2002，?37（3）：?44-47.

[7]?翟慶波，?蘇兆婧.?輕量化設(shè)計(jì)在工程機(jī)械中的應(yīng)用[J].?設(shè)計(jì)，?2016（17）：?120-121.

[8]?崔蓬勃.?500噸折彎?rùn)C(jī)的結(jié)構(gòu)與優(yōu)化[D].?青島：?青島大學(xué)，?2018.

[9]?劉祥.?WF67Y-160/3200數(shù)控折彎?rùn)C(jī)有限元分析[J].?科技與生活，?2012?（21）：?170-171.

[10]?胡亞民，?夏華，?陳善民，?等.?我國(guó)彎曲成形技術(shù)的發(fā)展?fàn)顩r[J].?鍛造與沖壓，?2005（4）：?26.

[11]?王宏，?劉翠.?折彎?rùn)C(jī)機(jī)架變形應(yīng)力的有限元分析[J].?重型機(jī)械，?2007（5）：?56-58.

[12]?張開(kāi)拓，?管殿柱，?白碩瑋，?等.?基于響應(yīng)面分析法的龍門式折彎?rùn)C(jī)輕量化設(shè)計(jì)[J].?制造業(yè)自動(dòng)化，?2019，?41（1）：?57-60.

[13]?范如明，?李恒，?吉桂生，?等.?模鍛壓力機(jī)機(jī)身結(jié)構(gòu)與剛度分析[J].?鍛壓裝備與制造技術(shù)，?2016，?51（4）：?19-22.

[14]?鞠家全，?邱自學(xué)，?崔德友，?等.?基于正交試驗(yàn)的機(jī)床移動(dòng)橫梁多目標(biāo)優(yōu)化設(shè)計(jì)[J].?機(jī)械強(qiáng)度，?2018，?40（2），?356-362.

[15]?宋明順，?黃佳，?張士朋，?等.?多指標(biāo)正交試驗(yàn)設(shè)計(jì)去量綱準(zhǔn)則及方法研究[J].?工業(yè)工程與管理，?2014，?19（1）：?41-46.

[16]?王巖，?隋思漣.?試驗(yàn)設(shè)計(jì)與MATLAB數(shù)據(jù)分析[M].?北京：?清華大學(xué)出版社，?2012.

[17]?趙韓，?陳興玉，?董玉德，?等.?參數(shù)化有限元分析的液壓機(jī)機(jī)身輕量化設(shè)計(jì)[J].?工程圖學(xué)學(xué)報(bào)，?2010，?31（1）：?20-25.

[18]?何曉群，?劉文卿.?應(yīng)用回歸分析[M].?3版.?北京：?中國(guó)人民大學(xué)出版社，?2011.

[19]?陳寶林.?最優(yōu)化理論與算法[M].?2版.?北京：?清華大學(xué)出版社，?2005.

Lightweight?Design?of?Gantry?Bender?Based?on?Orthogonal?Test?Method

LIN?Guanyi，?GUAN?Dianzhu，?ZHANG?Kaituo，?SONG?Zhanjie

（School?of?Electromechanic?Engineering，?Qingdao?University，?Qingdao?266071，?China）

Abstract：??Aiming?at?the?design?of?a?gantry?bending?machine，?computer?simulation?analysis?is?needed?before?it?is?put?into?production?and?work.?The?feasibility?of?the?design?scheme?is?verified?and?the?lightweight?design?is?carried?out?at?the?same?time.?Firstly，?the?threedimensional?model?of?the?bending?machine?is?built?by?SolidWorks，?which?is?imported?into?ANSYS?Workbench?15.0.?By?applying?load?and?pressure?on?the?main?parts?of?the?bending?machine?and?dividing?the?mesh，?the?analysis?results?of?the?overall?deformation?of?the?bending?machine?are?obtained?on?the?premise?that?the?maximum?stress?of?the?material?is?less?than?the?material?yield?limit?value.?According?to?the?deformations?of?the?bending?machine，?the?structural?optimization?scheme?is?designed.?Five?parameters，?namely，?the?thickness?of?the?side?plate?of?the?table，?the?thickness?of?the?supporting?plate?of?the?table，?the?thickness?of?the?roof?of?the?table，?the?diameter?of?the?pull?rod?of?the?bending?machine?and?the?thickness?of?the?stiffened?plate?of?the?upper?pressure?plate，?are?added.?Orthogonal?experiments?are?carried?out，?and?the?finite?element?analysis?is?done?according?to?the?test?group?numbers.?The?lightest?weight?of?bending?machine?is?the?objective?function.?The?boundary?conditions?of?state?variables?and?size?variables?are?set.?The?lightweight?model?of?bending?machine?is?established.?The?model?is?solved?by?using?Matlab?software.?The?parameters?are?rounded?and?remodeled?and?analyzed?by?finite?element?method.?Finally，?the?weight?of?the?whole?bending?machine?is?reduced?by?42?510?kg，?and?the?working?performance?of?the?bending?machine?is?guaranteed.?The?optimization?effect?is?remarkable.?This?research?has?certain?practical?application?value.

Key?words：??orthogonal?test?method;?lightweight?design;?gantry?bender;?ANSYS?Workbench;?finite?element?analysis