Sound Insulation of Porous Plates

Because of the still increasing noise pollution the numerical simulation of acoustic problems becomes more and more important. One essential aspect is the numerical treatment of noise insulation of solid walls. The main noise source is the bending vibration of separating components. In general, they consist of porous material, e.g., concrete or bricks. To take into account the porous structure as well as the damping effect of the porosity of these components a poroelastic plate theory is developed. The Finite Element implementation of this plate theory shows the importance of taking porous materials into account. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Poroelastodynamic Plate Formulation of Extendable Order

In order to model the sound insulation properties of porous plate-like structures, an adequate poroelastic plate theory is needed. The approach used here to obtain such a theory consists in expanding the primal variables into power series in thickness direction. This allows an analytical integration over the thickness, thus reducing the dimension of the problem from 3d to 2d. The resulting out-of-plane and in-plane problems are discussed. Numerical calculations show a good agreement between the proposed model and the three-dimensional system. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A symplectic analytical wave propagation model for damping and steady state forced vibration of orthotropic composite plate structure

An analytical wave propagation model is proposed in this paper for damping and steady state forced vibration of orthotropic composite plate structure by using the symplectic method. By solving an eigen-problem derived in the symplectic dual system of free bending vibration of orthotropic rectangular thin plates, the wave shape of plate is obtained in symplectic analytical form for any combination of simple boundary conditions along the plate edges. And then the specific damping capacity of wave mode is obtained symplectic analytically by using the strain energy theory. The steady state forced vibration of built-up plates structure is calculated by combining the wave propagation model and the finite element method. The vibration of the uniform plate domain of the built-up plates structure is described using symplectic analytical waves and the connector with discontinuous geometry or material is modeled using finite elements. In the numerical examples, the specific damping capacity of orthotropic rectangular thin plate with three different combinations of boundary condition is first calculated and analyzed. Comparisons of the present method results with respect to the results from the finite element method and from the Rayleigh–Ritz method validate the effectiveness of the present method. The relationship between the specific damping capacity of wave mode and that of modal mode is expounded. At last, the damped steady state forced vibration of a two plates system with a connector is calculated using the hybrid solution technique. The availability of the symplectic analytical wave propagation model is further validated by comparing the forced response from the present method with the results obtained using the finite element method.     

Modeling of amplitude-dependent damping characteristics of fiber reinforced composite thin plate

A damping model of fiber reinforced composite thin plate with consideration of amplitude-dependent property is established using the Jones–Nelson nonlinear theory in conjunction with the classical laminated plate theory, polynomial fitting method and strain energy method. In this model, the elastic moduli are expressed as the function of strain energy density and the loss factors in the longitudinal, transverse and shear directions are expressed as the functions of excitation amplitude. Moreover, three TC300 carbon/epoxy composite plates are taken as research objects to carry out a case study. One of them is used to determine the amplitude-dependent coefficients of loss factors in fiber reinforced composites by combining the least square method with polynomial fitting method, and the other two plates are used to verify the correctness of the theoretical model. The results of the developed model considering amplitude dependence and experimental test show a good consistency. It is discovered that the viscoelastic effect of epoxy resin materials will contribute to the increased degree of damping. So, if such plate structures exhibit more pronounced viscoelastic characteristics, there will be more significant in the nonlinear degree of amplitude-dependent damping phenomenon.     

3D nonlinear variable strain-rate-dependent-order fractional thermoviscoelastic dynamic stress investigation and vibration of thick transversely graded rotating annular plates/discs

In the present article, the idea of using the variable-order fractional-derivative thermoviscoelastic constitutive laws in dynamic stress and vibration analysis of the engineering structures, the required implementation backgrounds, and the relevant numerical solution procedures are investigated for the first time. In this regard, dynamic 3D stress and displacement fields and radial/transverse vibrations of transversely graded viscoelastic spinning thick plates/discs exposed to sudden thermoelastic loads are investigated. Instead of using the approximate plate theories, the exact thermoviscoelasticity theory is employed in the development of the governing equations. Since the variable fractional order is dependent on the localized deformation rates, the resulting thermoviscoelastic integro-differential equations are nonlinear. These equations are solved by utilizing a combination of the second-order backward/central/forward finite difference discretization of the spatial and time domains, numerical evaluation and updating of the Caputo-type fractional derivatives, updating the growing number of terms of the governing equations, and Picard's iterations. Various edge conditions are considered. Finally, comprehensive sensitivity analyses and various 3D plots are presented and discussed regarding the effects of the variable fractional order of the constitutive law, time variations of the nonuniformly distributed transverse loads, and edge conditions on the distributions and damping of the resulting displacement and stress components.     

不同理论下圆板特征值之间的解析关系

基于经典板理论(CPT)、一阶剪切变形板理论(FPT)以及Reddy三阶剪切变形板理论(RPT)之间,圆板轴对称特征值问题在数学上的相似性,研究了不同理论之间圆板特征值间的解析关系．将特征值问题的求解转化为代数方程的求解,并导出了不同理论之间圆板特征值的显式精确解析关系．从而,只要已知圆板特征值（临界屈曲载荷和固有频率）的经典结果,便很容易从这些解析关系中得到一阶和三阶理论下圆板特征值的相应结果,这便于工程应用,同时也可检验一阶和三阶理论下板特征值的数值结果的有效性、收敛性以及精确性等问题．     

不可压饱和多孔弹性杆动力响应的多辛方法

研究了不可压饱和多孔弹性杆的一维动力响应问题.基于多孔介质理论,在流相和固相微观不可压、固相骨架小变形的假定下,建立了不可压流体饱和多孔弹性杆一维轴向动力响应的数学模型.利用Hamilton空间体系的多辛理论,构造了不可压饱和多孔弹性杆轴向振动方程的多辛形式及其多种局部守恒律.采用中点Box离散方法得到轴向振动方程的多辛离散格式和局部能量守恒律以及局部动量守恒律的离散格式;数值模拟了不可压饱和多孔弹性杆的轴向振动过程,记录了每一时间步的局部能量数值误差和局部动量数值误差.结果表明,已构造的多辛离散格式具有很高的精确性和较长时间的数值稳定性,这为解决饱和多孔介质的动力响应问题提供了新的途径.     

Discrete mass conservation for porous media saturated flow

Global and local mass conservation for velocity fields associated with saturated porous media flow have long been recognized as integral components of any numerical scheme attempting to simulate these flows. In this work, we study finite element discretizations for saturated porous media flow that use Taylor–Hood (TH) and Scott–Vogelius (SV) finite elements. The governing equations are modified to include a stabilization term when using the TH elements, and we provide a theoretical result that shows convergence (with respect to the stabilization parameter) to pointwise mass‐conservative solutions. We also provide results using the SV approximation pair. These elements are pointwise divergence free, leading to optimal convergence rates and numerical solutions. We give numerical results to verify our theory and a comparison with standard mixed methods for saturated flow problems. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 625–640, 2014     

横观各向同性饱和弹性多孔介质非轴对称动力响应

应用Fourier展开和Hankel变换求解了简谐激励下横观各向同性饱和弹性多孔介质的非轴对称Biot波动方程,得到了一般解。用一般解给出了多孔介质总应力分量的表达式。最后对求解横观各向同性饱和弹性多孔介质非轴对称动力响应边值问题的方法作了系统说明,并且给出了数值分析特例。     

Formulação numérica de estruturas compósitas amortecidas utilizando as teorias da Deformação Cisalhante de Primeira Ordem e de Alta Ordem

Composite materials have been used in the design of the aircrafts structures because their low weight and high mechanical strength. However, structures made in composite material are exposed to dynamical and/or static loading environments. Therefore, a major research effort is undertaken in the development of tools numerical for analysis and design of composite structures. This paper presents a numerical formulation of the composite structures using the Finite Element Method (FEM). The damped composite structures, using inserted viscoelastic devices, and undamped composite structures are formulated by FEM. Viscoelastic materials are applied as continuous layers inserted on composite structures. The intrinsic damping of the composite material is included in the studies, too. The First‐order (FSDT) and Higher‐order Shear Deformation (HSDT) theories are formulated. They are distinguished by order of the approximation functions used in the mechanical displacements field. Both theories are computationally implemented using the Serendipity finite element. This is a rectangular finite element with 8 nodes, 5 or 11 degrees of freedom per node. The results are compared with papers predictions. The advantages and disadvantages of using each theory in the modeling of composite (thin or thick) and thick sandwiches structures, including the intrinsic and the viscoelastic damping, are discusses.     

An operator splitting approach for the interaction between a fluid and a multilayered poroelastic structure

We develop a loosely coupled fluid‐structure interaction finite element solver based on the Lie operator splitting scheme. The scheme is applied to the interaction between an incompressible, viscous, Newtonian fluid, and a multilayered structure, which consists of a thin elastic layer and a thick poroelastic material. The thin layer is modeled using the linearly elastic Koiter membrane model, while the thick poroelastic layer is modeled as a Biot system. We prove a conditional stability of the scheme and derive error estimates. Theoretical results are supported with numerical examples. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1054–1100, 2015     

A vibration analysis of composite laminated plates

A finite element formulation of the equations governing laminated anisotropic plates using Reddy's higher-order theory is presented. This simple higher-order shear deformable theory takes into account the parabolic distribution of the transverse shear deformation through the thickness of the plate and contains the same unknowns as in the first-order shear deformation theory. Finite element solutions are presented for rectangular plates of different layups, such as cross-ply, antisymmetric angle-ply, and sandwich plates with various material properties, boundaries, and plate aspect ratios. The numerical results are compared with the available closed-form results, the 3-D linear elasticity theory results, and the other available numerical results. A comparison is also made with test data from a laminated cantilever plate.     

On the application of genetic algorithms and gradient methods to the problem of the reconstruction of the initial stress field in a poroelastic inhomogeneous column

Initial stress fields play an important role in deformational processes in porous structures. Accurate modeling of dynamics behavior of the porous medium taking initial stresses into account requires an appropriate mathematical model. One of the common approaches for the construction of a mathematical model of the dynamic behavior of a medium taking the initial state into account is linearization. This paper presents equations of longitudinal vibrations of a poroelastic inhomogeneous body in the presence of the initial stress field. The influence of the initial stress rate and Biot’s modulus on the dynamic behavior of the environment is analyzed. The problem of reconstructing the initial stress field for the poroelastic inhomogeneous column based on the genetic algorithm and nonlinear optimization with the use of modern packages of numerical simulation is studied. A series of numerical experiments on the reconstruction of different initial stress field distribution laws using the described methods is carried out.     

中厚度复合材料夹芯层板变分渐近精细模型

为准确预测对中厚度复合材料夹芯层板分层开裂至关重要的沿厚向应力/应变分布,利用板固有小参数将原三维板分析严格拆分为沿厚向的一维分析和二维板非线性分析,并将原三维能量渐近扩展为系列二维近似能量泛函;通过对近似能量泛函中主导变分项(含翘曲项)的渐近修正,得到与原三维模型尽可能接近的近似能量,从而构建无需任何场变量假设的精细模型,并转换为工程常用的Reissner模型形式．通过4层复合材料夹芯板柱形弯曲算例表明:基于所构建模型重构的三维场精度较一阶剪切变形理论和经典层合理论更好,与精确解基本一致;由于所构建的变分渐近模型为等效单层板模型,在保证足够精度的前提下,相比三维有限元计算可减少2~3阶计算量,在精确性和有效性间取得较好的折衷．     

Convolution quadrature time-domain boundary element method for 2-D fluid-saturated porous media

For large-scale wave analyses of fluid-saturated porous media, a conventional time-domain boundary element method (BEM) cannot be applied because of the following reasons: (1) no time-domain fundamental solutions are known for some problems, (2) the method sometimes suffers from instability, and (3) the analyses require large amounts of computational time and memory. In this study, an innovative time-domain BEM is developed for a fluid-saturated porous medium. The formulation presented herein overcomes the above disadvantages using a convolution quadrature method (CQM), first proposed by Lubich, and hybrid-parallelization with both MPI and OpenMP. Problems involving the scattering of an incident plane wave by cavities in a 2-D poroelastic medium are solved as a means of validating the proposed method.     

Approximation of the vibration modes of a plate by Reissner-Mindlin equations

This paper deals with the approximation of the vibration modes of a plate modelled by the Reissner-Mindlin equations. It is well known that, in order to avoid locking, some kind of reduced integration or mixed interpolation has to be used when solving these equations by finite element methods. In particular, one of the most widely used procedures is the mixed interpolation tensorial components, based on the family of elements called MITC. We use the lowest order method of this family. Applying a general approximation theory for spectral problems, we obtain optimal order error estimates for the eigenvectors and the eigenvalues. Under mild assumptions, these estimates are valid with constants independent of the plate thickness. The optimal double order for the eigenvalues is derived from a corresponding -estimate for a load problem which is proven here. This optimal order -estimate is of interest in itself. Finally, we present several numerical examples showing the very good behavior of the numerical procedure even in some cases not covered by our theory.

     

Thermoelastic damping effect analysis in micro flexural resonator of atomic force microscopy

In the design of high-Q micro/nano-resonators, dissipation mechanisms may have damaging effects on the quality factor (Q). One of the major dissipation mechanisms is thermoelastic damping (TED) that needs an accurate consideration for prediction. Aim of this paper is to evaluate the effect of TED on the vibrations of thin beam resonators. In particular, we will focus on cantilever beam resonator used in atomic force microscopy (AFM). AFM resonator is actually a cantilever with a spring attached to its free end. The end spring is considered to capture the effect of surface stiffness between tip and sample surface. The coupled governing equations of motion of thin beam with consideration of TED effects are derived. In general, there are four elastic equations that are coupled with thermal conduction equation. Based on accurate assumptions, these equations are simplified and the various boundary conditions have been used in order to validate the computational procedure. In order to accurately determine TED effects, the coupled thermal conduction equation is solved for the temperature field by considering three-dimensional (3-D) heat conduction along the length, width and thickness of the beam. Weighted residual Galerkin technique is used to obtain frequency shift and the quality factor of the thin beam resonator. The obtained results for quality factor, frequency shift and sensitivity change due to thermo-elastic coupling are presented graphically. Furthermore, the effects of beam aspect ratio, stress-free temperature on the quality factor and the influence of the surface stiffness on the frequencies and modal sensitivity of the AFM cantilever with and without considering thermo-elastic damping effects are discussed.     

Higher-order theories for symmetric and unsymmetric fiber reinforced composite beams with C finite elements

A simple C⁰ isoparametric finite element formulation based on a set of higher-order displacement models for the analysis of symmetric and asymmetric multilayered composite and sandwich beams subjected to sinusoidal loading is presented. These theories do not require the usual shear correction coefficients which are generally associated with the Timoshenko theory. The four-noded Lagrangian cubic element with kinematic models having four, five and six degrees of freedom per node is used. A computer algorithm is developed which incorporates realistic prediction of transverse interlaminar stresses from equilibrium equations. By comparing the results obtained with the elasticity solution and the CPT (classical laminated plate theory) it is shown that the present higher-order theories give a much better approximation to the behaviour of laminated composite beams, both thick and thin. In addition numerical results for unsymmetric sandwich beams are presented which may serve as benchmark for future investigations.