Spin tests to determine brittle fracture under plane strain

The feasibility of using cyclic thermal stress or hydrostatic pressure to generate a fatigue crack in a large test-rotor blank is demonstrated. Test rotors, having test notches with fatigue-crack terminations, were spun to fracture to determine optimum test-notch design. It was found that, for the other test-notch dimensions held, it was necessary to extend the fatigue crack a minimum of 0.1 in. from the machined portion of the test notch to obtain a most effectively notched test rotor. In another series of tests, the influence of temperature on brittle-fracture strength of a Cr-Mo-V steel under plane strain was evaluated. It was found that, although there is a significant increase in fracture strength with increasing temperature, no knee in the curve is apparent in the vicinity of the conventionally measured transition temperatures ofNDT andT ₅₀. Also, design against brittle fracture is still required at temperatures aboveNDT andT ₅₀.Paper was presented at 1968 SESA Spring Meeting held in Albany, N. Y., on May 7–10.     

Plastic shearing of materials exhibiting strain hardening or strain softening

We consider the plastic shearing of a strain-rate dependent material exhibiting strain hardening or strain softening, subjected to steady shearing. We establish the existence of classical solutions and study the stability of uniform shearing. For materials exhibiting strain hardening or a moderate degree of strain softening we show that, as t , every solution approaches, at specific rates of convergence, uniform shearing; thus shear bands do not form.     

Dual-parameter elastic-plastic fracture criterion of metals and alloys in the plane strain case

In recent years, some investigators discussed the applicability of the HRR theory for engineering materials based on the results of numerical analyses and experimental studies. In the present paper, the finite element method is employed to analyze the crack tip fields of the engineering elastic-plastic material with a variety of geometry configurations of cracked specimens from elastic state to intensely general yielded state in the plane strain case. The results indicate that the HRR theory loses its validity of application for engineering elastic-plastic materials in the plane strain case. The reasons for this are analyzed. A dual-parameter fracture criterion is suggested for this case.     

Progressive failure constitutive model of fracture plane in geomaterial based on strain strength distribution

Progressive failure constitutive model of fracture plane in geomaterial based on strain strength distribution is proposed. The basic assumption is that strain strength of geomaterial comply with a certain distribution law in space. Failure of tensile fracture plane and shear fracture plane in representative volume element (RVE) with iso-strain are discussed, and generalized failure constitutive model of fracture plane in RVE is established considering combined effect of tension and shear. Fracture plane consists of elastic microplanes and fractured microplanes. Elastic microplanes are intact parts of the fracture plane, and fractured microplanes are the rest parts of the fracture plane whose strain have ever exceeded their strain strength. Interaction mode on elastic microplanes maintains linear elasticity, while on fractured microplanes it turns into contact and complies with Coulomb’s friction law. Intact factor and fracture factor are defined to describe damage state of the fracture plane which can be easily expressed with cumulative integration of distribution density function of strain strength. Strong nonlinear macroscopic behavior such as yielding and strain softening can be naturally obtained through statistical microstructural damage of fracture plane due to distribution of strain strength. Elastic–brittle fracture model and ideal elastoplastic model are special cases of this model when upper and lower limit of distribution interval are equal.     

Anisotropic work hardening of steel sheets under plane stress

This work is a follow-up of the previous report by Kim and Yin [Kim, K.H., Yin, J.J., 1997. Evolution of anisotropy under plane stress. J. Mech. Phys. Solids 45, 841–851] regarding the anisotropic work hardening of cold rolled steel sheets. Tensile prestrain has been applied at angles to the rolling direction and then tensile uniaxial yield stress and R-value distributions are measured. As reported earlier, the orientations of local maxima and minima in the yield stress are altered when the prestrain axis is not in the rolling direction. This led Kim and Yin [Kim and Yin (1997)] to suggest that the orientations of orthotropy axes are altered by the tensile prestrain at angles to the rolling direction. However, R-value distribution is found to be hardly affected by the prestrain. The unchanging R-value distribution shows that the material remembers the rolling direction even after the prestrain. An attempt is made to approximate the observed yield and flow behavior based upon isotropic-kinematic hardening with the quadratic yield function (Hill, 1948). The degree of approximation raises the issues of yield point definition, flexibility of yield function, non-associated flow rules, distortional hardening and others.     

A finite element analysis of strain localization for soft rock using a constitutive equation with strain softening

Summary It is well known that deformation and/or strain of geological materials localize when they come close to the failure state. In the present: study, we try to analyze the deformation problem using a constitutive relation with strain hardening and strain softening. The constitutive model of a soft rock and overconsolidated clay using an elasto-plastic constitutive theory with memory was originally developed by Adachi and Oka. This type of formulation is shown to be easily applied to analyze the material behavior of strain softening because there is a similarity to that in viscoplasticity. Using the proposed model, the loss of uniqueness of the solution to the initial value problem can be avoided and a special or complicated numerical technique, e.g., an arc length method, does not need to be used. When we use constitutive equations with strain softening in a finite element analysis, there is a problem of strong mesh size dependency of numerical results. To remedy the mesh size dependency, we generalize the Adachi-Oka model based on the concept of non-localization by Bazant. We apply the proposed constitutive model to the behavior of a sedimentary soft rock in the drained triaxial compression test. It is found that mesh size dependency becomes smaller using the non-localization of the constitutive model.

---

Eine Finite-Element-Analyse der lokalisierten Verformung von weichem Gestein bei einem Stoffgesetz mit Entfestigung

Übersicht Bekanntlich findet bei geologischen Stoffen in der Nähe des Versagenszustandes eine Lokalisierung der Verformung statt. Die Analyse dieses Problems wird in diesem Beitrag auf der Grundlage eines Stoffansatzes mit Ver- und Entfestigung unternommen. Entwickelt wurde das zugrunde gelegte elastischplastische Stoffmodell mit Gedächtnis von Adachi und Oka für weiches Gestein und übermäßig verdichteten Ton. Wegen der Ähnlichkeit zur Viskoplastizität läßt sich diese Formulierung des Stoffgesetzes leicht auf die Analyse des Verhaltens von entfestigendem Material anwenden, da der Eindeutigkeitsverlust der Lösung des Anfangswertproblems vermieden wird und besondere Rechenverfahren wie etwa die Bogenlängenmethode nicht benötigt werden. Bei der Benutzung von Stoffgesetzen in einer Finite-Element-Rechnung hängen die Ergebnisse stark von der Netzeinteilung ab. Um dies abzustellen, wird das Modell von Adachi und Oka auf der Grundlage von Bazant's Konzept der Nichtlokalisierung verallgemeinert. Anwendungsbeispiel ist das Verhalten von weichem Gestein im drainierten Triaxial-Test. Es zeigt sich, daß mit dem Konzept der Nichtlokalisierung im Stoffmodell der Einfluß der Netzeinteilung geringer wird.

---

---

Presented at the workshop on Numerical Methods for Localization and Bifurcation of Granular bodies, held at the Technical University of Gdask (Poland), September 25–30, 1989.     

Finite element solutions for plane strain mode I crack with strain gradient effects

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom ω_i are introduced in addition to the conventional three translational degrees of freedom u_i. ω_i is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l₁. Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale.     

Influence of damage on properties of strain localization in geomaterials at plane stress and plane strain

Summary By regarding geomaterials under loading as a mixture of intact and damaged parts, we investigate the influence of damage on the properties of strain localization in elastoplastic geomaterials at plane stress and plane strain. Conditions for the onset of strain localization including the effects of damage are derived for the cases of plane strain and plane stress. Discussed are the inclination of the localized band and the hardening modulus corresponding to the onset of strain localization. It is shown that the properties of the strain localization are dependent on the damage and the capacity of bearing hydrostatic pressure by the damaged part, and that damage may induce an earlier onset of strain localization and lead to instability of a geomaterial.accepted for publication 11 March 2004     

Modeling of plane strain ductile rupture

The formation of slanted fracture under plane strain conditions is studied using the finite Element (FE) method. Constitutive models proposed by Rousselier and by Gurson are used. Rice's condition for localization is checked at every point of the FE mesh for each time step. The role of mesh design (element size, element aspect ratio, symmetry) is first studied. The different constitutive models are then compared. It is in particular shown that the use of the f function in the Gurson model favors flat fracture.     

Plane strain pure bending of sheets with damage evolution at large strains

An analysis of plane-strain bending at large strains for the rigid/plastic incompressible material model including arbitrary strain-hardening and damage evolution laws is performed. The fracture criterion is based on a critical value of the damage parameter. Numerical treatment is reduced to the system of two partial differential equations written in characteristic coordinates. The through-thickness distribution of the principal stresses and damage parameter as well as the variation of the bending moment with the radius of curvature of the concave surface are found for Swift’s hardening law and one specific damage evolution law. General tendencies in solution behaviour are in agreement with physical expectations.     

应变软化对复合材料分层尖端场的影响

提出了适用于纤维增强复合材料应变软化分析的多标量连续损伤模型,并将其应用于含分层裂纹的复合材料层板的后屈曲损伤破坏有限元分析,研究了软化参数对损伤场和应力场的影响。计算结果表明:(1)应变软化使分层尖端的应力奇异性降低或消失。(2)应变软化参数影响极限损伤区的大小及其扩展速度,应力跌落使极限损伤加剧。(3)由于损伤的影响,裂尖的能量释放率受软化参数的影响出现波动,已难反映分层特征,可以根据损伤程度来判断分层的扩展。(4)分层上下表面的纤维方向影响损伤形式及其扩展方向。     

脆性固体损伤和局部软化的有限元分析

本文研究混凝土、岩石一类工程中常用的应变软化材料的有限元分析方法。在作者以往有关粘塑性损伤本构模型的工作基础上，给出了一组便于有限元计算的本构方程表达式。包括损伤弹性矩阵和局部损伤软化矩阵，分别运用于计算硬化和软化阶段的有限元刚度矩阵；对所提出的本构方程的实验验证计算和一些算例的有限元数值分析，表明文中给出的本构方程是可行的，相应的有限元算法能较好地对损伤固体的局部软化效应进行数值分析，并可成功地追踪应力应变响应的软化曲线     

Influence of softening on fracture propagation in large-scale shell structures

In this study, fracture propagation in large shell elements is modelled with the softening law. This law is given in a general form, enabling investigations of different softening behaviours to be conducted. The final fracture is simulated by removing elements. The softening parameters are derived using the energy-based representative volume element (RVE) approach. Tracing crack propagation through the RVE defines the physically justified softening parameters for the current model. The softening model is implemented into ABAQUS using VUMAT subroutines for the shell elements. A large-scale tearing experiment is simulated with the current model and RVE-based softening parameters. In addition, the softening laws from the literature have been used. The fracture propagation is assessed in terms of plastic energy dissipation in the RVE and the whole structure, load–displacement, and crack growth. The RVE-based model is shown to have better performance compared with other models from the literature.     

Prediction of the softening and damage effects with permanent set in fibrous biological materials

Deformation induced softening is an inelastic phenomenon frequently accompanying mechanical response of soft biological tissues. Inelastic phenomena which occur in mechanical testing of biological tissues are very likely to be associated with alterations in the internal structure of these materials.In this study, a novel structural constitutive model is formulated to describe the inelastic effects in soft biological tissues such as Mullins type behavior, damage and permanent set as a result of residual strains after unloading. Anisotropic softening is considered by evolution of internal variables governing the anisotropic properties of the material. We consider two weight factors w_i (softening) and s_k (discontinuous damage) as internal variables characterizing the structural state of the material. Numerical simulations of several soft tissues are used to demonstrate the performance of the model in reproducing the inelastic behavior of soft biological tissues.     

Plane-stress fracture testing of finite sheets under biaxial loads

A procedure which combines the Williams series-type stress- and displacement-field expressions at the crack-tip neighborhood with a suitable numerical scheme away from the crack-tip was employed in the determination of the plane-stress fracture properties of four finite 7076-T6 aluminum sheets containing cracks emanating from a circular hole under four biaxial loads. The compatibility of the analytical and numerical displacements at the nodal points along the boundary of the crack-tip neighborhood was utilized in formulating displacement-continuity expressions containing some undetermined constants which solution depends on the nature of the boundary loading conditions. By linear superposition of the displacement due to remote uniaxial load and the displacements due to remotely applied transverse load in the neighborhood of the crack-tip, biaxial-displacement-continuity expressions containing these important fracture properties—namely, the opening Mode I stress-intensity factorK, the nonsingular stress term associated with the stresses in the direction parallel to the plane of cracksA and the integration termB associated with the displacement in this direction—were evaluated. Because no known biaxial testing of this geometry had been reported prior to this research, the analytical procedure was used to select the optimum geometry required in a biaxial fracture test of a finite-sheet specimen containing cracks emanating from a circular hole. This geometric optimization of the specimen guaranteed uniformity of stress all over the volume of specimen and also made the alteration of the existing MTS test fixtures unnecessary. Four square sheets of 7075-T6 aluminum alloy containing a central hole with two collinear cracks emanating radially at the edge of the hole were then fabricated in accordance with the analytically determined geometric requirements. The biaxial fracture test was then conducted under four biaxial load factors (λ) of 0.0, 0.5, 1.0 and 1.5. The fracture toughness obtained in this research was compared with those reported for uniaxial loading of large panels. It was found that there is a good correlation between the reported fracture toughness and this work.     

Upper-bound approach to plane strain strip rolling with free deformation zones

Summary A recently developed extension of the upper-bound theorem is applied to the problem of rigid/plastic deformation in case where the surface configuration of the concerned body is previously unknown. By assuming a simple velocity field the strip rolling process is analysed by considering the free deformation zones in front of, and behind the roll gap. The results so obtained show that the transitional plastic deformation zone exists only in front of the roll gap without any bulge. Furthermore, we find that this zone becomes larger with increase of the strain sensitivity as well as the strain rate sensitivity of material.

---

Obere-Schranken-Abschätzung zum ebenen Bandwalzen mit freien Umformzonen

Übersicht Eine vor kurzem entwickelte Erweiterung des Satzes von der oberen Schranke für starr-plastische Umformvorgänge wird auf ein Problem angewendet, bei dem die Gestalt der Materialoberfläche im voraus unbekannt ist. Und zwar wird mittels eines elementaren Geschwindigkeitsansatzes der ebene Walzvorgang unter Berücksichtigung der freien Umformzonen vor und hinter dem Walzspalt ermittelt. Die berechneten Ergebnisse zeigen, daß nur vor dem Walzspalt eine Übergangszone der plastischen Umformung ohne Wulstbildung existiert und daß diese Zone größer wird mit anwachsender Verfestigungs- und Geschwindigkeitsabhängigkeit des Materials.

     

Plastic surface strain mapping of bent sheets by image correlation

A technique using a single CCD camera, a precision rotation/translation stage, a telecentric zoom lens, and digital image correlation software is described for measuring surface profiles and surface plastic strain distributions of a bent thin sheet. The measurement principles, based on both parallel and pinhole perspective projections, are outlined and the relevant mathematical equations for computing the profiles and displacement fields on a curved surface are presented. The typical optical setup as well as the experimental measurement and digital image correlation analysis procedure are described. The maximum errors in the in-plane and out-of-plane coordinates or displacements are about ±5 and ±25 μm, respectively, and the maximum errors in surface strain mapping are about 0.1% or less based on a series of evaluation tests on flat and curved sample surfaces over a physical field of view of 15.2 × 11.4 mm². As an application example, the shape and surface plastic strain distribution example, the shape and surface plastic strain distributions around a bent apex of a flat 2 mm thick automotive aluminum AA5182-O sheet, which underwent a 90° bend with three bend ratios of 2t, 1t, and 0.6t, are determined using the proposed technique.     

The effect of strain softening in the matrix material during void growth

A Finite element analysis has been employed to investigate the growth of an initially spherical void embedded in a cylinder of elastic-plastic material. The boundary displacement of this cylindrical cell is regulated by the value of a parameter α which controls the radial shrinkage of the cell as it elongates. A large strain analysis was used and results for both strain hardening and strain softening (after an appropriate amount of hardening has taken place) have been obtained. The effects of different mean tensile stresses, equivalent strains and initial void volume fractions have also been included. The numerical work shows relationships between the mechanical and geometrical variables that may favour ductile fracture by void coalescence or by shear decohesion.