使用焦散线法与光弹法测定三维裂纹混合型应力强度因子

本文介绍了将焦散线法与应力冻结，“解冻”技术相结合，使焦散线法用于测量三维体内部裂纹前缘应力强度因子的实验方法，并针对复杂应力状态下三维裂纹前缘的不同应变奇异场，合理地综合运用焦散线法与三维光弹法，实际测量与分离了三维裂纹前缘混合型应力强度因子。     

有限长压电层合简支板自由振动的三维精确解

基于三维弹性理论和压电理论，导出了有限长矩形压电层合简支板的动力学方程及相应的边界条件，给出了一种求解压电层合板自由振动三维精确解的方法；分析了正、逆向压电效应对层合板振动频率的影响．本文所述的方法和结果对于求解其他三维动态问题，验证、比较其他简化模型、有限元计算结果以及工程应用都有指导意义．     

The method of caustics for static plane elasticity problems

By use of the complex stress function analysis of Muskhelishvili-Kolosov and conformal mapping procedures the general governing equations of the method of caustics or shadow spot technique have been developed for optically isotropic and anisotropic materials in static plane elasticity theory. Special cases of caustics formed about cutouts, cracks, and various singular regions in static elastic stress fields are obtained upon specialization.     

Spatial decay estimates and upper bounds in elasticity for domains with unbounded cross-sections

In this paper we investigate spatial decay estimates and upper bounds for the solutions of elastic problems when the cross-sections of the three dimensional solid are semi-infinite strips. We obtain spatial decay estimates for the solutions of a static problem in the theories of homogeneous and isotropic linear elasticity and linearised elasticity. Energy bounds and some spatial decay estimates are obtained for the solutions of a dynamical problem in the case of anisotropic linear elasticity. For both kinds of problems we use the energy methods.     

Experimental determination of stress-intensity factors using caustics and photoelasticity

The optical method of caustics has been used with considerable success in recent years for determining stress-intensity factors in both static and dynamic problems. However the midplane analysis explaining the formation of transmission caustics has certain approximations that need to be examined. In this paper it is demonstrated that the midplane analysis is in good agreement with the numerical solution of the exact geometrical optical equations for mapping the light-ray path in a cracked plate. Since both these analyses are obtained by imposing a two-dimensional crack-tip stress field the sensitivity of the method to the deviations from the imposed stress field is examined next. The implications of this examination on the photoelastic technique are then discussed.     

The method of analysis of crack problem in three-dimensional non-local elasticity

ｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｈｅｏｒｅｔｉｃａｌｆｒａｍｅｏｆｎｏｎ＿ｌｏｃａｌｅｌａｓｔｉｃｉｔｙｈａｓｂｅｎｓｅｔｕｐ［１,２］．Ｍａｎｙｉｍｐｏｒｔａｎｔｒｅｓｕｌｔｓｈａｖｅｂｅｎｇｏｔｂｙｕｓｉｎｇｎｏｎ＿ｌｏｃａｌｔｈｅｏｒｙｔｏｓｏｌ...     

横向力作用下悬臂梁固定端应力分布研究

在弹性力学Hamilton体系中,利用解析法,考虑圣维南原理所覆盖的解,对横向力作用下悬臂梁固定端应力分布问题进行研究,并对计算结果进行分析。研究结果表明,辛体系解析法采用对偶的二类变量求解,能很好地处理各种复杂边界条件,并且对此类问题的分析具有优越性,计算精度较高。该方法对其他边界问题的研究也具有指导意义。     

Flow simulation on moving boundary‐fitted grids and application to fluid–structure interaction problems

We present a method for the parallel numerical simulation of transient three‐dimensional fluid–structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non‐overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time‐dependent domains. To this end, we present a technique to solve the incompressible Navier–Stokes equation in three‐dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time‐dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian–Eulerian formulation of the Navier–Stokes equations. Here the grid velocity is treated in such a way that the so‐called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well‐known MAC‐method to a staggered mesh in moving boundary‐fitted coordinates which uses grid‐dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second‐order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid–structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid–structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Circularly cylindrical layered media in plane elasticity

This paper presents an alternative efficient procedure to analyze plane elasticity problems of a circularly cylindrical layered media subject to an arbitrary point force. Based on the method of analytical continuation in conjunction with the alternating technique, the elastic fields of the three-phase media are derived. A rapidly convergent series solution which is expressed in terms of an explicit general term of the complex potential of the corresponding homogeneous problem is obtained in an elegant form. As a numerical illustration, the interfacial stresses are presented for different material combinations and for different positions of the point force.     

Some problems in boundary element programming of axisymmetric elasticity

Axisymmetric problems in elasticity can be reduced to two dimensional ones, but they are a little more complicated than plane problems. Therefore, some special problems will be encountered in the boundary element programming of axisymmetric elasticity. In this paper, the methods to treat these problems and some remarks are given according to our experience in programming. Numerical examples are presented for the checking of these treatments.     

New Solution of Axisymmetric Contact Problem of Elasticity

We consider two problems of elasticity, namely, the Boussinesq problem about the action of a lumped force on a half-space and the related problem about the interaction of the half-space with a cylindrical rigid punch with plane base. In the classical statement, these problems have singular solutions. In the Boussinesq problem, the displacement under the action of the force is infinitely large, and in the punch problem, the infinitely large variable is the pressure on the punch boundary. In the present paper, these problems are solved with the use of relations of generalized elasticity derived regarding a medium element of small but finite dimensions rather than a traditional infinitesimal element. The structure parameter of the medium contained in the solutions can be determined experimentally. The obtained generalized solutions of the problems under study are regular.     

带双裂纹的椭圆孔口问题的应力分析

利用复变方法,通过构造新的保角映射,研究了带双裂纹椭圆孔口的平面弹性问 题,得到了I型与II型裂纹问题应力强度因子的解析解.在极限情形下, 不仅可以还原经典的Griffith裂纹的结果,而且可模拟出十字裂纹和带双裂纹的圆 形孔口问题.     

Nonlinear equations of elastic deformation of plates

A method for constructing nonlinear equations of elastic deformation of plates with boundary conditions for stresses and displacements at the face surfaces in an arbitrary coordinate system is proposed. The initial three–dimensional problem of the nonlinear theory of elasticity is reduced to a one–parameter sequence of two–dimensional problems by approximating the unknown functions by truncated series in Legendre polynomials. The same unknowns are approximated by different truncated series. In each approximation, a linearized system of equations whose differential order does not depend on the boundary conditions at the face surfaces which can be formulated in terms of stresses or displacements is obtained.     

Axisymmetric bending for thick laminated circular plate under a concentrated load

Ｎｏｗａｄａｙｓ,ｔｈｅｃｕｒｒｅｎｔｔｈｅｏｒｉｅｓｏｆｐｌａｔｅｓａｎｄｓｈｅｌｌｓ,ｓｕｃｈａｓＫｉｒｃｈｈｏｆｆ’ｓｔｈｉｎｐｌａｔｅｔｈｅｏｒｙａｎｄＲｅｉｓｓｎｅｒ’ｓｍｏｄｅｒａｔｅｌｙｔｈｉｃｋｐｌａｔｅｔｈｅｏｒｙｅｔｃ．,ａｒｅｅ...     

Development of highly accurate interpolation method for mesh‐free flow simulations II. Application of CIVA method to incompressible flow simulations

This is the second report on the development of a highly accurate interpolation method, which is called cubic interpolation with volume/area (CIVA) co‐ordinates, for mesh‐free flow simulations. In this paper, the method of determining the c‐parameter of CIVA using a constant curvature condition is first considered for the two‐ and three‐dimensional cases. A computation of a three‐dimensional passive scalar advection problem is performed for accuracy verification and for comparison with widely used methods. Then, an application algorithm of the CIVA method respecting incompressible fluid simulation is presented. As the incompressible condition based on Lagrangian approaches causes problems, in this paper we consider the condition based on the conventional Eulerian approach. The CIVA‐based incompressible flow simulation algorithm enables a highly accurate simulation of many kinds of problems that have complicated geometries and involve complicated phenomena. To confirm the facts, numerical analyzes are executed for some benchmark problems, namely flow in a square cavity, free surface sloshing and moving boundary problems in complex geometries. The results show that the method achieves high accuracy and has high flexibility, even for the flows involving high Reynolds number, complicated geometries, moving boundaries and free surfaces. Copyright © 2000 John Wiley & Sons, Ltd.     

Resonance frequencies and stability of a current-carrying suspended nanobeam in a longitudinal magnetic field

The resonance frequencies and stability of a nanobeam in a longitudinal magnetic field are investigated. To this aim, a three dimensional beam model is used in which the small-scale effect is taken into account based on the nonlocal elasticity theory. The Lorentz forces are obtained in terms of the local elastic rotations of the beam and the thermal stress due to current is modeled as an axial compressive force. Using the Galerkin method, the governing equations of motion are solved and the stability boundary of the nanobeam is determined.     

解三维摩擦接触问题的一个二阶锥线性互补法

针对三维摩擦接触问题的求解,给出了一种基于参变量变分原理的二阶锥线性互补法. 首先,基于三维Coulomb摩擦锥在数学表述上属于二阶锥的事实,利用二阶锥规划对偶理论,建立了三维Coulomb摩擦接触条件的参变量二阶锥线性互补模型,它是二维Coulomb摩擦接触条件参变量线性互补模型在三维情形下的自然推广;随后,利用参变量变分原理与有限元方法,建立了求解三维摩擦接触问题的二阶锥线性互补法. 较之于将三维Coulomb摩擦锥进行显式线性化的线性互补法,该方法无需对三维Coulomb摩擦锥进行线性化,因而在保证精度的前提下所解问题的规模要小很多. 最后通过算例展示了该方法的特点.      

Experimental estimation of the load distribution in bearings by the method of caustics

The objective of this paper is the experimental estimation of the load distribution in roller-bearing by the metrhod of caustics. Contact problems have many practical and important applications¹. For the solution of such problems, besides mechanical analysis, the experimental method of caustics² can also be applied. The optical method of caustics is suitable for the experimental study of singularities in stress fields created either by discontinuities or by loading. Previously, caustics has been applied to the study of singular stress fields developed near concentrated or uniformly distributed loads which are applied along straight boundaries.^3,4 In this work, it is applied to study the load distribution in rollerbearings.     

Fourth‐order accurate compact‐difference discretization method for Helmholtz and incompressible Navier–Stokes equations

A fourth‐order accurate solution method for the three‐dimensional Helmholtz equations is described that is based on a compact finite‐difference stencil for the Laplace operator. Similar discretization methods for the Poisson equation have been presented by various researchers for Dirichlet boundary conditions. Here, the complicated issue of imposing Neumann boundary conditions is described in detail. The method is then applied to model Helmholtz problems to verify the accuracy of the discretization method. The implementation of the solution method is also described. The Helmholtz solver is used as the basis for a fourth‐order accurate solver for the incompressible Navier–Stokes equations. Numerical results obtained with this Navier–Stokes solver for the temporal evolution of a three‐dimensional instability in a counter‐rotating vortex pair are discussed. The time‐accurate Navier–Stokes simulations show the resolving properties of the developed discretization method and the correct prediction of the initial growth rate of the three‐dimensional instability in the vortex pair. Copyright © 2004 John Wiley & Sons, Ltd.     

随机杆系结构几何非线性分析的递推求解方法

建立了随机静力作用下考虑几何非线性的随机杆系结构的随机非线性平衡方程. 将和 位移耦合的随机割线弹性模量以及随机响应量表示为非正交多项式展开式，运用传统的摄动方法获 得了关于非正交多项式展式的待定系数的确定性的递推方程. 在求解了待定系数后，利用非 正交多项式展开式和正交多项式展开式的关系矩阵，可以很方便地得到未知响应量的二阶统计矩. 两杆结构和平面桁架拱的算例结果表明，当随机量涨落较大时，递推随机有限元方法比基于 二阶泰勒展开的摄动随机有限元方法更逼近蒙特卡洛模拟结果，显示了该方法对几何非线性 随机问题求解的有效性.