Analytical solution for motion of an elliptical inclusion in gradient stress field

We present a rigorous analytical solution for motion of an elliptical inclusion in isotropic matrix driven by gradient stress field. The interfacial diffusion is considered as the dominant mechanism for the motion. We demonstrate that normal stress gradient on the interface is the major driven force, while the strain energy density gradient is negligible. A key prediction of the solution is that for a given inclusion the motion velocity is proportional to stress gradient only, indicating that the solution is applicable for inclusion motion in nonuniform stress field of varying stress gradient, and that the inclusion tends to move towards the region of lower stress in nonuniform stressed materials.     

应用边界积分法求圆形夹杂问题的解析解

边界元方法作为一种数值方法, 在各种科学工程问题中得到了广泛的应用.本文参考了边界元法的求解思路, 从Somigliana等式出发, 利用格林函数性质,得到了一种边界积分法, 使之可以用来寻求弹性问题的解析解.此边界积分法也可以从Betti互易定理得到. 应用此新方法, 求解了圆形夹杂问题.首先设定夹杂与基体之间完美连接, 将界面处的位移与应力按照傅里叶级数展开,根据问题的对称性与三角函数的正交性来简化假设, 减少待定系数的个数.其次选择合适的试函数(试函数满足位移单值条件以及无体力的线弹性力学问题的控制方程),应用边界积分法, 求得界面处的位移与应力的值. 然后再求解域内位移与应力.得到了问题的精确解析解, 当夹杂弹性模量为零或趋向于无穷大时,退化为圆孔或刚性夹杂问题的解析解. 求解过程表明,若问题的求解区域包含无穷远处时, 所取的试函数应满足无穷远处的边界条件.若求解区域包含坐标原点, 试函数在原点处位移与应力应是有限的.结果表明了此方法的有效性.      

On a bound for deformations in terms of Strains in an elastic cylinder

Consider a long thin isotropic elastic cylinder with a self-equilibrated loading on each end face, but which is stress-free on the sides and which has no internal body forces. It is shown that if the displacement gradient is pointwise sufficiently small, then, in any subcylinder of length 1/4a, it is possible to add a rigid body motion such that the L ₂ norm of the resulting displacement gradient can be bounded by a constant times the L ₂ norm of the strain in a subcyclinder of length 2a centered at the same point. The parameter a depends upon the pointwise bound for the displacement gradient (the smaller the bound, the larger a can be) and the constant is independent of the length thickness ratio of the subcylinder.     

The energy-release rate and “self-force” of dynamically expanding spherical and plane inclusion boundaries with dilatational eigenstrain

In the context of the linear theory of elasticity with eigenstrains, the radiated field including inertia effects of a spherical inclusion with dilatational eigenstrain radially expanding is obtained on the basis of the dynamic Green's function, and one of the half-space inclusion boundary (with dilatational eigenstrain) moving from rest in general subsonic motion is obtained by a limiting process from the spherically expanding inclusion as the radius tends to infinity while the eigenstrain remains constrained, and this is the minimum energy solution. The global energy-release rate required to move the plane inclusion boundary and to create an incremental region of eigenstrain is defined analogously to the one for moving cracks and dislocations and represents the mechanical rate of work needed to be provide for the expansion of the inclusion. The calculated value, which is the “self-force” of the expanding inclusion, has a static component plus a dynamic one depending only on the current value of the velocity, while in the case of the spherical boundary, there is an additional contribution accounting for the jump in the strain at the farthest part at the back of the inclusion having the time to reach the front boundary, thus making the dynamic “self-force” history dependent.     

On unsteady unidirectional flows of a second grade fluid

Some properties of unsteady unidirectional flows of a fluid of second grade are considered for flows produced by the sudden application of a constant pressure gradient or by the impulsive motion of one or two boundaries. Exact analytical solutions for these flows are obtained and the results are compared with those of a Newtonian fluid. It is found that the stress at the initial time on the stationary boundary for flows generated by the impulsive motion of a boundary is infinite for a Newtonian fluid and is finite for a second grade fluid. Furthermore, it is shown that initially the stress on the stationary boundary, for flows started from rest by sudden application of a constant pressure gradient is zero for a Newtonian fluid and is not zero for a fluid of second grade. The required time to attain the asymptotic value of a second grade fluid is longer than that for a Newtonian fluid. It should be mentioned that the expressions for the flow properties, such as velocity, obtained by the Laplace transform method are exactly the same as the ones obtained for the Couette and Poiseuille flows and those which are constructed by the Fourier method. The solution of the governing equation for flows such as the flow over a plane wall and the Couette flow is in a series form which is slowly convergent for small values of time. To overcome the difficulty in the calculation of the value of the velocity for small values of time, a practical method is given. The other property of unsteady flows of a second grade fluid is that the no-slip boundary condition is sufficient for unsteady flows, but it is not sufficient for steady flows so that an additional condition is needed. In order to discuss the properties of unsteady unidirectional flows of a second grade fluid, some illustrative examples are given.     

含椭圆形刚性夹杂的压电材料平面问题

应用复变函数的Ｆａｂｅｒ级数展开方法，本文研究了含椭圆形刚性夹杂的压电材料平面问题，给出了问题的封闭解。解签表明，夹杂内的电场强度和电位移为常量。并通过算例分析，讨论了正，逆压电效应在基体孔周处的机电行为。     

Joint unidirectional motion of two viscous heat-conducting fluids in a tube

This paper studies an invariant solution of the problem of joint motion of two heat-conducting viscous immiscible fluids which have a common interface in a cylindrical tube under an unsteady pressure gradient. The problem reduces to a coupled initial-boundary-value problem for parabolic equations. A priori estimates of velocity and temperature perturbations are obtained. The steady state of the system is determined, and it is proved that if, in one of the fluids, the pressure gradient rapidly approaches zero, the perturbations of all quantities tend to zero. It is shown that if the pressure gradient has a nonzero limit, the solution reaches a steady state. In this case, the velocity field in the limit is the same as in conjugate Poiseuille flow, and the temperature is represented as a polynomial of the fourth order on the radial coordinate.     

Thermo-fluid dynamics of the unsteady channel flow

We present an exact analytical representation of the unsteady thermo-fluid dynamic field arising in a two-dimensional channel with parallel walls for a fluid with constant properties. We assume that the axial pressure gradient is an arbitrary function of time that can be expanded in Taylor series; a particular case is the impulsive motion generated by a sudden jump to a constant value; for large time values the flow reaches the well-known steady Poiseuille solution. As boundary conditions for the dynamic field we consider fixed and moving walls (unsteady Couette flow). The assigned temperature on the walls can be an arbitrary function of time. We also consider the coupling of the energy and momentum equations (i.e. Eckert number different from zero). The solution is obtained by series with simple expressions of the coefficients in terms of the error functions. The fundamental physical parameters, such as shear stress, mass flow and heat flux at the wall are obtained in explicit analytical form and discussed by means of their diagrams.     

Plane problems in piezoelectric media with elliptic inclusion

I.IntroductionPiezoelectricmedia,asa"ex\'typeoffullctionalmaterial.arex'idel}'appliedtomanytechnologicalfieldsduetoitselectronlechallicalcouplillgeffect.Defects.likethatofothermaterials.arenotlimitedtocracks.x'oidsandinclusionsillpiezoelectricmaterialsorelements.Yet,stressconcentrationsornoll-ullitbrllldistl-ibutionsofelectricfieldillducedbythosedefectsareoneofthehe}l'filctorswllicllwouldleadpiezoelectricstructurestonon-normalfailure.Therel'ore.itisofgrealimportancetostudythepropertiesofthos…     

Non-linear stochastic response of a shallow cable

The paper considers the stochastic response of geometrical non-linear shallow cables. Large rain-wind induced cable oscillations with non-linear interactions have been observed in many large cable stayed bridges during the last decades. The response of the cable is investigated for a reduced two-degrees-of-freedom system with one modal coordinate for the in-plane displacement and one for the out-of-plane displacement. At first harmonic varying chord elongation at excitation frequencies close to the corresponding eigenfrequencies of the cable is considered in order to identify stable modes of vibration. Depending on the initial conditions the system may enter one of two states of vibration in the static equilibrium plane with the out-of-plane displacement equal to zero, or a whirling state with the out-of-plane displacement different from zero. Possible solutions are found both analytically and numerically. Next, the chord elongation is modelled as a narrow-banded Gaussian stochastic process, and it is shown that all the indicated harmonic solutions now become instable with probability one. Instead, the cable jumps randomly back and forth between the two in-plane and the whirling mode of vibration. A theory for determining the probability of occupying either of these modes at a certain time is derived based on a homogeneous, continuous time three states Markov chain model. It is shown that the transitional probability rates can be determined by first-passage crossing rates of the envelope process of the chord wise component of the support point motion relative to a safe domain determined from the harmonic analysis of the problem.     

SH波对浅埋弹性圆柱及裂纹的散射与地震动

采用Green函数、复变函数和多极坐标等方法研究含圆柱形弹性夹杂的弹性半空间中任意位置、任意方位有限长度裂纹对SH波的散射与地震动. 构造了含圆柱形弹性夹杂的半空间对SH波的散射波,并求解了适合本问题Green函数,即含有圆柱形弹性夹杂的半空间内(表面)任意一点承受时间谐和的出平面线源载荷作用时位移函数的基本解答. 利用裂纹``切割'方法在任意位置构造任意方位的裂纹,可以得到基体中圆柱形弹性夹杂和裂纹同时存在条件下的位移场与应力场. 通过数值算例,讨论各种参数对夹杂上方地表位移的影响.      

Solution of Eshelby's inclusion problem with a bounded domain and Eshelby's tensor for a spherical inclusion in a finite spherical matrix based on a simplified strain gradient elasticity theory

A solution for Eshelby's inclusion problem of a finite homogeneous isotropic elastic body containing an inclusion prescribed with a uniform eigenstrain and a uniform eigenstrain gradient is derived in a general form using a simplified strain gradient elasticity theory (SSGET). An extended Betti's reciprocal theorem and an extended Somigliana's identity based on the SSGET are proposed and utilized to solve the finite-domain inclusion problem. The solution for the disturbed displacement field is expressed in terms of the Green's function for an infinite three-dimensional elastic body in the SSGET. It contains a volume integral term and a surface integral term. The former is the same as that for the infinite-domain inclusion problem based on the SSGET, while the latter represents the boundary effect. The solution reduces to that of the infinite-domain inclusion problem when the boundary effect is not considered. The problem of a spherical inclusion embedded concentrically in a finite spherical elastic body is analytically solved by applying the general solution, with the Eshelby tensor and its volume average obtained in closed forms. This Eshelby tensor depends on the position, inclusion size, matrix size, and material length scale parameter, and, as a result, can capture the inclusion size and boundary effects, unlike existing Eshelby tensors. It reduces to the classical Eshelby tensor for the spherical inclusion in an infinite matrix if both the strain gradient and boundary effects are suppressed. Numerical results quantitatively show that the inclusion size effect can be quite large when the inclusion is very small and that the boundary effect can dominate when the inclusion volume fraction is very high. However, the inclusion size effect is diminishing as the inclusion becomes large enough, and the boundary effect is vanishing as the inclusion volume fraction gets sufficiently low.     

Interface crack extension at any constant speed in orthotropic or transversely isotropic bimaterials––I. General exact solutions

A semi-infinite crack along the interface of two dissimilar half-spaces extends under in-plane loading. Each half-space belongs to a class of orthotropic or transversely isotropic elastic materials, the crack can extend at any constant speed, and all six possible relations between the four body wave speeds are considered. A steady dynamic situation is treated, and exact full displacement fields derived. A key step is a factorization that produces, despite anisotropy, simple solution forms and compact crack speed-dependent functions that exhibit the Rayleigh and Stoneley speeds as roots. These roots are calculated for various representative bimaterials.Closed-form crack opening displacement gradient and interface stress fields are also derived from a general set of coupled singular integral equations. The equation eigenvalues can, depending on crack speed, be complex/imaginary conjugates, purely real, or zero. This suggests possibilities observed in other studies: oscillations and square-root singular behavior at the crack edge, non-singular behavior, singular behavior not of square-root order, and the radiation of displacement gradient discontinuities at crack speeds beyond the purely sub-sonic range.These possibilities are explored further in terms of two important special cases in Part II of this study [Int. J. Solids Struct., 39, 1183–1198].     

Strain gradient solution for a finite-domain Eshelby-type anti-plane strain inclusion problem

A solution for the finite-domain Eshelby-type inclusion problem of a finite elastic body containing an anti-plane strain inclusion of arbitrary cross-sectional shape prescribed with a uniform eigenstrain and a uniform eigenstrain gradient is derived in a general form using a simplified strain gradient elasticity theory (SSGET). The formulation is facilitated by an extended Betti’s reciprocal theorem and an extended Somigliana’s identity based on the SSGET and suitable for anti-plane strain problems. The disturbed displacement field is obtained in terms of the SSGET-based Green’s function for an infinite anti-plane strain elastic body. The solution reduces to that of the infinite-domain anti-plane strain inclusion problem when the boundary effect is not considered. The problem of a circular cylindrical inclusion embedded concentrically in a finite cylindrical elastic matrix undergoing anti-plane strain deformations is analytically solved by applying the general solution, with the Eshelby tensor and its average over the circular cross section of the inclusion obtained in closed forms. This Eshelby tensor, being dependent on the position, inclusion size, matrix size, and a material length scale parameter, captures the inclusion size and boundary effects, unlike existing ones. It reduces to the classical linear elasticity-based Eshelby tensor for the circular cylindrical inclusion in an infinite matrix if both the strain gradient and boundary effects are suppressed. Numerical results quantitatively show that the inclusion size effect can be quite large when the inclusion is small and that the boundary effect can dominate when the inclusion volume fraction is high. However, the inclusion size effect is diminishing with the increase of the inclusion size, and the boundary effect is vanishing as the inclusion volume fraction becomes sufficiently low.     

Spatial Estimates for the Constrained Anisotropic Elastic Cylinder

This paper establishes spatial estimates in a prismatic (semi-infinite) cylinder occupied by an anisotropic homogeneous linear elastic material, whose elasticity tensor is strongly elliptic. The cylinder is maintained in equilibrium under zero body force, zero displacement on the lateral boundary and pointwise specified displacement over the base. The other plane end is subject to zero displacement (when the cylinder is finite, say). The limiting case of a semi-infinite cylinder is also considered and zero displacement on the remote end (at large distance) is not assumed in this case. A first approach is developed by considering two mean-square cross-sectional measures of the displacement vector whose spatial evolution with respect to the axial variable is studied by means of a technique based on a second-order differential inequality. Conditions on the elastic constants are derived that show the cross-sectional measures exhibit alternative behaviour and in particular for the semi-infinite cylinder that there is either at least exponential growth or at most exponential decay. A second approach considers cross-sectional integrals involving the displacement and its gradient and furnishes information upon the spatial evolution, without restricting the range of strongly elliptic elastic constants. Such models are principally based upon a first-order differential inequality as well as on one of second order. The general results are explicitly presented for transversely isotropic materials and graphically illustrated for a cortical bone.     

三相有限延迟型非Fourier介质热力耦合传播模式和群速度

给出了三相有限延迟型非Fourier介质的热力耦合传播的基本解,得到传播速度及其界限的结论。正值对数温升率相应的基本解表示温度和变形的指数分布形式较一定速度的传播,其传播速度与对数温升率有关;各个对数温升率相应的传播模式的总和,不存在恒定的传播速度,而是两个群速度。两个群速度都是有界数集。在热力无耦合情况下,一个群速度蜕化为具有恒定数值的膨胀波速,即第一声速;另一个群速度蜕化为温度波波速。后者也是一个群速度,其上确界为CV型介质温度波速的0.6065倍。     

Asymptotic of the flow upon shock incidence on a wedge cavity

The asymptotic of the motion originating because of shock incidence on a wedge cavity in a metal is investigated as the wave amplitude tends to zero. It has been shown in [1] that the flow is hence divided into two domains. The principal term governing the flow in the first domain agrees with the acoustic approximation. The flow in the second domain is described by incompressible fluid equations in the principal term. Determination of the flow in the second domain is reduced herein to the solution of a singular nonlinear integral equation. A numerical solution is found for a series of values of the cavity aperture.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 129–138, May–June, 1972.     

内聚力模型裂纹问题分析的解析奇异单元

内聚力模型已经被广泛应用于需要考虑断裂过程区的裂纹问题当中,然而常用的数值方法应用于分析内聚力模型裂纹问题时还存在着一些不足,比如不能准确的给出断裂过程区的长度、需要网格加密等。为了克服这些缺点,论文构造了一个新型的解析奇异单元,并将之应用于基于内聚力模型的裂纹分析当中。首先将虚拟裂纹表面处的内聚力用拉格拉日插值的方法近似表示为多项式的形式,而多项式表示的内聚力所对应的特解可以被解析地给出。然后利用一个简单的迭代分析,基于内聚力模型的裂纹问题就可以被模拟出来了。最后,给出二个数值算例来证明本文方法的有效性。