On the torsion of a cylinder with several cracks

In this paper, based on paper [1], the analytic expression of the torsion function for a cylinder containing arbitrary oriented cracks is obtained. The problem is reduced to solve a system of singular integral equations for the unknown dislocation density functions. Using the numerical method of the singular integral equations^[2,7] the torsional rigidities and stress intensity factors are evaluated for several multicracked cylinders. Next, the creak-cutting method^[5] is firstly extended to lve the torsion problem for a rectangular prism. The numerical results show that the method presented here is successful. Projects Supported by the Science Fund of the Chinese Academy of Sciences.     

Stress rate integral equations of elastoplasticity

The stress rate integral equations of elastoplaticity are deduced based on Ref. [1] by consistent methods. The point at which the stresses and/or displacements are calculated can be in the body or on the boundary, and in the plastic region or elastic one. The existence of the principal value integral in the plastic region is demonstrated strictly, and the theoretical basis is presented for the paticular solution method by unit initial stress fields. In the present method, programming is easy and general, and the numerical results are excellent. The project supported by the National Natural Science Foundation of China     

应力梯度不均匀时平板表面裂纹的应力强度因子

本文采用考虑裂纹面上具有任意分布载荷的线弹簧模型,在Kirchhoff板弯曲理论的假设下,将含半椭圆型表面裂纹的平板问题化为一组耦合的积分方程组进行求解,对均匀拉伸和纯弯曲两种载荷作用下的应力强度因子数值解,同经典线弹簧模型和有限元解进行了比较,并给出了经典线弹簧模型不能得到的、裂纹面上承受幂次不均匀应力分布时应力强度因子的数值解.     

Calculation of stresses and deformations of bellows by initial parameter method of numerical integration

By reducing the boundary value problem in stress analysis of bellows into initial value problem, this paper presents a numerical solution of stress distribution in semi-circular arc type bellows based upon the toroidal shell equation of V. V. Novozelov^[8]. Throughout the computation, S. Gill’s method^[1O] of extrapolation is used. The stresses and deformations of bellows under axial load and internal pressure are c-alculated, the results of which agree completely with those derived from the general solution of Prof. Chien Wei-zang^[1-4]. The extrapolation formula presented in this paper greatly promotes the accuracy of discrete calculation.The computer program in BASIC language of Wang 2200 VS computer is included in the appendix.     

Mode I crack tip field with strain gradient effects

A plane strain mode I crack tip field with strain gradient effects is investigated. A new strain gradient theory is used. An elastic-power law hardening strain gradient material is considered and two hardening laws, i. e. a separation law and an integration law are used respectively. As for the material with the separation law hardening, the angular distributions of stresses are consistent with the HRR field, which differs from the stress results^[19]; the angular distributions of couple stresses are the same as the couple stress results^[19]. For the material with the integration law hardening, the stress field and the couple stress field can not exist simultaneously, which is the same as the conclusion^[19], but for the stress dominated field, the angular distributions of stresses are consistent with the HRR field; for the couple stress dominated field, the angular distributions of couple stresses are consistent with those in Ref. [19]. However, the increase in stresses is not observed in strain gradient plasticity because the present theory is based on the rotation gradient of the deformation only, while the crack tip field of mode I is dominated by the tension gradient, which will be shown in another paper. Supported by the National Science Foundation of China (No. 19704100), Science Foundation of Chinese Academy of Sciences (Project KJ951-1-20), CAS K. C. Wong Post-doctoral Research Award Fund and the Post Doctoral Science Fund of China.     

高精度九节点非协调曲边三角形平面应力单元

本文利用[1]的方法,构造了一个九节点非协调三角形平面单元.与一般有限元相比可以提高一阶收敛精度,应力可直接在单元节点上得到.形成单刚矩阵时,不需要在单元域内进行数值积分,容易构造曲边单元.文末的算例表明,仅用很少的单元,位移和应力即可获得较高的精度.     

Calculation of aerodynamic characteristics of arrays of profiles of arbitrary shape

It is well known that the problem on nonseparating potential flow of an incompressible fluid about an array of profiles reduces to an integral equation for a certain real function, determined on the contours of the profiles of the array. As such a function one can take, as was done, for instance, in [1–5], the relative velocity of the fluid on the profiles of the array. For arrays of profiles of arbitrary shape it is necessary to solve the corresponding integral equation numerically. In the particular examples of the calculation of aerodynamic arrays that are available [1–3] the numerical methods used were based on the approximate evaluation of contour integrals by rectangle formulas. As investigations showed, sizeable errors arose thereby in the approximate solution obtained, these being especially significant in the case of curved profiles of relatively small bulk. In the present paper a method for the numerical solution of the integral equation obtained in [5] is proposed. The method is based on the replacement of a profile of the array with an inscribed N polygon, the length of whose sides is of the order N^–1 and whose internal angles are close to . Convergence with increasing N of the numerical solution to an exact solution of the integral equations at the reference points is demonstrated. Examples of the calculation are given.Novosibirsk. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 105–112, March–April, 1972.     

Dynamic response of elastic layer on stiff foundation under time harmonic surface vertical concentrated load

In this paper, the line-load integral equation method proposed in reference [1] is first used for solving the elastodynamic problems. A set of one-dimensional regular integral equation is derived for calculating the dynamic response of elastic layer on stiff foundation under time harmonic surface vertical concentrated load. And the numerical solution of the integral equation is obtained.     

压电介质表面圆形电极边缘的电-弹场分析

基于线性压电理论，研究了表面附有圆形电极板的半无限压电介质电极边缘的电弹场奇异分布特性。利用Hankel变换技术，将空间轴对称混合边值问题简化为一组对偶积分方程的求解，并得到了应力分量和电位移分量的解析表达式。算例表明在电极边缘的邻域内存在剧烈的应力集中和电位移增强，足以造成介质的脆性断裂或介质击穿。并对这种类裂纹奇异性质进行了讨论。     

垂直磁压电材料界面三维裂纹的超奇异积分法

应用有限部积分概念和广义位移基本解，垂直于磁压电双材料界面三维复合型裂纹问题被转 化为求解一组以裂纹表面广义位移间断为未知函数的超奇异积分方程问题. 进而，通过主部 分析法精确地求得裂纹尖端光滑点附近的奇性应力场解析表达式. 然后，通过将裂纹表面 位移间断未知函数表达为位移间断基本密度函数与多项式之积，使用有限部积分法对超奇异 积分方程组建立了数值方法. 最后，通过典型算例计算，讨论了广义应力强度因子的变化规 律.     

Crystallographic orientation and size dependence of tension-compression asymmetry in molybdenum nano-pillars

Uniaxial tension and compression experiments on [0 0 1] and [0 1 1] oriented molybdenum nano-pillars exhibit tension-compression asymmetry, a difference in attained stresses in compression vs. tension, which is found to depend on crystallographic orientation and sample size. We find that (1) flow stresses become higher at smaller diameters in both orientations and both loading directions, (2) compressive flow stresses are higher than tensile ones in [0 0 1] orientation, and visa versa in [0 1 1] orientation, and (3) this tension-compression asymmetry is in itself size dependent. We attribute these phenomena to the dependence of twinning vs. antitwinning deformation on loading direction, to the non-planarity of screw dislocation cores in Mo crystals, and to the possibly lesser role of screw dislocations in governing nano-scale plasticity compared with bulk Mo.     

Viscoelastic layer under the action of a moving load

Recently, problems concerning the dynamic behavior of imperfect continuous media under various types of actions have been widely investigated. The method of Laplace transformation is very convenient for describing physical processes concerning unsteady phenomena. In viscoelastic media two complications are added: the representation of the properties of a medium depending on time, and the inversion of the obtained solutions containing this additional complication. Certain approximate methods of inversion in the analysis of viscoelastic stresses are discussed in [1]. In [2, 3] a discussion is given for an effective method of constructing the solution of unsteady problems for finite and for infinite imperfect media using auxiliary functions, and a solution is presented for a half-space. Making use of the idea of the inversion of transforms, discussed in [4], in [5] a solution is obtained and a complete picture is presented for the dynamics of the variation of the stress field in a viscoelastic half-space. In the present study we consider the action of a normal moving load that is suddenly applied to the free surface of a viscoelastic layer. By Laplace and Fourier integral transformations we obtain a solution in the form of a uniformly converging series based on longitudinal and transverse waves reflected in the layer. By means of inverting the transforms by the method discussed in [4, 5], we obtain an exact solution for the stress field in the medium under investigation. We consider the special case of a viscoelastic medium of Boltzmann type, for which we obtain a numerical realization of the solution on a digital computer.     

The distributions of the interlaminar stress in a laminated shell with free end boundary

This paper studies the free edge effect yielded by interlaminar stress in a laminated cylindrical shell made up of fiber reinforced layer [0°], [90°] and the isotropic material layer under axisymmetric thermal load or radial pressure. Both ends of the shell are in free boundary condition. The exact solution of the problem can be obtained by using the three-dimensional theory of elasticity. For illustration, the numerical laminar stresses in a double-layer laminated shell under thermal load or radial pressure are calculated.     

A self-consistent solution for nonlinear theory of the positive column in a magnetic field

Only the electron and ion gases were taken into account in all previous theories of the positive column of intermediately-low-pressure arc discharge with or without the longitudinal magnetic field, while the motion of neutral gas was neglected. In 1982, the authors^[1] presented a nonlinear theory of a positive column which indicated that the rotating velocities of neutral gas and ion gas were nearly equal, and the motion of neutral gas could not be ignored. They further discussed the problem of validity of Bohm's criterion. However, some of the parameters with which the computation was worked out in Ref. [1] were not correlated to the initial discharge parameters. In the present paper, two integral relations are supplemented, so that a complete mathematical formation of the problem is given. A convergent numerical solution is obtained by iteration and the solution of Ref. [1] turns out to be the first iteration approximation. It is shown that both functions and parameters obtained by self-consistent solution differ significantly from those obtained in the first iteration approximation. According to this paper the computation can be conducted when the initial discharge parameters are given, so this method could have certain practical applications.     

一个新的八节点非协调曲边四边形平面单元

本文在文[1]和文[2]的基础上,提出构造非协调有限元的新方法。该方法不用一般的变分原理,可适用任意变系数正定和非正定偏微分方程。利用这一方法得到一个新的八节点四边形平面应力单元。与一般有限元相比,位移和应力可提高一阶收敛精度。形成单刚矩阵时,不需要进行数值积分。单元之间的协调条件容易满足,文中给出收敛性证明。文末给出数值算例,表明利用本文的方法,应力和位移均可获得满意的数值精度。     

磁电热弹耦合材料三维多裂纹超奇异积分法

用超奇异积分方程法将多场耦合载荷作用下磁电热弹耦合材料内含任意形状和位置三维多裂纹问题转化为求解一以广义位移间断为未知函数的超奇异积分方程组问题,退化得到内含任意形状平行三维多裂纹问题的超奇异积分方程组;推导出平行三维多裂纹问题的裂纹前沿广义奇异应力场解析表达式、定义了广义(应力、应变能)强度因子和广义能量释放率;应用有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以平行双裂纹为例,通过典型算例,研究了广义(应力、应变能)强度因子随裂纹位置、裂纹形状及材料参数变化规律,得到裂纹断裂评定准则. 最后,分析了裂纹间干扰、屏蔽作用及其在工程实际中的应用.      

Existence and comparison results for solutions of random integral and differential equations

This paper is the continuation of [1]. In this paper, we give another criteria of the existence of solutions for nonlinear random Volterra integral. A comparison theorem and the existence of random extremal solutions are also obtained by using the notion of ordering with respect to a cone. Our theorems generalize the corresponding results of Vaughan^[2,3] and Lakshmikantham^[4,5].     

A new analytical-numerical method for calculating interacting stresses of a multi-hole problem under both remote and arbitrary surface stresses

Based on the elementary solutions and new integral equations, a new analytical-numerical method is proposed to calculate the interacting stresses of multiple circular holes in an infinite elastic plate under both remote stresses and arbitrarily distributed stresses applied to the circular boundaries. The validity of this new analytical-numerical method is verified by the analytical solution of the bi-harmonic stress function method, the numerical solution of the finite element method, and the analytical-numerical solutions of the series expansion and Laurent series methods. Some numerical examples are presented to investigate the effects of the hole geometry parameters (radii and relative positions) and loading conditions (remote stresses and surface stresses) on the interacting tangential stresses and interacting stress concentration factors (SCFs). The results show that whether the interference effect is shielding (k <1) or amplifying (k> 1) depends on the relative orientation of holes (α) and remote stresses (σ_x^∞, σ_y^∞). When the maximum principal stress is aligned with the connecting line of two-hole centers and σ_y^∞ <0.5σ_x^∞, the plate containing two circular holes has greater stability than that containing one circular hole, and the smaller circular hole has greater stability than the bigger one. This new method not only has a simple formulation and high accuracy, but also has an advantage of wide applications over common analytical methods and analytical-numerical methods in calculating the interacting stresses of a multi-hole problem under both remote and arbitrary surface stresses.     

折线裂纹和两相材料界面的相互作用

利用螺位错基本解建立了和界面相交的折线裂纹的Cauchy型积分方程，根据奇异积分方程理论，得出了确定折线裂纹和界面交点处的奇性应力指数的特征方程，以及交点处各角形域内的奇性应力，利用所得的交点处的奇性应力定义了折线裂纹和界面交点处的应力强度因子，对所得积分方程进行数值求解，可得裂纹端点以及裂纹和界面交点处的应力强度因子。     

ELASTICITY SOLUTIONS OF SPHERICALLY ISOTROPIC CONES UNDER CONCENTRATED LOADS AT APEX

ＥＬＡＳＴＩＣＩＴＹＳＯＬＵＴＩＯＮＳＯＦＳＰＨＥＲＩＣＡＬＬＹＩＳＯＴＲＯＰＩＣＣＯＮＥＳＵＮＤＥＲＣＯＮＣＥＮＴＲＡＴＥＤＬＯＡＤＳＡＴＡＰＥＸＤｉｎｇＨａｏｊｉａｎｇ（丁皓江）；ＺｏｕＤａｏｑｉｎ（邹道勤）；ＲｅｎＹｏｎｇｊｉａｎ；（任永坚）（...