蜂窩夾層結(jié)構(gòu)中粘結(jié)劑層的建模與仿真

摘要： 為分析夾層結(jié)構(gòu)制造過程中粘結(jié)劑層對表面質(zhì)量的影響，以聚甲基丙烯酸甲脂（PMMA）為主要研究對象，在有限元分析軟件ABAQUS中進(jìn)行建模分析，通過加熱和加壓2個過程研究在夾層結(jié)構(gòu)的熱壓制造過程中粘結(jié)劑對表面質(zhì)量的影響，結(jié)果表明：當(dāng)僅考慮加熱過程時，由于粘結(jié)劑的粘彈性，隨著溫度的升高，PMMA變軟，導(dǎo)致表面波紋度增大；當(dāng)僅考慮加壓過程時，粘彈性對表面波紋度幾乎無影響；當(dāng)熱壓耦合時，表面波紋度比單獨考慮加熱或加壓的更大；粘結(jié)劑越軟，表面波紋度越大。通過調(diào)整粘結(jié)劑的特性，可以有效控制夾層結(jié)構(gòu)的表面質(zhì)量，從而提高產(chǎn)品的整體質(zhì)量。

關(guān)鍵詞：紙基蜂窩夾層結(jié)構(gòu) 粘結(jié)劑 有限元分析 粘彈性 表面波紋度

中圖分類號：U463.8 " 文獻(xiàn)標(biāo)志碼：B " DOI： 10.19710/J.cnki.1003-8817.20240379

Modeling and Simulation of the Adhesive Layer in Honeycomb Sandwich Structures

Liu Shaopeng1，2， Zhai Liming1， Li Zhiqiang1， Liang Fuxiang1， Rao Cong1， Wang Xin1

（1. FAW Jiefang Engine Division， Wuxi 214000; 2. Technische Universitaet Dresden Institut fuer Festkoepermechanik， Dresden 01069）

Abstract： In order to analyze the influence of the adhesive layer on the surface quality during the manufacturing process of the sandwich structure， modeling and analysis are carried out in the finite element analysis software ABAQUS with PMMA as the main research object. During the hot-pressing manufacturing process of the sandwich structure， the influence of the adhesive on the surface quality is studied through the 2 processes of heating and pressurization. The results show that when only the heating process is considered， due to the viscoelasticity of the adhesive， PMMA will become softer as the temperature rises， which will lead to an increase in the surface waviness; when only the pressurization process is considered， the viscoelasticity has almost no effect on the surface waviness. When the hot-pressing is coupled， the surface waviness will be greater than that when only heating or pressurization is considered alone; the softer the adhesive， the greater the surface waviness. By adjusting the characteristics of the adhesive， the surface quality of the sandwich structure can be effectively controlled， thereby improving the overall quality of the product.

Key words： Paper-based honeycomb sandwich structure， Adhesive， Finite element analysis， Viscoelasticity， Surface waviness

1 前言

汽車制造中越來越多地采用輕量化結(jié)構(gòu)，例如使用復(fù)合材料或紙基材料的結(jié)構(gòu)[1]，在不降低機械性能的同時實現(xiàn)輕量化。

夾層結(jié)構(gòu)是提高輕量化水平的典型策略[2-3]，通過熱壓成形技術(shù)可以制造車身所需的夾層結(jié)構(gòu)。蜂窩芯夾層結(jié)構(gòu)在飛機艙部件中得到了廣泛使用[4]。在高速客車中使用夾層蜂窩結(jié)構(gòu)能減重40%，同時提供降噪和減振等功能[5-6]。典型的夾層結(jié)構(gòu)由多層材料組成，包括上、下兩層和中間芯層。生產(chǎn)夾層結(jié)構(gòu)時，上、下層與中間芯層之間需要使用粘結(jié)劑，使各層間緊密結(jié)合，熱壓后形成一個整體結(jié)構(gòu)。

為了使產(chǎn)品能夠用作車身外殼，需保證紙基夾層結(jié)構(gòu)的表面質(zhì)量。目前，夾層結(jié)構(gòu)表面存在明顯的不平整，可以通過表面波紋度來衡量。為分析和改善夾層結(jié)構(gòu)的表面質(zhì)量，基于有限元分析軟件ABAQUS開發(fā)適用于紙基夾層結(jié)構(gòu)壓制過程的有限元模型。在實際生產(chǎn)中，粘結(jié)劑會對夾層表面的形狀偏差造成影響。因此，把粘結(jié)劑層納入模型中進(jìn)行分析，選擇聚甲基丙烯酸甲酯（Polymethyl Methacrylate， PMMA）作為主要粘結(jié)劑進(jìn)行模擬，考慮粘結(jié)劑的粘彈性對夾層結(jié)構(gòu)表面波紋度的影響，并分析其在熱壓過程中對夾層結(jié)構(gòu)表面質(zhì)量的影響，同時研究具有不同彈性模量的粘結(jié)劑對表面質(zhì)量的影響。

2 理論基礎(chǔ)

粘結(jié)劑的材料行為受溫度和時間的影響，無法通過傳統(tǒng)的彈塑性變形來精確描述。為了通過有限元方法（Finite Element Method，F(xiàn)EM）研究粘結(jié)劑的特性，Hahn [7] 和 Wang [8] 提供了一種描述粘結(jié)劑的策略：粘結(jié)劑的材料行為隨溫度升高表現(xiàn)出強度降低和蠕變增加的特點。因此，與線性彈性相比，粘彈性更能描述粘結(jié)劑的特性。粘彈性指的是材料部分具有彈性，部分具有粘性的特性。粘彈性材料的特性受時間和溫度的影響，在有限元模擬中，將粘結(jié)劑建模為具有溫度依賴性屬性的粘彈性單元，使粘結(jié)劑的建模成為可能[9]。

2.1 ABAQUS中的粘彈性

粘彈性材料的反應(yīng)包括彈性（瞬時）和粘性（與時間、溫度等相關(guān)）行為。材料的瞬時彈性反應(yīng)后，當(dāng)其受到固定應(yīng)力時發(fā)生蠕變，并且在受到固定應(yīng)變時發(fā)生應(yīng)力松弛。在ABAQUS中，材料的變形被分為剪切和體積行為。

線性各向同性粘彈性的基本遺傳積分方程如下[10]：

[σt=0t2Gτ-τ'edt'+Ι0tΚτ-τ'?dt'] （1）

式中：e為剪切變形，Φ為體積變形，G為剪切模量，K為壓縮模量，I為調(diào)整體積變形對應(yīng)力貢獻(xiàn)比例的系數(shù)，t′為積分變量，τ為壓縮時間，τ′與時間延遲或者記憶效應(yīng)有關(guān)。

τ-τ′表明材料的性質(zhì)（剪切模量和壓縮模量）與時間差有關(guān)。通過積分微分方程將壓縮時間與實際時間聯(lián)系起來：

[τ=0tdt'Aθ（θ（t'）） "dτdt=1Aθ（θ（t））] （2）

式中：θ為溫度，Aθ為位移函數(shù)（當(dāng)Aθ=1時，τ=t）。

由此便可將溫度與變形程度聯(lián)系起來，以實現(xiàn)材料的粘彈性。在ABAQUS中通過設(shè)置相對應(yīng)的參數(shù)來實現(xiàn)粘彈性的設(shè)置。

2.2 ABAQUS中使用的網(wǎng)格單元類型

在有限元分析中，網(wǎng)格類型的選擇對于模擬結(jié)果的準(zhǔn)確性、計算效率和經(jīng)濟(jì)性均有重要影響。因此，有必要根據(jù)模型的不同區(qū)域、受力情況和幾何形狀來選擇適當(dāng)?shù)木W(wǎng)格類型。SC8R網(wǎng)格在計算精度和計算時間上表現(xiàn)更好，因此，可作為粘結(jié)劑的網(wǎng)格類型之一。C3D8R網(wǎng)格是ABAQUS中常用的體積網(wǎng)格，它是帶有減少積分的線性六面體元。這2種網(wǎng)格均為三維建模中常用的網(wǎng)格類型，因此本文同時使用這兩種網(wǎng)格進(jìn)行對比，選出較優(yōu)的網(wǎng)格單元。

SC8R網(wǎng)格在厚度方向上沒有明確的層數(shù)規(guī)定。在早期研究中，使用SC8R網(wǎng)格時只采用了1層厚度方向，并取得了良好的模擬結(jié)果。因此，SC8R網(wǎng)格的模型被設(shè)定為只有1層厚度方向。而C3D8R網(wǎng)格中，厚度方向的層數(shù)應(yīng)至少設(shè)置為4，以控制沙漏值。SC8R網(wǎng)格屬于連續(xù)殼單元，僅具有位移自由度，這意味著厚度是根據(jù)元素節(jié)點幾何形狀確定的。正確對齊連續(xù)殼單元非常重要，因為厚度方向和平面方向之間的行為有所不同。在ABAQUS中，通過指定模型的上面和下面來確定厚度方向。如果未設(shè)置厚度方向，則標(biāo)準(zhǔn)方向通常與實際厚度方向不符，可能導(dǎo)致錯誤的計算結(jié)果。由于C3D8R網(wǎng)格本身為體積單元，已包含實際厚度，因此不需要額外的設(shè)置。

2.3 表面質(zhì)量的評價標(biāo)準(zhǔn)

表面質(zhì)量的評價通?；?個指標(biāo)：表面粗糙度和表面波紋度。表面粗糙度通常用來描述表面的微觀結(jié)構(gòu)，比如微小的凹凸、起伏、坑洼等。表面波紋度是指表面的波動或波形，它通常用來描述表面的整體形態(tài)、起伏和輪廓，可以是微觀或宏觀的。表面波紋度的特征可以包括波峰、波谷、波長、波高等。表面粗糙度更側(cè)重于描述表面的微小凹凸和不規(guī)則性，而表面波紋度更關(guān)注表面整體形態(tài)和波動特征。對于車身外殼的表面質(zhì)量，更看重宏觀上的形態(tài)和輪廓，因此本文選擇表面波紋度對表面質(zhì)量進(jìn)行評價。

為了評價模擬結(jié)果的收斂性，需要對表面數(shù)據(jù)進(jìn)行處理，僅通過簡單比較應(yīng)力云圖或特定路徑的坐標(biāo)不能準(zhǔn)確確定表面質(zhì)量。因此，需要使用濾波器對表面的坐標(biāo)數(shù)據(jù)進(jìn)行處理，通過計算經(jīng)過濾波處理的表面波紋度，可以分析結(jié)果的收斂性。

本文使用高斯濾波器對表面數(shù)據(jù)進(jìn)行處理。高斯濾波器是一種低通濾波器。低通濾波器允許低頻信號（即形狀輪廓）通過，同時濾除高頻信號（即表面粗糙度）。λf濾波器用于定義從形狀波紋到更長波長的過渡。λc濾波器用于定義從粗糙度到形狀波紋的過渡。處理后的表面波紋度是通過將λc和λf 濾波器應(yīng)用于原始輪廓而產(chǎn)生的輪廓，其中短波和長波形狀偏差被抑制，留下穩(wěn)定的部分。高斯濾波器所使用的權(quán)重函數(shù)為：

[sx=1α?λf?e（-π（xα?λf）2）] （3）

式中：x為權(quán)重函數(shù)中心（最大值）的距離，λf為從形狀波紋過渡到更長波長的邊界波長，α為提供在邊界波長λf處50%的傳輸特性而設(shè)定的常數(shù)。

使用以下2個特征參數(shù)對過濾后的表面波紋度進(jìn)行評價。

如圖1所示，Wa為表面形狀輪廓線縱坐標(biāo)絕對值的算術(shù)平均值，它表示了輪廓與中心線的平均偏差：

[Wa=1lr0lrZ（x）dx] （4）

如圖2所示，Wz為表面形狀輪廓線高度的平均值。它是沿著單個測量段lr的最大輪廓峰值高度與最大輪廓谷底深度之間的差值的算術(shù)平均值：

[Wz=1ni=1nZci] （5）

3 夾層結(jié)構(gòu)模型搭建

3.1 夾層結(jié)構(gòu)的壓制模型

本文在現(xiàn)有的壓制模型中添加粘結(jié)劑層來進(jìn)行夾層結(jié)構(gòu)的建模。壓制模型基于某公司設(shè)計的夾層蜂窩結(jié)構(gòu)（見圖3）。夾層結(jié)構(gòu)通常由2個外表面層和1個內(nèi)部芯層組成。紙基正弦波狀蜂窩是夾層結(jié)構(gòu)的中間芯層，玻璃纖維與聚氨酯（Polyurethane，PU）樹脂制造的復(fù)合材料為外表面層。此外，最外層通過噴漆涂裝等加工后用作隔離層。

夾層結(jié)構(gòu)的壓制和成形過程已在現(xiàn)有的壓制模型中進(jìn)行了模擬。通過FEM和試驗制造，研究了熱壓成形過程對紙基復(fù)合材料表面質(zhì)量的影響。正弦波蜂窩的形狀在壓制過程中導(dǎo)致夾層結(jié)構(gòu)表面形成凸起且不均勻的波紋。如圖4所示為沒有粘結(jié)劑層的壓制模型。該模型共由5層組成：蜂窩芯層、上模具、下模具、上面板和下面板。蜂窩層的形狀參數(shù)可以根據(jù)需求進(jìn)行調(diào)整，例如振幅、周期和相位角。

3.2 考慮粘結(jié)劑層的壓制模型

在實際制造夾層結(jié)構(gòu)時，上面板不直接與蜂窩芯層接觸，而是通過粘結(jié)劑連接。通過壓合過程，粘結(jié)劑層被壓入蜂窩和上板之間。為了簡化模型并減少仿真時間，可以假設(shè)蜂窩和上板是相互連接固定的。在上面板外層還有一層Senoplast塑料材料。在上面板和Senoplast 之間，也有一層粘結(jié)劑，其影響不可忽視。本文討論該粘結(jié)劑層對表面質(zhì)量的影響。為了有效模擬含粘結(jié)劑層的壓制過程，本文采用了簡化模型、單單元模型和多單元模型。在簡化模型中，只考慮了Senoplast、粘結(jié)劑和上面板，蜂窩芯層被視為剛性體。通過簡化模型，確定了所需的網(wǎng)格單元、接觸條件和邊界條件等，如圖5所示。

確定所需的接觸條件、邊界條件后，通過添加下面板和下模具得到完整的模型。同樣，假設(shè)下面板和蜂窩芯層之間的粘結(jié)劑是固定耦合的。如圖6所示為完整的單單元模型。在單單元模型模擬后，可以將模型從單單元擴展到六單元（即多單元模型，如圖7所示）。一些模擬結(jié)果可以與單單元模型的結(jié)果進(jìn)行比較：加熱、加壓過程，粘結(jié)劑為彈性或者粘彈性，使用不同彈性模量的粘結(jié)劑，考慮熱壓耦合過程。

3.3 邊界條件

材料間的接觸條件是影響計算結(jié)果的關(guān)鍵因素。在確定接觸條件時需確定接觸面的主從性。原則上，應(yīng)首先選擇較大的表面作為主導(dǎo)表面；然后選擇網(wǎng)格較粗糙的表面作為主導(dǎo)表面；如果網(wǎng)格密度大致相等，則應(yīng)選擇硬度較大的材料作為主導(dǎo)表面。詳細(xì)的接觸條件如表1所示，摩擦因數(shù)為0.3，全部為面對面接觸。

在簡化模型中，將蜂窩芯層設(shè)置為剛性體，在相互作用（Interaction）選項中進(jìn)行配置。蜂窩芯層的中心點被選為參考點（Reference Point，RP）。在建模過程中，蜂窩芯層被設(shè)置為剛性體，其RP在加熱或加壓過程中均被設(shè)定約束所有自由度（U1=U2=U3=UR1=UR2=UR3=0）。上模具被設(shè)定為只沿z軸方向移動（U3=0）。在單單元或多單元模型中，蜂窩芯層設(shè)置為可變形體，下模具被固定。施加在剛性體上的約束條件必須作用于RP，在整個分析過程中，剛性體內(nèi)的節(jié)點相對于RP保持不變。

在具有復(fù)雜接觸條件的復(fù)雜模型中，通常通過使用顯式方法來求解材料力學(xué)行為的本構(gòu)方程。在顯式方法中，通常選擇時間離散化的二次方法（如中心差分法）。顯式求解器非常快速，因為它們只需要將時間點t的狀態(tài)向量與一個矩陣相乘，而不需要在每個時間步長Δt中解完整的方程組。因此，它們被用于高動態(tài)和高非線性過程，例如碰撞或成形模擬，以及非常復(fù)雜的模型。并且由于蜂窩結(jié)構(gòu)的幾何形狀和變形復(fù)雜，當(dāng)在ABAQUS中考慮幾何非線性時，將為所有步驟啟用幾何非線性開關(guān)（Nlgeom）設(shè)置。

本文將熱壓過程分解為單獨加熱、單獨加壓或者同時加熱加壓。在加熱過程中，上模具由室溫25 ℃加熱到130 ℃，下模加熱到100 ℃，然后冷卻到室溫。由于溫度升高，材料會發(fā)生熱膨脹。與其他層相比，粘結(jié)劑層受溫度影響更大，因為它不僅會熱膨脹，還會由于溫度導(dǎo)致粘彈性變化。因此，在考慮溫度影響時，重點應(yīng)放在正常彈性粘結(jié)劑和粘彈性粘結(jié)劑之間的差異上。在加壓過程中，通過上模具的下壓和抬升模擬壓制過程。模具的運動速度與試驗中保持一致，為0.5 mm/s。

4 計算結(jié)果分析

4.1 不同網(wǎng)格類型的模擬結(jié)果比較

粘結(jié)劑層分別使用SC8R和C3D8R網(wǎng)格單元進(jìn)行建模。首先需要在簡化模型中確定網(wǎng)格的類型與大小。如圖8所示為在加熱模型中使用不同網(wǎng)格類型和尺寸的模擬結(jié)果。通過表面中某一剖線來展示表面變形，并可由此剖線計算表面波紋度。

隨著網(wǎng)格尺寸的變化，表面波紋度也在變化。由此可以得出，并非越細(xì)致的網(wǎng)格得到的模擬結(jié)果一定越好的，需要后續(xù)與試驗結(jié)果進(jìn)行比對，同時還需考慮計算成本，選擇合適的網(wǎng)格大小進(jìn)行模擬。從試驗結(jié)果和計算成本角度出發(fā)，網(wǎng)格尺寸選擇0.08 mm或0.1 mm。

比較結(jié)果表明，在加熱過程中，使用SC8R網(wǎng)格和C3D8R網(wǎng)格得到的模擬結(jié)果之間存在明顯差異。為了排除厚度方向?qū)訑?shù)對模擬結(jié)果的影響，將SC8R網(wǎng)格的厚度方向?qū)訑?shù)從1增加到5，結(jié)果如表2所示。

結(jié)果表明，厚度方向的層數(shù)幾乎不會對結(jié)果產(chǎn)生影響。這種差異可能是由于SC8R和C3D8R網(wǎng)格單元的本構(gòu)方程在熱膨脹方面的差異造成的。由于SC8R單元基于板理論，忽略了一些因素，這導(dǎo)致了結(jié)果的偏差。其次，可能是因為模型中只有一個單元，模型過小使得結(jié)果并不穩(wěn)定。如果在更大的模型中增加單元格的數(shù)量，處理更多的數(shù)據(jù)，結(jié)果的差異可能會減小。

因此，在比較中加入多單元模型的模擬結(jié)果。表3顯示，在擁有六單元的大模型中，網(wǎng)格單元類型造成的差異明顯減小。相對于只擁有一個單元的模型，更大的模型的穩(wěn)定性更好，但是相對應(yīng)的計算成本將大幅增加。在后續(xù)的模擬中，可以使用簡單模型確認(rèn)邊界條件等參數(shù)，待模擬能夠穩(wěn)定運行時，可以適當(dāng)擴充模型，增加模型穩(wěn)定性。

試驗測得僅受溫度影響的表面波紋度Wa約為0.02～0.05 μm，SC8R網(wǎng)格單元計算得到的結(jié)果可以接受。因此，在后續(xù)的建模中應(yīng)使用殼單元SC8R進(jìn)行建模，在確保準(zhǔn)確率的前提下，計算成本將大幅降低。

在加壓過程中，上模具的移動模擬了壓制過程。對于網(wǎng)格類型與尺寸的影響，壓制過程得出的結(jié)論印證了加熱過程的模擬結(jié)果，在此不再贅述，如圖9所示為加壓過程的模擬結(jié)果。與加熱過程的結(jié)果相比，在加壓過程中，SC8R和C3D8R的差異并未過大（表4）。由此可見，當(dāng)表面變形足夠大（表面波紋度Wa從0.019 μm提升到 5 μm）或模型足夠大時，網(wǎng)格差異對結(jié)果的影響越小。

4.2 粘彈性的影響

通過比較彈性粘結(jié)劑和粘彈性粘結(jié)劑的模擬結(jié)果，分析粘彈性在加熱過程中對表面質(zhì)量的影響。由圖10和表5可以得出以下結(jié)論：在加熱過程中，由于粘彈性導(dǎo)致粘結(jié)劑變軟，表面波紋度變大。這是由于溫度的升高加速了粘結(jié)劑在松弛和蠕變過程中的軟化，從而導(dǎo)致在溫度升高的影響下彈性模量降低，使得夾層結(jié)構(gòu)的表面變形變大。

而在加壓過程中，由于沒有溫度的影響，粘彈性對表面波紋度幾乎無影響，模擬結(jié)果如表6所示。

4.3 不同粘結(jié)劑對表面波紋度的影響

由于在不加熱的壓制過程中，粘彈性幾乎不產(chǎn)生影響，為簡化模型并節(jié)省模擬時間，在壓制模型中不考慮粘結(jié)劑的粘彈性。本文主要針對3種典型的粘結(jié)劑開展研究：PMMA（硬）、PU（軟）和EP（中）。從理論研究的角度來看，其他不同的粘結(jié)劑將直接根據(jù)PMMA的彈性模量的不同百分比近似替代，例如80%、60%、40%、20%和10%的PMMA彈性模量作為不同的粘結(jié)劑。如圖11所示為模擬結(jié)果。

從結(jié)果來看，所用的粘結(jié)劑越軟，表面波紋度越大，這與在實驗室中所得結(jié)論一致。因此，為減小表面波紋度，應(yīng)當(dāng)選擇相對較硬的粘結(jié)劑。

4.4 加熱加壓耦合的模擬結(jié)果

實際生產(chǎn)中，通過熱壓過程生產(chǎn)夾層結(jié)構(gòu)。因此，需要通過熱壓耦合模擬：模具下壓過程中，模具升溫達(dá)到預(yù)定溫度；壓制結(jié)束，模具抬起，溫度冷卻到室溫。結(jié)果如表7所示，熱壓耦合下，表面波紋度明顯加大。

5 結(jié)論

a. 通過完善和擴大模型可以減小網(wǎng)格單元類型對結(jié)果的影響，從計算經(jīng)濟(jì)性和準(zhǔn)確性角度考慮，SC8R網(wǎng)格單元應(yīng)作為首選；在確定網(wǎng)格單元尺寸時，考慮計算成本和準(zhǔn)確性的同時，應(yīng)對比試驗結(jié)果，選擇合適的網(wǎng)格單元尺寸。

b. 在加熱過程中，由于粘結(jié)劑的粘彈性，粘結(jié)劑變軟，表面波紋度增大。與單純的彈性材料相比，粘彈性材料在高溫下的表現(xiàn)需要特別關(guān)注。

c. 在加壓過程中，由于沒有溫度的影響，粘彈性對結(jié)果無明顯影響。這說明單獨加壓過程的表面質(zhì)量主要受物理壓力控制。

d. 在熱壓耦合過程中，相比單過程，表面波紋度明顯增大。在生產(chǎn)中需要注意耦合效應(yīng)對表面質(zhì)量的影響。

e. 當(dāng)使用較軟的粘結(jié)劑時，表面波紋度會增大。因此，選擇硬度較高的粘結(jié)劑有助于提高表面質(zhì)量。

參考文獻(xiàn)：

[1] STAMPB S. Material- und Prozesscharakterisierung zur Herstellung von Sandwichleichtbauteilen mit Papierwabenkern und Thermoplastischen Faserverst?rkten Deckschichten für die Anwendung im Fahrzeuginterieur[J]. Stuttgart： Universit?t Stuttgart， 2013.

[2] RATHBUN H， ZOK F， EVANS A. Strength Optimization of Metallic Sandwich Panels Subject to Bending[J]. International Journal of Solids and Structures， 2005， 42（26）： 6643-6661.

[3] XUE X， ZHANG C， CHEN W， et al. Study on the Impact Resistance of Honeycomb Sandwich Structures Under Low-Velocity/Heavy Mass[J]. Composite Structures， 2019， 226.

[4] HEIMBS S， MIDDENDORF P， MAIER M. Sandwichstrukturen in der Flugzeugkabine： Materialmodellierung und Dynamische Simulation[C]// IVW-Kolloquium. 2006： 45-56.

[5] WENNBERG D. Light-Weighting Methodology in Rail Vehicle Design Through Introduction of Load Carrying Sandwich Panels[D]. Stockholm： KTH Royal Institute of Technology， 2011.

[6] WENNBERG D， STICHEL S. Multi-Functional Design of a Composite High-Speed Train Body Structure[J]. Structural and Multidisciplinary Optimization， 2014， 50（3）： 475-488.

[7] HAHN J， JENDRY J， MAHNKEN R. Simulation der Fügeteilverformung w?hrend der Warmaush?rtung von Klebstoffen： Modellierung und Experimentelle Untersuchungen an einem Schiebedachdeckel[C]// 22-nd CADFEM User Meeting， International Congress on FEM Technology with ANSYS CFX amp; ICEM CFD Conference， 2004.

[8] WANG S， STRATFORD T，REYNOLDS T. Linear Creep of Bonded FRP-Strengthened Metallic Structures at Warm Service Temperatures[J]. Construction and Building Materials， 2021， 283.

[9] ROYLANCE D. Engineering Viscoelasticity. Massachusetts Institute of Technology[M]. Cambridge， MA： ， 2001： 1-37.

[10] SMITH M. ABAQUS/Standard User’s Manual[DB]. United States： Dassault Systèmes Simulia Corp.， 2020.