A boundary element formulation for analysis of elastoplastic plates with geometrical nonlinearity

In this paper a new boundary element method formulation for elastoplastic analysis of plates with geometrical nonlinearities is presented. The von Mises criterion with linear isotropic hardening is considered to evaluate the plastic zone. Large deflections are assumed but within the context of small strain. To derive the boundary integral equations the von Karman`s hypothesis is taken into account. An initial stress field is applied to correct the true stresses according to the adopted criterion. Isoparametric linear elements are used to approximate the boundary unknown values while triangular internal cells with linear shape function are adopted to evaluate the domain value influences. The nonlinear system of equations is solved by using an implicit scheme together with the consistent tangent operator derived along the paper. Numerical examples are presented to demonstrate the accuracy and the validity of the proposed formulation.

Non-linear boundary element formulation with tangent operator to analyse crack propagation in quasi-brittle materials

This work deals with analysis of cracked structures using BEM. Two formulations to analyse the crack growth process in quasi-brittle materials are discussed. They are based on the dual formulation of BEM where two different integral equations are employed along the opposite sides of the crack surface. The first presented formulation uses the concept of constant operator, in which the corrections of the nonlinear process are made only by applying appropriate tractions along the crack surfaces. The second presented BEM formulation to analyse crack growth problems is an implicit technique based on the use of a consistent tangent operator. This formulation is accurate, stable and always requires much less iterations to reach the equilibrium within a given load increment in comparison with the classical approach. Comparison examples of classical problem of crack growth are shown to illustrate the performance of the two formulations. (C) 2009 Elsevier Ltd. All rights reserved.

Analysis of stiffened plates by the boundary element method

In this paper, a formulation for representation of stiffeners in plane stress by the boundary elements method (BEM) in linear analysis is presented. The strategy is to adopt approximations for the displacements in the central line of the stiffener. With this simplification the Spurious oscillations in the stress along stiffeners with small thickness is prevented. Worked examples are analyzed to show the efficiency of these techniques, especially in the insertion of very narrow sub-regions, in which quasi-singular integrals are calculated, with stiffeners that are much stiffer than the main domain. The results obtained with this formulation are very close to those obtained with other formulations. (C) 2007 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.

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.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

Long-term stability of hydrogen and organic acids production in an anaerobic fluidized-bed reactor using heat treated anaerobic sludge inoculum

This study evaluates the stability of hydrogen and organic acids production in an anaerobic fluidized-bed reactor (AFBR) that contains expanded clay (2.8-3.35 mm in diameter) as a support medium and is operated on a long-term basis. The reactor was inoculated with thermally pre-treated anaerobic sludge and operated with decreasing hydraulic retention time (HRT), from 8 h to 1 h, at a controlled temperature of 30 degrees C and a pH of about 3.8. Glucose (2000 mg L(-1)) was used as the substrate, generating conversion rates of 92-98%. Decreasing the HRT from 8 h to 1 h led to an increase in average hydrogen-production rates, with a maximum value of 1.28 L h(-1) L(-1) for an HRT of 1 h. In general, hydrogen yield production increased as HRT decreased, reaching 2.29 mol of H(2)/mol glucose at an HRT of 2 h and yielding a maximum hydrogen content of 37% in the biogas. No methane was detected in the biogas throughout the period of operation. The main soluble metabolites (SMP) were acetic acid (46.94-53.84% of SMP) and butyric acid (34.51-42.16% of SMP), with less than 15.49% ethanol. The steady performance of the AFBR may be attributed to adequate thermal treatment of the inoculum, the selection of a suitable support medium for microbial adhesion, and the choice of satisfactory environmental conditions imposed on the system. The results show that stable hydrogen production and organic acids production were maintained in the AFBR over a period of 178 days. (C) 2009 International Association for Hydrogen Energy. Published by Elsevier Ltd. All rights reserved.

Experimental study of a Tollmien-Schlichting wave interacting with a shallow 3D roughness element

The paper is devoted to an experimental study of the effect of a shallow 3D roughness element on the evolution of a 2D Tollmien-Schlichting wave in a Blasius boundary layer. The experiments were carried out under controlled disturbance conditions on an airfoil section which could provide a long run with zero pressure gradient flow. A pneumatically driven slit source was used to introduce the Tollmien-Schilichting wave upstream of the lower branch of the neutral stability curve. A few wavelengths downstream, the T-S wave interacts with a cylindrical roughness element. The height of the roughness was slowly oscillating in time, which allows a continuous measurement of the T-S wave response downstream the roughness. The oscillation frequency was approximately 1500 times lower than the frequency of the studied Tollmien-Schlichting wave and therefore, behaved as a steady roughness with respect to the T-S wave. Hot wire anemometry was used to measure wall normal profiles and spanwise scans close to the maximum of the eigenfunction of the T-S wave. The oscillation of the roughness and the synchronization of all-equipments permitted the use of ensemble average techniques. Two different amplitudes of T-S waves with a non-dimensional frequency of F120E-06 were studied. They show a strong amplification of the disturbances in a small spanwise wave number range. The analysis of the wall normal T-S profiles suggests the growth of oblique modes.

A 2-D Delaunay refinement algorithm using an initial prerefinement from the boundary mesh

An improvement to the quality bidimensional Delaunay mesh generation algorithm, which combines the mesh refinement algorithms strategy of Ruppert and Shewchuk is proposed in this research. The developed technique uses diametral lenses criterion, introduced by L. P. Chew, with the purpose of eliminating the extremely obtuse triangles in the boundary mesh. This method splits the boundary segment and obtains an initial prerefinement, and thus reducing the number of necessary iterations to generate a high quality sequential triangulation. Moreover, it decreases the intensity of the communication and synchronization between subdomains in parallel mesh refinement.

Unsaturated Infinite Slope Stability Considering Surface Flux Conditions

A slope stability model is derived for an infinite slope subjected to unsaturated infiltration flow above a phreatic surface. Closed form steady state solutions are derived for the matric suction and degree of saturation profiles. Soil unit weight, consistent with the degree of saturation profile, is also directly calculated and introduced into the analyzes, resulting in closed-form solutions for typical soil parameters and an infinite series solution for arbitrary soil parameters. The solutions are coupled with the infinite slope stability equations to establish a fully realized safety factor function. In general, consideration of soil suction results in higher factor of safety. The increase in shear strength due to the inclusion of soil suction is analogous to making an addition to the cohesion, which, of course, increases the factor of safety against sliding. However, for cohesive soils, the results show lower safety factors for slip surfaces approaching the phreatic surface compared to those produced by common safety factor calculations. The lower factor of safety is due to the increased soil unit weight considered in the matric suction model but not usually accounted for in practice wherein the soil is treated as dry above the phreatic surface. The developed model is verified with a published case study, correctly predicting stability under dry conditions and correctly predicting failure for a particular storm.

Post-processing and visualization techniques in 2D boundary element analysis

Most post-processors for boundary element (BE) analysis use an auxiliary domain mesh to display domain results, working against the profitable modelling process of a pure boundary discretization. This paper introduces a novel visualization technique which preserves the basic properties of the boundary element methods. The proposed algorithm does not require any domain discretization and is based on the direct and automatic identification of isolines. Another critical aspect of the visualization of domain results in BE analysis is the effort required to evaluate results in interior points. In order to tackle this issue, the present article also provides a comparison between the performance of two different BE formulations (conventional and hybrid). In addition, this paper presents an overview of the most common post-processing and visualization techniques in BE analysis, such as the classical algorithms of scan line and the interpolation over a domain discretization. The results presented herein show that the proposed algorithm offers a very high performance compared with other visualization procedures.

An EFG method for the nonlinear analysis of plates undergoing arbitrarily large deformations

The applicability of a meshfree approximation method, namely the EFG method, on fully geometrically exact analysis of plates is investigated. Based on a unified nonlinear theory of plates, which allows for arbitrarily large rotations and displacements, a Galerkin approximation via MLS functions is settled. A hybrid method of analysis is proposed, where the solution is obtained by the independent approximation of the generalized internal displacement fields and the generalized boundary tractions. A consistent linearization procedure is performed, resulting in a semi-definite generalized tangent stiffness matrix which, for hyperelastic materials and conservative loadings, is always symmetric (even for configurations far from the generalized equilibrium trajectory). Besides the total Lagrangian formulation, an updated version is also presented, which enables the treatment of rotations beyond the parameterization limit. An extension of the arc-length method that includes the generalized domain displacement fields, the generalized boundary tractions and the load parameter in the constraint equation of the hyper-ellipsis is proposed to solve the resulting nonlinear problem. Extending the hybrid-displacement formulation, a multi-region decomposition is proposed to handle complex geometries. A criterium for the classification of the equilibrium`s stability, based on the Bordered-Hessian matrix analysis, is suggested. Several numerical examples are presented, illustrating the effectiveness of the method. Differently from the standard finite element methods (FEM), the resulting solutions are (arbitrary) smooth generalized displacement and stress fields. (c) 2007 Elsevier Ltd. All rights reserved.

Numerical prediction of gas concentrations and fluctuations above a triangular hill within a turbulent boundary layer

The objective of the present work is to propose a numerical and statistical approach, using computational fluid dynamics, for the study of the atmospheric pollutant dispersion. Modifications in the standard k-epsilon turbulence model and additional equations for the calculation of the variance of concentration are introduced to enhance the prediction of the flow field and scalar quantities. The flow field, the mean concentration and the variance of a flow over a two-dimensional triangular hill, with a finite-size point pollutant source, are calculated by a finite volume code and compared with published experimental results. A modified low Reynolds k-epsilon turbulence model was employed in this work, using the constant of the k-epsilon model C(mu)=0.03 to take into account the inactive atmospheric turbulence. The numerical results for the velocity profiles and the position of the reattachment point are in good agreement with the experimental results. The results for the mean and the variance of the concentration are also in good agreement with experimental results from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

Floquet stability analysis of the flow around an oscillating cylinder

This investigative work is concerned with the flow around a circular cylinder submitted to forced transverse oscillations. The goal is to investigate how the transition to turbulence is initiated in the wake for cases with different Reynolds numbers (Re) and displacement amplitudes (A). For each Re the motion frequency is kept constant, close to the Strouhal number of the flow around a fixed cylinder at the same Re. Stability analysis of two-dimensional periodic flows around a forced-oscillating cylinder is carried out with respect to three-dimensional infinitesimal perturbations. The procedure consists of performing a Floquet type analysis of time-periodic base flows, computed using the spectral/hp element method. With the results of the Floquet calculations, considerations regarding the stability of the system are drawn, and the form of the instability at its onset is obtained. The critical Reynolds number is observed to change with the amplitude of oscillation. With respect to instabilities, unstable modes with the same symmetry as mode A of a fixed cylinder are observed; however, they present different wavelengths. Also, the instabilities observed for the oscillating cylinder are distinctively stronger in the braid shear layers. Other unstable modes similar to mode B are found. Quasi-periodic modes are observed in the 2S wake, and subharmonic mode occurrences are reported in P + S wakes. (C) 2009 Elsevier Ltd. All rights reserved.