Exact dynamic solutions to piezoelectric smart beams including peel stresses Part II: Numerical results,comparison and approximate solution

This work presents numerical results for the exact dynamic solution of piezoelectric (PZT) smart beams including peel stresses, which was developed in Part I. Numerical results are presented in details for frequency spectra, natural frequencies, normal mode shapes, harmonic responses of the shear and peel stresses, and sensing electric charges for a cantilever beam with a bonded PZT patch to the clamped end. The exact dynamic solution can provide useful data for benchmarking other methods. The numerical results of the present model including peel stresses (PSM) are also compared with those obtained using the shear lag beam model and the shear lag rod model. On the basis of the equivalent forces derived in the static analysis, simple approximate dynamic solutions are obtained and compared with the exact solutions, and then the application and limitation of the simple approximate solutions are investigated. By comparing numerical results predicted by the present PSM model with the shear lag models and the approximate solutions based on the static equivalent forces, effects of the dynamic shear and peel stresses on natural frequencies and dynamic responses of the smart structures are examined.     

功能梯度梁与均匀梁静动态解间的相似转换

基于Euler-Bernoulli 梁理论, 研究了功能梯度材料梁的弯曲、屈曲和自由振动. 通过分析和比较功能梯度材料梁 和均匀梁的控制方程, 得到了功能梯度材料梁与均匀梁的解之间的相似转换关系. 在给定功 能梯度材料梁的材料性质在横向按幂函数分布的情况下, 导出了解之间的相似转换系数的解 析表达式. 该系数集中反映了功能梯度梁的材料非均匀性. 因此, 可将功能梯度材料梁的静 动态问题的求解转换为同样载荷和边界条件下均匀梁的静动态问题求解以及相似转换系数的 计算.     

Explicit solutions for shearing and radial stresses in curved beams

In this paper the formulae for the shearing and radial stresses in curved beams are derived analytically based on the solution for a Volterra integral equation of the second kind. These formulae satisfy both the equilibrium equations and the static boundary conditions on the surfaces of the beams. As some applications, the resulting solutions are used to calculate the shearing and radial stresses in a cantilevered curved beam subjected to a concentrated force at its free end. The numerical results are compared with other existing approximate solutions as well as the corresponding solutions based on the theory of elasticity. The calculations show a better agreement between the present solution and the one based on the theory of elasticity. The resulting formulae can be applied to more general cases of curved beams with arbitrary shapes of cross-sections.     

Two-dimensional elasticity solution for bending of functionally graded beams with variable thickness

This paper studies the stress and displacement distributions of functionally graded beam with continuously varying thickness, which is simply supported at two ends. The Young’s modulus is graded through the thickness following the exponential-law and the Poisson’s ratio keeps constant. On the basis of two-dimensional elasticity theory, the general expressions for the displacements and stresses of the beam under static loads, which exactly satisfy the governing differential equations and the simply supported boundary conditions at two ends, are analytically derived out. The unknown coefficients in the solutions are approximately determined by using the Fourier sinusoidal series expansions to the boundary conditions on the upper and lower surfaces of the beams. The effect of Young’s modulus varying rules on the displacements and stresses of functionally graded beams is investigated in detail. The two-dimensional elasticity solution obtained can be used to assess the validity of various approximate solutions and numerical methods for the aforementioned functionally graded beams.     

复合材料叠层梁和金属梁的固有振动特性

对根据三种梁理论得到的金属梁和复合材料叠层梁的固有振动特性进行了对比性的研究对常用的三种梁理论在弹性碰撞分析中的应用进行了分析和比较     

Solution of generalized coordinate for warping for naturally curved and twisted beams

A theoretical method for static analysis of naturally curved and twisted beams under complicated loads was presented, with special attention devoted to the solving process of governing equations which take into account the effects of torsion-related warping as well as transverse shear deformations. These governing equations, in special cases, can be readily solved and yield the solutions to the problem. The solutions can be used for the analysis of the beams, including the calculation of various internal forces, stresses, strains and displacements. The present theory will be used to investigate the stresses and displacements of a plane curved beam subjected to the action of horizontal and vertical distributed loads. The numerical results obtained by the present theory are found to be in very good agreement with the results of the FEM results. Besides, the present theory is not limited to the beams with a double symmetric cross-section, it can also be extended to those with arbitrary cross-sectional shape.     

基于一种广义梁理论的弹性地基FGM简支梁自由振动解析解

基于一种扩展的n阶广义剪切变形梁理论(n-GBT),应用Hamilton原理,建立了以轴向位移、横向位移及转角为未知函数的Winkler-Pasternak弹性地基功能梯度材料(FGM)梁的自由振动方程,采用Navier法获得了弹性地基FGM简支梁自由振动的精确解。与多种梁理论预测结果进行比较,讨论并给出了GBT阶次n的理想取值;分析了梯度指标、跨厚比及地基刚度对FGM梁频率的影响。结果表明:本文方法有效且适用范围广,若采用高阶剪切梁理论模型,宜取n≥3的奇数;FGM梁的自振频率随材料梯度指标的增大而减小;随跨厚比的增加而增大,但当跨厚比大于20,跨厚比增加对频率的影响很小;随地基刚度的增加而增大,地基刚度足够大时,频率趋于收敛。     

Photoelastic analysis with an interferometer simulating reinforced simple concrete beams

This investigation was undertaken to study the static stresses developed in reinforced models simulating reinforced-concrete beams. The models maintained the correct modular ratio of the reinforcing rod to the beam material. In addition, the percentage of reinforcing material was made to produce an elastically balanced beam. The stresses were obtained with a Favre interferometer. Comparison of the boundary stresses was made with a standard circular polariscope using monochromatic light. Samples of the stresses plotted are provided to show graphically the distribution of the stresses obtained in the analyses.     

Cohesive zone modeling of interfacial stresses in plated beams

The elastic analysis of interfacial stresses in plated beams has been the subject of several investigations. These studies provided both first-order and higher-order solutions for the distributions of interfacial shear and normal stresses close to the plate end in the elastic range. The notable attention devoted to this topic was driven by the need to develop predictive models for plate end debonding mechanisms, as the early models of this type adopted debonding criteria based on interfacial stresses. Currently, approaches based on fracture mechanics are becoming increasingly established. Cohesive zone modeling bridges the gap between the stress- and energy-based approaches. While several cohesive zone analyses of bonded joints subjected to mode-II loading are available, limited studies have been conducted on cohesive zone modeling of interfacial stresses in plated beams. Moreover, the few available studies present complex formulations for which no closed-form solutions can be found. This paper presents an analytical cohesive zone model for the determination of interfacial stresses in plated beams. A first-order analysis is conducted, leading to closed-form solutions for the interfacial shear stresses. The mode-II cohesive law is taken as bilinear, as this simple shape is able to capture the essential properties of the interface. A closed-form expression for the debonding load is proposed, and the comparison between cohesive zone modeling and linear-elastic fracture mechanics predictions is discussed. Analytical predictions are also compared with results of a numerical finite element model where the interface is described with zero-thickness contact elements, using the node-to-segment strategy and incorporating decohesion and contact within a unified framework.     

Higher-order impact modeling of sandwich structures with flexible core

In this study, a higher-order impact model is presented to simulate the response of a soft-core sandwich beam subjected to a foreign object impact. A free vibration problem of sandwich beams is first solved, and the results are validated by comparing with numerical finite element modeling results of ABAQUS and the solution by Frostig and Baruch [Frostig, Y., Baruch, M., 1994. Free vibration of sandwich beams with a transversely flexible core: a high order approach. Journal of Sound and Vibration 176(2), 195–208]. Then a foreign object impact process is incorporated in the higher-order model, and the contact force and deflection history as well as the propagation of transverse normal, shear, and axial stresses during the impact are analyzed and discussed. The validity of the model in the impact response predictions is demonstrated by comparing with finite element solutions of LS-DYNA. The calculated stresses caused by a foreign object impact are then used to assess failure locations, failure time, and failure modes in sandwich beams, which are shown to compare well with the available experimental results. The effects of impact mass, initial velocity, core stiffness, and core height on the impact stresses generated in the beams are discussed. The influences of impact mass and initial velocity on the contact force history are close to those by the linearized impact solution, but the proposed higher-order impact model captures the non-linear impact process and different generated stresses. Compared to the fully backed sandwich case, the core height shows a great influence over the impact process of a simply supported sandwich system, in which the global behavior of the sandwich is dominant; while the core stiffness shows minor effect over the impact process. The higher-order impact model of sandwich beams developed in the study provides accurate predictions of the generated stresses and impact process and can be used effectively in design analysis of anti-impact structures made of sandwich materials.     

冲击荷载作用下简支圆板的塑性动力响应统一解

采用统一强度理论求解了简支圆板在中等脉冲荷载作用下的动力响应问题,得出了统一的动力塑性极限荷载、内力场和速度场,并给出了上限解和下限解。讨论了静力许可条件和运动许可条件。利用本文的解还得出了简支圆板在静力荷载作用下的极限荷载、内力场和速度场。根据选择不同的拉压比参数,本文所给出的解可以适用于各种拉压异性和拉压同性材料。Tresca解、Mohr Coulomb解和双剪统一屈服准则解是本文的特例,Mises解是本文当=1和b=0.5时的线性逼近。研究结果表明,拉压比和强度理论参数b对动力解的影响要大于对静力解的影响,所以,根据材料的不同选择合适的强度理论,对于更好的发挥材料的强度潜力,减轻结构的重量具有重要的意义。     

Planar multiple-contact problems subject to unilateral and bilateral kinetic constraints with static Coulomb friction

This paper is concerned with contact problems. A planar multiple-contact problem subject to unilateral and bilateral kinetic constraints with static friction is studied using the complementarity method. First, this paper discusses the one-to-one correspondence of solutions of the contact problems of concern and of the corresponding complementarity models. An enhanced complementarity model is proposed by adding missed tangential acceleration constraints into previous complementarity models. Solutions of the proposed complementarity model and solutions of the contact problem are proven to exhibit one-to-one correspondence, which may not be guaranteed in the previous complementarity models. Then, this paper applies linear complementarity theory to investigate the properties of the solutions of the proposed complementarity model. For both unilaterally constrained contact problems and bilaterally constrained contact problems, the existence of solutions and boundedness of solutions are proven. Sufficient conditions for the uniqueness of solutions and finiteness of the number of solutions are also provided. Several numerical examples are given to show the non-uniqueness of solutions or the infiniteness of the number of solutions. Such phenomena demonstrate the non-smoothness of the contact problems discussed herein.     

A higher-order theory for static and dynamic analyses of functionally graded beams

The higher-order theory is extended to functionally graded beams (FGBs) with continuously varying material properties. For FGBs with shear deformation taken into account, a single governing equation for an auxiliary function F is derived from the basic equations of elasticity. It can be used to deal with forced and free vibrations as well as static behaviors of FGBs. A general solution is constructed, and all physical quantities including transverse deflection, longitudinal warping, bending moment, shear force, and internal stresses can be represented in terms of the derivatives of F. The static solution can be determined for different end conditions. Explicit expressions for cantilever, simply supported, and clamped-clamped FGBs for typical loading cases are given. A comparison of the present static solution with existing elasticity solutions indicates that the method is simple and efficient. Moreover, the gradient variation of Young’s modulus and Poisson’s ratio may be arbitrary functions of the thickness direction. Functionally graded Rayleigh and Euler–Bernoulli beams are two special cases when the shear modulus is sufficiently high. Moreover, the classical Levinson beam theory is recovered from the present theory when the material constants are unchanged. Numerical computations are performed for a functionally graded cantilever beam with a gradient index obeying power law and the results are displayed graphically to show the effects of the gradient index on the deflection and stress distribution, indicating that both stresses and deflection are sensitive to the gradient variation of material properties.     

POLYNOMIAL SOLUTIONS TO PIEZOELECTRIC BEAMS (Ⅱ)——ANALYTICAL SOLUTIONS TO TYPICAL PROBLEMS

IntroductionThis paper is a continuation of Ref.[1],in which a series of orthotropic piezoelectricplane problems was solved and the corresponding exact solutions were obtained with the trial-and-error method,on the basis of the general solution expressed …     

POLYNOMIAL SOLUTIONS TO PIEZOELECTRIC BEAMS(Ⅱ)--ANALYTICAL SOLUTIONS TO TYPICAL PROBLEMS

For the orthotropic piezoelectric plane problem, a series of piezoelectric beams is solved and the corresponding analytical solutions are obtained with the trialand-error method on the basis of the general solution in the case of three distinct eigenvalues, in which all displacements, electrical potential, stresses and electrical displacements are expressed by three displacement functions in terms of harmonic polynomials. These problems are cantilever beam with cross force and point charge at free end, cantilever beam and simply-supported beam subjected to uniform loads on the upper and lower surfaces, and cantilever beam subjected to linear electrical potential.     

Boundary element method for model analysis of 2-D composite structure

The boundary element method is used for the modal analysis of free vibration of 2-D composite structures in this paper. Since the particular solution method is used to treat the terms of body forces (inertial forces) in the equation of motion, only static fundamental solutions are needed in solving the problem. For an isotropic cantilever beam, the numerical results obtained by using the BEM presented in this paper are in good agreement, with, those of using FEM or other BEM, but this BEM can also be used to analyze problems for anisotropic materials. For simply supported composite laminated beams, the comparisons of the numerical reslts obtained by this method with the analytical results obtained by 1-D laminated beam theory indicate that if the ratio of length/thickness is greater than 20, the results of the two methods are in good agreement, but if the ratio of length/thickness is less than 20, big errors will occur for 1-D laminated beam theory.     

Dynamic simulation of frictional contacts of thin beams during large overall motions via absolute nodal coordinate formulation

The aim of this study is to develop an approach of simulating the frictional contact dynamics of thin beams with large deformations and continuous contact zones of large size during their large overall motions. For this purpose, the thin beams are meshed via initially straight and gradient deficient thin beam elements of the absolute nodal coordinate formulation (ANCF) degenerated from a curved beam element of ANCF. A detection strategy for contact zone is proposed based on the combination of the minimal distance criterion and master-slave approach. By making use of the minimal distance criterion, the closest points of two thin beams can be found efficiently. The master-slave approach is employed to determine the continuous contact zone. The generalized frictional contact forces and their Jacobians are derived based on the principle of virtual work. Gauss integration is used to integrate the contact forces over the continuous contact zone. The generalized-alpha method is used to solve the dynamic equations of contacting beams. Numerical simulations of four static and dynamic contact problems, including those with continuous contact zones of large size, are completed to validate the high performance of the approach.     

Homogenized and classical expressions for static bending solutions for functionally graded material Levinson beams

The exact relationship between the bending solutions of functionally graded material (FGM) beams based on the Levinson beam theory and those of the corresponding homogenous beams based on the classical beam theory is presented for the material properties of the FGM beams changing continuously in the thickness direction. The deflection, the rotational angle, the bending moment, and the shear force of FGM Levinson beams (FGMLBs) are given analytically in terms of the deflection of the reference homogenous Euler-Bernoulli beams (HEBBs) with the same loading, geometry, and end supports. Consequently, the solution of the bending of non-homogenous Levinson beams can be simplified to the calculation of transition coefficients, which can be easily determined by variation of the gradient of material properties and the geometry of beams. This is because the classical beam theory solutions of homogenous beams can be easily determined or are available in the textbook of material strength under a variety of boundary conditions. As examples, for different end constraints, particular solutions are given for the FGMLBs under specified loadings to illustrate validity of this approach. These analytical solutions can be used as benchmarks to check numerical results in the investigation of static bending of FGM beams based on higher-order shear deformation theories.     

基于初应力法的作大范围运动柔性梁动力学建模理论研究

研究了初应力法的作大范围运动柔性梁的建模理论.根据连续介质理论,考虑应变-位移中的非线性项,用一致质量有限元法对柔性梁进行离散,基于Jourdain速度变分原理导出定轴转动下大范围运动为自由的柔性梁刚-柔耦合动力学方程.从其刚柔耦合动力学方程出发,考虑在大范围运动已知情况下的结构动力学方程.通过引入准静态概念,把其结构动力学方程转化为准静态方程.对纵向和横向变形节点坐标进行坐标分离,解出与纵向变形相关的准静态方程,得到准静态时的纵向应力表达式,从而获得附加刚度项.并对此非惯性系下作大范围运动柔性梁的结构动力学方程进行数值仿真,对零次近似模型、一次近似模型、初应力法动力学模型的仿真结果进行分析,揭示三种模型的动力学性质的差异.     

In-plane free vibration and stability of loaded and shear-deformable circular arches

The first known equations governing vibrations of preloaded, shear-deformable circular arches are derived according to a variational principle for dynamic problems concerning an elastic body under equilibrium initial stresses. The equations are three partial differential equations with variable coefficients. The governing equations are solved for arches statically preloaded with a uniformly distributed vertical loading, by obtaining a static, closed-form solution and an analytical dynamic solution from series solutions and dynamic stiffness matrices. Convergence to accurate results is obtained by increasing the number of elements or by increasing both the number of terms in the series solution and the number of terms in the Taylor expansion of the variable coefficients. Graphs of non-dimensional frequencies and buckling loads are presented for preloaded clamped arches. They clarify the effects of opening angle and thickness-to-radius ratio on vibration frequencies and buckling loads. The effects of static deformations on vibration frequencies are also investigated. This work also compares the results obtained from the proposed governing equations with those obtained from the classical theory neglecting shear deformation.