Este trabalho de pesquisa tem por objetivo apresentar e investigar a viabilidade de um método numérico que contempla o paralelismo no tempo. Este método numérico está associado a problemas de condição inicial e de contorno para equações diferenciais parciais (evolutivas). Diferentemente do método proposto neste trabalho, a maioria dos métodos numéricos associados a equações diferencias parciais evolutivas e tradicionalmente encontrados, contemplam apenas o paralelismo no espaço. Daí, a motivação em realizar o presente trabalho de pesquisa, buscando não somente um método com paralelismo no tempo mas, sobretudo, um método viável do ponto de vista computacional. Para isso, a implementação do esquema numérico proposto está por conta de um algoritmo paralelo escrito na linguagem C e que utiliza a biblioteca MPI. A análise dos resultados obtidos com os testes de desempenho revelam um método numérico escalável e que exige pouco nível de comunicação entre processadores. / This research aims to present and investigate the feasibility of a numerical method that considers the parallelism in time. This numerical method is associated with problems of initial and boundary conditions for (evolutionary) partial differential equations. Unlike the method proposed in this work, most of the numerical methods associated with evolutionary partial diferential equations and traditionally found include only parallelism in space. Hence, the motivation for carrying out the present research work,seeking not only a method with parallelism in time but,above all, a viable method. The implementation of this proposed computational parallel algorithm was written with the language C and uses the MPI library. The results obtained from performance tests reveal a scalable and numerical method that requires little level of communication amount processors

Made available in DSpace on 2016-06-02T19:56:53Z (GMT). No. of bitstreams: 1 DissMF.pdf: 1677378 bytes, checksum: 5f2c8e8509c97ea1775dbbda2508174b (MD5) Previous issue date: 2005-03-23 / Universidade Federal de Minas Gerais / This work studied the horizontal-flow anaerobic immobilized biomass (HAIB) for the treatment of wastewater, developed in the Department of Hydraulics and Sanitation of the School of Engineering of São Carlos - USP. For the simulation and modeling of RAHLF, data of two HAIB were used, and; when the regime is the dynamic, one in bench scale treating wastewater containing poisonous substances (BTEX) and when the regime is the stationary, other in pilot scale treating sewer sanitarium. The study was divided in two parts: stationary regime and transient regime. In the first case the models were analyzed: pseudo - homogeneous with and without axial dispersion and the heterogeneous with axial dispersion; in the second case, the models were analyzed: pseudo - homogeneous and heterogeneous. The models of the first case were resolved respectively for the numeric methods of Runge - Kutta, of the finite differences and of the orthogonal collocation. In those models, when if it despised the axial dispersion the kinetic constant of first order was adjusted. Already in the other two models acted in this same case, they were adjusted the kinetic constant of first order and the coefficient of axial dispersion simultaneously. All of the methods numeric employees in the first case described the tendency of variation of the concentration well along the reactor and inside the particle, making possible the convergence of the solutions. The models without axial dispersion foresaw concentration values in the exit closest of the experimental values than the models with dispersion, and this, for his/her time foresaw closer values to the experimental ones in the intermediate points. In the second case, the models were solved through two methods, the method of the finite differences and the method of the sequence. The methods got to solve the equations that describe the behavior of the reactor satisfactorily and they show that the answer of the exit in the reactor is appropriately made calculations. / Este trabalho estudou o reator anaeróbio horizontal de leito fixo (RAHLF) para o tratamento de águas residuárias, desenvolvido no Departamento de Hidráulica e Saneamento da Escola de Engenharia de São Carlos USP. Para a simulação e modelagem do RAHLF, foram utilizados dados de dois RAHLFs, sendo que; quando o regime é o dinâmico, um em escala de bancada tratando águas residuárias contendo substâncias tóxicas (BTEX) e quando o regime é o estacionário, outro em escala piloto tratando esgoto sanitário. O estudo foi dividido em duas partes: regime estacionário e regime transiente. No primeiro caso foram analisados os modelos: pseudo homogêneo com e sem dispersão axial e o heterogêneo com dispersão axial; no segundo caso, foram analisados os modelos: pseudo homogêneo e heterogêneo. Os modelos do primeiro caso foram resolvidos respectivamente pelos métodos numéricos de Runge Kutta, das diferenças finitas e da colocação ortogonal. Nesses modelos, quando se desprezava a dispersão axial a constante cinética de primeira ordem foi ajustada. Já nos outros dois modelos representados neste mesmo caso, foram ajustados simultaneamente a constante cinética de primeira ordem e o coeficiente de dispersão axial. Todos os métodos numéricos empregados no primeiro caso descreveram bem a tendência de variação da concentração ao longo do reator e no interior da partícula, possibilitando a convergência das soluções. Os modelos sem dispersão axial previram valores de concentração na saída mais próximo dos valores experimentais do que os modelos com dispersão, e este, por sua vez previram valores mais próximos aos experimentais nos pontos intermediários. No segundo caso, os modelos foram solucionados através de dois métodos, o método das diferenças finitas e o método da seqüência. Os métodos conseguiram resolver satisfatoriamente as equações que descrevem o comportamento do reator e mostram que a resposta da saída no reator é adequadamente calculada.

Neste trabalho, é estudada a descrição do fluxo de um gás rarefeito em um duto cilíndrico de comprimento infinito. A formulação matemática do problema está baseada na forma integral de equações cinéticas derivadas da Equação de Boltzmann. Particularmente são estudados os modelos cinéticos conhecidos como BGK e S. Métodos espectrais são propostos para obtenção de soluções, em forma fechada, para quantidades de interesse como o perfil de velocidade do gás, bem como taxas de fluxo. As formulações espectrais são baseadas em duas abordagens: expansão clássica em termos de Polinômios de Legendre e expansão em termos de splines cúbicas de Hermite, neste caso, associada a um esquema de colocação. A implementação das propostas produz resultados computacionais satisfatórios do ponto de vista prático. Para obtenção de resultados com maior precisão, técnicas de tratamento da singularidade do núcleo da equação integral foram introduzidas, resultando em ganho computacional significativo. Finalmente, a proposta de solução espectral para problemas em geometria cilíndrica se mostrou adequada para problemas em que se admite reflexão especular na superfície do cilindro, situação onde outras abordagens clássicas disponíveis na literatura não podem ser utilizadas. / In this work, rarefied gas flows in cylindrical ducts are studied. The mathematical formulation of the problems are based on the integral form of kinetic equations derived from the Boltzmann equation. Particularly, the BGK and S models are studied. Spectral methods are proposed to obtain closed form solutions for quantities of interest as velocity profile of the gas as well as flow rates. The spectral formulations are based on two approaches: classical expansions in terms of Legendre Polynomials and Hermite cubic splines expansions. In this case, associated with a collocation scheme. The approaches provide good computational results, from the practical point of view. On the other hand, for obtaining higher accuracy, some techniques were introduced to deal with the inherent singularity of the integral kernel. In this context, a significant computational gain is achieved. Finally, this spectral approach has shown to be adequate to solve problems where specular reflection is assumed at the surface, in which cases, classical approaches available in the literature can not be used.

O trabalho apresenta uma modificação do modelo do número de onda cumulativo (CW – Cumulative Wavenumber) para atender ao balanço de energia na determinação da transferência radiante de gases participantes não uniformes. De um modo geral, o modelo CW fornece resultados que se mostram precisos para o divergente do fluxo de calor radiante (ou taxa volumétrica de geração de energia) quando comparados à solução benchmark (integração linhapor- linha). No entanto, conforme é apresentado nesta pesquisa, uma hipótese importante feita pelo modelo o torna não conservativo e, não garantindo o balanço de energia radiante quando um meio não uniforme é considerado. Como conseqüência, o fluxo de calor radiante pode apresentar desvios consideráveis em relação à solução exata. O modelo do número de onda cumulativo modificado (CWM) foi desenvolvido de forma a manter o mesmo valor para a taxa volumétrica de geração de energia, mas satisfazendo ao balanço de energia radiante. O modelo CWM é aplicado juntamente com o método das ordenadas discretas para resolver a transferência de calor radiante em sistemas unidimensionais formados por placas planas infinitas e paralelas, contendo camadas uniformes e não uniformes de gases típicos da combustão de metano ou óleo combustível. As paredes são negras para a radiação térmica. A aquisição de dados das linhas de absorção espectrais é feita pela utilização dos bancos de dados HITRAN e HITEMP. São apresentados os bancos de dados, bem como a forma em que é efetuada a extração destes e, posteriormente, a elaboração dos espectros de absorção das espécies químicas. São analisados meios uniformes (isotérmicos e homogêneos) primeiramente, evidenciando que o modelo original (CW) não apresenta problemas nesses casos. Posteriormente, meios não uniformes (nãoisotérmicos e homogêneos; ou não-isotérmicos e não homogêneos) são avaliados. Os resultados para ambos os modelos, CW e CWM, são comparados com a solução linha-por-linha (LBL). Casos com diferentes perfis de temperatura e concentração da espécie química são considerados. A solução a partir da nova metodologia se apresentou mais dispendiosa computacionalmente devido ao acréscimo de novos passos iterativos. Os resultados mostram que o CWM leva a resultados mais precisos para o fluxo de calor radiante, satisfazendo ao balanço de energia com um desvio que se deve principalmente à discretização espacial da equação de transporte radiante. / This work presents a modification of the cumulative wavenumber (CW) method to determine the radiative transfer in non-uniform participating gases to enforce the radiative energy balance to be satisfied. In particular, the CW model leads to results for the divergent radiative heat flux (or radiative volumetric heat source) that proved accurate in comparison to the benchmark solution ( LBL integration) for non-isothermal medium. However, as will be shown in this work, one important assumption of the method prevents it of satisfying the radiative energy balance when a non uniform medium is considered. As consequence, the radiative heat flux can present considerable deviation of the correct solution. The modified cumulative wavenumber (CWM) model was developed to keep the same value of the radiative volumetric heat source, but also to satisfy the radiative energy balance. The CWM model is applied together with the discrete ordinates method to solve the radiation heat transfer in a one-dimentional slab containing a uniform/non-uniform layer of typical gases from the methane or fuel oil combustion. The walls are black for the thermal radiation. The HITRAN and HITEMP are used to extract the spectral lines information that is required for modeling. The database are briefly presented as well as the procedure that is used to extract the data and the spectra drawing. Firstly, uniform media (isothermal and homogeneous) are analyzed proving that the original model (CW) presents accurate results in those cases. Next, non-uniform media (non-isothermal and homogeneous; or non-isothermal and non-homogeneous) are taken into account. The results of both the CW and the CWM modeling are compared to the benchmark line-by-line (LBL) integration. Different temperature and concentration profiles are considered. The solution obtained from the new model presents an increase in the computational time due to the insertion of new iterative loops. The results show that the CWM leads to accurate estimation of both the radiative heat flux and volumetric heat source, satisfying the radiative energy balance with an error that is mostly due to the spatial discretization of the radiative transfer equation.

Este trabalho apresenta algumas contribuições para métodos de controle baseados em dados obtidos em apenas um experimento, a fim de torná-los mais atrativos quanto à aplicação em processos industriais. A partir de dados obtidos em experimentos no processo, os métodos baseados em dados estimam os parâmetros de um controlador de estrutura fixa por meio da minimização do erro entre a saída do sistema real em malha fechada e uma saída desejada, dada por um modelo de referência. O método Virtual Reference Feedback Tuning - VRFT é o método mais expressivo na literatura que estima os parâmetros do controlador usando apenas uma batelada de dados, mas o mesmo apresenta alguns inconvenientes em sua formulação que limitam a sua aplicação. Neste trabalho, o método VRFT é modificado de forma que se obtém um método VRFT flexível, o qual minimiza um critério flexível, pelo qual são estimados tanto os parâmetros do controlador quanto os parâmetros do numerador do modelo de referência; assim, caso a planta que se deseja controlar seja de fase não-mínima, o critério é capaz de estimar esses zeros e os mesmos devem ser incluídos no modelo de referência que será utilizado no projeto do controlador. Além disso, em sistemas com ruído, o método VRFT necessita de uma variável instrumental para que a estimativa dos parâmetros do controlador seja não-polarizada. Para eliminar a necessidade de usar variáveis instrumentais, um novo método de controle baseado em dados é proposto, o qual é descrito sob a ótica de identificação. Este método pode ser visto como a identificação de um sistema, no qual a função de transferência do processo é reparametrizada em função do controlador ideal e do modelo de referência. Além disso, estende-se a teoria de projeto de experimento com solução baseada em restrições LMI para o caso da identificação do controlador ótimo. Todas essas contribuições são ilustradas através de simulações. / This work presents some contributions to data-based control methods where data is obtained in only one experiment, in order to make them more attractive to industrial process applications. Using data from experiments on the process, data-based methods estimate the parameters of fixed structure controller, through the minimization of the error between the closed loop response of the system and the desired response, given by a reference model. The Virtual Reference Feedback Tuning - VRFT method is the most expressive method in the literature that estimates the controller parameters using only one batch of data, but this method presents some inconveniences in its formulation which limit its application. In this work, the VRFT method is modified in order to obtain a flexible VRFT method, which minimizes a flexible criterion, and then obtains the controller parameters together with the parameters related to the reference model numerator; thus, if the plant which we want to control is non-minimum phase, then the criterion is able to estimate these zeros and they need to be included in the reference model that will be used in the control design. Besides, when dealing with noisy systems, the VRFT method needs an instrumental variable so the controller parameters estimate is unbiased. In order to eliminate this necessity, a new data-based control is proposed in this work, which is formulated using identification theory. This method can be seen as the identification of a system, where the process transfer function is reparameterized as a function of the ideal controller and the reference model. Besides, we extend the experimental design theory where the problem is solved using LMI constraints to the case of the optimal controller identification method. All these contributions are illustrated through simulations.

Métodos de controle baseado em dados utilizam dados de operação da planta para minimizar algum critério de desempenho, tipicamente quadrático, sem que seja necessário o conhecimento de um modelo da planta. O Iterative Feedback Tuning (IFT) é um destes métodos, utilizando dados obtidos da operação da planta em malha fechada para gerar uma estimativa não-polarizada do gradiente do critério de desempenho, que é então minimizado através de algum variante do método de Newton. Alguns resultados da literatura tratam da aplicação do IFT em sistemas não-lineares, sempre sob a condição de diferenciabilidade da não-linearidade (linearização em torno do ponto de operação). Um exemplo prático da aplicação do método na sintonia do controlador de um servomotor com folga é o único resultado encontrado. Este trabalho investiga, através de casos simulados e práticos, a aplicabilidade do método IFT em sistemas com não-linearidades não-diferenciáveis comumente encontradas em sistemas mecânicos, como a folga e a zona-morta, e a saturação, presente em todo e qualquer sistema físico. / Data based control methods use data collected from a plant’s “regular” operation to minimize a performance criterion, typically in quadratic form, without any information about a system’s model. Iterative Feedback Tuning (IFT) is one such method that uses closed-loop operating data to generate an unbiased estimate of the cost criterion gradient, which is then minimized with some form of the Newton method. Results in the literature for IFT applied to non-linear systems require that the non-linearity is differentiable (basically, linearizing the system around an operating point). One practical application of the algorithm to tune a controller for a servomotor system with backlash is the only result found. This work investigates, using simulated and practical examples, the suitability of the use of IFT in systems with non-differentiable non-linearities commonly found in mechanical systems, such as backlash and dead-zone, and saturation, found in every physical system.

Engenheiros geotécnicos, defrontam-se frequentemente com solos de baixa capacidade de suporte. Para viabilizar projetos nestes materiais, pode-se utilizar fundações com grandes dimensões e/ou melhorar a propriedades mecânicas destes solos. Fundações de grandes dimensões podem ser extremamente onerosas e gerar impactos ambientais indesejáveis provocados por grandes movimentações de solo. Em contrapartida, o melhoramento das propriedades mecânicas do solo pode ser obtido por meio de tratamento com agentes cimentantes. Esta técnica tem-se mostrado bastante promissora nas diversas subáreas da geotecnia, como por exemplo, leito e subleito de pavimentação, estabilidade de taludes e terrenos para assentamento de fundações à compressão. A utilização desta técnica em fundações escavadas submetidas à tração ainda é incipiente. Até meados dos anos 1950, os métodos de previsão de desempenho consideravam somente o peso do solo contido em uma superfície ruptura somado ao peso da fundação, como contribuindo na capacidade de carga ao arrancamento. Estudos recentes passaram a incluir uma terceira parcela de resistência devido à resistência ao cisalhamento do solo. A partir de então diversos autores passaram a estudar técnicas de melhoramento e reforços dos reaterros. Num primeiro momento, foi estudada a influência das técnicas de compactação do reaterro e a substituição do material por outro com melhores propriedades. Atualmente, os estudos concentram-se no reforço dos reaterros com geossintéticos e adição de agentes cimentícios. Neste contexto, o objetivo deste trabalho é realizar um estudo avaliando os benefícios gerados pelo aumento da capacidade de carga, de fundações escavadas e reaterradas com areia fina e homogênea (proveniente da cidade de Osório/RS), cimentada, avaliando sua influência em termos das propriedades geométricas (diâmetro das fundações, profundidade de assentamento e diâmetro de tratamento) e geotécnicas (coesão, ângulo de atrito e módulo de elasticidade, influenciados pelo teor de cimento). Para tanto foram realizadas provas de carga à tração, variando-se o diâmetro (30 e 45 cm), profundidades de embutimentos (entre 0,5 e 2,0), e teores de cimento (0%, 3% e 7%). Através das provas de carga, verificou-se que o aumento do teor de cimento e da profundidade de assentamento aumentam a capacidade de carga à tração. Os resultados experimentais foram reproduzidos pelo método de elementos finitos por meio de retroanálise. A partir do modelo numérico gerado pela retroanálise, foram executadas duas análises paramétricas. Na primeira foram definidos os parâmetros significativos – coesão e embutimento – e determinada uma equação adimensional para determinação da carga máxima para tratamento de uma camada infinita. Na segunda análise foram definidos ábacos para determinação da perda de carga em função da redução do diâmetro de tratamento, para três níveis de cimentação – baixa, média e alta. / Geotechnical engineers frequently have to deal with soils that have reduced strength. In order to carry out the design of footings in such materials, their base has to be quite large and/or mechanical properties of the soils have to be improved. Shallow foundations of large dimensions can be extremely expensive and generate undesirable environmental impacts due to large soil movements. As an alternative, the improvement of mechanical properties of the local soil can be obtained by treating it with cementitious agents. Such technique has been used with success in several earthworks such as the improvement of base and sub-base of pavements, slope stability and particularly as a soil-cement mixture of a compacted layer over a low bearing capacity soil. The use of such technique in footings subjected to pullout forces is still insipient. Until middle 1950’s, pullout design methodologies of shallow foundations embedded in soil backfills considered only the weight of the soil contained inside a specific failure surface plus the foundation self-weight as the pullout failure load. Further studies included the shear strength of the soil at the failure surface. Since then several authors considered several techniques to improve and/or reinforce the backfills. At first the influence of backfill compaction and material substitution to improve soil properties were considered. Nowadays, studies concentrate in the reinforcement of backfills with geosynthetics and in the addition of cementitious agents. In such context, the main aim of present work is to carry out a research evaluating the benefits towards the improvement of uplift performance of footings embedded in cemented stabilized backfill, analyzing the influence of the geometrical characteristics of the problem (diameter of the footings, depth of embedment and size of the improved area) and geotechnical properties of the backfill (cohesion intercept, friction angle and Young’s modulus), the latter influenced by the amount of cementitious agents inserted on the soil. For doing so, several pullout tests were carried out varying foundation diameter (0.30 and 0.45 m), embedment depth according foundation diameter (ranging from 0.5 e 2.0) and cement contents (0%, 3% and 7%) in the stabilized backfill. Results of pullout tests have shown that the increase of cement content and embedment depth of the backfill increased uplift capacity. Field results were reproduced through a back-analysis using the finite element technique. In the sequence, two distinct parametrical analyses were also carried out. In the first analysis the parameters that show significance were cohesion intercept and embedment depth – a non-dimensional equation determining the failure pullout load (considering infinite horizontal improvement) was determined. In the second analysis distinct abacus were produced to allow considering reduction of horizontal treatment area for three cementation levels – low, medium e high.

Este trabalho apresenta o resultado da aplicação do Método dos Volumes Finitos, adotado pelo software comercial Star-CCM+ na simulação para o regime permanente de escoamentos turbulentos incompressíveis e compressíveis sobre seções de aerofólios. Para o caso incompressível modelam-se seções do aerofólio NACA 0012 com ângulo de ataque zero. Para o caso compressível, uma seção do aerofólio supercrítico OAT15A em pequeno ângulo de ataque é modelada. Os domínios computacionais são discretizados por malhas não estruturadas de células poliédricas ou por malhas estruturadas de geração hiperbólica para diferentes topologias e parâmetros construtivos determinados pela estimativa de grandezas do fenômeno físico, como a altura da primeira camada de células quadrilaterais junto à parede. A qualidade e adequação dessas malhas para as simulações são verificadas por estudo de dependência quanto ao nível de refinamento e também quanto à posição do contorno onde o escoamento é livre no caso de escoamento compressível. Na metodologia de verificação, o índice de convergência de malha GCI e a ordem observada de convergência do método (dos Volumes Finitos) são obtidos para três níveis de refinamento com o propósito de selecionar uma malha de trabalho que concilie precisão e esforço computacional com os recursos disponíveis. As simulações são conduzidas para dois modelos de turbulência – o modelo Spalart-Allmaras e o modelo k-ω/SST. Os resultados obtidos pela aplicação desses modelos são interpretados sob o ponto de vista fenomenológico e comparados com os resultados experimentais disponíveis em literatura. / The Finite Volumes Method adopted by the commercial software Star-CCM+ is applied to the simulation of the steady state regime of incompressible and compressible turbulent flows over selected airfoil’s sections. The physical model used with the incompressible flow case is a NACA 0012 airfoil section at zero angle of attack. The ONERA’s OAT15A supercritical airfoil section at small angle of attack applies to the compressible flow case. The computational domains are discretized by structured and unstructured meshes for different topologies and far field configurations. The structured meshes are of the quadrilateral type with hyperbolic node distribution whilst the unstructured meshes use polyhedral cells. The grids are generated by applying a methodology where estimates of the flow variables are used as input for the grid’s constructive parameters like the near wall cell thickness. Grid dependency studies are carried out in order to verify the grid’s quality and suitability to represent the physical phenomena. The grid’s asymptotic convergence index GCI and its observed order of convergence are evaluated for three refinement levels and far field position for the compressible flow cases. The objective is to select the most suitable grid taking into account the accuracy requirements and the computational resources available. The one-equation Spalart-Allmaras turbulence model and the two-equation k-ω/SST turbulence models are used. The numerical results are discussed from the physical point of view and compared with the experimental ones available in literature.

Este trabalho de pesquisa tem por objetivo apresentar e investigar a viabilidade de um método numérico que contempla o paralelismo no tempo. Este método numérico está associado a problemas de condição inicial e de contorno para equações diferenciais parciais (evolutivas). Diferentemente do método proposto neste trabalho, a maioria dos métodos numéricos associados a equações diferencias parciais evolutivas e tradicionalmente encontrados, contemplam apenas o paralelismo no espaço. Daí, a motivação em realizar o presente trabalho de pesquisa, buscando não somente um método com paralelismo no tempo mas, sobretudo, um método viável do ponto de vista computacional. Para isso, a implementação do esquema numérico proposto está por conta de um algoritmo paralelo escrito na linguagem C e que utiliza a biblioteca MPI. A análise dos resultados obtidos com os testes de desempenho revelam um método numérico escalável e que exige pouco nível de comunicação entre processadores. / This research aims to present and investigate the feasibility of a numerical method that considers the parallelism in time. This numerical method is associated with problems of initial and boundary conditions for (evolutionary) partial differential equations. Unlike the method proposed in this work, most of the numerical methods associated with evolutionary partial diferential equations and traditionally found include only parallelism in space. Hence, the motivation for carrying out the present research work,seeking not only a method with parallelism in time but,above all, a viable method. The implementation of this proposed computational parallel algorithm was written with the language C and uses the MPI library. The results obtained from performance tests reveal a scalable and numerical method that requires little level of communication amount processors

objetivo principal do trabalho é determinar a capacidade última de tração (arrancamento) de uma fundação com solo tratado no topo a partir de uma modelagem numérica com o programa computacional Abaqus a partir dos ensaios feitos em solos tratados. Os modelos numéricos foram feitas a partir dos ensaios ao arrancamento em solos tratados por Ruver, 2011 e Girardello, 2014 para um aterro de areia (proveniente da cidade de Osório/RS) e areia com fibra (monofilamentos de polipropileno). Foi feito um modelo numérico inicial com as propriedades geotécnicas de ensaios feitos em laboratório e campo para depois fazer um ajuste do modelo com os ensaios ao arrancamento. Foi avaliado o modelo numérico de areia e areia-fibra para as tensões iniciais, deslocamentos inciais e a etapa geostatic porque não só é ter bons resultados nas figuras do ajuste, mas se ter um modelo numérico bem feito que apresente um comportamento próximo da realidade numericamente. Foram avaliados os benefícios gerados pelo aumento da capacidade de carga de fundações escavadas e reaterradas com areia fina e homogênea, areia-fibra, avaliando sua influência geotécnica (coesão, ângulo de atrito e módulo de elasticidade) a partir de uma análise paramétrica. Foi feito um método de dimensionamento e estimativa da capacidade ao arrancamento para um aterro com areia e areia-fibra a partir das metodologias internacionais e os modelos numéricos. / The main objective of the work is to determine the ultimate traction (pullout) capacity of a foundation with treated soil at the top from a numerical modeling with the Abaqus computational program from the tests on treated soils. The numerical models were made from the pullout tests on soils treated by Ruver, 2011 and Girardello, 2014 for a sand dump (from the city of Osório / RS) and sand with fiber (polypropylene monofilaments). An initial numerical model was made with the geotechnical properties of laboratory and field tests to later calibrate the model with the pullout tests. It was evaluated the numerical model of sand and sand-fiber for the initial stresses, displacements and the geostatic step because not only is it to have good results in the figures of the calibration, but it is necessary to have a well-designed numerical model that presents behavior close to reality numerically. It was evaluated the benefits generated by the increase of the load capacity of excavated foundations and recovered with fine sand and homogeneous, sand-fiber, evaluating its geotechnical influence (cohesion, friction angle and modulus of elasticity) from a parametric analysis. A method of dimensioning and estimation of the pullout capacity for a backfill with sand and sand-fiber was done from the international methodologies and the numerical models.