The contact problem of a rubber half-space dented by a rigid cone apex

Summary  The smooth contact of a rubber half-space dented by a rigid cone apex is analyzed based on the large deformation theory. The problem is treated as an axially symmetric case, and the material is assumed to be incompressible. The asymptotic equations for the domain near the apex are derived. They are solved analytically for the shrinking domain, while a numerical solution is found for the expanding domain in the vicinity of the stress singularity. The purpose of this paper is not only to solve a typical problem but also to provide an analytical method to solve a large-strain problem with a singular point. Received 10 July 2001; accepted for publication 24 January 2002     

On generalized Kelvin solutions in a multilayered elastic medium

This paper presents fundamental singular solutions for the generalized Kelvin problems of a multilayered elastic medium of infinite extent subjected to concentrated body force vectors. Classical integral transforms and a backward transfer matrix method are utilized in the analytical formulation of solutions in both Cartesian and cylindrical coordinates. The solution in the transform domain has no functions of exponential growth and is invariant with respect to the applied forces. The convergence of the solutions in the physical domain is rigorously and analytically verified. The solutions satisfy all required constraints including the basic equations and the interfacial conditions as well as the boundary conditions. In particular, singular terms of the generalized Kelvin solutions associated with the point and ring types of concentrated body force vectors are obtained in exact closed-forms via an asymptotic analysis. Numerical results presented in the paper illustrate that numerical evaluation of the solutions can be easily achieved with very high accuracy and efficiency and that the layering material inhomogeneity has a significant effect on the elastic field.The Canadian Government right to retain a non-exclusive, royalty-free licence in and to any copyright is acknowledged.     

Large deformation field near a crack tip in rubber-like material

This paper deals with strain field near a crack tip in a rubber-like material under plane strain condition. The constitutive relation adopted here is valid for both small and large strain. The asymptotic equations are derived for a shrinking sector and expanding sector. The closed mathematic solution is obtained for the latter while a numerical solution is found for the former. By connecting deformation of the two sectors, the crack tip field character is found.     

Thermodynamic driving force in ferroelectric crystals with a rank-2 laminated domain pattern, and a study of enhanced electrostriction

The thermodynamic driving force for domain growth in a rank-2 laminated ferroelectric crystal is derived in this article, and we used it, together with a homogenization theory, to study the issue of enhanced electrostrictive actuation recently reported by Burcsu et al. [2004. Large electrostrictive actuation of barium titanate single crystals. J. Mech. Phys. Solids 52, 823-846]. We derived this force from the reduction of Gibbs free energy with respect to the increase of domain concentration. It is shown that both the free energy and the thermodynamic force consist of three parts: the first arises from the difference in M₀ and M₁, the linear electromechanical compliances of the parent and product domains, respectively, at a given level of applied stress and electric field, the second stems from the electromechanical work associated with the change of spontaneous strain and spontaneous polarization during domain switch, and the third from the internal energy due to the distribution of polarizations strain and electric polarization inside the crystal. We prove that the first term is substantially lower than the second one, and the third one is identically zero with compatible domain pattern. The second one is, however, not exactly equal to the commonly written sum of the products of stress with strain, and electric field with polarization during switch, unless both domains have identical moduli in the common global axes. We also show that, with compatible domain patterns and when M₁=M₀, this driving force is identical to Eshelby's driving force acting on a flat interface due to the jump of energy-momentum tensor. Applications of the theory to a BaTiO₃ crystal subjected to a fixed axial compression and decreasing electric field from the [0 0 1] state reveal that the crystal undergoes a three-stage switching process: (i) the 0→90° switch to form a rank-1 laminate, (ii) the 0→180° switch inside the 0° domain to form a laminate I with a concurrent 90°→−90° switch inside the 90° domain to form laminate II, creating a rank-2-laminated domain pattern, and (iii) finally the 90→180° switch. It is the exchange of stability between the 0, 90°, and 180° domains under compression and electric field that is the origin of the enhanced actuation. We illustrate these intrinsic features by showing the evolution of these domains, and demonstrate how the reported large actuation strain can be attained with a rank-2 laminate.     

Large strain field near a crack tip in a rubber sheet

The distribution of stress-strain near a crack tip in a rubber sheet is investigated by employing the constitutive relation given by Gao (1997). It is shown that the crack tip field is composed of two shrinking sectors and one expanding sector. The stress state near the crack tip is in uniaxial tension. The analytical solutions are obtained for both expanding and shrinking sectors.     

A transverse crack tip field in bimaterial

The asymptotic field near an interface crack tip is analyzed with the fully nonlinear theory. By dividing the crack tip field into narrowing sectors and an expanding sector, the asymptotic equations for the crack tip field are derived and solved. The singular characters of stress and strain near the crack tip are revealed.     

Analysis of the Interface Crack for Rubber-like Materials

The interface crack between two dissimilar rubber materials under mixed load is asymptotically analyzed in this paper. It is shown that the crack tip field is composed of one expanding sector and two shrinking sectors. For the case considered, the interface is located in the expanding sector. Analytical solutions are obtained for both the expanding sector and the shrinking sectors. The structures of the crack tip field obtained in this paper is compared with those in previous works for different constitutive equations.     

Spheroidal wave solutions for sound radiation problems in the near-field of planar structures

Current research in active noise control and in the reconstruction of vibrating sources often requires knowledge of the independent source–field components that best represent the complex acoustical transfer paths observed between a radiating structure and a given control or observation domain. In this paper, closed-form solutions are provided for the singular value expansion of the radiation operator that maps the boundary velocity of a baffled rectangular structure onto the acoustic pressure observed in the half-space domain over a hemi-spheroidal surface located at an arbitrary separation distance from the radiator, including in the near-field zone. Independent contributions of the evanescent and propagating wave components to the complex power are examined for a baffled beam when varying the frequency and the source–field distance parameter. It is shown that the reactive-to-active power ratio induced by each singular mode follows an inverse power law that scales on the product between the reduced frequency and the source–field distance parameter. A transitional region is defined in the space-frequency domain within which the reactive power components are preponderant and should be accounted for when controlling or imaging the near-field zone of a planar radiator. The optimality of the singular source modes is found to be of interest to actively reduce the active and reactive power components in the near-field zone of a radiator with a limited number of independent control channels.     

Helmholtz边界积分方程中奇异积分间接求解方法

提出了间接求解传统Helmholtz边界积分方程CBIE的强奇异积分和自由项系数,以及Burton-Miller边界积分方程BMBIE中的超强奇异积分的特解法。对于声场的内域问题,给出了满足Helmholtz控制方程的特解,间接求出了CBIE中的强奇异积分和自由项系数。对于声场外域对应的BMBIE中的超强奇异积分,按Guiggiani方法计算其柯西主值积分需要进行泰勒级数展开的高阶近似,公式繁复,实施困难。本文给出了满足Helmholtz控制方程和Sommerfeld散射条件的特解,提出了间接求出超强奇异积分的方法。推导了轴对称结构外场问题的强奇异积分中的柯西主值积分表达式,并通过轴对称问题算例证明了本文方法的高效性。数值结果表明,对于内域问题,采用本文特解法的计算结果优于直接求解强奇异积分和自由项系数的结果,且本文的特解法可避免针对具体几何信息计算自由项系数,因而具有更好的适用性。对于外域问题,两者精度相当,但本文的特解法可避免对核函数进行高阶泰勒级数展开,更易于数值实施。     

On stability of elastic domain during isothermal solid-solid phase transformation in a tube configuration

Under isothermal quasi-static stretching the phasetransition of a superelastic NiTi tube involves the formation(during loading) and vanishing(in unloading) of a high strain(martensite) domain.The two events are accompanied by arapid stress drop/rise due to the formation/vanishing of domain fronts.From a thermodynamic point of view,both areinstability phenomena that occur once the system reaches itscritical state.This paper investigates the stability of a shrinking cylindrical domain in a tube configuration during unloading.The energetics and thermodynamic driving force of thecylindrical domain are quantified by using an elastic inclusion model.It is demonstrated that the two domain fronts exhibit strong interaction when they come close to each other,which brings a peak in the total energy and a sign changein the thermodynamic driving force.It is proved that suchdomain front interaction plays an important role in controlling the stability of the domain and in the occurrence of stressjumps during domain vanishing.It is also shown that the process is governed by two nondimensional length scales(thenormalized tube length and normalized wall-thickness) andthat the length scale dependence of the critical domain lengthand stress jump for the domain vanishing can be quantifiedby the elastic inclusion model.     

集中力作用下的两相饱和介质位移场Green函数

以复模两相饱和介质Biot动力学方程为基础,根据该方程D'Alembert解的Fourier变换所属的Homholtz方程特性,由Biot方程解的相容性条件及δ函数性质较好地解决了快、慢纵波位势的耦合问题.较为简便地得到了两相饱和介质在集中力作用下低频(ω＜ωc)时的频域和时域的Green函数.     

非牛顿流体在渐变管中压力和剪切应力的二次摄动解

本文利用双摄动方法求解缓慢变化管道中Johnson-Segalman(J-S)流体流动的渐近解．将管道的扩张（或收缩）角度和粘弹性参数分别作为双摄动的参数,由流函数和涡量函数的形式,推导出压力和壁面剪切应力的渐近解．在此基础上,分析了管道角度,粘弹性参数和雷诺数等参数对压力以及剪切力影响．主要结论如下：(1) 管道扩张角度增加时,流向同一位置处径向压力以及壁面剪切应力随扩张角度减小;(2) 在同一扩张管道中,径向压力随着流向位移减小,收缩管与之相反;(3) 扩张角度与雷诺数对流场起主导作用,粘弹性系数起次要作用．     

Prediction of thermal crack spacing

When the thermally induced stress in a shrinking pavement layer reaches the tensile strength of asphalt, regularly spaced thermal cracks form across the width of the pavement. A one-dimensional analytical solution for the stress distribution in a thermally shrinking elastic pavement layer placed on an elastoplastic, cohesive–frictional base is developed and validated by comparison with a 2D numerical solution. From the analytical model, a prediction of a length parameter that provides bounds on the thermal crack spacing is obtained. Predicted bounds on crack spacing are validated by comparison with field observations. It is demonstrated that the proposed formulation can also be applied to estimate the average crack density observed in thin ceramic films subjected to the application of an axial strain; in the latter system, the crack spacing is six decades smaller than that observed in pavement systems.     

膨胀环受力历史对应力应变关系的影响

利用两类实验装置开展了无氧铜TU1膨胀环实验研究,发现:电磁膨胀环在加载阶段,样品受体力作用,满足均匀变形的假定;而爆炸膨胀环在加载阶段,样品内壁受面力冲击作用,不满足均匀变形的假定。针对这个差异,发展了一种考虑冲击阶段变形不均匀性的新方法,利用回收样品几何变形,将冲击阶段试样环内轴向塑性应变、径向塑性应变纳入等效塑性应变的计算中,通过修正后的方法更准确地获得了材料的应力应变关系。     

A Normal Force on the Free Surface of a Coated Material

The three-dimensional theoretical solution of a concentrated normal force acting on the free surface of a coated material has been deduced by applying the reflection method. It is found that all stress functions defined in the local coordinate systems with their origins placed at each mirror point can be deduced from the fundamental solution of a concentrated normal force acting on the free surface of a semi-infinite homogeneous medium. The structure of the elastic solution has been illustrated by numerical analysis. It is found that only the stress functions corresponding to the first few mirror points are influential. It is also found that the effect of material combination on the stress field shall be described by three parameters, the two Dundurs' parameters and one additional parameter.     

The self-similarly expanding Eshelby ellipsoidal inclusion: I. Field solution

The solution of a self-similarly (subsonically) dynamically expanding ellipsoidal inclusion with general spatially uniform transformation strain temporally constant is obtained by the use of the Radon transform and the satisfaction of the zero initial conditions and the radiation condition at infinity. It constitutes the self-similar evolution of the inclusion singularity (jump discontinuity at the inclusion boundary) starting from zero dimension. The field solutions for the displacement gradient and particle velocity are presented. Due to the fact that for a self-similarly expanding subsonic motion the hyperbolic system of the partial differential equations of motion becomes elliptic (as proved in Ni and Markenscoff, 2015), it is shown here explicitly that the solution for the displacement gradient in the interior domain of the expanding ellipsoid is constant, thus extending the Eshelby property to the self-similarly expanding ellipsoids as pointed out by Burridge and Willis (1969). Also, the particle velocity is shown to be zero in the interior domain (lacuna) as the waves emitted by the self-similarly expanding inclusion cancel each other due to the symmetries of geometry and motion.     

A real-space phase field model for the domain evolution of ferromagnetic materials

A real-space phase field model based on the time-dependent Ginzburg–Landau (TDGL) equation is developed to predict the domain evolution of ferromagnetic materials. The phase field model stems from a thermodynamic theory of ferromagnetic materials which employs the strain and magnetization as independent variables. The phase field equations are shown to reduce to the common micromagnetic model when the magnetostriction is absent and the magnitude of magnetization is constant. The strain and magnetization in the equilibrium state are obtained simultaneously by solving the phase field equations via a nonlinear finite element method. The finite-element based phase field model is applicable for the domain evolution of ferromagnetic materials with arbitrary geometries and boundary conditions. The evolution of magnetization domains in ferromagnetic thin film subjected to external stresses and magnetic fields are simulated and the magnetoelastic coupling behavior is investigated. Phase field simulations show that the magnetization vectors form a single magnetic vortex in ferromagnetic disks and rings. The configuration and size of the simulated magnetization vortex are in agreement with the experimental observation, suggesting that the phase field model is a powerful tool for the domain evolution of ferromagnetic materials.     

Force at a point in the interior of a layered elastic half space

This study formulates, by the technique of integral transforms, the solution of a layered half space subjected to a concentrated force which may act either vertically or horizontally in the interior of the system. Accurate approximations of the reciprocals of the common denominators in the solution integrals are suggested in such a way that the latter are in standard closed forms and can be identified by two parts. The first part is the singular part of Mindlin's solution which is singular at the point of application of the force, and the second is non-singular. The solutions for plane problems are also obtained in closed forms by performing appropriate integrations of the solutions for the corresponding three-dimensional cases.     

Electroelastic interaction between a piezoelectric screw dislocation and a confocal elliptic blunt crack with elliptical inhomogeneity

The electroelastic interaction between a piezoelectric screw dislocation and an elliptical inhomogeneity containing a confocal blunt crack under infinite longitudinal shear and in-plane electric field is investigated. Using the sectionally holomorphic function theory, Cauchy singular integral, singularity analysis of complex functions and theory of Rieman boundary problem, the explicit series solution of stress field is obtained when the screw dislocation is located in inhomogeneity. The intervention law of the interaction between blunt crack and screw dislocation in inhomogeneity is discussed. The analytical expressions of generalized stress and strain field of inhomogeneity are calculated, while the image force, field intensity factors of blunt crack are also presented. Moreover, a new matrix expression of the energy release rate and generalized strain energy density (SED) are deduced. With the size variation of blunt crack, the results can be reduced to the case of the interaction between a piezoelectric screw dislocation and a line crack in inhomogeneity. Numerical analysis are then conducted to reveal the effects of the dislocation location, the size of inhomogeneity and blunt crack and the applied load on the image force, energy release rate and strain energy density. The influence of dislocation on energy release rate and strain energy density is also revealed.     

Poincaré’s Variational Problem in Potential Theory

One of the earliest attempts to rigorously prove the solvability of Dirichlet’s boundary value problem was based on seeking the solution in the form of a “potential of double layer”, and this leads to an integral equation whose kernel is (in general) both singular and non-symmetric. C. Neumann succeeded with this approach for smoothly bounded convex domains, and H. Poincaré, by a tremendous tour de force, showed how to push through the analysis for domains with sufficiently smooth boundaries but no hypothesis of convexity. In this work he was (according to his own account) guided by consideration of a variational problem involving the partition of energy of an electrostatic field induced by charges placed on the boundary of a domain, more precisely the charge distributions which render stationary the energy of the field inside the domain divided by the energy of the field outside the domain. Unfortunately, a rigorous treatment of this problem was not possible with the tools available at that time (as Poincaré was well aware). So far as we know, the only one to propose a rigorous treatment of Poincaré’s problem was T. Carleman (in the two-dimensional case) in his doctoral dissertation. Thanks to later developments (especially concerning Sobolev spaces, and spectral theory of operators on Hilbert space) we can now give a complete, general and rigorous account of Poincaré’s variational problem, and that is the main object of the present paper. As a by-product, we refine some technical aspects in the theory of symmetrizable operators and prove in any number of dimensions the basic properties of the analogue of the planar Bergman–Schiffer singular integral equation. We interpret Poincaré’s variational principle as a non-selfadjoint eigenvalue problem for the angle operator between two distinct pairs of subspaces of potentials. We also prove a series of novel spectral analysis facts (some of them conjectured by Poincaré) related to the Poincaré–Neumann integral operator.