Numerical evaluation of material model and thickness effect on LSP induced residual stresses

Laser Shock Peening (LSP) is a technological process used to improve mechanical properties in metallic components. When a short and intense laser pulse irradiates a metallic surface, high pressure plasma is generated on the treated surface; elasto-plastic waves, then, propagate inside the target and create plastic strain. This surface treatment induces a deep compressive residual stresses field on the treated area and through the thickness; such compressive residual stress is expected to increase the fatigue resistance, and reduce the detrimental effects of corrosion and stress corrosion cracking.

Evaluation of LS-DYNA wood material model 143 /

"HRDI-04/08-05(1M)E"--P. [4] of cover.

Multi-objective optimization of material model parameters of an adhesive layer by using SPEA2

The usage of multi material structures in industry, especially in the automotive industry are increasing. To overcome the difficulties in joining these structures, adhesives have several benefits over traditional joining methods. Therefore, accurate simulations of the entire process of fracture including the adhesive layer is crucial. In this paper, material parameters of a previously developed meso mechanical finite element (FE) model of a thin adhesive layer are optimized using the Strength Pareto Evolutionary Algorithm (SPEA2). Objective functions are defined as the error between experimental data and simulation data. The experimental data is provided by previously performed experiments where an adhesive layer was loaded in monotonically increasing peel and shear. Two objective functions are dependent on 9 model parameters (decision variables) in total and are evaluated by running two FEsimulations, one is loading the adhesive layer in peel and the other in shear. The original study converted the two objective functions into one function that resulted in one optimal solution. In this study, however, a Pareto frontis obtained by employing the SPEA2 algorithm. Thus, more insight into the material model, objective functions, optimal solutions and decision space is acquired using the Pareto front. We compare the results and show good agreement with the experimental data.

Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

Defining of the material models for stainless steels to be used in FE modeling

The aim of this work is to study the results of tensile tests for austenitic stainless steel type 304 and make accurate FE-models according to the results of the tests. Tensile tests were made at Central Research Institute of Structural Material, Prometey at Saint Petersburg and Mariyenburg in Russia. The test specimens for the tensile tests were produced at Lappeenranta University of Technology in a Laboratory of Steel Structures. In total 4 different tests were made, two with base material specimens and two with transverse butt weld specimens. Each kind of a specimen was tested at room temperature and at low temperature. By comparing the results of room and low temperature tests of similar test specimen we get to study the results of work hardening that affect the austenitic steels at below room temperature. The produced specimens are to be modeled accurately and then imported for nonlinear FEM- analyzing. Using the data gained from the tensile tests the aim is to get the models work like the specimens did during the tests. By using the analyzed results of the FE-models the aim is to calculate and get the stress-strain curves that correspond to the results acquired from the tensile tests.

Hygro-elasto-plastic behavior of planar orthotropic material

The mechanical and hygroscopic properties of paper and board are factors affecting the whole lifecycle of a product, including paper/board quality, production, converting, and material and energy savings. The progress of shrinkage profiles, loose edges of web, baggy web causing wrinkling and misregistration in printing are examples of factors affecting runnability and end product quality in the drying section and converting processes, where paper or board is treated as a moving web. The structural properties and internal stresses or plastic strain differences built up during production also cause the end-product defects related to distortion of the shape of the product such as sheet or box. The objective of this work was to construct a model capable of capturing the characteristic behavior of hygroscopic orthotropic material under moisture change, during different external in-plane stretch or stress conditions. Two independent experimental models were constructed: the elasto-plastic material model and the hygroexpansivity-shrinkage model. Both describe the structural properties of the sheet with a fiber orientation probability distribution, and both are functions of the dry solids content and fiber orientation anisotropy index. The anisotropy index, introduced in this work, simplifies the procedure of determining the constitutive parameters of the material model and the hygroexpansion coefficients in different in-plane directions of the orthotropic sheet. The mathematically consistent elasto-plastic material model and the dry solids content dependent hygroexpansivity have been constructed over the entire range from wet to dry. The presented elastoplastic and hygroexpansivity-shrinkage models can be used in an analytical approach to estimate the plastic strain and shrinkage in simple one-dimensional cases. For studies of the combined and more complicated effects of hygro-elasto-plastic behavior, both models were implemented in a finite element program for a numerical solution. The finite element approach also offered possibilities for studying different structural variations of orthotropic planar material, as well as local buckling behavior and internal stress situations of the sheet or web generated by local strain differences. A comparison of the simulation examples presented in this work to results published earlier confirms that the hygro-elasto-plastic model provides at least qualitatively reasonable estimates. The application potential of the hygro-elasto-plastic model is versatile, including several phenomena and defects appearing in the drying, converting and end-use conditions of the paper or board webs and products, or in other corresponding complex planar materials.

A new damage model for composite laminates

Aircraft composite structures must have high stiffness and strength with low weight, which can guarantee the increase of the pay-load for airplanes without losing airworthiness. However, the mechanical behavior of composite laminates is very complex due the inherent anisotropy and heterogeneity. Many researchers have developed different failure progressive analyses and damage models in order to predict the complex failure mechanisms. This work presents a damage model and progressive failure analysis that requires simple experimental tests and that achieves good accuracy. Firstly, the paper explains damage initiation and propagation criteria and a procedure to identify the material parameters. In the second stage, the model was implemented as a UMAT (User Material Subroutine), which is linked to finite element software, ABAQUS (TM), in order to predict the composite structures behavior. Afterwards, some case studies, mainly off-axis coupons under tensile or compression loads, with different types of stacking sequence were analyzed using the proposed material model. Finally, the computational results were compared to the experimental results, verifying the capability of the damage model in order to predict the composite structure behavior. (C) 2011 Elsevier Ltd. All rights reserved.

Designing piezoresistive plate-based sensors with distribution of piezoresistive material using topology optimization

Piezoresistive sensors are commonly made of a piezoresistive membrane attached to a flexible substrate, a plate. They have been widely studied and used in several applications. It has been found that the size, position and geometry of the piezoresistive membrane may affect the performance of the sensors. Based on this remark, in this work, a topology optimization methodology for the design of piezoresistive plate-based sensors, for which both the piezoresistive membrane and the flexible substrate disposition can be optimized, is evaluated. Perfect coupling conditions between the substrate and the membrane based on the `layerwise' theory for laminated plates, and a material model for the piezoresistive membrane based on the solid isotropic material with penalization model, are employed. The design goal is to obtain the configuration of material that maximizes the sensor sensitivity to external loading, as well as the stiffness of the sensor to particular loads, which depend on the case (application) studied. The proposed approach is evaluated by studying two distinct examples: the optimization of an atomic force microscope probe and a pressure sensor. The results suggest that the performance of the sensors can be improved by using the proposed approach.

The Cohesive Crack Model Applied To Notched PMMA Specimens Obeying a Non Linear Behaviour Under Torsion Loading

This article presents a new material model developed with the aim of analyzing failure of blunt notched components made of nonlinear brittle materials. The model, which combines the cohesive crack model with Hencky's theory of total deformations, is used to simulate an experimental benchmark carried out previously by the authors. Such combination is achieved through the embedded crack approach concept. In spite of the unavailability of precise material data, the numerical predictions obtained show good agreement with the experimental results.

A coupled elastoplastic-damage constitutive model with Lode angle dependent failure criterion

A coupled elastoplastic-damage constitutive model with Lode angle dependent failure criterion for high strain and ballistic applications is presented. A Lode angle dependent function is added to the equivalent plastic strain to failure definition of the Johnson–Cook failure criterion. The weakening in the elastic law and in the Johnson–Cook-like constitutive relation implicitly introduces the Lode angle dependency in the elastoplastic behaviour. The material model is calibrated for precipitation hardened Inconel 718 nickel-base superalloy. The combination of a Lode angle dependent failure criterion with weakened constitutive equations is proven to predict fracture patterns of the mechanical tests performed and provide reliable results. Additionally, the mesh size dependency on the prediction of the fracture patterns was studied, showing that was crucial to predict such patterns

Root cause analysis applied to a finite element model's refinement of a negative stiffness structure

Negative Stiffness Structures are mechanical systems that require a decrease in the applied force to generate an increase in displacement. They are structures that possess special characteristics such as snap-through and bi-stability. All of these features make them particularly suitable for different applications, such as shock-absorption, vibration isolation and damping. From this point of view, they have risen awareness of their characteristics and, in order to match them to the application needed, a numerical simulation is of great interest. In this regard, this thesis is a continuation of previous studies in a circular negative stiffness structure and aims at refine the numerical model by presenting a new solution. To that end, an investigation procedure is needed. Amongst all of the methods available, root cause analysis was the chosen one to perform the investigation since it provides a clear view of the problem under analysis and a categorization of all the causes behind it. As a result of the cause-effect analysis, the main causes that have influence on the numerical results were obtained. Once all of the causes were listed, solutions to them were proposed and it led to a new numerical model. The numerical model proposed was of nonlinear type of analysis with hexagonal elements and a hyperelastic material model. The results were analyzed through force-displacement curves, allowing for the visualization of the structure’s energy recovery. When compared to the results obtained from the experimental part, it is evident that the trend is similar and the negative stiffness behaviour is present.

Mechanical characterization of aortas using 2D ultrasound elastography

Rupture of aortic aneurysms (AA) is a major cause of death in the Western world. Currently, clinical decision upon surgical intervention is based on the diameter of the aneurysm. However, this method is not fully adequate. Noninvasive assessment of the elastic properties of the arterial wall can be a better predictor for AA growth and rupture risk. The purpose of this study is to estimate mechanical properties of the aortic wall using in vitro inflation testing and 2D ultrasound (US) elastography, and investigate the performance of the proposed methodology for physiological conditions. Two different inflation experiments were performed on twelve porcine aortas: 1) a static experiment for a large pressure range (0 – 140 mmHg); 2) a dynamic experiment closely mimicking the in vivo hemodynamics at physiological pressures (70 – 130 mmHg). 2D raw radiofrequency (RF) US datasets were acquired for one longitudinal and two cross-sectional imaging planes, for both experiments. The RF-data were manually segmented and a 2D vessel wall displacement tracking algorithm was applied to obtain the aortic diameter–time behavior. The shear modulus G was estimated assuming a Neo-Hookean material model. In addition, an incremental study based on the static data was performed to: 1) investigate the changes in G for increasing mean arterial pressure (MAP), for a certain pressure difference (30, 40, 50 and 60 mmHg); 2) compare the results with those from the dynamic experiment, for the same pressure range. The resulting shear modulus G was 94 ± 16 kPa for the static experiment, which is in agreement with literature. A linear dependency on MAP was found for G, yet the effect of the pressure difference was negligible. The dynamic data revealed a G of 250 ± 20 kPa. For the same pressure range, the incremental shear modulus (Ginc) was 240 ± 39 kPa, which is in agreement with the former. In general, for all experiments, no significant differences in the values of G were found between different image planes. This study shows that 2D US elastography of aortas during inflation testing is feasible under controlled and physiological circumstances. In future studies, the in vivo, dynamic experiment should be repeated for a range of MAPs and pathological vessels should be examined. Furthermore, the use of more complex material models needs to be considered to describe the non-linear behavior of the vascular tissue.

Numerical Assessment of the out-of-plane response of a brick masonry structure without box behaviour

This paper presents the assessment of the out-of-plane response due to seismic loading of a masonry structure without rigid diaphragm. This structure corresponds to real scale brick masonry specimen with a main façade connected to two return walls. Two modelling approaches were defined for this evaluation. The first one consisted on macro modelling, whereas the second one on simplified micro modelling. As a first step of this study, static nonlinear analyses were conducted to the macro model aiming at evaluating the out-of-plane response and failure mechanism of the masonry structure. A sensibility analyses was performed in order to assess the mesh size and material model dependency. In addition, the macro models were subjected to dynamic nonlinear analyses with time integration in order to assess the collapse mechanism. Finally, these analyses were also applied to a simplified micro model of the masonry structure. Furthermore, these results were compared to experimental response from shaking table tests. It was observed that these numerical techniques simulate correctly the in-plane behaviour of masonry structures. However, the

Simulation of the deformation behaviour of two-phase composites

Two-phase polycrystal, elasto-viscoplastic material model, heterogeneity, anisotropy, flow behaviour, crystallographic texture, strain field

Ab initio absorption spectrum of NiO combining molecular dynamics with the embedded cluster approach in a discrete reaction field

We developed a procedure that combines three complementary computational methodologies to improve the theoretical description of the electronic structure of nickel oxide. The starting point is a Car-Parrinello molecular dynamics simulation to incorporate vibrorotational degrees of freedom into the material model. By means ofcomplete active space self-consistent field second-order perturbation theory (CASPT2) calculations on embedded clusters extracted from the resulting trajectory, we describe localized spectroscopic phenomena on NiO with an efficient treatment of electron correlation. The inclusion of thermal motion into the theoretical description allowsus to study electronic transitions that, otherwise, would be dipole forbidden in the ideal structure and results in a natural reproduction of the band broadening. Moreover, we improved the embedded cluster model by incorporating self-consistently at the complete active space self-consistent field (CASSCF) level a discrete (or direct) reaction field (DRF) in the cluster surroundings. The DRF approach offers an efficient treatment ofelectric response effects of the crystalline embedding to the electronic transitions localized in the cluster. We offer accurate theoretical estimates of the absorption spectrum and the density of states around the Fermi level of NiO, and a comprehensive explanation of the source of the broadening and the relaxation of the charge transferstates due to the adaptation of the environment