Antiplane shear fields near a crack-tip in a brittle damaged material

In this paper, a simplified brittle damage model is proposed according to the Mazarz-Lemaitre damage model for concrete. A closed-form solution for a mode III crack is obtained based on the simplified model under small scale damage conditions, which allows for discontinuities of displacement-gradient and tangential stress on the damage boundary. It is pointed out that the discontinuities of field variables near the tip region exist for the brittle damaged material induced by the softening effect of the material. The preoject supported by the National Natural Science Foundation of China     

Antiplane deformation of a partially bonded elliptical inclusion

A new method that introduces two holomorphic potential functions (the two-phase potentials) is applied to analyze the antiplane deformation of an elliptical inhomogeneity partially-bonded to an infinite matrix. Elastic fields are obtained when either the matrix is subject to a uniform longitudinal shear or the inhomogeneity undergoes a uniform shear transformation. The stress field possesses the square-root singularity of a Mode III interface crack, which, in the special case of a rigid line inhomogeneity, changes in order, as the crack tip approaches the inhomogeneity end. In the latter situation the crack-tip elastic fields are linear in two real stress intensity factors related to a strong and a weak singularity of the stress field.     

Higher order fields for damaged nonlinear antiplane shear notch, crack and inclusion problems

Damaged nonlinear antiplane shear problems with a variety of singularities are studied analytically. A deformation plasticity theory coupled with damage is employed in analysis. The effect of microscopic damage is considered in terms of continuum damage mechanics approach. An exact solution for the general damaged nonlinear singular antiplane shear problem is derived in the stress plane by means of a hodograph transformation, then corresponding higher order asymptotic solutions are obtained by reversing the stress plane solution to the physical plane. As example, traction free sharp notch and crack, rigid sharp wedge and flat inclusion, and mixed boundary sharp notch problems are investigated, respectively. Consequently, higher order fields are obtained, in which analytical expressions of the dominant and second order singularity exponents and angular distribution functions of the near tip fields are derived. Effects of the damage and hardening exponents of materials and the geometric angle of notch/wedge on the near tip quantities are discussed in detail. It is found that damage leads to a weaker dominant singularity of stress, but to little stronger singularities of the dominant and second order terms of strain compared to that for undamaged material. It is also seen that damage has important effect on the angular distribution functions of the near tip stress and strain fields. As special cases, higher order analytical solutions of the crack and rigid flat inclusion tip fields are obtained, respectively, by reducing the notch/wedge tip solutions. Effects of damage and hardening exponents on the dominant and second order terms in the solutions of the crack and inclusion tip fields are discussed.     

Plasticity of damaged solids and shear band localization

Summary The main objective of the paper is the investigation of shear band localization conditions for finite elastic-plastic rate independent deformations of damaged solids. The first part of the paper is devoted to the formulation of the constitutive relations for elastic-plastic solids when isotropic and kinematic hardening effects and the micro-damage process are taken into consideration. The isotropic work-hardening effect is incorporated in the theory directly by defining the work-hardening-softening material function while the kinematic hardening effect and the softening effect generated by the micro-damage process are described by means of the internal state variable method. The second part of the paper aims at the investigation of the localization of plastic deformations. Different effects on the localization phenomenon are investigated. Particular attention is focused on kinematic hardening and micro-damage effects. It has been found that the influence of these two cooperative phenomena on the onset of localization within shear bands has synergetic nature. The results obtained are in good agreement with recent experimental observations.

---

Plastizität von geschädigten Feststoffen und Lokalisierung in Scherzonen

Übersicht Hauptgegenstand der Arbeit ist die Untersuchung der Bedingungen, die bei großen elastischplastischen Formänderungen von geschädigten, formänderungsgeschwindigkeitunabhängigen Feststoffen zur Lokalisierung in Scherzonen führen. Der erste Teil dient der Formulierung des Stoffgesetzes für elastischplastische Werkstoffe mit isotroper und kinematischer Verfestigung sowie Mikro-Schädigung. Die isotrope Verfestigung wird unmittelbar durch eine Verfestigungs-Entfestigungsfunktion berücksichtigt, während die kinematische Verfestigung und die Entfestigung infolge Mikro-Schädigung durch innere Zustandsgrößen beschrieben werden. Der zweite Teil befaßt sich mit der Lokalisierung der plastischen Formänderung, wobei verschiedene Einflüsse untersucht werden. Besondere Aufmerksamkeit wird auf die kinematische Verfestigung und Mikro-Schädigung gerichtet. Es stellt sich heraus, daß beide Erscheinungen bei der Lokalisierung in Scherzonen zusammenwirken. Die Ergebnisse stehen in guter Übereinstimmung mit neueren experimentellen Beobachtungen.

     

Antiplane shear crack in an infinite plate with a circular inclusion

Summary The effect of an elastic circular inclusion of a different material on the stress state of a single-cracked infinite sheet subjected to anti-plane shear is investigated. The proposed method makes use of complex variables in conjunction with Muskhelishvili's technique and has given interesting results for the stress field in the cracked plate. Numerical results, derived from this technique, investigate the influence of the regional inhomogeneity, produced by the inclusion in the structure, on the variation of stress-intensity factors at the crack tips.

---

Nichtebener Schubspannungsriß in einer unendlichen Platte mit einem kreisförmigen Einschluß

Übersicht Es wird der Einfluß untersucht, den ein elastischer, kreisförmiger Einschluß aus verschiedenem Material in einer unendlichen Platte mit einem Riß, welche einer longitudinalen Schubspannung unterworfen ist, auf den Spannungszustand der Platte ausübt. Die vorgeschlagene Methode benutzt komplexe Variablen in Zusammenhang mit der Muskhelishvilischen Technik und gibt für den Spannungszustand in der Platte interessante Resultate. Numerische Ergebnisse wurden gefunden, um den Einfluß der Inhomogenität auf die Spannungsintensitätsfaktoren des Risses zu bestimmen.

     

Three dimensional structures of the geometrically necessary dislocations in matrix-inclusion systems under uniaxial tensile loading

Plastic slip deformation in matrix-inclusion systems, in which a cuboidal or spherical shaped inclusion is embedded in a softer matrix, are numerically analyzed by a finite element technique. Edge and screw components of the geometrically necessary dislocations on slip systems are evaluated for each finite element from the spatial gradient of plastic shear strain. The character of the dislocation segments in each element is deduced from the data for edge and screw components and the directions of dislocation segments are determined. The aggregate of the dislocation segments in the whole specimen shows typical structures of dislocations, such as the Orowan loops around the inclusion and tilt boundaries that develop perpendicular to the primary slip plane. Stress state and shape of dislocations in deformable inclusions are discussed.     

Antiplane electro-mechanical field of a piezoelectric finite wedge under shear loading and at fixed-grounded boundary conditions

This paper presents the general solutions of antiplane electro-mechanical field solutions for a piezoelectric finite wedge subjected to a pair of concentrated forces and free charges. The boundary conditions on the circular segment are considered as fixed and grounded. Employing the finite Mellin transform method, the stress and electrical displacement at all fields of the piezoelectric finite wedge are derived analytically. In addition, the singularity orders and intensity factors of stress and electrical displacement can also be obtained. These parameters can be applied to examine the fracture behavior of the wedge structure. After being reduced to the problem of an antiplane edge crack or an infinite wedge in a piezoelectric medium, the results compare well with those of previous studies.     

Overall electromechanical properties of a binary composite with 622 symmetry constituents.: Antiplane shear piezoelectric state

A binary composite is studied here, where the electroelastic properties of the constituent materials belong to the crystal class 622. A square arrangement of long continuous circular cylinders, the fiber phase, embedded in a homogeneous medium is consider here. The composite is in a state of antiplane shear piezoelectricity, that is, a coupled state of out-of-plane mechanical displacement and in-plane electric field, which is characterized by three electroelastic parameters: longitudinal shear modulus, shear stress piezoelectric coefficient and transverse dielectric constant. Our interest here lies in the determination of its effective properties. They are derived by means of the method of two spatial scales. Closed-form expressions are obtained for them. Only one of the four local (or canonical) problems that arise is needed. Two properties are thus found. The Milgrom–Shtrikman compatibility relation is used to fix the remaining one. The local problem is solved using potential methods of a complex variable. The solution involves doubly periodic Weierstrass elliptic and related functions. The final formulae for the overall properties show explicitly the dependence on (i) the properties of the phases, (ii) the radius of the cylindrical fiber and (iii) the lattice sums associated with the square array. The shear modulus is shown to depend explicitly not only on the rigidity of the phases but also on their piezoelectric and dielectric coefficients. Some natural organic substances have the symmetry 622 like collagen. Recently Silva et al. measured its electroelastic properties. Their data is used to show some numerical results of the derived formulae as a function of the fiber volumetric fraction.     

Crack shielding and material deterioration in damaged materials: an antiplane shear fracture problem

Summary A simple damage evolution model is proposed for a quasibrittle material in the framework of continuum damage mechanics. The model is used to obtain a closed form solution for a mode-III stationary crack under small scale damage conditions. It is found that the crack tip stress intensity factor is reduced, i.e., the crack is shielded by the damage. However, this shielding effect is completely offset by the material deterioration caused by the damage. It also holds for steady state crack growth. When the most effective shielding is reached for the stationary crack, the zone dominated by the stress intensity factor shrinks to the crack tip. The results for the antiplane shear problem should shed some light on the in- plane fracture problem. Received 4 August 1997; accepted for publication 7 October 1997     

Complex dynamics of perfect discrete systems under partial follower forces

Equilibrium points, primary and secondary static bifurcation branches, and periodic orbits with their bifurcations of discrete systems under partial follower forces and no initial imperfections are examined. Equilibrium points are computed by solving sets of simultaneous, non-linear algebraic equations, whilst periodic orbits are determined numerically by solving 2- or 4-dimensional non-linear boundary value problems. A specific application is given with Ziegler's 2-DOF cantilever model. Numerous, complicated static bifurcation paths are computed as well as complicated series of periodic orbit bifurcations of relatively large periods. Numerical simulations indicate that chaotic-like transient motions of the system may appear when a forcing parameter increases above the divergence state. At these forcing parameter values, there co-exist numerous branches of bifurcating periodic orbits of the system; it is conjectured that sensitive dependence on initial conditions due to the large number of co-existing periodic orbits causes the chaotic-like transients observed in the numerical simulations.     

On equilibrium instability for conservative and partially dissipative systems

The investigation brings some contributions to the classical problem of inverting the Lagrange-Dirichlet stability theorem. First, an example is given of a conservative holonomic mechanical system with a stable equilibrium at the origin, although the potential function is strictly negative along some rays issuing from the origin. Then, one establishes a new instability result in the conservative case. Last, by means of a vector auxiliary function, one proves an instability theorem for holonomic systems with partial dissipation.     

Nonlinear analysis of shear deformable beam-columns partially supported on tensionless three-parameter foundation

In this paper, a boundary element method is developed for the nonlinear analysis of shear deformable beam-columns of arbitrary doubly symmetric simply or multiply connected constant cross-section, partially supported on tensionless three-parameter foundation, undergoing moderate large deflections under general boundary conditions. The beam-column is subjected to the combined action of arbitrarily distributed or concentrated transverse loading and bending moments in both directions as well as to axial loading. To account for shear deformations, the concept of shear deformation coefficients is used. Five boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to two stress functions and solved using the Analog Equation Method, a BEM-based method. Application of the boundary element technique yields a system of nonlinear equations from which the transverse and axial displacements are computed by an iterative process. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress functions using only boundary integration. The proposed model takes into account the coupling effects of bending and shear deformations along the member as well as the shear forces along the span induced by the applied axial loading. Numerical examples are worked out to illustrate the efficiency, wherever possible, the accuracy and the range of applications of the developed method.     

Superimposed steady and oscillatory shear in dispersed systems

Rheologica Acta - A viscometer has been designed to measure the behaviour of dispersed systems in oscillatory shear. The apparatus permits determination of the dynamic properties within the linear...     

用材料力学解的切应力求解弹性力学问题

分析了切应力与正应力的关系,讨论了导出切应力公式的条件,提出按切应力求解弹性力学问题的方法. 证明凡是$\sigma_y$与$x$无关的梁或者已知切应力为零的问题,按切应力求解都是可行的. 用这种方法求解比传统方法方便,运算简单. 该文给出了用切应力求解弹性力学平面问题的两个算例.     

Stochastic averaging and Lyapunov exponent of quasi partially integrable Hamiltonian systems

An n degree-of-freedom Hamiltonian system with r(1<r<n) independent first integrals which are in involution is called partially integrable Hamiltonian system and a partially integrable Hamiltonian system subject to light dampings and weak stochastic excitations is called quasi partially integrable Hamiltonian system. In the present paper, the averaged Itô and Fokker-Planck-Kolmogorov (FPK) equations for quasi partially integrable Hamiltonian systems in both cases of non-resonance and resonance are derived. It is shown that the number of averaged Itô equations and the dimension of the averaged FPK equation of a quasi partially integrable Hamiltonian system is equal to the number of independent first integrals in involution plus the number of resonant relations of the associated Hamiltonian system. The technique to obtain the exact stationary solution of the averaged FPK equation is presented. The largest Lyapunov exponent of the averaged system is formulated, based on which the stochastic stability and bifurcation of original quasi partially integrable Hamiltonian systems can be determined. Examples are given to illustrate the applications of the proposed stochastic averaging method for quasi partially integrable Hamiltonian systems in response prediction and stability decision and the results are verified by using digital simulation.     

Some problems in the antiplane shear deformation of bi-material wedges

The antiplane shear deformation of a bi-material wedge with finite radius is studied in this paper. Depending upon the boundary condition prescribed on the circular segment of the wedge, traction or displacement, two problems are analyzed. In each problem two different cases of boundary conditions on the radial edges of the composite wedge are considered. The radial boundary data are: traction–displacement and traction–traction. The solution of governing differential equations is accomplished by means of finite Mellin transforms. The closed form solutions are obtained for displacement and stress fields in the entire domain. The geometric singularities of stress fields are observed to be dependent on material property, in general. However, in the special case of equal apex angles in the traction–traction problem, this dependency ceases to exist and the geometric singularity shows dependency only upon the apex angle. A result which is in agreement with that cited in the literature for bi-material wedges with infinite radii. In part II of the paper, Antiplane shear deformation of bi-material circular media containing an interfacial edge crack is considered. As a special case of bi-material wedges studied in part I of the paper, explicit expressions are derived for the stress intensity factor at the tip of an edge crack lying at the interface of the bi-material media. It is seen that in general, the stress intensity factor is a function of material property. However, in special cases of traction–traction problem, i.e., similar materials and also equal apex angles, the stress intensity factor becomes independent of material property and the result coincides with the results in the literature.     

Fracture dynamics problems of orthotropic solids under anti-plane shear loading

By application of the approaches of the theory of complex functions, fracture dynamics problems of orthotropic solids under anti-plane shear loading were researched. Universal representation of analytical solutions was obtained by means of self-similar functions. The problems dealt with can be facilely transformed into Riemann–Hilbert problems by this technique, and analytical solutions of the stress, the displacement and dynamic stress intensity factor under the actions of moving increasing loads Px ²/t ² and Pt ³/x ² for the edges of asymmetrical mode III crack, respectively, were acquired. In the light of corresponding material properties, the variable rule of dynamic stress intensity factor was illustrated very well.     

Friedrichs systems for systems of wave equations and shear waves in a three-dimensional elastic medium

This paper considers a two-dimensional hydraulic model which describes gas or oil flow to a horizontal well with hydraulic fractures and takes into account the reservoir geometry and fluid flow between the reservoir and the well. A computational algorithm is proposed, and calculations for gas and oil reservoirs are performed. A comparison of the calculation results and the solutions of the corresponding problems in a three-dimensional formulation show that the calculations using the approximate hydraulic model yield reasonably accurate results.