Plane strain transient finite thermomechanical deformations of heat-conducting particulate composites comprised of circular tungsten particulates in nickel-iron matrix are analyzed using the finite element method to delineate the initiation and propagation of brittle/ductile failures by the nodal release technique. Each constituent and composites are modeled as strain hardening, strain-rate-hardening and thermally softening microporous materials. Values of material parameters of composites are derived by analyzing deformations of a representative volume element whose minimum dimensions are determined through numerical experiments. These values are found to be independent of sizes and random distributions of particulates, and are close to those obtained from either the rule of mixtures or micromechanics models. Brittle and ductile failures of composites are first studied by homogenizing their material properties; subsequently their ductile failure is analyzed by considering the microstructure. It is found that the continuously varying volume fraction of tungsten particulates strongly influences when and where adiabatic shear bands (ASB) initiate and their paths. Furthermore, an ASB initiates sooner in the composite than in either one of its constituents. We have studied the initiation and propagation of a brittle crack in a precracked plate deformed in plane strain tension, and a ductile crack in an infinitely long thin plate with a rather strong defect at its center and deformed in shear. The crack may propagate from the tungsten-rich region to nickel-iron-rich region or vice-a-versa. It is found that at the nominal strain-rate of 2000/s the brittle crack speed approaches Rayleigh's wave speed in the tungsten-plate, the nickel-iron-plate shatters after a small extension of the crack, and the composite plate does not shatter; the minimum nominal strain-rate for the nickel-iron-plate to shatter is 1130/s. The ductile crack speed from tungsten-rich to tungsten-poor regions is nearly one-tenth of that in the two homogeneous plates. The maximum speed of a ductile crack in tungsten and nickel-iron is found to be about 1.5 km/s. Meso and multiscale analyses have revealed that microstructural details strongly influence when and where ASBs initiate and their paths. ASB initiation criteria for particulate composites and their homogenized counterparts are different. / Ph. D.

The accurate and efficient numerical simulation of geophysical flows is of great interest in numerical weather prediction and climate modeling as well as in numerous critical areas and industries, such as agriculture, construction, tourism, transportation, weather-related disaster management, and sustainable energy technologies. Oceanic and atmospheric flows display an enormous range of temporal and spatial scales, from seconds to decades and from centimeters to thousands of kilometers, respectively. Scale interactions, both spatial and temporal, are the dominant feature of all aspects of general circulation models in geophysical fluid dynamics. In this thesis, to decrease the cost for these geophysical flow computations, several types of multiscale methods were systematically developed and tested for a variety of physical settings including barotropic and stratified wind-driven large scale ocean circulation models, decaying and forced two-dimensional turbulence simulations, as well as several benchmark incompressible flow problems in two and three dimensions. The new models proposed here are based on two classes of modern multiscale methods: (i) interpolation based approaches in the context of the multigrid/multiresolution methodologies, and (ii) deconvolution based spatial filtering approaches in the context of large eddy simulation techniques. In the first case, we developed a coarse-grid projection method that uses simple interpolation schemes to go between the two components of the problem, in which the solution algorithms have different levels of complexity. In the second case, the use of approximate deconvolution closure modeling strategies was implemented for large eddy simulations of large-scale turbulent geophysical flows. The numerical assessment of these approaches showed that both the coarse-grid projection and approximate deconvolution methods could represent viable tools for computing more realistic turbulent geophysical flows that provide significant increases in accuracy and computational efficiency over conventional methods. / Ph. D.

The goal of this research is developing a computational framework to study mechanical fatigue and fracture at different length scales for a broad range of materials. The developed multiscale framework is utilized to study the details of fracture and fatigue for the rolling contact in rails, additively manufactured alloys, and bio-inspired hierarchical materials. Rolling contact fatigue (RCF) is a major source of failure and a dominant cause of maintenance and replacements in many railways around the world. The highly-localized stress in a relatively small contact area at the wheel-rail interface promotes micro-crack initiation and propagation near the surface of the rail. 2D and 3D microstructural-based computational frameworks are developed for studying the rolling contact fatigue in rail materials. The method can predict RCF life and simulate crack initiation sites under various conditions. The results obtained from studying RCF behavior in different conditions will help better maintenance of the railways and increase the safety of trains. The developed framework is employed to study the fracture and fatigue behavior in 3D printed metallic alloys fabricated by selective laser melting (SLM) method. SLM method as a part of metal additive manufacturing (AM) technologies is revolutionizing the manufacturing sector and is being utilized across a diverse array of industries, including biomedical, automotive, aerospace, energy, consumer goods, and many others. Since experiments on 3D printed alloys are considerably time-consuming and expensive, computational analysis is a proper alternative to reduce cost and time. In this research, a computational framework is developed to study fracture and fatigue in different scales in 3D printed alloys fabricated by the SLM method. Our method for studying the fatigue at the microstructural level of 3D printed alloys is pioneering with no similar work being available in the literature. Our studies can be used as a first step toward establishing comprehensive numerical frameworks to investigate fracture and fatigue behavior of 3D metallic devices with complex geometries, fabricated by 3D printing. Composite materials are fabricated by combining the attractive mechanical properties of materials into one system. A combination of materials with different mechanical properties, size, geometry, and order of different phases can lead to fabricating a new material with a wide range of properties. A fundamental problem in engineering is how to find the design that exhibits the best combination of these properties. Biological composites like bone, nacre, and teeth attracted much attention among the researchers. These materials are constructed from simple building blocks and show an uncommon combination of high strength and toughness. By inspiring from simple building blocks in bio-inspired materials, we have simulated fracture behavior of a pre-designed composite material consisting of soft and stiff building blocks. The results show a better performance of bio-inspired composites compared to their building blocks. Furthermore, an optimization methodology is implemented into the designing the bio-inspired composites for the first time, which enables us to perform the bio-inspired material design with the target of finding the most efficient geometries that can resist defects in their structure. This study can be used as an effective reference for creating damage-tolerant structures with improved mechanical behavior. / Doctor of Philosophy / The goal of this research is developing a multiscale framework to study the details of fracture and fatigue for the rolling contact in rails, additively manufactured alloys, and bio-inspired hierarchical materials. Rolling contact fatigue (RCF) is a major source of failure and a dominant cause of maintenance and replacements in many railways around the world. Different computational models are developed for studying rolling contact fatigue in rail materials. The method can predict RCF life and simulate crack initiation sites under various conditions and the results will help better maintenance of the railways and increase the safety of trains. The developed model is employed to study the fracture and fatigue behavior in 3D printed metals created by the selective laser melting (SLM) method. SLM method as a part of metal additive manufacturing (AM) technologies is revolutionizing industries including biomedical, automotive, aerospace, energy, and many others. Since experiments on 3D printed metals are considerably time-consuming and expensive, computational analysis is a proper alternative to reduce cost and time. Our method for studying the fatigue at the microstructural level of 3D printed alloys can help to create more fatigue and fracture resistant materials. In the last section, we have studied fracture behavior in bio-inspired materials. A fundamental problem in engineering is how to find the design that exhibits the best combination of mechanical properties. Biological materials like bone, nacre, and teeth are constructed from simple building blocks and show a surprising combination of high strength and toughness. By inspiring from these materials, we have simulated fracture behavior of a pre-designed composite material consisting of soft and stiff building blocks. The results show a better performance of bio-inspired structure compared to its building blocks. Furthermore, an optimization method is implemented into the designing the bio-inspired structures for the first time, which enables us to perform the bio-inspired material design with the target of finding the most efficient geometries that can resist defects in their structure.

Financial incentives offered by payers to health care providers and patients have been identified as a key mechanism to lower costs while improving quality of care. How to effectively design incentive programs that can align the varying objectives of health care stakeholders, as well as predict programs' performance and stakeholders' decision response is an unresolved research challenge. The objective of this study is to establish a novel approach based on multiscale decision theory (MSDT) that can effectively model and efficiently analyze such incentive programs, and the complex health care system in general. The MSDT model captures the interdependencies of stakeholders, their decision processes, uncertainties, and how incentives impact decisions and outcomes at the payer, hospital, physician, and patient level. In the first part of this thesis, we study the decision processes of agents pertaining to the investment and utilization of imaging technologies. We analyze the payer-hospital-physician relationships and later extend the model to include radiologist and patient as major stakeholders in the second part of this thesis. We focus on a specific incentive program, the Medicare Shared Savings Program (MSSP) for Accountable Care Organizations (ACOs). The multi-level interactions between agents are mathematically formulated as a sequential non-cooperative game. We derive the equilibrium solutions using the subgame perfect Nash equilibrium (SPNE) concept and the backward induction principle, and determine the conditions under which the MSSP incentive leads to the desired outcomes of cost reduction and quality of care improvements. In the third part of this thesis, we study the multi-level decision making in chronic disease management. We model and analyze patients' and physicians' decision processes as a general-sum stochastic game with perfect information and switching control structure. We incorporate the Health Belief Model (HBM) as the theoretical foundation to capture the behavioral aspect of agents. We analyze how incentives and interdependencies affect patients' engagement in health-promoting activities and physicians' delivery of primary care services. We show that a re-alignment of incentives can improve the effectiveness of chronic disease management. / Ph. D.

In this research, a first-principles, non-equilibrium thermodynamic-ensemble approach applicable across all temporal and spatial scales is developed based on steepest-entropy-ascent quantum thermodynamics (SEAQT). The SEAQT framework provides an equation of motion consisting of both reversible mechanical dynamics and irreversible relaxation dynamics, which is able to describe the evolution of any state of any system, equilibrium or non-equilibrium. Its key feature is that the irreversible dynamics is based on a gradient dynamics in system state space instead of the microscopic mechanics of more traditional approaches. System energy eigenstructure and density operator (or ensemble probability distribution) describe the system and system thermodynamic state, respectively. Extensive properties (i.e., energy, entropy, and particle number) play a key role in formulating the equation of motion and in describing non-equilibrium state evolutions. All the concepts involved in this framework (i.e., eigentstructure, density operator, and extensive properties) are well defined at all temporal and spatial scales leading to the extremely broad applicability of SEAQT. The focus of the present research is that of developing non-equilibrium thermodynamic models based specifically on the irreversible part of the equation of motion of SEAQT and applying these to the study of pure relaxation processes of systems in non-equilibrium states undergoing chemical reactions and heat and mass diffusion. As part of the theoretical investigation, the new concept of hypo-equilibrium state is introduced and developed. It is able to describe any non-equilibrium state going through a pure relaxation process and is a generalization of the concept of stable equilibrium of equilibrium thermodynamics to the non-equilibrium realm. Using the concept of hypo-equilibrium state, it is shown that non-equilibrium intensive properties can be fundamentally defined throughout the relaxation process. The definition of non-equilibrium intensive properties also relies on various ensemble descriptions of system state. In this research, three SEAQT ensemble descriptions, i.e., the canonical, grand canonical, and isothermal-isobaric, are derived corresponding, respectively, to the definition of temperature, chemical potential, and pressure. To computationally and not just theoretically permit the application of the SEAQT framework across all scales, a density of states method is developed, which is applicable to solving the SEAQT equation of motion for all types of non-equilibrium relaxation processes. In addition, a heterogeneous multiscale method (HMM) algorithm is also applied to extend the application of the SEAQT framework to multiscale modeling. Applications of this framework are given for systems involving chemical kinetics, the heat and mass diffusion of indistinguishable particles, power cycles, and the complex, coupled reaction-diffusion pathways of a solid oxide fuel cell (SOFC) cathode. / Ph. D.

Most engineering systems have some degree of uncertainty in their input and operating parameters. The interaction of these parameters leads to the uncertain nature of the system performance and outputs. In order to quantify this uncertainty in a computational model, it is necessary to include the full range of uncertainty in the model. Currently, there are two major technical barriers to achieving this: (1) in many situations -particularly those involving multiscale phenomena-the stochastic nature of input parameters is not well defined, and is usually approximated by limited experimental data or heuristics; (2) incorporating the full range of uncertainty across all uncertain input and operating parameters via conventional techniques often results in an inordinate number of computational scenarios to be performed, thereby limiting uncertainty analysis to simple or approximate computational models. This first objective is addressed through combining molecular and macroscale modeling where the molecular modeling is used to quantify the stochastic distribution of parameters that are typically approximated. Specifically, an adsorption separation process is used to demonstrate this computational technique. In this demonstration, stochastic molecular modeling results are validated against a diverse range of experimental data sets. The stochastic molecular-level results are then shown to have a significant role on the macro-scale performance of adsorption systems. The second portion of this research is focused on reducing the computational burden of performing an uncertainty analysis on practical engineering systems. The state of the art for uncertainty analysis relies on the construction of a meta-model (also known as a surrogate model or reduced order model) which can then be sampled stochastically at a relatively minimal computational burden. Unfortunately these meta-models can be very computationally expensive to construct, and the complexity of construction can scale exponentially with the number of relevant uncertain input parameters. In an effort to dramatically reduce this effort, a novel methodology "QUICKER (Quantifying Uncertainty In Computational Knowledge Engineering Rapidly)" has been developed. Instead of building a meta-model, QUICKER focuses exclusively on the output distributions, which are always one-dimensional. By focusing on one-dimensional distributions instead of the multiple dimensions analyzed via meta-models, QUICKER is able to handle systems with far more uncertain inputs. / Ph. D.

This thesis presents a multiscale model of an industrial nylon-6 leacher. We develop several models at various spatial scales and implement them together in a simplistic, efficient way to develop an overall leacher model. We solve dynamic transport differential equations using the finite-volume method and method of lines in an in-house-developed FORTRAN program. We use the ODEPACK package of ordinary differential equation (ODE) solvers to solve our system of coupled ODEs. Our multiscale model performs transport, thermodynamic, physical property, and mass-transfer calculations at a finite-volume scale. We introduce two additional scales: a mesoscale, in which we perform computational fluid dynamic (CFD) simulations, and a molecular scale. Our CFD simulations solve for turbulent properties of fluid flowing over a packed bed. We incorporate the turbulent diffusivity of the fluid into our finite-volume leacher model. We perform molecular simulations and use the conductor-like screening model-segment activity coefficient (COSMO-SAC) model to generate solubility predictions of small, cyclic oligomers in water and ε-caprolactam. Additionally, we develop an extension of COSMO-SAC to model polymer species, which we refer to as Polymer-COSMO-SAC, and apply it to solve liquid-liquid equilibrium equations. We present a unique methodology to apply COSMO-based models to polymer species, which shows reasonable results for nylon-6. Because of the computational intensity of our Polymer-COSMO-SAC liquid-liquid equilibrium algorithm, we generate pre-computed tables of equilibrium predictions that we may import into our leacher model. Our integration of multiscale models maximizes efficiency and feasibility with accuracy. We are able to use our multiscale models to estimate necessary parameters, but we need to fit two mass-transfer related parameters to industrial data. We validate our model against the plant data and find average-absolute errors in the final mass percent of ε-caprolactam and cyclic dimer in polymer chips of 25.0% and 54.7%, respectively. Several plant data sets are suspected outliers and we believe an unforeseen equilibrium limitation may cause this discrepancy. If we remove these outlying data sets, we then find average-absolute errors of 7.5% and 19.3% for ε-caprolactam and cyclic dimer, respectively. We then use our validated model to perform application and sensitivity studies to gain critical insight into the leacher's operating conditions. / Master of Science

This dissertation focuses on two different types of responses of Earth; that is, seismic and electromagnetic, and aims to better understand Earth processes at a wider range of scales than those conventional approaches offer. Electromagnetic responses resulting from the subsurface diffusion of applied electromagnetic fields through heterogeneous geoelectrical structures are utilized to characterize the underlying geology. Geology exhibits multiscale hierarchical structure which brought about by almost all geological processes operating across multiple length scales and the relationship between multiscale electrical properties of underlying geology and the observed electromagnetic response has not yet been fully understood. To quantify this relationship, the electromagnetic responses of textured and spatially correlated, stochastic geologic media are herein presented. The modelling results demonstrate that the resulting electromagnetic responses present a power law distribution, rather than a smooth response polluted with random, incoherent noise as commonly assumed; moreover, they are examples of fractional Brownian motion. Furthermore, the results indicate that the fractal behavior of electromagnetic responses is correlated with the degree of the spatial correlation, the contrasts in ground electrical conductivity, and the preferred orientation of small-scale heterogeneity. In addition, these inferences are also supported by the observed electromagnetic responses from a fault zone comprising different lithological units and varying wavelengths of geologic heterogeneity. Seismic signals generated by aftershocks are generally recorded by local aftershock networks consisted of insufficient number of stations which result in strongly spatially-aliased aftershock data. This limits aftershock detections and locations at smaller magnitudes. Following the 23 August 2011 Mineral, Virginia earthquake, to drastically reduce spatial aliasing, a temporary dense array (AIDA) consisting of ~200 stations at 200-400 m spacing was deployed near the epicenter to record the 12 days of the aftershocks. The backprojection imaging method is applied to the entire AIDA dataset to detect and locate aftershocks. The method takes advantage of staking of many seismograms and improves the signal-to-noise ratio for detection. The catalog obtained from the co-deployed, unusually large temporal traditional network of 36 stations enabled a quantitative comparison. The aftershock catalog derived from the dense AIDA array and the backprojection indicates event detection an order of magnitude smaller including events as small as M–1.8. The catalog is complete to magnitude –1.0 while the traditional network catalog was complete to M–0.27 for the same time period. The AIDA backprojection catalog indicate the same major patterns of seismicity in the epicentral region, but additional details are revealed indicating a more complex fault zone and a new shallow cluster. The b-value or the temporal decay constant were not changed by inclusion of the small events; however, they are different for two completeness periods and are different at shallow depth than greater depth. / Ph. D. / This dissertation revolves around two ends of geophysics: seismology and electromagnetics. The electromagnetic method of exploration geophysics aims to characterize the underlying geology by evaluating the electromagnetic responses resulting from the interaction between the applied electromagnetic fields and the subsurface electrical properties. In case of rough geology comprising heterogeneity at every scale, the electromagnetic responses are more complicated than the response of a piecewise smooth Earth structure. Most analyses treat the responses of small-scale heterogeneities as random, uncorrelated noise. Here, more realistic geologic models comprising spatially-correlated, fine-scale heterogeneities are incorporated into electromagnetic modeling to better understand the relationship between the causative multiscale geoelectrical heterogeneities and the electromagnetic responses. The numerical results indicate that these electromagnetic responses are not random as commonly assumed, in contrast, they are repeatable and fractally distributed presenting spatial fluctuations that appear on all length scales. Moreover, the numerical results indicate that the fractal behavior of electromagnetic responses is correlated with the degree of the spatial correlation, the contrasts in ground electrical conductivity, and the preferred orientation of small-scale heterogeneity. In addition, the analysis of the observed electromagnetic responses from a fault zone comprising multiscale heterogeneity also support these inferences. Seismic signals generated by aftershocks are generally recorded by local aftershock networks consisted of insufficient number of stations which result in strongly spatially-aliased aftershock data. This limits aftershock detections and locations at smaller magnitudes. Following the 23 August 2011 Mineral, Virginia earthquake, to drastically reduce spatial aliasing, a temporary dense array (AIDA) consisting of ∼200 stations at 200-400 m spacing was deployed near the epicenter to record the 12 days of the aftershocks. The backprojection imaging method is applied to the entire AIDA dataset to detect and locate aftershocks. The method takes advantage of summing of many seismograms and improves the signal-to-noise ratio for detection. The catalog obtained from the co-deployed, unusually large temporal traditional network of 36 stations enabled a quantitative comparison. The aftershock catalog derived from the dense AIDA array and the backprojection indicates event detection an order of magnitude smaller. The AIDA backprojection catalog indicate the same major patterns of seismicity in the epicentral region, but additional details are revealed indicating a more complex fault zone and a new shallow cluster. The decay of aftershock rate and the distribution of earthquakes with respect to the magnitude do not show a significant change by inclusion of the small events; however, they differ at shallow and greater depth, and for different completeness periods.

Ce travail de recherche propose un modèle numérique pour simuler les procédés par infusion de résine en utilisant la méthode des éléments finis. Ce modèle permet de représenter l'écoulement d'une résine liquide dans des préformes poreuses subissant de grandes déformations. Dans cette étude, une modélisation macroscopique est utilisée. Au niveau du procédé, une zone de résine liquide est déposée sur les préformes. Ces dernières étant considérées comme un milieu poreux. Les équations de Stokes et de Darcy sont utilisées pour modéliser l'écoulement de la résine respectivement dans le drainant et dans les préformes. L'originalité du modèle réside dans le fait qu'un seul maillage est utilisé pour les deux milieux. La discrétisation est réalisée avec des éléments mixtes : dans Stokes, des éléments P1+/P1 sont utilisés et dans Darcy, des éléments P1/P1 stabilisés avec une formulation multi-échelle sont employés. Des fonctions distances signées sont utilisées pour représenter l'interface entre Stokes-Darcy et pour représenter le front de résine. Concernant la déformation des préformes, une formulation Lagrangienne réactualisée est utilisée. Dans cette formulation Lagrangienne, le comportement des préformes humides est représenté à l'aide du modèle de Terzaghi dans lequel les préformes sèches ont un comportement élastique non-linéaire. La perméabilité est reliée à la porosité via la relation de Carman-Kozeny. Celle-ci est déterminée à partir de l'équation de conservation de la masse. Ce modèle a été implémenté dans ZéBuLoN. Plusieurs simulations numériques d'infusion de résine sont présentées à la fin de ce manuscrit. / This work proposes a numerical model to simulate the manufacturing processes by resin infusion using the finite element method. This model allows to represent the resin flow into porous preforms, which are themselves subject to large deformations. In this study, a macroscopic description is used. The preforms are considered as a porous medium. The Stokes and the Darcy equations are used respectively to describe the resin flow into the liquid zone and into the preforms.The originality of the model consists in using one single unstructured mesh. The discretization is ensured by using a mixed velocity-pressure formulation. Indeed, a P1/P1 formulation is employed throughout the entire discretized domain, stabilized in the Darcy region with a multi-scale formulation and in the Stokes subdomain with a hierarchical-based bubble, i.e. a P1+/P1 finite element. Signed distance functions are used both to represent the Stokes-Darcy interface and to capture the moving flow front. Concerning the deformations of the preforms, an updated Lagrangian scheme is used. In the Lagrangian formulation, the behavior of the wet preforms is represented by using the Terzaghi model in which the dry preforms have a non-linear elastic behavior. The permeability depends on the porosity through the Carman-Kozeny relationship. This model has been implemented in Zset. Several numerical simulations of manufacturing processes by resin infusion are presented at the end of this manuscript.

Este trabalho propõe uma técnica de modelagem multiescala concorrente do concreto considerando duas escalas distintas: a mesoescala, onde o concreto é modelado como um material heterogêneo, e a macroescala, na qual o concreto é tratado como um material homogêneo. A heterogeneidade da estrutura mesoscópica do concreto é idealizada considerando três fases distintas, compostas pelos agregados graúdos e argamassa (matriz), estes considerados materiais homogêneos, e zona de transição interfacial (ZTI), tratada como a parte mais fraca entre as três fases. O agregado graúdo é gerado a partir de uma curva granulométrica e posicionado na matriz de forma aleatória. Seu comportamento mecânico é descrito por um modelo constitutivo elástico-linear, devido a sua maior resistência quando comparado com as outras duas fases do concreto. Elementos finitos contínuos com alta relação de aspecto em conjunto com um modelo constitutivo de dano são usados para representar o comportamento não linear do concreto, decorrente da iniciação de fissuras na ZTI e posterior propagação para a matriz, dando lugar à formação de macrofissuras. Os elementos finitos de interface com alta relação de aspecto são inseridos entre todos os elementos regulares da matriz e entre os da matriz e agregados, representando a ZTI, tornando-se potenciais caminhos de propagação de fissuras. No estado limite, quando a espessura do elemento de interface tende a zero (h ?0) e, consequentemente, a relação de aspecto tende a infinito, estes elementos apresentam a mesma cinemática da aproximação contínua de descontinuidades fortes (ACDF), sendo apropriados para representar a formação de descontinuidades associados a fissuras, similar aos modelos coesivos. Um modelo de dano à tração é proposto para representar o comportamento mecânico não linear das interfaces, associado à formação de fissuras, ou até mesmo ao eventual fechamento destas. A fim de contornar os problemas causados pela malha de elementos finitos de transição entre as malhas da macro e da mesoescala, que, em geral, apresentam diferenças expressivas 5 de refinamento, utiliza-se uma técnica recente de acoplamento de malhas não conformes. Esta técnica é baseada na definição de elementos finitos de acoplamento (EFAs), os quais são capazes de estabelecer a continuidade de deslocamento entre malhas geradas de forma completamente independentes, sem aumentar a quantidade total de graus de liberdade do problema, podendo ser utilizados tanto para acoplar malhas não sobrepostas quanto sobrepostas. Para tornar possível a análise em multiescala em casos nos quais a região de localização de deformações não pode ser definida a priori, propõe-se uma técnica multiescala adaptativa. Nesta abordagem, usa-se a distribuição de tensões da escala macroscópica como um indicador para alterar a modelagem das regiões críticas, substituindo-se a macroescala pela mesoescala durante a análise. Consequentemente, a malha macroscópica é automaticamente substituída por uma malha mesoscópica, onde o comportamento não linear está na iminência de ocorrer. Testes numéricos são desenvolvidos para mostrar a capacidade do modelo proposto de representar o processo de iniciação e propagação de fissuras na região tracionada do concreto. Os resultados numéricos são comparados com os resultados experimentais ou com aqueles obtidos através da simulação direta em mesoescala (SDM). / A concurrent multiscale analysis of concrete is presented, in which two distinct scales are considered: the mesoscale, where the concrete is modeled as a heterogeneous material and the macroscale that treats the concrete as a homogeneous material. The mesostructure heterogeneities are idealized as three phase materials composed of the coarse aggregates, mortar matrix and the interfacial transition zone (ITZ). The coarse aggregates are generated from a grading curve and placed into the mortar matrix randomly. Their behavior is described using an elastic-linear constitutive model due to their significant higher strength when compared with the other two phases of the concrete. Special continuum finite elements with a high aspect ratio and a damage constitutive model are used to describe the nonlinear behavior associated to the propagation of cracks, which initiates in the ITZ and then propagates to the mortar matrix given place to a macro-crack formation. These interface elements with a high aspect ratio are inserted in between all regular finite elements of the mortar matrix and in between the mortar matrix and aggregate elements, representing the ITZ. In the limit case, when the thickness of interface elements tends to zero (h ?0) and consequently the aspect ratio tends to infinite, these elements present the same kinematics as the continuous strong discontinuity approach (CSDA), so that they are suitable to represent the formation of discontinuities associated to cracks, similar to cohesive models. A tensile damage model is proposed to model the nonlinear mechanical behavior of the interfaces, associated to the crack formation and also to the possible crack closure. To avoid transition meshes between the macro and the mesoscale meshes, a new technique for coupling non-matching meshes is used. This technique is based on the definition of coupling finite elements (CFEs), which can ensure the continuity of displacement between independent meshes, without increasing the total number of degrees of freedom of the problem. This technique can be used to couple non-overlapping and overlapping meshes.To make possible the concurrent multiscale analysis, where the strain localization region cannot be defined a priori, an adaptive multiscale model is proposed. In this approach the macroscale stress distribution is used as an indicator to properly change from the macroscale to the mesoscale modeling in the critical regions during the analysis. Consequently, the macroscopic mesh is automatically replaced by a mesoscopic mesh where the nonlinear behavior is imminent. A variety of tests are performed to show the ability of the proposed methodology in predicting the behavior of initiation and propagation of cracks in the tensile region of the concrete. The numerical results are compared with the experimental ones or with those obtained by the direct simulation in mesoscale (DSM).