A solid-like FEM for geometrically non-linear 3D frames

This study presents a solid-like finite element formulation to solve geometric non-linear three-dimensional inhomogeneous frames. To achieve the desired representation, unconstrained vectors are used instead of the classic rigid director triad; as a consequence, the resulting formulation does not use finite rotation schemes. High order curved elements with any cross section are developed using a full three-dimensional constitutive elastic relation. Warping and variable thickness strain modes are introduced to avoid locking. The warping mode is solved numerically in FEM pre-processing computational code, which is coupled to the main program. The extra calculations are relatively small when the number of finite elements. with the same cross section, increases. The warping mode is based on a 2D free torsion (Saint-Venant) problem that considers inhomogeneous material. A scheme that automatically generates shape functions and its derivatives allow the use of any degree of approximation for the developed frame element. General examples are solved to check the objectivity, path independence, locking free behavior, generality and accuracy of the proposed formulation. (C) 2009 Elsevier B.V. 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.

Preloaded RC columns strengthened with high-strength concrete jackets under uniaxial compression

Although use of high-strength reinforced concrete (RC) jackets has become common practice worldwide, there are still two unresolved issues regarding the contribution of the original concrete and the effects of existing loads. Twelve RC-jacketed columns were tested with and without preloading under uniaxial compression. Tests showed the entire core to contribute to the capacity of the jacketed column, as long as adequate confinement is provided. Also, preloading does not adversely affect the capacity of the jacketed column, while it may increase its deformability, especially in square sections. Transverse reinforcement in the jacket directly improves ductility of the strengthened column, especially in circular sections.

Influence of physical and geometrical parameters on three-dimensional load transfer mechanism at tunnel face

Three-dimensional discretizations used in numerical analyses of tunnel construction normally include excavation step lengths much shorter than tunnel cross-section dimensions. Simulations have usually worked around this problem by using excavation steps that are much larger than the actual physical steps used in a real tunnel excavation. In contrast, the analyses performed in this study were based on finely discretized meshes capable of reproducing the excavation lengths actually used in tunnels, and the results obtained for internal forces are up to 100% greater than those found in other analyses available in the literature. Whereas most reports conclude that internal forces depend on support delay length alone, this study shows that geometric path dependency (reflected by excavation round length) is very strong, even considering linear elasticity. Moreover, many other solutions found in the literature have also neglected the importance of the relative stiffness between the ground mass and support structure, probably owing to the relatively coarse meshes used in these studies. The analyses presented here show that relative stiffness may account for internal force discrepancies in the order of 60%. A dimensionless expression that takes all these parameters into account is presented as a good approximation for the load transfer mechanism at the tunnel face.

Performance and molecular evaluation of an anaerobic system with suspended biomass for treating wastewater with high fat content after enzymatic hydrolysis

The effect of a lipase-rich fungal enzymatic preparation, produced by a Penicillium sp. during solid-state fermentation, was evaluated in an anaerobic digester treating dairy wastewater with 1200 mg of oil and grease/L The oil and grease hydrolysis step was carried out with 0.1% (w/v) of solid enzymatic preparation at 30 degrees C for 24 h, and resulted in a final free acid concentration eight times higher than the initial value. The digester operated in sequential batches of 48 h at 30 degrees C for 245 days, and had high chemical oxygen demand (COD) removal efficiencies (around 90%) when fed with pre-hydrolyzed wastewater. However, when the pre-hydrolysis step was removed, the anaerobic digester performed poorly (with an average COD removal of 32%), as the oil and grease accumulated in the biomass and effluent oil and grease concentration increased throughout the operational period. PCR-DGGE analysis of the Bacteria and Archaea domains revealed remarkable differences in the microbial profiles in trials conducted with and without the pre-hydrolysis step, indicating that differences observed in overall parameters were intrinsically related to the microbial diversity of the anaerobic sludge. (C) 2009 Elsevier Ltd. All rights reserved.

Porous alumina-spinel ceramics for high temperature applications

In order to reduce energy costs, high-temperature insulation porous refractory ceramics have been subjected to increasing demands. Among the techniques used to produce these materials (such as the addition of foaming agents and organic compounds), the pore generation via phase transformation presents key aspects, such as easy processing and the absence of toxic volatiles. In this study, this technique was applied to produce porous ceramics by decomposing an aluminum magnesium hydro-carbonate known as hydrotalcite (Mg(6)Al(2)(CO(3))(OH)(16)center dot 4H(2)O). It was found out that by using this complex compound, a large fraction of pores can be generated and kept at high temperatures (above 1300 degrees C) due to the in situ formation of spinel-like phases (MgAl(2)O(4)). (C) 2011 Elsevier Ltd and Techna Group S.r.l. All rights reserved.

Wavelet-Galerkin method for one-dimensional elastoplasticity and damage problems: Constitutive modeling and computational aspects

This work presents an analysis of the wavelet-Galerkin method for one-dimensional elastoplastic-damage problems. Time-stepping algorithm for non-linear dynamics is presented. Numerical treatment of the constitutive models is developed by the use of return-mapping algorithm. For spacial discretization we can use wavelet-Galerkin method instead of standard finite element method. This approach allows to locate singularities. The discrete formulation developed can be applied to the simulation of one-dimensional problems for elastic-plastic-damage models. (C) 2007 Elsevier Inc. All rights reserved.

Proposal of a High-Level Language for Writing Self Modifying Programs

This paper proposes a simple high-level programming language, endowed with resources that help encoding self-modifying programs. With this purpose, a conventional imperative language syntax (not explicitly stated in this paper) is incremented with special commands and statements forming an adaptive layer specially designed with focus on the dynamical changes to be applied to the code at run-time. The resulting language allows programmers to easily specify dynamic changes to their own program`s code. Such a language succeeds to allow programmers to effortless describe the dynamic logic of their adaptive applications. In this paper, we describe the most important aspects of the design and implementation of such a language. A small example is finally presented for illustration purposes.

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.

Suction equilibration time for a high capacity tensiometer

It is well established, the importance of the measurement of soil suction for the assessment of mechanical and hydraulic behavior of unsaturated soils. Among the methods to obtain the soil suction, the tensiometer is one of the most convenient and reliable. However conventional tensiometer has a limitation related to the maximum suction it is capable of measure. This limitation was overcome by Ridley and Burland (1993), with the development of a high capacity tensiometer, which is capable of measure suction well above 100 kPa. The equipment has a quick response time, allowing the determination of suction in minutes. This paper presents a study about the factors that affect the equilibrium time for high capacity tensiometers in the laboratory. Soil specimens were prepared at three different conditions, creating different soil structures. In addition to that an investigation about the characteristic of the interface that is required between the soil sample and the porous ceramic of the tensiometer was carried out; showing the role of the paste on the technique. The results also suggested that it is possible to infer the hydraulic conductivity function using the equilibrium curve obtained during the measurement of the soil suction using the high capacity tensiometer.

Hybrid and multi-field variational principles for geometrically exact three-dimensional beams

This paper addresses the development of several alternative novel hybrid/multi-field variational formulations of the geometrically exact three-dimensional elastostatic beam boundary-value problem. In the framework of the complementary energy-based formulations, a Legendre transformation is used to introduce the complementary energy density in the variational statements as a function of stresses only. The corresponding variational principles are shown to feature stationarity within the framework of the boundary-value problem. Both weak and linearized weak forms of the principles are presented. The main features of the principles are highlighted, giving special emphasis to their relationships from both theoretical and computational standpoints. (C) 2010 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.

Numerical solution of the two-phase expansion of a metastable flashing liquid jet using the dispersion-controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two-dimensional axisymmetric domain, while in the two-phase region. In general, the phenomenon starts off when a high-pressure, high-temperature liquid jet emerges from a small nozzle or orifice expanding into a low-pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two-phase high-speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation-of-state for the substance. A false-transient model is implemented using the shock capturing scheme: dispersion-controlled dissipative (DCD), which was used to calculate the flow conditions as the steady-state condition is reached. Numerical results with computational code DCD-2D vI have been analyzed. Copyright (C) 2009 John Wiley & Sons, Ltd.

A transient three-dimensional heat transfer model of the human body

The objective of this work is to develop an improved model of the human thermal system. The features included are important to solve real problems: 3D heat conduction, the use of elliptical cylinders to adequately approximate body geometry, the careful representation of tissues and important organs, and the flexibility of the computational implementation. Focus is on the passive system, which is composed by 15 cylindrical elements and it includes heat transfer between large arteries and veins. The results of thermal neutrality and transient simulations are in excellent agreement with experimental data, indicating that the model represents adequately the behavior of the human thermal system. (C) 2009 Elsevier Ltd. 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.