An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 1: Rods

A fully conserving algorithm is developed in this paper for the integration of the equations of motion in nonlinear rod dynamics. The starting point is a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, which results in an extremely simple update of the rotational variables. The weak form is constructed with a non-orthogonal projection corresponding to the application of the virtual power theorem. Together with an appropriate time-collocation, it ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that nonlinear hyperelastic materials (and not only materials with quadratic potentials) are permitted without any prejudice on the conservation properties. Spatial discretization is performed via the finite element method and the performance of the scheme is assessed by means of several numerical simulations.

An EFG method for the nonlinear analysis of plates undergoing arbitrarily large deformations

The applicability of a meshfree approximation method, namely the EFG method, on fully geometrically exact analysis of plates is investigated. Based on a unified nonlinear theory of plates, which allows for arbitrarily large rotations and displacements, a Galerkin approximation via MLS functions is settled. A hybrid method of analysis is proposed, where the solution is obtained by the independent approximation of the generalized internal displacement fields and the generalized boundary tractions. A consistent linearization procedure is performed, resulting in a semi-definite generalized tangent stiffness matrix which, for hyperelastic materials and conservative loadings, is always symmetric (even for configurations far from the generalized equilibrium trajectory). Besides the total Lagrangian formulation, an updated version is also presented, which enables the treatment of rotations beyond the parameterization limit. An extension of the arc-length method that includes the generalized domain displacement fields, the generalized boundary tractions and the load parameter in the constraint equation of the hyper-ellipsis is proposed to solve the resulting nonlinear problem. Extending the hybrid-displacement formulation, a multi-region decomposition is proposed to handle complex geometries. A criterium for the classification of the equilibrium`s stability, based on the Bordered-Hessian matrix analysis, is suggested. Several numerical examples are presented, illustrating the effectiveness of the method. Differently from the standard finite element methods (FEM), the resulting solutions are (arbitrary) smooth generalized displacement and stress fields. (c) 2007 Elsevier Ltd. All rights reserved.

Effects of elastic indenter deformation on spherical instrumented indentation tests: the reduced elastic modulus

Although the Hertz theory is not applicable in the analysis of the indentation of elastic-plastic materials, it is common practice to incorporate the concept of indenter/specimen combined modulus to consider indenter deformation. The appropriateness was assessed of the use of reduced modulus to incorporate the effect of indenter deformation in the analysis of the indentation with spherical indenters. The analysis based on finite element simulations considered four values of the ratio of the indented material elastic modulus to that of the diamond indenter, E/E(i) (0, 0.04, 0.19, 0.39), four values of the ratio of the elastic reduced modulus to the initial yield strength, E(r)/Y (0, 10, 20, 100), and two values of the ratio of the indenter radius to maximum total displacement, R/delta(max) (3, 10). Indenter deformation effects are better accounted for by the reduced modulus if the indented material behaves entirely elastically. In this case, identical load-displacement (P - delta) curves are obtained with rigid and elastic spherical indenters for the same elastic reduced modulus. Changes in the ratio E/E(i), from 0 to 0.39, resulted in variations lower than 5% for the load dimensionless functions, lower than 3% in the contact area, A(c), and lower than 5% in the ratio H/E(r). However, deformations of the elastic indenter made the actual radius of contact change, even in the indentation of elastic materials. Even though the load dimensionless functions showed only a little increase with the ratio E/E(i), the hardening coefficient and the yield strength could be slightly overestimated when algorithms based on rigid indenters are used. For the unloading curves, the ratio delta(e)/delta(max), where delta(e) is the point corresponding to zero load of a straight line with slope S from the point (P(max), delta(max)), varied less than 5% with the ratio E/E(i). Similarly, the relationship between reduced modulus and the unloading indentation curve, expressed by Sneddon`s equation, did not reveal the necessity of correction with the ratio E/E(i). The most affected parameter in the indentation curve, as a consequence of the indentation deformation, was the ratio between the residual indentation depth after complete unloading and the maximum indenter displacement, delta(r)/delta(max) (up to 26%), but this variation did not significantly decrease the capability to estimate hardness and elastic modulus based on the ratio of the residual indentation depth to maximum indentation depth, h(r)/h(max). In general, the results confirm the convenience of the use of the reduced modulus in the spherical instrumented indentation tests.

A critical reassessment of elastic unloading in sharp instrumented indentation experiments and the extraction of mechanical properties

This work examines the extraction of mechanical properties from instrumented indentation P-h(s) curves via extensive three-dimensional finite element analyses for pyramidal tips in a wide range of solids under frictional and frictionless contact conditions. Since the topography of the imprint changes with the level of pile-up or sink-in, a relationship is identified between correction factor beta in the elastic equation for the unloading indentation stage and the amount of surface deformation effects. It is shown that the presumption of a constant beta significantly affects mechanical property extractions. Consequently, a new best-fit function is found for the correlation between penetration depth ratios h(e)/h(max), h(r)/h(max) and n, circumventing the need for the assumption of a constant value for beta, made in our prior investigation [Acta Mater. 53 (2005) pp. 3545-3561]. Simulations under frictional contact conditions provide sensible boundaries for the influence of friction on both h(e)/h(max) and h(r)/h(max). Friction is essentially found to induce an overestimation in the inferred n. Instrumented indentation experiments are also performed in three archetypal metallic materials exhibiting distinctly different contact responses. Mechanical property extractions are finally demonstrated in each of these materials.

Influence of hardness of the harder body on wear regime transition in a sliding pair of steels

In the unlubricated sliding wear of steels the mild-severe and severe-mild wear transitions have long been investigated. The effect of system inputs such as normal load, sliding speed, environment humidity and temperature, material properties, among others, on those transitions have also been studied. Although transitions seem to be caused by microstructural changes, surfaces oxidation and work-hardening, some questions remain regarding the way each aspect is involved. Since the early studies in sliding wear, it has usually been assumed that only the material properties of the softer body influence the wear behavior of contacting surfaces. For example, the Archard equation involves only the hardness of the softer body, without considering the hardness of the harder body. This work aims to discuss the importance of the harder body hardness in determining the wear regime operation. For this, pin-on-disk wear tests were carried out, in which the disk material was always harder than the pin material. Variations of the friction force and vertical displacement of the pin were registered during the tests. A material characterization before and after tests was conducted using stereoscopy and scanning electron microscopy (SEM) methods, in addition to mass loss, surface roughness and microhardness measurements. The wear results confirmed the occurrence of a mild-severe wear transition when the disk hardness was decreased. The disk hardness to pin hardness ratio (H(d)/H(p)) was used as a criterion to establish the nature of surface contact deformation and to determine the wear regime transition. A predominantly elastic or plastic contact, characterized by H(d)/H(p) values higher or lower than one, results in a mild or severe wear regime operation, respectively. (c) 2009 Elsevier B.V. All rights reserved.

Friction behavior of lubricated zinc phosphate coatings

The aim of the present work is to elucidate the influence of lubricants on the friction behavior of zinc phosphated coatings and provide an explanation for the results in terms of physical-chemical interactions between lubricant and phosphate. The friction behavior was studied through a sliding wear test, with a conventional ball-on-disc configuration. Discs, made of AISI 1006 low carbon steel. uncoated and coated with zinc phosphate, were tested against bearing steel balls. A stearate sodium soap, paraffinic oil and both soap and oil were used as lubricants. The sodium stearate soap was found to have the best seizure resistance. The nature of the interfacial forces between the lubricant and surface has an important role in determining the friction behavior. (C) 2008 Elsevier B.V. All rights reserved.

Axial wedge effect in tilted hydrodynamic journal bearings

Hydrodynamic journal bearings are susceptible to static angular misalignment, resulting from improper assemblage, elastic and thermal distortion of the shaft and bearing housing, and also manufacturing errors. Several previous works on the theme, both theoretical and experimental, focused on the determination of the static properties of angular misaligned bearings. Although some reports show agreement between theoretical and experimental results, the increasingly severe operating conditions of hydrodynamic bearings (heavy loads and high rotational speeds) require more reliable theoretical formulations for the evaluation of the journal performance during the design process. The consideration of the angular misalignment in the derivation of the Reynolds equation is presented here in detail, showing that properly conducted geometric and magnitude-order analyses lead to the inclusion of an axial wedge effect term that influences the velocity and pressure fields in the lubricant film. Numerical results evidence that this axial wedge effect more significantly affects the hydrodynamic forces and static operational properties of tilted short journal bearings.

Identification of Elastic, Dielectric, and Piezoelectric Constants in Piezoceramic Disks

Three-dimensional modeling of piezoelectric devices requires a precise knowledge of piezoelectric material parameters. The commonly used piezoelectric materials belong to the 6mm symmetry class, which have ten independent constants. In this work, a methodology to obtain precise material constants over a wide frequency band through finite element analysis of a piezoceramic disk is presented. Given an experimental electrical impedance curve and a first estimate for the piezoelectric material properties, the objective is to find the material properties that minimize the difference between the electrical impedance calculated by the finite element method and that obtained experimentally by an electrical impedance analyzer. The methodology consists of four basic steps: experimental measurement, identification of vibration modes and their sensitivity to material constants, a preliminary identification algorithm, and final refinement of the material constants using an optimization algorithm. The application of the methodology is exemplified using a hard lead zirconate titanate piezoceramic. The same methodology is applied to a soft piezoceramic. The errors in the identification of each parameter are statistically estimated in both cases, and are less than 0.6% for elastic constants, and less than 6.3% for dielectric and piezoelectric constants.

Experiments on stress-triaxiality dependence of material behavior of aluminum alloys

The paper presents and discusses experimental procedures, visual observations and test results considered important to obtain data that can be used in validation of constitutive relations and failure criteria. The aim is to investigate the combined effects of stress intensity, stress-triaxiality and Lode parameter on the material response and failure behavior of aluminum alloys. Smooth and pre-notched tensile and shear specimens were manufactured from both very thin sheets and thicker plates to cover a wide range of stress triaxialities and Lode parameters. In addition, modified Arcan specimens were designed allowing investigation of the effect of sudden changes in stress states and deformation modes on the material behavior. (C) 2009 Elsevier Ltd. All rights reserved.

A ductile damage criterion at various stress triaxialities

The paper discusses the effect of stress triaxiality on the onset and evolution of damage in ductile metals. A series of tests including shear tests and experiments oil smooth and pre-notched tension specimens wits carried Out for it wide range of stress triaxialities. The underlying continuum damage model is based oil kinematic definition of damage tensors. The modular structure of the approach is accomplished by the decomposition of strain rates into elastic, plastic and damage parts. Free energy functions with respect to fictitious undamaged configurations as well as damaged ones are introduced separately leading to elastic material laws which are affected by increasing damage. In addition, a macroscopic yield condition and a flow rule are used to adequately describe the plastic behavior. Numerical simulations of the experiments are performed and good correlation of tests and numerical results is achieved. Based oil experimental and numerical data the damage criterion formulated in stress space is quantified. Different branches of this function are taken into account corresponding to different damage modes depending oil stress triaxiality and Lode parameter. In addition, identification of material parameters is discussed ill detail. (C) 2007 Elsevier Ltd. All rights reserved.

Effect of alpha prime due to 475 degrees C aging on fracture behavior and corrosion resistance of DIN 1.4575 and MA 956 high performance ferritic stainless steels

The 475 degrees C embrittlement in stainless steels is a well-known phenomenon associated to alpha prime (alpha`) formed by precipitation or spinodal decomposition. Many doubts still remain on the mechanism of alpha` formation and its consequence on deformation and fracture mechanisms and corrosion resistance. In this investigation, the fracture behavior and corrosion resistance of two high performance ferritic stainless steels were investigated: a superferritic DIN 1.4575 and MA 956 superalloy were evaluated. Samples of both stainless steels (SS) were aged at 475 degrees C for periods varying from 1 to 1,080 h. Their fracture surfaces were observed using scanning electron microscopy (SEM) and the cleavage planes were determined by electron backscattering diffraction (EBSD). Some samples were tested for corrosion resistance using electrochemical impedance spectroscopy (EIS) and potentiodynamic polarization. Brittle and ductile fractures were observed in both ferritic stainless steels after aging at 475 degrees C. For aging periods longer than 500 h, the ductile fracture regions completely disappeared. The cleavage plane in the DIN 1.4575 samples aged at 475 degrees C for 1,080 h was mainly {110}, however the {102}, {314}, and {131} families of planes were also detected. The pitting corrosion resistance decreased with aging at 475 degrees C. The effect of alpha prime on the corrosion resistance was more significant in the DIN 1.4575 SS comparatively to the Incoloy MA 956.

Annealing behavior of ferritic-martensitic 9%Cr-ODS-Eurofer steel

Oxide dispersion strengthened ferritic-martensitic steels are potential candidates for applications in future fusion power plants. High creep resistance, good oxidation resistance, reduced neutron activation and microstructural long-term stability at temperatures of about 650-700 degrees C are required in this context. In order to evaluate its thermal stability in the ferritic phase field, samples of the reduced activation ferritic-martensitic 9%Cr-ODS-Eurofer steel were cold rolled to 50% and 80% reductions and further annealed in vacuum from 300 to 800 degrees C for 1 h. The characterization in the annealed state was performed by scanning electron microscopy in the backscattered electron mode, high-resolution electron backscatter diffraction and transmission electron microscopy. Results show that the fine dispersion of Y-based particles (about 10 nm in size) is effective to prevent recrystallization. The low recrystallized volume fraction (<0.1) is associated to the nuclei found at prior grain boundaries and around large M(23)C(6) particles. Static recovery was found to be the predominant softening mechanism of this steel in the investigated temperature range. (c) 2010 Elsevier B.V. All rights reserved.

Influence of Water Content, Time, and Temperature on the Rheological Behavior of Polyethylene Terephtalate

In this work, the main factors affecting the rheological behavior of polyethylene terephtalate (PET) in the linear viscoelastic regime (water content, time delay before test, duration of experiment, and temperature) were accessed. Small amplitude oscillatory shear tests were performed after different time delays ranging from 300 to 5000 s for samples with water contents ranging from 0.02 to 0.45 wt %. Time sweep tests were carried out for different durations to explain the changes undergone by PET before and during small amplitude oscillatory shear measurements. Immediately after the time sweep tests, the PET samples were removed from the rheometer, analyzed by differential scanning calorimetry and their molar mass was obtained by viscometry analysis. It was shown that for all the samples, the delay before test and residence time within the rheometer (i.e. duration of experiment) result in structural changes of the PET samples, such as increase or decrease of molar mass, broadening of molar mass distribution, and branching phenomena. (C) 2010 Wiley Periodicals, Inc. J Appl Polym Sci 116: 3525-3533, 2010

Modeling of Plastic Deformation Effects in Ferromagnetic Thin Films

To explain the magnetic behavior of plastic deformation of thin magnetic films (Fe and permalloy) on an elastic substrate (nitinol), it is noted that unlike in the bulk, the dislocation density does not increase dramatically because of the dimensional constraint. As a result, the resulting residual stress, even though strain hardening is limited, dominates the observed magnetic behavior. Thus, with the field parallel to the stress axis, the compressive residual stress resulting from plastic deformation causes a decrease in remanence and an increase in coercivity; and with the field perpendicular to the stress axis, the resulting compressive residual stress causes an increase in remanence and a decrease in coercivity. These elements have been inserted into the model previously developed for plastic deformation in the bulk, producing the aforementioned behavior, which has been observed experimentally in the films.

High temperature flexural strength and fracture toughness of AlN with Y(2)O(3) ceramic

The effects of temperature on the fast fracture behavior of aluminum nitride with 5 wt% Y(2)O(3) ceramic were investigated. Four-point flexural strength and fracture toughness were measured in air at several temperatures (30-1,300 A degrees C). The flexural strength gradually decreased with the increase of temperature up to 1,000 A degrees C due to the change in the fracture mode from transgranular to intergranular, and then became almost constant up to 1,300 A degrees C. Two main flaw types as fracture origin were identified: small surface flaw and large pores. The volume fraction of the large pores was only 0.01%; however, they limited the strength on about 50% of the specimens. The fracture toughness decreased slightly up to 800 A degrees C controlled by the elastic modulus change, and then decreased significantly at 1,000 A degrees C due to the decrease in the grain-boundary toughness. Above 1,000 A degrees C, the fracture toughness increased significantly, and at 1,300 A degrees C, its value was close to that measured at room temperature.