Orientador : Prof. Dr. José Antonio Marques Carrer / Coorientador : Prof. Ph.D. Guilherme Bontorin Alves / Tese (doutorado) - Universidade Federal do Paraná, Setor de Tecnologia, Programa de Pós-Graduação em Métodos Numéricos em Engenharia. Defesa: Curitiba, 26/08/2016 / Inclui referências : f. 71-76 / Área de concentração : Mecânica Computacional / Resumo: O estudo em questão é baseado no desenvolvimento da formulação MEC-DT (com DT significando Domínio do Tempo) para a solução do problema de propagação de ondas unidimensionais, utilizando a solução fundamental dependente do tempo. Diferentes condições iniciais e de contorno são testadas na análise de 4 exemplos, utilizando uma barra unidimensional de comprimento L. Pode-se observar que o MEC-DT produz resultados eficazes em termos de precisão numérica e computacional, tendo em vista a exatidão das soluções numéricas quando comparadas com as soluções analíticas correspondentes em cada exemplo abordado. Palavras-chave: Método dos Elementos de Contorno; Ondas Unidimensionais; Condições Iniciais. / Abstract: The study in question is based on the development of DT-BEM (DT meaning Time Domain) formulation to solve the problem of one-dimensional wave propagation using a time-dependent fundamental solution. Di_erent initial and boundary conditions are tested in the analysis of 4 examples, using a one-dimensional bar of length L. It can be observed that the DT-BEM produces e_ective results in terms of numerical and computational precision, in view of the accuracy of the numerical solutions when compared with the corresponding analytical solutions in each example discussed. Key-words: Boundary Element Method; Wave Unidimensional; Initial Conditions.

Submitted by maria angelica Varella (angelica@sibi.ufrj.br) on 2018-01-31T14:04:10Z No. of bitstreams: 1 149229.pdf: 686174 bytes, checksum: f170031e3e4162fda482686a266b9293 (MD5) / Made available in DSpace on 2018-01-31T14:04:10Z (GMT). No. of bitstreams: 1 149229.pdf: 686174 bytes, checksum: f170031e3e4162fda482686a266b9293 (MD5) Previous issue date: 1978-11 / FINEP / CNPq / Neste trabalho de problema de contorno livre, é analisado um modelo transiente no sentido de representar os perfis de temperatura nas fases líquida e sólida, durante a fusão térmica. Os resultados experimentais obtidos em câmara especialmente: construída são confrontados com o modelo teórico condutivo, unidirecional, a fim de se estabelecer o comportamento dinâmico da frente de fusão. É também estudada a eficiência do processo térmico e sua maximização como uma função das propriedades físicas, incluindo o calor latente. / In this work a transient model of a free boundary value problem is analysed in arder to represent the temperature profile in liquid and solid phases during thermal fusion. The experimental results which were obtained in a special probe are compared with the results of the conductive unidirectional model in arder to determine the dynamical behavoir of the fusion front. The efficiency of the thermal process is also anklised with respect to its maximizing conditions as a function of the metal's physical properties, mainly its latent heat.

Made available in DSpace on 2015-05-08T14:59:44Z (GMT). No. of bitstreams: 1 Arquivototal.pdf: 3772690 bytes, checksum: e032df086ce22dd8ceaa668979d9ed25 (MD5) Previous issue date: 2012-10-30 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / This work deals with the analysis of soil-structure interaction modeling of pipeline problems in static behavior using the coupling between Finite Element Method (FEM) and Boundary Element Method (BEM).The representation of the pipe is made by MEF using two finite elements in the cylindrical panel: The first, formulated from the classical theory of cylindrical shells by assumptions of Donnel and initially proposed by Djoudi & Bahai, using the Assumed Strain Based Model. The second formulated from the theory of equivalent discrete layers (Layerwise Theory), proposed by J. N. Reddy. The soil is represented by elastic continum infinite or semi-infinite and modeled using boundary elements with special curved surface, associated with cylindrical panel, used to represent the soil-structure interaction within the soil, especially at the contact surface with the pipe. Besides the mathematical representation of the model, this work propose a library DLL (Dynamic Link Library) written using the object-oriented programming (OOP) for processing soil structure interaction problems and a user-friendly environment built using the ObjectARX library for customizing AutoCAD. / Neste trabalho aborda-se a análise da interação solo-estrutura para modelagem de problemas de dutos em regime estático através do acoplamento entre o método dos Elementos Finitos (MEF) e o método dos elementos de contorno (MEC). A representação do duto homogêneo ou laminado é feita pelo MEF utilizando-se dois elementos finitos na forma de um painel cilíndrico: O primeiro foi desenvolvido a partir do modelo isotrópico homogêneo proposto inicialmente por Djoudi e Bahai, que utiliza a filosofia do modelo de deformações assumidas (Assumed Strain Based Model). O segundo foi proposto por Reddy e é baseado na teoria de camadas equivalentes discretas (Layerwise Theory). O solo é admitido ser um contínuo elástico infinito ou semi-infinito e modelado pelo MEC onde elementos de contorno especiais são propostos com superfície curva ou bordos curvos para representar a interação solo-estrutura principalmente nas superfícies de contato com o duto. Além da representação matemática do modelo da interação solo-estrutura, é proposta uma biblioteca DLL (Dynamic Link Library) escrita utilizando-se a filosofia de programação orientada a objetos (POO) para o processamento de problemas de interação solo estrutura e um ambiente amigável ao usuário construído a partir da biblioteca ObjectARX para customização do AutoCAD.

Made available in DSpace on 2015-05-14T12:09:35Z (GMT). No. of bitstreams: 1 arquivototal.pdf: 9782170 bytes, checksum: e9aff4656bf80c325fe0eccf48ab0a1e (MD5) Previous issue date: 2014-04-03 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / Roads performed an important role in the design process of city, by facilitating long distance journeys. When necessary, contours were constructed that supposedly facilitate connectivity between the city and the highway. However, when absorbed into the urban fabric of cities, these stretches of road contours develop individual dynamics with intra-urban space in which it is inserted, with the change of accessibility, morphology and the use and occupation of urban land. These interfaces occur during the growth process of the city and during the subsequent stages. Within this context, this master's research aims to analyze the dynamics between the implementation of BR-230 road contour and the process of production of the urban space of the city of João Pessoa, Paraíba, between the years 1963 and 2013, in terms of morphology, use and occupation of urban land. To this end, we employed the method of evolutionary logic of the urban fabric of Panerai, which divides the urban evolution in three stages: overcoming limits, growth and combination and conflict. The method was associated with analytical tools such as Space Syntax and the use of Geographic Information System software. The results show that the highway boosted the process of spreading the city to the south, which produced a sprawling urban area, this being the fastest growing sector in that period, where they were also able to identify the pattern of horizontal and fragmented spatial growth, changes in accessibility based on morphology and pattern of use and occupation of spaces bordering the highway / As rodovias desempenham função importante no processo de estruturação das cidades, por facilitar os deslocamentos interurbanos. Quando necessário, eram construídos contornos que supostamente facilitam a conectividade entre a cidade e a rodovia. No entanto, quando absorvidas pelo tecido urbano das cidades, estes trechos de contornos rodoviários desenvolvem dinâmicas particulares com o espaço intraurbano no qual está inserido, com a alteração da acessibilidade, da morfologia e do uso e ocupação da terra urbana. Estas interfaces ocorrem durante o processo de crescimento da cidade e durante os estágios subsequentes. Com base nesse contexto, esta pesquisa de mestrado tem por objetivo analisar as dinâmicas entre a implantação do contorno rodoviário da BR-230 e o processo de produção do espaço intraurbano da cidade de João Pessoa, na Paraíba, entre os anos de 1963 e 2013, em termos de morfologia, uso e ocupação do solo urbano. Para tal, foi empregado o método da lógica evolutiva do tecido urbano de Panerai, que divide a evolução urbana em três estágios: superação de limites, crescimento e combinação e conflitos. O método foi associado a ferramentas analíticas, como a sintaxe espacial e o uso de softwares de Sistema de Informação Geográfica. Os resultados mostram que a rodovia impulsionou o processo de espraiamento da cidade para o sul, que produziu uma malha urbana tentacular, sendo esse setor o que mais cresceu naquele período, onde também se puderam identificar o padrão de crescimento espacial horizontal e fragmentado, alterações na acessibilidade com base na morfologia e no padrão de uso e ocupação dos espaços lindeiros à rodovia

Resumo: Este trabalho trata da inclusão de pilares numa formulação do Método dos Elementos de Contorno para análise linear de pavimentos de edifícios. Inicialmente é apresentada a formulação para análise do pavimento sujeito à flexão simples, onde o pavimento é modelado por uma placa composta por sub-regiões de diferentes rigidezes, sendo cada sub-região a representação de uma viga ou laje. Nesse modelo as forças de superfície são eliminadas das interfaces e todas as sub-regiões são representadas por suas superfícies médias. Então, a fim de reduzir o número de graus de liberdade do problema, aproximam-se os deslocamentos ao longo da largura da viga, obtendo-se um novo modelo em que as variáveis ficam definidas ao longo dos eixos das vigas ao invés de seus contornos. Em seguida, introduz-se o efeito de membrana nessa formulação para que se possa fazer a análise do pavimento sujeito à flexão composta, isto é, representando todas as sub-regiões por uma mesma superfície de referência. Nesse modelo considera-se aproximação linear para as forças de membrana ao longo da largura das vigas. Finalmente, estendem-se essas formulações para análise de pavimentos de edifícios com pilares, através da inclusão de pontos do domínio com carregamentos prescritos, que são escritos em termos das reações do pilar sobre a placa. Os esforços de flexão do pilar sobre a placa são escritos em termos da tensão uniformemente distribuída sobre a seção transversal do pilar e os esforços de membrana são escritos em termos de linhas de carga atuando na seção transversal do pilar. / Abstract: This work deals with the inclusion of columns in a formulation of the Boundary Element Method (BEM) for linear analysis of building floors structures. Initially it is presented the BEM formulation for analyzing building floors subjected to simple bending, where the floor is modeled by a plate divided into sub-regions of different rigidities, being each sub-region the representation of either a beam or slab. In this model the tractions are eliminated on the interfaces and all the sub-regions are represented by their middle surfaces. Then, in order to reduce the number of degrees of freedom, some approximations for the displacements are adopted along the beam width, leading to a model where the values are defined on the beam axis, instead of its boundaries. Then the membrane effects are introduced into this formulation in order to represent all sub-regions by a same reference surface. In this model are assumed linear approximation for the in-plane tractions along the beam cross section. Finally, those formulations are extended in order to consider columns in the plate domain. This is done by considering domain points with prescribed loadings which are written in terms of the column reactions on the plate. The columns bending reactions are written in terms of a stress uniformly distributed over the column cross section and the membrane forces are written in terms of loading lines acting in the column cross section. / Orientador: Gabriela Rezende Fernandes / Coorientador: Rogério de Oliveira Rodrigues / Banca: Luiz Carlos Facundo Sanches / Banca: Carlos Humberto Martins / Mestre

Orientador: Francisco José Teixeira-Neto / Resumo: Introdução e objetivos: Mensurações do débito cardíaco (DC) pela técnica termodiluição transpulmonar (DCTP) e pela análise de contorno de pulso com calibração pela técnica transpulmonar (DCACP) são alternativas menos invasivas em comparação ao DC fornecido pela técnica de termodiluição de artéria pulmonar (DCP). Entretanto, instabilidades hemodinâmicas podem alterar o desempenho destes métodos. A Fase I do estudo objetivou avaliar se a utilização de 10 mL de indicador térmico (solução fisiológica a ≤ 5oC) para o DCTP (artéria femoral) promoveria melhor concordância e habilidade em detectar alterações no DCP em comparação a 5 mL de indicador. Na Fase II, objetivou-se verificar se alterações na resistência vascular sistêmica (RVS) influenciariam a concordância e a habilidade em detectar tendências entre o DCACP e o DCP. Métodos: Em 8 cães adultos (20,8–31,5 kg), mensurações simultâneas em triplicata do DCTP e DCP foram obtidas utilizando 5 e 10 mL de indicador térmico durante anestesia com isoflurano associado ou não com a infusão contínua intravenosa de remifentanil (0,3 e 0,6 μg/kg/min) ou de dobutamina (2,5 e 5,0 μg/kg/min) (Fase I). Durante a Fase II, o DCACP e o DCP foram mensurados simultaneamente (triplicata) antes e durante alterações na RVS induzidas pela infusão contínua de fenilefrina (1,0 μg/kg/min) ou de nitroprussiato (1,0 μg/kg/min). A acurácia e a precisão da concordância entre métodos foram estudadas pela análise de Bland-Altman para medidas múltiplas (Fase I) ... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: Background and objectives: Cardiac output (CO) measurements by transpulmonary thermodilution (TPTDCO) and by pulse contour analysis calibrated with transpulmonary thermodilution (PCACO) are less invasive alternatives to pulmonary artery thermodilution (PATDCO). However, hemodynamic instability could affect the performance of these methods. The objective of Phase I of the study was to determine if the use of 10 mL of thermal indicator (physiological saline at ≤ 5oC) for TPTDCO (measured in the femoral artery) would improve the agreement and trending ability with PATDCO in comparison to 5 mL of indicator. During Phase II, the aim was to verify if changes in systemic vascular resistance (SVR) would alter the agreement and trending ability between PCACO and PATDCO. Methods: In eight adult dogs (20.8–31.5 kg), simultaneous TPTDCO and PATDCO measurements (averaged from 3 repetitions) using 5 and 10 mL of thermal indicator were obtained during isoflurane anesthesia combined or not with intravenous remifentanil (0.3 e 0.6 μg/kg/min) or dobutamine (2.5 e 5.0 μg/kg/min) (Phase-1). During Phase-2, triplicate PCACO and PATDCO measurements were recorded before and during phenylephrine (1.0 μg/kg/min) or nitroprusside (1.0 μg/kg/min) induced changes in SVR. The accuracy and precision of agreement was evaluated by the Bland-Altman method for multiple measurements (Phase I) and for single measurements per subject (Phase 2). The ability of the test methods (PCACO and TPTDCO) to detect changes... (Complete abstract click electronic access below) / Doutor

Nessa tese, mostramos que a Hamiltoniana H e o operador helicidade de uma partícula de Dirac que se movimenta em duas dimensões na presença de um tubo de fluxo magnético infinitamente fino na origem admitem, cada um, uma família de quatro parâmetros de extensões auto-adjuntas. Para cada extensão correspondem condições de contorno a serem satisfeitas pelas auto-fuções na origem. Apesar dos operadores H e formalmente comutarem antes da especificação das condições de contorno, para garantirmos a conservação da helicidade, não é suficiente obtermos as mesmas condições de contorno para ambos os operadores, ou seja, não é suficiente a determinação de um domínio comum a ambos. Mostramos que, para certas relações entre os parâmetros das extensões satisfeitas, é possível a determinação dos domínios mais gerais onde ambos os operadores H e são auto-adjuntos e onde a helicidade é conservada, simultaneamente com a preservação da simetria de Aharonov-Bohm ( + 1), onde é o fluxo magnético em unidades naturais. Nossos resultados implicam que, nem a conservação da helicidade nem a simetria de Aharonov-Bohn, resolvem o problema da escolha da condição de contorno fisicamente correta. / We show that both the Hamiltonian H and the helicity operator of a Dirac particle moving in two dimension in the presence of an infinitely thin magnetic flux tube admit each a four- parameter family of self-adjoint extensions. Each extension is in one-to-one correspondence with the boundary conditions (BC\'s) to be satisfied by the eigenfunctions at the origin. Althou- gh the actions af these two operators commute before specification of boundary conditions, to ensure helicity conservation it is not sufficient to take the same BC\'s for both operators. We show that, given certain relations between the parameters of the extensions it is possible to write down the most general domain where both operators H and are self-adjoint with heli- city conservation and also Aharonov-Bohm symmetry ( + 1) preserved, where is the magnetic flux in natural units. The continuity of the dynamics is also obtained. Our results im- ply that neither helicity conservation nor Aharonov-Bohm symmetry by themselves solves the problem of choosing the \"physical \"boundary conditions for this system.

O método adjunto tem sido extensivamente utilizado como ferramenta de síntese no projeto de aeronaves por permitir que se obtenham sensibilidades de distintas medidas de mérito com relação a parâmetros que controlam a geometria de superfícies aerodinâmicas. O presente trabalho visa uma ampliação das aplicações da formulação contínua do método, ao utilizar propriedades físicas do escoamento nas fronteiras permeáveis do domínio computacional como parâmetros de controle de uma particular medida de mérito. Desse modo é possível, entre muitas possibilidades, determinar a sensibilidade de integrais como sustentação ou arrasto de uma aeronave com relação às condições de cruzeiro, por exemplo. Mais do que isso, essa informação pode ser obtida com a mesma solução adjunta computada para realizar otimização de forma. Vale destacar, ainda, que para que se consiga obter essa informação a partir das equações adjuntas, é necessário que se implemente condições de contorno baseadas em equações diferenciais características, resolvendo o problema de Riemann completo nas fronteiras do domínio. A implementação das usuais condições de contorno homogêneas, vastamente difundidas na literatura, resultaria em gradientes nulos. Esta nova abordagem do método é então aplicada a escoamentos modelados pelas equações de Euler 2-D compressíveis em estado estacionário. Ambos os problemas, físico e adjunto, são resolvidos numericamente com um código computacional que utiliza o método dos volumes finitos com segunda ordem de precisão no espaço e discretização centrada com dissipação artificial. As soluções estacionárias são obtidas ao se postular um termo tempo-dependente e integra-lo com um esquema Runge-Kutta de 5 passos e 2a ordem de precisão. As simulações são realizadas em malhas não-estruturadas formadas por elementos triangulares em 4 geometrias distintas: um bocal divergente, um perfil diamante, um aerofólio simétrico (NACA 0012) e o outro assimétrico (RAE 2822). Os gradientes adjuntos são então validados por meio da comparação com os obtidos pelo método de diferenças finitas nos regimes de escoamento subsônico, supersônico e transônico. / The adjoint method has been extensively used as an aircraft design tool, since it enables one to obtain sensitivities of many different mesures of merit with respect to parameters that control the aerodynamic surface geometry. This works aims to open up the possibilities of the method\'s applications by using flow physical properties at the permeable boundaries of the computational domain as control parameters of a particular measure of merit. This way it is possible, among many possibilities, to compute lift or drag sensitivities of an aircraft with respect to cruise conditions, for instance. Moreover, this information can be obtained with the same adjoint solution used to perform shape optimization. It is also worth noting that in order to obtain this information from the adjoint equations it is necessary to implement characteristics-based boundary conditions, resolving the complete Riemann problem at the boundaries of the computational domain. The use of the traditional homogeneous boundary conditions, widely spread in the literature, would lead the gradient to vanish. This new approach of the method is, then, applied to flows modeled by the 2-D steady state compressible Euler equations. Both, physical and adjoint problems are numerically solved with a computational code that makes use of a 2nd order finite volume method and central differences with artifficial dissipation. The steady solutions are obtained by postulating a time-dependent term and integrating it with a 5-stage 2nd order Runge-Kutta scheme. The simulations are performed on unstructured triangular meshes to 4 different geometries: a divergent nozzle, a diamond profile, a symmetric airfoil (NACA 0012) and a assymmetric airfoil (RAE 2822). The adjoint gradients are then validated by comparison with those obtained by finite differences method in subsonic, supersonic and transonic flow regimes.

A busca por novos modelos matemáticos e técnicas inovadoras para análises numéricas tem sido tema de muitas pesquisas. Em análises de modelos que possuem domínios infinitos e semi-infinitos, o Método dos Elementos de Contorno (MEC) sobressai-se como uma das mais eficientes ferramentas numéricas. Por outro lado, em análises não-lineares o MEC requer a avaliação de integrais de domínio, diminuindo as vantagens de uma discretização apenas do contorno do modelo analisado. Neste trabalho apresenta-se uma técnica inovadora que trata as integrais de domínio, não adequadas para uma representação pura do contorno, em análises de modelos com materiais viscoplásticos. Na abordagem proposta, utiliza-se um novo algoritmo de visualização proposto por Noronha & Pereira para detectar as regiões de plastificação automaticamente. Este procedimento de detecção é realizado de forma incremental por meio de predições (gradiente como direção de busca) e iterações (Newton-Raphson). Uma vez que as regiões sejam obtidas, torna-se possível transformar as integrais de domínio em integrais de contorno de forma direta. Obtém-se assim uma abordagem baseada apenas na discretização do contorno dos modelos, mantendo uma das principais vantagens da utilização do MEC. Foram realizados neste trabalho alguns exemplos numéricos que apresentaram excelentes resultados em comparação com o Método dos Elementos Finitos (MEF) e com resultados encontrados na literatura. / The search for new mathematical models and innovative techniques for numerical analyses has been subject of many research studies. In analysis of models with infinite and semi-infinite domains, the Boundary Element Method (BEM) has been proved to be one of the most efficient numerical tools. On the other hand, in nonlinear analyses the BEM requires the evaluation of domain integrals, diminishing the advantages of a discretization only of the boundary of the model. This work presents an innovative technique that treats the domain integrals, not suitable for pure boundary representations, in analyses of models with viscoplastic materials. The proposed approach is based on a new post-processing algorithm developed by Noronha & Pereira to detect the plastic regions automatically. The detection procedure herein proposed is an incremental technique that uses prediction (along the gradient direction) and iteration (Newton-Raphson) loops. Once the plastic regions are found, it becomes possible to transform the domain integrals in boundary integrals in a straightforward manner. The proposed approach results in a pure boundary discretization, preserving the main advantage of the BEM. The numerical examples presented in this work are in good agreement with the Finite Element Method (FEM) and with results found in the literature.

Um dos principais problemas que o Método de Elementos de Contorno (MEC) apresenta encontra-se na avaliação de integrais singulares e quase-singulares que envolvem as soluções fundamentais de Kelvin em deslocamento e força. O processo de integração numérica em MEC tem sido o objetivo de inúmeras pesquisas, pois dele depende a qualidade das respostas quando se deseja obter uma excelente precisão numérica em uma análise. Este trabalho apresenta uma nova proposta de integração numérica para análise tridimensional com MEC. Esta técnica possui três características importantes. A primeira é a determinação da parcela efetiva de singularidade que ocorre na função raio, distância entre o ponto fonte e o elemento de contorno bidimensional. A correta obtenção desta parcela permite representar sem aproximações o comportamento da singularidade da função raio, que é a verdadeira fonte de singularidade e quase-singularidade nas soluções fundamentais. A segunda característica da técnica proposta é que ela baseia-se em um Método Semi-Analítico de avaliação de integrais, onde, para cada parcela efetiva de singularidade, utiliza-se uma quadratura numérica cujos pesos específicos são calculados analiticamente. A terceira característica da técnica proposta é apresentar um tratamento unificado para todos os tipos de integrais singulares, quasesingulares e regulares. Esta técnica foi implementada na plataforma computacional desenvolvida pelo grupo GoBEM, utilizando o conceito de Programação Orientada a Objetos através da Linguagem de programação Java. Com a implementação da nova técnica de integração na plataforma computacional torna-se possível realizar o desenvolvimento de vários tipos de pesquisa sobre análise tridimensional com o MEC como, por exemplo: visualização de isosuperfícies em análise tridimensional sem discretização do domínio, automatização do cálculo elasto-plástico, modelagem de problemas geotécnicos de forma precisa, etc. Para a validação da técnica proposta três procedimentos foram considerados: análise da eficiência da parcela efetiva de singularidade, testes de convergência da integração numérica específica e exemplos numéricos utilizando o MEC em problemas de engenharia. / One of the main problems with the Boundary Elements Method (BEM) is the evaluation of singular and quasi-singular integrals due to the Kelvin\'s Fundamental Solutions in displacement and traction. Today there is an increasing body of research work that focuses on numerical evaluation of BEM integrals, for this is a crucial issue in order to achieve highly accurate results. This work presents an innovative numerical integration procedure for threedimensional analysis with BEMs. The proposed technique encompasses three important features. First, it corresponds to an accurate representation of the effective term of singularity in the radius function, which measures the distance between the source point and a twodimensional boundary element. The correct evaluation of this term leads without approximation to the actual singularity behavior of the radius function, which is the true source of singularity and quasi-singularity in the fundamental solutions. Second, the proposed technique is based on a semi-analytical procedure, i.e, for each singularity effective term it uses a quadrature scheme with specific weights analytically evaluated. Last, the proposed technique represents a unified procedure to singular, quasi-singular and regular integrals. This technique was implemented in the computational platform developed by the group GoBEM, using Object Oriented Programming and the Java programming language. The implementation of the proposed technique into this computational plataform opens new possibilities for future researches on three-dimensional BEM, e.g.: visualization of isosurfaces in three-dimensional analysis without any domain discretization, automatic elasto-plastic analysis, accurate modeling of geomechanical problems, etc. Validation of the proposed technique was carried out using three procedures: efficiency analysis of the singularity effective term, convergence tests for the specific numerical integration and numerical examples using BEM in engineering problems.