197 research outputs found
Investigation of a gradient enriched Gurson-Tvergaard model for porous strain hardening materials
Size effects in a strain hardening porous solid are investigated using the Gurson-Tvergaard (GT) model enriched by a constitutive length parameter, as proposed by Niordson and Tvergaard [C.F. Niordson, V. Tvergaard, A homogenised model for size effects in porous metals, J. Mech. Phys. Solids (2019)]. The results are compared with unit cell calculations of regularly distributed voids embedded in a strain gradient enhanced matrix material. The strain gradient plasticity theory proposed by Fleck and Willis [N.A. Fleck, J.R. Willis, A mathematical basis for strain gradient plasticity theory. Part II: tensorial plastic multiplier, J. Mech. Phys. Solids 57 (2009) 1045–1057], extended to finite strains, is adopted for the cell model, consistent with the gradient enriched Gurson model. The gradient model allows for a material length parameter to enter the constitutive framework for dimensional consistency, while the enriched GT model has the same length parameter introduced through prefactors of the usual and factors. The continuum model featuring size-dependent Tvergaard-constants is used to investigate a strain hardening material with the strain gradient plasticity enriched cell model as reference. The two models are compared for three triaxialities, three initial void volume fractions, and three hardening exponents. The enriched GT model captures the effect of elevated yield point and suppressed void growth with increasing length parameter for all the cases investigated. The agreement between the models is good until severe void distortion or plastic flow localisation between neighbouring voids. The response curves and void growth curves for the enriched GT model deviate from those of the cell model at high axial strains. Void shape plots, which are only available for the cell model, show that the length parameter influences the shape of the void which in turn has impact on the material response curves and the void evolution. This is not captured by the enriched GT model as the voids are accounted for solely through a volume fraction parameter
Numerical simulation of grain-size effects on creep crack growth by means of grain elements
The effect of grain size on creep crack growth is investigated by means of a numerical technique in which the actual crack growth process is simulated in a discrete manner by grain elements and grain boundary elements. The grain elements account for the creep deformation of individual grains, while grain boundary cavitation and sliding are accounted for by grain boundary elements between the grains. This grain-element technique allows for an independent study of multiple grain size effects: a (direct) size effect related to the specimen size/grain size ratio or an (indirect) effect related to the effect of grain size on nucleation rate and creep resistance. Preliminary numerical results are presented concerning the direct effect of grain size, which predict that the crack growth rate and brittleness increase with grain size.
Molecular dynamics simulations of phase transformations in niti bicrystals
The influence of grain boundaries and grain misorientation on the nucleation and growth of martensite in an equi-atomic nickeltitanium (NiTi) shape memory alloy (SMA) is investigated by performing molecular dynamics (MD) simulations on bicrystals with a modified embedded atom method (MEAM) interatomic potential. Stress-induced martensitic transformations are simulated in bicrystals with mixed grain boundaries and the behavior of the bicrystal is compared to that of individual single crystals. Here, a particular bicrystal with < 110 > and < 111 > oriented austenite grains is chosen as an example. Results indicate that the mixed grain boundary in the austenite bicrystal acts as a nucleation site for stress-induced martensitic transformation in the grains. The deformation behavior and the transformation strain of the bicrystal fall in between those of the two corresponding single crystals.Applied Mechanics(OLD) MSE-
Multi-material topology optimization of viscoelastically damped structures: 21-26 August, Montreal, Canada
The design of high performance instruments often involves the attenuation of poorly damped resonant modes. Current design methods typically rely on informed trial and error based modifications to improve dynamic performance. In this contribution, we present a multi-material topology optimization as an alternative, systematic methodology to design structures with optimized damping characteristics. A parametric, level set-based topology optimization is employed to simultaneously distribute structural and viscoelastic material to optimize the structure’s damping characteristics. To model the viscoelastic behavior a complex-valued material modulus is applied. The structural loss factor is determined from the complex-valued eigensolutions and its value is maximized during the optimization. We demonstrate the performance of the optimization by maximizing the damping of a cantilever beam.Computational Design and Mechanic
Auto-generation by interaction of weak eddies
For channel flow, we explore how the interaction of weak eddies produces additional eddies by means of auto-generation. This is done by DNS of two eddies with different initial strengths, initial sizes and initial stream-wise spacing between them. The numerical procedure followed is similar to Zhou et al[1]. The two eddies merge into a single stronger eddy when a larger upstream and a smaller downstream eddy are placed within a certain initial stream-wise separation distance. Subsequently, the resulting stronger eddy is observed to auto-generate new eddies. The non-merging cases with small initial stream wise separation also auto-generate. The auto-generation is characterized by a rapid lift-up of an initial eddy, which blocks the incoming flow and leads to shear- layer roll-up and formation of a new eddy. The same sequence of events is observed in a fully developed turbulent boundary layer[2].Fluid MechanicsMulti Phase System
Use of abrupt strain path change for determining subsequent yield surface::Experimental study with metal sheets
A basic idea for a method for determining the subsequent yield surface in the vicinity of a current loading point by using an abrupt strain path change has been proposed recently by Kuroda and Tvergaard (Acta mater., 1999, 47, 3879). The proposed method is applied to real experimental studies. In a biaxial tensile testing apparatus, a cruciform specimen is used, with the strains measured by a biaxial-strain gauge. Then, with the hydraulic pressure of two sets of opposing hydraulic cylinders servo-controlled independently, the testing apparatus can be used to prescribe an abrupt change of the strain path. Both a cold-rolled steel sheet and an aluminum alloy sheet are investigated. The differences between the yield surface shapes found by the strain path change procedure and the shapes found by probing the yield points from the elastic region are shown and discussed for different cases. (C) 2000 Acta Metallurgica Inc. Published by Elsevier Science Ltd. All rights reserved
A consistent use of the Gurson-Tvergaard-Needleman damage model for the R-curve calculation
The scope of the present work is to point out a consistent simulation procedure for the quasi-static fracture processes, starting from the micro-structural characteristics of the material. To this aim, a local nineparameters
Gurson-Tvergaard-Needleman (GTN) damage law has been used. The damage parameters depend on the micro-structural characteristics and must be calculated, measured or opportunely tuned. This can be done, as proposed by the author, by using an opportunely tuned GTN model for the representative volume element simulations, in order to enrich the original damage model by considering also the defect size distribution. Once determined all the material parameters, an MT fracture test has been simulated by a FE code, to calculate the R-curve in an aeronautical Al-based alloy. The simulation procedure produced results in a very good agreement with the experimental data
A consistent use of the Gurson-Tvergaard-Needleman damage model for the R-curve calculation
The scope of the present work is to point out a consistent simulation procedure for the quasi-static fracture processes, starting from the micro-structural characteristics of the material. To this aim, a local nine-parameters Gurson-Tvergaard-Needleman (GTN) damage law has been used. The damage parameters depend on the micro-structural characteristics and must be calculated, measured or opportunely tuned. This can be done, as proposed by the author, by using an opportunely tuned GTN model for the representative volume element simulations, in order to enrich the original damage model by considering also the defect size distribution. Once determined all the material parameters, an MT fracture test has been simulated by a FE code, to calculate the R-curve in an aeronautical Al-based alloy. The simulation procedure produced results in a very good agreement with the experimental data
Bending effects on instabilities of internally pressurized tubes modelled by a nonlocal membrane theory
Multiscale Modeling and Simulation of Damage by Void Nucleation and Growth
Voids are observed to be generated under sufficient loading in many materials, ranging from polymers and metals to biological tissues. The presence of these voids can have drastic implications at the macroscopic level including strong material softening and more incipient fracture. Developing tools to appropriately account for these effects is therefore very desirable.
This thesis is concerned with both, the appearance of voids (nucleation process) and the modeling and simulation of materials in the presence of voids. A particular nucleation mechanism based on vacancy aggregation in high purity metallic single crystals is analyzed. A multiscale model is developed in order to obtain an approximate value of the time required for vacancies to form sufficiently large clusters for further growth by plastic deformation. It is based on quantum mechanical results, kinetic Monte Carlo methods and continuum mechanics estimates calibrated with quasi-continuum results. The ultimate goal of these simulations is to determine the feasibility of this nucleation mechanism under shock loading conditions, where the temperature and tensions are high and vacancy diffusion is promoted.
On the other hand, the effective behavior of materials with pre-existent voids is analyzed within the general framework of continuum mechanics and is therefore applicable to any material. The overall properties of the heterogeneous material are obtained through a two-level characterization: a representative volume element consisting of a hollow sphere is used to describe the "microscopic" fields, and an equivalent homogeneous material is used for the "macroscopic" behavior. A variational formulation of this two-scale model is presented. It provides a consistent definition of the macro-variables under general loading conditions, extending the well-known static averaging results so as to include microdynamic effects under finite deformations. This variational framework also provides a suitable starting point for time discretization and consistent definitions within discrete time. The spatial boundary value problem resulting from this multiscale model is solved with a particular spherical shell element specially developed for this problem. The approximation space is based on spherical harmonics, which respects the symmetries of the porous material and allows the representation of the fields on the sphere with very few degrees of freedom. Numerical tools, such as the exact representation of the boundary conditions and an exact quadrature rule, are also provided. The resulting numerical model is verified extensively, demonstrating good convergence results, and its applicability is shown through several material point calculations and a full two-scale finite element implementation.</p
- …
