基于新型裂尖杂交元的压电材料断裂力学研究

提出了一种裂尖邻域杂交元模型，将其与标准杂交应力元结合来求解压电材料裂纹尖 端的奇性电弹场和断裂参数的数值解．裂纹尖端杂交元的建立步骤为：1) 利用高次内插有限元特征法求解特征问题，得到反映裂尖奇异性电弹场状况的特 征值和特征角分布函数；2) 利用广义Hellinger-Reissner变分泛函以及特征问题的解来建立裂尖邻域杂交元模型．该 方法求解电弹场时，摒弃了传统有限元方法中裂尖奇异性场需要借助解析解的做法，也避免 了单纯有限元方法中需要在裂尖端部进行高密度单元划分．采用PZT5板中心裂纹问题 作为考核例，数值结果显示了良好的精确性．作为进一步应用，求解了含中心界面裂纹 的PZT4-PZT5两相压电材料的应力强度因子和电位移强度因子．所有的算例都考虑 了3种裂纹面电边界条件．     

Spline sectorial elements

In this paper, the quadratic and cubic splines local interpolation on a sectorial element in polar coordinates is discussed and a class of spline sectorial elements for analyses of plane and thin, problems are presented. A reasonable treatment of the assumed displacement fields for elements with nodes at the origin, (r=0) has been made so that the elements can not only characterize the geometrical properties at the origin but also remove the singularity of strains and stresses there. Some numerical examples are given to show the efficiency of the proposed elements.     

Comment on “A theorem for a shear stress-free plane boundary in Stokes flow”

The results proved in a recent paper by Akhtar and Sen [Mech. Res. Commun. 38 (2011) 529–531] are not new and are already known in the literature. In fact, the main expressions – presented in the form of a theorem – in the cited article follow as a special case of the more general results for a two-fluid Stokes flow with a planar interface proved earlier by Palaniappan [Acta Mechanica 139 (2000) 1–13]. We briefly demonstrate the reduction of particular cases of the solution given by Palaniappan which also reveals the duplication of another special solution elsewhere.     

Invariant solutions of the equations of a multicomponent charged-particle beam

The concept of the invariant-group solution (H-solution) was introduced and a general method for obtaining it was developed in [1–3]. The group properties of the equations of a monoenergetic charged-particle beam with the same value and sign of the specific charge, assuming univalency of the velocity vector V, were studied in [4–6], where all essentially different H-solutions were also constructed. Below, the results of [4–6] are extended to the case of a beam in the presence of a fixed background of density ₀ (§1), and also to the case of multivelocity (V is an s-valued function) and multicomponent beams (i.e., beams formed by particles of several kinds) (§2). A number of analytic solutions that describe some nonstationary processes in devices with plane, cylindrical, and spherical geometry —among them a continuous periodic solution for a plane diode with a period determined by the background density -are obtained in §1. A transformation that contains arbitrary functions of time and preserves Vlasov's equations is given (§2). The equations studied can be treated as the equations of a rarefied plasma in the magnetohydrodynamic approximation, when the pressure gradients are negligible as compared with forces of electromagnetic origin.     

Asymptotic theory of the flow around an obstacle by a sonic flow

The first investigation of the problem of the flow around an obstacle by a gas flow whose velocity is equal to the speed of sound at infinity was carried out in [1, 2], where it is shown in particular that the principal term of the appropriate asymptotic expansion is a self-similar solution of Tricomi's equation, to which the problem reduces in the first approximation upon a hodographic investigation. The requirement that the stream function be analytic as a function of the hodographic variables on the limiting characteristic was an important condition determining the selection of the self-similarity exponent n (xy^–n is an invariant of the self-similar solution). The analytic nature of the velocity field everywhere in the flow above the shock waves, which arise from necessity upon flow around an obstacle, follows from this condition. The latter was found in [3], where one of the branches of the solution obtained in [1] was used in the region behind the shock waves. The principal and subsequent terms of the asymptotic expansion describing a sonic flow far from an obstacle were discussed in [4], where the author restricted himself to Tricomi's equation. Each term of the series constructed in [4] contains an arbitrary coefficient (we will call it a shape parameter) which is not determined within the framework of a local investigation, and consideration of the problem of flow around a given obstacle as a whole is necessary in order to determine these shape parameters. It follows from the results of [4] that the problem of higher approximations to the solution of [1] coincides with the problem, of constructing a flow in the neighborhood of the center of a Laval nozzle with an analytic velocity distribution along the longitudinal axis (a Meyer-type flow). Along with the Meyer-type flow in the vicinity of the nozzle center, which corresponds to a self-similarity exponent n=2, two other types of flow are asymptotically possible with n=3 and 11, given in [5]. The appropriate solutions are written out in algebraic functions in [6]. The results of [5] show that the condition that the velocity vector be analytic on the limiting characteristic in the flow plane is broader than the condition that the stream function be analytic as a function of the hodographic variables, which is employed in [1, 2, 4]. Therefore, the necessity has arisen of reconsidering the problem of higher approximations for the obstacle solution of F. I. Frankl'. It has proved possible for the region in front of the shock waves to use a series which is more general than in [4], which implies the inclusion of an additional set of shape parameters. The solution is given in the hodograph plane in the form of the sum of two terms; the series discussed in [4] corresponds to the first one, and the series generated by the self-similar solution with n=3 or with n=11 corresponds to the second one.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 99–107, May–June, 1979.The authors thank S. V. Fal'kovich for a useful discussion.     

THE SOLUTIONS OF STEADY－STATE CONVECTION EQUATIONS IN THE SPACES THAT POSSESS RESTORI NGNUCLEUS

ＴＨＥＳＯＬＵＴＩＯＮＳＯＦＳＴＥＡＤＹ－ＳＴＡＴＥＣＯＮＶＥＣＴＩＯＮＥＱＵＡＴＩＯＮＳＩＮＴＨＥＳＰＡＣＥＳＴＨＡＴＰＯＳＳＥＳＳＲＥＳＴＯＲＩＮＧＮＵＣＬＥＵＳＺｈａｎｇＣｈｉ－ｐｉｎｇ（张池平）;ＣｕｉＭｉｎｇ－ｇｅｎ（崔明根）（Ｈａｒｂｉｎ...     

Calculating the flows of a viscoplastic medium in tubes using the local variation method

This article considers the problem of the motion of a visco-plastic medium in tubes and channels. The results of [1, 2] are used, which present the variational formulation and the qualitative analysis of this problem. The method of local variations suggested in [3] is used for the numerical solution of the variational problem. A more detailed presentation of the algorithm of this method in application to boundary-value and variational problems is given in [4]. Results of calculations of certain concrete problems on an electronic computer are presented.The author wishes to express his sincere gratitude to F. L. Chernous'ko for the problem formulation and helpful counsel, and G. I. Barenblatt and S. S. Grigoryan for useful discussions.     

Independence of the Presence of Cracks or Inclusions of the Net Interaction Force Acting on a Dislocation by a Skewed Line Singularity

Orlov and Indenbom [1] have shown that the net (integrated) interaction force F between two skew dislocations with Burgers vectors separated by a distance h in an infinite anisotropic elastic medium is independent of h. Nix [2] computed numerically the net interaction force F between two skew dislocations that are parallel to the traction-free surface X2=0 of an isotropic elastic half-space. His numerical results showed that F was independent of h; a partial result of what Barnett [3] called Nix"s theorem. The separation-independence portion of Nix"s theorem has been proved to hold for a general anisotropic elastic half-space with a traction-free, rigid, or slippery surface, and for bimaterials [3-5]. In this paper, we show that the net interaction force is independent of the presence of inclusions. We will consider the case in which the line dislocation b is a more general line singularity which can include a coincident line force with strength f per unit length of the line singularity. An inclusion is an infinitely long dissimilar anisotropic elastic cylinder of an arbitrary cross-section whose axis is parallel to the line singularity (f, b). The (skew) line dislocation does not intersect the inclusion. The special cases of an inclusion are a void, crack, or rigid inclusion. There can be more than one inclusion of different cross sections and different materials. The line singularity (f, b) can be outside the inclusions or inside one of the inclusions. The inclusions and the matrix need not have a perfect bonding. One can have a debonding with or without friction. This revised version was published online in August 2006 with corrections to the Cover Date.     

Hydrodynamic methods in the inverse problem of gravitational prospecting

In [1, 2], a dynamical method is proposed for solving stationary inverse problems of potential theory, including the inverse problem of gravitational prospecting. It is based on analogy with the problem of establishing the interface of two immiscible fluids flowing in a porous medium. In the present paper, a system of two functional equations is derived from which one can obtain, as special cases, an equation corresponding to the method of [1, 2], and also a system of equations that enables one to propose a new and different method for solving the inverse problem of gravitational prospecting. Equations are derived in polar coordinates for plane Cauchy problems corresponding to both methods, and the results are also given of the solution of some model problems by these methods. Finally, ways of generating new methods of solution of the inverse problem of gravitational prospecting are considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 63–71, July–August, 1980.     

Proper solutions and limit boundary value problems of nonlinear second-order systems of differential equations

For the system of differential equations x=r(t)y,y=-a(t)f(x)g(y) where a(t)>0, r(t)>0 for t≥t; f(x) >0 and is decreasing for x>0 g(y)>0, we give necessary and sufficient condition of the existence of a proper solution, a bounded proper solution or solutions of two kinds of boundary value problems on an infinite interval [c,∞] c≥t_g. Several examples are given to illustrate the conditions of these results.     

Acoustic resonance associated with vibrations of an array of plates in a subsonic flow

The acoustic resonance effect occurring in the vibrations of an array of profiles (cascade) in a gas flow has been studied by a number of authors [1–8]. Relying on assumptions of a heuristic nature [1, 2, 7, 8] and using rather crude models [4, 6], they have derived criteria governing the acoustic resonance regimes and given the effect a certain physical interpretation. However, many problems of a physical bearing with regard to the quantitative and qualitative principles of the effect have been left unresolved. For a more complete and rigorous solution of the problem the author has previously [9, 10] analyzed the natural modes of a gas flowing past an array of plates It was determined that in the array domain the vibrational modes of the gas are localized in the vicinity, of the array and the eigenvalues are determined by the characteristic dimensions of the interstitial channel (as an open resonator). Also, the eigenfrequences were determined for the gas in the flow plane with the array absent [9]. Under spatial periodic conditions, such that the flow in the plane can be considered as a certain model of flow in an annular duct, these eigenfrequencies concurred with those obtained earlier in [1, 2, 4–6]. The results of [9, 10] are used in the present study as a basis for investigating certain laws and relations governing the unsteady aerodynamic characteristics of arrayed profiles in or close to regimes such that the gas can execute natural vibrations in the array domain and in an annular duct in the absence of the array.Novosibirsk. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 139–144, September–October, 1972.     

Self-similar solution to a problem of axisymmetric motion of gas through a porous medium in the case of a quadratic law of resistance

The results of solution of the self-similar problem of planar flow of gas through a porous medium in the case of a quadratic law of resistance [1] are generalized to the case of axisymmetric motion. The equation in similarity variables for the velocity of isothermal gas flow is reduced to an equation having cylindrical functions as solution. Analytic dependences of the pressure and the gas velocity on the coordinate and time are obtained for a given flow rate of the gas at the coordinate origin and for zero Initial gas pressure in the porous medium.Translated from Izvestlya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza., No. 4, pp. 168–171, July–August, 1982.     

Supersonic flow of a combustible gas mixture past a sphere

In recent years considerable interest has developed in the problems of steady-state supersonic flow of a mixture of gases about bodies with the formation of detonation waves and slow combustion fronts. This is due in particular to the problem of fuel combustion in a supersonic air stream.In [1] the problem of supersonic flow past a wedge with a detonation wave attached to the wedge apex is solved. This solution is based on using the equation of the detonation polar obtained in [2]-the analog of the shock polar for the case of an exothermic discontinuity. In [3] a solution is given of the problem of cone flow with an attached detonation wave, and [4] presents solutions of the problems of supersonic flow past the wedge and cone with the formation of attached adiabatic shocks with subsequent combustion of the mixture in slow combustion fronts. In the two latter studies two different solutions were also found for the problem of flow past a point ignition source, one solution with gas combustion in the detonation wave, the other with gas combustion in the slow combustion front following the adiabatic shock. These solutions describe two different asymptotic pictures of flow of a combustible gas mixture past bodies.In an experimental study of the motion of a sphere in a combustible gas mixture [5] it was found that the detonation wave formed ahead of the sphere splits at some distance from the body into an ordinary (adiabatic) shock and a slow combustion front. Arguments are presented in [6] which make it possible to explain this phenomenon and in certain cases to predict its occurrence.The present paper presents examples of the calculation of flow of a combustible gas mixture past a sphere with a detonation wave in the case when the wave does not split. In addition, the flow near the point at which the detonation wave splits is analyzed for the case when splitting occurs where the gas velocity behind the wave is greater than the speed of sound. This analysis shows that in the given case the flow calculation may be carried out without any particular difficulties. On the other hand, the calculation of the flow for the case when the point of splitting is located in the subsonic portion of the flow behind the wave (or in the region of influence of the subsonic portion of the flow) presents difficulties. This flow case is similar to the problem of the supersonic jet of finite width impacting on an obstacle.     

Transition from Planar to Whirling Oscillations in a Certain Nonlinear System

A single-mass two-degrees-of-freedom system is considered, witha radially oriented nonlinear restoring force. The latter is smooth andbecomes infinite at a certain value of a radial displacement. Stabilityanalysis is made for planar oscillation, or motion along a givendirection. As long as this motion is periodic, the nonlinearity in therestoring force provides a periodic parametric excitation in thetransverse direction. The linearized stability analysis is reduced tostudy of the Mathieu equation for the (infinitesimal) motions in thetransverse direction. For the case of free oscillations in the givendirection an exact solution is obtained, since a specific analyticalform is used for the (strongly nonlinear) restoring force, which permitsexplicit integration of the equation of motion. Stability of the planarmotion in this case is shown to be very sensitive to even slightdeviations from polar symmetry in the restoring force (as well as to theamplitude of oscillations in the given direction). Numerical integrationof the original equations of motion shows the resulting motion to be awhirling type indeed in case of the transversal instability. For thecase of a sinusoidal forcing in the given direction solution for the(periodic) response is obtained by Krylov–Bogoliubov averaging. Thisresults in the transmitted Ince–Strutt chart – namely, stabilitychart for transverse direction on the amplitude-frequency plane of theexcitation in the original direction.     

II型动态扩展裂纹尖端场的仿真分析

采用弹粘塑性力学模型,对II型动态扩展裂纹尖端的对数奇异性,进行了数值仿真计算。详细地分析了粘性系数α、马赫数M2对裂纹尖端的应力场影响。指出了文献[1]中对数奇异性区域存在的问题,解释了文献[1]中过度区的成因,对过度区尖端场解的形式和求解方法做了合理的推测。     

Plane stationary flows with ionizing shock waves in a magnetic field

The formulation of stationary, plane, and self-similar problems is considered when the flow parameters depend only on the polar angle, and the magnetic field lies in the flow plane. The case in which the magnetic field is perpendicular to the flow plane has been examined in [1]. The conditions are found under which the solution depends on an arbitrary parameter and the reasons for this nonuniqueness are explained. Self-similar solutions are constructed to describe the flow around an insulating wedge and a wall.     

极坐标系下Fourier-Legendre谱元方法与有限差分法数值扩散的比较

提出一种Fourier-Legendre谱元方法用于求解极坐标系下的Navier-Stokes方程,其中极点所在单元的径向采用Gauss-Radau积分点,避免了r=0处的1/r坐标奇异性。时间离散采用时间分裂法,引入数值同位素模型跟踪同位素的输运过程验证数值模拟的精度,分别利用谱元法和有限差分法的迎风差分格式求解匀速和加速坩埚旋转流动中的同位素方程。计算结果表明,有限差分法中的一阶迎风差分格式存在严重的数值假扩散,二阶迎风差分格式的数值结果较精确,增加节点可以有效地缓解数值扩散。然而,谱元法具有以较少节点得到高精度解的优势。     

Fluid flow in an arbitrary cascade on an axisymmetric stream surface in a variable thickness layer

A solution is given for the problem of flow past a cascade on an axisymmetric stream surface in a layer of variable thickness, which is a component part of the approximate solution of the three-dimensional problem for a three-dimensional cascade. Generalized analytic functions are used to obtain the integral equation for the potential function, which is solved via iteration method by reduction to a system of linear algebraic equations. An algorithm and a program for the Minsk-2 computer are formulated. The precision of the algorithm is evaluated and results are presented of the calculation of an example cascade.In the formulation of [1, 3] the problem of flow past a three-dimensional turbomachine cascade is reduced approximately to the joint solution of two-dimensional problems of the averaged axisymmetric flow and the flow on an axisymmetric stream surface in an elementary layer of variable thickness.In the following we solve the second problem for an arbitrary cascade with finite thickness rotating with constant angular velocity in ideal fluid flow: the solution is carried out on a Minsk-2 computer.Many studies have been devoted to this problem. A method for solving the direct problem for a cascade of flat plates in a hyperbolic layer was presented in [2]. Methods were developed in [1, 3] for constructing the flow for the case of a channel with variable thickness; these methods are approximately applicable for dense cascades but yield considerable error for small-load turbomachine cascades. The solution developed in [4], somewhat reminiscent of that of [2], is applicable for thin, slightly curved profiles in a layer with monotonically varying thickness. A solution has been given for a circular cascade for layers varying logarithmically [5] and linearly [6]. Approximate methods for slightly curved profiles in a monotonically varying layer with account for layer variability only in the discharge component were examined in [7–9]. A solution is given in [10] for an arbitrary layer by means of the relaxation method, which yields a roughly approximate flow pattern. The general solution of the problem by means of potential theory and the method of singularities presented in [11] is in error because of neglect of the crossflow through the skeletal line. The computer solution of [12] contains an unassessed error for the calculations in an arbitrary layer. The finite difference method is used in [13] to solve the differential equation of flow, which is illustrated by numerical examples for monotonie layers of axial turbomachines. The numerical solution of [13] is very complex.The solution presented below is found in the general formulation with respect to the geometric parameters of the cascade and the axisymmetric surface and also in terms of the layer thickness variation law.The numerical solution requires about 15 minutes of machine time on the Minsk-2 computer.     

Calculation of supersonic viscid flow near stagnation line of a blunt body

A numerical study is made of supersonic flow of a viscous gas in the vicinity of the stagnation line of plane and axisymmetric blunt bodies (cylinder, sphere). As in [1–5], which consider the compressed layer of a viscous gas in the vicinity of the stagnation point, use is made of the locally self-similar approximation, which is used to transform the Navier-Stokes equations into a system of ordinary differential equations. In the present paper the solution is sought with the simplifications of [5] and with more general conditions, which makes it possible to study a broad class of flows. The proposed numerical algorithm permits obtaining the structure of the compressed layer near the stagnation line, including the shock wave and the boundary layer. The calculations made on a computer for different flow conditions are illustrated by graphs.The author wishes to thank G. I. Petrov, G. F. Telenin, and L. A., Chudov for their interest in the study and for their helpful discussions. discussions.     

On several new methods to polar decomposition computation

In this paper, the polar decomposition of a deformation gradient tensor is analyzed in detail. The four new methods for polar decomposition computation are given: (1) the iterated method, (2) the principal invariant's method, (3) the principal rotation axis's method, (4) the coordinate transformation's method. The iterated method makes it possible to establish the nonlinear finite element method based on polar decomposition. Furthermore, the material time derivatives of the stretch tensor and the rotation tensor are obtained by explicit and simple expressions. The authors gratefully acknowledge the support rendered by the National Natural Science Foundation of China and the Natural Science Foundation of Jiangxi of China in 1998.