BEM formulation for von Karman plates

This work deals with nonlinear geometric plates in the context of von Karman`s theory. The formulation is written such that only the boundary in-plane displacement and deflection integral equations for boundary collocations are required. At internal points, only out-of-plane rotation, curvature and in-plane internal force representations are used. Thus, only integral representations of these values are derived. The nonlinear system of equations is derived by approximating all densities in the domain integrals as single values, which therefore reduces the computational effort needed to evaluate the domain value influences. Hyper-singular equations are avoided by approximating the domain values using only internal nodes. The solution is obtained using a Newton scheme for which a consistent tangent operator was derived. (C) 2009 Elsevier Ltd. All rights reserved.

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 extended BEM formulation for plates reinforced by rectangular beams

This article presents a BEM formulation developed particularly for analysis of plates reinforced by rectangular beams. This is an extended version of a Previous paper that only took into account bending effects. The problem is now re-formulated to consider bending and membrane force effects. The effects of the reinforcements are taken into account by using a simplified scheme that requires application of ail initial stress field to locally correct the bending and stretching stiffness of the reinforcement regions. The domain integrals due to the presence of the reinforcements are then transformed to the reinforcement/plate interface. To reduce the number of degrees of freedom related to the presence of the reinforcement, the proposed model was simplified to consider only bending and stretching rigidities in the direction of the beams. The complete model can be recovered by applying all six internal force correctors, corresponding to six degrees of freedom per node. Examples are presented to confirm the accuracy of the formulation and to illustrate the level of simplification introduced by this strong reduction in the number of degrees of freedom. (C) 2008 Elsevier Ltd. All rights reserved.

BEM formulation for reinforced plates

This article presents a BEM formulation developed to analyse reinforced plate bending. The reinforcements are formulated using a simplified scheme based on applying an initial moment field adopted to locally correct the stiffness of the reinforcement regions. The domain integrals due to the presence of the reinforcements are then transformed to the reinforcement/plate interface. The increase in system stiffness due to the reinforcements can be taken into account independently for each coefficient. Thus, one can conveniently reduce the number of degrees of freedom required in considering the reinforcement. Only one degree-of-freedom is required at each internal node when taking into account only the flexural stiffness of beams. Examples are presented to confirm the accuracy of the formulation. (C) 2009 Elsevier Ltd. All rights reserved.

Multiple random crack propagation using a boundary element formulation

This paper proposes a boundary element method (BEM) model that is used for the analysis of multiple random crack growth by considering linear elastic fracture mechanics problems and structures subjected to fatigue. The formulation presented in this paper is based on the dual boundary element method, in which singular and hyper-singular integral equations are used. This technique avoids singularities of the resulting algebraic system of equations, despite the fact that the collocation points coincide for the two opposite crack faces. In fracture mechanics analyses, the displacement correlation technique is applied to evaluate stress intensity factors. The maximum circumferential stress theory is used to evaluate the propagation angle and the effective stress intensity factor. The fatigue model uses Paris` law to predict structural life. Examples of simple and multi-fractured structures loaded until rupture are considered. These analyses demonstrate the robustness of the proposed model. In addition, the results indicate that this formulation is accurate and can model localisation and coalescence phenomena. (C) 2010 Elsevier Ltd. All rights reserved.

On the analysis of notched concrete beams: From measurement with digital image correlation to identification with boundary element method of a cohesive model

A way of coupling digital image correlation (to measure displacement fields) and boundary element method (to compute displacements and tractions along a crack surface) is presented herein. It allows for the identification of Young`s modulus and fracture parameters associated with a cohesive model. This procedure is illustrated to analyze the latter for an ordinary concrete in a three-point bend test on a notched beam. In view of measurement uncertainties, the results are deemed trustworthy thanks to the fact that numerous measurement points are accessible and used as entries to the identification procedure. (C) 2010 Elsevier Ltd. All rights reserved.

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.

Constitutive model for long term municipal solid waste mechanical behavior

This paper presents some improvements in the model proposed by Machado et al. [Machado SL, Carvalho MF, Vilar OM. Constitutive model for municipal solid waste. J Geotech Geoenviron Eng ASCE 2002; 128(11):940-51] now considering the influence of biodegradation of organic matter in the mechanical behavior of municipal solid waste. The original framework considers waste as composed of two component groups; fibers and organic paste. The particular laws of behavior are assessed for each component group and then coupled to represent waste behavior. The improvements introduced in this paper take into account the changes in the properties of fibers and mass loss due to organic matter depletion over time. Mass loss is indirectly calculated considering the MSW gas generation potential through a first order decay model. It is shown that as the biodegradation process occurs the proportion of fibers increases, however, they also undergo a degradation process which tends to reduce their ultimate tensile stress and Young modulus. The way these changes influence the behavior of MSW is incorporated in the final framework which captures the main features of the MSW stress-strain behavior under different loading conditions. (C) 2007 Elsevier Ltd. All rights reserved.

Internal Residual Stresses in Sintered and Commercial Low Expansion Li(2)O-Al(2)O(3)-SiO(2) Glass-Ceramics

We performed Synchrotron X-ray diffraction (XRD) analyses of internal residual stresses in monolithic samples of a newly developed Li(2)O-Al(2)O(3)-SiO(2) (LAS) glass-ceramic produced by sintering and in a commercial LAS glass-ceramic, CERAN (R), produced by the traditional crystal nucleation and growth treatments. The elastic constants were measured by instrumented indentation and a pulse-echo technique. The thermal expansion coefficient of virgilite was determined by high temperature XRD and dilatometry. The c-axis contracts with the increasing temperature whereas the a-axis does not vary significantly. Microcracking of the microstructure affects the thermal expansion coefficients measured by dilatometry and thermal expansion hysteresis is observed for the sintered glass-ceramic as well as for CERAN (R). The measured internal stress is quite low for both glass-ceramics and can be explained by theoretical modeling if the high volume fraction of the crystalline phase (virgilite) is considered. Using a modified Green model, the calculated critical (glass) island diameter for spontaneous cracking agreed with experimental observations. The experimental data collected also allowed the calculation of the critical crystal grain diameters for grain-boundary microcracking due to the anisotropy of thermal expansion of virgilite and for microcracking in the residual glass phase surrounding the virgilite particles. All these parameters are important for the successful microstructural design of sintered glass-ceramics.

Analytical solutions for geometrically nonlinear trusses

This paper presents an analytical method for analyzing trusses with severe geometrically nonlinear behavior. The main objective is to find analytical solutions for trusses with different axial forces in the bars. The methodology is based on truss kinematics, elastic constitutive laws and equilibrium of nodal forces. The proposed formulation can be applied to hyper elastic materials, such as rubber and elastic foams. A Von Mises truss with two bars made by different materials is analyzed to show the accuracy of this methodology.

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.

Hybrid Reinforcement of Sisal and Polypropylene Fibers in Cement-Based Composites

Several studies using vegetable fibers as the exclusive reinforcement in fiber-cement composites have shown acceptable mechanical performance at the first ages. However, after the exposure to accelerated aging tests, these composites have shown significant reduction in the toughness or increase in embrittlement. This was mainly attributed to the improved fiber-matrix adhesion and fiber mineralization after aging process. The objective of the present research was to evaluate composites produced by the slurry dewatering technique followed by pressing and air curing, reinforced with combinations of polypropylene fibers and sisal kraft pulp at different pulp freeness. The physical properties, mechanical performance, and microstructural characteristics of the composites were evaluated before and after accelerated and natural aging. Results showed the great contribution of pulp refinement on the improvement of the mechanical strength in the composites. Higher intensities of refinement resulted in higher modulus of rupture for the composites with hybrid reinforcement after accelerated and natural aging. The more compact microstructure was due to the improved packing of the mineral particles with refined sisal pulp. The toughness of the composites after aging was maintained in relation to the composites at 28 days of cure.

Exploring the potential of functionally graded materials concept for the development of fiber cement

In this study we establish the concept of functionally graded fiber cement. We discuss the use of statistical mixture designs to choose formulations and present ideas for the production of functionally graded fiber cement components for Hatschek machines. The feasibility of producing functionally graded fiber cement by grading PVA fiber content has been experimentally evaluated. Thermogravimetric analysis (TG) was employed to assess fiber distribution profiles and four-point bending tests were applied to evaluate the mechanical performance of both conventional and graded composites. The results show that grading PVA fiber content is an effective way to produce functionally graded fiber cement, which allows for a reduction of the total fiber volume without a significant reduction on modulus of rupture of composite. TG tests were found adequate to assess the fiber content at different points in functionally graded fiber cements. (C) 2009 Elsevier Ltd. All rights reserved.

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.