有限元法在折疊紙盒強度分析中的應(yīng)用

時間：2023-03-17 06:36:06 理工畢業(yè)論文我要投稿

相關(guān)推薦

　　摘要：介紹了有限元分析的方法，討論了有限元法在折疊紙盒強度分析中的應(yīng)用，并且結(jié)合實際研究提出了有限元法在分析紙盒強度時要考慮的系列問題，同時為折疊紙盒強度的研究提供了一種新思路，推動了計算機輔助分析在包裝工程方面的應(yīng)用。

　　關(guān)鍵詞：折疊紙盒強度等效模擬有限元法

　　折疊紙盒憑借其原材料的環(huán)保性、可裝潢性、可印刷性等其他包裝不可替代的特性[1]，在物資流通和包裝行業(yè)深受人們的青睞，廣泛應(yīng)用于食品、藥品、化妝品、煙酒等日常生活用品的包裝。盡管折疊紙盒不是運輸包裝容器，但作為一種“剛”性包裝容器而言，折疊紙盒也應(yīng)具有足夠的強度來抵抗變形和破壞的力量，才能對商品提供足夠的保護和便于銷售時對商品進行陳列和展示[2]。

　　近年來，國內(nèi)外學者對折疊紙盒的抗壓強度進行了不少研究，王金林[3]論證了紙盒抗壓強度和白板紙挺度的關(guān)系;Grangard[4]通過大量的試驗確定經(jīng)驗公式BRDA式具有較高的相關(guān)性;余本農(nóng)[2]根據(jù)瓦楞紙箱抗壓強度公式推導出折疊紙盒的經(jīng)驗公式;劉慧和張新昌[5-6]研究了管式折疊紙盒的結(jié)構(gòu)和幾何尺寸對抗壓強度的影響，并且探討折疊紙盒抗壓強度經(jīng)驗計算公式;袁毅等[7]研究鎖底結(jié)構(gòu)形式對紙盒承載力的影響。但以上學者的研究主要還停留在試驗研究的基礎(chǔ)上，由于試驗的不穩(wěn)定性以及各種尺寸折疊紙盒試驗樣品獲取困難，要想更好地掌握折疊紙盒承載能力和抗壓變形能力不能只停留在單純的試驗研究上，應(yīng)該從原材料自身入手，結(jié)合理論與試驗結(jié)論，借助于有限元分析技術(shù)來對折疊紙盒承載能力和抗壓變形能力進行仿真模擬，不僅可以縮短實驗周期，節(jié)約試驗費用，還能夠?qū)埡懈鱾€部位作研究和討論，更為直觀地得到折疊紙盒各個部位的力學狀態(tài)對包裝性能的影響。

　　一、有限元法的概念

　　有限單元法最初是20世紀50年代作為處理固體力學問題的方法出現(xiàn)的，并隨著電子計算機的發(fā)展而迅速發(fā)展起來的一種現(xiàn)代計算方法。作為計算機輔助工程CAE的一種，使用它大大提高了產(chǎn)品開發(fā)、設(shè)計、分析和制造的效率和產(chǎn)品性能[8]。

　　有限元法是一種基于變分法而發(fā)展起來的求解微分方程的數(shù)值計算方法，其基本思想是把一個連續(xù)的介質(zhì)(或構(gòu)件)看成是由有限數(shù)目的單元組成的集合體，并在每一個單元中設(shè)定有限數(shù)量的節(jié)點，在各單元內(nèi)假定具有一定的理想化的位移和應(yīng)力分布模式，各單元間通過節(jié)點相連接，并藉以實現(xiàn)應(yīng)力的傳遞，各單元之間的交接面要求位移協(xié)調(diào)，利用彈性力學、固體力學、結(jié)構(gòu)力學等力學中的變分原理去建立一套線性方程組，從而將一個連續(xù)域中的無限自由度問題轉(zhuǎn)化為離散域中的有限自由度問題，求解這些方程組，便可得到各單元和結(jié)點的位移、應(yīng)力。簡言之，就是用較簡單的問題代替復(fù)雜問題后再求解，即化整為零分析，積零為整研究。對于不同物理性質(zhì)和數(shù)學模型的問題，有限元求解法的基本步驟是相同的，分為前處理、處理和后處理3個階段，只是具體公式推導和運算求解不同。

　　二、有限元法在分析紙盒強度時要考慮的問題

　　2.1 研究參數(shù)的選取問題

　　根據(jù)理論分析可以得出，折疊紙盒的抗壓強度受到紙板原材料性能以及紙盒的結(jié)構(gòu)形式和結(jié)構(gòu)參數(shù)的影響。紙板原材料性能包括紙板的環(huán)壓強度、挺度、厚度和定量等，紙盒結(jié)構(gòu)形式包括管式、盤式、管盤式和非管非盤式等，紙盒結(jié)構(gòu)參數(shù)主要指的是紙盒的幾何尺寸長、寬和高。對于研究紙盒強度問題，在單純的試驗研究時首先是確定紙板原材料和紙盒結(jié)構(gòu)形式，然后對改變不同幾何結(jié)構(gòu)參數(shù)的紙盒進行測試，看紙盒是否滿足強度和剛度的要求。

　　其實，紙盒的結(jié)構(gòu)形式和原材料性能在很大程度上決定了紙盒的特點和工藝效果，提高紙板的挺度和環(huán)壓強度，優(yōu)化紙盒的結(jié)構(gòu)和幾何尺寸都可以提高紙盒的抗壓強度，若是僅僅靠改變結(jié)構(gòu)參數(shù)來研究紙盒強度性能，很難做到對紙盒抗壓強度的系統(tǒng)化研究。而且目前對紙盒強度的研究主要是基于“設(shè)計―盒樣―試驗―修改”不斷循環(huán)的優(yōu)化過程模式，這樣既麻煩又浪費大量成本。所以在基于有限元法研究紙盒各個參數(shù)的同時還應(yīng)該考慮紙盒結(jié)構(gòu)形式、原材料屬性以及還未進行過研究的參數(shù)。另外，有限元模擬分析既可以針對某一參數(shù)變量的研究，也可以對多參數(shù)的交互作用進行研究，我們可以根據(jù)需要靈活運用。

　　2.2 有限元等效模擬分析問題

　　有限元法對紙盒壓縮過程的數(shù)值模擬就是利用一定的方法，將結(jié)構(gòu)和應(yīng)力變化的復(fù)雜過程近似模擬出來，它包括結(jié)構(gòu)靜力學的計算和結(jié)構(gòu)非線性的計算兩部分。在有限元仿真模擬分析中，幾何實體模型的建立及簡化、單元類型的選取、有限元網(wǎng)格模型的生成、約束條件的設(shè)置和載荷的施加都顯得尤為重要，它們都將直接影響到最終仿真計算結(jié)果的精確性。

　　(1)幾何模型建立及簡化問題

　　折疊紙盒的受力壓縮過程是一個實際問題，考慮到壓縮的復(fù)雜性和多樣性對有限元分析計算結(jié)果的影響，必須對紙盒結(jié)構(gòu)進行合理簡化，把實際結(jié)構(gòu)簡化成計算模型，然后再利用有限元法來分析解決問題。為了方便有限元建模及計算，在不失真的前提下，對紙盒結(jié)構(gòu)可做如下的簡化處理。

　�、勹b于紙板原材料的特性，將紙板等效成一塊正交各向異性的復(fù)合材料板;②忽略非承載件的影響，將紙盒的六個面板都看作是四邊簡支板，忽略紙盒上的小結(jié)構(gòu)，如側(cè)板上的襟片以及插舌等;③認為結(jié)構(gòu)中的粘合是理想粘合，忽略襟片粘合不牢所引起的強度弱化問題。

　　由于紙盒結(jié)構(gòu)規(guī)則，對紙盒結(jié)構(gòu)進行合理簡化之后，可以選擇直接在有限元分析軟件ANSYS中采用“自底向上”法建立模型。需要注意的是，簡化一些無關(guān)緊要的細節(jié)能使分析求解盡可能地高效，但是隨著襟片和其它一些細節(jié)被簡化，在它們鄰近區(qū)域內(nèi)仿真模擬計算出的應(yīng)力值可能不準確，導致該區(qū)域計算的結(jié)果不能反映真實應(yīng)力。因此，可以考慮采用在模型中添加該細節(jié)重新計算、子模型法和外插值法3種方案對局部應(yīng)力進行考察。

　　(2)單元類型選擇問題

　　在有限元分析過程中，對于不同的問題，需要應(yīng)用到不同特性的單元，ANSYS單元庫中提供了超過250種的不同單元類型，每一種單元都是專門為有限元問題而設(shè)計的。因此，單元選擇不當，直接影響到計算能否進行和計算結(jié)果的精度。根據(jù)紙板原材料的結(jié)構(gòu)，由于紙板的厚度相對較薄，理論上可以將其簡化為正交各向異性薄板。在薄板構(gòu)件的有限元分析中，一般采用殼體單元或者三維實體單元進行分析。根據(jù)國內(nèi)外一些研究資料和有限元分析數(shù)據(jù)，采用殼體單元對紙板結(jié)構(gòu)進行分析可以獲得足夠精確的結(jié)果，并且采用殼體單元分析可以比采用實體單元分析節(jié)省大量的時間。

　　殼單元是為有效建立薄的結(jié)構(gòu)而特別設(shè)計的單元類型，它用于主尺寸不小于10倍厚度的結(jié)構(gòu)。ANSYS中提供了一系列殼體單元，常見的有：SHELL63、SHELL43、SHELL93、SHELL181等。其中SHELL181是一個4節(jié)點的四邊形薄板殼結(jié)構(gòu)，如圖1所示，它含有I、J、K、L四個節(jié)點，每個節(jié)點含有6個自由度，即沿X、Y、Z方向上的三個平動自由度和繞X、Y、Z軸的三個轉(zhuǎn)動自由度。它可以支持材料塑性、應(yīng)力剛化、大應(yīng)變，適用于從薄到中等厚度的薄板殼結(jié)構(gòu)。因此，根據(jù)紙板本身特性和紙盒壓縮的實際過程，為了獲得更精確的模擬數(shù)據(jù)，建議選擇SHELL181單元類型進行分析。

　　(3)網(wǎng)格劃分問題

　　ANSYS的運算最終是通過有限單元法進行，因此對建立好的幾何模型進行網(wǎng)格劃分是ANSYS分析的一個重要內(nèi)容。網(wǎng)格劃分直接影響到計算結(jié)果的精度和計算速度，甚至會因為網(wǎng)格劃分不合理而導致計算不收斂。目前對于網(wǎng)格劃分，有兩種劃分方式：自由網(wǎng)格劃分和映射網(wǎng)格劃分，這兩種網(wǎng)格劃分方法各有優(yōu)缺點。自由網(wǎng)格劃分適合于初學者使用，劃分耗時較少，而映射網(wǎng)格計算結(jié)果的精度一般要高于自由網(wǎng)格的計算精度。充分考慮到紙盒各部位的結(jié)構(gòu)，可以采用自由網(wǎng)格劃分和映射網(wǎng)格劃分相結(jié)合的方法。對紙盒面板結(jié)構(gòu)使用自由網(wǎng)格劃分，在可能出現(xiàn)應(yīng)力集中的位置或擬定的重點考察位置采用映射網(wǎng)格劃分，使之能以較高的精度得到這些位置的應(yīng)力，從而保證計算結(jié)果的精度。

　　另外，網(wǎng)格數(shù)量的多少影響計算結(jié)果的精度和計算分析時間，但并非網(wǎng)格劃分的越細得到的結(jié)果就越好，因為網(wǎng)格太密太細，會占用大量的分析時間。為了考察網(wǎng)格劃分的合理程度，在計算模型建立后進行試算，以考察模型精度，特別是重點考察區(qū)域的精度。對于紙板類結(jié)構(gòu)，保證單元邊長在lmm-2mm之間時基本都能取得理想的結(jié)果。

　　(4)約束與加載問題

　　建立紙盒有限元分析模型應(yīng)當在每個坐標軸方向上至少給定一個約束，受約束節(jié)點的數(shù)量根據(jù)實際情況具體確定。定義約束條件是一個容易產(chǎn)生較大誤差的地方，紙盒的約束是根據(jù)其具體結(jié)構(gòu)確定的。紙盒的壓縮試驗是在上下兩塊光滑的壓盤上進行的，顯然壓盤的剛度比紙盒大的多，上壓盤以緩慢的速度加載，最后將紙盒壓潰。當紙盒在水平壓盤上加載受壓時由于摩擦力的作用與上壓盤接觸的紙盒面板只會沿豎直方向被壓縮，而與下壓盤接觸的紙盒面板相當于被固定在水平的壓盤上，不會有水平的滑動也沒豎直方向的壓縮。故對紙盒施加約束時，應(yīng)將紙盒受壓時的邊界約束條件等效化，根據(jù)紙盒的受力以及實際破壞情況，紙盒面板與上下壓盤接觸的各節(jié)點豎向位移耦合，與上壓盤接觸的各節(jié)點豎向位移耦合，施加共同向下的位移，下壓盤固定不動，對與下壓盤接觸的各節(jié)點的所有自由度施加約束，即對紙盒模型的四個側(cè)板的底部進行約束，保證這些部位的節(jié)點不發(fā)生平動或轉(zhuǎn)動。

　　在ANSYS分析中，載荷包括邊界條件和內(nèi)外部的作用力函數(shù)，在不同的分析領(lǐng)域中具有不同的表征，為了真實地反映實際物理情況，將載荷分為自由度約束(DOF)、集中力載荷、面載荷、體載荷、慣性載荷以及耦合場載荷等。載荷可以直接在實體模型上加載，如將載荷施加到關(guān)鍵點、線或者面上;也可以在有限元模型上加載，即直接將載荷施加到節(jié)點和單元上。這兩種方法各有利弊，但無論選用何種途徑加載，程序都會自動將這些載荷轉(zhuǎn)移到有限元模型上。在紙盒壓縮過程中，紙盒面板主要承受壓盤施加的作用力，為了模擬實際壓縮情況，應(yīng)該選擇在紙盒上加一個剛體，最后在剛體上施加一定載荷。

　　(5)非線性求解問題

　　紙盒的受力壓縮過程是一個囊括了材料、狀態(tài)、幾何三種非線性問題的復(fù)雜受力過程。材料非線性是指紙板原材料的應(yīng)力-應(yīng)變關(guān)系引起的非線性行為，而且是不可逆轉(zhuǎn)的;狀態(tài)非線性是指紙板壓縮后引起的接觸問題;幾何非線性是指紙板在壓縮的過程中其變化的幾何形狀引起的非線性響應(yīng)，屬于大應(yīng)變類型。作為復(fù)雜的非線性分析，涉及到求解器的選擇、載荷步長的選取、子步數(shù)的設(shè)置、迭代收斂準則、自動時間步長等許多非線性因素，求解設(shè)置比較復(fù)雜，若處理不當就會造成分析結(jié)果不收斂或者分析中斷等情況出現(xiàn)。因此，為了保證模擬計算的順利實現(xiàn)并獲得準確的計算結(jié)果，在前期求解設(shè)置時可以針對不同情況進行多次試算，以獲取最佳的求解設(shè)置參數(shù)，特別是控制好平衡迭代及收斂準則之間的關(guān)系。

　　三、結(jié)語

　　近年來，隨著計算機技術(shù)、現(xiàn)代力學和計算數(shù)學的發(fā)展，有限元法作為一個具有堅實理論基礎(chǔ)和廣泛應(yīng)用效力的數(shù)值模擬分析工具，在包裝工程領(lǐng)域方面的應(yīng)用越來越多。包裝設(shè)計及研究人員將這種技術(shù)引入到研究包裝材料的系統(tǒng)中，不僅提供了一種新的研究思路，而且有效地提高了工作效率，降低了包裝設(shè)計成本。但目前情況下，有限元法在紙類容器強度研究方面還不夠成熟，尤其是折疊紙盒抗壓強度的研究。由于受紙板原材料特性的制約，尚未有用于紙板壓縮應(yīng)力分析評定的系統(tǒng)方法，這樣就不能區(qū)分出應(yīng)力的類型，造成有限元法不能對應(yīng)力進行進一步分析，得出的結(jié)論可能對紙盒結(jié)構(gòu)的優(yōu)化提供不了更加可信的科學依據(jù)，因此，這些方面還有待進一步的研究與發(fā)展。

　　總之，有限元法為包裝設(shè)計及研究人員提供了更為直觀更加快捷的研究手段，通過使用有限元分析的方法我們就能夠建立起紙盒的有限元模型并研究和分析它的力學特性�？梢詷酚^地估計，在不遠的將來，有限元方法將非常廣泛地應(yīng)用到紙盒強度的研究領(lǐng)域，它將推動計算機輔助分析在包裝工程方面的應(yīng)用向更高的水平發(fā)展。

　　參考文獻：

　　[1]段瑞俠，劉雪瑩，陳金周.紙盒結(jié)構(gòu)設(shè)計中要考慮的幾個要素[J].包裝工程，2008(12).

　　[2]余本農(nóng)，平幼妹等.關(guān)于折疊紙盒抗壓強度的討論[J].包裝工程，2002(1).

　　[3]王金林.紙盒抗壓強度與白板紙挺度[J].包裝工程，1989(1).

　　[4]HAKAN Grangard，JOSEF Kubat. Some Aspects of the Compressive Strength of Cartons[J].Svenak Papperstiding，1969(15).

　　[5]劉慧，張新昌.折疊紙盒結(jié)構(gòu)及其尺寸對抗壓強度的影響[J].包裝工程，2008(5).

　　[6]張新昌，劉慧.折疊紙盒抗壓強度經(jīng)驗計算公式的探討[J].包裝工程，2008(10).

　　[7]袁毅，肖穎喆.折疊紙盒鎖底結(jié)構(gòu)承載力實驗研究[J].包裝工程，2014(7).

　　[8]李曉麗.蜂窩紙板靜態(tài)壓縮力學性能仿真研究[D].淮南：安徽理工大學，2011.

　　[9]孫誠.包裝結(jié)構(gòu)設(shè)計[M].第4版.北京：中國輕工業(yè)出版社，2014.

　　[10]張朝暉.ANSYS12.0結(jié)構(gòu)分析工程應(yīng)用實例解析[M].北京：機械工業(yè)出版社，2010.

【有限元法在折疊紙盒強度分析中的應(yīng)用】相關(guān)文章：