Exact frequency response of two-node coupled bending-torsional beam element with attachments

This paper addresses the frequency response of coupled bending-torsional beams carrying an arbitrary number of in-span viscoelastic dampers and attached masses. Using the elementary coupled bending-torsion theory, along with appropriate generalized functions to treat the discontinuities of the response variables at the application points of dampers/masses, exact analytical expressions are derived for the frequency response of the beam under harmonically-varying, arbitrarily-placed point/polynomial loads. On this basis, the exact 6 × 6 dynamic stiffness matrix and 6 × 1 load vector of a two-node coupled bending-torsional beam finite element, with any number of in-span dampers/masses and harmonic loads, are obtained in a closed analytical form. Finally, the modal frequency response functions of the beam are built by a complex modal analysis approach, upon deriving pertinent orthogonality conditions for the modes. In this context, the modal impulse response functions are also obtained for time-domain analysis under arbitrary loads.     

Alternative Kelvin viscoelastic procedure for finite elements

In this work, an alternative Kelvin viscoelastic formulation for the finite element method (FEM) is described. This formulation performs spatial approximations before considering time integration and makes use of differential viscoelastic relations. A matrix time differential equation arises from the proposed formulation. It is solved numerically by a time marching procedure. It is shown that, after a small simplification, this methodology can be employed together with existent dynamic FEM packages. The methodology is extended to dynamic analysis leading to a rheological explanation for the first order modal decomposition stiffness proportional damping matrix. Plates and shells applications are shown in order to demonstrate the proposed formulation accuracy and stability.     

Damping sources interactions in impact of viscoelastic composite plates with damping treated SMA wires,using a hyperbolic plate theory

In the present paper, a comprehensive study is made on effects of the viscoelastic and the phase-transformation-based dissipations and their interactions on impact responses of viscoelastic composite plates with damping treated (structural hierarchy) shape memory alloy (SMA) wires, for the first time. In contrast to almost all of the available researches, a high-order hyperbolic plate theory that includes not only odd but also even functions of the transverse coordinate, is proposed and employed here. While a hierarchical viscoelastic constitutive law is employed for both the orthotropic and SMA materials, Brinson's constitutive law is refined to include the loading fluctuations and structural hierarchy of the SMA wire, simultaneously. The traditional Hertz and Yang-Sun contact laws are modified accordingly. The resulting highly nonlinear piecewise-defined integro-differential finite element governing equations are solved by an iterative algorithm within each time step. The presented discussions show that in contrast to the common belief, the zero-shear traction condition on the top and bottom surfaces of the viscoelastic orthotropic plate cannot be satisfied by the available plate theories, even for the symmetric lamination schemes. Results show that the viscoelasticity and phase-transformation effects on the resulting dynamic responses are more pronounced for the low and high energy impacts, respectively.     

The variability of dynamic responses of beams resting on elastic foundation subjected to vehicle with random system parameters

This paper investigates the variability of dynamic responses of a beam resting on an elastic foundation, which is subjected to a vehicle with uncertain parameters, such as random mass, stiffness, damping of the vehicle and random fields of mass density, and the elastic modulus of the beam and stiffness of elastic foundation. The vehicle is modeled as a two-degree-of-freedom spring-damper-mass system. The equations of motion of the beam was constructed using a finite element method. The mass and elastic properties of the beam, and the stiffness of foundation are assumed to be Gaussian random fields and were simulated by the spectral represent method. Masses, stiffness of the spring, and the damping coefficient of the vehicle are assumed as Gaussian random variables. The numerical analyses were performed using the finite element method (FEM) in conjunction with the Monte Carlo simulation (MCS). The variability of dynamic responses of the beam were investigated with various cases of random parameters. For each sample, the equations of motions were solved with the Wilson-q integral method to find dynamic responses. The influence of random system parameters and their correlation on the response variability is discussed in detail.     

A comparison between elastohydrodynamic and dry hertzian line contacts with regard to their nonlinear vibration behaviour

The nonlinear dynamic behaviour of a simple system employing two different types of contact is compared: elastohydrodynamic line contact and Hertzian line contact with appropriate damping. The elastohydrodynamic problem is solved with finite elements. For the Hertzian line contact problem an analytical solution is known. It is found that for highly loaded cases Hertzian and elastohydrodynamic contacts lead to very similar frequency response curves. This means that instead of solving the computational expensive elastohydrodynamic problem the Hertzian contact can be used. However, considering weakly loaded contacts the two types of contact show very different behaviour. In particular the natural frequencies are different and the elastohydrodynamic contact exhibits rich nonlinear dynamic behaviour including jump phenomena as well as sub- and superharmonic bifurcations. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling of viscoelastic contact-impact problems

An accurate finite element (FE) model for analyzing the response of viscoelastic structure under low-velocity impact is presented. Generalized standard linear solid (Wiechert) model is adopted to simulate the internal damping of the structure, because its capability of describing both creep and relaxation phenomena adequately. Newmark time integration scheme is proposed to transfer the problem into a static one for each time increment. The incremental convex programming method is modified to accommodate viscoelastic dynamic-contact problems. The Lagrange multiplier technique is selected to incorporate the contact condition. One, two and three-dimensional finite element model is presented to compare between the elastic and viscoelastic materials.     

非线性阻尼非线性刚度隔振系统随机动力学特性研究

针对随机激励环境,同时引入刚度和阻尼非线性来提高隔振系统的隔振性能.刚度和阻尼非线性分别是由水平弹簧和水平阻尼的几何布置获得.通过求解Fokker-Planck-Kolmogorov(FPK)方程等效非线性随机振动方程来研究非线性隔振系统在随机激励下的隔振性能,并使用路径积分和Monte-Carlo数值方法进行验证.在此基础上研究刚度非线性和阻尼非线性对隔振系统在随机激励下力传递率及其概率分布的影响.研究表明随着噪声强度的增加,非线性阻尼抑制振动的能力增强,但是在较小的随机激励下线性阻尼优于非线性阻尼.     

Analysis of a dynamic viscoelastic unilateral contact problem with normal damped response

In this paper we deal with a viscoelastic unilateral contact problem with normal damped response. The process is assumed to be dynamic and frictionless. Normal damping function is modeled by the Clarke subdifferential of a nonconvex and nonsmooth function. First, the variational formulation of this problem is provided in the form of a nonlinear first order variational–hemivariational inequality for the velocity field. Then, based on the surjectivity results for pseudomonotone and maximal monotone operators, we obtain the unique solvability for a new class of abstract evolutionary variational-hemivariational inequalities. Finally, we apply our abstract results to prove the existence of a unique weak solution to the corresponding contact problem.     

Optimal constrained layer damping with partial coverage

This paper deals with the optimal damping of beams constrained by viscoelastic layers when only one or several portions of the beam are covered. An efficient finite element model for dynamic analysis of such beams is used. The design variables are the dimensions and prescribed locations of the viscoelastic layers and the objective is the maximum viscoelastic damping factor. The method for non-linear programming in structural optimization is the so-called method of moving asymptotes.     

Numerical modeling of the active damping of forced thermomechanical resonance vibrations of viscoelastic shells of revolution with the help of piezoelectric inclusions

We consider the problem of active damping of forced resonance vibrations of viscoelastic shells of revolution with the help of piezoelectric sensors and actuators. Here, the interaction of electromechanical and thermal fields is taken into account. For modeling of vibrations, we use the Kirchhoff–Love hypotheses as well as hypotheses adequate to them and describing the distribution of temperature and electric field quantities. The shell temperature increases as a result of dissipative heating. For the active damping of vibrations, piezoelectric sensors and actuators are used. It is supposed that the electromechanical characteristics of materials depend on the temperature. The solution of this complex nonlinear problem has been obtained by the iterative method and finite element method. We have investigated the influence of temperature of dissipative heating on the efficiency of active damping of vibrations of a viscoelastic cylindrical panel with rigid restraint of its edges.     

Hybrid modeling approach for vehicle frame coupled with nonlinear dampers

The vehicle frame system comprises frame structure and nonlinear dampers. In order to investigate the effects of frame flexibility and nonlinear hysteresis, a hybrid modeling approach for vehicle frame coupled with nonlinear dampers will be proposed. Before that, a complex model for nonlinear damper is developed consisting of knowledge-based model and support vector machine (SVM) model. The frame structure is modeled by FEM where the SVM complex model of damper is embedded in. Thus a hybrid model for vehicle frame system is established and successfully validated via a dummy vehicle riding in different conditions. The results show that the hybrid model can capture the nonlinear dynamic characteristics accurately. The hybrid model can also provide a basis for structural design with the existing of FEM model.     

Numerical Analysis of a Dynamic Contact Problem with History-Dependent Operators

In this paper, we study a dynamic contact model with long memory which allows both the convex potential and nonconvex superpotentials to depend on history-dependent operators. The deformable body consists of a viscoelastic material with long memory and the process is assumed to be dynamic. The contact involves a nonmonotone Clarke subdifferential boundary condition and the friction is modeled by a version of the Coulomb's law of dry friction with the friction bound depending on the total slip. We introduce and study a fully discrete scheme of the problem, and derive error estimates for numerical solutions. Under appropriate solution regularity assumptions, an optimal order error estimate is derived for the linear finite element method. This theoretical result is illustrated numerically.     

A three-step Oseen-linearized finite element method for incompressible flows with damping

This paper introduces a three-step Oseen-linearized finite element method for the 2D/3D steady incompressible Navier–Stokes equations with nonlinear damping term. Within this method, we solve a nonlinear problem over a coarse grid followed by solving two Oseen-linearized problems over a fine grid, which possess the same stiffness matrices with only various right-hand sides. We theoretically analyze the stability of the present method, and derive optimal error estimates of the finite element solutions. We conduct a series of numerical experiments which support the theoretical analysis and test the effectiveness of the proposed method. We demonstrate numerically that there is a significant improvement in the accuracy of the approximate solutions over those for the standard two-level method.     

On a system of nonlinear wave equations with Balakrishnan–Taylor damping

In this paper, we study the initial-boundary value problem for a coupled system of nonlinear viscoelastic wave equations of Kirchhoff type with Balakrishnan–Taylor damping terms. For certain class of relaxation functions and certain initial data, we prove that the decay rate of the solution energy is similar to that of relaxation functions which is not necessarily of exponential or polynomial type. Also, we show that nonlinear source of polynomial type is able to force solutions to blow up in finite time even in presence of stronger damping.     

频域有限元计算的扩展面向目标误差估计

研究了针对频域有限元直接动态分析的面向目标误差估计以及误差范围估计计算方法．面向目标的误差估计方法就是专门针对如何准确和经济地估算特定值误差的一种方法,利用原问题的共轭偶问题进行计算．频域有限元的直接动态分析是模拟频域扫描实验的一种计算方法,专门针对谐振激励的线性动态响应问题,利用将原自由度分解为实部和虚部描述频率的变化,从而计算变形体的动态响应．利用扩展针对有限元的面向目标误差估计的自由度,将该方法应用到直接动态分析中进行误差估计．通过建立同时包含实部和虚部自由度的能量弱形式及偶问题,并将其数值实现,估算频域直接动态分析有限元解的误差及误差范围,并通过悬臂梁的激振算例进行了验证．     

A posteriori error analysis for the normal compliance problem

In this work, a contact problem between a linear elastic material and a deformable obstacle is numerically analyzed. The contact is modeled using the well-known normal compliance contact condition. The weak formulation leads to a nonlinear variational equation which is approximated by using the finite element method. A priori error estimates are recalled. Then, we define an a posteriori error estimator of residual type to evaluate the accuracy of the finite element approximation of the problem. Upper and lower bounds of the discretization error are proved for this estimator.     

Application of fuzzy sets to transient analysis of space structures

An application of fuzzy sets, in conjunction with finite elements, to the transient analysis of a precision-deployable space structure is presented. The structural members are modeled by using beam finite elements, and the structure's latch joint is modeled by using a spring–damper–Coulomb friction element. Two types of transient response simulations are performed: slow transient load–deflection response and transient impulse response. The first simulation is used to evaluate the stiffness and buckling loads at the structure's tip. The second simulation is used to evaluate the structure's natural frequencies, mode shapes and the precision of the final shape. For each simulation the possibility distributions of various response quantities are obtained. Fuzzy sets are used to represent three beam properties, namely: damping coefficient, bending stiffness, and axial stiffness; as well as two joint parameters: Coulomb friction force and damping coefficient. Fuzzy set techniques provide an insight into the range of possible responses associated with the combined selected variations in the system parameters.     

Vibrational analysis of electrostatically actuated microstructures considering nonlinear effects

The vibrational behavior of electrostatically actuated microstructures subjected to nonlinear squeeze film damping and in-plane forces is investigated. First-Order Shear Deformation Theory is used to model dynamical system by means of finite element method, while finite difference method is applied to solve the nonlinear Reynolds equation of squeeze film damping simultaneously. Vibrational analysis of microplates is performed by solving eigenvalue problem, after validating the model by pull-in phenomenon and transient behavior. In addition, considering nonlinear squeeze film damping and step-input actuations, response frequencies of microplates are calculated. Effect of ambient pressure and in-plane forces on dynamic pull-in phenomenon is also studied. Results for simplified models are verified and are in good agreement with the published literature. This investigation can reveal nonlinear vibrational behavior of microstructures.     

The Vibrational Behavior of Bladed Disks in Consideration of Friction Damping and Contact Elasticity

Rotating turbine blading is subjected to fluctuating gas forces during operation that cause blade vibrations. One of the main tasks in the design of turbomachinery blading is the reduction of the vibration amplitudes of the blades to avoid high resonance stresses that could damage the blading. The vibration amplitudes of the blades can be reduced significantly to a reasonable amount by means of friction damping devices such as underplatform dampers. In the case of blade vibrations, relative displacements between the friction damping devices and the neighboring blades occur and friction forces are generated that provide additional damping to the structure due to the dry friction energy dissipation. In real turbomachinery applications, spatial blade vibrations caused by a complex blade geometry and distributed excitation forces acting on the airfoil accur. Therefore, a three dimensional model including an appropriate spatial contact model to predict the generalized contact forces is necessary to describe the vibrational behavior of the blading with sufficient accuracy, see [1] and [2]. In this paper the contact model presented in [2] is extended to include also local deformations in the contacts between underplatform dampers and the contact surfaces of the adjacent blades. The additional elasticity in the contact influences the resonance frequency of the coupled bladed disk assembly. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)