A method for reconstruction of residual stress fields from measurements made in an incompatible region

A method is introduced by which the complete state of residual stress in an elastic body may be inferred from a limited set of experimental measurements. Two techniques for carrying out this reconstruction using finite element analysis are compared and it is shown that for exact reconstruction of the stress field via this method, the stress field must be measured over all eigenstrain-containing regions of the object. The effects of error and incompleteness in the measured part of the stress field on the subsequent analysis are investigated in a series of numerical experiments using synthetic measurement data based on the NeT TG1 round-robin weld specimen. It is hence shown that accurate residual stress field reconstruction is possible using measurement data of a quality achievable using current experimental techniques.     

3D dynamic stress intensity factors around two parallel square cracks in an infinite elastic medium under impact load

Summary  Transient stresses around two parallel cracks in an infinite elastic medium are investigated in the present paper. The shape of the cracks is assumed to be square. Incoming shock stress waves impinge upon the two cracks normal to tzheir surfaces. The mixed boundary value equations with respect to stresses and displacements are reduced to two sets of dual integral equations in the Laplace transform domain using the Fourier transform technique. These equations are solved by expanding the differences in the crack surface displacements in a double series of a function that is equal to zero outside the cracks. Unknown coefficients in the series are calculated using the Schmidt method. Stress intensity factors defined in the Laplace transform domain are inverted numerically to the physical space. Numerical calculations are carried out for transient dynamic stress intensity factors under the assumption that the shape of the upper crack is identical to that of the lower crack. Received 2 February 2000; accepted for publication 10 May 2000     

3D dynamic stress intensity factors at three rectangular cracks in an infinite elastic medium subjected to a time-harmonic stress wave

Summary Dynamic stresses around three coplanar cracks in an infinite elastic medium are determined in the paper. Two of the cracks are equal, rectangular and symmetrically situated on either side of the centrally located rectangular crack. Time-harmonic normal traction acts on each surface of the three cracks. To solve the problem, two kind of solutions are superposed: one is a solution for a rectangular crack in an infinite elastic medium, and the other one is that for two rectangular cracks in an infinite elastic medium. The unknown coefficients in the combined solution are determined by applying the boundary conditions at the surfaces of the cracks. Finally, stress intensity factors are calculated numerically for several crack configurations. Received 14 July 1998; accepted for publication 2 December 1998     

Elastic field in an infinite medium of one-dimensional hexagonal quasicrystal with a planar crack

This present work is concerned with planar cracks embedded in an infinite space of one-dimensional hexagonal quasicrystals. The potential theory method together with the general solutions is used to develop the framework of solving the crack problems in question. The mode I problems of three common planar cracks (a penny-shaped crack, an external circular crack and a half-infinite crack) are solved in a systematic manner. The phonon and phason elastic fundamental fields along with some important parameters in crack analysis are explicitly presented in terms of elementary functions. Several examples are given to show the applications of the present fundamental solutions. The validity of the present solutions is discussed both analytically and numerically. The derived analytical solutions of crack will not only play an important role in understanding the phonon–phason coupling behavior in quasicrystals, but also serve as benchmarks for future numerical studies and simplified analyses.     

Non-axisymmetric thermal stress of a functionally graded coated circular inclusion in an infinite matrix

This paper is to study the non-axisymmetric two-dimensional problem of thermal stresses in an infinite matrix with a functionally graded coated circular inclusion based on complex variable method. With using the method of piece-wise homogeneous layers, the general solution for the functionally graded coating having radial arbitrary elastic properties is derived when the matrix is subjected to uniform heat flux at infinity, and then numerical results are presented for several special examples. It is found that the existence of the functionally graded coating can change interfacial thermal stresses, and choosing proper change ways of the radial elastic properties in the coating can obviously reduce the thermal stresses.     

Elastic and elastic–plastic singular fields around a crack-tip in particulate-reinforced composites with progressive debonding damage

This paper deals with elastic and elastic–plastic singular fields around a crack-tip in particulate-reinforced composites with debonding damage of particle-matrix interface. Numerical analyses are carried out on a crack-tip field in elastic-matrix and elastic–plastic-matrix composites reinforced with elastic particles, using a finite element method developed based on an incremental damage theory of particulate-reinforced composites. A particle volume fraction and interfacial strength between particles and matrix of the composites are parametrically changed. In the elastic-matrix composites, a unique elastic singular field is created on the complete damage zone in the vicinity of a crack-tip in addition to the conventional elastic singular field on the no damage zone. The macroscopic stress level around a crack-tip is reduced by the debonding damage while the microscopic stress level of the matrix remains unchanged. In the elastic–plastic-matrix composites, the damage zone develops in addition to the plastic zone due to matrix plasticity, and both the macroscopic and microscopic stress revels around a crack-tip are reduced by the debonding damage. It is concluded from the numerical results that the toughening due to damage could be expected in the elastic–plastic-matrix composites, while it is questionable in the elastic-matrix composites.     

Lagrangian velocity and acceleration statistics of fluid and inertial particles measured in pipe flow with 3D particle tracking velocimetry

Three-dimensional particle tracking velocimetry (3D-PTV) has been applied to particle-laden pipe flow at Reynolds number 10,300, based on the bulk velocity and the pipe diameter. The volume fraction of the inertial particles was equal to 1.4 × 10⁻⁵. Lagrangian velocity and acceleration statistics were determined both for tracers and for inertial particles with Stokes number equal to 2.3, based on the particle relaxation time and the viscous time scale. The decay of Lagrangian velocity and acceleration correlation functions was measured both for the fluid and for the dispersed phase at various radial positions. The decay of Lagrangian velocity correlations is faster for inertial particles than for flow tracers, whereas the decay of Lagrangian acceleration correlations is about 25% slower for inertial particles than for flow tracers. Further differences between inertial and tracer particles are found in velocity fluctuations evaluated for both positive and negative time lags. The asymmetry in time of velocity cross-correlations is more pronounced for inertial particles. Quadrant analysis revealed another difference still near the wall: ejection and sweep events are less frequent for inertial particles than for tracers.     

弹性SHPB试件内部纵向应力波传播及所产生的横向约束应力

采用考虑横向惯性效应的Rayleigh-Love杆理论分析了一个弹性试件在分离式霍普金森压杆(SHPB)加载过程中的内部弹性波传播过程,运用Laplace变换和反变换方法,得到了试件内部各点的变形、速度、应变和应力解析解.通过数值计算,得到梯形入射波加载情况下,纵向应力在试件内部的连续变化过程,以及波传播所伴随的横向附加应力.计算表明：在试件/入射杆界面附近,初次加载所产生的横向附加应力最大,可达入射波平台的12%;在大部分试件区域,纵向应力波传播将造成入射波平台4%～6%的横向附加应力;材料的泊松比越大,或者杆/试件声阻抗比越小,所伴随的横向附加应力越大;梯形波的上升时间和试件长径比对横向附加应力影响不大.     

2-3型压电复合材料有效特性及其静水压性能的研究

压电复合材料是由至少一种压电材料与非压电材料按照一定的连通方式组合在一起而构成的具有压电效应的材料。采用均匀场理论,提出了2-3型压电复合材料的代表性体积单元,通过建立串联、并联等模型预估了压电复合材料的弹性刚度常数:C11、C12、C13、C22、C23、C33,压电常数:e31e、32e、33,介电常数:3ε3。然后用有限元方法(使用了有限元分析软件ANSYS)分析了压电复合材料的材料常数,并用有限元方法验证了均匀场理论求解2-3型压电复合材料有效特性参数的准确性。根据求得的有效特性参数进一步研究了2-3型压电复合材料的静水压性能,由于求得的静水压性能值很低,又改良了2-3型压电复合材料的模型,获得了较好的静水压性能。并将改良前后2-3型压电复合材料与3-2型压电复合材料的静水压性能进行了比较。     

On Cauchy singular integrals and stress intensity factors for 2D mode I cracks in an infinite solid

This paper comments on some of the different numerical techniques commonly employed in evaluating Cauchy singular integrals of the first kind; e.g. as pertaining to 2D through cracks in a brittle material undergoing Mode I loading. In addition, a different more direct method is proposed here. Also, two different ways to calculate the stress intensity factor (K_I) are contrasted. The accuracy attained by the different methods in calculating K_I, and the factors affecting the calculation, are compared. Finally, comments on calculating the stress field of a 2D crack and important considerations are presented.     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

The symmetry and loading-independency of multiple inclusions enclosing uniform stresses in an infinite elastic plane

The identification of multiple interacting inclusions with uniform internal stresses in an infinite elastic matrix subjected to a uniform remote loading is of fundamental importance in the mechanics and design of particulate composite materials. In anti-plane shear and plane deformations, certain sufficient conditions have been established in the literature which guarantee uniform internal stresses inside multiple interacting inclusions displaying various symmetries when the matrix is subjected ...     

A rigid inclusion in an elastic space under the action of a uniform heat flow in the inclusion plane

A solution is presented for the three dimensional static thermoelastic problem of an absolutely rigid inclusion (anticrack) in the case when a uniform heat flow is directed along the inclusion plane. By using the potential method and the Fourier transform technique, the problem is reduced to a system of coupled two-dimensional singular integral equations for the shear stress jumps across the inclusion. As an illustration, a typical application to the circular anticrack is presented. Explicit expressions for the thermal stresses in the inclusion plane are obtained and discussed from the point of view of material failure.     

二维弹性接触问题中随机边界元法与可靠性分析

发展了二维弹性接触问题中的随机边界元法，推导并建立了相应的随机边界元基本方程，并将所发展的方法用于静强度的可靠性分析，讨论了其数值解技术。通过算例分析表明，本文发展的方法是可行的。     

A shock tube study on the thermal decomposition of CS2 based on S(3P) and S(1D) concentration measurements

The thermal decomposition of CS₂ highly diluted in Ar was studied behind reflected shock waves by monitoring time-dependent absorption profiles of S(³P) and S(¹D) using atomic resonance absorption spectroscopy (ARAS). The rate coefficient of the reaction:

Influence of local cubic anisotropy on the transition towards an equipartition regime in a 3D texture-less random elastic medium

At long lapse times in randomly fluctuating media with macroscopic isotropy (texture-less media), the energy of elastic waves is equipartitioned between compressional (P) and shear (S) waves. This property is independent of the local isotropy or anisotropy of the heterogeneous constitutive tensor and of the type of source. However the local symmetry of the constitutive tensor does influence the rate of convergence to equipartition and this paper discusses the precise influence of local anisotropy on the time required to reach equipartition. More particularly, a randomly-fluctuating medium is considered, whose behavior is statistically isotropic, and locally cubic. After calculating all the differential and total scattering cross-sections in that case, an analytical formula is derived for the rate of convergence to the equipartition regime, function of the second-order statistics of the mechanical parameter fields (bulk and shear moduli and anisotropy parameter). The local anisotropy is shown to influence strongly that transition rate, with a faster transition when the fluctuations of the anisotropy parameter are positively correlated to those of the shear modulus. A numerical model is constructed to illustrate numerically these results. Since the asymptotic regime of equipartition cannot be simulated directly because it would require too large a computational domain, boundaries are introduced and mechanical properties are chosen so as to minimize their influence on equipartition.     

拉压不同弹性模量薄板上圆孔的应力分析

采用弹性理论研究了拉压不同弹性模量薄板上圆孔的孔边应力集中问题．采用广义虎克定律推导出了拉压不同弹性模量薄板上圆孔边的应力平衡方程,并联合利用应力函数及边界条件得到了拉压不同弹性模量薄板上圆孔边的应力表达式．算例分析表明,当薄板材料的拉压弹性模量相差较大时,采用经典弹性理论研究薄板上圆孔的孔边应力是不合适的,当经典弹性理论与拉压不同弹性模量弹性理论的计算结果间的差别超过工程允许误差5%时,应该采用拉压不同弹性模量弹性理论进行计算．     

Numerical simulation of cavitating flow in 2D and 3D inducer geometries

A computational method is proposed to simulate 3D unsteady cavitating flows in spatial turbopump inducers. It is based on the code FineTurbo, adapted to take into account two‐phase flow phenomena. The initial model is a time‐marching algorithm devoted to compressible flow, associated with a low‐speed preconditioner to treat low Mach number flows. The presented work covers the 3D implementation of a physical model developed in LEGI for several years to simulate 2D unsteady cavitating flows. It is based on a barotropic state law that relates the fluid density to the pressure variations. A modification of the preconditioner is proposed to treat efficiently as well highly compressible two‐phase flow areas as weakly compressible single‐phase flow conditions. The numerical model is applied to time‐accurate simulations of cavitating flow in spatial turbopump inducers. The first geometry is a 2D Venturi type section designed to simulate an inducer blade suction side. Results obtained with this simple test case, including the study of its general cavitating behaviour, numerical tests, and precise comparisons with previous experimental measurements inside the cavity, lead to a satisfactory validation of the model. A complete three‐dimensional rotating inducer geometry is then considered, and its quasi‐static behaviour in cavitating conditions is investigated. Numerical results are compared to experimental measurements and visualizations, and a promising agreement is obtained. Copyright © 2004 John Wiley & Sons, Ltd.     

Symmetrical Representation of Stresses in the Stroh Formalism and its Application to a Dislocation and a Dislocation Dipole in an Anissotropic Elastic Medium

Symmetrical stress representation in the Stroh formalism for anisotropic elastic bodies is introduced and the range of its applicability is analysed. By making use of this stress representation new formulae for influence functions giving stresses in an infinite anisotropic medium subjected to a straight dislocation and a straight dislocation dipole are derived. The advantage of the new formulae is that they explicitly show the symmetrical structure of these influence functions not referred to previously. Relations of these influence functions to influence functions giving stresses and Airy stress function due to a straight wedge disclination, whose explicit expressions are also introduced, are derived. Application of these results in computation of stresses by the hypersingular and regularized Somigliana stress identities is discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

共和干热岩地热井GR2井区三维地应力场反演分析

GR2 地热井位于共和恰卜恰镇，干热岩地热井井深达到4000 米，大深度井区极大增加了反演的难度．井区附近20 千米内无大型构造，附近地应力场的主要影响因素为板块构造应力和自重应力．由于井中未进行三维地应力场测井工作，所以我们以地面地应力调查结果和震源机制解反演结果为依据，基于MIDAS GTS 和大比例尺高清卫星云图建立有限元数值模型，使用改进的边界法进行反演计算，最终预估了井中地应力场分布特征．对于后续的井网布置以及干热岩井场开采具有重要意义．