Tension–torsion fracture experiments – Part II: Simulations with the extended Gurson model and a ductile fracture criterion based on plastic strain

An extension of the Gurson model that incorporates damage development in shear is used to simulate the tension–torsion test fracture data presented in Faleskog and Barsoum (2013) (Part I) for two steels, Weldox 420 and 960. Two parameters characterize damage in the constitutive model: the effective void volume fraction and a shear damage coefficient. For each of the steels, the initial effective void volume fraction is calibrated against data for fracture of notched round tensile bars and the shear damage coefficient is calibrated against fracture in shear. The calibrated constitutive model reproduces the full range of data in the tension–torsion tests thereby providing a convincing demonstration of the effectiveness of the extended Gurson model. The model reinforces the experiments by highlighting that for ductile alloys the effective plastic strain at fracture cannot be based solely on stress triaxiality. For nominally isotropic alloys, a ductile fracture criterion is proposed for engineering purposes that depends on stress triaxiality and a second stress invariant that discriminates between axisymmetric stressing and shear dominated stressing.     

An equivalent von Mises stress and corresponding equivalent plastic strain for elastic–plastic ordinary peridynamics

Meccanica - Simple formulas for calculating equivalent von Mises stress and von Mises effective plastic strain in an elastic–plastic ordinary peridynamic analysis are proposed. The equivalent...     

Ratchetting under tension–torsion loadings: experiments and modelling

This paper is concerned with the mechanical behaviour of 316 austenitic stainless steel under multiaxial loadings and particular attention is paid to ratchetting under tension–torsion non-proportional loadings. First, a series of uniaxial tests and biaxial tests has been carried out in order to calibrate five different cyclic plasticity models based on an isotropic hardening rule and a non-linear kinematic hardening rule. It is shown that this class of models gives quite good agreement between the experimental and numerical results. Second, another series of ratchetting tests has been carried out under tension–torsion loadings in order to test the prediction capacities of the previous models. It is shown that whereas the models have been calibrated with similar loading paths, four of the five selected models give poor predictions.     

Dynamic Separation of Resistance Spot Welded Joints: Part I—Experiments

An analysis for an impact system is presented. The results are then used to interpret the test data from dynamic separation of resistance spot welded joints. In this Part I of the investigation emphasis is placed on the design consideration, development of a test system and verification of the design from actual test data obtained from the test system. In addition, the inertia effect of a generic dynamic system is analyzed using the principle of rigid body dynamics. It is shown that the load recorded by a load cell could include both the load experienced by the test specimen and the inertia force generated from the mass and acceleration between the specimen and the load cell, when the load cell is placed on the fixed side of the test specimen. Impact fixtures designed for spot weld strength testing are then studied for the inertia effect.     

Interaction between a submerged evacuated cylindrical shell and a shock wave—Part I: Diffraction–radiation analysis

The interaction between a submerged elastic circular cylindrical shell and an external shock wave is addressed. A linear, two-dimensional formulation of the problem is considered. A semi-analytical solution is obtained using a combination of the classical analytical approach based on the use of the Laplace transform and separation of variables, and finite difference methodology. The study consists of two parts. Part I focuses on the simulation and analysis of the acoustic fields induced during the interaction. Both the diffraction (absolutely rigid cylinder) and complete diffraction–radiation (elastic shell) are considered. Special attention is paid to the lower-magnitude shell-induced waves representing radiation by the elastic waves circumnavigating the shell. The focus of Part II is on the numerical analysis of the solution. The convergence of the series solution and finite-difference scheme is analysed. The computation of the response functions of the problem is discussed as well, as is the effect of the bending stiffness on the acoustic field. The membrane model of the shell is considered to analyse such an effect, which, in combination with the models addressed in Part I, allows for the analysis of the evolution of the acoustic field around the structure as its elastic properties change from an absolutely rigid cylinder to a membrane. The results of the numerical simulations are compared to available experimental data, and a good agreement is observed.     

A micromechanics-based strain gradient damage model for fracture prediction of brittle materials – Part I: Homogenization methodology and constitutive relations

In this paper, we first describe a homogenization methodology with the aim of establishing strain gradient constitutive relations for heterogeneous materials. The methodology presented in this work includes two main steps. The first one is the construction of the average strain-energy density for a well-chosen RVE by using a homogenization technique. The second one is the transformation of the obtained average strain-energy density to that for the continuum. An important characteristic of this method is its self-consistency with respect to the choice of the RVE: the strain gradient constitutive law built by using the present method is independent of the size and the form of the RVE. In the frame of this homogenization procedure, we have constructed a strain gradient constitutive relation for a two-dimensional elastic material with many microcracks by adopting the self-consistent scheme. It was shown that the effective behavior of cracked solids depends not only on the crack density but also on the average crack size with which the strain gradient is associated. The proposed constitutive relation provides a starting point for the development of an evolution law of damage including strain gradient effect, which will be presented in the second part of this work.     

Ductile–brittle transitions in the fracture of plastically-deforming,adhesively-bonded structures. Part I: Experimental studies

Rate effects for adhesively-bonded joints in steel sheets failing by mode-I fracture and plastic deformation were examined. Three types of test geometries were used to provide a range of crack velocities between 0.1 and 5000 mm/s: a DCB geometry under displacement control, a wedge geometry under displacement control, and a wedge geometry loaded under impact conditions. Two fracture modes were observed: quasi-static crack growth and dynamic crack growth. The quasi-static crack growth was associated with a toughened mode of failure; the dynamic crack growth was associated with a more brittle mode of failure. The experiments indicated that the fracture parameters for the quasi-static crack growth were rate independent, and that quasi-static crack growth could occur even at the highest crack velocities. Effects of rate appeared to be limited to the ease with which a transition to dynamic fracture could be triggered. This transition appeared to be stochastic in nature, it did not appear to be associated with the attainment of any critical value for crack velocity or loading rate. While the mode-I quasi-static fracture behavior appeared to be rate independent, an increase in the tendency for dynamic fracture to be triggered as the crack velocity increased did have the effect of decreasing the average energy dissipated during fracture at higher loading rates.     

On the effect of Lüders bands on the bending of steel tubes. Part I: Experiments

In several practical applications hot-finished steel pipe that exhibits Lüders bands is bent to strains of 2–3%. Lüders banding is a material instability that leads to inhomogeneous plastic deformation in the range of 1–4%. This work investigates the influence of Lüders banding on the inelastic response and stability of tubes under rotation controlled pure bending. Part I presents the results of an experimental study involving tubes of several diameter-to-thickness ratios in the range of 33.2–14.7 and Lüders strains of 1.8–2.7%. In all cases the initial elastic regime terminates at a local moment maximum and the local nucleation of narrow angled Lüders bands of higher strain on the tension and compression sides of the tube. As the rotation continues the bands multiply and spread axially causing the affected zone to bend to a higher curvature while the rest of the tube is still at the curvature corresponding to the initial moment maximum. With further rotation of the ends the higher curvature zone(s) gradually spreads while the moment remains essentially unchanged. For relatively low D/t tubes and/or short Lüders strains, the whole tube eventually is deformed to the higher curvature entering the usual hardening regime. Subsequently it continues to deform uniformly until the usual limit moment instability is reached. For high D/t tubes and/or materials with longer Lüders strains, the propagation of the larger curvature is interrupted by collapse when a critical length is Lüders deformed leaving behind part of the structure essentially undeformed. The higher the D/t and/or the longer the Lüders strain is, the shorter the critical length. Part II presents a numerical modeling framework for simulating this behavior.     

Fluid velocity based simulation of hydraulic fracture: a penny shaped model—part I: the numerical algorithm

In the first part of this paper, a universal fluid velocity based algorithm for simulating hydraulic fracture with leak-off, previously demonstrated for the PKN and KGD models, is extended to obtain solutions for a penny-shaped crack. The numerical scheme is capable of dealing with both the viscosity and toughness dominated regimes, with the fracture being driven by a power-law fluid. The computational approach utilizes two dependent variables; the fracture aperture and the reduced fluid velocity. The latter allows for the application of a local condition of the Stefan type (the speed equation) to trace the fracture front. The obtained numerical solutions are carefully tested using various methods, and are shown to achieve a high level of accuracy.     

The characteristics of plastic flow and a physically-based model for 3003 Al–Mn alloy upon a wide range of strain rates and temperatures

The uniaxial compressive responses of 3003 Al–Mn alloy upon strain rates ranging from 0.001/s to about 10⁴/s with initial temperatures from 77 K to 800 K were investigated. Instron servohydraulic testing machine and enhanced split Hopkinson bar facilities have been employed in such uniaxial compressive loading tests. The maximum true strain up to 80% has been achieved. The following observations have been obtained from the experimental results: 1) 3003 Al–Mn alloy presents remarkable ductility and plasticity at low temperatures and high strain rates; 2) its plastic flow stress strongly depends on the applied temperatures and strain rates; 3) the temperature history during deformation strongly affects the microstructure evolution within the material. Finally, paralleled with the systematic experimental investigations, a physically-based model was developed based on the mechanism of dislocation kinetics. The model predictions are compared with the experimental results, and a good agreement has been observed.     

A micromechanics-based strain gradient damage model for fracture prediction of brittle materials – Part II: Damage modeling and numerical simulations

In this paper, we established a strain-gradient damage model based on microcrack analysis for brittle materials. In order to construct a damage-evolution law including the strain-gradient effect, we proposed a resistance curve for microcrack growth before damage localization. By introducing this resistance curve into the strain-gradient constitutive law established in the first part of this work (Li, 2011), we obtained an energy potential that is capable to describe the evolution of damage during the loading. This damage model was furthermore implemented into a finite element code. By using this numerical tool, we carried out detailed numerical simulations on different specimens in order to assess the fracture process in brittle materials. The numerical results were compared with previous experimental results. From these studies, we can conclude that the strain gradient plays an important role in predicting fractures due to singular or non-singular stress concentrations and in assessing the size effect observed in experimental studies. Moreover, the self-regularization characteristic of the present damage model makes the numerical simulations insensitive to finite-element meshing. We believe that it can be utilized in fracture predictions for brittle or quasi-brittle materials in engineering applications.     

Mixed mode fracture of a piezoelectric–piezomagnetic bi-layer structure with two un-coaxial cracks parallel to the interface and each in a layer

The purpose of the present work is to study the mixed mode fracture of a piezoelectric–piezomagnetic composite with two un-coaxial cracks parallel to the interface and each in a layer. Methods of generalized dislocation simulation, Green’s function, Cauchy singular integral equation and Lobatto–Chebyshev collocation are combined together to get the numerical results of mechanical strain energy release rate (MSERR). Three kinds of effects are revealed by parametric studies, i.e., the free-surface effect, the shielding effect and the interference effect, and they are used to interpret the characteristics of COD and MSERR curves. In addition, the effects of shear loading, magnetic loading and electric loading on MSERR are also disclosed, respectively, by varying the corresponding loading factor.     

Collaborative optimization of screening efficiency and screen surface load — Part I: Modeling and evaluation

The screen surface load (SSL) caused by granular materials is an important factor affecting the structural performance of vibrating screen. Based on virtual experiment, a multi-objective collaborative optimization method is proposed to control the SSL under high screening efficiency (SE) in this work. Firstly, a DEM model was established to study the influence of process parameters on SE and SSL. Secondly, the NSGA-II (Non-dominated Sorting Genetic Algorithm) was employed to optimize the screening parameters with both SE and SSL as targets. The optimization method proves to be effective implementing on a linear vibrating screening. With SE equals to 98.5%, the SSL optimizable range is 39.2%. While compromising the SE to 88.7%, the SSL optimizable range improves to 48.6%. The result shows that the collaborative optimization could effectively control the SSL while maintaining a high SE, which is of great significance to improve the service life of screen surface and screen body.     

Analytical study of the nonlinear behavior of a shape memory oscillator: Part I—primary resonance and free response at low temperatures

In this work, the response of a single-degree-of-freedom shape memory oscillator subjected to the excitation harmonic has been investigated. Equation of motion is formulated assuming a polynomial constitutive model to describe the restitution force of the oscillator. Here the method of multiple scales is used to obtain an approximate solution to the equations of the motion describing the modulation equations of amplitude and phase, and to investigate theoretically its stability. This work is presented in two parts. In Part I of this study we showed the modeling of the problem where the free vibration of the oscillator at low temperature is analyzed, where martensitic phase is stable. Part I also presents the investigation dynamics of the primary resonance of the pseudoelastic oscillator. Part II of the work is focused on the study in the secondary resonance of a pseudoelastic oscillator using the model developed in Part I. The analysis of the system in Part I as well as in Part II is accomplished numerically by means of phase portraits, Lyapunov exponents, power spectrum and Poincare maps. Frequency-response curves are constructed for shape memory oscillators for various excitation levels and detuning parameter. A rich class of solutions and bifurcations, including jump phenomena and saddle-node bifurcations, is found.     

A second strain gradient elasticity theory with second velocity gradient inertia – Part I: Constitutive equations and quasi-static behavior

A multi-cell homogenization procedure with four geometrically different groups of cell elements (respectively for the bulk, the boundary surface, the edge lines and the corner points of a body) is envisioned, which is able not only to extract the effective constitutive properties of a material, but also to assess the “surface effects” produced by the boundary surface on the near bulk material. Applied to an unbounded material in combination with the thermodynamics energy balance principles, this procedure leads to an equivalent continuum constitutively characterized by (ordinary, double and triple) generalized stresses and momenta. Also, applying this procedure to a (finite) body suitably modelled as a simple material cell system, in association with the principle of the virtual power (PVP) for quasi-static actions, an equivalent structural system is derived, featured by a (macro-scale) PVP having the typical format as for a second strain gradient material model. Due to the surface effects, the latter model does work as a combination of two subsystems, i.e. the bulk material behaving as a Cauchy continuum, and the boundary surface operating as a membrane-like boundary layer, each subsystem being in (local and global) equilibrium by its own. Further, the applied (ordinary) boundary traction splits into two (response-dependent) parts, i.e. the “Cauchy traction” transmitted to the bulk material and the “Gurtin–Murdoch traction” acting, together with all other boundary tractions, upon the boundary layer. The role of the boundary layer as a two-dimensional manifold enclosing a Cauchy continuum is elucidated, also with the aid of a discrete model. A strain gradient elasticity theory is proposed which includes a minimum total potential energy principle featuring the relevant boundary-value problem for quasi-static loads and its (unique) solution. A simple application is presented. Two appendices are included, one reports the proof of the global equilibrium of the boundary layer, the other is concerned with double and triple stresses. The paper is complemented by a companion Part II one on dynamics. Previous findings by the author [Polizzotto, C., 2012. A gradient elasticity theory for second-grade materials and higher order inertia. Int. J. Solids Struct. 49, 2121–2137] are improved and extended.     

A thermodynamic approach to the rheology of highly interactive filler-polymer mixtures: Part I – Theory

The rheological behavior of highly interactive filler-polymer mixtures is simulated utilizing a double network created by the entangled polymer matrix and the adsorbed polymer. Both networks are represented by a nonlinear viscoelastic constitutive equation. The dependence of rheological properties on filler concentration is taken into account through the bridging density resulting from polymer-filler interactions and a hydrodynamic reinforcement. The relative contribution of both networks is computed through the energy balance consistent with the thermodynamics of the polymer-filler chemical interactions and fluid mechanics. This self-consistent approach allows one to calculate the strain dependence of dynamic properties under oscillatory flow and shear rate dependence of stresses under steady simple shear flow and upon start up and cessation of shear flow. Received: 11 May 2000 Accepted: 8 March 2001     

Elastic–plastic stresses in a rotating disc of transversely isotropic material fitted with a shaft and subjected to thermal gradient

The elastic–plastic stresses in a rotating disc of transversely isotropic material fitted with a shaft and subjected to thermal gradient has been analyzed by using Seth’s transition theory and generalized strain measure. It has been observed that disc made of beryl and magnesium materials requires higher angular speed to yield at the inner surface in comparison to the disc made of brass material. The radial stress has a maximum at the internal surface of the disc made of beryl, magnesium and brass materials, but circumferential stress neither maximum nor minimum at this surface. With the introduction of thermal effect, the value of circumferential stress has a maximum at the external surface of the disc made of the beryl and magnesium, but the reverse results are obtained for the disc made of brass material. The combined impacts of temperature and angular speed have been displayed numerically and depicted graphically.