A hybrid-mixed finite element formulation for the geometrically exact analysis of three-dimensional framed structures

This paper addresses the development of a hybrid-mixed finite element formulation for the quasi-static geometrically exact analysis of three-dimensional framed structures with linear elastic behavior. The formulation is based on a modified principle of stationary total complementary energy, involving, as independent variables, the generalized vectors of stress-resultants and displacements and, in addition, a set of Lagrange multipliers defined on the element boundaries. The finite element discretization scheme adopted within the framework of the proposed formulation leads to numerical solutions that strongly satisfy the equilibrium differential equations in the elements, as well as the equilibrium boundary conditions. This formulation consists, therefore, in a true equilibrium formulation for large displacements and rotations in space. Furthermore, this formulation is objective, as it ensures invariance of the strain measures under superposed rigid body rotations, and is not affected by the so-called shear-locking phenomenon. Also, the proposed formulation produces numerical solutions which are independent of the path of deformation. To validate and assess the accuracy of the proposed formulation, some benchmark problems are analyzed and their solutions compared with those obtained using the standard two-node displacement/ rotation-based formulation.

On the finite element modeling of fatigue crack growth in pressurized cylindrical shells

A methodology for the computational modeling of the fatigue crack growth in pressurized shell structures, based on the finite element method and concepts of Linear Elastic Fracture Mechanics, is presented. This methodology is based on that developed by Potyondy [Potyondy D, Wawrzynek PA, Ingraffea, AR. Discrete crack growth analysis methodology for through crack in pressurized fuselage structures. Int J Numer Methods Eng 1995;38:1633-1644], which consists of using four stress intensity factors, computed from the modified crack integral method, to predict the fatigue propagation life as well as the crack trajectory, which is computed as part of the numerical simulation. Some issues not presented in the study of Potyondy are investigated herein such as the influence of the crack increment size and the number of nodes per element (4 or 9 nodes) on the simulation results by means of a fatigue crack propagation simulation of a Boeing 737 airplane fuselage. The results of this simulation are compared with experimental results and those obtained by Potyondy [1]. (C) 2008 Elsevier Ltd. All rights reserved.

Buckling and post-buckling of extensible rods revisited: A multiple-scale solution

An exact non-linear formulation of the equilibrium of elastic prismatic rods subjected to compression and planar bending is presented, electing as primary displacement variable the cross-section rotations and taking into account the axis extensibility. Such a formulation proves to be sufficiently general to encompass any boundary condition. The evaluation of critical loads for the five classical Euler buckling cases is pursued, allowing for the assessment of the axis extensibility effect. From the quantitative viewpoint, it is seen that such an influence is negligible for very slender bars, but it dramatically increases as the slenderness ratio decreases. From the qualitative viewpoint, its effect is that there are not infinite critical loads, as foreseen by the classical inextensible theory. The method of multiple (spatial) scales is used to survey the post-buckling regime for the five classical Euler buckling cases, with remarkable success, since very small deviations were observed with respect to results obtained via numerical integration of the exact equation of equilibrium, even when loads much higher than the critical ones were considered. Although known beforehand that such classical Euler buckling cases are imperfection insensitive, the effect of load offsets were also looked at, thus showing that the formulation is sufficiently general to accommodate this sort of analysis. (c) 2008 Elsevier Ltd. All rights reserved.

Finite element and dimensional analysis algorithm for the prediction of mechanical properties of bulk materials and thin films

In this work, the applicability of a new algorithm for the estimation of mechanical properties from instrumented indentation data was studied for thin films. The applicability was analyzed with the aid of both three-dimensional finite element simulations and experimental indentation tests. The numerical approach allowed studying the effect of the substrate on the estimation of mechanical properties of the film, which was conducted based on the ratio h(max)/l between maximum indentation depth and film thickness. For the experimental analysis, indentation tests were conducted on AISI H13 tool steel specimens, plasma nitrated and coated with TiN thin films. Results have indicated that, for the conditions analyzed in this work, the elastic deformation of the substrate limited the extraction of mechanical properties of the film/substrate system. This limitation occurred even at low h(max)/l ratios and especially for the estimation of the values of yield strength and strain hardening exponent. At indentation depths lower than 4% of the film thickness, the proposed algorithm estimated the mechanical properties of the film with accuracy. Particularly for hardness, precise values were estimated at h(max)/l lower than 0.1, i.e. 10% of film thickness. (C) 2010 Published by Elsevier B.V.

3D CFD Simulation of Vortex-induced Vibration of Cylinder

The 3D flow around a circular cylinder free to oscillate transversely to the free stream was simulated using Computational Fluid Dynamics (CFD) and the Spalart-Allmaras Detached Eddy Simulation (DES) turbulence model for a Reynolds number Re = 10(4). Simulations were carried out for a small mass-damping parameter m*zeta = 0.00858, where m* = 3.3 and zeta = 0.0026. We found good agreement between the numerical results and experimental data. The simulations predicted the high observed amplitudes of the upper branch of vortex-induced vibrations for low mass-damping parameters.

Residual stresses in titanium nitride thin films obtained with step variation of substrate bias voltage during deposition

In this work, a series of depositions of titanium nitride (TiN) films on M2 and D2 steel substrates were conducted in a Triode Magnetron Sputtering chamber. The temperature; gas flow and pressure were kept constant during each run. The substrate bias was either decreased or increased in a sequence of steps. Residual stress measurements were later conducted through the grazing X-ray diffraction method. Different incident angles were used in order to change the penetration depth and to obtain values of residual stress at different film depths. A model described by Dolle was adapted as an attempt to calculate the values of residual stress at each incident angle as a function of the value from each individual layer. Stress results indicated that the decrease in bias voltage during the deposition has produced compressive residual stress gradients through the film thickness. On the other hand, much less pronounced gradients were found in one of the films deposited with increasing bias voltage. (C) 2010 Elsevier B.V. All rights reserved.

Experimental investigation of flow-induced vibration on isolated and tandem circular cylinders fitted with strakes

The effect of varying the geometric parameters of helical strakes on vortex-induced vibration (VIV) is investigated in this paper. The degree of oscillation attenuation or even suppression is analysed for isolated circular cylinder cases. How a cylinder fitted with strakes behaves when immersed in the wake of another cylinder in tandem arrangement is also investigated and these results are compared to those with a single straked cylinder. The experimental tests are conducted at a circulating water channel facility and the cylindrical models are mounted on a low-damping air bearing elastic base with one degree-of-freedom, restricted to oscillate in the transverse direction to the channel flow. Three strake pitches (p) and heights (h) are tested: p = 5, 10, 15d, and h = 0.1, 0.2, 0.25d. The mass ratio is 1.8 for all models. The Reynolds number range is from 1000 to 10000, and the reduced velocity varies up to 21. The cases with h = 0.1d strakes reduce the amplitude response when compared to the isolated plain cylinder, however the oscillation still persists. On the other hand, the cases with h = 0.2, 0.25d strakes almost completely suppress VIV. Spanwise vorticity fields, obtained through stereoscopic digital particle image velocimetry (SDPIV), show an alternating vortex wake for the p = 10d and h = 0.1d straked cylinder. The p = 10d and h = 0.2d cylinder wake has separated shear layers with constant width and no roll-up close to the body. The strakes do not increase the magnitude of the out-of-plane velocity compared to the isolated plain cylinder. However, they deflect the flow in the out-of-plane direction in a controlled way, which can prevent the vortex shedding correlation along the span. In order to investigate the wake interference effect on the strake efficiency, an experimental arrangement with two cylinders in tandem is employed. The centre-to-centre distance for the tandem arrangement varies from 2 to 6. When the downstream p = 10d and h = 0.2d cylinder is immersed in the wake of an upstream fixed plain cylinder, it loses its effectiveness compared with the isolated case. Although the oscillations have significant amplitude, they are limited, which is a different behaviour from that of a tandem configuration with two plain cylinders. For this particular case, the amplitude response monotonically increases for all gaps, except one, a trait usually found in galloping-like oscillations. SDPIV results for the tandem arrangements show alternating vortex shedding and oscillatory wake. (C) 2010 Elsevier Ltd. All rights reserved.

Residual stresses in TiN, DLC and MoS(2) coated surfaces with regard to their tribological fracture behaviour

Thin hard coatings on components and tools are used increasingly due to the rapid development in deposition techniques, tribological performance and application skills. The residual stresses in a coated surface are crucial for its tribological performance. Compressive residual stresses in PVD deposited TiN and DLC coatings were measured to be in the range of 0.03-4 GPa on steel substrate and 0.1-1.3 GPa on silicon. MoS(2) coatings had tensional stresses in the range of 0.8-1.3 on steel and 0.16 GPa compressive stresses on silicon. The fracture pattern of coatings deposited on steel substrate were analysed both in bend testing and scratch testing. A micro-scale finite element method (FEM) modelling and stress simulation of a 2 mu m TiN-coated steel surface was carried out and showed a reduction of the generated tensile buckling stresses in front of the sliding tip when compressive residual stresses of 1 GPa were included in the model. However, this reduction is not similarly observed in the scratch groove behind the tip, possibly due to sliding contact-induced stress relaxation. Scratch and bending tests allowed calculation of the fracture toughness of the three coated surfaces, based on both empirical crack pattern observations and FEM stress calculation, which resulted in highest values for TiN coating followed by MoS(2) and DLC coatings, being K(C) = 4-11, about 2, and 1-2 MPa M(1/2), respectively. Higher compressive residual stresses in the coating and higher elastic modulus of the coating correlated to increased fracture toughness of the coated surface. (C) 2009 Elsevier B.V. All rights reserved.

Energy dissipation in depth-sensing indentation as a characteristic of the nanoscratch behavior of coatings

Wear behavior of coatings has usually been described in terms of mechanical properties such as hardness (H) and effective elastic modulus (E*). Alternatively, an energy approach appears as a promising analysis taking into account the influence of those properties. In a nanoindentation test, the dissipated energy depends not only on the hardness and elastic modulus, but also on the elastic recovery (W(e)). This work aims to establish a relation between plastic deformation energy (E(p)) during depth-sensing indentation method and the grooving resistance of coatings in nanoscratch tests. An energy dissipation coefficient (K(d)) was defined, calculated as the ratio of the plastic to the total deformation energy (E(p)/E(t)), which represents the energy dissipation of materials. Reactive depositions using titanium as the target and nitrogen and methane as reactive gases were obtained by triode magnetron sputtering, in order to assess wear and nanoindentation data. A topographical, chemical and microstructural characterization has been conducted using X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), wave dispersion spectroscopy (WDS), scanning electron (SEM) and atomic force microscopy (AFM) techniques. Nanoscratch results showed that the groove depth was well correlated to the energy dissipation coefficient of the coatings. On the other hand, a reduction in the coefficient was found when the elastic recovery was increased. (C) 2009 Elsevier B.V. All rights reserved.

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.

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.

Ultrasonic material characterization using large-aperture PVDF receivers

This work describes the use of a large-aperture PVDF receiver in the measurement of liquid density and composite material elastic constants. The density measurement of several liquids is obtained with accuracy of 0.2% using a conventional NDE emitter transducer and a 70-mm-diameter, 52-mu m P(VDF-TrFE) membrane with gold electrodes. The determination of the elastic constants is based on the phase velocity measurement. Diffraction can lead to errors around 1% in velocity measurement when using alternatively the conventional pair of ultrasonic transducers (1-MHz frequency and 19-mm-diameter) operating in through-transmission mode, separated by a distance of 100 mm. This effect is negligible when using a pair of 10-MHz, 19-mm-diameter transducers. Nevertheless, the dispersion at 10 MHz can result in errors of about 0.5%, when measuring the velocity in composite materials. The use of an 80-mm diameter, 52-mu m-thick PVDF membrane receiver practically eliminates the diffraction effects in phase velocity measurement. The elastic constants of a carbon fiber reinforced polymer were determined and compared with the values obtained by a tensile test. (C) 2009 Elsevier B. V. All rights reserved.

Toward Optimal Design of Piezoelectric Transducers Based on Multifunctional and Smoothly Graded Hybrid Material Systems

This work explores the design of piezoelectric transducers based on functional material gradation, here named functionally graded piezoelectric transducer (FGPT). Depending on the applications, FGPTs must achieve several goals, which are essentially related to the transducer resonance frequency, vibration modes, and excitation strength at specific resonance frequencies. Several approaches can be used to achieve these goals; however, this work focuses on finding the optimal material gradation of FGPTs by means of topology optimization. Three objective functions are proposed: (i) to obtain the FGPT optimal material gradation for maximizing specified resonance frequencies; (ii) to design piezoelectric resonators, thus, the optimal material gradation is found for achieving desirable eigenvalues and eigenmodes; and (iii) to find the optimal material distribution of FGPTs, which maximizes specified excitation strength. To track the desirable vibration mode, a mode-tracking method utilizing the `modal assurance criterion` is applied. The continuous change of piezoelectric, dielectric, and elastic properties is achieved by using the graded finite element concept. The optimization algorithm is constructed based on sequential linear programming, and the concept of continuum approximation of material distribution. To illustrate the method, 2D FGPTs are designed for each objective function. In addition, the FGPT performance is compared with the non-FGPT one.

Multi-actuated functionally graded piezoelectric micro-tools design: A multiphysics topology optimization approach

Micro-tools offer significant promise in a wide range of applications Such as cell Manipulation, microsurgery, and micro/nanotechnology processes. Such special micro-tools consist of multi-flexible structures actuated by two or more piezoceramic devices that must generate output displacements and forces lit different specified points of the domain and at different directions. The micro-tool Structure acts as a mechanical transformer by amplifying and changing the direction of the piezoceramics Output displacements. The design of these micro-tools involves minimization of the coupling among movements generated by various piezoceramics. To obtain enhanced micro-tool performance, the concept of multifunctional and functionally graded materials is extended by, tailoring elastic and piezoelectric properties Of the piezoceramics while simultaneously optimizing the multi-flexible structural configuration using multiphysics topology optimization. The design process considers the influence of piezoceramic property gradation and also its polarization sign. The method is implemented considering continuum material distribution with special interpolation of fictitious densities in the design domain. As examples, designs of a single piezoactuator, an XY nano-positioner actuated by two graded piezoceramics, and a micro-gripper actuated by three graded piezoceramics are considered. The results show that material gradation plays an important role to improve actuator performance, which may also lead to optimal displacements and coupling ratios with reduced amount of piezoelectric material. The present examples are limited to two-dimensional models because many of the applications for Such micro-tools are planar devices. Copyright (c) 2008 John Wiley & Sons, Ltd.

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.