曹敏娜，楊國(guó)超，張求慧

紙木蜂窩板面外靜態(tài)壓縮有限元仿真分析

曹敏娜a，楊國(guó)超b，張求慧b

（北京林業(yè)大學(xué) a.木質(zhì)材料科學(xué)與應(yīng)用教育部重點(diǎn)實(shí)驗(yàn)室 b.木材科學(xué)與工程北京重點(diǎn)實(shí)驗(yàn)室，北京 100083）

探究紙木蜂窩板的壓潰模式和單個(gè)蜂窩單元的變形失效模式，為紙木蜂窩板內(nèi)部失效機(jī)制提供一定的理論基礎(chǔ)。通過(guò)有限元仿真模擬，分析面外靜態(tài)壓縮得到變形云圖、應(yīng)力–應(yīng)變?cè)茍D和應(yīng)力–應(yīng)變曲線(xiàn)。紙木蜂窩板的受壓過(guò)程分為4個(gè)階段，應(yīng)力集中分布在蜂窩芯層，最外側(cè)的蜂窩最早出現(xiàn)變形且變形較大，隨后中間區(qū)域的蜂窩逐步變形增大。蜂窩芯層在受壓時(shí)存在至多1個(gè)應(yīng)力峰，蜂窩壁彎曲折疊了1～2次，產(chǎn)生了1～2個(gè)塑性鉸，并在此過(guò)程中吸收能量。紙木蜂窩板的壓潰失效趨勢(shì)為由外側(cè)蜂窩向內(nèi)部蜂窩，由單層壁向雙層壁，由最弱的壓潰點(diǎn)向其余位置。

紙木蜂窩板；蜂窩芯層；面外靜態(tài)壓縮；有限元

蜂窩夾層板是一種環(huán)境友好的仿生材料，多孔的蜂窩結(jié)構(gòu)使其具備比強(qiáng)度高、質(zhì)量小、抗沖擊性高、緩沖隔振好等優(yōu)點(diǎn)，目前已廣泛應(yīng)用在緩沖包裝、運(yùn)輸和家具領(lǐng)域[1-3]。在運(yùn)輸傳遞過(guò)程中，蜂窩板會(huì)受到靜載荷、沖擊振動(dòng)等作用力，有限元仿真技術(shù)已經(jīng)應(yīng)用到蜂窩板平面承載性能和緩沖振動(dòng)的模擬分析中[4-6]。

目前，國(guó)內(nèi)外學(xué)者對(duì)蜂窩板的研究集中在蜂窩芯層和面板為同質(zhì)材料的分析，尤其集中在對(duì)蜂窩紙板和金屬蜂窩板的力學(xué)性能研究和有限元分析，對(duì)關(guān)于采用2種材料的復(fù)合蜂窩板的研究較為少見(jiàn)[7-8]。將紙蜂窩和木板2種基材制備成紙木蜂窩板，其力學(xué)性能和抗彎性能要優(yōu)于傳統(tǒng)蜂窩紙板的，強(qiáng)重比又高于普通板材的。已有學(xué)者對(duì)紙木蜂窩板進(jìn)行了力學(xué)性能研究，但有關(guān)有限元仿真的分析還未見(jiàn)報(bào)道。在紙木蜂窩板的各項(xiàng)力學(xué)性能中，面外靜態(tài)壓縮機(jī)理復(fù)雜，涉及到蜂窩芯層的局部失穩(wěn)。由于測(cè)試手段的限制，試驗(yàn)只能觀察到最外側(cè)的部分蜂窩壁的變形失效，無(wú)法測(cè)得蜂窩芯層內(nèi)部的應(yīng)力、應(yīng)變變化和單個(gè)蜂窩單元的失效模式?；诖耍疚囊哉呅渭埛涓C作為芯層，以三層楊木膠合板作為結(jié)構(gòu)面板，提出一種紙木復(fù)合的蜂窩板結(jié)構(gòu)。應(yīng)用試驗(yàn)和有限元分析相結(jié)合的方式分析紙木蜂窩板的面外靜態(tài)壓縮過(guò)程，探究紙木蜂窩板的壓潰模式和蜂窩單元的變形失效模式，為紙木蜂窩板的平面力學(xué)性能和后續(xù)應(yīng)用提供理論基礎(chǔ)。

1 紙木蜂窩板仿真模型及材料參數(shù)

1.1 幾何結(jié)構(gòu)

紙木蜂窩板的正六邊形蜂窩芯層由瓦楞紙沿縱向施膠黏接后拉伸制成，故蜂窩單元在厚度上會(huì)分為單層壁厚和雙層壁厚。規(guī)定蜂窩芯層拉伸成型前的長(zhǎng)度方向?yàn)榭v方向，簡(jiǎn)稱(chēng)，對(duì)應(yīng)空間坐標(biāo)方向；蜂窩芯層的拉伸方向?yàn)闄M方向，簡(jiǎn)稱(chēng)，對(duì)應(yīng)空間坐標(biāo)方向；蜂窩芯層的厚度方向?yàn)?，?duì)應(yīng)空間坐標(biāo)方向，如圖1所示。除最外層蜂窩單元外壁外，所有垂直于方向的蜂窩單元壁均為雙層壁厚。在Solidworks中對(duì)紙蜂窩芯層和膠合板面板分別建模，采用自下而上的方式進(jìn)行裝配體組合，導(dǎo)入Ansys Workbench進(jìn)行有限元仿真分析。由于瓦楞紙芯紙?jiān)诤穸确较蛏系某叽邕h(yuǎn)遠(yuǎn)小于其他2個(gè)方向，在壓縮過(guò)程中會(huì)產(chǎn)生失穩(wěn)導(dǎo)致大變形。為確保仿真精度，蜂窩芯層采用殼單元Shell181建模更為合理，膠合板面板采用實(shí)體單元Solid186建模。根據(jù)GB/T 1453—2022《夾層結(jié)構(gòu)或芯子平壓性能試驗(yàn)方法》[9]中的規(guī)定，蜂窩芯層尺寸為60 mm×60 mm，厚度為15 mm。在切割制樣過(guò)程中為保證紙蜂窩尺寸符合國(guó)家標(biāo)準(zhǔn)，蜂窩芯層在2個(gè)側(cè)面會(huì)產(chǎn)生非封閉的蜂窩單元，在受力后會(huì)更容易出現(xiàn)局部應(yīng)力集中現(xiàn)象[10]。為更好地等效紙木蜂窩板受力情況，將兩側(cè)未封閉的蜂窩單元進(jìn)行還原建模，即紙木蜂窩板的4個(gè)側(cè)面形狀各不相同。將紙蜂窩芯層全封閉的2個(gè)側(cè)面記作Ⅰ，未封閉的2個(gè)側(cè)面記作Ⅱ，紙木蜂窩板的模型如圖2所示。

1.2 材料參數(shù)

紙木蜂窩板的面外靜態(tài)壓縮過(guò)程涉及材料非線(xiàn)性、幾何非線(xiàn)性和狀態(tài)非線(xiàn)性，是一個(gè)復(fù)雜的非線(xiàn)性響應(yīng)過(guò)程，還存在壓縮后的大變形和接觸問(wèn)題[11]。采用雙線(xiàn)性正交各向異性彈塑性材料模型定義瓦楞原紙和膠合板的材料本構(gòu)關(guān)系。紙蜂窩芯層由定量為95.2 g/m2的高強(qiáng)瓦楞紙制備，瓦楞原紙的厚度根據(jù)GB/T 451.3—2002《紙和紙板的厚度測(cè)定》[12]進(jìn)行測(cè)定。由于瓦楞紙?jiān)埖牟牧咸匦?，壓縮試驗(yàn)難度極大，根據(jù)GB/T 12914—2018《紙和紙板抗張強(qiáng)度的測(cè)定恒速拉伸法（20 mm/min）》[13]，通過(guò)萬(wàn)能力學(xué)試驗(yàn)機(jī)對(duì)瓦楞紙進(jìn)行拉伸試驗(yàn)，獲得其力學(xué)參數(shù)。參考Mann等[14]、Baum等[15]和Persson等[16]的研究成果，根據(jù)經(jīng)驗(yàn)公式（1）獲得其他參數(shù)。面外靜態(tài)壓縮方向?yàn)橥呃慵埖目v向，對(duì)應(yīng)材料參數(shù)中的方向，試驗(yàn)過(guò)程中瓦楞紙?jiān)诖朔较虬l(fā)生屈曲變形。瓦楞紙屬于屈服現(xiàn)象不明顯的材料，規(guī)定屈服強(qiáng)度為材料產(chǎn)生0.2%的應(yīng)變時(shí)的應(yīng)力值，測(cè)得試驗(yàn)數(shù)值為2.28 MPa。瓦楞紙?jiān)诤穸确较蛏系某叽邕h(yuǎn)遠(yuǎn)小于其他2個(gè)方向的。泊松比μ取0.3，μ和μ在計(jì)算中通常取0.01。按照GB/T 17657—2013《人造板及飾面人造板理化性能試驗(yàn)方法》[17]對(duì)膠合板的密度進(jìn)行測(cè)定，通過(guò)萬(wàn)能力學(xué)試驗(yàn)機(jī)進(jìn)行抗壓力學(xué)試驗(yàn)。根據(jù)GB/T 1453—2022《夾層結(jié)構(gòu)或芯子平壓性能試驗(yàn)方法》[9]對(duì)紙木蜂窩板進(jìn)行面外靜態(tài)壓縮試驗(yàn)。最后獲得的材料參數(shù)如表1所示。

圖1 紙蜂窩芯層結(jié)構(gòu)

圖2 紙木蜂窩板模型

表1 材料參數(shù)

Tab.1 Material parameters

1.3 分析設(shè)定

將紙木蜂窩板上下面板和蜂窩芯層分別建立部件，部件表示具有相同單元類(lèi)型的單元集合，可實(shí)現(xiàn)共享網(wǎng)格、節(jié)點(diǎn)耦合，每個(gè)部件可以單獨(dú)配置材料屬性，自成一體。紙木蜂窩板的上下面板和蜂窩芯層是黏接而成的，在試驗(yàn)過(guò)程中未發(fā)生脫膠現(xiàn)象，在有限元分析的接觸定義中選擇綁定約束，以便在仿真模擬過(guò)程中實(shí)現(xiàn)整體效應(yīng)，不發(fā)生分離。部件內(nèi)部已默認(rèn)接觸無(wú)須再重復(fù)設(shè)置，只需對(duì)部件間設(shè)置綁定接觸，即將蜂窩芯層上表面的所有邊與上面板的下表面進(jìn)行綁定，將蜂窩芯層下表面的所有邊與下面板的上表面進(jìn)行綁定。網(wǎng)格劃分精度對(duì)仿真模擬結(jié)果會(huì)產(chǎn)生較大影響，高質(zhì)量的網(wǎng)格有更好的邊界網(wǎng)格形態(tài)，能夠獲取更精確的解。對(duì)膠合板面板采用Hex六面體網(wǎng)格劃分，對(duì)紙蜂窩芯層采用Qaud四邊形網(wǎng)格劃分。劃分后單元質(zhì)量高達(dá)0.999 6，共計(jì)82 903個(gè)節(jié)點(diǎn)、20 460個(gè)單元。

對(duì)紙木蜂窩板施加方向的位移載荷，在求解結(jié)果中分別提取支反力與方向定向變形值隨時(shí)間變化的曲線(xiàn)，轉(zhuǎn)化整合成仿真的應(yīng)力–應(yīng)變曲線(xiàn)。為了更好地對(duì)比試驗(yàn)與仿真結(jié)果，本文挑選一組試驗(yàn)的應(yīng)力–應(yīng)變曲線(xiàn)與仿真曲線(xiàn)進(jìn)行對(duì)比，結(jié)果見(jiàn)圖3。試驗(yàn)獲得的平壓強(qiáng)度均值為0.153 MPa，有限元模擬值為0.143 MPa，數(shù)值非常接近，如表2所示。仿真曲線(xiàn)相比試驗(yàn)曲線(xiàn)更平滑，尤其是在平臺(tái)階段，這主要是由于實(shí)際試驗(yàn)時(shí)蜂窩芯層在坍塌堆疊后可能會(huì)與周邊蜂窩接觸，蜂窩間存在支撐作用，讓已經(jīng)到達(dá)平臺(tái)階段的蜂窩，仍有可能產(chǎn)生應(yīng)力的二次上升?？偟膩?lái)看，仿真結(jié)果與試驗(yàn)結(jié)果呈現(xiàn)高度一致性，驗(yàn)證了有限元仿真的準(zhǔn)確性。

2 結(jié)果與分析

2.1 整體變形

紙木蜂窩板整體壓縮變形過(guò)程和應(yīng)力云圖如圖4所示。從應(yīng)力分布可以清晰看到，應(yīng)力集中分布在蜂窩芯層，面板受到的應(yīng)力遠(yuǎn)遠(yuǎn)小于紙蜂窩芯層的應(yīng)力。變形失效破壞的位置發(fā)生在紙蜂窩芯層部位，面板的變形值很小，這也與試驗(yàn)結(jié)果相符合。在一定范圍內(nèi)，紙木蜂窩板的承載能力主要由芯層決定，面板的強(qiáng)度對(duì)蜂窩板性能影響不大[18]。

圖3 紙木蜂窩板試驗(yàn)與仿真的應(yīng)力–應(yīng)變曲線(xiàn)

表2 紙木蜂窩板試驗(yàn)與仿真的平壓強(qiáng)度

Tab.2 Test and simulation flatwise compression properties of paper-wood honeycomb board

紙木蜂窩板的靜態(tài)壓縮過(guò)程可以分成4個(gè)階段。第1階段是彈性階段，這個(gè)階段分為2個(gè)部分：在壓縮初期，蜂窩壁在軸向發(fā)生彈性變形，此時(shí)應(yīng)力隨應(yīng)變線(xiàn)性增加；若繼續(xù)加載超過(guò)彈性極限后，就變成非線(xiàn)性的彈性變形，直至達(dá)到蜂窩芯層的臨界最大屈曲載荷，試件發(fā)生破壞。第2階段是彈塑性屈曲階段，應(yīng)力隨應(yīng)變的增加而減少，此時(shí)由彈性變形導(dǎo)致的局部坍塌轉(zhuǎn)變?yōu)橐运苄糟q為特征的屈曲；屈曲首先發(fā)生在單層壁上，隨后發(fā)生在雙層壁上。第3階段是塑性坍塌階段，應(yīng)變持續(xù)增加，應(yīng)力基本保持不變，呈現(xiàn)出近似平臺(tái)的區(qū)域，平臺(tái)的寬度代表了蜂窩板緩沖吸能的能力。此階段表現(xiàn)為以蜂窩壁塑性坍塌為主的持續(xù)壓潰過(guò)程，壓潰載荷隨應(yīng)變的增加呈現(xiàn)一定的上下波動(dòng)，直至完全被壓塌。第4階段是密實(shí)化階段，蜂窩壁充分坍塌變形發(fā)生接觸，直至蜂窩芯層完全被壓實(shí)，此時(shí)微小應(yīng)變將引起應(yīng)力急劇增大[19]。紙木蜂窩板的壓縮過(guò)程以蜂窩芯層單層壁的失穩(wěn)為開(kāi)始，以雙層壁達(dá)到最大承載力發(fā)生破壞為結(jié)束。

為更清楚地了解紙蜂窩芯層的壓潰過(guò)程，將上下面板進(jìn)行隱藏，觀察蜂窩芯層在持續(xù)載荷作用下的應(yīng)力云圖。如圖5所示，蜂窩芯層的變形不同步，呈現(xiàn)一定的先后順序。從總體來(lái)看，應(yīng)力集中最早出現(xiàn)在最外側(cè)的蜂窩單層壁的中心區(qū)域，而后逐漸向中間區(qū)域的蜂窩單元擴(kuò)展，外側(cè)相對(duì)自由的蜂窩單層壁應(yīng)力較大。有限元模擬的變形過(guò)程與試驗(yàn)圖相對(duì)應(yīng)，可見(jiàn)試驗(yàn)過(guò)程中蜂窩芯層最外側(cè)首先發(fā)生變形，尤其是4個(gè)角落的處的蜂窩變形明顯。這是由于最外側(cè)存在未封閉的單層壁蜂窩，其抗壓能力要弱于全封閉的蜂窩，也弱于雙層壁的蜂窩，因此首先出現(xiàn)應(yīng)力集中和大變形。隨著位移載荷的增大，中間區(qū)域的單層壁蜂窩也發(fā)生較大的變形，主要載荷由中間區(qū)域的雙層壁蜂窩承擔(dān)，出現(xiàn)大應(yīng)力，最終全部蜂窩都被屈服壓潰[20]。對(duì)應(yīng)試驗(yàn)圖同樣可以看到變形逐步由最外側(cè)向中心處的蜂窩移動(dòng)，隨著位移載荷的持續(xù)增大，中間區(qū)域蜂窩變形值也逐步增加。對(duì)每個(gè)獨(dú)立的蜂窩單元而言，應(yīng)力集中最早出現(xiàn)在蜂窩的上下節(jié)點(diǎn)連線(xiàn)的中心區(qū)域，隨著位移載荷的增加先沿著豎直方向由中心區(qū)域向上下節(jié)點(diǎn)擴(kuò)展；然后在水平方向上由各節(jié)點(diǎn)和節(jié)點(diǎn)連線(xiàn)向周?chē)鷶U(kuò)展，直至全部被壓潰變形。對(duì)于單個(gè)獨(dú)立的蜂窩的變形，實(shí)際試驗(yàn)不能達(dá)到有限元模擬的理想狀態(tài)，其屈服并非一定發(fā)生在蜂窩上下節(jié)點(diǎn)連線(xiàn)的中心區(qū)域，也有可能在蜂窩芯層最上端或最下端首先發(fā)生屈服變形。試驗(yàn)過(guò)程中屈服壓潰順序取決于對(duì)應(yīng)區(qū)域壓潰點(diǎn)的強(qiáng)弱，最弱的壓潰點(diǎn)所在區(qū)域首先被壓潰。

2.2 單個(gè)蜂窩變形

為進(jìn)一步探究受壓時(shí)處于不同位置的單個(gè)蜂窩單元在直立壁高度方向上的變形情況，分別在紙蜂窩芯層的中心蜂窩上取一條平行軸的路徑，4個(gè)側(cè)面Ⅰ、Ⅱ、Ⅰ、Ⅱ的蜂窩上同樣各取一條平行軸的同名路徑。路徑方向由節(jié)點(diǎn)1指向節(jié)點(diǎn)2，節(jié)點(diǎn)1代表起始位置0，分析5條異面路徑的等效應(yīng)力和應(yīng)變的分布情況。其中等效總應(yīng)變是包括彈性應(yīng)變和塑性應(yīng)變?cè)趦?nèi)的應(yīng)變總和。

如圖6所示，單個(gè)蜂窩單元的受壓變形主要表現(xiàn)為在近似節(jié)點(diǎn)1和節(jié)點(diǎn)2的連線(xiàn)中心處對(duì)稱(chēng)面發(fā)生彎曲、局部的變形折疊。根據(jù)Wierzbicki[21]的研究表明，蜂窩壁每彎曲折疊一次都會(huì)產(chǎn)生一個(gè)塑性鉸，按照波長(zhǎng)進(jìn)行漸進(jìn)折疊，而波長(zhǎng)約等于蜂窩單元邊長(zhǎng)。因此蜂窩單元高度越大，邊長(zhǎng)越小，產(chǎn)生的塑性鉸會(huì)越多，能量吸收就越多。本研究中紙蜂窩芯層的蜂窩單元邊長(zhǎng)為6 mm，故在壓縮過(guò)程中將按照=6 mm進(jìn)行漸進(jìn)折疊，蜂窩單元高度為15 mm，蜂窩芯層將折疊1～2次，如圖6可見(jiàn)蜂窩壁產(chǎn)生了1～2個(gè)塑性鉸。外側(cè)的封閉蜂窩單元折疊彎曲對(duì)稱(chēng)度更高，發(fā)生折疊的區(qū)域相對(duì)小而集中；外側(cè)的未封閉蜂窩單元彎曲扭曲度更高，發(fā)生折疊的區(qū)域更大，并呈現(xiàn)出一定的非對(duì)稱(chēng)性，這可能與較弱的承載能力所導(dǎo)致的應(yīng)力分布不均有關(guān)；處在中心區(qū)域的蜂窩單元雖承受較大應(yīng)力，但整體變形較小，未見(jiàn)其在高度方向發(fā)生彎曲折疊。外側(cè)的蜂窩壁比中心處的蜂窩壁有更大的變形，首先發(fā)生彎曲失效，這也與試驗(yàn)結(jié)果相匹配。

圖4 紙木蜂窩板靜態(tài)壓縮過(guò)程

圖5 紙蜂窩靜態(tài)壓縮過(guò)程

圖7為5條路徑由節(jié)點(diǎn)1至節(jié)點(diǎn)2的應(yīng)力和應(yīng)變曲線(xiàn)，處在封閉蜂窩單元的路徑Ⅰ和Ⅰ分別在靠近上下面板6 mm左右的位置，即在距離為6 mm和9 mm處產(chǎn)生了應(yīng)力峰值，這與上文提及的塑性鉸產(chǎn)生的位置對(duì)應(yīng)。應(yīng)力峰處對(duì)應(yīng)的位置代表可能壓潰的點(diǎn)，雖然在7.5 mm處，即1/2處的應(yīng)力較小，但由于塑性鉸的存在，在此產(chǎn)生了最大變形和應(yīng)變。未封閉的蜂窩單元Ⅱ和Ⅱ在靠近上下面板的位置和1/2處均產(chǎn)生了應(yīng)力峰，這些均為可能的壓潰點(diǎn)。壓潰點(diǎn)并非同時(shí)被壓潰，往往是較弱的首先被壓潰，然后逐步向其他的壓潰點(diǎn)擴(kuò)展[15]。由應(yīng)變曲線(xiàn)可見(jiàn)，未封閉的蜂窩單元也在約1/2處壓潰點(diǎn)發(fā)生了大變形，可見(jiàn)1/2處的壓潰點(diǎn)相對(duì)較弱。處在中心處的蜂窩單元的壓潰點(diǎn)更多，應(yīng)力–應(yīng)變分布更為均勻。處在外側(cè)的封閉的蜂窩單元可承受的應(yīng)力峰值最大，約為外側(cè)未封閉的蜂窩單元可承受應(yīng)力的2倍，中心的蜂窩單元應(yīng)力承受能力次之。

圖6 單個(gè)蜂窩的應(yīng)力應(yīng)變?cè)茍D

圖7 單個(gè)蜂窩的應(yīng)力–應(yīng)變曲線(xiàn)

3 結(jié)語(yǔ)

通過(guò)試驗(yàn)測(cè)試和有限元仿真模擬，研究了紙木蜂窩板在面外靜態(tài)壓縮下的整體壓潰機(jī)制和單個(gè)蜂窩單元的變形失效模式。紙木蜂窩板的受壓過(guò)程分為4個(gè)階段，應(yīng)力集中分布在蜂窩芯層，面板受到的應(yīng)力很小，承載能力主要由蜂窩芯層決定。蜂窩芯層的變形最早出現(xiàn)在最外側(cè)的蜂窩，隨著位移載荷的增大中間區(qū)域的蜂窩變形逐步增大，直至全部蜂窩屈服壓潰。處在外側(cè)的蜂窩相比于中心的蜂窩會(huì)發(fā)生更大的變形，蜂窩在受壓時(shí)會(huì)按照波長(zhǎng)約等于蜂窩單元邊長(zhǎng)進(jìn)行漸進(jìn)折疊，每折疊一次產(chǎn)生一個(gè)塑性鉸，在此過(guò)程中吸收能量。壓縮過(guò)程中蜂窩壁存在1到多個(gè)應(yīng)力峰，這些應(yīng)力峰代表著可能的壓潰點(diǎn)，較弱的首先被壓潰，文中多為1/2處的壓潰點(diǎn)首先被壓潰，然后逐步向其他的壓潰點(diǎn)擴(kuò)展。綜上，紙木蜂窩板的壓潰失效規(guī)律總的趨勢(shì)為由外側(cè)向內(nèi)部，由單層壁向雙層壁，由最弱的壓潰點(diǎn)向其余位置。本文探究了蜂窩芯層內(nèi)部的應(yīng)力–應(yīng)變和單個(gè)蜂窩的失效模式，為紙木蜂窩板的平壓性能探究和后續(xù)應(yīng)用提供一定的理論和技術(shù)支持。

[1] OZEN I, AVA K, GEDIKLI H, et al. Low-energy Impact Response of Composite Sandwich Panels with Thermoplastic Honeycomb and Reen-Trant Cores[J]. Thin-Walled Structures, 2020, 156: 106-989.

[2] ZHANG Xiao-ming, XIE Juan, CHEN Jin-xiang, et al. The Beetle Elytron Plate: A Lightweight, High-Strength and Buffering Functional-Structural Bionic Material[J]. Scientific Reports, 2017, 7(1): 4440.

[3] 王強(qiáng), 王寧, 王鈴偉, 等. 輕質(zhì)蜂窩形仿生材料在家具設(shè)計(jì)中的應(yīng)用與發(fā)展研究[J]. 家具與室內(nèi)裝飾, 2020(10): 82-84.

WANG Qiang, WANG Ning, WANG Ling-wei, et al. Study on Application and Development of Lightweight Honeycomb Bionic Materials in Furniture Design[J]. Furniture & Interior Design, 2020(10): 82-84.

[4] 李娜, 尤青文. 蜂窩夾層復(fù)合材料小質(zhì)量沖擊接觸力的分析[J]. 科學(xué)技術(shù)與工程, 2020, 20(14): 5558-5565.

LI Na, YOU Qing-wen. Analysis on Contact Force of Small-Mass-Impact on Honeycomb Sandwich Composites[J]. Science Technology and Engineering, 2020, 20(14): 5558-5565.

[5] 楊震, 李賜生. 實(shí)木家具包裝箱的設(shè)計(jì)與探討[J]. 家具與室內(nèi)裝飾, 2018(12): 110-111.

YANG Zhen, LI Ci-sheng. Design and Discussion of Solid Wood Furniture Packaging Cases[J]. Furniture & Interior Design, 2018(12): 110-111.

[6] 楊穩(wěn), 張勝蘭, 李瑩. 蜂窩夾層結(jié)構(gòu)等效模型研究進(jìn)展[J]. 復(fù)合材料科學(xué)與工程, 2020(10): 122-128.

YANG Wen, ZHANG Sheng-lan, LI Ying. Research Progress of Equivalent Model of Honeycomb Sandwich Structure[J]. Composites Science and Engineering, 2020(10): 122-128.

[7] 王宏磊. 蜂窩夾芯復(fù)合材料的力學(xué)性能研究[D]. 吉林: 吉林大學(xué), 2019: 1-3.

Wang Hong-lei. Mechanical Properties of Honeycomb Sandwich Composites[D]. Jilin: Jilin University, 2019: 1-3.

[8] 楊雨恒, 張周鎖, 史文博, 等. 鋁蜂窩夾層板動(dòng)力學(xué)建模仿真與試驗(yàn)驗(yàn)證[J]. 機(jī)械設(shè)計(jì), 2021, 38(8): 9-16.

YANG Yu-heng, ZHANG Zhou-suo, SHI Wen-bo, et al. Dynamic Modeling Simulation and Experimental Verification of the Aluminum Honeycomb Sandwich Panel[J]. Journal of Machine Design, 2021, 38(8): 9-16.

[9] GB/T 1453—2022, 夾層結(jié)構(gòu)或芯子平壓性能試驗(yàn)方法[S].

GB/T 1453—2022, Test Method for Flatwise Compression Properties of Sandwich Constructions or Cores[S].

[10] 潘鳳麗, 李光. 基于Ansys Workbench的塑料蜂窩板仿真分析及優(yōu)化[J]. 包裝工程, 2016, 37(23): 44-49.

PAN Feng-li, LI Guang. Ansys Workbench-Based Simulation Analysis and Optimization of Plastic Perforated Board[J]. Packaging Engineering, 2016, 37(23): 44-49.

[11] 湯金若, 朱瑾, 于江祥. 曲面蜂窩板的彎曲特性仿真分析研究[J]. 現(xiàn)代防御技術(shù), 2016, 44(5): 186-193.

TANG Jin-ruo, ZHU Jin, YU Jiang-xiang. Simulations on Buckling Characteristics of Curving Honeycomb[J]. Modern Defense Technology, 2016, 44(5): 186-193.

[12] GB/T 451.3—2002, 紙和紙板的厚度測(cè)定[S].

GB/T 451.3—2002, Paper and Board-Determination of Thickness[S].

[13] GB/T 12914—2018, 紙和紙板抗張強(qiáng)度的測(cè)定恒速拉伸法(20 mm/min)[S].

GB/T 12914—2008, Paper and Board-Determination of Tensile Properties-Constant Rate of Elongation Method(20 mm/min)[S].

[14] MANN R W, BAUM G A, HABEGER C C, et al. Determination of All Nine Orthotropic Elastic Constants for Machine-made Paper[J]. Appleton Wisconsin the Institute, 1980, 63(2): 163-166.

[15] BAUM G A, BRENNAN D C, HABEGER C C. Orthotropic Elastic-constants of Paper[J]. Tappi Journal, 1981, 64(8): 97-101.

[16] PERSSON K. Material Model for Paper: Experimental and Theoretical Aspects[D]. Lund: Lund University, 1991.

[17] GB/T 17657—2013, 人造板及飾面人造板理化性能試驗(yàn)方法[S].

GB/T 17657—2013, Test Methods of Evaluating the Properties of Wood-based Panels and Surface Decorated Wood-based Panels[S].

[18] 李厚民. 蜂窩紙板受壓時(shí)的有限元分析[J]. 包裝工程, 2006, 27(1): 34-36.

LI Hou-min. FEM Analysis of the Compression Properties of Honeycomb Paperboard[J]. Packaging Engineering, 2006, 27(1): 34-36.

[19] 陳岳鈴, 王軍, 盧立新, 等. 凹六邊形蜂窩紙板面內(nèi)承載性能有限元仿真分析[J]. 包裝工程, 2021, 42(13): 174-179.

CHEN Yue-ling, WANG Jun, LU Li-xin, et al. Finite Element Simulation Analysis of In-Plane Load-Bearing Performance of the Paperboard with re-Entrant Hexagonal Honeycomb Core[J]. Packaging Engineering, 2021, 42(13): 174-179.

[20] 李楊, 賁小軍, 袁建東, 等. 熱塑性蜂窩板的平壓性能分析[J]. 包裝學(xué)報(bào), 2022, 14(2): 93-97.

LI Yang, BEN Xiao-jun, YUAN Jian-dong, et al. Analysis of Flat Compression Performance of Thermoplastic Honeycomb Panel[J]. Packaging Journal, 2022, 14(2): 93-97.

[21] WIERZBICKI T. Crushing Analysis of Metal Honeycombs[J]. International Journal of Impact Engineering, 1983, 1(2): 157-174.

Finite Element Simulation Analysis on Out of Plane Static Compression of Paper-wood Honeycomb Board

CAO Min-naa, YANG Guo-chaob, ZHANG Qiu-huib

(a. Key Laboratory of Wood Material Science and Application, Ministry of Education b. Beijing Key Laboratory of Wood Science and Engineering, Beijing Forestry University, Beijing 100083, China)

The work aims to explore the collapse mode of paper-wood honeycomb boards and the deformation failure mode of a single honeycomb unit, so as to provide a theoretical basis for the internal failure mechanism of paper-wood honeycomb boards. Through finite element simulation, the out of plane static compression was analyzed and its deformation pattern, stress-strain pattern and stress-strain curve were obtained. The results showed that the compression process of paper-wood honeycomb boards was divided into four stages. The stress was concentrated in the honeycomb core layer. The outer most honeycomb deformed greatly first, and then the deformation of the middle honeycomb gradually increased. There was one or more stress peaks in the honeycomb core layer under compression. The honeycomb wall bent and folded for 1-2 times, generating 1-2 plastic hinges and absorbing energy in this process. In a word, the collapse failure trend of paper-wood honeycomb boards is from the outer honeycomb to the inner honeycomb, from the single wall to the double wall, and from the weakest collapse point to the rest.

paper-wood honeycomb board; honeycomb core layer; out of plane static compression; finite element

TS653.8、O242.21

A

1001-3563(2023)09-0199-07

10.19554/j.cnki.1001-3563.2023.09.024

2022?10?24

曹敏娜（1998—），女，碩士生，主攻紙木蜂窩板有限元仿真分析。

張求慧（1960—），女，博士，教授，研究方向主要為家具材料與包裝材料。

責(zé)任編輯：曾鈺嬋