Vibration analysis of linear coupled thermoviscoelastic thin plates by a variational approach

According to the integral type constitutive relation of linear coupled thermoviscoelasticity, a mathematical model of thin plates is set up by the introduction of “structural functions” and “thermal functions” in the sense of the Kirchhoff’s hypothesis. The corresponding integral type variational formulations are presented by means of modern convolution bilinear forms as well as classical Cartesian bilinear forms. The Ritz method in the spatial domain and the differentiating method in the temporal domain are used to approximate the mathematical model in a system of rectangular Cartesian coordinates. By properties of inequality and parabola, the structure of dynamic solution to vibration of a thermoviscoelastic thin plate under a harmonic thermal load is studied in the space splayed by material parameter and loading parameter. The influences of thermal excitation frequency, mechanical relaxation time and thermal relaxation time on amplitude and phase difference of steady-state vibration of a square plate are investigated by amplitude-frequency analysis and phase-frequency analysis. Double-peak resonance vibration of thermoviscoelastic plates exists for given parameters.     

QUASI-STATIC ANALYSIS FOR VISCOELASTIC TIMOSHENKO BEAMS WITH DAMAGE

Based on convolution-type constitutive equations for linear viscoelastic materials with damage and the hypotheses of Timoshenko beams, the equations governing quasi-static and dynamical behavior of Timoshenko beams with damage were first derived. The quasi-static behavior of the viscoelastic Timoshenko beam under step loading was analyzed and the analytical solution was obtained in the Laplace transformation domain. The deflection and damage curves at different time were obtained by using the numerical inverse transform and the influences of material parameters on the quasi-static behavior of the beam were investigated in detail.     

热黏弹性梁拟静态弯曲问题的解析分析

基于平面应力假设和热黏弹性材料的积分型本构关系,建立了以位移分量为未知量的热黏弹性梁静动力学分析的二维数学模型。针对拟静态弯曲问题,首先,在Laplace变换域,引入位移势函数,将控制方程解耦;其次,根据给定的平面温度场和边界条件,采用分离变量法,引入热应力函数,得到了热黏弹性梁的热应力分布;最后,利用Laplace逆变换,获得了热黏弹性梁拟静态弯曲热应力响应的解析解,考察了热载荷作用下几何、黏弹性等参数对梁应力和位移的影响。     

Transient analysis of viscoelastic helical rods subject to time-dependent loads

The transient analysis of viscoelastic helical rods subject to time-dependent loads are examined in the Laplace domain. The governing equations for naturally twisted and curved spatial rods obtained using the Timoshenko beam theory are rewritten for cylindrical helical rods. The curvature of the rod axis, effect of rotary inertia and, shear and axial deformations are considered in the formulation. The material of the rod is assumed to be homogeneous, isotropic and linear viscoelastic. The viscoelastic constitutive equations are written in the Boltzmann–Volterra form. Ordinary differential equations in canonical form obtained in the Laplace domain are solved numerically using the complementary functions method to calculate the dynamic stiffness matrix of the problem. The solutions obtained are transformed to the real space using an appropriate numerical inverse Laplace transform method. Numerical results for quasi-static and dynamic response of viscoelastic models are presented in the form of graphics.     

Pull-in analysis of electrically actuated viscoelastic microbeams based on a modified couple stress theory

A new beam model is developed for the viscoelastic microbeam based on a modified couple stress model which contains only one material length scale parameter. The governing equations of equilibrium together with initial conditions and boundary conditions are obtained by a combination of the basic equations of modified couple stress theory and Hamilton’s principle. This new beam model is then used for an electrically actuated microbeam-based MEMS structure. The dynamic and quasi-static governing equations of an electrically actuated viscoelastic microbeam are firstly given where the axial force created by the midplane stretching effect is also considered. Galerkin method is used to solve above equation and this method is also validated by the finite element method (FEM) when our model is reduced into an elastic case. The numerical results show that the instantaneous pull-in voltage, durable pull-in voltage and pull-in delay time predicted by this newly developed model is larger (longer) than that predicted by the classical beam model. A comparison between the quasi-static model results and the dynamic model results is also given.     

Nonlinear dynamic stability analysis of Euler–Bernoulli beam–columns with damping effects under thermal environment

In this study, a unified nonlinear dynamic buckling analysis for Euler–Bernoulli beam–columns subjected to constant loading rates is proposed with the incorporation of mercurial damping effects under thermal environment. Two generalized methods are developed which are competent to incorporate various beam geometries, material properties, boundary conditions, compression rates, and especially, the damping and thermal effects. The Galerkin–Force method is developed by implementing Galerkin method into force equilibrium equations. Then for solving differential equations, different buckled shape functions were introduced into force equilibrium equations in nonlinear dynamic buckling analysis. On the other hand, regarding the developed energy method, the governing partial differential equation for dynamic buckling of beams is also derived by meticulously implementing Hamilton’s principles into Lagrange’s equations. Consequently, the dynamic buckling analysis with damping effects under thermal environment can be adequately formulated as ordinary differential equations. The validity and accuracy of the results obtained by the two proposed methods are rigorously verified by the finite element method. Furthermore, comprehensive investigations on the structural dynamic buckling behavior in the presence of damping effects under thermal environment are conducted.     

功能梯度变曲率曲梁的几何非线性模型及其数值解

基于弹性曲梁平面问题的精确几何非线性理论,建立了功能梯度变曲率曲梁在机械和热载荷共同作用下的无量纲控制方程和边界条件,其中基本未知量均被表示为变形前的轴线坐标的函数。以椭圆弧曲梁为例,采用打靶法求解非线性常微分方程的两点边值问题,获得了两端固定功能梯度椭圆弧曲梁在横向非均匀升温下的热弯曲变形数值解,分析了材料梯度指数、温度参数、结构几何参数等对曲梁受力及变形的影响。     

Differential quadrature method for analyzing nonlinear dynamic characteristics of viscoelastic plates with shear effects

The integro-partial differential equations governing the dynamic behavior of viscoelastic plates taking account of higher-order shear effects and finite deformations are presented. From the matrix formulas of differential quadrature, the special matrix product and the domain decoupled technique presented in this work, the nonlinear governing equations are converted into an explicit matrix form in the spatial domain. The dynamic behaviors of viscoelastic plates are numerically analyzed by introducing new variables in the time domain. The methods in nonlinear dynamics are synthetically applied to reveal plenty and complex dynamical phenomena of viscoelastic plates. The numerical convergence and comparison studies are carried out to validate the present solutions. At the same time, the influences of load and material parameters on dynamic behaviors are investigated. One can see that the system will enter into the chaotic state with a paroxysm form or quasi-periodic bifurcation with changing of parameters.     

Stability of curvilinear Euler-Bernoulli beams in temperature fields

In this work, stability of thin flexible Bernoulli-Euler beams is investigated taking into account the geometric non-linearity as well as a type and intensity of the temperature field. The applied temperature field T(x,z) is yielded by a solution to the 2D Laplace equation solved for five kinds of thermal boundary conditions, and there are no restrictions put on the temperature distribution along the beam thickness. Action of the temperature field on the beam dynamics is studied with the help of the Duhamel theory, whereas the motion of the beam subjected to the thermal load is yielded employing the variational principles.The heat transfer (Laplace equation) is solved with the use of the finite difference method (FDM) of the third-order accuracy, while the integrals along the beam thickness defining the thermal stress and moments are computed using Simpson's method. Partial differential equations governing the beam motion are reduced to the Cauchy problem by means of application of FDM of the second-order accuracy. The obtained ordinary differential equations are solved with the use of the fourth-order Runge-Kutta method.The problem of numerical results convergence versus a number of beam partitions is investigated. A static solution for a flexible Bernoulli-Euler beam is obtained using the dynamic approach based on employment of the relaxation/set-up method.Novel stability loss phenomena of a beam under the thermal field are reported for different beam geometric parameters, boundary conditions, and the temperature intensity. In particular, it has been shown that stability of the flexible beam during heating the beam surface essentially depends on the beam thickness.     

Dynamic response to a moving load of a Timoshenko beam resting on a nonlinear viscoelastic foundation

The present paper investigates the dynamic response of finite Timoshenko beams resting on a sixparameter foundation subjected to a moving load. It is for the first time that the Galerkin method and its convergence are studied for the response of a Timoshenko beam supported by a nonlinear foundation. The nonlinear Pasternak foundation is assumed to be cubic. Therefore, the efects of the shear deformable beams and the shear deformation of foundations are considered at the same time. The Galerkin method is utilized for discretizing the nonlinear partial differential governing equations of the forced vibration. The dynamic responses of Timoshenko beams are determined via the fourth-order Runge–Kutta method. Moreover, the efects of diferent truncation terms on the dynamic responses of a Timoshenko beam resting on a complex foundation are discussed. The numerical investigations shows that the dynamic response of Timoshenko beams supported by elastic foundations needs super high-order modes. Furthermore, the system parameters are compared to determine the dependence of the convergences of the Galerkin method.     

Quasi-static response of linear viscoelastic cantilever beams subject to a concentrated harmonic end load

This study evaluates the response of a uniform cantilever beam with a symmetric cross-section fixed at one end, and submitted to a lateral concentrated sinusoidal load at the free extremity. The beam material is assumed to be homogeneous, isotropic and linear viscoelastic. Due to the nature of the loading and the beam slenderness, large displacements are developed but the strains are considered small. Consequently, the mathematical formulation only involves geometrical non-linearity. It is also assumed that the beam is inextensible (neutral axis length is constant) and that inertial forces are negligible, i.e., dynamic effects are insignificant and the system can thus be modeled quasi-statically. The beam is therefore subject to oscillations caused by the sinusoidal time-dependent load, leading to a transient response until the material stabilizes and the system exhibits a periodic response, which can be conveniently described in the frequency domain. The time domain solution of this problem is elaborated by considering the quasi-static response for each time interval. The mathematical equations are presented in dimensional and dimensionless forms, and for the latter case, a numerical solution is generated and several case studies are presented. The problem is governed by a set of non-linear ordinary differential equations encompassing functions of space and time that relate the curvature, rotation angle, bending moment and geometrical coordinates. In this study, an elegant solution is deduced using perturbation theory, yielding a precise steady-state solution in the frequency domain with considerable computational economy. The solutions for both time and frequency domain methods are developed and compared using a case study for a series of dimensionless parameters that influence the response of the system.     

Effects of the dynamic wheel–rail interaction on the ground vibration generated by a moving train

Based on Biot’s fully dynamic poroelastic theory, the dynamic responses of the poroelastic half-space soil medium due to quasi-static and dynamic loads from a moving train are investigated semi-analytically. The dynamic loads are assumed to be generated from the rail surface irregularities. The vehicle dynamics model is used to simulate the axle loads (quasi-static loads) and the dynamic loads from a moving train. The compatibility of the displacements at wheel–rail contact points couple the vehicle and the track–ground subsystem, and yield equations for the dynamic wheel–rail loads. A linearized Hertzian contact spring between the wheel and rail is introduced to calculate the dynamic loads. Using the Fourier transform, the governing equations for the poroelastic half-space are then solved in the frequency–wavenumber domain. The time domain responses are evaluated by the fast inverse Fourier transform. Numerical results show that the dynamic loads can make important contribution to dynamic response of the poroelastic half-space for different train speed, and the dynamically induced responses lie in a higher frequency range. The ground vibrations caused by the moving train can be intensified as the primary suspension stiffness of the vehicle increases.     

Nonlinear vibrations of rotating thin circular cylindrical shell

Nonlinear vibrations of thin circular cylindrical shells are investigated in this paper. Based on Love thin shell theory, the governing partial differential equations of motion for the rotating circular cylindrical shell are formulated using Hamilton principle. Taking into account the clamped-free boundary conditions, the partial differential system is truncated by using the Galerkin method. Sequentially, the effects of temperature, geometric parameters, circumferential wave number, axial half wave number and rotating speed on the nature frequency of the rotating circular cylindrical shell are studied. The dynamic responses of the rotating circular cylindrical shell are also investigated in time domain and frequency domain. Then, the effects of nonlinearity, excitation and damping on frequency responses of steady solution are investigated.     

Thermoviscoelastic deformations of functionally graded thin plates

By considering the influence of the mid-plane strains due to the inhomogeneous property of FGMs across the thickness, a constitutive equation of thermoviscoelastic functionally graded thin plates is reduced on the basis of the Kirchhoff's hypothesis for the classical plate theory. The corresponding simplified Gurtin's type variational principle of functionally graded thin plates is presented by means of the modern convolution bilinear forms. By using the Navier analytic method or combining the Ritz method in the spatial domain and the Legendre interpolation method in the temporal domain, the static thermoelastic deformations or quasi-static thermoviscoelastic deformations of functionally graded thin plates under mechanical or thermal loads are studied.     

Aerothermoelastic behavior of supersonic rotating thin-walled beams made of functionally graded materials

In this study, a thin-walled beam made of functionally graded material (FGM) which is used as rotating blades in turbomachinery under aerothermoelastic loading is investigated. The governing equations, which are based on first-order shear deformation theory, include the effects of the presetting angle, the secondary warping, temperature gradient through the wall thickness of the beam and also the rotational speed. Moreover, quasi-steady aerodynamic pressure loadings are determined using first-order piston theory, and steady beam surface temperature is obtained from gas dynamics theory. Then, the blade partial differential equations are transformed into a set of ordinary differential equations using the extended Galerkin method. Finally, having solved the resulting structural–fluid–thermal eigenvalue system of equations, the effects of Mach number and geometric parameters on natural frequencies are presented. The results demonstrate that the natural frequencies decrease under aerothermoelastic loading at high Mach numbers.     

Dynamic response of an infinite Timoshenko beam on a nonlinear viscoelastic foundation to a moving load

The present paper investigates the dynamic response of infinite Timoshenko beams supported by nonlinear viscoelastic foundations subjected to a moving concentrated force. Nonlinear foundation is assumed to be cubic. The nonlinear governing equations of motion are developed by considering the effects of the shear deformable beams and the shear modulus of foundations at the same time. The differential equations are, respectively, solved using the Adomian decomposition method and a perturbation method in conjunction with complex Fourier transformation. An approximate closed form solution is derived in an integral form based on the presented Green function and the theorem of residues, which is used for the calculation of the integral. The dynamic response distribution along the length of the beam is obtained from the closed form solution. The derivation process demonstrates that two methods for the dynamic response of infinite beams on nonlinear foundations with a moving force give the consistent result. The numerical results investigate the influences of the shear deformable beam and the shear modulus of foundations on dynamic responses. Moreover, the influences on the dynamic response are numerically studied for nonlinearity, viscoelasticity and other system parameters.     

基于微分求积法的钢筋混凝土梁静力分析

采用广义微分求积法(GDQM)对钢筋混凝土(RC)梁进行了准静力分析,得到了其抗静载的强度特性．首先,基于虚功原理导出了考虑钢筋和混凝土材料的非线性的RC梁准静力分析控制微分方程,并根据广义微分求积法对其离散,从而得到有限自由度的非线性代数方程组,进而采用Newton-Raphson迭代求解格式,建立了荷载增量法数值分析模型．其次,通过本文GDQM与有限元法分析结果比较,表明了新建算法的正确性;与有限元法的收敛性对比表明本文算法较有限元法有优越性．     

Dynamic Analysis of a Composite Hollow Sphere Composed of Elastic and Piezoelectric Layers

Dynamic analysis of a two-layered elasto-piezoelectric composite hollow sphere under spherically symmetric deformation is developed. An unknown function of time is first introduced in terms of the charge equation of electrostatics and then the governing equations of piezoelectric layer, in which the unknown function of time is involved, are derived. By the method of superposition, the dynamic solution for elastic and piezoelectric layers is divided into quasi-static and dynamic parts. The quasi-static part is treated independently by the state space method and the dynamic part is obtained by the separation of variables method. By virtue of the obtained quasi-static and dynamic parts, a Volterra integral equation of the second kind with respect to the unknown function of time is derived by using the electric boundary conditions for piezoelectric layer. Interpolation method is employed to solve the integral equation efficiently. The transient responses for elastic and electric fields are finally determined. Numerical results are presented and discussed.     

热载荷作用下梯度多孔材料梁弯曲及过屈曲力学行为的研究

基于Bernoulli-Euler梁理论,引入物理中面解耦了复合材料结构的面内变形与横向弯曲特性,研究了梯度多孔材料矩形截面梁在热载荷作用下的弯曲及过屈曲力学行为．假设沿梁厚度方向材料的性质是连续变化的,利用能量法推导了矩形截面梁的控制微分方程和边界条件,并用打靶法对无量纲化的控制方程进行数值求解．利用计算得到的结果分析了材料的性质、热载荷、边界条件对矩形截面梁非线性力学行为的影响．结果表明,对称材料模型下,固支梁与简支梁均显示出了典型的分支屈曲行为特征,而其临界屈曲热载荷值均会随着孔隙率系数的增加而单调增加．非对称材料模型下,固支梁仍显示出分支屈曲行为特征,但其临界屈曲热载荷不再随着孔隙率系数的变化而单调变化;而对于两端简支梁,发生了弯曲变形,弯曲挠度随载荷的增大而增大．     

A mixed modal-differential quadrature method for free and forced vibration of beams in contact with fluid

A simple and accurate mixed modal-differential quadrature formulation is proposed to study the dynamic behavior of beams in contact with fluid. Both free and forced vibration problems are considered. The proposed mixed methodology uses the modal technique for the structural domain while it applies the differential quadrature method (DQM) to the fluid domain. Thus, the governing partial differential equations of the beam and fluid are reduced to a set of ordinary differential equations in time. In the case of forced vibration, the Newmark time integration scheme is employed to solve the resulting system of ordinary differential equations. The proposed formulation, in general, combines the simplicity of the modal method and high accuracy and efficiency of the DQM. Its application is shown by solving some beam-fluid interaction problems. Comparisons with analytical solutions show that the present method is very accurate and reliable. To demonstrate its efficiency, the test problems are also solved using the finite element method (FEM). It is found that the proposed method can produce better accuracy than the FEM using less computational time. The technique presented in this investigation is general and can be used to solve various fluid-structure interaction problems.