轴向力作用任意不连续梁自由振动分析的广义函数解

摘要：首先运用分布理论建立了轴向力作用下含多个不连续点的欧拉梁的自由振动的统一微分方程。不连续点的影响由广义函数（Dirac delta函数）引入梁的振动方程。微分方程运用Laplace变换方法求解;与传统方法不同的是,本文方法适用于含任意类型的不连续点和多种不连续点组合情况的梁,求得的模态函数为整个不连续梁的一般解。由于模态函数的统一化以及连续条件的退化,特征值的求解得到了极大的简化。最后,以轴向力作用下多跨梁—弹簧质量块系统模型为例验证了本文方法的正确性与有效性。     

On Euler–Bernoulli discontinuous beam solutions via uniform-beam Green’s functions

The bending problem of Euler–Bernoulli discontinuous beams is dealt with. The purpose is to show that uniform-beam Green’s functions can be used to build efficient solutions for beams with internal discontinuities due to along-axis constraints and flexural-stiffness jumps. Specifically, upon deriving the equilibrium equation in the space of generalized functions, first it is seen that the original bending problem may be recast as linear superposition of a principal and an auxiliary bending problem, both involving a uniform reference beam and homogeneous boundary conditions. Then, based on the Green’s functions of the reference beam, closed-form solutions are developed for the principal beam response, while the auxiliary beam response is obtained by solving, in general, (r + 2s) algebraic equations written at the discontinuity locations, being r the number of discontinuities due to along-axis constraints, and s the number of flexural-stiffness jumps. In this manner, an appreciable reduction of computational effort is achieved as compared to alternative analytical solutions in the literature.     

Euler–Bernoulli beams with multiple singularities in the flexural stiffness

Euler–Bernoulli beams under static loads in presence of discontinuities in the curvature and in the slope functions are the object of this study. Both types of discontinuities are modelled as singularities, superimposed to a uniform flexural stiffness, by making use of distributions such as unit step and Dirac's delta functions. A non-trivial generalisation to multiple different singularities of an integration procedure recently proposed by the authors for a single singularity is presented in this paper. The proposed integration procedure leads to closed form solutions, dependent on boundary conditions only, which do not require enforcement of continuity conditions along the beam span. It is however shown how, from the solution of the clamped-clamped beam, by considering suitable singularities at boundaries in the flexural stiffness model, responses concerning several boundary conditions can be recovered. Furthermore, solutions in terms of deflection of the beam are obtained for imposed displacements at boundaries providing the so called shape functions. The above mentioned shape functions can be adopted to insert beams with singularities as frame elements in a finite element discretisation of a frame structure. Explicit expressions of the element stiffness matrix are provided for beam elements with multiple singularities and the reduction of degrees of freedom with respect to classical finite element procedures is shown. Extension of the proposed procedure to beams with axial displacement and vertical deflection discontinuities is also presented.     

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.     

Interaction between elliptical and ellipsoidal inclusions under bending stress fields

Summary  This paper deals with interaction problems of elliptical and ellipsoidal inclusions under bending, using singular integral equations of the body force method. The problems are formulated as a system of singular integral equations with Cauchy-type or logarithmic-type singularities, where unknown functions are densities of body forces distributed in the x,y and r,θ,z directions in infinite bodies having the same elastic constants as those of the matrix and inclusions. In order to satisfy the boundary conditions along the elliptical and the ellipsoidal boundaries, the unknown functions are approximated by a linear combination of fundamental density functions and polynomials. The present method is found to yield the exact solutions for a single elliptical or spherical inclusion under a bending stress field. It yields rapidly converging numerical results for interface stresses in the interaction of inclusions. Received 9 September 1999; accepted for publication 15 January 2000     

Stress-strain state of non-thin plates and shells. Generalized theory (survey)

Conclusions Thus, this part of our survey has presented the main approaches that have been taken to the construction of two-dimensional (in terms of the space coordinates) equations of a generalized theory of plates and shells. The solutions of these equations represent a certain approximation of the solution of the initial three-dimensional problem. They are based on expansion of the sought functions into Fourier series in Legendre polynomials of the thickness coordinate. Studies completed on the basis of the given variants of plate and shell theory were systematized and analyzed. In terms of the method of its construction, the theory involves a regular process of replacing the solution of the three-dimensional problem by the solution (or sequence of solutions) of two-dimensional boundary-value problems or initial-boundary-value problems. Numerical results illustrating the convergence of the successive approximation were presented. It should be noted that to make comparison with the results of classical or applied theories, several of the studies cited here presented solutions of problems for thin plates and shells with allowance only for the initial terms of expansions of the stress and displacement components into base functions (Legendre polynomials).S. P. Timoshenko Institute of Mechanics, Ukrainian Academy of Sciences, Kiev. Translated from Prikladnaya Mekhanika, Vol. 29, No. 11, pp. 3–34, November, 1993.     

Problems of the Flamant–Boussinesq and Kelvin Type in Dipolar Gradient Elasticity

This work studies the response of bodies governed by dipolar gradient elasticity to concentrated loads. Two-dimensional configurations in the form of either a half-space (Flamant–Boussinesq type problem) or a full-space (Kelvin type problem) are treated and the concentrated loads are taken as line forces. Our main concern is to determine possible deviations from the predictions of plane-strain/plane-stress classical linear elastostatics when a more refined theory is employed to attack the problems. Of special importance is the behavior of the new solutions near to the point of application of the loads where pathological singularities and discontinuities exist in the classical solutions. The use of the theory of gradient elasticity is intended here to model material microstructure and incorporate size effects into stress analysis in a manner that the classical theory cannot afford. A simple but yet rigorous version of the generalized elasticity theories of Toupin (Arch. Ration. Mech. Anal. 11:385–414, 1962) and Mindlin (Arch. Ration. Mech. Anal. 16:51–78, 1964) is employed that involves an isotropic linear response and only one material constant (the so-called gradient coefficient) additional to the standard Lamé constants (Georgiadis et al., J. Elast. 74:17–45, 2004). This theory, which can be viewed as a first-step extension of the classical elasticity theory, assumes a strain-energy density function, which besides its dependence upon the standard strain terms, depends also on strain gradients. The solution method is based on integral transforms and is exact. The present results show departure from the ones of the classical elasticity solutions (Flamant–Boussinesq and Kelvin plane-strain solutions). Indeed, continuous and bounded displacements are predicted at the points of application of the loads. Such a behavior of the displacement fields is, of course, more natural than the singular behavior present in the classical solutions.      

Exact solutions of solid-liquid phase-change heat transfer when subjected to convective boundary conditions

A closed-form expression is presented to find the location of solid-liquid interface motion in convection-dominated solidification and melting problems. In this regard, the solutions are expressed in terms of the generalized representations of error functions,E (u, v) andF (u, v), which are useful to heat-conduction problems with convective-type boundary conditions. It is demonstrated that for constant surface temperature, the interface solution reduces to the classical Neumann solution.     

POLYNOMIAL SOLUTIONS TO PIEZOELECTRIC BEAMS(Ⅰ)--SEVERAL EXACT SOLUTIONS

For the orthotropic piezoelectric plane problem, a series of piezoelectric beams is solved and the corresponding exact solutions are obtained with the trial-anderror 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 rectangular beams having rigid body displacements and identical electrical potential, rectangular beams under uniform tension and electric displacement as well as pure shearing and pure bending, beams of two free ends subjected to uniform electrical potential on the upper and lower surfaces.     

POLYNOMIAL SOLUTIONS TO PIEZOELECTRIC BEAMS (Ⅰ)——SEVERAL EXACT SOLUTIONS

IntroductionDue to its excellent piezoelectric properties,composites made of piezoelectric materialsare found widespread applications and attracted more attentions[1-10].Because of materialanisotropy and couplingbetween mechanical deformation and electric…     

Discontinuous fold bifurcations in mechanical systems

Summary  This paper treats discontinuous fold bifurcations of periodic solutions of discontinuous systems. It is shown how jumps in the fundamental solution matrix lead to jumps of the Floquet multipliers of periodic solutions. A Floquet multiplier of a discontinuous system can jump through the unit circle, causing a discontinuous bifurcation. Numerical examples are treated, which show discontinuous fold bifurcations. A discontinuous fold bifurcation can connect stable branches to branches with infinitely unstable solutions. Received 20 September 2000; accepted for publication 26 June 2001     

A novel beam finite element with singularities for the dynamic analysis of discontinuous frames

In this work, a model of the stepped Timoshenko beam in presence of deflection and rotation discontinuities along the span is presented. The proposed model relies on the adoption of Heaviside’s and Dirac’s delta distributions to model abrupt and concentrated, both flexural and shear, stiffness discontinuities of the beam that lead to exact closed-form solutions of the elastic response in presence of static loads. Based on the latter solutions, a novel beam element for the analysis of frame structures with an arbitrary distribution of singularities is here proposed. In particular, the presented closed-form solutions are exploited to formulate the displacement shape functions of the beam element and the relevant explicit form of the stiffness matrix. The proposed beam element is adopted for a finite element discretization of discontinuous framed structures. In particular, by means of the introduction of a mass matrix consistent with the adopted shape functions, the presented model allows also the dynamic analysis of framed structures in presence of deflection and rotation discontinuities and abrupt variations of the cross-section. The presented formulation can also be easily employed to conduct a dynamic analysis of damaged frame structures in which the distributed and concentrated damage distributions are modelled by means of equivalent discontinuities. As an example, a simple portal frame, under different damage scenarios, has been analysed and the results in terms of frequency and vibration modes have been compared with exact results to show the accuracy of the presented discontinuous beam element.     

General finite element description for non-uniform and discontinuous beam elements

The theory of generalized functions is used to address the static equilibrium problem of Euler–Bernoulli non-uniform and discontinuous 2-D beams. It is shown that if simple integration rules are applied, the full set of response variables due to end nodal displacements and to in-span loads can be derived, in a closed form, for most common beam profiles and arbitrary discontinuity parameters. On this basis, for finite element analysis purposes, a non-uniform and discontinuous beam element is implemented, for which the exact stiffness matrix and the fixed-end load vector are derived. Upon computing the nodal response, no numerical integration is required to build the response variables along the beam element.     

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.     

Stability of nonlinear periodic vibrations of 3D beams

The geometrically nonlinear periodic vibrations of beams with rectangular cross section under harmonic forces are investigated using a p-version finite element method. The beams vibrate in space; hence they experience longitudinal, torsional, and nonplanar bending deformations. The model is based on Timoshenko’s theory for bending and assumes that, under torsion, the cross section rotates as a rigid body and is free to warp in the longitudinal direction, as in Saint-Venant’s theory. The theory employed is valid for moderate rotations and displacements, and physical phenomena like internal resonances and change of the stability of the solutions can be investigated. Green’s nonlinear strain tensor and Hooke’s law are considered and isotropic and elastic beams are investigated. The equation of motion is derived by the principle of virtual work. The differential equations of motion are converted into a nonlinear algebraic form employing the harmonic balance method, and then solved by the arc-length continuation method. The variation of the amplitude of vibration in space with the excitation frequency of vibration is determined and presented in the form of response curves. The stability of the solution is investigated by Floquet’s theory.     

非傅里叶热弹性的时域间断迦辽金有限元方法

基于Lord-Shulman非傅里叶热弹性模型,提出了采用修正的时域间断迦辽金有限元方法(time discontinuousGalerkin finite element method, DGFEM)求解方法. DGFEM对温度场、位移场基本未知向量及其时间导数向量在时域中分别插值;在最终的求解公式中,引入了人工阻尼. 数值结果显示所发展的DGFEM 较好地捕捉了波的间断并消除了热冲击作用下虚假的数值振荡,能够良好地模拟热弹性问题并具有较高的精度.     

Generalized thermoelastic analysis of the stress-focusing effect in a solid cylinder under rotationally asymmetrical instantaneous thermal loading

Summary  This paper investigates the stress-focusing effect in an infinitely long cylinder under rotationally asymmetrical instantaneous thermal loading on the basis of the generalized thermoelastic Lord–Shulman (L-S) and Green–Lindsay (G-L) theories. Combined forms of the governing equations of both theories are given in a cylindrical coordinate system. The two-dimensional generalized thermoelastic problems are solved by numerical inversion of Laplace transform. Calculations have been performed to find distributions of thermal stresses on the basis of the L-S theory. Stress-focusing phenomena under different heating conditions are presented. The effects of thermomechanical coupling and relaxation time on the stress-focusing phenomena as well as the singularity of stresses are discussed. Received 15 November 2000; accepted for publication 15 November 2001     

Accurate modelling of piezoelectric plates: single-layered plate

Summary  The paper presents an efficient two-dimensional approach to piezoelectric plates in the framework of linear theory of piezoelectricity. The approximation of the through-the-thickness variations accounts for the shear effects and a refinement of the electric potential. Using a variational formalism, electromechanically coupled plate equations are obtained for the generalized stress resultants as well as for the generalized electric inductions. The latter are deduced from the conservative electric charge equation, which plays a crucial role in the present model. Emphasis is placed on the boundary conditions at the plate faces. The model is used to examine some problems for cylindrical bending of a single simply supported plate. Number of situations are examined for a piezoelectric plate subject to (i) an applied electric potential, (ii) a surface density of force, and (iii) a surface density of electric charge. The through-thickness distributions of electromechanical quantities (displacements, stresses, electric potential and displacement) are obtained, and compared with results provided by finite element simulations and by a simplified plate model without shear effects. A good agreement is observed between the results coming from the present plate model and finite element computations, which ascertains the effectiveness of the proposed approach to piezoelectric plates. Received 17 July 2000; accepted for publication 26 September 2000     

A function theory for thin elastic shells

It is well known in the theory of elastic shells that a first order approximation using the shell thickness as an expansion parameter leads to the membrane theory of shells. The membrane equations have as solutions thegeneralized analytic functions. These functions have been exhaustively studied by Ilya N. Vekua [6], [7] and his students. R.P. Gilbert and J. Hile [3] introduced an extension of these systems to include elliptic systems of 2n equations in the plane and named the solutions of these systemsgeneralized hyperanalytic functions.It is shown in this paper that the next order approximation to the shell, which permits, moreover, the introduction of bending, may be described in terms of the generalized hyperanalytic functions. It is strongly suspected that the higher order approximations may also be described in terms of corresponding hypercomplex systems.     

近场动力学方法及其应用

近场动力学(peridynamics,PD)是一种新兴的基于非局部作用思想建立模型并通过求解空间积分方程描述物质力学行为的方法.它兼有分子动力学方法和无网格方法的优点,避免了基于连续性假设建模和求解空间微分方程的传统宏观方法在面临不连续问题时的奇异性,又突破了经典分子动力学方法在计算尺度上的局限,在宏/微观不连续力学问题分析中均表现出很高的求解精度和效率.首先概述了PD方法的理论基础、建模思路和计算体系;进而介绍了PD方法在不同尺度不连续力学问题中的应用,包括均匀与非均匀材料和结构的大变形、损伤、断裂、冲击、穿透和失稳问题,结晶相变动力学问题以及纳米材料和结构的破坏问题;最后讨论了PD方法在理论、计算和应用等方面值得进一步研究的问题.