力学学报 2018 , 50 (4): 847-852 https://doi.org/10.6052/0459-1879-18-083

Orginal Article

多层复合壳体三维振动分析的谱--微分求积混合法

叶天贵^*^,, 靳国永, 刘志刚

哈尔滨工程大学动力与能源工程学院,哈尔滨 150001

A SPECTRAL-DIFFERENTIAL QUADRATURE METHOD FOR 3-D VIBRATION ANALYSIS OF MULTILAYERED SHELLS

Ye Tiangui^*^,, Jin Guoyong, Liu Zhigang

College of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China

中图分类号: TH113,TB33

通讯作者: ^*叶天贵,副教授,主要研究方向：振动噪声控制. E-mail:yetiangui@hrbeu.edu.cn^*叶天贵,副教授,主要研究方向：振动噪声控制. E-mail:yetiangui@hrbeu.edu.cn

收稿日期: 2018-03-20

接受日期: 2018-03-20

网络出版日期: 2018-07-18

基金资助: 国家自然科学基金(51709066, 51775125)和中国博士后科学基金(2017M621252)资助项目

展开

摘要

对于较厚的多层复合壳体,其振动位移沿厚度方向呈锯齿形变化且层间剪切和拉、压应力呈三维耦合状态,采用传统的等效单层理论分析已不能满足精度要求. 建立不受结构厚度、铺层材料性质和铺层方式限制的三维分析方法具有重要的研究价值. 本文以独立铺层为建模对象,结合广义谱方法与微分求积技术建立了一种适用一般边界条件和铺层方式的多层复合壳体三维分析新方法——谱--微分求积混合法. 该方法应用三维弹性理论对独立铺层进行精确建模,有效克服了二维简化理论对横向变形以及层间应力估计不确切的缺点;引入微分求积技术对铺层进行数值离散,将三维偏微分问题转化为二维偏微分问题,降低了求解维度和难度;应用广义谱方法近似地表述离散计算面上的场变量,将获取的二维偏微分方程转化为以场变量谱展开系数为未知量的线性代数方程组,避免了对超越方程的求解. 数值验证结果表明该方法收敛性好,计算精度高.

关键词： 多层复合壳体 ; 三维振动 ; 谱--微分求积混合法 ; 一般边界条件

Abstract

The equivalent single layer (ESL) theories can be grossly in error for predicting vibration characteristics of thick multilayered shells because the vibration displacement and stress field of such shells under vibration are in full 3-D coupling condition. It is necessary to develop more accurate and efficient methods which are capable of dealing with multilayered structures with different boundary conditions, general laminations as well as arbitrary thickness universally. In order to overcome the drawback of the existing three-dimensional methods that are only confined for very limited cases such as cross-ply laminated rectangular plates under simply-supported boundary conditions, a general spectral-differential quadrature method is proposed. This method is undertaken by the exact 3-D elasticity theory so that it’s able to study very well the dynamic behavior of thick multilayered structures which cannot be provided by the 2-D ESL theories. In each individual layer, the transverse domain is discretized by the differential quadrature technique. The displacement fields of the discretized surfaces are selected as fundamental unknowns. Then, each fundamental unknown is invariantly expanded by the general spectral method as a series of complete, orthogonal polynomials. The problems are stated in a variational form by the aid of penalty parameters which provides complete flexibilities to describe any prescribed boundary conditions. The current method can successfully avoid solving a highly nonlinear transcendental equation that is rely on roots-locating numerical method and all the modal information can be obtained just by solving linear algebraic equation systems. Numerical verification shows that the proposed method has high calculation precision. The method can be directly extend to the static and dynamic analysis of multilayered shells as well.

Keywords： multilayered shells ; 3-D vibration ; spectral-differential quadrature method ; general boundary restraints

PDF (438KB) 元数据多维度评价相关文章收藏文章

本文引用格式导出 EndNote Ris Bibtex

叶天贵, 靳国永, 刘志刚. 多层复合壳体三维振动分析的谱--微分求积混合法[J]. 力学学报, 2018, 50(4): 847-852 https://doi.org/10.6052/0459-1879-18-083

Ye Tiangui, Jin Guoyong, Liu Zhigang. A SPECTRAL-DIFFERENTIAL QUADRATURE METHOD FOR 3-D VIBRATION ANALYSIS OF MULTILAYERED SHELLS[J]. Acta Mechanica Sinica, 2018, 50(4): 847-852 https://doi.org/10.6052/0459-1879-18-083

引言

多层复合结构广泛应用于航空航天、高速列车、防卫装备等高科技行业,其振动特性分析一直是相关领域的研究焦点之一^[1,2,3]. 常用的多层复合结构振动建模理论主要有：经典等效单层理论、剪切变形理论、锯齿理论和三维层合理论四大类^{[4,5,6,7,8,9,10,11,12,13]}. 经典等效单层理论和剪切变形理论都是传统的板壳理论对多层复合结构的直接推广应用. 这两种理论先在厚度方向上把多层复合结构等效成各向异性的单一薄层,然后用结构中面的位移或者应力分量描述结构的整体振动行为. 由于组成材料的各向异性及呈层性等特点,多层复合结构受载后的应力状态和振动行为十分复杂,其振动位移沿厚度方向呈锯齿形分布且层间剪切和拉、压应力呈三维耦合状态. 因此,当结构厚度比较大或者各层材料差异较为明显时,位移连续但层间应力不连续的等效单层理论的建模误差急剧增大. 锯齿理论是以上述两类理论为基础,结合层间应力连续条件而得到的理论. 这类理论一般通过在结构的面内位移分量添加锯齿函数从而使横向剪切应力满足层间连续条件. 但该理论一般不满足上下表面自由边界条件也无法得到准确的层间应力^[12]. 直接应用三维层合理论可以满足位移和横向应力层间连续条件,这一理论和等效单层理论及锯齿理论最大的区别在于结构的每一个子层都看作是用独立的函数来描述的三维实体. 众所周知,只有边界条件和铺层方式简单的少数规则结构才能得到解析解,对其他较为复杂的边界条件、铺层方式以及结构形式,只能借助数值方法或者其他近似解法来求解^[14].

针对已有的大多数三维求解方法仅能求解少数形状规则,边界条件简单且铺层方式正交铺设的多层复合结构的限制,本文以独立铺层为建模对象,结合广义谱方法与微分求积技术建立了一种适用一般边界条件和铺层方式的多层复合壳体三维分析半解析方法——谱--微分求积混合法,以期为多层复合壳体振动计算分析提供高效可靠的参考方法.

1 多层复合壳体三维振动力学模型

1.1 分析模型

考虑如图1所示的多层复合壳体. 假设参考坐标系为正交曲线坐标系 $αβζ$ . 其中, $α$ 和 $β$ 沿底面的主曲率线方向, $ζ$ 沿厚度方向; $u$ , $v$ 和 $w$ 为壳体上任意一点沿 $α$ , $β$ 和 $ζ$ 方向的位移分量; $R_{α}$ 和 $R_{β}$ 代表 $α$ 和 $β$ 方向的主曲率半径.

显示原图| 下载原图ZIP| 生成PPT

图1 多层复合壳体模型及参考坐标系.

Fig. 1 Schematic illustration of a multilayered shell and the curvilinear coordinate system

1.2 几何方程及物理方程

当壳体因振动产生微小位移时,壳体上任意点的6个应变分量为

其中, $ε_{α}$ , $ε_{β}$ , $ε_{ζ}$ 分别表示主应变, $γ_{αβ}$ , $γ_{αζ}$ , $γ_{βζ}$ 为面内及横向剪切应变, $A_{α}$ 和 $A_{β}$ 为Lamè常数^[15].

1.3 计算面选取及应变应变分量计算

本文方法先从独立子层入手,沿子层厚度方向配置若干计算面并引入微分求积技术将结构的三维偏微分问题转化为二维偏微分问题进行求解.

假设第 $k$ 子层上下表面横向坐标分别为 $ζ_{k + 1}$ 和 $ζ_{k}$ . 为了提高方法收敛性和计算速度,本文所选取的各计算面满足切比雪夫零点分布形式,各计算面横向坐标为

$\begin{array}{l} ζ_{k}^{1} = ζ_{k} ζ_{k}^{J_{k}} = ζ_{k + 1} ζ_{k}^{j} = \frac{(ζ_{k}^{1} + ζ_{k}^{J_{k}})}{2} - \frac{h_{k}}{2} \cos (\frac{2 j π - 3 π}{2 J_{k} - 4}), j = 2,3, \dots \end{array}\} (2)$

式中, $J_{k}$ 表示计算面总数, $h_{k}$ 为第 $k$ 子层厚度.

将第 $k$ 子层第 $j$ 计算面的位移分量、Lamè常数和应变分量定义为

$\begin{array}{l} u (ζ_{k}^{j}) = u_{k}^{j}, & A_{α} (ζ_{k}^{j}) = A_{α k}^{j} v (ζ_{k}^{j}) = v_{k}^{j}, \\ A_{β} (ζ_{k}^{j}) = A_{β k}^{j} w (ζ_{k}^{j}) = w_{k}^{j}, & A_{ζ} (ζ_{k}^{j}) = A_{ζ k}^{j} \end{array}\} (3)$

$\begin{array}{l} ε_{α} (ζ_{k}^{j}) = ε_{α k}^{j}, & γ_{βζ} (ζ_{k}^{j}) = γ_{βζ k}^{j} ε_{β} (ζ_{k}^{j}) = ε_{β k}^{j}, \\ γ_{αζ} (ζ_{k}^{j}) = γ_{αζ k}^{j} ε_{ζ} (ζ_{k}^{j}) = ε_{ζ k}^{j}, & γ_{αβ} (ζ_{k}^{j}) = γ_{αβ k}^{j} \end{array}\} (4)$

因此,将式 (3)和式(4)代入式(1)并进行微分求积可得 $\begin{array}{l} ε_{α k}^{j} = \frac{1}{A_{α k}^{j}} \frac{\partial u_{k}^{j}}{\partial α} + \frac{v_{k}^{j}}{A_{α k}^{j} A_{β k}^{j}} \frac{\partial A_{α k}^{j}}{\partial β} + \frac{w_{k}^{j}}{A_{α k}^{j}} \sum_{r} ‍ a_{kjr} A_{α k}^{r} ε_{β k}^{j} = \frac{u_{k}^{j}}{A_{α k}^{j} A_{β k}^{j}} \frac{\partial A_{β k}^{j}}{\partial α} + \frac{1}{A_{β k}^{j}} \frac{\partial v_{k}^{j}}{\partial β} + \frac{w_{k}^{j}}{A_{β k}^{j}} \sum_{r} ‍ a_{kjr} A_{β k}^{r} ε_{ζ k}^{j} = \sum_{r} ‍ a_{kjr} w_{k}^{r} γ_{βζ k}^{j} \\ = \frac{1}{A_{β k}^{j}} \frac{\partial w_{k}^{j}}{\partial β} + A_{β k}^{j} \sum_{r} ‍ a_{kjr} \frac{v_{k}^{r}}{A_{k β}^{r *}} γ_{αζ k}^{j} = \frac{1}{A_{α k}^{j}} \frac{\partial w_{k}^{j}}{\partial α} + A_{α k}^{j} \sum_{r} ‍ a_{kjr} \frac{u_{k}^{r}}{A_{k α}^{r *}} γ_{αβ k}^{j} = \frac{A_{β k}^{j}}{A_{α k}^{j}} \frac{\partial}{\partial α} (\frac{v_{k}^{j}}{A_{β k}^{j}}) + \frac{A_{α k}^{j}}{A_{β k}^{j}} \frac{\partial}{\partial β} (\frac{u_{k}^{j}}{A_{α k}^{j}}) \end{array}\} (5)$

式中,权重系数 $a_{kjr}$ 表达式如下所示

$a_{kjr} = \{\begin{array}{l} \frac{1}{ζ_{k}^{r} - ζ_{k}^{j}} \prod_{s}^{\neq} ‍ \frac{ζ_{k}^{j} - ζ_{k}^{s}}{ζ_{k}^{r} - ζ_{k}^{s}}, & r \neq j - \sum_{r} ‍ a_{kjr}, \\ r = j \end{array} (6)$

根据拉格朗日插值定理,结构第 $k$ 子层的整体位移分布、几何分量分布及应变分布可写为

$\begin{array}{l} \{u, v, w\} = L_{k}^{j} (ζ) \sum_{j} ‍ \{u_{k}^{j}, v_{k}^{j}, w_{k}^{j}\} {A_{α}, A_{β}, A_{ζ}} = L_{k}^{j} (ζ) \sum_{j} ‍ \{A_{α k}^{j}, A_{β k}^{j}, A_{ζ k}^{j}\} \{ε_{α}, ε_{β}, ε_{ζ}\} \\ = L_{k}^{j} (ζ) \sum_{j} ‍ \{ε_{α k}^{j}, ε_{β k}^{j}, ε_{ζ k}^{j}\} \{γ_{βζ}, γ_{αζ}, γ_{αβ}\} = L_{k}^{j} (ζ) \sum_{j} ‍ \{γ_{βζ k}^{j}, γ_{αζ k}^{j}, γ_{αβ k}^{j}\} \end{array}\} (7)$

其中 $L_{k}^{j} (ζ)$ 为拉格朗日插值多项式,定义如下

$L_{k}^{j} (ζ) = \prod_{r}^{\neq} ‍ \frac{ζ - ζ_{k}^{r}}{ζ_{k}^{j} - ζ_{k}^{r}}, ζ_{k}^{1} # x 2264; ζ # x 2264; ζ_{k}^{J_{k}} (8)$

根据胡克定律,第 $k$ 子层上任一点的应力状态可由下列方程直接求得 $\begin{array}{l} σ_{k} = T_{k} C_{k} T_{k}^{T} ε_{k} σ_{k} = [σ_{α}, σ_{β}, σ_{ζ}, τ_{βζ}, τ_{αζ}, τ_{αβ}]^{T} ε_{k} = [ε_{α}, ε_{β}, ε_{ζ}, γ_{βζ}, γ_{αζ}, γ_{αβ}]^{T} \end{array}\} (9)$

式中, $σ_{α}$ , $σ_{β}$ , $σ_{ζ}$ 分别表示沿 $α$ , $β$ , $ζ$ 方向上的正应力; $τ_{αβ}$ , $τ_{αζ}$ , $τ_{βζ}$ 为面内和横向剪切应力, $C_{k}$ 为材料主方向的材料弹性常数,转换矩阵 $T_{k}$ 与材料铺层方式有关,一般为材料主坐标与结构主坐标夹角的函数^[20].

第 $k$ 子层的振动应变能及动能量方程可写为

$\begin{array}{l} U_{sk} = \frac{1}{2} \underset{α}{\int ‍} \underset{β}{\int ‍} \int_{ζ_{k}}^{ζ_{k + 1}} ‍ ε_{k}^{T} T_{k} C_{k} T_{k}^{T} v_{k} A_{α} A_{β} A_{ζ} d ζ d β d α T_{k} = \frac{1}{2} \underset{α}{\int ‍} \underset{β}{\int ‍} \int_{ζ_{k}}^{ζ_{k + 1}} ‍ ρ_{k} ({\dot{u}}^{2} + {\dot{v}}^{2} + {\dot{w}}^{2}) A_{α} A_{β} A_{ζ} d ζ d β d α \end{array}\} (10)$

借助哈密尔顿变分原理,式(10)可推导得到到结构的振动控制微分方程,具体形式参考文献[17].

1.4 子层特征矩阵及整体方程求取

多层复合壳体的边界条件可以分为3大类：几何边界条件、自然边界条件和混合边界条件. 在实际求解中,只有几何形状规则、至少一对边简支、铺层对称的极少数多层复合结构才能找到既满足控制微分方程又满足边界约束条件的解析解. 对于其他情况,可将结构的边界条件方程和控制微分方程转化成能量泛函,然后通过求解泛函的极值问题进行求解. 本文考虑几何边界条件下层合壳体的振动问题,为了适用一般边界条件,通过引入罚函数,将边界约束条件改写为如下泛函

$\begin{array}{l} Π_{pk} = \frac{1}{2} \int ‍ \int_{ζ_{k}}^{ζ_{k + 1}} ‍ \sum_{l}^{=} ‍ (k_{ku}^{α_{l}} u^{2} + k_{kv}^{α_{l}} v^{2} + k_{kw}^{α_{l}} w^{2}) |_{α = α_{l}} A_{β} A_{ζ} d ζ d β + \frac{1}{2} \int ‍ \int_{ζ_{k}}^{ζ_{k + 1}} ‍ \sum_{l}^{=} ‍ (k_{ku}^{β_{l}} u^{2} + k_{kv}^{β_{l}} v^{2} + k_{kw}^{β_{l}} w^{2}) |_{β = β_{l}} A_{α} A_{ζ} d ζ d α \end{array} (11)$

式中 $k_{ku}^{α_{l}}, k_{kv}^{α_{l}}, k_{kw}^{α_{l}}, k_{ku}^{β_{l}}, k_{kv}^{β_{l}}, k_{kv}^{β_{l}}$ 为罚参数. 上述处理方法克服了传统解析方法难以找到满足边界条件展开式的困难,极大拓宽了展开式的选择范围. 表1给出了不同边界条件下罚参数的取值情况.

表1 不同边界条件下罚参数的取值情况( $α = α l$ 边)

Table 1 Penalty parameters for various classical boundaries (at the end of $α = α l$ )

Boundary conditions	Penalty parameters
	$k_{ku}^{α_{l}}$	$k_{kv}^{α_{l}}$	$k_{kw}^{α_{l}}$
free ( $σ_{α} = 0$ ; $τ_{αβ} = 0$ ; $τ_{αζ} = 0$ )	0	0	0
simply-supported ( $σ_{α} = 0$ ; $v = 0$ ; $w = 0$ )	0	103E_1	103E_1
clamped ( $u = 0$ ; $v = 0$ ; $w = 0$ )	103E_1	103E_1	103E_1

新窗口打开

考虑到Chebyshev多项式良好的稳定性和收敛性,本文借助其来表述计算面上的场变量,从而将获取的二维偏微分方程转化为以场变量展开系数为未知量的线性代数方程组,避免了对超越方程的求解. 第 $k$ 子层第 $j$ 计算面上的各位移分量可延展为

$\begin{array}{l} u_{k}^{j} (α, β, t) = \sum_{n}^{=} ‍ \sum_{m}^{=} ‍ u_{k}^{jmn} (t) T_{m} (α) T_{n} (β) v_{k}^{j} (α, β, t) = \sum_{n}^{=} ‍ \sum_{m}^{=} ‍ v_{k}^{jmn} (t) T_{m} (α) T_{n} (β) w_{k}^{j} (α, β, t) = \sum_{n}^{=} ‍ \sum_{m}^{=} ‍ w_{k}^{jmn} (t) T_{m} (α) T_{n} (β) \end{array}\} (12)$

其中, $u_{k}^{jmn}$ , $v_{k}^{jmn}$ , $w_{k}^{jmn}$ 为展开系数. 非负整 $T_{n}$ 和 $T_{m}$ 分别为第 $n$ 和 $m$ 阶Chebyshev多项式. $M$ 和 $N$ 为截断量. 必须说明的是, $α$ 和 $β$ 的定义域为 $[- 1, + 1]$ ,在实际计算中需要进行坐标变换.

由式(10)和式(11)可知,子层 $k$ 的最终能量泛函可以写为

$Π_{k t otal} = T_{k} - U_{sk} - Π_{pk} (13)$

将计算面位移表达式(12)代入式 (13),经变分得到第 $k$ 子层的振动特征矩阵为

$(K_{k} - ω^{2} M_{k}) G_{k} = 0 (14)$

式中, $K_{k}$ 和 $M_{k}$ 为刚度和质量矩阵, $G_{k}$ 为系数向量.

层合结构层间界面需要满足位移连续和横向应力连续两个条件. 本文解法属于基于能量变分原理的位移解法,当所求得的解满足层间位移连续条件,层间横向应力连续条件将自动满足^[17]. 为了便于后续层与层之间特征矩阵的装配,在式(2)中,每一子层的上下表面分别选为该层的第一和最后一个计算面. 因此,式(14)中刚度/质量矩阵的第一行和最后一行分别为计算面位移场 $u_{k}^{1}$ , $v_{k}^{1}$ , $w_{k}^{1}$ 和 $u_{k}^{J_{k}}$ , $v_{k}^{J_{k}}$ , $w_{k}^{J_{k}}$ 的刚度/质量矩阵元素. 从物理角度来说,第 $k$ 层的第一个计算面与第 $k - 1$ 层的最后一个计算面为同一个计算面,第 $k$ 层的最后一个计算面与第 $k + 1$ 层的第一个计算面为同一个计算面. 因此,将所有子层的刚度/质量矩阵中共面的元素相叠加即可得到多层复合壳体的整体矩阵方程,组装方式如图2所示(质量矩阵和系数向量同理可得). 通过对整体矩阵方程进行特征值分解即可得到壳体振动特性信息.

显示原图| 下载原图ZIP| 生成PPT

图2 整体刚度矩阵组装方式

Fig. 2 Assemblage of the global stiffness matrix

2 多层复合壳体振动算例与讨论

在前面推导的基础上,文中通过多层复合板和圆柱壳算例来验证方法的可靠性. 下文采用字母F,S,C来表示自由、简支和固支边界条件.

2.1 多层复合板算例分析

表2给出了 $[0^{\circ} / 9 0^{\circ}]_{2}$ 和 $[0^{\circ} / 9 0^{\circ}]_{4}$ 多层复合板的无量纲基频 $Ω = ω h \sqrt[]{ρ / E_{2}}$ 与文献[18,19]的对比情况. 文献[18]采用的是状态空间与微分求积混合架构的半解析解法,其边界条件采用混合变量的方式直接描述,而文献[19]采用的是基于广义变分的位移解法. 复合板的材料和几何参数分别为： $E_{1} = 40 E_{2}$ , $G_{12} = G_{13} = 0.6 E_{2}$ , $G_{23} = 0.5 E_{2}$ , $ν_{12} = ν_{13} = 0.25$ , $ν_{23} = 0.49$ ; $a = b = 1$ , $h / a = 0.1$ . 计算面数 $J_{k} = 5$ ,面内截断量取 $M = N = 15$ . 可以看出,各种边界条件下,本文方法得到的解与文献[18,19]的数据对比吻合良好,最大偏差分别小于0.26%和0.25%, 从而说明了本文方法的正确性.

表2 $[0 ∘ / 9 0 ∘] 2$ 和 $[0 ∘ / 9 0 ∘] 4$ 复合板无量纲基频

Table 2 The fundamental frequency parameter $Ω$ of $[0 ∘ / 9 0 ∘] 2$ and $[0 ∘ / 9 0 ∘] 4$ laminated plates

Layout	Boundary	Frequency parameters $Ω$
		Ref.[18]	Ref.[19]	This paper
$[0 / 9 0^{\circ}]_{2}$	S-S-S-S	0.145 0	0.145 0	0.145 0
	C-S-S-S	0.165 5	0.165 5	0.165 5
	C-S-C-S	0.188 0	0.187 5	0.187 8
$[0 / 9 0^{\circ}]_{4}$	S-S-S-S	0.155 7	0.155 7	0.155 7
	C-S-S-S	0.177 9	0.177 8	0.177 8
	C-S-C-S	0.201 8	0.201 3	0.201 7

新窗口打开

表3给出了[0/90/0]复合板前四阶无量纲频率的收敛情况. 复合板的材料和几何参数与表2所用一致,面内截断量取 $M = N = 15$ . 可以看出,四边自由和四边简支边界条件下,计算面数取 $J_{k} = 4$ 与 $J_{k} = 6$ 的两组解相等,从而说明了本文方法收敛性良好.

表 3 [0/90/0]三层复合板无量纲频率收敛情况

Table 3 Convergency of frequency parameters $Ω$ of [0/90/0] laminated plates

Boundary	$J_{k}$	Frequency parameters $Ω$
		1	2	3	4
F-F-F-F	3	0.048 6	0.095 6	0.132 2	0.237 6
	4	0.048 6	0.095 2	0.131 9	0.235 1
	5	0.048 6	0.095 2	0.131 9	0.235 1
	6	0.048 6	0.095 2	0.131 9	0.235 1
S-S-S-S	3	0.147 1	0.218 3	0.243 3	0.243 3
	4	0.147 0	0.217 1	0.243 3	0.243 3
	5	0.147 0	0.217 1	0.243 3	0.243 3
	6	0.147 0	0.217 1	0.243 3	0.243 3

新窗口打开

表4给出的是四边简支软芯夹层方板( $h / a = 0.1$ )的无量纲频率参数与文献[20]的对比情况,其中, $m_{λ}$ 和 $n_{λ}$ 分别表示长度和宽度方向的振动波数. LW和ESL分别代表文献中基于分层理论(LW)和单层等效理论(ESL)给出的解. 该夹层板由复合表层与软芯中间层复合而成, 即[0°/90°/core/0°/90°],各层厚度比为1:1:20:1:1. 各层材料参数[ $E_{1}$ , $E_{2}$ , $E_{3}$ , $ν_{12}$ , $ν_{13}$ , $ν_{23}$ , $G_{12}$ , $G_{13}$ , $G_{23}$ , $ρ$ ]为：表层-[131 GPa, 10.34 GPa, 10.34 GPa, 0.22, 0.22, 0.49, 6.895 GPa, 6.205 GPa, 6.895 GPa, 1 627 kg/m3], 芯-[6.89 GPa, 6.89 GPa, 6.89 GPa, 0, 0, 0, 3.45 GPa, 3.45 GPa, 3.45 GPa, 97 kg/m3]. 计算面数取 $J_{k} = 5$ ,面内截断量取 $M = N = 15$ . 对比结果一方面说明了等效单层理论计算误差较大.

表 4 四边简支[0°/90 °/core/0°/90 °]软芯夹层板无量纲频率参数

Table 4 Frequency parameters Ω of [0°/90°/core/0°/90 °] sandwich plates of S-S-S-S restraint

$\begin{array}{l} M odes \\ (m_{λ}, n_{λ}) \end{array}$	$J_{f ace} - J_{c ore} - J_{f ace}$		Ref.[20]
	4-7-4	5-9-5	LW	ESL
, 1	1.849	1.849	1.848	4.962
, 2	3.222	3.222	3.220	8.193
, 2	4.292	4.292	4.289	10.30
, 3	5.228	5.227	5.223	11.99
, 3	6.099	6.099	6.094	13.75
, 3	7.683	7.683	7.676	16.45

新窗口打开

2.2 多层复合圆柱壳算例分析

表5给出了不同边界条件和周向波数下一软芯夹层圆柱壳结构的第一阶固有频率与文献[21]的对比情况. 其中, $n$ 表示周向波数. 该圆柱壳表层为金属铝,夹层的材料参数与表3所用的芯体材料相同. 计算面数取 $J_{k} = 5$ ,面内截断量取 $M = 15$ . 可以看出,本文方法具有良好的计算精度.

表 5 不同边界条件下软芯夹层圆柱壳第一阶频率

Table 5 The first frequencies of sandwich circular cylindrical shell for different boundary conditions

F-F		S-S		C-C
Pre.	3D^[21]	Pre.	3D^[21]	Pre.	3D^[21]
296.1	295.9	445.2	444.8	455.5	455.0
16.87	16.89	251.7	251.5	293.4	292.6
35.90	35.88	152.6	152.5	204.8	203.8
57.57	57.54	112.3	112.2	157.9	156.9
82.88	82.85	108.7	108.6	140.2	139.3
112.4	112.3	126.5	126.4	145.2	144.6
146.3	146.3	156.0	155.9	166.6	166.2
184.8	184.8	192.7	192.7	199.0	198.7

新窗口打开

3 结论

本文以独立铺层为建模对象,结合广义谱方法与微分求积技术建立了一种多层复合壳体三维分析新方法——谱--微分求积混合法. 该方法的特点为：(1)应用三维弹性理论对独立铺层进行精确建模;(2)引入微分求积技术对铺层进行数值离散,从而将三维偏微分问题转化为二维偏微分问题;(3)应用广义谱方法近似地表述离散计算面上的场变量,从而将获取的二维偏微分方程转化为以场变量谱展开系数为未知量的线性代数方程组. 与传统三维求解方法相比,谱--微分求积混合法在处理不同边界条件、铺层方式时具有很好的普适性和灵活性. 针对给定的不同边界条件和铺层方式,仅需要通过改变相关参数设置来完成边界约束的施加和铺层方式的配置而不需要逐一重新编程处理. 数值验证结果表明该方法收敛性好,计算精度高.

The authors have declared that no competing interests exist.

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]	陈建霖, 励争, 储鹏程. 大开口复合材料层合板强度破坏研究 . 力学学报, 2016, 48(7):1326-1333 [本文引用: 1] (Chen Jianlin, Li Zheng, Chu Pengcheng. Strength analysis of fiber reinforced composite laminates with big cutouts . Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(7): 1326-1333 (in Chinese)) [本文引用: 1]
[2]	马强, 周凤玺, 刘杰. 梯度波阻板的地基振动控制研究 . 力学学报,2017,49(7):1360-1369 [本文引用: 1] (Ma Qiang, Zhou Fengxi, Liu Jie. Analysis of ground vibration control by graded wave impeding block . Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(7): 1360-1369 (in Chinese)) [本文引用: 1]
[3]	郭洪宝, 王波, 贾普荣等. 平纹编织陶瓷基复合材料面内剪切细观损伤行为研究 . 力学学报, 2016, 48(3): 361-368 [本文引用: 1] (Guo Hongbao, Wang Bo, Jia Purong, et al. Mesoscopic damage behaviors of plain woven ceramic composite under in-plane shear loading . Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(3): 361-368 (in Chinese)) [本文引用: 1]
[4]	Reddy JN.Mechanics of Laminated Composite Plates and Shells: Theory and Analysis. Florida: CRC Press, 2003 [本文引用: 1]
[5]	Jin G, Ye T, Wang X, et al. A unified solution for the vibration analysis of FGM doubly-curved shells of revolution with arbitrary boundary conditions . Composites Part B Engineering, 2016, 89: 230-252 [本文引用: 1]
[6]	Ye T, Jin G. Elasticity solution for vibration of generally laminated beams by a modified Fourier expansion-based sampling surface method . Computers & Structures, 2016, 167: 115-130 [本文引用: 1]
[7]	瞿叶高, 华宏星, 谌勇等. 复合材料旋转壳自由振动分析的新方法 . 力学学报, 2013, 45(1): 139-143 [本文引用: 1] (Qu Yegao, Hua Hongxing, Chen Yong, et al. Vibroacoustic performance of simply supported honeycomb sandwich panels . Chinese Journal of Theoretical and Applied Mechanics, 2013, 45(1): 139-143 (in Chinese)) [本文引用: 1]
[8]	Carrera E. Historical review of zig-zag theories for multilayered plates and shells . Applied Mechanics Reviews, 2003, 56(3): 287-308 [本文引用: 1]
[9]	任树伟, 辛锋先, 卢天健. 蜂窝层芯夹层板结构振动与传声特性研究 . 力学学报, 2013, 45(3): 349-358 [本文引用: 1] (Ren Shuwei, Xin Fengxian, Lu Tianjian. Vibroacoustic performance of simply supported honeycomb sandwich panels . Chinese Journal of Theoretical and Applied Mechanics, 2013, 45(3): 349-358 (in Chinese)) [本文引用: 1]
[10]	刘人怀, 薛江红. 复合材料层合板壳非线性力学的研究进展 . 力学学报, 2017, 49(3): 487-506 [本文引用: 1] (Liu Renhuai, Xue Jianghong. Development of non-linear theories for laminated composite plates and shells . Chinese Journal of Theoretical and Applied Mechanics, 2013, 45(1): 139-143 (in Chinese)) [本文引用: 1]
[11]	Liu B, Xing Y. Exact solutions for free vibrations of orthotropic rectangular Mindlin plates . Composite Structures, 2011, 93: 1664-1672 [本文引用: 1]
[12]	吴振. 高性能整体--局部高阶理论及高阶层合板单元. [博士论文]. 大连:大连理工大学, 2007 [本文引用: 2] (Wu Zhen.Higher-performance global-local higher-order theory and higher-order laminated plate elements. [PhD Thesis]. Dalian: Dalian University of Technology, 2001 (in Chinese)) [本文引用: 2]
[13]	Lo SH, Zhen W, Cheung YK, et al. Hygrothermal effects on multilayered composite plates using a refined higher order theory . Composite Structures, 2010, 92(3): 633-646 [本文引用: 1]
[14]	吕朝锋. 基于状态空间架构的微分求积法及其应用. [博士论文]. 杭州:浙江大学, 2006 [本文引用: 1] (Lü Chaofeng.State-space-based differential quadrature method and its applications. [PhD Thesis]. Hangzhou: Zhejiang University, 2006 (in Chinese)) [本文引用: 1]
[15]	Saada AS.Elasticity: Theory and Applications, second ed. Florida: Ross Publishing Inc, 2009 [本文引用: 1]
[16]	Ye T, Jin G, Su Z. Three-dimensional vibration analysis of sandwich and multilayered plates with general ply stacking sequences by a spectral-sampling surface method . Composite Structures, 2017, 176: 1124-1142
[17]	Reddy JN.Energy Principles and Variational Methods in Applied Mechanic. New York: John Wiley & Sons, 1984 [本文引用: 2]
[18]	Chen WQ, Lü CF. 3D free vibration analysis of cross-ply laminated plates with one pair of opposite edges simply supported . Composite Structures, 2005, 69(1): 77-87 [本文引用: 3]
[19]	Qu Y, Wu S, Li H, et al. Three-dimensional free and transient vibration analysis of composite laminated and sandwich rectangular parallelepipeds: Beams, plates and solids . Composites Part B: Engineering, 2015, 73: 96-110 [本文引用: 3]
[20]	Rao MK, Desai YM. Analytical solutions for vibrations of laminated and sandwich plates using mixed theory . Composite Structures, 2004, 63(3): 361-373 [本文引用: 2]
[21]	Qu Y, Meng G. Dynamic analysis of composite laminated and sandwich hollow bodies of revolution based on three-dimensional elasticity theory . Composite Structures, 2014, 112: 378-396 [本文引用: 4]

大开口复合材料层合板强度破坏研究

2016

... 多层复合结构广泛应用于航空航天、高速列车、防卫装备等高科技行业,其振动特性分析一直是相关领域的研究焦点之一^[1,2,3]. 常用的多层复合结构振动建模理论主要有：经典等效单层理论、剪切变形理论、锯齿理论和三维层合理论四大类^{[4,5,6,7,8,9,10,11,12,13]}. 经典等效单层理论和剪切变形理论都是传统的板壳理论对多层复合结构的直接推广应用. 这两种理论先在厚度方向上把多层复合结构等效成各向异性的单一薄层,然后用结构中面的位移或者应力分量描述结构的整体振动行为. 由于组成材料的各向异性及呈层性等特点,多层复合结构受载后的应力状态和振动行为十分复杂,其振动位移沿厚度方向呈锯齿形分布且层间剪切和拉、压应力呈三维耦合状态. 因此,当结构厚度比较大或者各层材料差异较为明显时,位移连续但层间应力不连续的等效单层理论的建模误差急剧增大. 锯齿理论是以上述两类理论为基础,结合层间应力连续条件而得到的理论. 这类理论一般通过在结构的面内位移分量添加锯齿函数从而使横向剪切应力满足层间连续条件. 但该理论一般不满足上下表面自由边界条件也无法得到准确的层间应力^[12]. 直接应用三维层合理论可以满足位移和横向应力层间连续条件,这一理论和等效单层理论及锯齿理论最大的区别在于结构的每一个子层都看作是用独立的函数来描述的三维实体. 众所周知,只有边界条件和铺层方式简单的少数规则结构才能得到解析解,对其他较为复杂的边界条件、铺层方式以及结构形式,只能借助数值方法或者其他近似解法来求解^[14]. ...

大开口复合材料层合板强度破坏研究

2016

梯度波阻板的地基振动控制研究

2017

梯度波阻板的地基振动控制研究

2017

平纹编织陶瓷基复合材料面内剪切细观损伤行为研究

2016

平纹编织陶瓷基复合材料面内剪切细观损伤行为研究

2016

2003

A unified solution for the vibration analysis of FGM doubly-curved shells of revolution with arbitrary boundary conditions

2016

Elasticity solution for vibration of generally laminated beams by a modified Fourier expansion-based sampling surface method

2016

复合材料旋转壳自由振动分析的新方法

2013

复合材料旋转壳自由振动分析的新方法

2013

Historical review of zig-zag theories for multilayered plates and shells

2003

蜂窝层芯夹层板结构振动与传声特性研究

2013

蜂窝层芯夹层板结构振动与传声特性研究

2013

复合材料层合板壳非线性力学的研究进展

2013

复合材料层合板壳非线性力学的研究进展

2013

Exact solutions for free vibrations of orthotropic rectangular Mindlin plates

2011

2001

... [12]. 直接应用三维层合理论可以满足位移和横向应力层间连续条件,这一理论和等效单层理论及锯齿理论最大的区别在于结构的每一个子层都看作是用独立的函数来描述的三维实体. 众所周知,只有边界条件和铺层方式简单的少数规则结构才能得到解析解,对其他较为复杂的边界条件、铺层方式以及结构形式,只能借助数值方法或者其他近似解法来求解^[14]. ...

2001

Hygrothermal effects on multilayered composite plates using a refined higher order theory

2010

2006

2009

... 其中,

ε_{α}

ε_{β}

ε_{ζ}

分别表示主应变,

γ_{αβ}

γ_{αζ}

γ_{βζ}

为面内及横向剪切应变,

A_{α}

和

A_{β}

为Lamè常数^[15]. ...

Three-dimensional vibration analysis of sandwich and multilayered plates with general ply stacking sequences by a spectral-sampling surface method

2017

1984

... 借助哈密尔顿变分原理,式(10)可推导得到到结构的振动控制微分方程,具体形式参考文献[17]. ...

... 层合结构层间界面需要满足位移连续和横向应力连续两个条件. 本文解法属于基于能量变分原理的位移解法,当所求得的解满足层间位移连续条件,层间横向应力连续条件将自动满足^[17]. 为了便于后续层与层之间特征矩阵的装配,在式(2)中,每一子层的上下表面分别选为该层的第一和最后一个计算面. 因此,式(14)中刚度/质量矩阵的第一行和最后一行分别为计算面位移场

u_{k}^{1}

v_{k}^{1}

w_{k}^{1}

和

u_{k}^{J_{k}}

v_{k}^{J_{k}}

w_{k}^{J_{k}}

的刚度/质量矩阵元素. 从物理角度来说,第

k

层的第一个计算面与第

k - 1

层的最后一个计算面为同一个计算面,第

k

层的最后一个计算面与第

k + 1

层的第一个计算面为同一个计算面. 因此,将所有子层的刚度/质量矩阵中共面的元素相叠加即可得到多层复合壳体的整体矩阵方程,组装方式如图2所示(质量矩阵和系数向量同理可得). 通过对整体矩阵方程进行特征值分解即可得到壳体振动特性信息. ...

3D free vibration analysis of cross-ply laminated plates with one pair of opposite edges simply supported

2005

... 表2给出了

[0^{\circ} / 9 0^{\circ}]_{2}

和

[0^{\circ} / 9 0^{\circ}]_{4}

多层复合板的无量纲基频

Ω = ω h \sqrt[]{ρ / E_{2}}

与文献[18,19]的对比情况. 文献[18]采用的是状态空间与微分求积混合架构的半解析解法,其边界条件采用混合变量的方式直接描述,而文献[19]采用的是基于广义变分的位移解法. 复合板的材料和几何参数分别为：

E_{1} = 40 E_{2}

G_{12} = G_{13} = 0.6 E_{2}

G_{23} = 0.5 E_{2}

ν_{12} = ν_{13} = 0.25

ν_{23} = 0.49

;

a = b = 1

h / a = 0.1

. 计算面数

J_{k} = 5

,面内截断量取

M = N = 15

. 可以看出,各种边界条件下,本文方法得到的解与文献[18,19]的数据对比吻合良好,最大偏差分别小于0.26%和0.25%, 从而说明了本文方法的正确性. ...

... ]的对比情况. 文献[18]采用的是状态空间与微分求积混合架构的半解析解法,其边界条件采用混合变量的方式直接描述,而文献[19]采用的是基于广义变分的位移解法. 复合板的材料和几何参数分别为：

E_{1} = 40 E_{2}

G_{12} = G_{13} = 0.6 E_{2}

G_{23} = 0.5 E_{2}

ν_{12} = ν_{13} = 0.25

ν_{23} = 0.49

;

a = b = 1

h / a = 0.1

. 计算面数

J_{k} = 5

,面内截断量取

M = N = 15

. 可以看出,各种边界条件下,本文方法得到的解与文献[18,19]的数据对比吻合良好,最大偏差分别小于0.26%和0.25%, 从而说明了本文方法的正确性. ...

... . 可以看出,各种边界条件下,本文方法得到的解与文献[18,19]的数据对比吻合良好,最大偏差分别小于0.26%和0.25%, 从而说明了本文方法的正确性. ...

Three-dimensional free and transient vibration analysis of composite laminated and sandwich rectangular parallelepipeds: Beams, plates and solids

2015

... 表2给出了

[0^{\circ} / 9 0^{\circ}]_{2}

和

[0^{\circ} / 9 0^{\circ}]_{4}

多层复合板的无量纲基频

Ω = ω h \sqrt[]{ρ / E_{2}}

E_{1} = 40 E_{2}

G_{12} = G_{13} = 0.6 E_{2}

G_{23} = 0.5 E_{2}

ν_{12} = ν_{13} = 0.25

ν_{23} = 0.49

;

a = b = 1

h / a = 0.1

. 计算面数

J_{k} = 5

,面内截断量取

M = N = 15

. 可以看出,各种边界条件下,本文方法得到的解与文献[18,19]的数据对比吻合良好,最大偏差分别小于0.26%和0.25%, 从而说明了本文方法的正确性. ...

... ]采用的是状态空间与微分求积混合架构的半解析解法,其边界条件采用混合变量的方式直接描述,而文献[19]采用的是基于广义变分的位移解法. 复合板的材料和几何参数分别为：

E_{1} = 40 E_{2}

G_{12} = G_{13} = 0.6 E_{2}

G_{23} = 0.5 E_{2}

ν_{12} = ν_{13} = 0.25

ν_{23} = 0.49

;

a = b = 1

h / a = 0.1

. 计算面数

J_{k} = 5

,面内截断量取

M = N = 15

. 可以看出,各种边界条件下,本文方法得到的解与文献[18,19]的数据对比吻合良好,最大偏差分别小于0.26%和0.25%, 从而说明了本文方法的正确性. ...

... ,19]的数据对比吻合良好,最大偏差分别小于0.26%和0.25%, 从而说明了本文方法的正确性. ...

Analytical solutions for vibrations of laminated and sandwich plates using mixed theory

2004

... 式中,

σ_{α}

σ_{β}

σ_{ζ}

分别表示沿

α

β

ζ

方向上的正应力;

τ_{αβ}

τ_{αζ}

τ_{βζ}

为面内和横向剪切应力,

C_{k}

为材料主方向的材料弹性常数,转换矩阵

T_{k}

与材料铺层方式有关,一般为材料主坐标与结构主坐标夹角的函数^[20]. ...

... 表4给出的是四边简支软芯夹层方板(

h / a = 0.1

)的无量纲频率参数

与文献[20]的对比情况,其中,

m_{λ}

和

n_{λ}

分别表示长度和宽度方向的振动波数. LW和ESL分别代表文献中基于分层理论(LW)和单层等效理论(ESL)给出的解. 该夹层板由复合表层与软芯中间层复合而成, 即[0°/90°/core/0°/90°],各层厚度比为1:1:20:1:1. 各层材料参数[

E_{1}

E_{2}

E_{3}

ν_{12}

ν_{13}

ν_{23}

G_{12}

G_{13}

G_{23}

ρ

]为：表层-[131 GPa, 10.34 GPa, 10.34 GPa, 0.22, 0.22, 0.49, 6.895 GPa, 6.205 GPa, 6.895 GPa, 1 627 kg/m3], 芯-[6.89 GPa, 6.89 GPa, 6.89 GPa, 0, 0, 0, 3.45 GPa, 3.45 GPa, 3.45 GPa, 97 kg/m3]. 计算面数取

J_{k} = 5

,面内截断量取

M = N = 15

. 对比结果一方面说明了等效单层理论计算误差较大. ...

Dynamic analysis of composite laminated and sandwich hollow bodies of revolution based on three-dimensional elasticity theory

2014

... 表5给出了不同边界条件和周向波数下一软芯夹层圆柱壳结构的第一阶固有频率与文献[21]的对比情况. 其中,

n

表示周向波数. 该圆柱壳表层为金属铝,夹层的材料参数与表3所用的芯体材料相同. 计算面数取

J_{k} = 5

,面内截断量取

M = 15

. 可以看出,本文方法具有良好的计算精度. ...

... The first frequencies of sandwich circular cylindrical shell for different boundary conditions

F-F		S-S		C-C
Pre.	3D^[21]	Pre.	3D^[21]	Pre.	3D^[21]
296.1	295.9	445.2	444.8	455.5	455.0
16.87	16.89	251.7	251.5	293.4	292.6
35.90	35.88	152.6	152.5	204.8	203.8
57.57	57.54	112.3	112.2	157.9	156.9
82.88	82.85	108.7	108.6	140.2	139.3
112.4	112.3	126.5	126.4	145.2	144.6
146.3	146.3	156.0	155.9	166.6	166.2
184.8	184.8	192.7	192.7	199.0	198.7

3 结论

... [21] Pre. 3D^[21] 296.1 295.9 445.2 444.8 455.5 455.0 16.87 16.89 251.7 251.5 293.4 292.6 35.90 35.88 152.6 152.5 204.8 203.8 57.57 57.54 112.3 112.2 157.9 156.9 82.88 82.85 108.7 108.6 140.2 139.3 112.4 112.3 126.5 126.4 145.2 144.6 146.3 146.3 156.0 155.9 166.6 166.2 184.8 184.8 192.7 192.7 199.0 198.7 3 结论

... [21] 296.1 295.9 445.2 444.8 455.5 455.0 16.87 16.89 251.7 251.5 293.4 292.6 35.90 35.88 152.6 152.5 204.8 203.8 57.57 57.54 112.3 112.2 157.9 156.9 82.88 82.85 108.7 108.6 140.2 139.3 112.4 112.3 126.5 126.4 145.2 144.6 146.3 146.3 156.0 155.9 166.6 166.2 184.8 184.8 192.7 192.7 199.0 198.7 3 结论

〈

〉

多层复合壳体三维振动分析的谱--微分求积混合法

A SPECTRAL-DIFFERENTIAL QUADRATURE METHOD FOR 3-D VIBRATION ANALYSIS OF MULTILAYERED SHELLS

引言

1 多层复合壳体三维振动力学模型

1.1 分析模型

1.2 几何方程及物理方程

1.3 计算面选取及应变应变分量计算

1.4 子层特征矩阵及整体方程求取

2 多层复合壳体振动算例与讨论

2.1 多层复合板算例分析

2.2 多层复合圆柱壳算例分析

3 结 论

参考文献 文献选项 原文顺序 文献年度倒序 文中引用次数倒序 被引期刊影响因子

3 结论

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子