Desenvolvem-se a formulação do Método dos Elementos de Contorno e correspondente algoritmo (para implementação em microcomputador) para a análise de propagação de fraturas em domínios bidimensionais. São utilizados elementos lineares isoparamétricos, tanto para discretizar o contorno quanto para simular a fratura. Os elementos de fratura são descontínuos. A formulação é baseada em equações integrais de tensões e de deslocamentos, onde o termo que considera tensões iniciais concentradas na linha de fratura é formulado a partir da definição de dipolos. O critério adotado é o modelo de fratura coesiva. Os termos singulares e hiper-singulares da formulação são tratados analiticamente e os termos quase-singulares são calculados através de um esquema numérico baseado na utilização de sub-elementos. Os valores dos dipolos são estimados ponto a ponto. Ao longo das fraturas, o valor máximo da tensão normal de tração permite definir novos elementos. As tensões de cisalhamento são removidas para manter a direção principal durante o processo. / The Boundary Element Method formulation and corresponding algorithrn (for microcomputer implementation) are developed for crack growth analysis in two-dimensional domains. Linear isoparametric elements are used to discretize both boundary and crack path. Fracture elements are assumed to be discontinuous. The formulation is based on stress and displacement integral equations, where the term that takes into account initial stresses concentrated along fracture line is formulated from dipoles definition. The coesive fracture rnodel is the criterium adopted. Singular and hipersingular formulation terms are anallitically treated and quasi-singular terms are computed by a numerical scheme based on element subdivision. Dipole values are estirnated point by point. Along fractures, the maximum normal tensile strenght is used to define new elements. Shear stresses are also removed to maintain the principal direction during the process.

Neste trabalho, estudam-se problemas bidimensionais de propagação de ondas acústicas e elásticas, no domínio da freqüência, formulados através do Método dos Elementos de Contorno. A formulação é baseada nas representações integrais das equações diferenciais que governam os fenômenos de propagação de ondas acústicas num meio fluido e de ondas elásticas numa estrutura elástica. Analisa-se também a interação entre o fluido e a estrutura com o uso de sistemas de equações acoplados. As soluções fundamentais utilizadas são expressões exatas e não há necessidade de subdivisão dos domínios em células de integração. São aplicadas técnicas de integração alternativas na escolha das equações algébricas no domínio do fluido, visando a melhora das respostas globais do conjunto. Apresentam-se ainda exemplos numéricos, com o objetivo de possibilitar a modelagem numérica de problemas de acoplamento fluido-estrutura e de radiação e difusão acústica. / In this work, acoustic and elastic wave propagation problems in 2D, in frequency domain, are studied and formulated with the Boundary Element Methods. The formulation is based on the integral representations derived from the differential equations that govern the phenomena of acoustic wave propagation in a fluid medium and elastic wave propagation inside an elastic domain. The fluid-structure interaction is also formulated by coupling appropriately the corresponding systems of equations. The fundamental solutions adopted in this work are conveniently chosen to avoid the mass integral terms in the elastic wave integral representation and the equivalent terms in the acoustic integral equation. Thus, the algebraic representations of both problems are written only in terms of boundary values. Subdivisions of the domain to perform integrals over cells are not required. In an attempt to improve the global answers of the fluid problem, several integration techniques have been experimented to build alternative algebraic matrix equations. Numerical examples are presented in order to shown the accuracy of the studied acoustic radiation and scattering problems and also to verify the proposed fluid-structure coupling.

O objetivo deste trabalho é apresentar uma nova formulação direta do Método dos Elementos de Contorno (M.E.C.) para análise de placas, utilizando a teoria de Kirchhoff, admitindo três parâmetros nodais de deslocamentos para sua representação integral: deslocamento transversal e suas derivadas nas direções normal e tangencial ao contorno. Dois valores nodais são usados para os esforços, momento fletor normal mn e força cortante equivalente Vn. Desta forma são escritas três equações integrais de contorno por nó, obtidas a partir da discretização da placa, segundo a forma usual do método. A vantagem mais perceptível desta formulação é a possibilidade de se fazer a ligação da placa analisada pelo M.E.C. com elementos lineares, representados por três valores nodais de deslocamentos que passam a ser compatibilizados diretamente, para a análise de edifícios. São apresentados exemplos numéricos da formulação e das ligações para comprovação da formulação. / The aim of this work is to present an alternative formulation for plate bending analysis, using Kirchhoff\'s theory, in wich the boundary equation for displacements and its derivative in tangential and normal directions to the boundary for each boundary node are used. The efforts, according to Kirchhoff\'s theory, are the normal bending mn and the equivalent shear force Vn. This formulation is adequate for the analysis of plates coupled with flexible colunms and beams because these structural elements have three nodal displacement values at its nodes. Many examples of single plates and buildings slab are presented using the formulation proposed in this work.

Large silencers used in the power generation industry usually have a very large cross section at the inlet and outlet. Higher-order modes will populate the inlet and outlet even at very low frequencies. Although the silencer itself is often modeled by a three-dimensional analysis tool such as the boundary element method (BEM) or finite element method (FEM), a direct computation of the transmission loss (TL) from the BEM or FEM model can be challenging without incorporating certain forms of modal expansion. A so-called “impedance-to-scattering matrix method” is proposed to extract the modes at the inlet and outlet from the BEM impedance matrix based on the point collocation method. The BEM impedance matrix relates the sound pressures at the inlet and outlet to the corresponding particle velocities, while the scattering matrix relates the modes at the inlet and outlet. Normally there are more boundary elements than the total number of modes at the inlet and outlet, and a least-squares procedure is used to condense the element-based impedance matrix to the mode-based scattering matrix. The TL computation will follow if a certain form of the incident wave is assumed and the outlet is non-reflective. Several commonly used inlet/outlet configurations are considered in this dissertation, which include axisymmetric, non-axisymmetric circular, and rectangular inlet/outlet shapes. In addition to the single inlet and outlet silencers, large multi-inlet and multi-outlet silencers are also investigated. Besides the collocation-based impedance-to-scattering matrix method, an integral-based impedance-to-scattering matrix method based on the reciprocal identity is also proposed for large silencer analysis. Although it may be more time-consuming to perform the additional numerical integration, an integral-based method is free of any uncertainties associated with collocation points. The computational efficiency, accuracy and stability are compared between two proposed methods. One bonus effect of producing the scattering matrix is that it can also be used to combine subsystems in series connection. The Redheffer’s star product is introduced to combine scattering matrices of subsystems. In the design stage, rapid assessment of the silencer performance is always preferred. However, the existing analytical approaches are only suitable for simple dissipative silencers such as straight lined ducts. A two-dimensional first-mode semi-analytical solution is developed to quickly evaluate the performance of tuned dissipative silencers below the cut-off frequency. The semi-analytical solution can also serve as a validation tool for the BEM.

This thesis presents the development of a reliable and efficient technique for the numerical simulation of dynamic soil-structure interaction problems in anisotropic and nonhomogeneous unbounded soils of arbitrary geometry. Such a technique is indispensable in the seismic analysis of large-scale engineering constructions and, to my best knowledge, does not exist at present. The theoretical framework of the research is based on the scaled boundary finite-element method. The following advances are achieved: The scaled boundary finite-element method is extended to simulate the dynamic response of non-homogeneous unbounded domains. The scaled boundary finite element equations in the frequency and time domains are derived for power-type non-homogeneity frequently employed in geotechnical engineering. A high-frequency asymptotic expansion of the dynamic-stiffness matrix is developed. The frequency domain analysis is performed by integrating the scaled boundary finite-element equation in dynamic stiffness. In the time domain, the scaled boundary finite-element equation including convolution integrals is solved for the unit-impulse response at discrete time stations. A Pad?? series solution for the scaled boundary finite-element equation in dynamic stiffness is developed. It converges over the whole frequency range as the order of the approximation increases. The computationally expensive task of numerically integrating the scaled boundary finite-element equation is circumvented. Exploiting the sparsity of the coefficientmatrices in the scaled boundary finite-element equation leads to a significant reduction in computer time and memory requirements for solving large-scale problems. Furthermore, lumped coefficient matrices are obtained by adopting the auss-Lobatto-Legendre shape functions with nodal quadrature, which avoids the eigenvalue problem in determining the asymptotic expansion. A high-order local transmitting boundary constructed from a continued-fraction solution of the dynamic-stiffness matrix is developed. An equation of motion as occurring in standard structural dynamics with symmetric and frequency-independent coefficient matrices is obtained. This transmitting boundary condition can be coupled seamlessly with standard finite elements. Transient responses are evaluated by using a standard timeintegration scheme. The expensive task of evaluating convolution integrals is circumvented. The advances developed in this thesis are applicable in other disciplines of engineering and science to the analysis of scalar and vector waves in unbounded media.

No description available.

In this dissertation, we introduce and analyse a scheme for the fast evaluation of integrals stemming from boundary element methods including discretisations of the classical single and double layer potential operators. Our method is based on the parametrisation of boundary elements in terms of a d-dimensional parameter tuple. We interpret the integral as a real-valued function f depending on d parameters and show that f is smooth in a d-dimensional box. A standard interpolation of f by polynomials leads to a d-dimensional tensor which is given by the values of f at the interpolation points. This tensor may be approximated in a low rank tensor format like the canonical format or the hierarchical format. The tensor approximation has to be done only once and allows us to evaluate interpolants in O(dr(m+1)) operations in the canonical format, or O(dk³ + dk(m + 1)) operations in the hierarchical format, where m denotes the interpolation order and the ranks r, k are small integers. In particular, we apply an efficient black box scheme in the hierarchical tensor format in order to adaptively approximate tensors even in high dimensions d with a prescribed (but heuristic) target accuracy. By means of detailed numerical experiments, we demonstrate that highly accurate integral values can be obtained at very moderate costs.

Efficiency is one of the key issues in numerical simulation of large-scale problems with complex 3-D geometry. Traditional domain based methods, such as finite element methods, may not be suitable for these problems due to, for example, the complexity of mesh generation. The Boundary Element Method (BEM), based on boundary integral formulations (BIE), offers one possible solution to this issue by discretizing only the surface of the domain. However, to date, successful applications of the BEM are mostly limited to linear and continuum problems. The challenges in the extension of the BEM to nonlinear problems or problems with non-continuum boundary conditions (BC) include, but are not limited to, the lack of appropriate BIE and the difficulties in the treatment of the volume integrals that result from the nonlinear terms. In this thesis work, new approaches and techniques based on the BEM have been developed for 3-D nonlinear problems and Stokes problems with slip BC. For nonlinear problems, a major difficulty in applying the BEM is the treatment of the volume integrals in the BIE. An efficient approach, based on the precorrected-FFT technique, is developed to evaluate the volume integrals. In this approach, the 3-D uniform grid constructed initially to accelerate surface integration is used as the baseline mesh to evaluate volume integrals. The cubes enclosing part of the boundary are partitioned using surface panels. No volume discretization of the interior cubes is necessary. This grid is also used to accelerate volume integration. Based on this approach, accelerated BEM solvers for non-homogeneous and nonlinear problems are developed and tested. Good agreement is achieved between simulation results and analytical results. Qualitative comparison is made with current approaches. Stokes problems with slip BC are of particular importance in micro gas flows such as those encountered in MEMS devices. An efficient approach based on the BEM combined with the precorrected-FFT technique has been proposed and various techniques have been developed to solve these problems. As the applications of the developed method, drag forces on oscillating objects immersed in an unbounded slip flow are calculated and validated with either analytic solutions or experimental results.

Since the micro-electro mechanical system microphone, MEMS microphone, has the advantages of superior sound quality, low power consumption, higher temperature resistance and anti-noise ability in used. The researchers therefore have studied the functions of MEMS microphone since 1980s. The MEMS microphones is applied as the part of 3G mobile phone in the market. Though the functions of microphone are improved by manufacturing process technique and new material designed, this study tends to provide a new, low-cost and rapid design idea to gain the performance in chamber of microphone. Taguchi method and BEASY software, which is boundary element method, are combined to evaluate the results of the design in sound field. Taguchi method is a famous method in industrial design to find out relations between system parameters and chamber size. BEASY is a tool for sound field analysis in the research. The result from Taguchi method appears the sound pressure level gain about 2.2 dB to 2.4 dB due to the change of microphone chamber size only. It is also interested in studying the optimization design for position of microphone. It is displayed that the location of port is closer to the boundary of chip will also increase about 0.3 dB to 0.6dB sound pressure level in sound field. The higher frequency of sound source will also create larger sound pressure level at two corners on the port.

This work studied the productivity evaluation and well test analysis of horizontal wells. The major components of this work consist of a 3D coupled reservoir/wellbore model, a productivity evaluation, a deconvolution technique, and a nonlinear regression technique improving horizontal well test interpretation. A 3D coupled reservoir/wellbore model was developed using the boundary element method for realistic description of the performance behavior of horizontal wells. The model is able to flexibly handle multiple types of inner and outer boundary conditions, and can accurately simulate transient tests and long-term production of horizontal wells. Thus, it can serve as a powerful tool in productivity evaluation and analysis of well tests for horizontal wells. Uncertainty of productivity prediction was preliminarily explored. It was demonstrated that the productivity estimates can be distributed in a broad range because of the uncertainties of reservoir/well parameters. A new deconvolution method based on a fast-Fourier-transform algorithm is presented. This new technique can denoise "noisy" pressure and rate data, and can deconvolve pressure drawdown and buildup test data distorted by wellbore storage. For cases with no rate measurements, a "blind" deconvolution method was developed to restore the pressure response free of wellbore storage distortion, and to detect the afterflow/unloading rate function using Fourier analysis of the observed pressure data. This new deconvolution method can unveil the early time behavior of a reservoir system masked by variable-wellbore-storage distortion, and thus provides a powerful tool to improve pressure transient test interpretation. The applicability of the method is demonstrated with a variety of synthetic and actual field cases for both oil and gas wells. A practical nonlinear regression technique for analysis of horizontal well testing is presented. This technique can provide accurate and reliable estimation of well-reservoir parameters if the downhole flow rate data are available. In the situation without flow rate measurement, reasonably reliable parameter estimation can be achieved by using the detected flow rate from blind deconvolution. It has the advantages of eliminating the need for estimation of the wellbore storage coefficient and providing reasonable estimates of effective wellbore length. This technique provides a practical tool for enhancement of horizontal well test interpretation, and its practical significance is illustrated by synthetic and actual field cases.