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.

Buckling and post-buckling of extensible rods revisited: A multiple-scale solution

An exact non-linear formulation of the equilibrium of elastic prismatic rods subjected to compression and planar bending is presented, electing as primary displacement variable the cross-section rotations and taking into account the axis extensibility. Such a formulation proves to be sufficiently general to encompass any boundary condition. The evaluation of critical loads for the five classical Euler buckling cases is pursued, allowing for the assessment of the axis extensibility effect. From the quantitative viewpoint, it is seen that such an influence is negligible for very slender bars, but it dramatically increases as the slenderness ratio decreases. From the qualitative viewpoint, its effect is that there are not infinite critical loads, as foreseen by the classical inextensible theory. The method of multiple (spatial) scales is used to survey the post-buckling regime for the five classical Euler buckling cases, with remarkable success, since very small deviations were observed with respect to results obtained via numerical integration of the exact equation of equilibrium, even when loads much higher than the critical ones were considered. Although known beforehand that such classical Euler buckling cases are imperfection insensitive, the effect of load offsets were also looked at, thus showing that the formulation is sufficiently general to accommodate this sort of analysis. (c) 2008 Elsevier Ltd. All rights reserved.

Computational evaluation of flexural toughness of FRC and fracture properties of plain concrete

In this paper, a computational tool concerning the computation of flexural and fracture toughness of cement based composites is presented. Firstly, RILEM`s (Reunion Internationale des Laboratoires d`Essais de Materiaux) recommendations related to the analysis of FRC in three-point bend tests are discussed in their relevant aspects regarding the computational implementations. The determination of other mechanical properties such as the Young modulus has been added to the program. Taking this into account, a new formulation based on displacements is used. In the second part of the paper, the determination of fracture properties of concrete, such as the fracture energy, G(F) , and the fracture toughness, K-IC(S), is discussed regarding the computational strategies used in the implementations. Several features whereby anterior data can be reanalyzed, obtained from other standards and recommendations, have been incorporated into the program, therefore allowing comparative studies and back analysis activities.

Short-term experimental data of drying shrinkage of ground granulated blast-furnace slag cement concrete

This paper presents results of laboratory testing of unrestrained drying shrinkage during a period of 154 days of different concrete mixtures from the Brazilian production line that utilize ground granulated blast-furnace slag in their compositions. Three concrete mixtures with water/cement ratio of 0.78(M1), 0.41(M2), and 0.37(M3) were studied. The obtained experimental data were compared with the analytical results from prediction models available in the literature: the ACI 209 model (ACI), the B3 model (B3), the Eurocode 2 model (EC2), the GL 2000 model (GL), and the Brazilian NBR 6118 model (NBR), and an analysis of the efficacy of these models was conducted utilizing these experimental data. In addition, the development of the mechanical properties (compressive strength and modulus of elasticity) of the studied concrete mixtures was also measured in the laboratory until 126 days. From this study, it could be concluded that the ACI and the GL were the models that most approximated the experimental drying shrinkage data measured during the analyzed period of time.

3D CFD Simulation of Vortex-induced Vibration of Cylinder

The 3D flow around a circular cylinder free to oscillate transversely to the free stream was simulated using Computational Fluid Dynamics (CFD) and the Spalart-Allmaras Detached Eddy Simulation (DES) turbulence model for a Reynolds number Re = 10(4). Simulations were carried out for a small mass-damping parameter m*zeta = 0.00858, where m* = 3.3 and zeta = 0.0026. We found good agreement between the numerical results and experimental data. The simulations predicted the high observed amplitudes of the upper branch of vortex-induced vibrations for low mass-damping parameters.