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.

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.

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.

Analysis of a boundary protein delimiting chromatin domains in yeast

SUMMARY The expression state of a eukaryotic gene depends in part on its location in the chromosome. This position effect results from the organization of eukaryotic genomes into discrete functional domains, defined by local differences in chromatin structure. The expression of genes within each domain appears to be defined and maintained by the concerted action of regulatory elements such as promoters, enhancers, silencers and locus control regions. Individual domains may be bordered by boundary elements that separate regions of permissive and silent chromatin. When located next to chromosomal elements such as telomeres, genes can be subjected to epigenetic silencing. In yeast, this is mediated by the propagation of the SIR proteins from telomeres towards more centromeric regions. Particular transcription factors can protect downstream genes from silencing when tethered between the gene and the telomere, and they may thus act as chromatin domain boundaries. Here we have studied one of these transcription factors, CTF-1, that binds directly histone H3. A deletion mutagenesis localized the barrier activity to CTF-1 histone-binding domain. A saturating point mutagenesis of this domain identified several amino-acid substitutions that similarly inhibited the boundary and histone-binding activities. Chromatin immunoprecipitation experiments indicated that the barrier protein efficiently prevents the spreading of SIR proteins, and that it separates domains of hypoacetylated and hyperacetylated histones. Together, these results suggest a mechanism by which proteins such as CTF-1 may interact directly with histone H3 to prevent the propagation of a silent chromatin structure, thereby defining boundaries of permissive and silent chromatin domains. RESUME L'expression des gènes eucaryotes dépend en partie de leur localisation sur les chromosomes. Cet effet de position résulte de l'organisation des génomes eucaryotes en domaines fonctionnels, définis par des changements locaux au niveau de la structure de la chromatine. Dans chacun de ces domaines, l'expression des gènes est définie et maintenue par l'action concertée de différents éléments régulateurs tels que les promoteurs, les amplificateurs, les silenceurs et les locus control régions. Ces domaines peuvent être entourés par des éléments barrière, séparant les régions de chromatine répressive des régions permissive pour l'expression des gènes. Lorsqu'ils se situent à proximité d'éléments chromosomiques comme les telomères, les gènes peuvent être réprimés de manière épigénétique. Chez la levure, cette répression est établie par la propagation des protéines SIR depuis les télomères vers les régions centromériques. Certains facteurs de transcription peuvent empêcher la répression d'un gène, lorsqu'ils sont placés entre ce gène et le télomère. Nous avons étudié un de ces facteurs, CTF-1, qui a la particularité de lier directement l'histone H3. La délétion de certaines parties de CTF-1 a permis de déterminer que la région responsable de l'activité barrière correspond au domaine d'interaction avec H3. Plusieurs mutations points effectuées dans ce domaine inhibent à la fois l'activité barrière et la capacité de lier H3. Des expériences d'immuno-précipitation de la chromatine indiquent que la protéine barrière CTF-1 prévient efficacement la propagation des protéines SIR et sépare des domaines contenant des histones hypo-acétylées de ceux constitués d'histones hyper-acétylées. Ces résultats suggèrent que CTF-1 interagit directement avec l'histone H3 pour empêcher la propagation de la chromatine répressive, délimitant ainsi des domaines de chromatine permissive et des domaines de chromatine silencieuse.

MAR elements as tools to increase protein production by CHO cells

One of the major hurdles of isolating stable, inducible or constitutive high-level producer cell lines is the time-consuming selection, analysis and characterization of the numerous clones required to identify one with the desired characteristics. Various boundary elements, matrix attachment regions, and locus control regions were screened for for their ability to augment the expression of heterologous genes in CHO and other cells. The 5'-matrix-attachment region (MAR) of the chicken lysozyme gene was found to significantly increase stable gene expression, in culture dishes and in bioreactors. These MAR elements can be easily combined with various existing expression systems, as they can be added in trans (i.e. on a separate plasmid) in co-transfections with previously constructed expression vectors. Using cell population analysis, we found that the use of the MAR increases the proportion of high-producing CHO cell clones, thus reducing the number of cell lines that need to be screened while increasing maximal productivity. Random cDNA cloning and sequencing indicated that over 12% of the ESTs correspond to the transgene. Thus, productivity is no longer limited by transcriptional events in such MAR-containing cell lines. The identification of small and more convenient active MAR portions will also be summarized. Finally, we will show examples of how MAR elements can be combined with short term expression to increase the simultaneous synthesis of many proteins in parallel by CHO cells. Overall, we conclude that the MAR sequence is a versatile tool to increase protein expression in short and long term production processes.

Sparse space-time boundary element methods for the heat equation

The goal of this work is the efficient solution of the heat equation with Dirichlet or Neumann boundary conditions using the Boundary Elements Method (BEM). Efficiently solving the heat equation is useful, as it is a simple model problem for other types of parabolic problems. In complicated spatial domains as often found in engineering, BEM can be beneficial since only the boundary of the domain has to be discretised. This makes BEM easier than domain methods such as finite elements and finite differences, conventionally combined with time-stepping schemes to solve this problem. The contribution of this work is to further decrease the complexity of solving the heat equation, leading both to speed gains (in CPU time) as well as requiring smaller amounts of memory to solve the same problem. To do this we will combine the complexity gains of boundary reduction by integral equation formulations with a discretisation using wavelet bases. This reduces the total work to O(h

Non-linear boundary element analysis of floor slabs reinforced with rectangular beams

In this work, a numerical model to perform non-linear analysis of building floor structures is proposed. The presented model is derived from the Kirchhoff-s plate bending formulation of the boundary element method (BENI) for zoned domains, in which the plate stiffness is modified by the presence of membrane effects. In this model, no approximation of the generalized forces along the interface is required and the compatibility and equilibrium conditions along interfaces are imposed at the integral equation level. In order to reduce the number of degrees of freedom, the Navier Bernoulli hypothesis is assumed to simplify the strain field for the thin sub-regions (rectangular beams). The non-linear formulation is obtained from the linear formulation by incorporating initial internal force fields, which are approximated by using the well-known cell sub-division. Then, the non-linear solution of algebraic equations is obtained by using the concept of the consistent tangent operator. The Von Mises criterion is adopted to govern the elasto-plastic material behaviour checked at points along the plate thickness and along the rectangular beam element axes. The numerical representations are accurately obtained by either computing analytically the element integrals or performing the numerical integration accurately using an appropriate sub-elementation scheme. (C) 2007 Elsevier Ltd. All rights reserved.

ANALYSIS OF CURVED EDGE PLATES BY THE BOUNDARY-ELEMENT METHOD

The aim of this work is to present a formulation of the boundary element method to analyse elastic and isotropic plates with curved boundaries. In this study the plate boundary is approximated, along each element, by a second degree polynomial relation or by a circular arch, in order to better represent the real boundary. The numerical integration is performed by the self-adaptive coordinate transformation proposed by Telles. The effective shear forces are approximated by concentrated reactions applied at the boundary element nodes, according to the alternative formulation introduced by Paiva. Some examples are presented to demonstrate the better accuracy obtained with the proposed elements.

Non-linear analysis of reinforced concrete plates by the boundary element method and continuum damage mechanics

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

An anisotropic linear elastic boundary element formulation for masonry wall analysis

This work presents an application of a Boundary Element Method (BEM) formulation for anisotropic body analysis using isotropic fundamental solution. The anisotropy is considered by expressing a residual elastic tensor as the difference of the anisotropic and isotropic elastic tensors. Internal variables and cell discretization of the domain are considered. Masonry is a composite material consisting of bricks (masonry units), mortar and the bond between them and it is necessary to take account of anisotropy in this type of structure. The paper presents the formulation, the elastic tensor of the anisotropic medium properties and the algebraic procedure. Two examples are shown to validate the formulation and good agreement was obtained when comparing analytical and numerical results. Two further examples in which masonry walls were simulated, are used to demonstrate that the presented formulation shows close agreement between BE numerical results and different Finite Element (FE) models. © 2012 Elsevier Ltd.

Probabilistic crack growth analyses using a boundary element model: Applications in linear elastic fracture and fatigue problems

This paper addresses the numerical solution of random crack propagation problems using the coupling boundary element method (BEM) and reliability algorithms. Crack propagation phenomenon is efficiently modelled using BEM, due to its mesh reduction features. The BEM model is based on the dual BEM formulation, in which singular and hyper-singular integral equations are adopted to construct the system of algebraic equations. Two reliability algorithms are coupled with BEM model. The first is the well known response surface method, in which local, adaptive polynomial approximations of the mechanical response are constructed in search of the design point. Different experiment designs and adaptive schemes are considered. The alternative approach direct coupling, in which the limit state function remains implicit and its gradients are calculated directly from the numerical mechanical response, is also considered. The performance of both coupling methods is compared in application to some crack propagation problems. The investigation shows that direct coupling scheme converged for all problems studied, irrespective of the problem nonlinearity. The computational cost of direct coupling has shown to be a fraction of the cost of response surface solutions, regardless of experiment design or adaptive scheme considered. (C) 2012 Elsevier Ltd. All rights reserved.

Application of evolutionary algorithms and the boundary element method in the optimization of noise barrier profiles

[EN]This Ph.D. thesis presents a general, robust methodology that may cover any type of 2D acoustic optimization problem. A procedure involving the coupling of Boundary Elements (BE) and Evolutionary Algorithms is proposed for systematic geometric modifications of road barriers that lead to designs with ever-increasing screening performance. Numerical simulations involving single- and multi-objective optimizations of noise barriers of varied nature are included in this document. results disclosed justify the implementation of this methodology by leading to optimal solutions of previously defined topologies that, in general, greatly outperform the acoustic efficiency of classical, widely used barrier designs normally erected near roads.

An efficient ‘a priori’ model reduction for boundary element models

The Boundary Element Method (BEM) is a discretisation technique for solving partial differential equations, which offers, for certain problems, important advantages over domain techniques. Despite the high CPU time reduction that can be achieved, some 3D problems remain today untreatable because the extremely large number of degrees of freedom—dof—involved in the boundary description. Model reduction seems to be an appealing choice for both, accurate and efficient numerical simulations. However, in the BEM the reduction in the number of degrees of freedom does not imply a significant reduction in the CPU time, because in this technique the more important part of the computing time is spent in the construction of the discrete system of equations. In this way, a reduction also in the number of weighting functions, seems to be a key point to render efficient boundary element simulations.