Analysis of nonlinear vibration for symmetric angle-ply laminated viscoelastic plates with damage

The behavior of nonlinear vibration for symmetric angle-ply laminated plates including the material viscoelasticity and damage evolution is investigated. By employing the von Karman's nonlinear theory, strain energy equivalence principle and Boltzmann superposition principle, a set of governing equations of nonlinear integro-differential type are derived. By applying the finite difference method, Newmark method and iterative procedure, the governing equations are solved. The effects of loading amplitudes, exciting frequencies and different ply orientations on the critical time to failure initiation and nonlinear vibration amplitudes of the structures are discussed. Numerical results are presented for the different parameters and compared with the available data. The project supported by the National Natural Science Foundation of China (10272042) and the Special Science Fund of the Doctoral Discipline of the Ministry of Education, China (20020532018) The English text was polished by Ron Marshall.     

Four-point combined DE/FE algorithm for brittle fracture analysis of laminated glass

A four-point combined DE/FE algorithm is proposed to constrain the rotation of a discrete element about its linked point and analyze the cracks propagation of laminated glass. In this approach, four linked points on a discrete element are combined with four nodes of the corresponding surface of a finite element. The penalty method is implemented to calculate the interface force between the two subdomains, the finite element (FE) and the discrete element (DE) subdomains. The sequential procedure of brittle fracture is described by an extrinsic cohesive fracture model only in the DE subdomain. An averaged stress tensor for granular media, which is automatically symmetrical and invariant by translations, is used to an accurate calculation of the averaged stress of the DE. Two simple cases in the elastic range are given to certify the effectiveness of the combined algorithm and the averaged stress tensor by comparing with the finite element method and the mesh-size dependency of the combined algorithm and the cohesive model is also investigated. Finally, the impact fracture behavior of a laminated glass beam is simulated, and the cracks propagation is compared with experimental results showing that the theory in this work can be used to predict some fracture characteristics of laminated glass.     

Vibration analysis of a stepped laminated composite Timoshenko beam

The goal of this study is to investigate the vibration characteristics of a stepped laminated composite Timoshenko beam. Based on the first order shear deformation theory, flexural rigidity and transverse shearing rigidity of a laminated beam are determined. In order to account for the effect of shear deformation and rotary inertia of the stepped beam, Timoshenko beam theory is then used to deduce the frequency function. Graphs of the natural frequencies and mode shapes of a T300/970 laminated stepped beam are given, in order to illustrate the influence of step location parameter exerts on the dynamic behavior of the beam.     

Vibration and stability of an axially moving viscoelastic beam with hybrid supports

Vibration and stability are investigated for an axially moving beam constrained by simple supports with torsion springs. A scheme is proposed to derive natural frequencies and modal functions from given boundary conditions of an elastic beam moving at a constant speed. For a beam constituted by the Kelvin model, effects of viscoelasticity on the free vibration are analyzed via the method of multiple scales and demonstrated via numerical simulations. When the axial speed is characterized as a simple harmonic variation about the constant mean speed, the instability conditions are presented for axially accelerating viscoelastic beams in parametric resonance. Numerical examples show the effects of the constraint stiffness, the mean axial speed, and the viscoelasticity.     

Large amplitude free vibration of symmetrically laminated magneto-electro-elastic rectangular plates on Pasternak type foundation

Nonlinear free vibration of symmetrically laminated magneto-electro-elastic rectangular plate resting on an elastic foundation is studied analytically. The plate is considered to be simply supported on all edges. It is also assumed that the magneto-electro-elastic body is poled along the z direction and subjected to electric and magnetic potentials between the upper and lower surfaces. To model the motion of the plate, the first order shear deformation theory along with the Gauss's equations for electrostatics and magnetostatics are used. Then equations of motion are reduced to a single nonlinear ordinary differential equation which is solved analytically by multiple scales method. The results are compared with the published results and good agreement is found. Some numerical examples are presented to investigate the effects of several parameters on the linear and nonlinear behavior of these plates.     

Three-dimensional exact solutions for the free vibration of laminated transversely isotropic circular,annular and sectorial plates with unusual boundary conditions

New state space formulations for the free vibration of circular, annular and sectorial plates are established by introducing two displacement functions and two stress functions. The state variables can be separated into two independent catalogues and two kinds of vibrations can be readily found. Expanding the displacements and stresses in terms of Bessel functions in the radial direction and trigonometric functions in the circumferential direction, we obtained the exact frequency equation for the free vibration for some uncommon boundary conditions. Numerical results are presented and compared with those of FEM to demonstrate the reliability of the proposed method. A parametric investigation is also performed.     

Fracture mechanics of laminated glass subjected to blast loading

A failure criterion based on energy balance approach is introduced for the laminated glass panel subjected to blast loading. Based on this failure criterion, a damage factor is developed to assess the failure of the laminated glass panel. If the damage factor is less than one, the plate is safe otherwise unsafe. Trigonometric function is employed to express the transverse deflection and the Airy’s stress function in von Karman’s large deflection equations of a thin plate. The nonlinear ordinary differential equation of motion obtained using the Galerkin method is solved using Runge–Kutta method. The predicted results indicate that the breakages of the laminated glass may be caused by the negative phase of the blast load if the positive phase blast load is not violent enough to cause failure. Also, the size of glass shards the laminated glass plies breaks in to is predicted using the surface energy based failure model.     

Experimental and numerical study of the rotation and the erosion of fillers suspended in viscoelastic fluids under simple shear flow

When a porous agglomerate immersed in a fluid is submitted to a shear flow, hydrodynamic stresses acting on its surface may cause a size reduction if they exceed the cohesive stress of the agglomerate. The aggregates forming the agglomerate are slowly removed from the agglomerate surface. Such a behaviour is known when the suspending fluid is Newtonian but unknown if the fluid is viscoelastic. By using rheo-optical tools, model fluids, carbon black agglomerates and particles of various shapes, we found that the particles had a rotational motion around the vorticity axis with a period which is independent on shape (flat particles not considered), but which is exponentially increasing with the elasticity of the medium expressed by the Weissenberg number (We). Spherical particles are always rotating for We up to 2.6 (largest investigated We in this study) but elongated particles stop rotating for We>0.9 while orienting along the flow direction. Erosion is strongly reduced by elasticity. Since finite element numerical simulation shows that elasticity increases the local stress around a particle, the origin of the erosion reduction is interpreted as an increase of cohesiveness of the porous agglomerate due to the infiltration of a viscoelastic fluid.     

Mode III fracture behavior of laminated composite with edge crack in torsion

A three-dimensional (3-D) finite element analysis was performed on a [90,(+45/−45)_n,(−45/+45)_n,90]_s class of laminated composites under the edge crack torsion (ECT) test configuration. Finite element delamination models were established and formulas for calculating the Mode III fracture toughness from 3-D finite element models were developed. The relations between the interlaminar fracture behavior and various configuration parameters were investigated and the effects of point loads, ends, geometry, Mode II interference, and friction were evaluated. Results showed that with proper selection of ECT specimen configuration and layup, the delamination could grow in pure Mode III in the middle region of the specimen. Specimen end effect played an important role in the ECT test. A Mode II component occurred in the end regions but it did not interfere significantly with the Mode III delamination state. Specimen dimension ratio, layup, and crack length exhibited significant effect on the interlaminar fracture behavior and the calculated strain energy release rates. However, friction between crackfaces was found to have negligible effect on the interlaminar properties.     

Vibration of a variable cross-section beam

Vibration of an isotropic beam which has a variable cross-section is investigated. Governing equation is reduced to an ordinary differential equation in spatial coordinate for a family of cross-section geometries with exponentially varying width. Analytical solutions of the vibration of the beam are obtained for three different types of boundary conditions associated with simply supported, clamped and free ends. Natural frequencies and mode shapes are determined for each set of boundary conditions. Results show that the non-uniformity in the cross-section influences the natural frequencies and the mode shapes. Amplitude of vibrations is increased for widening beams while it is decreased for narrowing beams.     

Large amplitude free flexural vibration of rings using finite element approach

Here, the large amplitude free flexural vibrations of isotropic/laminated orthotropic rings are investigated, using a shear flexible curved beam element based on field consistency principle. A laminated refined beam theory is introduced for developing the element, which satisfies the interface transverse shear stress and displacement continuity, and has a vanishing shear stress on the inner and outer surfaces of the beam. The formulation includes in-plane and rotary inertia effects, and the non-linearity due to the finite deformation of the ring. The governing equations obtained using Lagrange's equations of motion are solved through the direct integration technique. Amplitude-frequency relationships evaluated from the dynamic response history are examined. Detailed numerical results are presented considering various parameters such as radius-to-thickness ratio, circumferential wave number and ovality for isotropic and laminated orthotropic rings. The nature and degree of the participation of various modes in non-linear asymmetric vibration of oval ring brought out through the present study are useful for accurate modelling of the closed non-circular structures.     

Geometrically nonlinear,steady state vibration of viscoelastic beams

The problem of geometrically non-linear steady state vibrations of beams excited by harmonic forces is considered in this paper. The beams are made of a viscoelastic material defined by the classic Zener rheological model - the simplest model that takes into account all the basic properties of real viscoelastic materials. The constitutive stress-strain relationship for this type of material is given as a differential equation containing derivatives of both stress and strain. This significantly complicates the solution to the problem. The von Karman theory is applied to describe the effects of geometric nonlinearities of beam deformations. The equations of motions are derived using the finite element methodology. A polynomial approximation of bending moments is used. The order of basis functions is set so as to obtain a coherent approximation of moments and displacements. In the steady-state solution of equations of motion, only one harmonic is taken into account. The matrix equations of amplitudes are derived using the harmonic balance method and the continuation method is applied for solving them. The tangent matrix of equations of amplitudes is determined in an explicit form. The stability of steady-state solution is also examined. The resonance curves for beams supported in a different way are shown and the results of calculation are briefly discussed.     

The viscoelastic behavior of melts of virgin and recycled polycarbonate reinforced with short glass fibers

Three series of shear oscillatory tests are performed on polycarbonate melts reinforced with short glass fibers at the temperatures T₁=250 and T₂=290 °C. The content of glass fibers ranges from 0 to 20 wt.%. In the first series, virgin polycarbonate is used, in the other series, dynamic tests are performed on recycled polymer, whereas in the third series, a mixture of virgin with recycled polycarbonates is employed. Constitutive equations are derived for the viscoelastic behavior of a polymer melt at isothermal deformations with small strains. A polymer is treated as an equivalent transient network of strands that rearrange at random times as they are agitated by thermal fluctuations. The time-dependent response of a network is determined by four adjustable parameters that are found by fitting the experimental data. Excellent agreement is demonstrated between the observations and the results of numerical simulation. The study focuses on the effects of temperature and filler content on the material constants in the stress–strain relations.     

Vibration and buckling of laminated sandwich plates having interfacial imperfections

The vibration and buckling characteristics of sandwich plates having laminated stiff layers are studied for different degrees of imperfections at the layer interfaces using a refined plate theory. With this plate theory, the through thickness variation of transverse shear stresses is represented by piece-wise parabolic functions where the continuity of these stresses is satisfied at the layer interfaces by taking jumps in the transverse shear strains at the interfaces. The transverse shear stresses free condition at the plate top and bottom surfaces is also satisfied. The inter-laminar imperfections are represented by in-plane displacement jumps at the layer interfaces and characterized by a linear spring layer model. It is quite interesting to note that this plate model having all these refined features requires unknowns only at the reference plane. To have generality in the analysis, finite element technique is adopted and it is carried out with a new triangular element developed for this purpose, as any existing element cannot model this plate model. As there is no published result on imperfect sandwich plates, the problems of perfect sandwich plates and imperfect ordinary laminates are used for validation.     

Rate effect on sharing of passive lumbar motion segment under load-controlled sagittal flexion: viscoelastic finite element analysis

Industrial epidemiological studies have shown that jobs requiring a higher speed of trunk motion contribute to a higher risk of industrial low back disorders. Consideration of the loading dynamic characteristics, such as lifting at different speeds, requires modeling of the viscoelastic behavior of passive tissues. Detailed systematic analysis of the effects of loading rate has been lacking in the literature. A validated viscoelastic finite element model of a L2–L3 motion segment was used to identify the load sharing among the passive elements at different loading rates. Force controlled complex flexion movement was simulated by applying load at the top of the upper vertebra without constraining any coupled sagittal rotation, whereas the lower vertebra was fixed at the bottom. The load reached its maximum values of 2000 N compression, 400 N anterior shear, and 20 Nm flexion in three different durations of 0.3, 3 and 30 s to represent fast, medium and slow movement. The global force–displacement response of the motion segment, forces in facet joints and ligaments, stresses and strains in anulus fibrosus, and intradiscal pressure were compared across different rates. The higher rate of loading while reaching a prescribed complex forward flexion loading increased the intradiscal pressure and the stress in the anulus fibers at the posterolateral innermost layers, but reduced the global displacements, ligament forces and facet joint forces. The distribution of stress and strain was markedly affected by the loading rate. Consideration of the time-dependent material properties of passive elements is essential to improve our understanding of the responses of the motion segment to dynamic loading conditions. Speed of the manual materials handling (MMH) tasks should be included as a risk factor in the biomechanical and epidemiological studies and guidelines for safe lifting.     

基于边光滑有限元法的加筋板静力和自由振动分析

运用边光滑有限元法,研究分析了加筋板结构的静力和自由振动问题。在边光滑有限元法中,将基于边的应变光滑技术用于对原来的应变场进行光滑操作;由于应变光滑技术能够适当地软化原来过刚的有限元模型,从而能够得到更加接近于系统准确刚度的光滑有限元模型;鉴于三角形单元良好的适用性,选用三角形单元对模型进行网格划分;同时,为了解决低阶Reissner-Mindlin板单元弯曲过程中的横向剪切自锁问题,采用了一种新型的离散剪切间隙技术。算例的数值计算结果表明,与传统的有限元法相比,边光滑有限元法能够得到精度更高的计算结果,且收敛更快,计算效率更佳。     

Size-dependent free vibration analysis of composite laminated Timoshenko beam based on new modified couple stress theory

A micro-scale free vibration analysis of composite laminated Timoshenko beam (CLTB) model is developed based on the new modified couple stress theory. In this theory, a new anisotropic constitutive relation is defined for modeling the CLTB. This theory uses rotation–displacement as dependent variable and contains only one material length scale parameter. Hamilton’s principle is employed to derive the governing equations of motion and boundary conditions. This new model can be reduced to composite laminated Bernoulli–Euler beam model of the couple stress theory. An example analysis of free vibration of the cross-ply simply supported CLTB model is adopted, and an explicit expression of analysis solution is given. Additionally, the numerical results show that the present beam models can capture the scale effects of the natural frequencies of the micro-structure.     

重新认识横力弯曲梁的切应力强度效应

通过多层叠合梁的横力弯曲问题分析, 证明弯曲切应力不仅直接导致多层叠合梁的脱层, 而且这种脱层还会进一步引发弯曲正应力的急剧增加. 这种由弯曲切应力引发的连锁效应和潜在危害, 理应引起同学们的足够重视.     

Nonlinear free vibration analysis of piezoelastic laminated plates with interface damage

This paper presents a nonlinear model for piezoelastic laminated plates with damage effect of the intra-layers and inter-laminar interfaces. Discontinuity of displacement and electric potential on the interfaces are depicted by three shape functions. By using the Hamilton variation principle, the three-dimensional nonlinear dynamic equations of piezoelastic laminated plates with damage effect are derived. Then, by using the Galerkin method, a mathematical solution is presented. In the numerical studies, effects of various factors on the natural frequencies and nonlinear amplitude-frequency response of the simply-supported peizoelastic laminated plates with interfacial imperfections are discussed. These factors include different damage models, thickness of the piezoelectric layer, side-to-thickness ratio, and length-to-width ratio.     

Parameter determination and response analysis of viscoelastic material

The parameter determination of viscoelastic material is a multi-variable, multi-aim nonlinear optimization problem, which made the optimization process very complicated. In this paper a hybrid optimal algorithm was proposed to determine the viscoelastic parameters in the constitutive relation according to the experimentally obtained mechanical properties. This algorithm merges the Broydon–Fletcher–Goldfarb–Shanno search into a genetic algorithm framework as a basic operator in order to enhance the local search capability. The proposed hybrid algorithm not only can reduce the iterative times greatly but can abolish the limitation of initial parameter values. Nonlinear material characteristic curve-fitting was carried out using the proposed algorithm and other existing approaches. And the comparison results show this algorithm is accurate and effective. The numerical simulation and experimental study of viscoelastic cantilever beam also indicates that the finite element formulation and the calculative viscoelastic model parameters are reliable. The proposed optimization method can be extended to further complex parameter estimation researches.