Closed-form solutions for stress singularities at plane bi- and trimaterial junctions

Stress singularities at two-dimensional bi- and trimaterial junctions, consisting of dissimilar, homogeneous, isotropic and linear-elastic wedges under a plane strain state, are considered. The stresses formed at the vertex of this composite situation are analyzed by the complex variable method, based on an appropriate choice of the Kolosov potentials that are applicable in the vicinity of the vertex. In doing so, the identification of the singularity exponent is performed. On the basis of a couple of bi- and trimaterial configurations, it is demonstrated how to derive some closed-form analytical solutions for the orders of the stress singularities from the characteristic equation.     

Stress singularities resulting from various boundary conditions in angular corners of plates of arbitrary thickness in extension

The stress singularities in angular corners of plates of arbitrary thickness with various boundary conditions subjected to in-plane loading are studied within the first-order plate theory. By adapting an eigenfunction expansion approach a set of characteristic equations for determining the structure and orders of singularities of the stress resultants in the vicinity of the vertex is developed. The characteristic equations derived in this paper incorporate that obtained within the classical plane theory of elasticity (M.L. Williams’ solution) and also describe the possible singular behaviour of the out-of-plane shear stress resultants induced by various boundary conditions.     

Libration points of a rotating complexified triangle

Similarities and differences between the force fields of a classical real dipole and a complex dipole are analyzed. The complex dipole is a pair of points equipped with complex conjugate masses and situated in a complex domain. The results of this analysis are used in the problem of motion of a material point in the field of attraction of a triangle uniformly rotating in its plane about its center of mass. It is assumed that a complex dipole is assigned to each vertex of the triangle. The existence and stability of libration points are studied. In particular, it is shown that there exist libration points outside the plane of the triangle.     

Singular rigid/plastic solutions in anisotropic plasticity under plane strain conditions

Assuming a rigid plastic material model with arbitrary smooth yield criterion, it is shown that the plane strain solutions are singular in the vicinity of maximum friction surfaces. In particular, some components of the strain rate tensor and thus the equivalent strain rate approach infinity. It is also shown that the exact asymptotic representation of the solution near maximum friction surfaces depends on the shape of the yield contour in the Mohr stress plane.     

理想正交各向异性弹塑性材料平面应力条件下Ⅰ型静止裂纹尖端场

在本文中,以 Hill 的塑性理论为基础,详细地讨论了理想正交各向异性弹塑性材料,平面应力条件下Ⅰ型静止裂纹尖端场解。裂纹尖端应力场不包含应力间断线,但包含弹性区。分析的结果表明(i)对于平面应力静止裂纹问题,应力场解不是唯一的,场解中的自由参数必须由远场条件来确定。(ii)裂纹尖端的应力、应变的奇异性,无论是各向异性材料还是各向同性材料,都是相同的。但在各向异性材料中,各向异性参数影响着应力、应变的幅度和分布。     

Corner stress singularities in an FGM thin plate

Describing the behaviors of stress singularities correctly is essential for obtaining accurate numerical solutions of complicated problems with stress singularities. This analysis derives asymptotic solutions for functionally graded material (FGM) thin plates with geometrically induced stress singularities. The classical thin plate theory is used to establish the equilibrium equations for FGM thin plates. It is assumed that the Young’s modulus varies along the thickness and Poisson’s ratio is constant. The eigenfunction expansion method is employed to the equilibrium equations in terms of displacement components for an asymptotic analysis in the vicinity of a sharp corner. The characteristic equations for determining the stress singularity order at the corner vertex and the corresponding corner functions are explicitly given for different combinations of boundary conditions along the radial edges forming the sharp corner. The non-homogeneous elasticity properties are present only in the characteristic equations corresponding to boundary conditions involving simple support. Finally, the effects of material non-homogeneity following a power law on the stress singularity orders are thoroughly examined by showing the minimum real values of the roots of the characteristic equations varying with the material properties and vertex angle.     

Fracture in plates of finite thickness

Application of the plane theory of elasticity to planar crack or angular corner geometries leads to the concept of stress singularity and stress intensity factor, which are the cornerstone of contemporary fracture mechanics. However, the stress state near an actual crack tip or corner vertex is always three-dimensional, and the meaning of the results obtained within the plane theory of elasticity and their relation to the actual 3D problems is still not fully understood. In particular, it is not clear whether the same stress field as found from the well-known 2D solutions of the theory of elasticity do describe the corresponding stress components in a plate made of a sufficiently brittle material and subjected to in-plane loading, and what effect the plate thickness has. In the present study we adopt, so called, first order plate theory to attempt to answer these questions. New features of the elastic solutions obtained within this theory are discussed and compared with 2D analytical results and experimental studies as well as with 3D numerical simulations.     

复合材料切口应力奇性指数计算

提出一种计算广义平面应交状态下复合材料切口应力奇性指数的新方法.在切口尖端的位移幂级数渐近展开式被引入正交各向异性材料的物理方程后,将用位移表示的应力分量代入切口端部柱状邻域的线弹性理论控制方程,切口应力奇性指数的计算被转化为常微分方程组特征值的求解.采用插值矩阵法求解该常微分方程组,可一次性地获取切口尖端多阶应力奇性指数.本法适合平面和反平面应力场耦合或解耦的情形,并可退化计算裂纹或各向同性材料切口的应力奇性指数.算例表明,所提方法对分析复合材料切口应力奇性指数是一种准确有效的手段.     

Localization in elasto-plastic materials: Influence of the plasticity yield surface in biaxial loading conditions

The influence of the plasticity yield surface on the development of instabilities in plane plates in biaxial loading is analyzed in order to understand and simulate the localization pattern observed in an expanding hemisphere experiment. First, a criterion for the activation of slip bands is formulated in the form of a critical hardening coefficient: it is particularized to the Von Mises and Tresca surfaces. In the Von Mises case, the criterion gives a strongly negative hardening coefficient in biaxial loading conditions different from the ones of plane strain. In the Tresca case, the criterion is fulfilled for a perfectly plastic material in uniaxial and biaxial loading; besides, in equi-biaxial loading, two possible orientations for slip bands are exhibited; this can be understood, with a few approximations, by the existence of a vertex point on the Tresca yield surface which give additive degrees of freedom for the direction of the plastic strain rate. Second, the development of localization in the loading conditions met in an expanding hemisphere experiment is simulated using both plasticity yield surfaces; whereas the Von Mises simulation does not localize, the Tresca simulation exhibits a pattern composed of a network of shear bands of different orientations; this pattern is not far from the pattern observed experimentally.     

The strain rate intensity factor in the plane strain compression of thin anisotropic metal strip

Using an available analytic solution for instantaneous plane strain compression of a plastically anisotropic strip between two parallel plates the strain rate intensity factor is found assuming Hill’s quadratic yield criterion. The distribution of material properties is uniform. The effect of parameters characterizing plastic anisotropy of the strip on the magnitude of the strain rate intensity factor is demonstrated. A possibility to replace the strain rate intensity factor with the plastic work rate intensity factor is discussed. Singular behavior of the plastic spin in the vicinity of the friction surface is revealed and discussed.     

Viscoelastic properties in plane squeeze-film flows

In this paper, we consider the problem of plane squeeze-film flows in a kinematics based on the formalism of convected (moving and deforming) coordinates. The flows discussed are treated as instantaneous motions with proportional stretch histories (cf. [16]). Certain simplifications in the constitutive equations of an incompressible simple fluid (cf. [10]) have been achieved for moderately low Deborah numbers.Approximate solutions of plane flows are obtained either for slightly viscoelastic fluids or in a form valid in the vicinity of any arbitrarily chosen instant of time. The conditions of improved lubrication, leading to inequalities imposed on material constants or kinematic quantities, are discussed in detail. Also, the necessary conditions are discussed under which the time-dependent distance between the plates may decrease non-monotonically, showing some “bounces”.     

Dynamically induced fracture

Expressions are presented for the stresses generated in the vicinity of a crack tip when a plane transient tension-stress wave, which is incident under an arbitrary angle, is diffracted by the crack. The effect of ductility effects on the combined Mode I and Mode II fractures is examined by assuming that yielding is restricted to a thin strip in the plane of the crack. Within this zone of yielding the normal stress in a certain defined direction is set equal to a yield stress of the material. It is shown that for a step-stress wave the leading edge of the zone of yielding initially moves at a constant speed, which is computed. The time for rupture at the trailing edge of the zone of yielding is also determined.     

Boundary Integral Analysis of the Symmetric Dynamic Problem for an Infinite Bimaterial Solid with an Embedded Crack

The symmetric frequency domain problem for two ideally bonded elastic half-spaces with a perpendicular plane crack is considered. It is reduced to the boundary integral equation (BIE) with integration over the limited crack region. The contact conditions on the bimaterial interface are satisfied identically in the initial stage of obtaining the equation. After boundary element solution of the equation, the stress concentration in the vicinity of a penny-shaped crack under time-harmonic loading of constant amplitude is studied. The mode I stress intensity factors as functions of angular coordinate of a crack front point and wave number for various relations between the material parameters are computed. The crack depth relative to the bimaterial interface is determined, when the effect of the material dissimilarity on the crack can be neglected.     

Influence of the Longitudinal Shear Crack Tip Radius on the Feathering Structure

The initial phase of feather joint development in the vicinity of a turnpike longitudinal shear vertex is analyzed. Experiments with model materials demonstrate that the crack parameters and the distance between the cracks along the shear front in the primary echelon brittle fracture structure linearly depend on the shear radius. A model for the development of the primary echelon structure along the longitudinal shear front is proposed.     

A special type of codimension two bifurcation and unusual dynamics in a phase-modulated system with switched strategy

We investigate a new type of codimension two bifurcation and related dynamics in a phase-modulated system with switched strategy. Two curves intersect at a point and are called the crisis–Hopf bifurcation. At the critical crisis–Hopf vertex, a boundary crisis and Hopf bifurcation coincide. Metamorphosis of coexisting attractors can be observed in the vicinity of the vertex. Another novelty is that we discover some sets of measure zero with riddled holes in the neighborhood of the bifurcation point. These sets display some qualitative and quantitative features of riddled basins but they are essentially different from the riddled basin. It may provide a more comprehensive picture for unusual dynamical features.     

Contact problems for a two-layer elastic wedge

We obtain integral equations for plane contact problems for a two-layer wedge (composite) under three types of boundary conditions on one of its sides (absence of stresses, sliding, or rigid fixation). The composite consists of two wedges completely linked with each other, which have different opening angles and elasticity parameters. Using the symbols (Mellin transforms) of the kernels of integral equations for the two-layer wedge, one can derive the symbols of the kernels of integral equations for symmetric problems about a crack in a three-layer wedge or a three-layer strip and for contact problems for a two-layer strip (by passing to the limit in a special way). The complex zeros of the Mellin transform determine the asymptotics of the normal contact pressure at the corner point of the composite as the contact region approaches this point. It is important that this asymptotics is also preserved in three-dimensional contact problems as the die enters the edge of a two-layer wedge (outside the corner points of the die itself). Taking into account this asymptotics, we obtain solutions of the contact problems as the die enters the vertex of the composite. We show that by appropriately choosing the materials and the internal angle of the two-layer wedge one can avoid contact pressure oscillations at the vertex, which occur in the case of a homogeneous wedge and result in loss of contact. The contact pressure at the wedge vertex can be made zero for a composite, while for a homogeneous wedge with the same opening angle it increases unboundedly. We construct asymptotic solutions of the contact problems for a plane die located relatively close or to the vertex of a two-layer wedge or relatively far from the vertex. The asymptotic and other methods were earlier used to solve similar plane contact problems for a homogeneous wedge [1, 2]. In the case of sliding fixation of one of the sides of a plane homogeneous wedge, the closed solution of the contact problem is known for a die entering the corner point [3, p. 131]. Two-dimensional contact problems were studied for a truncated wedge [4] and for a wedge supported by a rod of equal resistance [5]. The out-of-plane shear vibrations were studied for wedge-shaped composites [6, 7]. The spatial contact problems were considered for a homogeneous wedge [8]. The plane contact problem was analyzed for a continuously inhomogeneous wedge one of whose sides was rigidly fixed (the shear modulus continuously depends on the angular coordinate and the Poisson coefficient is constant). For a two-layer composite, which is studied in the present paper, the kernel symbol has different asymptotic properties, which are used in asymptotic methods for solving the problem. A similar distinction of the symbol properties takes place in contact problems for a continuously inhomogeneous layer and a layered packet.     

On dissolution driven crack growth

The formation and growth of a crack in a body subjected to stress driven material dissolution is studied. The rate of material dissolution is proportional to strain energy and curvature of the body surface. The formation of a crack from a plane surface is preceded by an evolving surface roughness. The continued dissolution enhances roughness amplitude resulting in pit formation. As the pit grows deeper into the material, it assumes the shape of a crack. The sharpness of the crack reaches its maximum during this transition from a pit to a crack. As the crack grows, a self-similar state is gradually assumed. During this phase characteristic lengths of the crack shape scale with the crack length. In line with this the crack progressively becomes blunt. The widest part of the crack when unloaded is in the vicinity of the crack tip. A consequence of the model is that no criterion is needed for crack growth. Neither is a criterion needed for determination of the crack path. It also follows that the crack growth rate is almost independent of the remote load. Further, spontaneous crack branching is anticipated. A motivation for this is given.     

基于单位分解的新型有限元方法研究

本文对基于单位分解概念的新型有限元方法具有一些新的特点进行了研究。这些新的特点包括自由度全部定义在单元顶点上,不同于传统的非协调和协调P型有限元方法自由度主要定义在单位区域内,因此很便于P型方法的实施;网格畸变度对解的不敏感性较传统有限元有了很大改进。此外,Dirichlet边界条件的处理与有限元是完全相同的,非常便利。论文分析了几个平面弹性问题：R.L.Tayor高阶分片试验和受拉弯剪切的平面悬臂梁,Cook悬臂T形板等,得到了满意的数值结果。     

Singularities in the Density Gradient

The compressible steady state Navier–Stokes system with viscosity is considered on a plane polygonal domain. It is shown that if a solution satisfies certain reasonable regularity conditions, and if the velocity is non-zero and points into the domain at a vertex of the polygon, then the gradient of the density either becomes infinite or has a jump discontinuity at each point of the streamline emanating from the vertex.     

中厚扁壳断裂问题的特征根及特征函数

证明了含任意切口、任意切口边界的多材料中厚扁壳问题的特征根等于相应的平面切口问题和平板弯曲切口问题两部分的特征根组合。进而证明了中厚扁壳切口问题的特征根等于相应反平面切口问题和平面切口问题的组合。中厚扁壳切口问题的特征根及其对应的0级特征函数均可直接按相应的两类基本问题（反平面切口问题和平面切口问题）进行求解。