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

Following the approach developed for rods in Part 1 of this paper (Pimenta et al. in Comput. Mech. 42:715-732, 2008), this work presents a fully conserving algorithm for the integration of the equations of motion in nonlinear shell dynamics. We begin with a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, allowing for an extremely simple update of the rotational variables within the scheme. The weak form is constructed via non-orthogonal projection, the time-collocation of which ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that general hyperelastic materials (and not only materials with quadratic potentials) are permitted in a totally consistent way. Spatial discretization is performed using the finite element method and the robust performance of the scheme is demonstrated by means of numerical examples.

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.

Punching Strength of Flat Slabs with Unbraced Shear Reinforcement

This paper analyzes the punching strength of concrete flat slabs with shear reinforcement that does not embrace flexural reinforcement. This paper also reports the results of tests of slabs without shear reinforcement. Finally, this paper shows some comparisons of tests of similar slabs without shear reinforcement and slabs with different types of shear reinforcement. The obtained results show that the use of shear reinforcement elements without embracement in the flexural reinforcement improves the punching strength of reinforced concrete flat slabs.

Planning and operation of large-scale water distribution systems with preemptive priorities

This paper describes the development of an optimization model for the management and operation of a large-scale, multireservoir water supply distribution system with preemptive priorities. The model considers multiobjectives and hedging rules. During periods of drought, when water supply is insufficient to meet the planned demand, appropriate rationing factors are applied to reduce water supply. In this paper, a water distribution system is formulated as a network and solved by the GAMS modeling system for mathematical programming and optimization. A user-friendly interface is developed to facilitate the manipulation of data and to generate graphs and tables for decision makers. The optimization model and its interface form a decision support system (DSS), which can be used to configure a water distribution system to facilitate capacity expansion and reliability studies. Several examples are presented to demonstrate the utility and versatility of the developed DSS under different supply and demand scenarios, including applications to one of the largest water supply systems in the world, the Sao Paulo Metropolitan Area Water Supply Distribution System in Brazil.

Volumetric reconstruction of the mean flow around circular cylinders fitted with strakes

The volumetric reconstruction technique presented in this paper employs a two-camera stereoscopic particle image velocimetry (SPIV) system in order to reconstruct the mean flow behind a fixed cylinder fitted with helical strakes, which are commonly used to suppress vortex-induced vibrations (VIV). The technique is based on the measurement of velocity fields at equivalent adjacent planes that results in pseudo volumetric fields. The main advantage over proper volumetric techniques is the avoidance of additional equipment and complexity. The averaged velocity fields behind the straked cylinders and the geometrical periodicity of the three-start configuration are used to further simplify the reconstruction process. Two straked cylindrical models with the same pitch (p = 10d) and two different heights (h = 0.1 and 0.2d) are tested. The reconstructed flow shows that the strakes introduce in the wake flow a well-defined wavelength of one-third of the pitch. Measurements of hydrodynamic forces, fluctuating velocity, vortex formation length, and vortex shedding frequency show the interdependence of the wake parameters. The vortex formation length is increased by the strakes, which is an important effect for the suppression of vortex-induced vibrations. The results presented complement previous investigations concerning the effectiveness of strakes as VIV suppressors and provide a basis of comparison to numerical simulations.

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.

Dynamic Imaging in Electrical Impedance Tomography of the Human Chest With Online Transition Matrix Identification

One of the electrical impedance tomography objectives is to estimate the electrical resistivity distribution in a domain based only on electrical potential measurements at its boundary generated by an imposed electrical current distribution into the boundary. One of the methods used in dynamic estimation is the Kalman filter. In biomedical applications, the random walk model is frequently used as evolution model and, under this conditions, poor tracking ability of the extended Kalman filter (EKF) is achieved. An analytically developed evolution model is not feasible at this moment. The paper investigates the identification of the evolution model in parallel to the EKF and updating the evolution model with certain periodicity. The evolution model transition matrix is identified using the history of the estimated resistivity distribution obtained by a sensitivity matrix based algorithm and a Newton-Raphson algorithm. To numerically identify the linear evolution model, the Ibrahim time-domain method is used. The investigation is performed by numerical simulations of a domain with time-varying resistivity and by experimental data collected from the boundary of a human chest during normal breathing. The obtained dynamic resistivity values lie within the expected values for the tissues of a human chest. The EKF results suggest that the tracking ability is significantly improved with this approach.

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.

Particle size effect on abrasion resistance of mottled cast iron with different retained austenite contents

Effects of particle abrasive sizes on wear resistance of mottled cast iron with different retained austenite contents were studied. Abrasive wear tests using a pin test on alumina paper were carried out, using abrasive sizes between 16 mu m and 192 mu m. Retained austenite content of the matrix was determined by X-ray diffraction. The wear surface of samples and the alumina paper were examined by scanning electron microscopy for identifying the wear micromechanism. The results show that at lower abrasive sizes the mass loss was similar for the iron with different austenite contents. However, at higher abrasive sizes the samples with higher retained austenite content presented higher abrasion resistance. For lower abrasive sizes tested, samples with higher and lower retained austenite content both presented microcutting. On the other hand, the main wear micromechanism for the samples with higher retained austenite content and higher abrasive sizes was microploughing. The samples with lower retained austenite content presented microcutting and wedge formation at higher abrasive sizes. Higher abrasive size induced more microcutting in samples with lower retained austenite. The iron with lower retained austenite content presented wider grooves for the different abrasive sizes measured. SEM on the abrasive paper used on samples with higher retained austenite showed continuous and discontinuous microchips and the samples with lower retained austenite showed discontinuous microchips at 66 and 141 mu m. This research demonstrates the relation between abrasive size, wear resistance, groove width and wear micromechanism for mottled cast iron with different retained austenite contents. (C) 2009 Elsevier B.V. All rights reserved.

Abrasive wear study of white cast iron with different solidification rates

The abrasive wear resistance of white cast iron was studied. The iron was solidified using two solidification rates of 1.5 and 15 degrees C/s. Mass loss was evaluated with tests of the type pin on abrasive disc using alumina of different sizes. Two matrices were tested: one predominantly austenitic and the other predominantly martensitic, containing M(3)C carbides. Samples with cooling rate of 15 degrees C/s showed higher hardness and more refined microstructure compared with those solidified at 1.5 degrees C/s. During the test, the movement of successive abrasives gave rise to the strain hardening of the austenite phase, leading to the attainment of similar levels of surface hardness, which explains why the wear rate showed no difference compared to the austenite samples with different solidification rates. For the austenitic matrix the wear rate seems to depend on the hardness of the worn surface and not on the hardness of the material without deformation. The austenitic samples showed cracking and fracture of M(3)C carbides. For the predominantly martensitic matrix, the wear rate was higher at the solidification rate of 1.5 degrees C/s, for grain size of 66 and 93 mu m. Higher abrasive sizes were found to produce greater penetration and strain hardening of austenitic matrices. However, martensitic iron produces more microcutting, increasing the wear rate of the material. The analysis of the worn surface by scanning electron microscopy indicated abrasive wear mechanisms such as: microcutting, microfatigue and microploughing. Yet, for the iron of austenitic matrix, the microploughing mechanism was more severe. (C) 2009 Elsevier B.V. All rights reserved.

Load effect in abrasive wear mechanism of cast iron with graphite and cementite

In this study four irons were casted with different chromium and vanadium contents: 2.66% Cr, 5.01% Cr, 2.51% V and 5.19% V. Their microstructure is composed of: ledeburite, graphite and M(3)C carbides (cementite). Pin-abrasion tests were carried out using fixed alumina abrasive grains at different loads: 1, 2, 4.6 and 10 N. The wear surface and the abrasive paper were examined by scanning electron microscopy for identifying the wear micromechanism. The results reveal that the mass loss increased with the load increase, and the effect of the percentage of chromium on mass loss is inverted when the load is increased from 4.6 to 10 N; for 4.6 N the mass loss decreased when the chromium percentage was increased from 2.66% to 5.01%. Nevertheless, for 10 N the mass loss increased when the chromium percentage was increased. The worn surfaces of the materials tested at 1 N show microcutting caused by the abrasive tip that produces continuous microchips. The worn surfaces and the abrasive paper tested at 10 N show continuous microchips and brittle debris. The results show that high pressures produce a brittle wear mechanism and low pressures produce a more ductile wear micromechanism, for this, the applied pressure defines the dependence between the wear resistance and wear micromechanism. (C) 2009 Elsevier B.V. All rights reserved.

Safety of a railway wheelset - derailment simulation with increasing lateral force

The motivation for this research is to make a comparison between dynamic results of a free railway wheelset derailment and safety limits. For this purpose, a numerical simulation of a wheelset derailment submitted to increasing lateral force is used to compare with the safety limit, using different criteria. A simplified wheelset model is used to simulate derailments with different adhesion conditions. The contact force components, including the longitudinal and spin effects, are identified in a steady-state condition on the verge of a derailment. The contact force ratios are used in a three-dimensional (3D) analytical formula to calculate the safety limits. Simulation results obtained with two contact methods were compared with the published results and the safety limit was identified with the two criteria. Results confirm Nadal`s conservative aspect and show that safety 3D analytical formula presents slightly higher safety limits for lower friction coefficients and smaller limits for high friction, in comparison with the simulation results with Fastsim.

Simulation of the Barkhausen Noise using random field Ising model with long-range interaction

We present a method to simulate the Magnetic Barkhausen Noise using the Random Field Ising Model with magnetic long-range interaction. The method allows calculating the magnetic flux density behavior in particular sections of the lattice reticule. The results show an internal demagnetizing effect that proceeds from the magnetic long-range interactions. This demagnetizing effect induces the appearing of a magnetic pattern in the region of magnetic avalanches. When compared with the traditional method, the proposed numerical procedure neatly reduces computational costs of simulation. (c) 2008 Published by Elsevier B.V.

Integral Piezoactuator System with Optimum Placement of Functionally Graded Material - A Topology Optimization Paradigm

Piezoactuators consist of compliant mechanisms actuated by two or more piezoceramic devices. During the assembling process, such flexible structures are usually bonded to the piezoceramics. The thin bonding layer(s) between the compliant mechanism and the piezoceramic may induce undesirable behavior, including unusual interfacial nonlinearities. This constitutes a drawback of piezoelectric actuators and, in some applications, such as those associated to vibration control and structural health monitoring (e. g., aircraft industry), their use may become either unfeasible or at least limited. A possible solution to this standing problem can be achieved through the functionally graded material concept and consists of developing `integral piezoactuators`, that is those with no bonding layer(s) and whose performance can be improved by tailoring their structural topology and material gradation. Thus, a topology optimization formulation is developed, which allows simultaneous distribution of void and functionally graded piezoelectric materials (including both piezo and non-piezoelectric materials) in the design domain in order to achieve certain specified actuation movements. Two concurrent design problems are considered, that is the optimum design of the piezoceramic property gradation, and the design of the functionally graded structural topology. Two-dimensional piezoactuator designs are investigated because the applications of interest consist of planar devices. Moreover, material gradation is considered in only one direction in order to account for manufacturability issues. To broaden the range of such devices in the field of smart structures, the design of integral Moonie-type functionally graded piezoactuators is provided according to specified performance requirements.