Free vibration characteristics of corrugated sandwich plates under different boundary conditions
-
摘要: 波纹夹芯板作为一种特殊的复合材料结构，边界条件对其振动特性有重要影响。根据不同剪切方式下的剪切变形理论和基尔霍夫经典板理论(CLPT)，利用Hamilton原理建立波纹夹芯板的动力学方程。其中，波纹芯层等效成各向异性均质体。根据四边简支、四边固支、对边简支和固支、一边固支三边简支的边界条件，推导出位移形式的偏微分动力学方程。求解得到波纹夹芯板在不同边界条件下自由振动的固有频率，与有限元仿真结果进行对比，验证了理论结果的正确性。在此基础上，基于指数剪切变形理论(ESDT)，分析了不同边界条件下波纹夹芯板的基频随材料参数和结构几何参数的变化规律。结果表明，材料和几何参数对不同边界条件下波纹夹芯板的振动特性有重要影响。相关研究结果将对波纹夹芯板在工程应用中的减振设计及优化分析提供一定的理论依据。Abstract: As a special composite structure, the vibration characteristics of corrugated sandwich panel are greatly influenced by the boundary conditions. According to the shear deformation theory of different shear modes and Kirchhoff's classical plate theory(CLPT), the dynamic equation of corrugated sandwich plates was established by Hamilton principle. Among them, the corrugated core layer was equivalent to an anisotropic homogeneous body. According to the boundary conditions of four sides simply supported, four sides clamped, opposite sides simply supported and clamped, one side fixed and three edges clamped, the partial differential dynamic equation relative to the displacements was derived. By solving the equation, the natural frequencies of the corrugated sandwich plates under different boundary conditions were obtained. Compared with the finite element simulation results, the correctness of the theoretical results was verified. On this basis, based on the exponential shear deformation theory(ESDT), the variation of fundamental frequency of the corrugated sandwich plate with material parameters and structural geometric parameters under different boundary conditions was analyzed. The results show that the material and structural geometric parameters have an important influence on the vibration characteristics of the corrugated sandwich plates under different boundary conditions. Relevant research results will provide a theoretical basis for the vibration reduction design and optimization analysis of corrugated sandwich plates in engineering application.
-
近年来由于航空航天、高铁的飞速发展，航空航天器、高速列车等在高速运行时会产生剧烈的振动，该结构的减振设计将成为技术的关键。波纹夹芯结构作为一种新型的复合材料组合结构，具有高比强度、高比刚度、质量轻、吸能能力强、承载效率高等良好特性，在航空航天、高速列车、汽车、建筑、机械、船舶等领域有着广泛的应用[1-4]。与传统的复合结构材料相比，波纹夹芯板结构很好地克服了夹芯板不连续的问题，具有良好的承载性能。此外，波纹夹芯板芯层间隙具备一定的流通性能，在间隙中穿插泡沫、多孔纤维、玻璃纤维等材料还能起到一定的隔声和隔热作用。因此，对波纹夹芯板的研究具有重要的意义。
由于夹层板拥有广阔的应用前景，近年来国内外学者针对波纹夹芯板各方面特性展开了广泛的研究。Peng等[5] 提出一种无网格伽辽金法来研究非加筋和加筋波纹夹芯板的弹性弯曲问题，将波纹夹芯板视为在两个垂直方向上具有不同挠曲力的正交各项异性板，并推导了梯形、正弦型波纹夹芯板的等效参数。Semenyuk等[6]研究了三种用于纵向波纹圆柱壳稳定性分析的设计方法。Briassoulis等[7]基于正交各向异性壳的拉伸刚度和弯曲刚度的一组解析表达式，对波纹壳进行了数值模拟。Sohrab等[8]提出一种新型的分层波纹复合芯，并通过实验进行了测试。Tian等[9]对波纹夹芯板进行了优化分析。Samanta等[10]首次提出了梯形波纹夹芯板的非线性几何分析，将波纹夹芯板等效成各向异性模型，对梯形波纹夹芯板进行了自由振动分析。Xia等[11]提出了一种基于均质化的波纹夹芯板模型，该模型可用于任何波纹形状。Li 等[12]针对铝制波纹夹芯板开展了空中爆炸试验研究和有限元数值分析，验证了有限元分析技术的可行性。王红霞等[13]和王青伟等[14]推导了考虑波纹拉伸变形的三角形波纹夹芯板的等效弹性参数，后来进行了修正。吴建均等[15]推导了泡沫填充下的三角形波纹夹芯梁的等效弹性常数公式，并且对泡沫波纹夹芯梁的模态进行实验研究和数值计算。还有一些学者对波纹夹芯板的结构性能、力学性能和冲击性能进行了研究[16-18]。
在夹芯板或板结构的研究中，发展了不同的理论，如Reissner板理论、基尔霍夫经典板理论(CLPT)及后来发展的各种剪切变形理论[19-21]等。Talha等[22]基于高阶剪切理论研究了梯度板的静态响应和自由振动。WU等[23]提出了基于RMVT的三阶剪切变形理论(TSDT)。THAI等[24]用正弦剪切变形理论(SSDT)分析了功能梯度板的弯曲、屈曲和振动特性。曹源等[25]用SSDT和修正的欧应力理论研究了功能梯度三明治微梁的静态弯曲和自由振动。但是，用这些理论来研究波纹夹芯板性能的还比较少。
综上所述，波纹夹芯板具有很强的各向异性特性，目前多数文献研究四边简支边界条件下波纹夹芯板的动力学特性，但波纹夹芯板在实际应用中还有其它的边界条件，并且边界条件的影响较显著。鉴于此，本文将采用不同的夹层板理论研究波纹夹芯板在四边简支(SSSS)、四边固支(CCCC)、对边简支和固支(CCSS)、一边固支三边简支(CSSS)四种边界条件下的自由振动，与ABAQUS有限元仿真结果进行对比，验证理论模型的正确性，分析边界条件对波纹夹芯板振动特性的影响。此外，基于指数剪切变形理论（ESDT），讨论波纹夹芯板的材料参数和结构几何参数对系统振动特性的影响。所得结果为轻质波纹夹芯结构的优化设计提供了必要的理论指导。
1. 波纹夹芯板模型
波纹夹芯板的模型如图1(a)所示，由上、下面板及中间的波纹芯层组成。在板中面建立x-y坐标系，
z 轴垂直于x-y面，z>0 的一侧称为下表面，z<0 的一侧称为上表面。波纹夹芯板的长度为a，宽度为b，高度为h，上、下板厚为hf 。图1(b)为波纹芯层的胞元，单胞的底边长为2p ，波纹壁厚为tc ，波纹与面板的夹角为θ，斜边长lc=p/cosθ 芯层厚度hc= ptanθ 。在波纹夹芯板中，上、下面板采用各向同性均质体，中间波纹芯层等效为各向异性均质体，等效示意图如图2所示。仿照文献[13-15]，考虑波纹夹芯板的伸缩变形，波纹芯层的等效参数表达式为
E1(2)=EstcpsinθE2(2)=Estc3cosθhc3[1+(tchc)2cos2θ]G12(2)=Gsptcsinθhc2G13(2)=GstcsinθpG23(2)=Estcsin2θcosθhcv12(2)=vstc2cos2θ(hc2+tc2cos2θ)v21(2)=vsρ(2)=2ρstclcsin2θ (1) 其中：
E1(2) 、E2(2) 、G12(2) 、G13(2) 、G23(2) 、v12(2) 、v21(2) 、ρ(2) 分别表示波纹芯层在各方向上的等效弹性模量、剪切模量、泊松比和密度；Es、Gs和ρs分别为基体材料的弹性模量、剪切模量和密度。![]()
图 1 波纹夹芯板的模型(a)和波纹单胞示意图(b)Figure 1. Model diagrams of corrugated sandwich panel (a) and corrugated cell (b)a—Length of corrugated sandwich panel; b—Width of corrugated sandwich panel; lc—Length of the hypotenuse; hc—Height of core layer; tc—Wall thickness; θ—Corrugation angle; p—Length of bottom side2. 波纹夹芯板振动理论模型
考虑到波纹夹芯板的横向剪切变形效应，位移场可表示为如下形式：
u(x,y,z;t)=u0−z∂w0∂x+f(z)ϕxv(x,y,z;t)=v0−z∂w0∂y+f(z)ϕyw(x,y,z;t)=w0 (2) 其中：
u0 、v0 、w0 分别为波纹夹芯板中面上任意一点的位移；ϕx 和ϕy 分别为波纹夹芯板的直法线沿x轴和y轴的转角；f(z) 是决定横向剪切应力和应变沿厚度分布的形状函数。表1为几种不同剪切形状函数[19, 23-26]。CLPT不考虑横向剪切变形，低估了挠度，高估了固有频率和屈曲载荷，只适用于薄板，而一阶剪切变形理论(FSDT)在此基础上考虑了横向剪切变形，但是在上下表面处不满足面力自由的条件，需要剪切修正因子；ESDT、SSDT和TSDT都考虑了横向剪切变形，且在夹层板的上下表面处满足面力自由的条件，在计算中不需要横向剪切修正因子，很好地克服了CLPT和FSDT的局限性，三者的横向剪切应力和应变沿厚度的形状分别为指数函数型、正弦函数型和曲线型。表 1 不同板理论对应的剪切形状函数Table 1. Shear shape functions corresponding to different plate theoriesShear theory Function f (z) CLPT f(z)=0 FSDT f(z)=z SSDT f(z)=hπsin(πzh) TSDT f(z)=z[1−43(zh)2] ESDT f(z)=ze−2(zh)2 Notes: CLPT—Classical plate theory; FSDT—First-order shear plate theory; SSDT—Sinusoidal shear deformation theory; TSDT—Third-order shear deformation theory; ESDT—Exponential shear deformation theory. 根据小变形假设，位移-应变关系有：
εxx=∂u∂x,εyy=∂v∂y,εzz=∂w∂z,γxz=∂u∂z+∂w∂x,γxy=∂u∂y+∂v∂x,γyz=∂v∂z+∂w∂y (3) 将式(2)代入式(3)得到应变分量为
{εxxεyyγxy}={εxx(0)εyy(0)γxy(0)}+z{εxx(1)εyy(1)γxy(1)}+f{εxx(3)εyy(3)γxy(3)},{γyzγxz}=f′{γyz(2)γxz(2)},εzz=0 (4) 其中，各应变分量的具体表达式为
εxx(0)=∂u0∂x,εxx(1)=−∂2w0∂x2,εxx(3)=∂ϕx∂x,εyy(0)=∂v0∂y,εyy(1)=−∂2w0∂y2,εyy(3)=∂ϕy∂y,γxy(0)=∂u0∂y+∂v0∂x,γxy(1)=−2∂2w0∂x∂y,γxy(3)=∂ϕx∂y+∂ϕy∂x,γyz(2)=ϕy,γxz(2)=ϕx (5) 波纹夹芯板的本构关系可表示为
{σ(k)xxσ(k)yyτ(k)yzτ(k)xzτ(k)xy}=[Q(k)11Q(k)12000Q(k)21Q(k)2200000Q(k)4400000Q(k)5500000Q(k)66]{εxxεyyγyzγxzγxy} (6) 其中，
σ(k)xx 、σ(k)yy 和τ(k)yz 、τ(k)xz 、τ(k)xy 为波纹夹芯板的应力和剪应力。各刚度系数可以表示为
\begin{split} &Q_{11}^{(k)} = \frac{{E_1^{(k)}}}{{1 - \nu_{12}^{(k)}\nu_{21}^{(k)}}},Q_{12}^{(k)} = \frac{{E_1^{(k)}\nu_{21}^{(k)}}}{{1 - \nu_{12}^{(k)}ν_{21}^{(k)}}},Q_{22}^{(k)} = \frac{{E_2^{(k)}}}{{1 - \nu_{12}^{(k)}\nu_{21}^{(k)}}}, \\ &{Q_{66}} = {G_{12}},{Q_{44}} = {G_{23}},{Q_{55}} = {G_{13}},{Q_{21}} = {Q_{12}} \\[-12pt] \end{split} (7) 其中：k=1、2、3 分别代表波纹夹芯板的上、中、下层；E1(k)、E2(k)、G12(k)、G13(k)、G23(k)、ν12(k)、ν21(k)分别表示上、下面板和芯层的弹性模量、剪切模量和泊松比。
波纹夹芯板系统的势能、动能及外力势可以表示为
\begin{split} {\rm{\delta}}U = &\int_0^a {\int_0^b {\int_{ - \frac{h}{2}}^{\frac{h}{2}} {({\sigma _{xx}}{\rm{\delta}}{\varepsilon _{xx}}} } } + {\sigma _{yy}}{\rm{\delta}}{\varepsilon _{yy}} + {\sigma _{{\textit{z}}{\textit{z}}}}{\rm{\delta}}{\varepsilon _{{\textit{z}}{\textit{z}}}} + {\tau _{y{\textit{z}}}}{\rm{\delta}}{\gamma _{y{\textit{z}}}} + \\ &{\tau _{x{\textit{z}}}}{\rm{\delta}}{\gamma _{x{\textit{z}}}} + {\tau _{xy}}{\rm{\delta}}{\gamma _{xy}}){\rm{d}}x{\rm{d}}y{\rm{d}}{\textit{z}} \\[-12pt] \end{split} (8) {\rm{\delta}}K = \int_0^a {\int_0^b {\int_{ - \frac{h}{2}}^{\frac{h}{2}} {\rho (\dot u{\rm{\delta}}\dot u + \dot v{\rm{\delta}}\dot v + \dot w{\rm{\delta}}\dot w)} } } {\rm{d}}x{\rm{d}}y{\rm{d}}{\textit{z}} (9) {\rm{\delta}}V = \int_0^a {\int_0^b {\int_{ - \frac{h}{2}}^{\frac{h}{2}} {q{\rm{\delta}}w{\rm{d}}x{\rm{d}}y{\rm{d}}{\textit{z}}} } } (10) 其中：
\dot u 、\dot v 、\dot w 为波纹夹芯板在x 、y 、{\textit{z}} 方向的运动速度；q为外载荷。根据哈密顿变分原理[27]，得到波纹夹芯板系统的动力学方程为
\begin{split} &\frac{{\partial {N_{xx}}}}{{\partial x}} + \frac{{\partial {N_{xy}}}}{{\partial y}} = {I_0}{{\ddot u}_0} - {I_1}\frac{{\partial {{\ddot w}_0}}}{{\partial x}}{\rm{ + }}{J_1}{{\ddot \phi }_x} \\ &\frac{{\partial {N_{yy}}}}{{\partial y}} + \frac{{\partial {N_{xy}}}}{{\partial x}} = {I_0}{{\ddot v}_0} - {I_1}\frac{{\partial {{\ddot w}_0}}}{{\partial y}}{\rm{ + }}{J_1}{{\ddot \phi }_y} \\ &\frac{{{\partial ^2}{M_{xx}}}}{{\partial {x^2}}} + 2\frac{{{\partial ^2}{M_{xy}}}}{{\partial x\partial y}} + \frac{{{\partial ^2}{M_{yy}}}}{{\partial {y^2}}} + q = {I_0}{{\ddot w}_0}{\rm{ + }}{I_1}\frac{{\partial {{\ddot u}_0}}}{{\partial x}}- \\ &{I_2}\frac{{{\partial ^2}{{\ddot w}_0}}}{{\partial {x^2}}} + {K_2}\frac{{\partial {{\ddot \phi }_x}}}{{\partial x}} + {I_1}\frac{{\partial {{\ddot v}_0}}}{{\partial y}} - {I_2}\frac{{{\partial ^2}{{\ddot w}_0}}}{{\partial {y^2}}} + {K_2}\frac{{\partial {{\ddot \phi }_y}}}{{\partial y}} \\ &\frac{{\partial {H_{xx}}}}{{\partial x}} + \frac{{\partial {H_{xy}}}}{{\partial y}} - {Q_x} = {J_1}{{\ddot u}_0} - {K_2}\frac{{\partial {{\ddot w}_0}}}{{\partial x}} + {J_2}{{\ddot \phi }_x} \\ &\frac{{\partial {H_{yy}}}}{{\partial y}} + \frac{{\partial {H_{xy}}}}{{\partial x}} - {Q_y} = {J_1}{{\ddot v}_0} - {K_2}\frac{{\partial {{\ddot w}_0}}}{{\partial y}} + {J_2}{{\ddot \phi }_y} \\ \end{split} (11) 其中：
\begin{array}{*{20}{l}} {\left[ {\begin{array}{*{20}{c}} {{N_{\xi \eta }}}\\ {{M_{\xi \eta }}}\\ {{H_{\xi \eta }}} \end{array}} \right] = \displaystyle\int_{ - \frac{h}{2}}^{\frac{h}{2}} {{\sigma _{\xi \eta }}} \left[ {\begin{array}{*{20}{c}} 1\\ {\textit{z}}\\ f \end{array}} \right]{\rm{d}}{\textit{z}},{Q_\xi } = \displaystyle\int_{ - \frac{h}{2}}^{\frac{h}{2}} {{\sigma _{\xi {\textit{z}}}}} f'{\rm{d}}{\textit{z}}}\\ {[{I_0},{I_1},{I_2},{J_1},{J_2},{K_2}] = \displaystyle\sum\limits_{k = 1}^3 {\int_{{\zeta _k}}^{{\zeta _{k + 1}}} {{\rho ^{(k)}}[1,{\textit{z}},{{\textit{z}}^2},f,{f^2},{\textit{z}}f]} } {\rm{d}}{\textit{z}}} \end{array} (12) 其中，
\xi 、\eta 可以用x或y表示。将式(4)~(6)、式(12)代入式(11)得到位移形式的运动控制方程为
{E_{i1}}{u_0} + {E_{i2}}{v_0} + {E_{i3}}{w_0} + {E_{i4}}{\phi _x} + {E_{i5}}{\phi _y} = 0,\; i = 1\sim5 (13) 其中：
{E_{11}} = {A_{11}}\frac{{{\partial ^2}}}{{\partial {x^2}}} + {A_{66}}\frac{{{\partial ^2}}}{{\partial {y^2}}} - {I_0}\frac{{{\partial ^2}}}{{\partial {t^2}}} {E_{22}} = {A_{22}}\frac{{{\partial ^2}}}{{\partial {y^2}}} + {A_{66}}\frac{{{\partial ^2}}}{{\partial {x^2}}} - {I_0}\frac{{{\partial ^2}}}{{\partial {t^2}}} {E_{12}}{\rm{ = }}{E_{21}} = ({A_{21}} + {A_{66}})\frac{{{\partial ^2}}}{{\partial x\partial y}} \begin{split} {E_{33}} =& - {G_{11}}\frac{{{\partial ^4}}}{{{x^4}}} - ({G_{12}} + {G_{21}} + 4{G_{66}})\frac{{{\partial ^4}}}{{\partial {x^2}\partial {y^2}}} - \\ &{G_{22}}\frac{{{\partial ^4}}}{{{y^4}}} + {I_2}{\nabla ^2}(\frac{{{\partial ^2}}}{{\partial {t^2}}}) - {I_0}\frac{{{\partial ^2}}}{{\partial {t^2}}} \\ \end{split} {E_{44}} = {R_{11}}\frac{{{\partial ^2}}}{{\partial {x^2}}} + {R_{66}}\frac{{{\partial ^2}}}{{\partial {y^2}}} + {V_{55}} - {J_2}\frac{{{\partial ^2}}}{{{t^2}}} {E_{55}} = {R_{22}}\frac{{{\partial ^2}}}{{\partial {y^2}}} + {R_{66}}\frac{{{\partial ^2}}}{{\partial {x^2}}} + {V_{44}} - {J_2}\frac{{{\partial ^2}}}{{\partial {t^2}}} {E_{13}} = {E_{31}} = {B_{11}}\frac{{{\partial ^3}}}{{{x^3}}} + ({B_{21}} + 2{B_{66}})\frac{{{\partial ^3}}}{{\partial x\partial {y^2}}} - {I_1}\frac{{{\partial ^3}}}{{\partial x\partial {t^2}}} {E_{14}} = {E_{41}} = {D_{11}}\frac{{{\partial ^2}}}{{\partial {x^2}}} + {D_{66}}\frac{{{\partial ^2}}}{{\partial {y^2}}} - {J_1}\frac{{{\partial ^2}}}{{\partial {t^2}}} {E_{15}} = {E_{51}} = {E_{24}} = {E_{42}} = ({D_{21}} + {D_{66}})\frac{{{\partial ^2}}}{{\partial x\partial y}} {E_{23}} = {E_{32}} = - {B_{22}}\frac{{{\partial ^3}}}{{\partial {y^3}}} - ({B_{21}} + 2{B_{66}})\frac{{{\partial ^2}}}{{\partial {x^2}\partial y}} + {I_1}\frac{{{\partial ^3}}}{{\partial y\partial {t^2}}} {E_{25}} = {E_{52}} = {D_{22}}\frac{{{\partial ^2}}}{{\partial {y^2}}} + {D_{66}}\frac{{{\partial ^2}}}{{\partial {x^2}}} - {J_1}\frac{{{\partial ^2}}}{{\partial {t^2}}} {E_{34}} = {E_{43}} = {L_{11}}\frac{{{\partial ^3}}}{{{x^3}}} + ({L_{21}} + 2{B_{66}})\frac{{{\partial ^3}}}{{\partial x\partial {y^2}}} - {K_2}\frac{{{\partial ^3}}}{{\partial x\partial {t^2}}} {E_{35}} = {E_{53}} = ({L_{12}} + 2{L_{66}})\frac{{{\partial ^3}}}{{\partial {x^2}\partial y}} + {L_{22}}\frac{{{\partial ^3}}}{{{y^3}}} - {K_2}\frac{{{\partial ^3}}}{{\partial y\partial {t^2}}} {E_{45}} = {E_{54}} = ({R_{11}} + {R_{66}})\frac{{{\partial ^2}}}{{\partial x\partial y}} \begin{split} &[{A_{ij}},{B_{ij}},{D_{ij}},{G_{ij}},{L_{ij}},{R_{ij}},{V_{ij}}] \\ &= \sum\limits_{k = 1}^3 {\int_{{\zeta _k}}^{{\zeta _{k + 1}}} {{Q^{(k)}}_{ij}[1,{\textit{z}},f,{{\textit{z}}^2},{\textit{z}}f} ,} {f^2},{(f')^2})]{\rm{d}}{\textit{z}} \end{split} 假设位移分量为双三角函数，则
{u_0} 、{v_0} 、{w_0} 、{\phi _x} 、{\phi _y} 可表示为\begin{split} &\left\{ {{u_0},{\phi _x}} \right\} = \displaystyle\sum\limits_{m = 1}^\infty {\displaystyle\sum\limits_{n = 1}^\infty {\left\{ {{U_{mn}},{\Phi _{xmn}}} \right\}\frac{{\partial {X_m}(x)}}{{\partial x}}{Y_n}(y)} } {{\rm{e}}^{{\rm{i}}\omega t}}, \\ &\left\{ {{v_0},{\phi _y}} \right\} = \displaystyle\sum\limits_{m = 1}^\infty {\displaystyle\sum\limits_{n = 1}^\infty {\left\{ {{V_{mn}},{\Phi _{ymn}}} \right\}} } {X_m}(x)\frac{{\partial {Y_n}(y)}}{{\partial y}}{{\rm{e}}^{{\rm{i}}\omega t}}, \\ &\left\{ {{w_0}} \right\} = \displaystyle\sum\limits_{m = 1}^\infty {\displaystyle\sum\limits_{n = 1}^\infty {\left\{ {{W_{0mn}}} \right\}} } {X_m}(x){Y_n}(y){{\rm{e}}^{{\rm{i}}\omega t}} \end{split} (14) 其中：
{\rm{i}} = \sqrt { - 1} ；\omega 为波纹夹芯板系统自由振动下的固有频率。满足波纹夹芯板在4种边界条件下的位移分量形式可以用表2的函数表示。表 2 不同边界条件下的函数Xm(x) 和 Yn(y)Table 2. Functions Xm(x) and Yn(y) for different boundary conditionsBoundary conditions Function Xm(x) Function Yn(y) SSSS \sin \alpha x \sin \beta x CCCC 1 - \cos 2\alpha x 1 - \cos 2\beta x CCSS 1 - \cos 2\alpha x \sin \beta x CSSS \sin \alpha x(\cos \alpha x - 1) \sin \beta x Notes: SSSS—Four sides simply supported; CCCC—Four sides clamped; CCSS—Opposite sides simply supported and clamped; CSSS—One side fixed and three edges clamped; { {\alpha = m{\text{π}}} / a}; { {\beta = n{\text{π}} } / b}; m, n—Half-wave numbers in two orthogonal coordinate directions respectively. 将式(14)代入式(13)可得系统的特征方程，简写为
\left\{ {{{K}} - {\omega ^2}{{M}}} \right\}\left\{ {{\delta}} \right\} = 0 (15) 其中：M和K分别为波纹夹芯板系统的质量矩阵和刚度矩阵；
{\left\{ {\bf{\delta }} \right\}^{\bf{T}}} = \left\{ {{U_{mn}},{V_{mn}},{\Phi _x}_{mn},{\Phi _y}_{mn},{W_{0mn}}} \right\} 为系统的振动幅值。由于它们的取值是任意性的，因此令该特征方程的系数行列式为0，求解该代数方程，即可得到波纹夹芯板自由振动时的固有频率。3. 数值算例及讨论
根据以上理论模型，编写程序计算波纹夹芯板在不同板理论下的固有频率，并且与ABAQUS有限元仿真结果进行对比，验证理论模型的正确性。此外，计算了不同边界条件下波纹夹芯板的基频，分析边界条件对波纹夹芯板振动特性的影响。同时，研究在四种不同边界条件下波纹夹芯板材料参数及结构几何参数的变化对波纹夹芯板基频的影响。
3.1 波纹夹芯板的固有频率
设波纹夹芯板的长a=240 mm，宽b=188.3 mm，上、下面板厚度hf=1 mm，波纹与面板的夹角θ=45°，波纹壁厚tc=1 mm，芯层厚度hc=8 mm。波纹夹芯板三层所采用的基体均为铝材料，其材料参数为：杨氏模量Es=71 GPa，泊松比νs=0.3，剪切模量Gs=Es/2(1+νs)，密度ρs=2 810 kg/m3。
对波纹夹芯板系统进行无量纲化，无量纲化的频率可以定义为
\varpi = \frac{{\omega {a^2}}}{h}\sqrt {\frac{{12{\rho _{\rm{s}}}(1 - \nu_{\rm{s}}^2)}}{{{E_{\rm{s}}}}}} (16) 表3为5种不同板理论所求得的波纹夹芯板在四边简支边界条件下的前五阶固有频率。可以看出，采用SSDT、TSDT、ESDT理论所求得的固有频率与有限元误差较小，全部小于3%，FSDT由于需要剪切修正因子，所求结果误差比以上三个理论略偏大，CLPT由于未考虑横向剪切变形，高估了固有频率，所求结果相比其余4种剪切变形理论，误差要大很多。其中，在基频处，波纹夹芯板的运动和变形相对稳定，剪切应力对板变形的影响相对较小，因此在五种理论中，所求基频与ABAQUS有限元仿真相比误差较小。随着固有频率的增加，波纹夹芯板的变形越来越大，运动也愈加剧烈，因此剪切应力影响变大，考虑剪切变形且满足上下表面处面力自由的SSDT、TSDT和ESDT相对于CLPT和FSDT所求结果更准确。
表 3 四边简支波纹夹芯板在不同理论下的固有频率理论解与有限元仿真结果Table 3. Theoretical solutions and finite element simulation results of natural frequency of simply supported corrugated sandwich plates using different theoriesMode ABAQUS CLPT FSDT SSDT TSDT ESDT Result Error/% Result Error/% Result Error/% Result Error/% Result Error/% (1,1) 31.18 31.72 1.75 31.41 0.77 31.00 −0.57 31.01 −0.54 30.99 −0.58 (2,1) 65.39 69.52 6.32 68.09 4.13 65.82 0.66 65.88 0.75 65.79 0.61 (1,2) 86.48 88.51 2.36 86.23 −0.29 83.74 −3.16 83.78 −3.12 83.74 −3.17 (2,2) 116.13 125.70 8.24 121.21 4.38 115.60 −0.45 115.72 −0.35 115.57 −0.48 (3,1) 118.46 132.06 11.48 127.11 7.30 119.24 0.66 119.45 0.84 119.12 0.56 在四边简支边界条件下，波纹夹芯板前五阶振型模态如图3所示。
![]()
图 3 波纹夹芯板前五阶振型图Figure 3. The first five mode shapes of the corrugated sandwich plates在四种不同边界条件下，由不同板理论计算的波纹夹芯板在自由振动时的基频与有限元仿真结果如表4所示。可以看到，理论解与有限元仿真结果相比，CLPT由于忽略横向剪切应力导致所求结果稍微偏大，考虑横向剪切变形且满足上下表面处面力自由的SSDT、TSDT、ESDT所求结果比FSDT和CLPT的误差小。因为基频处波纹夹芯板的变形相对高频小，剪应力的影响有限，所以5种板理论所求结果误差都在工程允许的误差范围之内，由此也可验证四种边界条件下所假设的位移函数的正确性。比较四种边界条件下所得到的固有频率可以发现，波纹夹芯板在CCCC边界条件下的基频最大，SSSS边界条件下的基频最小，CCCC边界条件下的基频基本可以达到SSSS的1倍左右，CCSS和CSSS边界条件下的固有频率介于上述二者之间，并且两者相差较小，即四种边界条件下波纹夹芯板的基频的关系为CCCC>CCSS>CSSS>SSSS，由此可以得知，固支边界条件会使波纹夹芯板系统在自由振动时的基频增大。同样这也说明了波纹夹芯板边界的约束越多，系统结构整体刚度越大，导致整个系统频率特征值增大。因此，可以通过设置不同的边界条件来调整波纹夹芯板的固有振动频率。
表 4 不同边界条件和板理论下波纹夹芯板的基频理论解与有限元仿真结果Table 4. Theoretical solutions and finite element simulation results of fundamental frequency of corrugated sandwich plates with different boundary conditions and plate theoriesBoundary conditons ABAQUS CLPT FSDT SSDT TSDT ESDT Result Error/% Result Error/% Result Error/% Result Error/% Result Error/% SSSS 31.18 31.72 1.75 31.41 0.77 31.00 −0.57 31.01 −0.54 30.99 −0.58 CCCC 56.42 60.46 7.17 58.90 4.40 56.99 1.01 57.03 1.08 56.98 0.99 CCSS 40.47 43.19 6.71 42.44 4.87 41.23 1.88 41.26 1.96 41.22 1.84 CSSS 39.50 42.14 6.70 41.56 5.22 40.66 2.94 40.68 2.99 40.65 2.91 4种不同边界条件下的波纹夹芯板的一阶振型模态如图4所示。
![]()
图 4 不同边界条件下波纹夹芯板一阶振型图Figure 4. The first mode shapes of corrugated sandwich plates under different boundary conditions3.2 参数变化对波纹夹芯板振动特性的影响
波纹夹芯板的材料参数和结构几何参数对其振动有着重要的影响，通过进一步研究波纹夹芯板的振动特性，为其在工程应用方面提供足够的依据。本节将基于ESDT，研究四种不同边界条件下，波纹夹芯板材料参数和结构几何参数的变化对系统基频的影响。
3.2.1 材料参数变化
波纹夹芯板的上面板、下面板、波纹芯层选用不同的基体材料：Ti和Al，研究不同材料组合对波纹夹芯板固有频率的影响。Al的物理参数在前文已给出，Ti的材料参数为弹性模量Et=177 GPa，泊松比vt = 0.32，密度ρt=4 540 kg/m3。波纹夹芯板的组合方式为Ti-Al-Ti、Ti-Ti-Ti、Al-Al-Al和Al-Ti-Al。板的其它尺寸与上述算例一致，得到四种边界条件下不同材料组合的波纹夹芯板的基频如表5所示。可知，对于每一种边界条件，基于ESDT理论计算所得的四种材料组合下波纹夹芯板的基频大小依次为Ti-Al-Ti> Ti-Ti-Ti> Al-Al-Al> Al-Ti-Al，其中不同材料组合的频率在CCCC边界条件下变化稍大一些，但几种边界条件下的变化趋势基本一致。由式(15)可知，增大材料的弹性模量，即抗变形能力增强，波纹夹芯板结构的刚度增大，系统整体的频率会增大；增大材料的密度，波纹夹芯板的质量增大，波纹夹芯板的频率会降低，但对于波纹夹芯板的上、下面板，弹性模量的影响更显著；对于芯层，密度影响起主导作用。在工程应用中可以选取适当的材料组合以提高波纹夹芯板的固有频率。
表 5 不同材料组合下波纹夹芯板的基频Table 5. Fundamental frequencies of corrugated sandwich panels with different material combinationsBoundary conditon Al-Al-Al Ti-Ti-Ti Al-Ti-Al Ti-Al-Ti SSSS 30.99 38.72 29.07 40.53 CCCC 56.98 71.16 54.31 71.99 CCSS 41.22 51.46 39.96 52.38 CSSS 40.65 50.76 39.03 52.17 3.2.2 结构几何参数变化
作为一种复合结构材料，波纹夹芯板的结构参数对其振动同样有着重要的影响，为了更好地观察波纹夹芯板的基频随结构几何参数的变化趋势，以下计算不再对系统固有频率进行无量纲化，波纹夹芯板各层采用的材料均为铝。
保持波纹夹芯板的基本尺寸和芯层高度及总高度不变，波纹夹芯板在四种边界条件下的基频随波纹与面板夹角θ的变化如图5所示。可以看出，随着波纹与面板夹角θ的增大，在四种边界条件下，波纹夹芯板基频的变化趋势基本接近，都随着夹角θ的增大呈缓慢下降的趋势。由式(1)可知，波纹与面板夹角从
30^\circ 变化到80^\circ 度时，芯层弹性模量E1和剪切模量G13增大，等效密度增大，其它等效参数变化相对较小，波纹夹芯板系统的刚度虽有一定的增大，但等效密度增大更明显，导致波纹夹芯板的基频呈下降趋势。另外，四种边界条件下，CCCC和SSSS边界条件下基频相差较大，且CCCC>SSSS；CSSS和CCSS边界条件下所得到的基频位于CCCC和SSSS之间，CCSS稍大于CSSS，但随着夹角的增大两者对应的波纹夹芯板的基频差值增大。波纹与面板的夹角θ越大，即波纹折皱密度越大，铝制波纹夹芯板越接近实体铝板。由此也可以得知，铝制波纹夹芯板的固有频率大于相同尺寸的实体铝板。![]()
图 5 波纹与面板夹角对波纹夹芯板基频的影响Figure 5. Influence of the angle between corrugation and panel on the fundamental frequency of corrugated sandwich panel波纹芯层高度对波纹夹芯板基频的影响如图6所示。可知，随着芯层高度占比hc/h的增大，即芯层厚度增大，上、下面板厚度减小，4种边界条件下波纹夹芯板基频的变化趋势基本接近。由式(1)可知，随着波纹芯层高度hc的增大，波纹的结构参数(斜边长lc、半底边长p)增大，波纹芯层各个方向的弹性模量和剪切模量减小，导致系统的刚度减小，从而使基频降低，但是随着芯层高度占比hc/h的增大，波纹芯层的等效密度减小，使系统基频增大。因为在hc/h≤0.8时芯层等效密度的影响起主导作用，之后弹性模量和剪切模量的影响起主导作用，所以导致波纹夹芯板的基频先增大后减小。在hc/h=0.8附近，基频达到最大值，在hc/h>0.8时，波纹夹芯板的基频迅速下降。因此可以得知，在实际工程应用中选取适当的芯层占比可以提高波纹夹芯板的固有频率。
![]()
图 6 波纹芯层高度占比对波纹夹芯板基频的影响Figure 6. Influence of corrugated core height ratio on the fundamental frequency of corrugated sandwich panel波纹夹芯板的基频随波纹壁厚的变化如图7所示。由等效参数计算式(1)，随着壁厚的增加，波纹芯层的剪切模量和弹性模量增大，且E12呈直线增大，系统刚度增大，同样芯层的等效密度也呈直线增大，导致系统的基频减小。因为等效密度的影响起主导作用，所以波纹夹芯板的基频持续减小。另外，CCSS和CSSS两种边界条件下的基频的差值略微增大，且4种边界条件下，CCCC边界条件下波纹夹芯板的基频相对于其他3种边界条件在壁厚tc≤3 mm时下降速度略快一些，其余边界条件下，波纹夹芯板基频的下降趋势基本一致。
![]()
图 7 波纹壁厚度对波纹夹芯板基频的影响Figure 7. Influence of corrugated wall thickness on the fundamental frequency of corrugated sandwich panel波纹夹芯板厚度 h对波纹夹芯板基频的影响如图8所示。可见，波纹夹芯板芯层在各方向的弹性模量和剪切模量在h≤30时迅速减小，芯层刚度随之减小，之后几乎保持不变，导致系统的固有频率减小。但是由式(1)可知，芯层的等效密度同样在h≤30时迅速减小，之后几乎保持不变，且等效密度的影响起着主导作用。当波纹夹芯板的总厚度h<30 mm时，四种边界条件下波纹夹芯板系统的基频呈现明显增大的趋势，CCCC边界条件下的基频增长速率最快，CSSS和CCSS边界条件所对应的基频曲线接近重合，在h=15 mm附近CSSS条件下所求的基频开始高于CCSS边界条件下所求的基频。当h>30 mm时，四条曲线都趋于平稳，波纹夹芯板基频变化波动很小，此时板的宽厚比b/h<6，波纹夹芯板由薄板逐渐变为厚板。根据以上分析可以得知，在工程应用中适当提高波纹夹芯板的厚度可以提高系统的固有频率。
![]()
图 8 波纹夹芯板厚度h 对波纹夹芯板基频的影响Figure 8. Influence of corrugated plate thickness on the fundamental frequency of corrugated sandwich panel设k为变化比例因子，将波纹夹芯板结构的所有几何尺寸分别乘以比例因子k，即波纹夹芯板结构整体变大或者缩小，得到的波纹夹芯板基频的变化曲线如图9所示。随着k的变化，波纹夹芯板芯层的弹性模量、剪切模量和密度保持不变。由于波纹夹芯板长和宽的增大，导致整个系统质量增加，因此系统的固有频率将减小。图中所取变化比例因子k为0.2~3，四种边界下波纹夹芯板的基频在
k\!\leqslant\!1 时急剧下降，CCSS和CSSS曲线基本重合；在k \!\geqslant\!1 时，四种边界条件下的波纹夹芯板系统的基频越来越接近。由此可知，当波纹夹芯板尺寸无限大时，四种边界条件下波纹夹芯板系统的基频会接近相等，且波纹夹芯板尺寸越大系统的基频就越小。![]()
图 9 比例因子k 对波纹夹芯板基频的影响Figure 9. Influence of scale factor on the fundamental frequency of corrugated sandwich panel4. 结 论
研究了波纹夹芯板在四种边界条件下的自由振动特性，对比了不同板理论所求得的固有频率，并与有限元仿真结果进行了对比。此外，基于指数剪切变形理论(ESDT)分析了波纹夹芯板材料参数和结构几何参数的变化对系统固有频率的影响。
(1)对于基尔霍夫经典板理论(CLPT)、一阶剪切变形理论(FSDT)、正弦剪切变形理论(SSDT)、三阶剪切变形理论(TSDT)、ESDT五种板理论，其中CLPT由于不考虑板的横向剪切变形，因此求得系统自由振动的固有频率误差最大，FSDT所求固有频率次之，SSDT、TSDT、ESDT所求固有频率相近，且与有限元仿真结果相比，误差最小。
(2)四种边界条件下波纹夹芯板自由振动的基频大小依次为：四边固支(CCCC)>对边简支和固支(CCSS)>一边固支三边简支(CSSS)>四边简支(SSSS)，且CCCC的基频比SSSS的基频约大1倍, 即波纹夹芯板结构的边界约束越多，系统的固有频率越大。
(3)增大波纹夹芯板材料的弹性模量会导致其刚度变大，进而使系统的基频增大；增大材料密度会导致波纹夹芯板的质量增大，进而使系统的基频减小，对于上、下面板，弹性模量的影响较显著，而对于夹芯层，密度的影响起主导作用。
(4)波纹夹芯板的结构几何参数对系统的振动有着重要影响。随着波纹与面板的夹角
\theta 或波纹壁厚tc的增大，波纹夹芯板的基频变小；随着芯层占比hc/h的增大，波纹夹芯板的基频先增大后减小；保持hc/h、面板占比hf/h不变，随着波纹夹芯板厚度h的增大，波纹夹芯板的基频增大明显；随着波纹夹芯板尺寸的增大其基频下降，且四种边界条件下波纹夹芯板的基频趋近相等。 -
![]()
图 1 波纹夹芯板的模型(a)和波纹单胞示意图(b)
Figure 1. Model diagrams of corrugated sandwich panel (a) and corrugated cell (b)
a—Length of corrugated sandwich panel; b—Width of corrugated sandwich panel; lc—Length of the hypotenuse; hc—Height of core layer; tc—Wall thickness; θ—Corrugation angle; p—Length of bottom side
![]()
图 3 波纹夹芯板前五阶振型图
Figure 3. The first five mode shapes of the corrugated sandwich plates
![]()
图 4 不同边界条件下波纹夹芯板一阶振型图
Figure 4. The first mode shapes of corrugated sandwich plates under different boundary conditions
![]()
图 5 波纹与面板夹角对波纹夹芯板基频的影响
Figure 5. Influence of the angle between corrugation and panel on the fundamental frequency of corrugated sandwich panel
![]()
图 6 波纹芯层高度占比对波纹夹芯板基频的影响
Figure 6. Influence of corrugated core height ratio on the fundamental frequency of corrugated sandwich panel
![]()
图 7 波纹壁厚度对波纹夹芯板基频的影响
Figure 7. Influence of corrugated wall thickness on the fundamental frequency of corrugated sandwich panel
![]()
图 8 波纹夹芯板厚度
h 对波纹夹芯板基频的影响Figure 8. Influence of corrugated plate thickness on the fundamental frequency of corrugated sandwich panel
![]()
图 9 比例因子
k 对波纹夹芯板基频的影响Figure 9. Influence of scale factor on the fundamental frequency of corrugated sandwich panel
表 1 不同板理论对应的剪切形状函数
Table 1 Shear shape functions corresponding to different plate theories
Shear theory Function f ({\textit{z}} ) CLPT f({\textit{z}}) = 0 FSDT f({\textit{z} }) = {\textit{z}} SSDT f({\textit{z}}) = \dfrac{h}{{\text{π}}}\sin \left(\dfrac{{{\text{π}}{\textit{z}}} }{h}\right) TSDT f({\textit{z}}) = {\textit{z}}\left[1 - \dfrac{4}{3}{\left(\dfrac{{\textit{z}}}{h}\right)^2}\right] ESDT f({\textit{z}}) = {\textit{z}}{ {\rm{e} }^{ - 2{ {\left( {\dfrac{{\textit{z}}}{h} } \right)}^2} } } Notes: CLPT—Classical plate theory; FSDT—First-order shear plate theory; SSDT—Sinusoidal shear deformation theory; TSDT—Third-order shear deformation theory; ESDT—Exponential shear deformation theory. 
下载: 导出CSV
表 2 不同边界条件下的函数Xm(x) 和 Yn(y)
Table 2 Functions Xm(x) and Yn(y) for different boundary conditions
Boundary conditions Function Xm(x) Function Yn(y) SSSS \sin \alpha x \sin \beta x CCCC 1 - \cos 2\alpha x 1 - \cos 2\beta x CCSS 1 - \cos 2\alpha x \sin \beta x CSSS \sin \alpha x(\cos \alpha x - 1) \sin \beta x Notes: SSSS—Four sides simply supported; CCCC—Four sides clamped; CCSS—Opposite sides simply supported and clamped; CSSS—One side fixed and three edges clamped; { {\alpha = m{\text{π}}} / a}; { {\beta = n{\text{π}} } / b}; m, n—Half-wave numbers in two orthogonal coordinate directions respectively.
下载: 导出CSV
表 3 四边简支波纹夹芯板在不同理论下的固有频率理论解与有限元仿真结果
Table 3 Theoretical solutions and finite element simulation results of natural frequency of simply supported corrugated sandwich plates using different theories
Mode ABAQUS CLPT FSDT SSDT TSDT ESDT Result Error/% Result Error/% Result Error/% Result Error/% Result Error/% (1,1) 31.18 31.72 1.75 31.41 0.77 31.00 −0.57 31.01 −0.54 30.99 −0.58 (2,1) 65.39 69.52 6.32 68.09 4.13 65.82 0.66 65.88 0.75 65.79 0.61 (1,2) 86.48 88.51 2.36 86.23 −0.29 83.74 −3.16 83.78 −3.12 83.74 −3.17 (2,2) 116.13 125.70 8.24 121.21 4.38 115.60 −0.45 115.72 −0.35 115.57 −0.48 (3,1) 118.46 132.06 11.48 127.11 7.30 119.24 0.66 119.45 0.84 119.12 0.56
下载: 导出CSV
表 4 不同边界条件和板理论下波纹夹芯板的基频理论解与有限元仿真结果
Table 4 Theoretical solutions and finite element simulation results of fundamental frequency of corrugated sandwich plates with different boundary conditions and plate theories
Boundary conditons ABAQUS CLPT FSDT SSDT TSDT ESDT Result Error/% Result Error/% Result Error/% Result Error/% Result Error/% SSSS 31.18 31.72 1.75 31.41 0.77 31.00 −0.57 31.01 −0.54 30.99 −0.58 CCCC 56.42 60.46 7.17 58.90 4.40 56.99 1.01 57.03 1.08 56.98 0.99 CCSS 40.47 43.19 6.71 42.44 4.87 41.23 1.88 41.26 1.96 41.22 1.84 CSSS 39.50 42.14 6.70 41.56 5.22 40.66 2.94 40.68 2.99 40.65 2.91
下载: 导出CSV
表 5 不同材料组合下波纹夹芯板的基频
Table 5 Fundamental frequencies of corrugated sandwich panels with different material combinations
Boundary conditon Al-Al-Al Ti-Ti-Ti Al-Ti-Al Ti-Al-Ti SSSS 30.99 38.72 29.07 40.53 CCCC 56.98 71.16 54.31 71.99 CCSS 41.22 51.46 39.96 52.38 CSSS 40.65 50.76 39.03 52.17
下载: 导出CSV
-
[1] 岳灿甫, 吴始栋. 国外船用激光焊接波纹夹芯板的开发与应用[J]. 鱼雷技术, 2007, 15(4):1-5. YUE C F, WU S D. Introduction to development and applications of marine laser welded lightweight sandwich panel in the US navy and european countries[J]. Torpedo Technology,2007,15(4):1-5(in Chinese).
[2] 贺梦豪, 吴昊, 程远胜. 波纹夹芯杂交夹层板入水砰击动力响应特性研究[J]. 舰船科学技术, 2016, 38(8):11-17. HE M H, WU H, CHENG Y S. Dynamic response analysis of hybrid sandwich plate with corrugated-cores subjected to slam-ming impact[J]. Ship Science and Technology,2016,38(8):11-17(in Chinese).
[3] 孙小曼. 车用波纹夹芯板弯曲性能研究[D]. 大连: 大连理工大学, 2013. SUN X M. Analysis of bending performance of vehicle corrugated core sandwich structure[D]. Dalian: Dalian University of Technology, 2013(in Chinese).
[4] KUJALA P, KLANAC A. Steel sandwich panels in marine applications[J]. Brodogradnja,2005,56(4):305-314.
[5] PENG L X, LIEW K M, KITIPORNCHAI S. Analysis of stiffened corrugated plates based on the FSDT via the mesh-free method[J]. International Journal of Mechanical Sciences,2006,49(3):364-378.
[6] SEMENYUK N P, NESKHODOVSKAYA N A. On design models in stability problems for corrugated cylindrical shells[J]. International Applied Mechanics,2002,38(10):1245-1252. DOI: 10.1023/A:1022214630581
[7] BRIASSOULIS D. Equivalent orthotropic properties of corrugated sheets[J]. Computers & Structures,1986,23(2):129-138.
[8] SOHRAB K, DAN Z. Corrugated all-composite sandwich structures. Part 2: Failure mechanisms and experimental programme[J]. Composites Science and Technology,2009,69(7):920-925.
[9] TIAN Y S, LU T J. Optimal design of compression corrugated panels[J]. Thin-Walled Structures,2005,43(3):477-498. DOI: 10.1016/j.tws.2004.07.014
[10] SAMANTA A, MUKHOPADHYAY M. Finite element static and dynamic analyses of folded plates[J]. Engineering Structures,1999,21(3):277-287. DOI: 10.1016/S0141-0296(97)90172-3
[11] XIA Y, FRISWELL M I, FLORES E I S. Equivalent models of corrugated panels[J]. International Journal of Solids and Structures,2012,49(13):1453-1462. DOI: 10.1016/j.ijsolstr.2012.02.023
[12] LI X, WANG Z, ZHU F, et al. Response of aluminium corrugated sandwich panels under air blast loadings: Experiment and numerical simulation[J]. International Journal of Impact Engineering,2014,65:79-88. DOI: 10.1016/j.ijimpeng.2013.11.002
[13] 王红霞, 王德禹, 李喆. 三角形夹芯板夹心层的等效弹性常数[J]. 固体力学学报, 2007, 28(2):178-182. DOI: 10.3969/j.issn.0254-7805.2007.02.011 WANG H X, WANG D Y. LI Z. Equivalent elastic constants of truss-core sandwich plates[J]. Acta Mechanica Solida Sinica,2007,28(2):178-182(in Chinese). DOI: 10.3969/j.issn.0254-7805.2007.02.011
[14] 王青伟, 赵才其. 三角形桁架夹芯层等效弹性常数研究和夹芯板参数优化设计[J]. 特种结构, 2010, 27(5):61-66. DOI: 10.3969/j.issn.1001-3598.2010.05.015 WANG Q W, ZHAO C Q. Research of equivalent elastic constants of triangle truss-core sandwich panels and optimization design of constants of core sandwich panels[J]. Special Structures,2010,27(5):61-66(in Chinese). DOI: 10.3969/j.issn.1001-3598.2010.05.015
[15] 吴建均, 王永静, 谢敬尧, 等. 泡沫波纹复合夹心梁等效弹性常数的推导及其振动特性研究[J]. 西安交通大学学报, 2018, 52(11):73-80. WU J J, WANG Y J, XIE J Y, et al. Derivation of equivalent elastic constant and vibration analysis of foam-corrugated composite sandwich beam[J]. Journal of Xi'an Jiaotong University,2018,52(11):73-80(in Chinese).
[16] 崔岸, 刘芳芳, 张晗, 等. 车身泡沫填充铝合金波纹夹芯板结构性能分析与优化[J]. 汽车工程, 2019, 41(10):1221-1227. CUI A, LIU F F, ZHANG H, et al. Performance analysis and optimization of foam-filled aluminum-alloy corrugated sandwich panel structure for vehicle body[J]. Automotive Engineering,2019,41(10):1221-1227(in Chinese).
[17] 李静雯, 张博明, 孙义亮, 等. 不同铺层方式下连续玻璃纤维/聚丙烯复合材料波纹夹芯板的力学性能[J]. 复合材料学报, 2019, 36(5):1074-1082. LI J W, ZHANG B M, SUN Y L, et al. Mechanical properties of continuous glass fiber/polypropylene corrugated sandwich boards under different laminates[J]. Acta Materiae Compositae Sinica,2019,36(5):1074-1082(in Chinese).
[18] 刘昆, 王哲, 王自力. 波纹夹层板冲击响应理论计算方法研究[J]. 振动与冲击, 2019, 38(2):90-97. LIU K, WANG Z, WANG Z L. Theoretical calculation method for the impact responses of corrugated sandwich panels[J]. Journal of Vibration and Shock,2019,38(2):90-97(in Chinese).
[19] REDDY J N. Mechanics of laminated composite plates and shells: Theory and analysis[M]. New York: CRC Press, 2004.
[20] CROCE L D, VENINI P. Finite elements for functionally graded Reissner-Mindlin plates[J]. Computer Methods in Applied Mechanics and Engineering,2004,193(9–11):705-725.
[21] MATSUNAGA H. Free vibration and stability of functionally graded plates according to a 2-D higher-order deformation theory[J]. Composite Structures,2008,82(4):499-512. DOI: 10.1016/j.compstruct.2007.01.030
[22] TALHA M, SINGH B N. Static response and free vibration analysis of FGM plates using higher order shear deformation theory[J]. Applied Mathematical Modelling,2010,34(12):3991-4011. DOI: 10.1016/j.apm.2010.03.034
[23] WU C P, LI H Y. An RMVT-based third-order shear deformation theory of multilayered functionally graded material plates[J]. Composite Structures,2010,92:2591-2605. DOI: 10.1016/j.compstruct.2010.01.022
[24] THAI H T, VO T P. A new sinusoidal shear deformation theory for bending, buckling, and vibration of functionally graded plates[J]. Applied Mathematical Modelling,2013,37:3269-3281. DOI: 10.1016/j.apm.2012.08.008
[25] 曹源, 雷剑. 基于正弦剪切变形理论的功能梯度材料三明治微梁的静动态特性[J]. 复合材料学报, 2020, 37(1):223-235. CAO Y, LEI J. Static and dynamic properties of functionally graded materials sandwich microbeams based on sinusoidal shear deformation theory[J]. Acta Materiae Compositae Sinica,2020,37(1):223-235(in Chinese).
[26] KARAMA M, AFAQ K S, MISTOU S. Mechanical behaviour of laminated composite beam by the new multi-layered laminated composite structures model with transverse shear stress continuity[J]. International Journal of Solids and Structures,2003,40(6):1525-1546. DOI: 10.1016/S0020-7683(02)00647-9
[27] REDDY J N. Energy principles and variational methods in applied mechanics[M]. New York: John Wiley & Sons Inc, 2002.
-
期刊类型引用(4)
1. 于丹. 聚氯乙烯/碳纳米管复合材料的制备和性能研究. 塑料科技. 2024(01): 36-39 . 
百度学术
2. 余澎，涂操，郭博森，王闻达，赵航，彭玉婷，罗卫华. 木质素基碳纳米管/炭复合材料的制备及电化学性能研究. 现代化工. 2023(02): 92-97 . 百度学术
3. 冀佳帅，杜佳琪，陈俊琳，张新民，刘伟，宋朝霞. Co-Fe普鲁士蓝/多壁碳纳米管复合材料的超电容性能. 材料科学与工艺. 2023(04): 1-8 . 百度学术
4. 张开砚. 电感耦合等离子体发射光谱法测定普鲁士蓝类正极材料中铁和钠. 化学分析计量. 2022(08): 26-30 . 百度学术
其他类型引用(5)
-
计量
- 文章访问数: 1476
- HTML全文浏览量: 522
- PDF下载量: 105
- 被引次数: 9
