三向应力Mohr圆的真实构成及剪应力作用方向的确定

三向应力Mohr圆的构成在传统上是借助公式推证而得，并以平面图形来表示，缺乏三维的真实感和直观性。在应力应变分析中，对于平面应力问题，可以通过平面应力摩尔圆确定过一点不同斜面上的应力分量及其作用方向。对于三维问题，利用摩尔圆图解法可以确定某一斜面上的正应力和剪应力的数值，但不能表示剪应力的作用方向。剪应力的作用方向需要通过另外的图解方法来确定。本文分别从坐标系旋转和数值计算的角度解释了三向应力Mohr圆的构成过程，形象地说明了Mohr圆的物理本质。针对三向应力Mohr圆不能表示剪应力作用方向的问题，通过矢量运算，给出了剪应力作用方向的确定公式。     

内压椭球壳塑性变形的发生部位与扩展过程分析

对内压椭球的就发布作了图形描述，指出椭球壳上一点的应力状态只与壳体的轴比λ和该点在壳体上的位置有关，并进一步分析了内压椭球壳的塑性变形的发生部位及其扩展过程。     

单搭接胶接接头弹性应力的二维分析方法

提出一种能得到单搭接胶接接头二维弹性应力解析解的的一般分析方法.该方法从纵向正应力沿厚度线性分布的理论假设出发,用完整的几何方程与本构方程,在考虑了被粘物与胶层的剪应力和剥离正应力沿厚度的变化,而且严格满足包括搭接区端头胶层剪应力为零的所有边界条件下,能够准确地预测接头任意点的应力状态,并适用于被粘物/胶层几何和材料属性的任何组合.通过与以往的解析解和有限元分析结果对比表明:不仅胶层平均剪应力最大值发生在距胶层端头很近的位置处,而且平均剥离正应力最大值位置也距离接头端头很近.这与以往绝大多数理论解预测"胶层剪应力和剥离正应力最大值都发生在端头处"的结论不同;该解析解有很好的精度和较广的适用范围.     

绕径旋转横观各向同性圆球壳的变形和应力分析

本文用三维弹性力学理论研究绕径旋转横观各向同性圆球壳的变形和应力分布，给出任意厚度圆球壳的位移和应务的封闭形式解以及数值结果。研究表明球壳上的位移和应力大小与旋转角速度的平方成正比，球面内的正应力远大于径向正应力和剪应力。     

一种预测层合板层间剪应力的新后处理方法

层合板是航空航天领域典型的承力构件,过大的层间应力是导致其分层失效的主要原因.准确的层间应力预测往往依赖于三维平衡方程后处理方法(TPM).然而,该方法需要计算面内应力的一阶导,使得基于C0型板理论构造的线性单元无法使用TPM计算横向剪应力.本文在三维平衡方程后处理方法的基础上,提出了一种新后处理方法(NPM).新后处理方法通过虚功等效法消除了三维平衡方程后处理方法中产生的位移参数的高阶导.基于提出的新后处理方法和C0型板理论,仅需使用线性单元就可以预测层合板的横向剪应力.为了验证所提方法的有效性,本文基于修正锯齿理论(RZT)和所提方法构造了一种C0连续的三节点三角形线性板单元.数值算例表明,所提方法和三维平衡方程后处理方法具有相同的计算精度,提出的板单元能够准确高效地预测层合板的横向剪应力.此外,所提方法便于结合现有的有限元商用软件使用,基于商用软件中板壳单元获得的节点位移,使用新后处理方法极易获得准确的层间剪应力.     

裂纹端部细短纤维的应力分析

基于裂纹端部存在与其裂纹面相垂直的二相细短纤维分析模型，采用叠加原理推导了求解纤维表面应力分布函数的积分方程，通过简化得到了该方程的解析表达显式，该积分方程的特征值方程是纤维几何参数，材料常数以及纤维相对于裂纹位置的相关函数，当材料参数不满足特征方程时，积分方程将具有唯一解；并借助数值方法，给出了纤维剪应力分布算例，和纤维对应力强度因子的影响。     

谷粒的热——湿应力分析

针对谷物烘干问题,本文研究弹性及粘弹性圆球和椭球的热,湿应力。对圆球分析表明,热应力峰值远比湿应力峰值出现时间早,在假定材料为热,湿流变简单的Maxwell粘弹性体条件下,本文导出了圆球的热,湿应力分布公式,并应用有限元法研究了椭球的热,湿应力,文中以无量纲形成给出了圆球中心和表面应力随时间的 变化情况,以及椭球中心轴向应力峰值随椭球长,短半轴之比的变化情况。     

CONDITIONS FOR THE APPLICATION OF THE STRESS TRANSFORMATION FORMULA 1)

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

斜面应力公式的适用条件

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

曲梁剪应力和径向应力的积分方程解

直接对曲梁剪应力的积分方程进行求解。得到剪应力和径向应力计算的一般公式．最后给出计算实例．     

DYNAMIC STRESS FIELD AROUND THE MODE Ⅲ CRACK TIP IN AN ORTHOTROPIC FUNCTIONALLY GRADED MATERIAL

The problem of a Griffith crack in an unbounded orthotropic functionally graded material subjected to antipole shear impact was studied. The shear moduli in two directions of the functionally graded material were assumed to vary proportionately as definite gradient. By using integral transforms and dual integral equations, the local dynamic stress field was obtained. The results of dynamic stress intensity factor show that increasing shear moduli’s gradient of FGM or increasing the shear modulus in direction perpendicular to crack surface can restrain the magnitude of dynamic stress intensity factor.     

Nonlinear vibration of viscoelastic sandwich plates under narrow-band random excitations

The effect of the narrow-band random excitation on the non-linear response of sandwich plates with an incompressible viscoelastic core is investigated. To model the core, both the transverse shear strains and rotations are assumed to be moderate and the displacement field in the thickness direction is assumed to be linear for the in-plane components and quadratic for the out-of-plane components. In connection to the moderate shear strains considered for the core, a non-linear single-integral viscoelastic model is also used for constitutive modeling of the core. The fifth-order perturbation method is used together with the Galerkin method to transform the nine partial differential equations to a single ordinary integro-differential equation. Converting the lower-order viscoelastic integral term to the differential form, the fifth-order method of multiple scale is applied together with the method of reconstitution to obtain the stochastic phase-amplitude equations. The Fokker–Planck–Kolmogorov equation corresponding to these equations is then solved by the finite difference method, to determine the probability density of the response. The variation of root mean square and marginal probability density of the response amplitude with excitation deterministic frequency and magnitudes are investigated and the bimodal distribution is recognized in narrow ranges of excitation frequency and magnitude.     

Measurement of Torsionally Induced Shear Stresses in Shrink-Fit Assemblies

Shear stresses along the shaft/hub interface in shrink-fit components, generated by torsional loads, can drive premature failure through fretting mechanisms. It is difficult to numerically predict these shear stresses, and the associated circumferential slip along the shaft/hub interface, due to uncertainties in frictional behaviour and the presence of steep stress gradients which can cause meshing problems. This paper attempts to provide validation of a numerical modelling methodology, based on finite element analysis, so the procedure may be used with confidence in fitness-for-purpose cases. Very few experimental techniques offer the potential to make measurements of stress and residual stress interior to metallic components. Even fewer techniques provide the possibility of measuring shear stresses. This paper reports the results of neutron diffraction measurements of shear stress and residual shear stress in a bespoke test specimen containing a shrink-fit. One set of measurements was made with a torsional load ‘locked-in’. A second set of measurements was made to determine the residual shear stress when the torsional load had been applied and removed. Overall, measurement results were consistent with numerical models, but the necessity for a small test specimen to allow penetration of the neutron beam to the measurement locations meant the magnitude of shear stresses was at the limits of what could be measured experimentally.     

DYNAMIC STRESS FIELD AROUND THE MODE Ⅲ CRACK TIP IN AN ORTHOTROPIC FUNCTIONALLY GRADED MATERIAL

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ｇｒｅａｔａｔｔｅｎｔｉｏｎｓｈａｖｅｂｅｅｎｐａｉｄｔｏｔｈｅｒｅｓｅａｒｃｈｏｆＦｕｎｃｔｉｏｎａｌｌｙＧｒａｄｅｄＭａｔｅｒｉａｌｓ(ＦＧＭ).Ｆｒｏｍｔｈｅｖｉｅｗｐｏｉｎｔｓｏｆａｐｐｌｉｅｄｍｅｃｈａｎｉｃｓ,ＦＧＭａｒｅｎｏｎ＿ｈｏｍｏｇｅｎｅｏｕｓｓｏｌｉｄｓ.Ｔｈｅｎｏｎ＿ｈｏｍｏｇｅｎｅｉｔｙｏｆＦＧＭｈａｓａｇｒｅａｔｉｎｆｌｕｅｎｃｅｏｎｔｈｅｉｒｍｅｃｈａｎｉｃａｌｂｅｈａｖｉｏｒ,ｅｓｐｅｃｉａｌｌｙｗｈｅｎｔｈｅｃｏｍｐ…     

用材料力学公式计算任意截面短梁正应力时的修正项研究

对材料力学中梁的弯曲应力公式增加一修正项,以反映短梁弯剪翘曲变形对应力分布的影响。提出一种根据短梁横截面边界形状及艾瑞应力函数求解应力修正项的方法,应用弹性力学空间问题的一般理论,通过应力平衡方程、应变相容方程及应力边界条件,建立了关于任意截面短梁的应力修正项及剪应力的基本方程。在所建立的基本方程基础上,导出了矩形截面和圆形截面短梁修正应力的具体计算公式,该修正应力与均布荷载大小及弹性模量与剪切模量之比均成正比,但与截面惯性矩成反比。数值算例表明,本文方法计算的应力与通用有限元软件ANSYS计算的结果吻合良好,从而验证了本文方法及其基本公式的正确性。     

Reynolds stress modelling of rectangular open‐channel flow

A Reynolds stress model for the numerical simulation of uniform 3D turbulent open‐channel flows is described. The finite volume method is used for the numerical solution of the flow equations and transport equations of the Reynolds stress components. The overall solution strategy is the SIMPLER algorithm, and the power‐law scheme is used to discretize the convection and diffusion terms in the governing equations. The developed model is applied to a flow at a Reynolds number of 77000 in a rectangular channel with a width to depth ratio of 2. The simulated mean flow and turbulence structures are compared with measured and computed data from the literature. The computed flow vectors in the plane normal to the streamwise direction show a small vortex, called inner secondary currents, located at the juncture of the sidewall and the free surface as well as the free surface and bottom vortices. This small vortex causes a significant increase in the wall shear stress in the vicinity of the free surface. A budget analysis of the streamwise vorticity is carried out. It is found that both production terms by anisotropy of Reynolds normal stress and by Reynolds shear stress contribute to the generation of secondary currents. Copyright © 2006 John Wiley & Sons, Ltd.     

A solution for the vertical variation of stress,rather than velocity,in a three-dimensional circulation model

A simple technique is presented that allows a numerical solution to be sought for the vertical variation of shear stress as a substitute for the vertical variation of velocity in a three-dimensional hydrodynamic model. In its most general form the direct stress solution (DSS) method depends only upon the validity of an eddy viscosity relation between the shear stress and the vertical gradient of velocity. The rationale for preferring a numerical solution for shear stress to one for velocity is that shear stress tends to vary more slowly over the vertical than velocity, particularly near boundaries. Consequently, a numerical solution can be obtained much more efficiently for shear stress than for velocity. When needed, the velocity profile can be recovered from the stress profile by solving a one-dimensional integral equation over the vertical. For most practical problems this equation can be solved in closed form. Comparisons are presented between the DSS technique, the standard velocity solution technique and analytical solutions for wind-driven circulation in an unstratified, closed, rectangular channel governed by the linear equations of motion. In no case was the computational effort required by the velocity solution competitive with the DSS when a physically realistic boundary layer was included. The DSS technique should be particularly beneficial in numerical models of relatively shallow water bodies in which the bottom and surface boundary layers occupy a significant portion of the water column.     

刚性圆管中血液周期振荡流的切应力分布

本文通过求解圆管内血液振荡流的基本方程，求得圆管内血液流的压力梯度与切应力之间的关系式。在此基础上，详细讲座了圆管中轴向流速和切变率谐波的变化规律，指出流速谐波和切变率谐波的幅值都将随着谐波次数的增大而逐渐减小。为了使所得结果便于应用。文章通过管轴向中心线流速与压力梯度之间的关系式，进一步给出一种利用管轴向中心线流速计算管内切应力分布的简便方法。该方法用于检测活体血管内血液振荡流的切应力分布，具有操作简单，精度较高的优点。最后，以人体颈动脉为例，讨论血液周期振荡流的切应力的分布特性。发现在任意时刻，除了邻近管壁处切应力急剧增大到一定数值之外，沿管截面切应力分布相当均匀且接近于零，呈现出与定常流不同的切应力分布特征。     

Fast evaluation of friction forces in traction-drive contacts including thermal effects

Summary  In automotive traction drives, power is transmitted by friction forces. The friction forces result from the shear stresses developed in lubricated and highly loaded contacts between rolling bodies. Due to the kinematics of a traction drive, shear velocities occur in both the rolling direction and perpendicular to it. Due to these shear velocities and by normal pressure, the lubricant is forced to build up shear stresses. The increase of the shear stresses may be modelled by a nonlinear viscous element. The describing differential equations are coupled by the equivalent shear stress, which defines the nonlinear behaviour of the element. A fast method is described to evaluate the coupled differential equations. By using a known analytical approximation for the equivalent shear stress, the differential equations are decoupled and can be solved analytically. In an iterative procedure the equivalent shear stress is updated, and the complete solution is found. The iterative method is extended to account for thermal effects in the contact. Received 17 June 1999; accepted for publication 26 October 1999