A hybrid-Trefftz formulation for plane elasticity with selective enrichment of the approximations

This work is related to the so-called non-conventional finite element formulations. Essentially, a methodology for the enrichment of the initial approximation which is typical of the meshless methods and based on the clouds concept is introduced in the hybrid-Trefftz formulation for plane elasticity. The formulation presented allows for the approximation and direct enrichment of two independent fields: stresses in the domains and displacements on the boundaries of the elements. Defined by a set of elements and interior boundaries sharing a common node, the cloud notion is employed to select the enrichment support for the approximation fields. The numerical analysis performed reveals an excellent performance of the resulting formulation, characterized by the good approximation ability and a reduced computational effort. Copyright (C) 2009 John Wiley & Sons, Ltd.

A STUDY OF PENALTY FORMULATIONS USED IN THE NUMERICAL APPROXIMATION OF A RADIALLY SYMMETRIC ELASTICITY PROBLEM

We consider a class of two-dimensional problems in classical linear elasticity for which material overlapping occurs in the absence of singularities. Of course, material overlapping is not physically realistic, and one possible way to prevent it uses a constrained minimization theory. In this theory, a minimization problem consists of minimizing the total potential energy of a linear elastic body subject to the constraint that the deformation field must be locally invertible. Here, we use an interior and an exterior penalty formulation of the minimization problem together with both a standard finite element method and classical nonlinear programming techniques to compute the minimizers. We compare both formulations by solving a plane problem numerically in the context of the constrained minimization theory. The problem has a closed-form solution, which is used to validate the numerical results. This solution is regular everywhere, including the boundary. In particular, we show numerical results which indicate that, for a fixed finite element mesh, the sequences of numerical solutions obtained with both the interior and the exterior penalty formulations converge to the same limit function as the penalization is enforced. This limit function yields an approximate deformation field to the plane problem that is locally invertible at all points in the domain. As the mesh is refined, this field converges to the exact solution of the plane problem.

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.

Anomalous compressibility of ferropericlase throughout the iron spin cross-over

The thermoelastic properties of ferropericlase Mg(1-x)Fe(x)O (x = 0.1875) throughout the iron high-to-low spin cross-over have been investigated by first principles at Earth`s lower mantle conditions. This cross-over has important consequences for elasticity such as an anomalous bulk modulus (K(S)) reduction. At room temperature the anomaly is somewhat sharp in pressure but broadens with increasing temperature. Along a typical geotherm it occurs across most of the lower mantle with a more significant K(S) reduction at approximate to 1,400-1,600 km depth. This anomaly might also cause a reduction in the effective activation energy for diffusion creep and lead to a viscosity minimum in the mid-lower mantle, in apparent agreement with results from inversion of data related with mantle convection and postglacial rebound.

Real exchange rate and elasticity of labour supply in a balance-of-payments-constrained macrodynamics

A macrodynamic model is proposed in which the real exchange rate and the elasticity of labour supply interact defining different trajectories of growth and income distribution in a developing economy. Growth depends on imports of capital goods which are paid with exports (there are no capital flows) and hence is constrained by equilibrium in current account. The role of the elasticity of labour supply is to prevent the real exchange rate from appreciating as the economy grows, thereby sustaining international competitiveness. The model allows for endogenous technological change and considers the impact of migration from the subsistence to the modern sector on the cumulative (Kaldor-Verdoorn) process of learning.

Dibucaine effects on structural and elastic properties of lipid bilayers

In this work we report the interaction effects of the local anesthetic dibucaine (DBC) with lipid patches in model membranes by Atomic Force Microscopy (AFM). Supported lipid bilayers (egg phosphatidylcholine, EPC and dimyristoylphosphatidylcholine, DMPQ were prepared by fusion of unilamellar vesicles on mica and imaged in aqueous media. The AFM images show irregularly distributed and sized EPC patches on mica. On the other hand DMPC formation presents extensive bilayer regions on top of which multibilayer patches are formed. In the presence of DBC we observed a progressive disruption of these patches, but for DMPC bilayers this process occurred more slowly than for EPC. In both cases, phase images show the formation of small structures on the bilayer surface suggesting an effect on the elastic properties of the bilayers when DBC is present. Dynamic surface tension and dilatational surface elasticity measurements of EPC and DMPC monolayers in the presence of DBC by the pendant drop technique were also performed, in order to elucidate these results. The curve of lipid monolayer elasticity versus DBC concentration, for both EPC and DMPC cases, shows a maximum for the surface elasticity modulus at the same concentration where we observed the disruption of the bilayer by AFM. Our results suggest that changes in the local curvature of the bilayer induced by DBC could explain the anesthetic action in membranes. (C) 2008 Elsevier B.V. All rights reserved.

The association of Na,K-ATPase subunits studied by circular dichroism, surface tension and dilatational elasticity

Different stoichiometries are observed between alpha and beta subunits of Na,K-ATPase that depend on the method employed to solubilize and purify the enzyme. It is not known whether this variability is due to loss of protein-protein association, or is a result of the replacement of essential phospholipids by detergent molecules. With the aim of understanding the effect of enzyme/surfactant ratio on both the catalytic activity and the enzyme structure, we have investigated the bulk and surface properties of the enzyme. The circular dichroism (CD) spectra, surface tension and dilatational surface elasticity results were compared with the residual ATPase activity of the Na,K-ATPase in different surfactant and protein concentrations. Na,K-ATPase in the (alpha beta)(2) form dissociated to the alpha beta form on dilution, and associated to the (alpha beta)(4) form when concentrated. These different stoichiometries have similar ATPase activities and are in equilibrium at C(12)E(8) concentrations below the CIVIC (0.053 mg mL(-1)). At detergent concentrations above the CIVIC the ATPase activity of all forms was abolished, which is concomitant with the dissociation of the a and subunits. (C) 2008 Elsevier Inc. All rights reserved.

Influence of curing protocol on selected properties of light-curing polymers: Degree of conversion, volume contraction, elastic modulus, and glass transition temperature

Objectives. The purpose of this study was to investigate the effect of light-curing protocol on degree of conversion (DC), volume contraction (C), elastic modulus (E), and glass transition temperature (T(g)) as measured on a model polymer. It was a further aim to correlate the measured values with each other. Methods. Different light-curing protocols were used in order to investigate the influence of energy density (ED), power density (PD), and mode of cure on the properties. The modes of cure were continuous, pulse-delay, and stepped irradiation. DC was measured by Raman micro-spectroscopy. C was determined by pycnometry and a density column. E was measured by a dynamic mechanical analyzer (DMA), and T(g) was measured by differential scanning calorimetry (DSC). Data were submitted to two-and three-way ANOVA, and linear regression analyses. Results. ED, PD, and mode of cure influenced DC, C, E, and T(g) of the polymer. A significant positive correlation was found between ED and DC (r = 0.58), ED and E (r = 0.51), and ED and T(g) (r = 0.44). Taken together, ED and PD were significantly related to DC and E. The regression coefficient was positive for ED and negative for PD. Significant positive correlations were detected between DC and C (r = 0.54), DC and E (r = 0.61), and DC and T(g) (r = 0.53). Comparisons between continuous and pulse-delay modes of cure showed significant influence of mode of cure: pulse-delay curing resulted in decreased DC, decreased C, and decreased T(g). Influence of mode of cure, when comparing continuous and step modes of cure, was more ambiguous. A complex relationship exists between curing protocol, microstructure of the resin and the investigated properties. The overall performance of a composite is thus indirectly affected by the curing protocol adopted, and the desired reduction of C may be in fact a consequence of the decrease in DC. (C) 2009 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

The impact of maternal cortisol concentrations on child arterial elasticity

Epidemiological studies suggest that glucocorticoid excess in the fetus may contribute to the pathophysiology of cardiovascular diseases in adulthood. However, the impact of maternal glucocorticoid on the cardiovascular system of the offspring has not been much explored in studies involving humans, especially in childhood. The objective of this study was to assess the influence of maternal cortisol concentrations on child arterial elasticity. One hundred and thirty pregnant women followed from 1997 to 2000, and respective children 5-7 years of age followed from 2004 to 2006 were included in the study. Maternal cortisol was determined in saliva by an enzyme immunoassay utilizing the mean concentration of nine samples of saliva. Arterial elasticity was assessed by the large artery elasticity index (LAEI; the capacitive elasticity of large arteries) by recording radial artery pulse wave, utilizing the equipment HDI/PulseWave CR-2000 Cardiovascular Profiling System (R). The nutritional status of the children was determined by the body mass index (BMI). Insulin concentration was assessed by chemiluminescence, and insulin resistance by the homeostasis model assessment. Blood glucose, total cholesterol and fractions (LDL-c and HDL-c) and triglyceride concentrations were determined by automated enzymatic methods. The association between maternal cortisol and child arterial elasticity was assessed by multivariate linear regression analysis. There was a statistically significant association between maternal cortisol and LAEI (P=0.02), controlling for birth weight, age, BMI and HDL-c of the children. This study suggests that exposure to higher glucocorticoid concentrations in the prenatal period is associated to lower arterial elasticity in childhood, an earlier cardiovascular risk marker.

Mechanical properties and porosity of dental glass-ceramics hot-pressed at different temperatures

The objective of this work was to evaluate biaxial-flexural-strength (σf), Vickers hardness (HV), fracture toughness (K Ic), Young's modulus (E), Poisson's ratio (ν) and porosity (P) of two commercial glass-ceramics, Empress (E1) and Empress 2 (E2), as a function of the hot-pressing temperature. Ten disks were hot-pressed at 1065, 1070, 1075 and 1080 °C for E1; and at 910, 915, 920 and 925 °C for E2. The porosity was measured by an image analyzer software and s f was determined using the piston-on-three-balls method. K Ic and HV were determined by an indentation method. Elastic constants were determined by the pulse-echo method. For E1 samples treated at different temperatures, there were no statistical differences among the values of all evaluated properties. For E2 samples treated at different temperatures, there were no statistical differences among the values of σf, E, and ν, however HV and K Ic were significantly higher for 910 and 915 °C, respectively. Regarding P, the mean value obtained for E2 for 925 °C was significantly higher compared to other temperatures.

Keratoconus prediction using a finite element model of the cornea with local biomechanical properties

PURPOSE: The ability to predict and understand which biomechanical properties of the cornea are responsible for the stability or progression of keratoconus may be an important clinical and surgical tool for the eye-care professional. We have developed a finite element model of the cornea, that tries to predicts keratoconus-like behavior and its evolution based on material properties of the corneal tissue. METHODS: Corneal material properties were modeled using bibliographic data and corneal topography was based on literature values from a schematic eye model. Commercial software was used to simulate mechanical and surface properties when the cornea was subject to different local parameters, such as elasticity. RESULTS: The simulation has shown that, depending on the corneal initial surface shape, changes in local material properties and also different intraocular pressures values induce a localized protuberance and increase in curvature when compared to the remaining portion of the cornea. CONCLUSIONS: This technique provides a quantitative and accurate approach to the problem of understanding the biomechanical nature of keratoconus. The implemented model has shown that changes in local material properties of the cornea and intraocular pressure are intrinsically related to keratoconus pathology and its shape/curvature.

Relationship between birthweight and arterial elasticity in childhood

There is a considerable debate about the potential influence of fetal programming on cardiovascular diseases in adulthood. In the present prospective epidemiological cohort study, the relationship between birthweight and arterial elasticity in 472 children between 5 and 8 years of age was assessed. LAEI (large artery elasticity index), SAEI (small artery elasticity index) and BP (blood pressure) were assessed using the HDI/PulseWave CR-2000 CardioVascular Profiling System. Blood concentrations of glucose, total cholesterol and its fractions [LDL (low-density lipoprotein)-cholesterol and HDL (high-density lipoprotein)-cholesterol] and triacylglycerols (triglycerides) were determined by automated enzymatic methods. Insulin was assessed by a chemiluminescent method, insulin resistance by HOMA (homoeostasis model assessment) and CRP (C-reactive protein) by immunonephelometry. Two linear regression models were applied to investigate the relationship between the outcomes, LAEI and SAEI, and the following variables: birthweight, gestational age, glucose, LDL-cholesterol, HDL-cholesterol, triacylglycerols, insulin, CRP, HOMA, age, gender, waist circumference, per capita income, SBP (systolic BP) and DBP (diastolic BP). LAEI was positively associated with birthweight (P=0.036), waist circumference (P<0.001) and age (P<0.001), and negatively associated with CRP (P=0.024) and SBP (P<0.001). SAEI was positively associated with birthweight (P=0.04), waist circumference (P=0.001) and age (P<0.001), and negatively associated with DBP (P<0.001). Arterial elasticity was decreased in apparently healthy children who had lower birthweights, indicating an earlier atherogenetic susceptibility to cardiovascular diseases in adolescence and adult life. Possible explanations for the results include changes in angiogenesis during critical phases of intrauterine life caused by periods of fetal growth inhibition and local haemodynamic anomalies

Parametric model for the analysis of concrete expansion due to alkali-aggregate reaction

The alkali-aggregate reaction (AAR) is a chemical reaction that provokes a heterogeneous expansion of concrete and reduces important properties such as Young's modulus, leading to a reduction in the structure's useful life. In this study, a parametric model is employed to determine the spatial distribution of the concrete expansion, combining normalized factors that influence the reaction through an AAR expansion law. Optimization techniques were employed to adjust the numerical results and observations in a real structure. A three-dimensional version of the model has been implemented in a finite element commercial package (ANSYS(C)) and verified in the analysis of an accelerated mortar test. Comparisons were made between two AAR mathematical descriptions for the mechanical phenomenon, using the same methodology, and an expansion curve obtained from experiment. Some parametric studies are also presented. The numerical results compared very well with the experimental data validating the proposed method.

Analysis of the tip roundness effects on the micro- and macroindentation response of elastic-plastic materials

In this work, the effects of indenter tip roundness oil the load-depth indentation curves were analyzed using finite element modeling. The tip roundness level was Studied based on the ratio between tip radius and maximum penetration depth (R/h(max)), which varied from 0.02 to 1. The proportional Curvature constant (C), the exponent of depth during loading (alpha), the initial unloading slope (S), the correction factor (beta), the level of piling-up or sinking-in (h(c)/h(max)), and the ratio h(max)/h(f) are shown to be strongly influenced by the ratio R/h(max). The hardness (H) was found to be independent of R/h(max) in the range studied. The Oliver and Pharr method was successful in following the variation of h(c)/h(max) with the ratio R/h(max) through the variation of S with the ratio R/h(max). However, this work confirmed the differences between the hardness values calculated using the Oliver-Pharr method and those obtained directly from finite element calculations; differences which derive from the error in area calculation that Occurs when given combinations of indented material properties are present. The ratio of plastic work to total work (W(p)/W(t)) was found to be independent of the ratio R/h(max), which demonstrates that the methods for the Calculation of mechanical properties based on the *indentation energy are potentially not Susceptible to errors caused by tip roundness.

Analysis of the effects of conical indentation variables on the indentation response of elastic-plastic materials through factorial design of experiment

In this work, the effects of conical indentation variables on the load-depth indentation curves were analyzed using finite element modeling and dimensional analysis. A factorial design 2(6) was used with the aim of quantifying the effects of the mechanical properties of the indented material and of the indenter geometry. Analysis was based on the input variables Y/E, R/h(max), n, theta, E, and h(max). The dimensional variables E and h(max) were used such that each value of dimensionless Y/E was obtained with two different values of E and each value of dimensionless R/h(max) was obtained with two different h(max) values. A set of dimensionless functions was defined to analyze the effect of the input variables: Pi(1) = P(1)/Eh(2), Pi(2) = h(c)/h, Pi(3) = H/Y, Pi(4) = S/Eh(max), Pi(6) = h(max)/h(f) and Pi(7) = W(P)/W(T). These six functions were found to depend only on the dimensionless variables studied (Y/E, R/h(max), n, theta). Another dimension less function, Pi(5) = beta, was not well defined for most of the dimensionless variables and the only variable that provided a significant effect on beta was theta. However, beta showed a strong dependence on the fraction of the data selected to fit the unloading curve, which means that beta is especially Susceptible to the error in the Calculation of the initial unloading slope.