Analyzing static three-dimensional elastic domains with a new infinite boundary element formulation

A new two-dimensionally mapped infinite boundary element (IBE) is presented. The formulation is based on a triangular boundary element (BE) with linear shape functions instead of the quadrilateral IBEs usually found in the literature. The infinite solids analyzed are assumed to be three-dimensional, linear-elastic and isotropic, and Kelvin fundamental solutions are employed. One advantage of the proposed formulation over quadratic or higher order elements is that no additional degrees of freedom are added to the original BE mesh by the presence of the IBEs. Thus, the IBEs allow the mesh to be reduced without compromising the accuracy of the result. Two examples are presented, in which the numerical results show good agreement with authors using quadrilateral IBEs and analytical solutions. (C) 2010 Elsevier Ltd. All rights reserved.

Eucalyptus pulp fibres as alternative reinforcement to engineered cement-based composites

This paper evaluates the advantages of using hardwood short fibre pulp (eucalyptus) as alternative to softwood long fibre pulp (pinus) and polymer fibres, traditionally used in reinforcement of cement-based materials. The effects of cellulose fibre length on microstructure and on mechanical performance of fibre-cement composites were evaluated before and after accelerated ageing cycles. Hardwood pulp fibres were better dispersed in the cement matrix and provided higher number of fibres per unitary weight or volume, in relation to softwood long fibre pulp. The short reinforcing elements lead to an effective crack bridging of the fragile matrix, which contributes to the improvement of the mechanical performance of the composite after ageing. These promising results show the potential of eucalyptus short fibres for reducing costs by both the partial replacement of expensive synthetic fibres in air curing process and the energy savings during pulp refining. (C) 2009 Elsevier B.V. All rights reserved.

Cellulose modified fibres in cement based composites

The objective of the present work is to evaluate the effect of surface modification of cellulose pulp fibres on the mechanical and microstructure of fibre-cement composites. Surface modification of the cellulose pulps was performed with Methacryloxypropyltri-methoxysilane (MPTS) and Aminopropyltri-ethoxysilane (APTS) in an attempt to improve their durability into fibre-cement composites. The surface modification showed significant influence on the microstructure of the composites on the fibre-matrix interface and in the mineralization of the fibre lumen as seen by scanning electron microscopy (SEM) with back-scattered electron (BSE) detector. Accelerated ageing cycles decreased modulus of rupture (MOR) and toughness (TE) of the composites. Composites reinforced with MPTS-modified fibres presented fibres free from cement hydration products, while APTS-modified fibres presented accelerated mineralization. Higher mineralization of the fibres led to higher embrittlement of the composite after accelerated ageing cycles. These observations are therefore very useful for understanding the mechanisms of degradation of fibre-cement composites. (C) 2009 Elsevier Ltd. All rights reserved.

An alternative multi-region BEM technique for three-dimensional elastic problems

The main objective of this work is to present an alternative boundary element method (BEM) formulation for the static analysis of three-dimensional non-homogeneous isotropic solids. These problems can be solved using the classical boundary element formulation, analyzing each subregion separately and then joining them together by introducing equilibrium and displacements compatibility. Establishing relations between the displacement fundamental solutions of the different domains, the alternative technique proposed in this paper allows analyzing all the domains as one unique solid, not requiring equilibrium or compatibility equations. This formulation also leads to a smaller system of equations when compared to the usual subregion technique, and the results obtained are even more accurate. (C) 2008 Elsevier Ltd. All rights reserved.

Consolidation of elastic-plastic saturated porous media by the boundary element method

This paper presents a domain boundary element formulation for inelastic saturated porous media with rate-independent behavior for the solid skeleton. The formulation is then applied to elastic-plastic behavior for the solid. Biot`s consolidation theory, extended to include irreversible phenomena is considered and the direct boundary element technique is used for the numerical solution after time discretization by the implicit Euler backward algorithm. The associated nonlinear algebraic problem is solved by the Newton-Raphson procedure whereby the loading/unloading conditions are fully taken into account and the consistent tangent operator defined. Only domain nodes (nodes defined inside the domain) are used to represent all domain values and the corresponding integrals are computed by using an accurate sub-elementation scheme. The developments are illustrated through the Drucker-Prager elastic-plastic model for the solid skeleton and various examples are analyzed with the proposed algorithms. (c) 2008 Elsevier B.V. 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.

Characterization of the elastic-plastic region in AISI/SAE 1070 steel by the magnetic barkhausen noise

The magnetic Barkhausen energy in the rolling and transversal directions of AISI/SAE 1070 annealed surfaces is studied. The measurements were made in the samples under applied tension in the elastic-plastic region for different angular directions. The outcomes evidence that the magnetic anisotropy coefficient can be used to characterize the linear and nonlinear elastic limits of the material tinder tensile tresses. The results also show that the area of the curve corresponding to the angular dependence of the number of Barkhausen jumps with average energy presents a maximum value that corresponds to the elastic limit of the sample. (C) 2008 Elsevier Ltd. All rights reserved.

MEASUREMENT OF ELASTIC PROPERTIES OF MATERIALS BY THE ULTRASONIC THROUGH-TRANSMISSION TECHNIQUE

The elastic mechanical behavior of elastic materials is modeled by a pair of independent constants (Young`s modulus and Poisson`s coefficient). A precise measurement for both constants is necessary in some applications, such as the quality control of mechanical elements and standard materials used for the calibration of some equipment. Ultrasonic techniques have been used because wave velocity depends on the elastic properties of the propagation medium. The ultrasonic test shows better repeatability and accuracy than the tensile and indentation test. In this work, the theoretical and experimental aspects related to the ultrasonic through-transmission technique for the characterization of elastic solids is presented. Furthermore, an amorphous material and some polycrystalline materials were tested. Results have shown an excellent repeatability and numerical errors that are less than 3% in high-purity samples.

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.

Numerical study of tensile tests conducted on systems with elastic-plastic films deposited onto elastic-plastic substrates

In this work, a series of two-dimensional plane-strain finite element analyses was conducted to further understand the stress distribution during tensile tests on coated systems. Besides the film and the substrate, the finite element model also considered a number of cracks perpendicular to the film/substrate interface. Different from analyses commonly found in the literature, the mechanical behavior of both film and substrate was considered elastic-perfectly plastic in part of the analyses. Together with the film yield stress and the number of film cracks, other variables that were considered were crack tip geometry, the distance between two consecutive cracks and the presence of an interlayer. The analysis was based on the normal stresses parallel to the loading axis (sigma(xx)), which are responsible for cohesive failures that are observed in the film during this type of test. Results indicated that some configurations studied in this work have significantly reduced the value of sigma(xx) at the film/substrate interface and close to the pre-defined crack tips. Furthermore, in all the cases studied the values of sigma(xx) were systematically larger at the film/substrate interface than at the film surface. (C) 2010 Elsevier B.V. All rights reserved.

Effects of Mg2+ on Ca2+ release from sarcoplasmic reticulum of skeletal muscle fibres from yabby (crustacean) and rat

1. The role of myoplasmic [Mg2+] on Ca2+ release from the sarcoplasmic reticulum (SR) was examined in the two major types of crustacean muscle fibres, the tonic, long sarcomere fibres and the phasic, short sarcomere fibres of the fresh mater decapod crustacean Cherax: destructor (yabby) and in the fast-twitch rat muscle fibres using the mechanically skinned muscle fibre preparation. 2. A robust Ca2+-induced Ca2+-release (CICR) mechanism was present in both long and short sarcomere fibres and 1 mM Mg2+ exerted a strong inhibitory action on the XR Ca2+ release in both fibre types. 3. The XR displayed different properties with respect to Ca2+ loading in the long and the short sarcomere fibres and marked functional differences were identified with respect to Mg2+ inhibition between the two crustacean fibre types. Thus, in long sarcomere fibres, the submaximally loaded XR was able to release Ca2+ when [Mg2+] was lowered from 1 to 0.01 mw in the presence of 8 mM ATP(total) and in the virtual absence of Ca2+ (< 5 nM) even when the CICR was suppressed. In contrast, negligible Ca2+ was released from the submaximally loaded SR of short sarcomere yabby fibres when [Mg2+] was lowered from 1. to 0.01 mM under the same conditions as for the long sarcomere fibres. Nevertheless, the rate of XR Ca2+ release in short sarcomere fibres increased markedly when [Mg2+] was lowered in the presence of [Ca2+] approaching the normal resting levels (50-100 nM). 4. Rat fibres were able to release SR Ca2+ at a faster rate than the long sarcomere yabby fibres when [Mg2+] was lowered from 1 to 0.01 mM in the virtual absence of Ca2+ but, unlike with yabby fibres, the net rate of Ca2+ release was actually increased for conditions that were considerably less favourable to CICR. 5. In summary it is concluded that crustacean skeletal muscles have more that one functional type of Ca2+-release channels, that these channels display properties that are intermediate between those of mammalian skeletal and cardiac isoforms, that the inhibition exerted by Mg2+ at rest on the crustacean SR Ca2+-release channels must be removed during excitation-contraction coupling and that, unlike in crustacean fibres, CICR cannot play the major role in the activation of XR Ca2+-release channels in the rat skeletal muscle.

Tubular system volume changes in twitch fibres from toad and rat skeletal muscle assessed by confocal microscopy

The volume of the extracellular compartment (tubular system) within intact muscle fibres from cane toad and rat was measured under various conditions using confocal microscopy. Under physiological conditions at rest, the fractional volume of the tubular system (t-sys(Vol)) was 1.38 +/- 0.09% (n = 17),1.41 +/- 0.09% (n = 12) and 0.83 +/- 0.07% (n = 12) of the total fibre volume in the twitch fibres from toad iliofibularis muscle, rat extensor digitorum longus muscle and rat soleus muscle, respectively. In toad muscle fibres, the t-sys(Vol) decreased by 30% when the tubular system was fully depolarized and decreased by 15% when membrane cholesterol was depleted from the tubular system with methyl-beta-cyclodextrin but did not change as the sarcomere length was changed from 1.93 to 3.30 mum. There was also an increase by 30% and a decrease by 25% in t-sys(Vol) when toad fibres were equilibrated in solutions that were 2.5-fold hypertonic and 50% hypotonic, respectively. When the changes in total fibre volume were taken into consideration, the t-sys(Vol) expressed as a percentage of the isotonic fibre volume did actually decrease as tonicity increased, revealing that the tubular system in intact fibres cannot be compressed below 0.9% of the isotonic fibre volume. The results can be explained in terms of forces acting at the level of the tubular wall. These observations have important physiological implications showing that the tubular system is a dynamic membrane structure capable of changing its volume in response to the membrane potential, cholesterol depletion and osmotic stress but not when the sarcomere length is changed in resting muscle.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.