坐标变换在应力分析中的应用

<正> 一点应力状态分析是弹性力学中的重要内容,国内外不少教材都借助微分体的平衡推导任一法向的斜面中应力分量和三坐标面内应力分量面的关系,推导较繁,学生亦难于记住.若利用坐标变换则可简化推     

用材料力学解的切应力求解弹性力学问题

分析了切应力与正应力的关系,讨论了导出切应力公式的条件,提出按切应力求解弹性力学问题的方法.证明凡是σ_y与x无关的梁或者已知切应力为零的问题,按切应力求解都是可行的.用这种方法求解比传统方法方便,运算简单.该文给出了用切应力求解弹性力学平面问题的两个算例.     

弹性力学平面问题例说解的唯一性定理

解的唯一性定理是用逆解法或半逆解法求解弹性力学问题的理论依据,在此用应力函数法、应力法、应力和函数法求解弹性力学平面问题,让学生切实、深入地理解解的唯一性定理的内在含义,丰富和扩大弹性力学的解题方法和应用范围。     

EXAMPLES ARE GIVEN TO ILLUSTRATE THE UNIQUENESS LAW OF SOLUTIONS FOR PLANE PROBLEMS OF ELASTICITY

解的唯一性定理是用逆解法或半逆解法求解弹性力学问题的理论依据，在此用应力函数法、应力法、应力和函数法求解弹性力学平面问题，让学生切实、深入地理解解的唯一性定理的内在含义，丰富和扩大弹性力学的解题方法和应用范围。     

斜面应力公式的适用条件

斜面应力公式,即一点应力的坐标旋转变换公式,是材料力学和弹性力学里最常用的公式之一,并广泛用于固体力学和工程设计中。一个让学生感觉例外的典型例子是含V形切口的薄板在切口尖端的应力情况,如果利用斜面应力公式和切口面的自由面条件,就会得出切口尖点处于零应力状态的结果,而这与线弹性断裂力学给出切口处应力趋于无限大的结果不符。为消除这一疑虑,考察了尖端应力的特性,指出:只有在过一点的各个斜面上的应力是单值连续的情况下,斜面应力公式才能适用,此时,该点的全部应力分量组成应力张量。在V形切口的尖端、裂纹尖端,自由面与顺其延伸至介质内侧的面上的应力不同,应力在该面上就不是单值连续的,该点的应力状态就不能用张量表示,斜面应力公式在切口或裂纹尖端就不适用了。     

CONDITIONS FOR THE APPLICATION OF THE STRESS TRANSFORMATION FORMULA 1)

斜面应力公式,即一点应力的坐标旋转变换公式,是材料力学和弹性力学里最常用的公式之一,并广泛用于固体力学和工程设计中。一个让学生感觉例外的典型例子是含V 形切口的薄板在切口尖端的应力情况,如果利用斜面应力公式和切口面的自由面条件,就会得出切口尖点处于零应力状态的结果,而这与线弹性断裂力学给出切口处应力趋于无限大的结果不符。为消除这一疑虑,考察了尖端应力的特性,指出：只有在过一点的各个斜面上的应力是单值连续的情况下,斜面应力公式才能适用,此时,该点的全部应力分量组成应力张量。在V形切口的尖端、裂纹尖端,自由面与顺其延伸至介质内侧的面上的应力不同,应力在该面上就不是单值连续的,该点的应力状态就不能用张量表示,斜面应力公式在切口或裂纹尖端就不适用了。     

双层球壳结构准静态线性化学弹性问题的位移势函数解

本文研究了各向同性固体的化学-力学耦合问题,在传统化学弹性理论描述的扩散-变形耦合关系基础上,进一步考虑了化学反应与固体变形的相互作用关系,发展了等温状态下固体-扩散-反应-变形耦合的线性化学弹性理论,拓展了化学弹性力学的应用范围．该理论能够同时描述固体内介质扩散和固体与介质之间化学反应两个不同时间尺度的化学过程,并给出由此引起的弹性范围内的应变和应力．为应用该模型求解具体化学弹性问题,本文通过构造扩散-反应位移势函数来获得位移特解形式,再与齐次Lamé方程通解叠加获得完整解;针对反应控制问题,引入化学弹性准静态假设,将反应-扩散-变形全耦合的瞬态过程分解为两个可解耦的相继过程,从而获得相应位移解．基于此解法,本文获得了反应控制的双层球壳结构化学弹性问题的解析解,并分析了化学反应、几何结构和弹性模量对应力分布的影响．     

平面弹性问题的位移-应力混合重心插值配点法

提出数值分析平面弹性问题的位移-应力混合重心插值配点法。将弹性力学控制方程表达为位移和应力的耦合偏微分方程组,采用重心插值近似未知量,利用重心插值微分矩阵得到平面问题控制方程的矩阵形式离散表达式。使用重心插值离散位移和应力边界条件,采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,应用最小二乘法求解过约束方程组,得到平面弹性问题位移和应力数值解。数值算例结果表明,重心Lagrange插值方法的计算精度可达到10~(-10)量级。位移-应力混合重心插值配点法的计算公式简单、程序实施方便,是一种高精度的无网格数值分析方法。     

平面弹性力学问题的离散元法

根据离散元的基本原理,基于变形体的理论提出了适用于平面弹性力学问题的界面位移、应变和应力模式,建立了求解平面弹性力学问题的离散元方程和相应的迭代求解方法.通过界面位移可以简洁地将位移和力的边界条件引入离散系统的控制方程,也可以方便地求解节点位移.数值算例表明,与具有相同网格的有限元结果相比,离散元能同时给出精度相对较高的应力解和精度相当的位移解.     

考虑流固耦合效应的某飞行器力学性能分析

考虑流固耦合效应,研究了飞行器在服役状态下的力学特性。本文采用了一种基于CFD/CSD(Computational Fluid Dynamics/Computational Structure Dynamics)耦合的高精度通用静气弹求解器,该求解器可以同时分析飞行器结构的表面气动性能和结构的力学性能,对某一翼身组合体模型进行了静气弹分析,得出弹性飞行器和刚性飞行器的气动压力变化以及结构的应力分布。建模时流体和结构模型都使用翼身组合模型,且都考虑机身的影响。结果表明:发生静气弹变形后,飞行器的升力系数、阻力系数、力矩系数都相应下降;机翼翼根附近会产生较大的应力分布;当飞行器尤其是机翼表面受到弯矩、扭矩以及气动压力的共同作用时,翼根附近会出现复杂的应力分布状态。这说明在静气弹性分析中,机身的建模也是很有必要的。     

Size-dependent effect on biaxial and shear nonlinear buckling analysis of nonlocal isotropic and orthotropic micro-plate based on surface stress and modified couple stress theories using differential quadrature method

The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress,the modified couple stress theory(MCST),and the nonlocal elasticity theories using the differential quadrature method(DQM)is presented.Main advantages of the MCST over the classical theory(CT)are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter.Based on the nonlinear von K′arm′an assumption,the governing equations of equilibrium for the micro-classical plate considering midplane displacements are derived based on the minimum principle of potential energy.Using the DQM,the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained.Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature.A parametric study is conducted to show the effects of the aspect ratio,the side-to-thickness ratio,Eringen’s nonlocal parameter,the material length scale parameter,Young’s modulus of the surface layer,the surface residual stress,the polymer matrix coefficients,and various boundary conditions on the dimensionless uniaxial,biaxial,and shear critical buckling loads.The results indicate that the critical buckling loads are strongly sensitive to Eringen’s nonlocal parameter,the material length scale parameter,and the surface residual stress effects,while the effect of Young’s modulus of the surface layer on the critical buckling load is negligible.Also,considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate.The results show that the critical biaxial buckling load increases with an increase in G12/E2and vice versa for E1/E2.It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude.Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios,it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.     

Near field behavior of in-plane crack extension in nonlinear incompressible material

The nonlinear theory of finite elasticity is applied to obtain the in-plane displacement and stresses in the immediate vicinity of the crack tip. Incompressibility, homogeneity, elasticity and isotropy are assumed for the material while the resultant shear stress and shear strain are assumed to follow a nonlinear hardening/softening behavior. The system of governing differential equations becomes nonelliptical when the strains are sufficiently large.     

宽翼缘梁非线性温差自应力级数解

为了客观反映宽翼缘梁在非线性温差作用下的自应力分布特点,克服按梁理论计算温度自应力时的不足,本文从弹性力学的分析方法入手,按平面应力状态分析翼缘板和腹板的温度应力。通过翼缘板与腹板连接处的变形协调条件及平衡条件建立补充方程,求解艾瑞应力函数中的积分常数,推导了翼缘板与腹板的纵、横向正应力、剪应力及位移分量解析式。导出的翼缘板纵向正应力公式可自动考虑宽翼缘梁沿横截面宽度自应力不均匀分布的特点。对一宽翼缘T梁的计算表明,当翼缘板相对于腹板发生温差变化时,沿其宽度的纵向温度自应力分布很不均匀,在翼缘板根部自应力较大而在悬臂端则显著减小。按通常基于梁理论的温度自应力计算方法,无法反映这种应力分布规律。     

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 in plates of finite thickness

Application of the plane theory of elasticity to planar crack or angular corner geometries leads to the concept of stress singularity and stress intensity factor, which are the cornerstone of contemporary fracture mechanics. However, the stress state near an actual crack tip or corner vertex is always three-dimensional, and the meaning of the results obtained within the plane theory of elasticity and their relation to the actual 3D problems is still not fully understood. In particular, it is not clear whether the same stress field as found from the well-known 2D solutions of the theory of elasticity do describe the corresponding stress components in a plate made of a sufficiently brittle material and subjected to in-plane loading, and what effect the plate thickness has. In the present study we adopt, so called, first order plate theory to attempt to answer these questions. New features of the elastic solutions obtained within this theory are discussed and compared with 2D analytical results and experimental studies as well as with 3D numerical simulations.     

Nonlinear resonance responses of geometrically imperfect shear deformable nanobeams including surface stress effects

In this work, a thorough investigation is presented into the nonlinear resonant dynamics of geometrically imperfect shear deformable nanobeams subjected to harmonic external excitation force in the transverse direction. To this end, the Gurtin–Murdoch surface elasticity theory together with Reddy’s third-order shear deformation beam theory is utilized to take into account the size-dependent behavior of nanobeams and the effects of transverse shear deformation and rotary inertia, respectively. The kinematic nonlinearity is considered using the von Kármán kinematic hypothesis. The geometric imperfection as a slight curvature is assumed as the mode shape associated with the first vibration mode. The weak form of geometrically nonlinear governing equations of motion is derived using the variational differential quadrature (VDQ) technique and Lagrange equations. Then, a multistep numerical scheme is employed to solve the obtained governing equations in order to study the nonlinear frequency–response and force–response curves of nanobeams. Comprehensive studies into the effects of initial imperfection and boundary condition as well as geometric parameters on the nonlinear dynamic characteristics of imperfect shear deformable nanobeams are carried out through numerical results. Finally, the importance of incorporating the surface stress effects via the Gurtin–Murdoch elasticity theory, is emphasized by comparing the nonlinear dynamic responses of the nanobeams with different thicknesses.     

Automated Stress Separation Along Stress Trajectories

A procedure for the separation of principal stresses in automated photoelasticity is presented. It is based on the integration of indefinite equations of equilibrium along stress trajectories, also known as Lamè–Maxwell equations. A new algorithm for precise and reliable stress trajectory calculation, which is an essential feature of the procedure, has also been developed. Automated stress separation is carried out along stress trajectories starting from free boundaries. Experimental tests were performed on a disc in diametral compression and on a ring with internally applied pressure. Full-field principal stress values were obtained and results were compared with those from the theory of elasticity and with those obtained from the classical shear difference method. It was shown that the proposed method is more accurate and less affected by the presence of residual stresses or experimental errors at the boundaries than the shear difference method. In addition, the method requires little human interaction and is therefore well-suited for automated photoelasticity.     

Axisymmetric bending analysis of thick functionally graded circular plates using fourth-order shear deformation theory

In the present article, axisymmetric bending and stretching of functionally graded (FG) circular plates subjected to uniform transverse loading based on fourth-order shear deformation plate theory (FOST) have been studied. Using a fourth-order shear deformation theory, the solutions for deflection and rotation functions of FG plates are presented in terms of the corresponding quantities for a homogeneous plate using the classical plate theory (CPT), from which solutions one can easily obtain the FOST solutions for axisymmetric bending of FG circular plates. It is assumed that the effective mechanical properties of the functionally graded plates through the thickness are continuous functions of the volume fractions of the constituent parts which are themselves defined by a power-law function. Numerical results for maximum deflection and shear stress are presented for various percentages of ceramic–metal volume fractions. These results are also compared with those obtained from the first-order shear deformation plate theory of Mindlin (FST), the third-order shear deformation plate theory of Reddy (TST) as well as the exact three-dimensional elasticity solution. It is found that although the maximum deflections obtained using FOST and TST are close to each other, the through-thickness shear stress is predicted more accurately by the FOST formulation than by the TST.     

A BEM solution to transverse shear loading of composite beams

In this paper a boundary element method is developed for the solution of the general transverse shear loading problem of composite beams of arbitrary constant cross-section. The composite beam consists of materials in contact, each of which can surround a finite number of inclusions. The materials have different elasticity and shear moduli with same Poisson’s ratio and are firmly bonded together. The analysis of the beam is accomplished with respect to a coordinate system that has its origin at the centroid of the cross-section, while its axes are not necessarily the principal ones. The transverse shear loading is applied at the shear centre of the cross-section, avoiding in this way the induction of a twisting moment. Two boundary value problems that take into account the effect of Poisson’s ratio are formulated with respect to stress functions and solved employing a pure BEM approach, that is only boundary discretization is used. The evaluation of the transverse shear stresses is accomplished by direct differentiation of these stress functions, while both the coordinates of the shear center and the shear deformation coefficients are obtained from these functions using only boundary integration. Numerical examples with great practical interest are worked out to illustrate the efficiency, the accuracy and the range of applications of the developed method. The accuracy of the proposed shear deformation coefficients compared with those obtained from a 3-D FEM solution of the ‘exact’ elastic beam theory is remarkable.