蜂窩胞壁厚度與邊長對等效模型計(jì)算精度的影響

康洪軍，江蘭馨，王明猛

(1.中車唐山機(jī)車車輛有限公司 產(chǎn)品研發(fā)中心，河北 唐山 063000;2.西南交通大學(xué) 牽引動(dòng)力國家重點(diǎn)實(shí)驗(yàn)室，四川 成都 610031)

0 引言

近年來,隨著車輛輕量化技術(shù)的發(fā)展,鋁蜂窩板結(jié)構(gòu)以其特輕的質(zhì)量、良好的強(qiáng)度與剛度以及降噪等特性在國外軌道車輛尤其是高速列車(如磁浮列車)上開始得到應(yīng)用,并顯示出廣闊前景,在國內(nèi)也受到了較大的關(guān)注。

鋁蜂窩夾層板結(jié)構(gòu)多樣，可通過設(shè)計(jì)其面板厚度、芯層高度、蜂窩胞壁厚度和蜂窩胞體邊長等參數(shù)并進(jìn)行實(shí)驗(yàn)以提高其力學(xué)性能。當(dāng)前的工程軟件如ANSYS、NASTRAN中沒有專門的蜂窩結(jié)構(gòu)單元，若使用蜂窩的精確模型仿真車體，則工作量和計(jì)算量過大。所以，在進(jìn)行鋁蜂窩車體仿真計(jì)算時(shí)只能采用三維實(shí)體模型或等效板模型[1]。其中三維實(shí)體模型計(jì)算精度較高，但是計(jì)算量大；等效板模型計(jì)算量小，應(yīng)用更為廣泛。當(dāng)前蜂窩夾層板的等效板模型有Reissner理論模型[2-3]、蜂窩板理論模型和層合板低階剪切理論模型。分析以上等效理論可發(fā)現(xiàn)，正六面體蜂窩結(jié)構(gòu)的胞壁厚度與邊長的比值對等效模量、等效密度和等效泊松比均有影響。以上理論研究并沒有考慮蜂窩胞壁厚度與邊長比值對不同等效板模型計(jì)算精度的影響，但在實(shí)際應(yīng)用中胞壁厚度與邊長之比對等效板模型計(jì)算精度的影響規(guī)律是不一致的。研究蜂窩夾層板長厚比對等效板模型計(jì)算精度的影響規(guī)律，對于含蜂窩夾層板用于車體計(jì)算時(shí)等效模型的選擇方面有著重要的影響。

本文首先歸納了芯層均勻壁厚正六邊形鋁蜂窩夾層板的等效板模型，包括Reissner理論模型、蜂窩板理論模型和層合板低階剪切理論模型；然后針對自由邊界下的不同胞壁厚度的蜂窩夾層板，將等效板模型的模態(tài)分析結(jié)果與精細(xì)化模型進(jìn)行對比，分析胞壁厚度與邊長之比值比對模態(tài)頻率計(jì)算精度的影響；最后，總結(jié)了胞壁厚度與邊長之比對等效板模型動(dòng)力學(xué)計(jì)算精度的影響規(guī)律，為含蜂窩夾層板構(gòu)件的計(jì)算提供參考。

1 正六邊形鋁蜂窩夾層板的等效模型

在鋁蜂窩夾層板結(jié)構(gòu)中，蜂窩與面板之間膠層厚度小、質(zhì)量輕且不是主要承載件，故在進(jìn)行蜂窩精確建模時(shí)忽略膠層的影響。本文主要研究芯層均勻壁厚的正六邊形鋁蜂窩夾層板的簡化板模型。圖1為蜂窩夾層板的幾何參數(shù)：hf為上、下面板厚度；hc為芯層厚度；H為總厚度；t為蜂窩胞壁厚度；l為蜂窩胞體邊長。

圖1 蜂窩夾層板簡易模型

本文研究鋁蜂窩胞壁厚度與邊長之比對蜂窩夾層板等效模型計(jì)算精度的影響時(shí)，保持蜂窩胞體邊長不變，通過改變胞壁厚度來探究兩者之比對等效模型計(jì)算精度的影響。胞壁厚度與邊長之比為λ=t/l。圖2為蜂窩夾層板的精細(xì)化有限元模型。

圖2 蜂窩夾層板精細(xì)化有限元模型

蜂窩夾層板中，面板和蜂窩芯材均采用牌號為6A02的鋁材。材料6A02常規(guī)力學(xué)性能參數(shù)如表1所示，蜂窩夾層板仿真結(jié)構(gòu)組合類型如表2所示。

表1 材料6A02常規(guī)力學(xué)性能參數(shù)

表2 蜂窩夾層板結(jié)構(gòu)參數(shù)仿真組合類型 單位：mm

1.1 芯材均質(zhì)化

芯層均質(zhì)化是將精細(xì)化模型的蜂窩芯層等效為等厚度實(shí)體芯層。假設(shè)芯層結(jié)構(gòu)極軟，僅僅能夠抵抗橫向切應(yīng)力，忽略其彎曲剛度和面內(nèi)剛度；因其上、下面板極薄，假設(shè)上、下面板服從Kirchhoff 假設(shè)，只能承受面內(nèi)應(yīng)力，忽略其抵抗橫向切應(yīng)力的能力[4]。根據(jù)Gibson提出的胞元材料理論，在進(jìn)行分析時(shí)將其簡化為線彈性Bernoulli-Euler梁，則可將蜂窩芯層等效為正交各向異性的均質(zhì)化材料，材料參數(shù)計(jì)算公式如下：

(1)

式中：Es、Gs、μs、ρs分別為蜂窩芯層結(jié)構(gòu)材料的彈性模量、切變模量、泊松比和密度；Ecx、Ecy、Ecz分別為蜂窩芯層均質(zhì)化后等效的3個(gè)主方向的彈性模量；Gcxy、Gcxz、Gcyz為蜂窩芯層均質(zhì)化后等效的面內(nèi)切變模量；μcxy、μcxz、μcyz為蜂窩芯層均質(zhì)化后等效的泊松比；ρc為蜂窩芯層均質(zhì)化后等效的芯層密度。

1.2 Reissner理論模型

Reissner理論等效模型是將芯層均質(zhì)化后的蜂窩夾層板等效為與原模型不等厚的單層板。假設(shè)表層處于薄膜應(yīng)力狀態(tài)，忽略芯層平行于xy面內(nèi)應(yīng)力分布，由于在芯層和表面中z方向的應(yīng)力極小，故忽略z方向的應(yīng)力[4]。根據(jù)以上假設(shè)，Reissner理論等效模型的基本參數(shù)計(jì)算公式如下：

(2)

式中：tr、ρr、Er、μr分別為Reissner理論等效模型等效后的厚度、密度、彈性模量和泊松比；Ef、μf、ρf分別為蜂窩夾層板面板材料的彈性模量、泊松比和密度。

1.3 蜂窩板理論模型

蜂窩板理論將整個(gè)蜂窩夾層板等效成等剛度、同尺寸、等厚度的正交各向異性板。利用哈密爾頓原理可得出蜂窩夾層結(jié)構(gòu)的等效參數(shù)[5]，根據(jù)以上原理，蜂窩板理論模型的基本參數(shù)計(jì)算公式如下：

(3)

式中：Ex、Ey、Gxy、Gxz、Gyz、μxy、ρr分別為蜂窩板理論模型等效后的彈性模量、切變模量、泊松比和密度。其中中間變量計(jì)算公式如下：

其中

1.4 層合板低階剪切理論模型

層合板低階剪切理論是將芯層均質(zhì)化后的蜂窩夾層板等效為與原板相同尺寸且層間位移低階連續(xù)的3層剪切板，等效板的第1層和第3層為各向同性，中間層為正交各向異性[6]。上、下面板參數(shù)為已知鋁板參數(shù)、芯材參數(shù)為均質(zhì)化后芯層參數(shù)，利用ANSYS中復(fù)合層合板理論對91號單元進(jìn)行建模，低層合板理論模型的基本參數(shù)計(jì)算公式如下：

(4)

式中：E′、μ′、ρ′分別為鋪層建模時(shí)第1層和第3層材料的彈性模量、泊松比和密度；E″x、E″y、E″z為鋪層建模時(shí)中間層材料的3個(gè)方向的彈性模量；G″xy、G″xz、G″yz為中間層材料的面內(nèi)切變模量；μ″xy、μ″xz、μ″yz為中間層材料的3個(gè)主方向的泊松比；ρ?為中間層材料的等效密度。

2 蜂窩結(jié)構(gòu)胞壁厚度與邊長之比對等效模型的影響分析

運(yùn)用有限元計(jì)算軟件ANSYS對精細(xì)化模型和相應(yīng)等效模型進(jìn)行自由邊界下的模態(tài)計(jì)算，將等效模型的模態(tài)計(jì)算結(jié)果與精細(xì)化模型計(jì)算結(jié)果做誤差分析。各模型建模過程中，面板和芯層厚度、胞體邊長不變，改變蜂窩胞壁厚度，分析λ對等效模型計(jì)算精度的影響。

進(jìn)行精細(xì)化模型建模時(shí)，蜂窩夾芯層與上、下面板均采用殼單元建模。進(jìn)行芯層均質(zhì)化后的理論模型建模時(shí)，芯層采用實(shí)體單元(Solid185)建模，上、下面板采用殼(Shell181)單元建模。Reissner理論模型采用均質(zhì)板模，蜂窩板理論模型采用lamina板建模，層合板低階剪切理論模型采用復(fù)合材料層合板模塊建模。仿真計(jì)算時(shí)采用的蜂窩夾層板幾何尺寸為500 mm×500 mm，蜂窩胞體邊長l=6 mm,蜂窩胞壁厚度t=0.1 mm、0.15 mm、0.2 mm、0.25 mm、0.3 mm、0.35 mm。

通過式(1)計(jì)算不同胞壁厚度下芯層均質(zhì)化的材料參數(shù)，并且建立不同胞壁厚度下與精細(xì)模型對應(yīng)的均質(zhì)化模型。在自由邊界下，計(jì)算精細(xì)化模型與均質(zhì)化模型的前4階模態(tài)。表3為精細(xì)化模型自由邊界下前4階模態(tài)頻率值，圖3為t=0.2 mm的精細(xì)化模型的陣型圖。分析不同λ下，均質(zhì)化等效模型與精細(xì)化模型的誤差值，結(jié)果如表4所示。

表3 精細(xì)化模型自由邊界下的模態(tài)頻率

圖3 t=0.2 mm的精細(xì)化模型陣型圖

表4 不同t下均質(zhì)化模型的模態(tài)頻率

通過式(2)計(jì)算不同胞壁厚度下Reissner理論模型的材料參數(shù)，在自由邊界下，計(jì)算Reissner理論模型的前4階模態(tài)頻率，并將其結(jié)果與精細(xì)化模型的模態(tài)結(jié)果相對比，分析不同λ下Reissner等效模型與精細(xì)化模型的誤差值，結(jié)果如表5所示。

表5 不同t下Reissner模型的模態(tài)頻率

通過式(3)計(jì)算不同胞壁厚度下蜂窩板理論等效模型的材料參數(shù)。在自由邊界下，計(jì)算蜂窩板理論等效模型的前4階模態(tài)，并將其結(jié)果與精細(xì)化模型的模態(tài)結(jié)果相對比，分析不同λ下蜂窩板等效模型與精細(xì)化模型的誤差值，結(jié)果如表6所示。

表6 不同t下蜂窩板模型的模態(tài)頻率

通過式(4)計(jì)算不同胞壁厚度下層合板低階剪切理論模型的材料參數(shù)。在自由邊界下，計(jì)算層合板低階剪切理論模型的前4階模態(tài)，并將其結(jié)果與精細(xì)化模型的模態(tài)結(jié)果相對比，分析不同λ下層合板等效模型與精細(xì)化模型的誤差值，結(jié)果如表7所示。

表7 不同t下層合板低階剪切模型的模態(tài)頻率

將胞壁厚度不等的各等效模型的計(jì)算結(jié)果與相應(yīng)精細(xì)化模型進(jìn)行對比，將其誤差絕對值繪制成折線圖，如圖4所示。

圖4 各等效模型誤差絕對值

由圖4可以觀察得到蜂窩胞壁的厚度變化對均質(zhì)化等效模型的模態(tài)計(jì)算精度基本無影響，前4階的計(jì)算誤差絕對值維持在1%以下，誤差主要來源為模型簡化時(shí)計(jì)算量的簡化。胞壁厚度的變化對層合板低階剪切理論模型的模態(tài)計(jì)算精度影響不大，前4階的計(jì)算誤差絕對值維持在3%～6% 的范圍內(nèi)。在t≤0.2 mm時(shí)，層合板低階剪切理論模型的計(jì)算優(yōu)于蜂窩板理論等效模型；在t≤0.15 mm時(shí)，層合板低階剪切理論模型的計(jì)算優(yōu)于Reissner和蜂窩板理論等效模型。Reissner理論等效模型的模態(tài)計(jì)算精度隨著胞壁厚度的增大而減小，并隨著胞壁厚度的增加計(jì)算精度趨于穩(wěn)定，在t=0.1 mm時(shí)，Reissner模型的精度最低，當(dāng)t≥0.15 mm時(shí)，Reissner理論等效模型有較好的計(jì)算精度。同樣，蜂窩板理論等效模型的模態(tài)計(jì)算精度也隨著胞壁厚度的增加而減小，并隨著胞壁厚度的增加計(jì)算精度趨于穩(wěn)定，在t≥0.15 mm時(shí)，其計(jì)算精度維持在5% 左右，在胞壁厚度在0.2 mm≤t≤0.35 mm時(shí)，Reissner等效模型精度為板殼類模型中最高的，精度誤差<2%。

3 結(jié)語

對蜂窩不同胞壁厚度與邊長之比的各個(gè)理論等效模型與精細(xì)化模型模態(tài)計(jì)算結(jié)果比較可得出以下結(jié)論：

1)不同胞壁厚度與邊長的比值下，均質(zhì)化理論等效模型模態(tài)計(jì)算結(jié)果精度最高，與精細(xì)化模型對比誤差<1%。

2)Reissner理論等效模型和蜂窩板理論等效模型模態(tài)的計(jì)算精度隨著胞壁厚度與邊長比例的增大而減小，在λ≥0.025 mm時(shí)，兩種理論模型前4階模態(tài)的誤差<10%。蜂窩板理論等效模型的精度低于Reissner理論等效模型。

3)層合板低階剪切理論模型的模態(tài)計(jì)算精度隨胞壁厚度與邊長的比值變化基本不變，與精細(xì)化模型的計(jì)算誤差維持在5%左右。