Este documento apresenta o desenvolvimento de um modelo computacional para prover informações aos tomadores de decisão frente as alternativas de investimentos na expansão das redes de distribuição. Estas redes, por estarem localizadas em áreas urbanas e rurais, apresentam variadas topologias e número elevado de equipamentos na sua composição, dificultando a análise técnica e econômica diante do número elevado de variáveis envolvidas. O método proposto inclui a avaliação dos principais fatores que tem influência na qualidade da energia, incluindo os índices de confiabilidade. Adicionalmente, aspectos relacionados com os custos financeiros e eficiência são considerados para diferentes condições de operação. As decisões são feitas considerando o desempenho dos transformadores e requisitos exigidos em normas de confiabilidade da agência reguladora do sistema elétrico. Perda de vida em transformadores, perdas técnicas, níveis de tensão e fator de utilização são avaliados em conjunto com a confiabilidade no fornecimento de energia, a fim de combinar os resultados destes critérios com indicadores financeiros. O resultado é a obtenção de indicadores individuais e coletivos que facilitam a tomada de decisão no direcionamento de investimentos a curto e médio prazo, proporcionando de maneira prática e fácil a comparação entre alimentadores, subestações ou mesmo regiões que compõem um sistema de distribuição. Por fim, para verificação do método, realiza-se um estudo de caso num sistema de distribuição real onde são analisados e discutidos os resultados. / This document aims to present the development of a computational model to provide information for decision-making toward the alternatives of investments in distribution system planning. These networks, as they are located in urban and rural areas, can present varied topologies and large numbers of equipment in its arrangement, making it difficult in performing technical and economic analysis due to the large number of variables involved. The proposed method includes the evaluation of main factors that have influence on power quality, including the reliability indexes. Additionally, aspects related to economic costs and energy efficiency are considered for different operating conditions. The decision is made taking into account the performance of medium voltage transformers and distribution system reliability standards of Regulatory Agency requirements. Loss of life in transformers, technical losses, voltage levels, and utilization factor are assessment together with continuity and availability of power supply in order to combine the results of these criteria with financial indicators. The result is the achievement of individual and collective indicators that facilitate decision making in directing investments in the short and medium term, providing a practical and easy way to compare feeders, substations or even regions that comprise a distribution system. The methodology is applied in a real distribution system, and the result are analyzed and discussed.

No presente trabalho, a conformabilidade do aço inoxinoxidável AISI 201 foi comparada com o aço inoxidável AISI 304. O aço inoxidável AISI 201 é uma liga baixo níquel ligado com manganês e nitrogênio. Nesse estudo a conformabilidade dos dois materiais foi examinada por meio de ensaios tecnológicos como ensaio de tração, determinação da curva de escoamento, determinação do índice de anisotropia, ensaio Erichsen, dureza, composição química, simulação computacional do processo de estampagem e estampagem dos blanks. O aço inoxidável AISI 201 possui propriedades como limite de escoamento e tensão de ruptura superior as do inoxidável AISI 304. Os dois aços possuem uma similaridade na anisotropia. A máxima altura alcançada no momento da fratura pelo ensaio Erichsen também é semelhante para os dois materiais. O aço inox AISI 201 apresenta uma dureza maior que o inoxidável AISI 304. Na composição química os dois aços apresentam elementos fora do especificado, caracterizando um problema de qualidade na fabricação desses aços. A simulação computacional do processo de estampagem apresentou uma redução de espessura na região mais critica, sem comprometer a estampagem do componente. Uma vez estampadas, obteve-se peças sem indícios de trincas, conforme previsto pela simulação computacional. Os resultados desse trabalho sugerem a possibilidade de utilização do aço inoxidável AISI 201 como opção para substituição ao inoxidável AISI 304, tendo uma observação a ser feita quanto aos cuidados na qualidade durante a fabricação do aço para que atenda os padrões exigidos. / In this study, the formability of stainless steel AISI 201 was compared to stainless steel AISI 304 stainless steel AISI 201 is a low alloy nickel alloyed with manganese and nitrogen. In this study, the formability of the two materials was examined by means of technological tests such as tensile test, determination of the flow curve, determining the anisotropy index, Erichsen test, hardness, chemical composition, computer simulation of the stamping process and stamping of the blanks. Stainless steel AISI 201 has properties such as yield strength and higher breakdown voltage of the stainless steel AISI 304. The two steels have a similarity in anisotropy. The maximum height reached at the time of fracture by Erichsen test is also similar for the two materials. Stainless steel AISI 201 has a hardness greater than the stainless steel AISI 304. In chemistry the two steels have elements outside the specified, featuring a quality problem in manufacturing these steels. A computer simulation of the printing process showed a reduction in thickness in the most critical region, without compromising the component stamping. Once stamped, gave no broken pieces of evidence as provided by the computer simulation. The findings suggest the possibility of use of stainless steel AISI 201 as an option to replace the stainless steel AISI 304, with a point to be made about the care as during the manufacture of steel that meets the required standards.

Mapas simpléticos são bem conhecidos por preservarem o volume do espaço de fase em dinâmica Hamiltoniana e são particularmente apropriados para problemas que requerem longos tempos de integração. Nesta tese nós desenvolvemos abordagens baseadas em mapas simpléticos para o acoplamento de multi sub-sistemas/domínios astrofísicos/códigos de simulação, para integração eficiente de sistemas de N-corpos auto-gravitantes com grandes variações nas escalas de tempo características. Nós estabelecemos uma família de 48 novos mapas simpléticos baseados numa separação Hamiltoniana recursiva, que permite que o acoplamento ocorra de uma maneira hierárquica, contemplando assim todas as escalas de tempo das interações envolvidas. Nossa formulação é geral o suficiente para permitir que tal método seja utilizado como receita para combinar diferentes fenômenos físicos, que podem ser modelados independentemente por códigos especializados. Nós introduzimos também uma separação Hamiltoniana baseada em Hamiltonianos de Kepler, para resolver o problema gravitacional geral de N-corpos como uma composição de N2 problemas de 2-corpos. O método resultante é exato para cada problema de 2-corpos individual e produz resultados rápidos e precisos para sistemas de N-corpos quase- Keplerianos, como sistemas planetários ou um aglomerado de estrelas que orbita um buraco-negro supermassivo. O método é também apropriado para integração de sistemas de N-corpos com hierarquias intrínsecas, como um aglomerados de estrelas com binárias compactas. Nós apresentamos a implementação dos algoritmos mencionados e descrevemos o nosso código tupan, que está publicamente disponível na seguinte url: https://github.com/ggf84/tupan. / Symplectic maps are well know for preserving the phase space volume in Hamiltonian dynamics and are particularly suited for problems that require long integration times. In this thesis we develop approaches based on symplectic maps for the coupling of multi sub-systems/astrophysics domains/simulation codes for efficient integration of self-gravitating N-body systems with large variation in characteristic time-scales. We establish a family of 48 new symplectic maps based on a recursive Hamiltonian splitting, which allow the coupling to occur in a hierarchical manner, thus contemplating all time-scales of the involved interactions. Our formulation is general enough to allow that such method be used as a recipe to combine different physical phenomena which can be modeled independently by specialized simulation codes. We also introduce a Keplerian-based Hamiltonian splitting for solving the general gravitational Nbody problem as a composition of N2 2-body problems. The resulting method is precise for each individual 2-body solution and produces quick and accurate results for near-Keplerian N-body systems, like planetary systems or a cluster of stars that orbit a supermassive black-hole. The method is also suitable for integration of N-body systems with intrinsic hierarchies, like a star cluster with compact binaries. We present the implementation of the mentioned algorithms and describe our code tupan, which is publicly available on the following url: https://github.com/ggf84/tupan.

Esta pesquisa foi realizada em um Instituto Federal de Educação, Ciência e Tecnologia e dedicou-se ao estudo de alguns conceitos de Física Quântica na formação de professores de Física, sob a perspectiva do referencial sociocultural, tendo o foco dirigido à complementaridade de segundo tipo, uma das teses da interpretação de Copenhagen, e ao emaranhamento quântico. O marco teóricometodológico adotado foi a teoria da mediação de Vygotsky e a filosofia translinguística de Bakhtin. Foram realizados três estudos qualitativos com o objetivo de investigar quais são os significados socialmente construídos sobre conceitos de Física Quântica por professores de Física em formação, sendo um estudo piloto e duas intervenções didáticas junto a uma disciplina de Física Moderna e Contemporânea em nível de graduação. As atividades de ensino foram conduzidas a partir de exposições teóricas e simulação computacional com o Interferômetro Virtual de Mach-Zehnder, orientada por roteiros exploratórios. A partir do estudo dos fenômenos da interferência e emaranhamento quânticos nesse dispositivo semiótico, buscou-se compreender em que medida as ações mediadas pelo uso do software como recurso didático atuaram no sentido de promover entre os professores uma construção conceitual colaborativa sobre esses conceitos da Física Quântica. A análise do discurso, fundamentada em conceitos bakhtinianos, possibilitou uma melhor compreensão sobre o processo de construção de conceitos culturalmente válidos acerca da interferência quântica, da complementaridade e do emaranhamento. No contexto explorado, verificamos que se estabeleceu o dilema da manutenção de visões antagônicas, o que pode ser evidenciado a partir dos modelos explicativos elaborados pelos estudantes. Para aqueles que atribuíram um caráter corpuscular aos fótons, os fenômenos intermediários permanecem como um conflito ainda não resolvido, enquanto para os demais a complementaridade onda-partícula passou a compor um nova classe de conceitos. A noção de não-localidade, necessária para compreender o emaranhamento quântico, não foi reconhecida pelos estudantes, evidenciando a forte presença de noções clássicas nos modelos explicativos utilizados na tentativa de obter uma explicação plausível para o fenômeno. / The research presented in this thesis was developed in a Federal Institute of Education, Science and Technology is directed to research about an approach of the some foundations of Quantum Physics in pre-service teachers’ training, based on the sociocultural perspective and this work focuses on the second kind of complementarity, one of the thesis of the Copenhagen interpretation, and the quantum entanglement. The theoretical and methodological framework adopted in this research is the Vygotsky’s theory of mediation and Bakhtin’s translinguistics philosophy. Three qualitative researches were performed in order to investigate which are the socially constructed meanings on QP foundations in the context of pre-service teachers’ training: one in an introductory course with physics teachers and two along a didactic intervention during a pre-service teachers’ training course on Modern and Contemporary Physics. The teaching activities were conducted mediated by theoretical lectures and computer simulation using the Virtual Mach- Zehnder Interferometer, conducted by the exploratory guide. From the study of interference phenomena and quantum entanglement in this semiotic tool, one tries to understand how the pupils actions are mediated by the software in the sense of promoting a collaborative conceptual construction among teachers about these foundations of QP. Discourse analysis grounded on Bakhtin’s concepts allowed us to evaluate the potentialities of the software and a better understanding of the process of conceptual construction that are culturally valid about quantum interference, complementarity and entanglement. In this context we perceive that was created the dilemma of maintaining competing visions, which can be evidenced from the explanatory models designed by the students. For the students that admitted a corpuscular character to photons so that intermediate phenomena remain as a unresolved cognitive conflict, while for others the wave - particle complementarity has became a new class of concepts. The non-locality notion, pre-requisite to understand the quantum entanglement, wasn’t learned by the students, demonstrating the strong presence of classical notions in explanatory models used to obtain a plausible explanation for the phenomenon.

Trocadores de Calor Solo-Ar (TCSA) são dispositivos usados para melhorar o condicionamento térmico de ambientes construídos. Eles funcionam através da ventilação do ar por um ou mais dutos enterrados, utilizando o solo como fonte ou sumidouro de calor. Em virtude da defasagem entre as temperaturas do ar e das camadas superficiais da Terra, é possível resfriar o ar no verão e aquecê-lo no inverno. Seus princípios de operação baseiam-se na mecânica dos fluidos e transferência de calor, áreas onde a teoria Construtal tem sido usada para melhorar o desempenho, ou reduzir imperfeições, de vários sistemas térmicos, após uma estratégica análise de suas geometrias (ou de seu design), também chamada de método Design Construtal que se fundamenta na lei construtal. Dessa lei, projetos em engenharia devem começar, por exemplo, descobrindo as arquiteturas que facilitam o escoamento entre uma fonte pontual e um volume, ou vice-versa. Explorando esses conceitos para TCSA, onde a corrente é o calor que flui entre os dutos e o solo, esta tese objetiva centralmente avaliar o desempenho térmico desses dispositivos em função de possíveis desenhos de suas estruturas. Particularmente, isso foi feito: (1) usando um duto e variando seu diâmetro e vazão de ar; (2) inserindo novos dutos, até um total de cinco, mantendo a vazão de ar constante, e estudando diversas disposições geométricas. Numericamente, foram empregados dois modelos testados e validados a partir de dados experimentais, que foram simulados no código comercial de dinâmica dos fluidos computacional, FLUENT, o qual baseia-se em métodos de volumes finitos. As geometrias e malhas foram geradas no software GAMBIT. Dos resultados numéricos, esta tese também avançou desenvolvendo novos modelos para o conceito de potencial térmico instantâneo, que é uma medida das diferenças de temperatura entre a saída e a entrada dos dutos de TCSA. Descobriu-se que estes modelos podiam ser simplificados em termos de funções senoidais, facilitando análises e definições analíticas subsequentes, criando um quadro resumo para grandes volumes de dados simulados. Das avaliações dos resultados nos estudos com um duto, foram encontradas faixas para a obtenção de rendimentos térmicos anuais estimados em pelo menos 80%, com mudanças no diâmetro e/ou vazão do ar. Analisando múltiplos dutos, foram reveladas melhorias significativas de rendimento com: (1) o aumento no número de dutos; (2) a redução, até certos limites, da razão entre os espaçamentos verticais e horizontais entre eles; (3) a ampliação da razão entre o volume da instalação e o volume do domínio computacional. Destaca-se, ainda, que em todos os casos foram calculados balanços energéticos, com a estimativa das taxas (e quantidade) de calor trocado, bem como da energia elétrica consumida pelos ventiladores, mostrando sua viabilidade no condicionamento de ar com redução do consumo de energia elétrica. / The Earth-Air Heat Exchangers (EAHE) are devices used to improve the thermal conditions of built environments. They work by blowing the air inside buried ducts to use the soil as heat source or sink. Due to a phase difference between the air and ground temperatures, it is possible to cool the air in the summer and heat it in the winter. These operating principles are based on the areas of fluid mechanics and heat transfer, where the Constructal theory has been used to improve the performance, or reduce imperfections, of several thermal systems, after a strategic evaluation of their geometries (or their design), also called Constructal Design method, which is anchored by the constructal law. From this law, engineering projects should begin, for example, finding the architectures which ease the flow between a source point and a volume, or vice versa. Bringing these concepts to EAHE, where the current is the heat which flows between the ducts and the soil, this thesis mainly aims to assess the thermal performance of EAHE due to possible designs of their layouts. Particularly, this was done by: (1) adopting one duct and varying its diameter and/or the air flow; (2) using up to five ducts, keeping the air flow constant, and studying various geometric shapes. Numerically, two models were used, which have been tested and validated from experimental data, and the simulations were done in the CFD software FLUENT, which is based in the finite volume methods. The geometries and meshes were generated by the software GAMBIT. From the numerical results, this thesis also pursued developing new models for the concept of the instantaneous thermal potential, which is a measure of the temperatures differences between the ducts outlet and inlet. It was found that these models can be simplified in terms of sinusoidal functions, helping analysis, definitions, and creating a framework to summarize large volumes of simulated data. From the evaluation of the data in the studies with one duct, ranges for changes in their diameter and/or air flow were found to obtain at least 80% of estimated annual efficiency. Analyzing multiple ducts, significant improvements in efficiency were also obtained by: (1) increasing the number of ducts; (2) reducing, to some extent, the ratio between their vertical and horizontal spacings; (3) increasing the ratio between the installation volume and the computational domain. It should also be highlighted that in all cases the energy balances were computed, to find the rates (and amounts) of exchanged heat, as well as the electric energy consumed by the fans, showing the economic viability of using EAHE as devices for air conditioning.

A análise e projeto de grandes estruturas de engenharia requerem, geralmente, a utilização de fundações profundas baseadas em estacas de aço ou concreto armado. Sendo o problema de natureza tridimensional pelas condições de contorno e a variação espacial das propriedades do solo, necessita-se de uma ferramenta computacional capaz de simular esse tipo de problema em situações de interesse, tais como carregamento sísmico ou quase-estático por consolidação. Neste trabalho, propõe-se um modelo numérico tridimensional baseado no método dos elementos finitos sob pequenas deformações para a modelagem do problema de interação solo-estaca sob carregamento estático, quase-estático e dinâmico. Elementos finitos hexaédricos são utilizados para representar o meio poroso saturado e as estacas de concreto armado. Considera-se a interação parcial entre ambos meios, mediante a inclusão de elementos de interface capazes de simular separação e escorregamento. A não-linearidade física de todos os materiais envolvidos é considerada mediante a utilização da teoria de plasticidade, onde esquemas de integração explícita são utilizados. Um modelo constitutivo baseado na teoria de plasticidade generalizada é utilizado para simular o fenômeno de liquefação em areias. No caso da estaca de concreto armado, utilizou-se um modelo de fissuras distribuídas para representar o processo de fissuração, enquanto a armadura de reforço pode ser representada por uma camada distribuída equivalente ou por um modelo discreto incorporado. Exemplos numéricos são apresentados para validar a correta implementação do modelo numérico. / Deep foundations using steel or reinforced concrete piles are conmmonly used in large civil engineering structures. Due to the three-dimensional nature of the problem, its boundary conditions and spatial variability of soil properties, it is necessary to employ numerical models that must be able to simulate this problem under seismic loading as well as under quasi static consolidation processes. A three-dimensional numerical model based on the finite element method, for the static, quase-static and dynamic analysis of the soil-pile interaction problem under the small strain assumption is proposed in this work. For this purpose, coupled hexahedral finite elements are used to represent soil and concrete piles. Interface elements are used to simulate slip, bonding and opening processes at the soil-pile interface. Material nonlinear behaviour of the soil and concrete are considered by using the theory of plasticity, where explicit integration schemes are used. A suitable constitutive model based on the generalized theory of plasticity is employed to represent the liquefaction phenomena in sands under cyclic loading. For the concrete pile, a smeared approach is used to represent concrete cracking due to traction stresses, while both an equivalent smeared layer representation or a discrete unidimensional representation can be used for modeling the reinforcing steel. Numerical examples are presented in order to validate the implementation of the numerical model.

The opening of energy markets for the production of energy from renewable sources has been positively influencing the economic investments in the sector. The growing investment from the financial market is based on the reduction of market risks inherent in projects related to the generation of clean energy. The emergence of new methodologies for the evaluation of projects concomitant with the improvement of the traditional techniques applied to the renewable energy sector are fundamental for the reduction of market uncertainties and the liquidity of investment projects in the energy sector. The valuation of investment options by the financial market occurs through the use of traditional techniques such as Discounted Cash Flow (CDF) or real options. The real option approach is considered by many investors to be complex, making discounted cash flow the most widely used approach. This thesis presents a contribution to the improvement of the use of the discounted cash flow technique in the evaluation of investment projects applied to the renewable energy sector considering the use of perpetuity or terminal value. This study is based on the evaluation of two hybrid water systems already in the literature with methods used in the evaluation of financial investments and discusses the results seeking to lay the foundations for a new way to evaluate hybrid systems based on renewable energy resources. The two systems under study were simulated with the Homer software, and in one, the choice for the final configuration did not result in the ideal configuration suggested by the software. The applied method allowed to justify this originally divergent choice of the result presented by the software Homer. The final results presented are the optimization spaces (one of the traditional ways in which the software Homer presents its results) constructed according to the presented evaluations. / A abertura de mercados energéticos para a produção de energia proveniente de fontes renováveis vem influenciando de forma positiva os investimentos econômicos no setor. O crescente investimento oriundo do mercado financeiro se baseia na redução de riscos de mercado inerentes aos projetos ligados à geração de energia limpa. O surgimento de novas metodologias de avaliação de projetos concomitante ao aprimoramento de técnicas consideradas tradicionais aplicadas ao setor de energias renováveis são fundamentais para a diminuição das incertezas de mercado e a liquidez dos projetos de investimento no setor energético. A avaliação de opções de investimento pelo mercado financeiro ocorre através do uso de técnicas tradicionais como o Fluxo de Caixa Descontado (FCD) ou as opções reais. A abordagem por opções reais é considerada por muitos investidores como complexa, tornando o fluxo de caixa descontado como a abordagem utilizada de forma mais ampla. Esta tese apresenta uma contribuição para o aprimoramento do uso da técnica do fluxo de caixa descontado na avaliação de projetos de investimentos aplicada ao setor de energias renováveis considerando a utilização da perpetuidade ou valor terminal. Este estudo se baseia na avaliação de dois sistemas híbridos hídricos já na literatura com métodos utilizados na avaliação de investimentos financeiros e discute os resultados buscando construir as bases para uma nova maneira de avaliar sistemas híbridos com base em recursos energéticos renováveis. Os dois sistemas em estudo foram simulados com o software Homer e em um deles a escolha para a configuração final não resultou na configuração ótima sugerida pelo software. O método aplicado permitiu justificar essa escolha originalmente divergente do resultado apresentado pelo software Homer. Os resultados finais apresentados são os espaços de otimização (uma das formas tradicionais na qual o software Homer apresenta seus resultados) construídos de acordo com as avaliações apresentadas.

Conselho Nacional de Desenvolvimento Científico e Tecnológico / Nesta dissertação é modelado computacionalmente o processo de corrosão por pites em aço inoxidável. Foi implementado um programa computacional que modela a distribuição da concentração local de espécies agressivas, o grau de dano do filme passivo, a suscetibilidade a pites e a distribuição dos pites, com base no modelo de PUNCKT et al. [2004]. Como uma contribuição desta dissertação foi considerada uma dependência entre a concentração local de espécies agressivas e a corrente anódica simulada. Também foi considerada, como outra contribuição do presente trabalho, a existência de pites estáveis e o efeito determinístico decorrente das regiões de inclusões de MnS em aços inoxidáveis. A geração de pites sobre uma superfície metálica é modelada usando um processo randômico de nucleação de pites, e a solução do modelo matemático permite a visualização da superfície corroída. / In the present work a computational model of the pitting corrosion process in stainless steal is developed. A computational program was implemented for the determination of the distribution of the local concentration of aggressive species, the degree of the passive film damage, the susceptibility of pitting and the pitting distribution, based on the PUNCKT et al. [2004] model. As a contribution of the present work dependence between the local concentration of aggressive species and the simulated anodic current is considered. It was also considered, as another contribution, the existence of stable pits and the deterministic effect of regions of MnS inclusions in stainless steel. The generation of pits over a metallic surface is modeled using a random process of pits nucleation, and the mathematical model solution gives the corroded surface visualization.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Apresentam-se aqui os resultados da simulação do Espectrômetro de Dipolo do Forward Próton Detector (FPD) do experimento DØ usando o Geant4 O FPD consiste de um conjunto de espectrômetros de momentum localizados no tubo do feixe em ambos os lados do detector DØ com o objetivo de detectar prótons e/ou antiprótons produzidos em eventos difrativos resultantes das colisões próton-antipróton no centro do DØ Para determinar quais os fatores influenciam a resposta do FPD, utilizamos eventos com = 2 TeV gerados pelo programa PHOJET Nesta dissertação realizou-se a simulação de um detector bem conhecido, neste caso os espectrômetros de dipolo do FPD a fim de entender como desenvolver simulações utilizando o Geant4 para futuras aplicações no LHC / We present here the results from the simulation of the Dipole Spectrometer of the Forward Proton Detector (FPD), using the Geant4 simulation package. The FPD consists of momentum spectrometers placed in the beam pipe in both sides of DØ experiment. It was built to tag protons and/or antiprotons produced in diffractive events as a result of the collisions at the center of DØ. In order to estimate the response of our simulation we used events at = 2 TeV generated with PHOJET. Our motivation to simulate the dipole spectrometer of FPD is to learn and understand how Geant4 operates in order to use it in future applications at LHC

O presente estudo busca desenvolver um modelo baseado em agentes para estudar em ambiente simulado baseado em sistemas computacionais as propriedades emergentes decorrentes da aplicação e execução da lei penal. Os modelos baseados em agentes têm-se revelado uma alternativa bastante promissora no estudo de sistemas complexos onde a grande quantidade de variáveis dinâmicas e interagentes tornam impossível o uso de técnicas puramente analíticas para a sua abordagem. Especialmente considerando que normas legais penais feitas para regular a conduta das pessoas em sociedade devem durar, a dimensão tempo acaba recebendo uma considerável relevância, o que incentiva ainda mais o uso de sistemas multiagentes, particularmente propícios nestes casos. Buscou-se construir um modelo no qual fosse possível analisar qualquer lei de natureza penal, quaisquer que fossem suas características e a quem estivesse destinada. Chegou-se, para tanto, no construto Cidadão- Estado-Oportunidade-Ambiente, que definem precisamente os contornos de qualquer conduta criminosa e que foram modelados em ambiente mutiagente. Para construir os fundamentos do comportamento dos agentes foram utilizados especialmente os referenciais teóricos de Shavell e Polinsky (2000, 2004), Becker (1968), Cohen e Felson (1979) e Clarke(1995), sem deixar de agregar outras contribuições teóricas, ainda que em escala menor. O construto, no entanto, permite que quaisquer contribuições e análises teóricas de atributos do comportamento dos agentes propostos possam ser utilizadas, uma vez que a sua generalidade permite que qualquer forma de conduta possa ser modelada. Para criar o simulador, foi utilizado o ambiente NetLogo 3.1.4 (WILENSKY,1999), cuja especial capacidade de tratar fenômenos que se desenvolvem ao longo do tempo se revelou de extrema valia e praticidade, demonstrando o potencial surgimento de propriedades imprevisíveis sob outras condições. Os testes realizados para aferir a aderência do modelo, usando distribuições de freqüência onde as médias e os desvios-padrão eram controlados, mostraram que o comportamento dos agentes em ambiente dinâmico segue perfeitamente os conceitos teóricos utilizados como referencial, especialmente no que diz respeito ao comportamento do agente cidadão, modelado segundo os conceitos estabelecidos por Becker (1968). Restou demonstrado pela simulação que a percepção de possibilidade de punição tem impacto considerável no comportamento dos agentes diante da possibilidade de delito, assim como o somatório dos custos envolvidos é grandemente influenciado pelo tipo de delito cometido, ou seja, a dimensão monetária do dano, pela estrutura de custos da justiça e pelos meios de arrecadação destes recursos utilizados pelo estado, conforme afirmam Shavell e Polinsky(2000). Para realizar os testes foram construídas tabelas de valores utilizando a técnica 2k fatorial, descrita em detalhes em Law e Kelton (1982), onde cada um dos fatores de impacto foi desdobrado em valores dinâmicos compatíveis com o sistema real no qual o modelo foi baseado. O teste 2k fatorial é especialmente prático no sentido de aferir a aderência do modelo de simulação aos referenciais teóricos, pois permite grande flexibilidade ao analista na formulação de cenários, pois quaisquer combinações possíveis podem ser testadas, dependendo apenas do foco de interesse da pesquisa. Foram realizados também experimentos de simulação com o objetivo de verificar o comportamento do modelo ao atingir estado estacionário. Os testes de validação e os experimentos de simulação apresentaram resultados semelhantes. / Agent based models are considered a milestone in social sciences research mainly because it´s random behavior allows a variety of analysis that were impossible to make due to complexity constraints. The term complexity usually derives from a large number of interacting and dynamic variables which cannot be considered singularly, in order to keep the whole system attributes intact. For that matter, a systemic approach of analysis is required not only to keep this information within the framework of the system, but also to get the emergent properties that come along with the interaction, not observable singling out one or another factor. Considering these very specific characteristics an agent based model, or multiagent system turns to be a more suitable analysis than analytic counterparts due to its very particular capability to deal simultaneously with a large number of random variables. This research work develops an agent based model to study the emergent properties due to public enforcement of criminal law. In order to build a model of criminal behavior suitable to the project, well known theories about crime were taken in consideration, including Shavell and Polinsky (2000), Shavell (2004), Becker (1968), Cohen and Felson(1979) and Clarke (1995). The construct derived from the analysis is a model based on three interacting agents and the environment, namely: citizen-opportunity-state-environment (offender-victim-public enforcement of law-environment). These interacting agents have the ability to define the boundaries of any given offence, regardless the kind, being capable to cope with from a simple misdemeanor (e.g. a traffic violation) to a felony (e.g. a capital murder). Particularly, the citizen-agent behavior has been modeled according to Becker’s (1968) Crime and Punishment seminal work. These remarks explains why originally the term citizen has been used in first place, and not offender. Rational Choice theory says that anyone can be an offender, depending on certain conditions. The State performance regarding its economic impact in terms of the public enforcement of criminal law is seen through the lenses of Shavell and Polinsky’s theory (2000). The opportunity and environment agents’ behavior are mainly based on Cohen and Felson’s Routine Activity Theory (1979) and Clarke’s Rational Choice Theory (1995). The model was implemented using NetLogo 3.1.4 platform (WILENSKY,1999). The simulation program provides a wide range of dynamic options making it very easy to perform any kind of test in order to assess the behavior of a given criminal rule in its dynamic operation. Validation and experimental tests were performed. The resulting responses were very consistent with the theoretical basis on which de model was based on.