Accurate and straightforward symplectic approach for fracture analysis of fractional viscoelastic media

An accurate and straightforward symplectic method is presented for the fracture analysis of fractional two-dimensional(2D) viscoelastic media. The fractional Kelvin-Zener constitutive model is used to describe the time-dependent behavior of viscoelastic materials. Within the framework of symplectic elasticity, the governing equations in the Hamiltonian form for the frequency domain(s-domain) can be directly and rigorously calculated. In the s-domain, the analytical solutions of the displacement ...     

辛体系下平面热黏弹性圣维南问题的解析解

借助积分变换,将辛体系引入平面热黏弹性问题,建立了基本问题的对偶方程,并将全 部圣维南问题归结为满足共轭辛正交关系的零本征值本征解问题. 同时,利用变量代换和本 征解展开等技术给出了一套求解边界条件问题的具体方法. 算例讨论了几种典型边界条件问 题,描述了热黏弹性材料的蠕变和松弛特征,体现了这种辛方法的有效性.     

A reciprocity theorem in linear gradient elasticity and the corresponding Saint-Venant principle

In many practical applications of nanotechnology and in microelectromechanical devices, typical structural components are in the form of beams, plates, shells and membranes. When the scale of such components is very small, the material microstructural lengths become important and strain gradient elasticity can provide useful material modelling. In addition, small scale beams and bars can be used as test specimens for measuring the lengths that enter the constitutive equations of gradient elasticity. It is then useful to be able to apply approximate solutions for the extension, shear and flexure of slender bodies. Such approach requires the existence of some form of the Saint-Venant principle. The present work presents a statement of the Saint-Venant principle in the context of linear strain gradient elasticity. A reciprocity theorem analogous to Betti’s theorem in classic elasticity is provided first, together with necessary restrictions on the constitutive equations and the body forces. It is shown that the order of magnitude of displacements are in accord with the Sternberg’s statement of the Saint-Venant principle. The cases of stretching, shearing and bending of a beam were examined in detail, using two-dimensional finite elements. The numerical examples confirmed the theoretical results.     

Closed form stress distribution in 2D elasticity for all boundary conditions

This paper applies a Hamiltonian method to study analytically the stress dis- tributions of orthotropic two-dimensional elasticity in(x,z)plane for arbitrary boundary conditions without beam assumptions.It is a method of separable variables for partial differential equations using displacements and their conjugate stresses as unknowns.Since coordinates(x,z)can not be easily separated,an alternative symplectic expansion is used. Similar to the Hamiltonian formulation in classical dynamics,we treat the x coordinate as time variable so that z becomes the only independent coordinate in the Hamiltonian ma- trix differential operator.The exponential of the Hamiltonian matrix is symplectic.There are homogenous solutions with constants to be determined by the boundary conditions and particular integrals satisfying the loading conditions.The homogenous solutions consist of the eigen-solutions of the derogatory zero eigenvalues(zero eigen-solutions) and that of the well-behaved nonzero eigenvalues(nonzero eigen-solutions).The Jordan chains at zero eigenvalues give the classical Saint-Venant solutions associated with aver- aged global behaviors such as rigid-body translation,rigid-body rotation or bending.On the other hand,the nonzero eigen-solutions describe the exponentially decaying localized solutions usually ignored by Saint-Venant's principle.Completed numerical examples are newly given to compare with established results.     

Stress singularity analysis of anisotropic multi-material wedges under antiplane shear deformation using the symplectic approach

Symplectic approach has emerged a popular tool in dealing with elasticity problems especially for those with stress singularities. However, anisotropic material problem under polar coordinate system is still a bottleneck. This paper presents a subfield method coupled with the symplectic approach to study the anisotropic material under antiplane shear deformation. Anisotropic material around wedge tip is considered to be consisted of many subfields with constant material properties which can be handled by the symplectic approach individually. In this way, approximate solutions of the stress and displacement can be obtained. Numerical examples show that the present method is very accurate and efficient for such wedge problems. Besides, this paper has extended the application of the symplectic approach and provides a new idea for wedge problems of anisotropic material.     

On new symplectic elasticity approach for exact bending solutions of rectangular thin plates with two opposite sides simply supported

This paper presents a bridging research between a modeling methodology in quantum mechanics/relativity and elasticity. Using the symplectic method commonly applied in quantum mechanics and relativity, a new symplectic elasticity approach is developed for deriving exact analytical solutions to some basic problems in solid mechanics and elasticity which have long been bottlenecks in the history of elasticity. In specific, it is applied to bending of rectangular thin plates where exact solutions are hitherto unavailable. It employs the Hamiltonian principle with Legendre’s transformation. Analytical bending solutions could be obtained by eigenvalue analysis and expansion of eigenfunctions. Here, bending analysis requires the solving of an eigenvalue equation unlike in classical mechanics where eigenvalue analysis is only required in vibration and buckling problems. Furthermore, unlike the semi-inverse approaches in classical plate analysis employed by Timoshenko and others such as Navier’s solution, Levy’s solution, Rayleigh–Ritz method, etc. where a trial deflection function is pre-determined, this new symplectic plate analysis is completely rational without any guess functions and yet it renders exact solutions beyond the scope of applicability of the semi-inverse approaches. In short, the symplectic plate analysis developed in this paper presents a breakthrough in analytical mechanics in which an area previously unaccountable by Timoshenko’s plate theory and the likes has been trespassed. Here, examples for plates with selected boundary conditions are solved and the exact solutions discussed. Comparison with the classical solutions shows excellent agreement. As the derivation of this new approach is fundamental, further research can be conducted not only on other types of boundary conditions, but also for thick plates as well as vibration, buckling, wave propagation, etc.     

基于辛体系的正交各向异性地基梁解析解

基于二维弹性理论, 利用Hellinger-Reissner变分原理, 通过引入对偶变量, 推导 了双参数地基上正交各向异性梁平面应力问题的辛对偶方程组; 采用分离变量法和本征展 开方法, 将原问题归结为求解零本征值本征解和非零本征值本征解, 得到了适用于任意横纵 比的梁的解析解. 由于在求解过程中不需要事先人为地选取试函数, 而是从梁的基本方程出 发, 直接利用数学方法求出问题的解, 使得问题的求解更加合理化. 其中, 地基对梁的力学 行为的影响看作是侧边边界条件, 类似于外载, 可通过零本征解的线性展开来评价, 非零本 征值本征解对应圣维南原理覆盖的部分. 还利用哈密顿变分原理, 给出了两端固支梁的 一种新的改进边界条件. 编程计算了细梁和深梁等算例, 研究了地基上梁的变形沿着厚度方 向的变化特性, 验证了辛方法的有效性.     

End effects for anti-plane shear deformations of sandwich structures

The purpose of this research is to further investigate the effects of material inhomogeneity and the combined effects of material inhomogeneity and anisotropy on the decay of Saint-Venant end effects. Saint-Venant decay rates for self-equilibrated edge loads in symmetric sandwich structures are examined in the context of anti-plane shear for linear anisotropic elasticity. The problem is governed by a second-order, linear, elliptic, partial differential equation with discontinuous coefficients. The most general anisotropy consistent with a state of anti-plane shear is considered, as well as a variety of boundary conditions. Anti-plane or longitudinal shear deformations are one of the simplest classes of deformations in solid mechanics. The resulting deformations are completely characterized by a single out-of-plane displacement which depends only on the in-plane coordinates. They can be thought of as complementary deformations to those of plane elasticity. While these deformations have received little attention compared with the plane problems of linear elasticity, they have recently been investigated for anisotropic and inhomogeneous linear elasticity. In the context of linear elasticity, Saint-Venant's principle is used to show that self-equilibrated loads generate local stress effects that quickly decay away from the loaded end of a structure. For homogeneous isotropic linear elastic materials this is well-documented. Self-equilibrated loads are a class of load distributions that are statically equivalent to zero, i.e., have zero resultant force and moment. When Saint-Venant's principle is valid, pointwise boundary conditions can be replaced by more tractable resultant conditions. It is shown in the present study that material inhomogeneity significantly affects the practical application of Saint-Venant's principle to sandwich structures.     

Fracture analysis of mode III crack problems for the piezoelectric bimorph

In this paper, a symplectic method based on the Hamiltonian system is proposed to analyze the interfacial fracture in the piezoelectric bimorph under anti-plane deformation. A set of Hamiltonian governing equations is derived from the Hamiltonian function by introducing dual variables of generalized displacements and stresses which can be expanded in series in terms of the symplectic eigensolutions. With the aid of the adjoint symplectic orthogonality, coefficients of the series are determined by the boundary conditions along the crack faces and along the external geometry. The stress\electric displacement intensity factors and energy release rates (G) directly relate to the first few terms of the nonzero eigenvalue solutions. The two ideal crack boundary conditions, namely the electrically impermeable and permeable crack assumptions, are considered. Numerical examples including the complex mixed boundary conditions are considered to show fracture behaviors of the interface crack and discuss the influencing factors.     

Exact Analysis of Curved Beams and Arches with Arbitrary End Conditions: A Hamiltonian State Space Approach

An exact analysis of stress and displacement fields in curved beams and arches subjected to inplane loads is conducted, with emphasis on the end effects. The material considered is cylindrically orthotropic, including transverse isotropy and isotropy as special cases. On the basis of the Hamiltonian state space approach, exact solutions that satisfy any combination of the fixed, free, and sliding-contact end conditions in a pointwise fashion are determined through symplectic eigenfunction expansion. The study allows us to evaluate the conventional solutions based on the elementary theory of curved beams and plane elasticity under simplifying assumptions, thereby, to assess the St. Venant principle as applied to this class of problems.     

On the role of symmetry in Saint-Venant’s problem

In the relaxed Saint-Venant’s elastic problem, in virtue of Saint-Venant’s Postulate, the pointwise assignments of tractions at cylinder plane ends are replaced by the assignments of the corresponding resultant forces and moments. The solution indeterminacy so introduced is traditionally overcome by postulating that some specific features characterize the elastic state. In this work a relaxed incremental equilibrium problem is posed for a heterogenous anisotropic cylinder, whose tangent elasticity tensor field possesses the usual major and minor symmetries, is positive definite, independent from the axial coordinate and endowed with a plane of elastic symmetry orthogonal to the cylinder axis. Symmetry has been consistently employed to formulate the basic problems of extension, bending, torsion and flexure as symmetric and antisymmetric problems respectively. It is shown that Saint-Venant’s postulate, momentum balance and symmetry are sufficient, without resorting to any a priori assumption, to determine the general form of the displacement field and to remove the solution indeterminacy.     

Corner stress singularities in an FGM thin plate

Describing the behaviors of stress singularities correctly is essential for obtaining accurate numerical solutions of complicated problems with stress singularities. This analysis derives asymptotic solutions for functionally graded material (FGM) thin plates with geometrically induced stress singularities. The classical thin plate theory is used to establish the equilibrium equations for FGM thin plates. It is assumed that the Young’s modulus varies along the thickness and Poisson’s ratio is constant. The eigenfunction expansion method is employed to the equilibrium equations in terms of displacement components for an asymptotic analysis in the vicinity of a sharp corner. The characteristic equations for determining the stress singularity order at the corner vertex and the corresponding corner functions are explicitly given for different combinations of boundary conditions along the radial edges forming the sharp corner. The non-homogeneous elasticity properties are present only in the characteristic equations corresponding to boundary conditions involving simple support. Finally, the effects of material non-homogeneity following a power law on the stress singularity orders are thoroughly examined by showing the minimum real values of the roots of the characteristic equations varying with the material properties and vertex angle.     

Variational formulation of a simplified strain gradient elasticity theory and its application to a pressurized thick-walled cylinder problem

A detailed variational formulation is provided for a simplified strain gradient elasticity theory by using the principle of minimum total potential energy. This leads to the simultaneous determination of the equilibrium equations and the complete boundary conditions of the theory for the first time. To supplement the stress-based formulation, the coordinate-invariant displacement form of the simplified strain gradient elasticity theory is also derived anew. In view of the lack of a consistent and complete formulation, derivation details are included for the tutorial purpose. It is shown that both the stress and displacement forms of the simplified strain gradient elasticity theory obtained reduce to their counterparts in classical elasticity when the strain gradient effect (a measure of the underlying material microstructure) is not considered. As a direct application of the newly obtained displacement form of the theory, the problem of a pressurized thick-walled cylinder is analytically solved. The solution contains a material length scale parameter and can account for microstructural effects, which is qualitatively different from Lamé’s solution in classical elasticity. In the absence of the strain gradient effect, this strain gradient elasticity solution reduces to Lamé’s solution. The numerical results reveal that microstructural effects can be large and Lamé’s solution may not be accurate for materials exhibiting significant microstructure dependence.     

New assessment on the Saint-Venant solutions for functionally graded beams

A symplectic approach is proposed to investigate the Saint-Venant problem of functionally graded beams with Young's modulus varying exponentially in the axial direction and constant Poisson radio. A matrix state equation is derived with a shift-Hamiltonian operator matrix whose particular eigenvalues are proved to compose the basic solutions of the Saint-Venant problem. The present analyses demonstrate that the Saint-Venant solutions under simple extension and pure bending can be derived using either the direct expansion method or the rigid motion removing method.     

含裂纹纳米板振动问题的哈密顿体系方法

基于裂纹处范德华力效应,采用非局部弹性理论构造纳米板模型,并通过导入哈密顿体系建立含裂纹纳米板振动问题的对偶正则控制方程组。在全状态向量表示的哈密顿体系下,将含裂纹纳米板的固有频率和振型问题归结为广义辛本征值和本征解问题。利用哈密顿体系具有的辛共轭正交关系,得到问题解的级数解析表达式。结合边界条件,得到固有频率与辛本征值的代数方程关系式,进而直接给出固有频率的表达式。数值结果表明,非局部尺寸参数和裂纹长度对纳米板振动的各阶固有频率有直接的影响。对比表明,辛方法是准确且可靠的,可为工程应用提供依据。     

Symplectic solution system for reissner plate bending

Based on the Hellinger-Reissner variatonal principle for Reissner plate bendingand introducing dual variables, Hamiltonian dual equations for Reissner plate bending werepresented. Therefore Hamiltonian solution system can also be applied to Reissner platebending problem, and the transformation from Euclidian space to symplectic space and fromLagrangian system to Hamiltonian system was realized. So in the symplectic space whichconsists of the original variables and their dual variables, the problem can be solved viaeffective mathematical physics methods such as the method of separation of variables andeigenfunction-vector expansion. All the eigensolutions and Jordan canonical formeigensolutions for zero eigenvalue of the Hamiltonian operator matrix are solved in detail, and their physical meanings are showed clearly. The adjoint symplectic orthonormal relation of the eigenfunction vectors for zero eigenvalue are formed. It is showed that the alleigensolutions for zero eigenvalue are basic solutions of the Saint-Venant problem and theyform a perfect symplectic subspace for zero eigenvalue. And the eigensolutions for nonzeroeigenvalue are covered by the Saint-Venant theorem. The symplectic solution method is notthe same as the classical semi-inverse method and breaks through the limit of the traditional semi-inverse solution. The symplectic solution method will have vast application.     

哈密顿体系下机电耦合问题的基本方程

利用多变量变分原理,针对具有机电耦合特性的压电材料,在哈密顿体系下推导出了机电耦合问题的对偶方程．求出了问题的零本征解和非零本征解的一般表达式。并就“二维压电平板对边受均布载荷”这一具体问题进行了完整的分析计算,得出了应力和电位移的具体解答,结果与一般弹性力学所得的结论相吻合,而且精度极高。表明哈密顿体系适用于压电材料的力学分析。     

Saint-Venant problem of three-dimensional linear viscoelasticity in the Hamiltonian system

The traditional Saint-Venant problem of three-dimensional viscoelasticity is discussed under the Hamiltonia system with the use of the Laplace integral transformation, and the original problem is transformed into finding eigenvalues and eigenvectors of the Hamiltonia operator matrix. Since local effect near the boundary is usually neglected, all solutions of Saint-Venant problems can be obtained directly by the combinations of zero eigenvectors. Moreover, the adjoint relationships of the symplectic orthogonality of zero eigenvectors in the Laplace domain are generalized to the time domain. Therefore the problem can be discussed directly in the eigenvector space of the time domain, and the iterative application of Laplace transformation is not needed. Simply by applying the adjoint relationships of the symplectic orthogonality, an effective method for boundary condition is given. Based on this method, some typical examples are discussed, in which the whole character of total creep and relaxation of viscoelasticity is clearly revealed.     

Saint-Venant problem for solids with helical anisotropy

We discuss the solution of Saint-Venant’s problem for solids with helical anisotropy. Here the governing relations of the theory of elasticity in terms of displacements were presented using the helical coordinate system. We proposed an approach to construct elementary Saint-Venant solutions using integration of ordinary differential equations with variable coefficients in the case of a circular cylinder with helical anisotropy. Elementary solutions correspond to problems of extension, of torsion, of pure bending and of bending of shear force. The solution of the problem is obtained using small parameter method for small values of twist angle and numerically for arbitrary values. Numeric results correspond to problems of extension–torsion. Dependencies of the stiffness matrix (in dimensionless form) on angle between the tangent to the helical coil and the axis of the cylinder corresponding to stiffness on stretching and torsion are illustrated graphically for different values of material and geometrical parameters.     

Torsion of cylindrically orthotropic elastic circular bars with radial inhomogeneity: some exact solutions and end effects

Torsion of elastic circular bars of radially inhomogeneous, cylindrically orthotropic materials is studied with emphasis on the end effects. To examine the conjecture of Saint-Venant’s torsion, we consider torsion of circular bars with one end fixed and the other end free on which tractions that results in a pure torque are prescribed arbitrarily over the free end surface. Exact solutions that satisfy the prescribed boundary conditions point by point over the entire boundary surfaces are derived in a unified manner for cylindrically orthotropic bars with or without radial inhomogeneity and for their counterparts of Saint-Venant’s torsion. Stress diffusion due to the end effect is examined in the light of the exact solutions.