Esta dissertação trata da perspectiva sobre a avaliação de investimentos em condições de risco, especificamente o caso da infraestrutura urbana, buscando-se um aperfeiçoamento dos métodos tradicionais de análise de investimento, através da aplicação do Método de Monte Carlo (MMC) e da criação de planilha eletrônica, possibilitando o refino e agilidade no processo de tomada de decisão quanto ao investimento em infraestrutura urbana. O objetivo principal deste trabalho consiste na promoção de critérios objetivos para a aceitação de projetos de investimento em infraestrutura urbana, cenário que apresenta inúmeras condições de risco. Para isto será desenvolvida uma ferramenta computacional, utilizando-se do MMC para refinar os métodos tradicionais. Para testar o modelo desenvolvido será utilizado o caso da ampliação física do Hospital de Clínicas de Porto Alegre. Assume-se a premissa que o emprego de métodos tradicionais como análise do VPL, payback e TIR não são suficientes para embasar as tomadas de decisão quando os investimentos estão em condições de risco. A partir das contribuições teóricas dos campos de estudos sobre Estatística, Engenharia Econômica, Análise de Investimentos e Tecnologia Computacional, propõe-se a elaboração de planilha eletrônica através de implementos em Visual Basic for Applications (VBA) e o uso do MMC objetivando expandir as possibilidades de análise existentes de investimento em condições de risco; situação esta encontrada em investimentos de infraestrutura, atingindo um aperfeiçoamento dos métodos tradicionais. Para isto, será realizada simulação em cenário real: a ampliação física do Hospital de Clínicas de Porto Alegre. Conclui-se que, enquanto a análise de investimento do caso em estudo pelo método do VPL embasou a viabilidade financeira de execução do investimento através de sua aceitabilidade, quando o modelo é aperfeiçoado através do MMC e princípios estatísticos, como por exemplo o Princípio de Pareto, o mesmo investimento demonstrou-se instável e concluiu-se que sua viabilidade não poderia ser garantida. Isto se deve ao fato de que, em projetos de infraestrutura, os riscos devem ser apreciados, cenário não identificado no uso de métodos tradicionais. / This dissertation deals with the perspective on the evaluation of investments in risky conditions, specifically the case of urban infrastructure, aiming at an improvement of the traditional methods of investment analysis, through the application of the Monte Carlo Method (MCM) and the creation of spreadsheet, enabling refining and agility in the decision-making process regarding the investment in urban infrastructure. The main objective of this work is the promotion of objective criteria for the acceptance of investment projects in urban infrastructure, a scenario that presents risk conditions. For such, a computational tool using the MCM to refine the traditional methods will be developed. To test the developed model, the case of the physical enlargement of the Hospital de Clínicas of Porto Alegre will be used. It is assumed that the use of traditional methods such as analysis of NPV, payback and IRR are not enough to support decision making when investments are in risky conditions. From the theoretical contributions of the field of studies on Statistics, Economic Engineering, Analysis of Investments and Computational Technology, it is proposed the elaboration of electronic spreadsheets through Visual Basic for Applications (VBA) tools and the use of MCM in order to expand the existing possibilities of investment in risky conditions, a situation found in investments in infrastructure, reaching an improvement of the traditional methods. For such, a simulation in a real scenario - the physical expansion of the Hospital of Clinics of Porto Alegre - will be done. It is concluded that, while the investment analysis of the case under review by the VPL method supported the feasibility of execution of the investment through its acceptability, when the model is improved through MCM and statistical principles such as the Pareto Principle, the same investment was proved unstable and it was concluded that its viability could not be guaranteed. This is due to the fact that in infrastructure projects the risks should be appreciated, an unidentified scenario in the use of traditional methods.

A água é uma das substâncias mais importantes em nosso planeta. São conhecidas atualmente 70 anomalias para a água. Uma das principais anomalias da água é a anomalia na densidade, através da qual a densidade da água aumenta com a diminuição da temperatura e atinge um máximo. É devido a ela que o gelo flutua em água e assim animais aquáticos sobrevivem quando a temperatura ambiente é muito baixa. Podemos citar também a anomalia na difusão, na qual a molécula de água se difunde mais rapidamente quando o sistema está mais denso. Acredita-se que as anomalias da água surgem devido a possibilidade de formação de ligações de hidrogênio entre suas moléculas. A combinação de ligações de hidrogênio e de interações de van der Walls leva à formação de dois tipos de estrutura: uma em que as moléculas formam quatro ligações de hidrogênio e que implica em moléculas mais distantes, e outra estrutura mais compacta sem tantas ligações de hidrogênio. Acredita-se que a competição entre estas duas estruturas seja responsável pela existência de anomalias na água. Nos modelos atomísticos, a direcionalidade das ligações de hidrogênio é considerada ingrediente importante para presença de anomalias. Além disso, as interações atrativas consideradas são igualmente importantes. Neste trabalho verificamos se a direcionalidade e as interações atrativas são fundamentais para a presença de anomalia. Estudou-se um modelo puramente repulsivo de duas escalas. O modelo estudado possui 3 diferentes fases para T = 0: duas fases ordenadas, uma na qual 1/4 da rede encontra-se ocupada e outra de 1/3 de ocupação, e a fase gasosa com rede vazia. Para temperatura não nula o sistema apresenta as duas estruturas ordenadas já observadas em temperatura nula e uma fase fluida desordenada. Foi obtido o diagrama de fases a potencial químico vs temperatura para o modelo. As transições de fase foram analisadas. Além disso foi observada a transição entre os comportamentos forte e fraco para a difusão, além das anomalias na densidade e no coeficiente de difusão para uma certa região do diagrama de fases. Observou-se que a linha de temperatura de máxima densidade engloba a região de anomalia na difusão no diagrama de fases, e estas anomalias encontram-se na região reentrante da linha de coexistência, próxima ao ponto bicrítico. / Water is one of the most important substances in our world. Nowadays 70 anomalies are known for water. One of the major anomalies of water is the density anomaly, in which water density increases with the decrease of temperature and reaches a maximum. Because of this anomaly ice floats on liquid water and aquatic animals can survive when room temperature is too low. We can also cite the anomaly in the diffusion, where a water molecule can spread faster when the system has higher density. It is believed that water anomalies rise from the possibility of hydrogen bonds interaction between two water molecules. The combination of hydrogen bonds and the van der Waals interactions leads to the formation of two kinds of structures: one of them where water molecules form four hydrogen bonds and implies on more distant molecules; the other structure is more compact and doesn’t have so much hydrogen bonds. It is believed that the two structures competition is responsible for the existence of water anomalies. In the atomistic models the directionality of hydrogen bonds is considered an important ingredient for the presence of anomalies. Furthermore, the attractive interactions considered are equally important. In this work we verified if the directionality and the attractive interactions are fundamental for the presence of anomalies. We studied a purely repulsive model with two scales. The studied model has three different phases when temperature T = 0: two ordered phases, one with a 1/4 of the lattice occupied and other with density 1/3 occupied sites, and a gas phase that is all empty. For temperature greater then zero the system presents the same two ordered structures and a fluid disordered phase. We obtained the temperature vs chemical potential phase diagram for this model. The phase transitions were analyzed. Moreover, we have observed the transition between the fragile and strong behaviors for diffusion, besides the density and diffusion coefficient anomalies for a certain region of the phase diagram. We have observed that the temperature of maximum density encompasses the diffusion anomaly region in phase diagram, and these anomalies are in the reentrant coexistence line region, near to the bicritical point.

Esta tese tem por objetivo avaliar a incerteza de medição no ensaio de fadiga do aço inoxidável AISI 316L em corpos de prova instrumentados ensaiados em temperatura ambiente em três níveis de rugosidade. Para isso, foi utilizado o método matemático de Monte Carlo para calcular o valor da incerteza com base nas variáveis de entrada. Para atingir esse objetivo foram realizados diversos ensaios de fadiga utilizando um número de corpos de prova suficientes para garantir confiabilidade estatística conforme norma ASTM E739-10 (ASTM, 2010), a uma frequência de ensaio de 20 Hz e razão de carregamento de 0,1. Nesta etapa 36 corpos de prova foram divididos em três grupos com valores de tensão 480, 500 e 520 MPa com três níveis de rugosidade distintos. Para o desenvolvimento do modelo matemático, cada fonte de incerteza teve distribuição de probabilidade determinada com base no comportamento estatístico característico de cada variável. Sendo assim, a partir de um modelo matemático desenvolvido foi possível determinar incerteza de medição do valor da tensão corrigida e do n° de ciclos, com uso do software Matlab®. Os resultados de incerteza de medição foram comparados aos valores obtidos pelo método Kragten, mostrando valores semelhantes entre os dois métodos matemáticos. / This thesis aims to evaluate the measurement uncertainty in the fatigue test of AISI 316L stainless steel in instrumented specimens tested at room temperature in three levels of roughness. For this, the mathematical method of Monte Carlo was used to evaluate the value of uncertainty based on the input variables. To achieve this goal several fatigue tests using a number of specimens sufficient to ensure statistical reliability according to ASTM E739-10 (ASTM, 2010) at test frequency of 20 Hz and stress ratio of 0.1 were tested. In this step 36 specimens were divided into three stress values of 480, 500 and 520 MPa with three different levels of roughness. For the mathematical development, each source of uncertainty probability distribution was initially determined based on the statistical behavior characteristic of each variable. A mathematical model was derived to determine the uncertainty of the stress value and number of cycles, using the Matlab® software. The combined uncertainty of results were compared to values obtained by Kragten method, showing similar values between the two mathematical methods.

A nova norma NBR 6118:2014, Projeto de estruturas de concreto – Procedimento, traz recomendações específicas para o dimensionamento de estruturas de concreto com resistência característica à compressão (fck) superior a 50 MPa, às quais não faziam parte de sua versão anterior NBR 6118:2007. Os pilares de concreto armado são as peças estruturais que mais se beneficiam com a utilização de concretos de elevada resistência, pois o seu emprego torna possível a redução da seção transversal, o emprego de menores taxas de armadura e um melhor desempenho estrutural. Então, o objetivo deste trabalho é avaliar a confiabilidade de pilares de concreto armado projetados de acordo com estas novas especificações para concreto de alta resistência. Inicialmente foi feita a validação de um modelo numérico de elementos finitos para representar o comportamento de pilares esbeltos de concreto armado. Após, considerando-se a variabilidade das resistências dos materiais e das dimensões da seção transversal, empregando-se o método de Monte Carlo é feita a determinação do valor médio, desvio padrão e distribuição de probabilidade da capacidade resistente destes pilares. Finalmente, combinando-se as propriedades da distribuição de probabilidade da capacidade resistente destes pilares, as propriedades estatísticas do carregamento, empregando-se o Método de Monte Carlo, foram determinados os índices de confiabilidade alcançados no projeto de diversos tipos de pilares, variando-se a resistência à compressão do concreto, a taxa de armadura, o índice de esbeltez e a relação entre as parcelas variável e permanente do carregamento. Desta forma, através de modelos matemáticos implementados computacionalmente, é possível analisar a confiabilidade de pilares de concreto projetados pela NBR-6118:2014 para diferentes situações. Os resultados obtidos indicam que o índice de confiabilidade está normalmente em níveis adequados, sendo que os casos mais críticos são os de pilares mais esbeltos e cuja parcela da carga acidental seja muito significativa. / The new Brazilian design code NBR 6118:2014, Design of concrete structures – Procedure, introduces new recommendations to the design of concrete structures with a concrete characteristic compressive strength (fck) higher than 50 MPa, which were not part of the previous version of this code NBR 6118:2007. The use of high-strength concrete is most advantageous in the design of concrete columns because it allows a reduction of the concrete cross sectional area, the use of lower steel ratio and a better structural performance. The objective of this study is to evaluate, through numerical modeling, the reliability reached in the design of reinforced concrete columns according to these new specifications for high strength concrete present in the new code. At first, a finite element model for slender reinforced concrete columns is validated. Then, once the variability of the materials strength and column’s dimensions are known, by using the Monte Carlo method, the mean value, standard deviation and the probability distribution of the column resistance are determined. Finally, by combining the statistical properties of the column resistance with the loading probability distribution and using the Monte Carlo method the reliability index for each designed column is calculated. Different types of columns are studied by varying the concrete characteristic compressive strength, the slenderness ratio, the first order eccentricity, the steel ratio and the live to dead load ratio. Therefore, through numerical models, it is possible to analyze the reliability for concrete columns in different situations designed by the NBR-6118:2014. In conclusion, the results achieved are normally in adequate reliability levels, and the critical cases are the long columns under condition of high live to dead load ratio.

A presente dissertação devolveu uma metodologia para determinar a probabilidade de escape de nêutrons conforme a energia e posição no reator. Para tanto, simulou-se um reator qualitativo semi-infinito de secção quadrada composto por três regiões homogêneas distintas. O domínio do reator foi subdividido em cem subcamadas concêntricas uniformes para a análise da fuga de nêutrons. Desta maneira, o nascimento de cada nêutron em cada camada foi registrado, bem como sua energia inicial. Os cálculos das trajetórias dos nêutrons foram efetuados usando o Método de Monte Carlo Físico. Assim, o código gerou a história paralela de 4x106 nêutrons, armazenando a energia final, posição final e fluxo angular na superfície. Desta forma, foi possível atribuir a probabilidade de escape de nêutrons provenientes de diferentes camadas conforme suas energias e posições iniciais. O método foi capaz de estabelecer o espectro de fuga, relações de dependência entre energia inicial e probabilidade de escape, além de observar que, sob certas condições, a probabilidade de escape possui crescimento exponencial ao longo do domínio. / This work presents a methodology to determine the neutron escape probability according to its energy and start position in the reactor. A semi-infinite qualitative reactor was simulated by a C++ program. This reactor has three distinct homogeneous regions. It was subdivided into hundred uniform concentric layers for a statistical analysis, allowing to record the birth and initial energy of each neutron in each layer. The neutron's path calculation was performed using Monte Carlo. The program generated 4x106 parallel neutron stories and has stored the final energy, position and angular flux. Thus, we determined the neutron escape probability from different layers. The method was able to estimate the leakage dependency with initial energy and position and it showed that the escape probability has a exponential growth tendency along the domain in certain conditions.

A aplicação de modelos matemáticos baseados em técnicas numéricas aumentou com o avanço da informática através da criação de microprocessadores mais velozes e periféricos de armazenamento de dados com grande capacidade. Buscando uma maior produtividade e a melhoria da qualidade final do produto solidificado, propõe-se neste trabalho, desenvolver um modelo computacional de origem físico para o crescimento de grão no Lingotamento Contínuo. O modelo é simulado utilizando o Método de Monte Carlo juntamente com o Método das Diferenças Finitas, com o objetivo de obter e caracterizar a transição colunar-equiaxial através desse método. As simulações foram realizadas utilizando a programação Fortran 90/95 no ambiente Linux através do software Developer Studio e aplicados nos aços SAE 1015 e 1020. A seguir, foram realizadas comparações entre as macroestruturas simuladas e as macroestruturas das amostras dos aços obtidas pelo Laboratório de Fundição (LAFUN) da UFRGS. Destas simulações observou-se que o modelo oferece a possibilidade da simulação de diferentes condições operacionais para prever a evolução macroestrutural de lingotes. / The application of mathematical models based on numerical techniques increased with the advancement of information technology by creating faster microprocessors and peripherals with high data storage capacity. Seeking for a major productivity and the improvement of the final quality of the solidified product. In this assignment, the development of a computational model with a physical origen for the growth of the continuous casting, is proposed. The model is simulated by using the Monte Carlo Method with the Finite Differences Method, the goal is to obtain and characterize the columnar-equiaxed transition through this method. The simulations were performed using the Fortran 90/95 programming in Linux environment using the Developer Studio software and applied to the steel SAE 1015 and 1020. Following it, comparisons were made between the macro and simulated macrostructures of the samples obtained by the Laboratory of Steel Casting - Laboratório de Fundição (LAFUN) at UFRGS. From these simulations, it was observed that the modified model offers the possibility of simulating different operating conditions to predict the evolution of macrostructural ingots.

Nesta tese realizamos um estudo de campo médio dos efeitos da desordem do campo cristalino sobre os diagramas de fases do modelo ANNNI com spin 1. Resultados de campo médio anteriores obtidos para o modelo com campo cristalino constante mostraram fases parcialmente desordenada em seus diagramas de fases. No entanto, a existência de tais fases não é garantida, devido às flutuações não levadas em conta na aproximação de campo médio. Sendo assim, estudamos o modelo com campo cristalino constante por meio de simulações de Monte Carlo. Nosso objetivo foi fornecer maior credibilidade aos resultados de campo médio obtidos anteriormente, e com isso adquirirmos maior confiança nos resultados de campo médio apresentados nesta tese para um campo cristalino aleatório. Os resultados das simulações mostraram evidências de que a configuração com planos desordenados da fase comensurável de período 6, observada em baixas temperaturas, pode sofrer uma transição interna para uma configuração na qual tais planos estão ausentes. Esse resultado está em concordância com os estudos de campo médio. Para o caso de um campo cristalino aleatório, investigamos os diagramas de fases do modelo na aproximação de campo médio. Para uma distribuição de probabilidades genérica do campo cristalino, obtivemos o funcional da energia livre e as equações de campo médio, assim como as expressões das linhas críticas e dos pontos tricríticos do modelo. Para investigarmos os efeitos da desordem sobre os diagramas de fases, consideramos as distribuições de probabilidade do campo cristalino delta-bimodal e gaussiana. / In this thesis we present a mean-field study of the effects of the disorder of the single-ion anisotropy on the phase diagrams of the spin-l ANNNI model. Previous mean-field results for the model with constant single-ion anisotropy showed partially disordered phases in the phase diagrams of the model. H owever , the existence of such phases cannot be taken as granted, due to the fiuctuations not taken into account in the mean-field approximation. Thus, we studied the model with constant single-ion anisotropy by means of Monte Carlo simulations. Qur aim was to give credit to previous mean-field results and to get more confidence in the mean-field results showed in this thesis obtained for a random single-ion anisotropy. The results obtained from the simulations showed evidences that the configuration of the commensurate phase with period 6 displaying disordered planes, found at low temperatures, may undergo an internaI phase transition to a phase in which such planes are absent. This result is in agreement with previous mean-field studies. For a random single-ion anisotropy, we investigated the phase diagrams of the model within mean-field approximation. For a generic probability distribution of the anisotropy, we obtained the free-energy functional and the mean-field equations, as well as the expressions for the criticallines and tricritical points of the model. In order to investigate disorder effects, we obtained phase diagrams of the model by considering a double-delta and a Gaussian probability distribution of the single-ion anisotropy.

O estudo da polarização eletrostática em líquidos homogêneos é crucial para compreender as propriedades termodinâmicas e de solvatação. Estes efeitos de polarização são também importantes para estudar a dinâmica molecular de reações químicas, conformações (tais como proteínas), e do meio solvente (solvatocromismo e RMN). Nesta dissertação, estudou-se a variação do momento de dipolo da fase gasosa para a fase líquida das moléculas: acetona, acetonitrila, amônia e água. Os momentos de dipolo em fase gasosa foram obtidos com cálculos quânticos ab initio para a molécula isolada. Para a fase líquida, cálculos equivalentes foram feitos sobre estruturas geradas de simulação computacional clássica (Monte Carlo). Tal procedimento foi feito segundo a metodologia \"Monte Carlo / Mecânica Quântica - Seqüencial\". Este método foi aplicado iterativamente, tornando possível estudar a variação do momento de dipolo e das propriedades termodinâmicas devido à mútua polarização das moléculas do solvente. Usando os resultados calculados, uma discussão foi realizada sobre a polarização eletrostática e o dipolo induzido em fase líquida. / The study of the electrostatic polarization in homogeneous liquids is crucial to the understanding of thermodynamic and solvation properties. These polarization effects are also important to study molecular dynamics of chemical reactions, conformations (such as proteins), and solvent environment (solvatocromism and NMR). In this dissertation, the variation of dipole moments from the gas phase to the liquid phase of acetone, acetonitrile, ammonia and water, were studied. The dipole moments in the gas phase were obtained with quantum mechanical ab initio calculation on the isolated molecule. For the liquid phase, equivalent quantum mechanical calculations were made on structures generated from classical computer simulations (Monte Carlo). Such procedure was made according to the methodology \"Sequential - Monte Carlo/Quantum Mechanics\". This method was applied iteratively, becoming possible to study the variation of the dipole moment and the thermodynamic properties due to the mutual polarization of the solvent molecules. Using the calculated results, a discussion was made of the electrostatic polarization and the induced dipole moment in the liquid phase.

O presente trabalho teve por objetivo avaliar a influência dos custos de transação em operações de hedge. O estudo se diferencia das tradicionais pesquisas relacionadas ao cálculo da razão ótima de hedge à medida que define os custos de transação em mercados futuros de maneira mais ampla. Considerou-se não apenas as taxas cobradas pela bolsa e corretora, mas também o custo de oportunidade do dinheiro destinado aos ajustes diários e a tributação incidente em mercados futuros. Como consequência dessa abordagem mais ampla, argumenta-se que os custos podem assumir caráter estocástico por dependerem da trajetória de preços a ser verificada na operação de hedge. Tal caracterização faz com que haja uma incerteza sobre seu verdadeiro valor, antes de iniciada a operação. Os cálculos dos custos de transação foram realizados através de simulação Monte Carlo, em que foram geradas diversas trajetórias de preços futuros e, para cada uma delas, calculados os respectivos custos de transação. Dessa forma, a incerteza em relação ao custo tornou-se passível de mensuração, sendo representada pelo desvio-padrão de sua distribuição de probabilidade e denominada \"risco de custo\". Com base nas estatísticas relacionadas aos custos de transação, adotou-se um modelo de decisão baseado no paradigma da utilidade esperada, de forma a gerar as razões ótimas de hedge para diferentes cenários. Os resultados foram consistentes em relação a estudos anteriores, ao verificar um impacto negativo dos custos de transação na razão ótima de hedge. Quanto à influência do risco de custo na decisão de um produtor rural sobre quanto negociar em bolsa, a influência mostrou-se pouco significativa. / The objective of the current research is to evaluate the influence of transaction costs on hedging decisions. This study expands on previous research by considering a broader and more realistic range of costs involved in hedging operations, such as income tax on both spot and futures positions. In addition to the traditional approach of investigating how the magnitude of transaction costs affect optimal hedge ratio, this study also explores the uncertainty related to those costs. The exact amount of some transaction costs in futures markets are not known by the time the hedge is placed, because they depend on the random trajectory of futures prices during the hedging period. Thus the hedger will only know the exact cost of his hedge at the end of the hedging period, which creates a \"cost risk\". This uncertain portion of transaction costs has never been addressed in the hedging literature, which typically assumes all costs are known to the hedger. This study incorporates uncertainty in transaction costs with Monte Carlo simulation, which generates a distribution of price series that is used to calculate expected values and variances of transaction costs. These values are then used in an expected utility framework to generate optimal hedge ratios in different scenarios. Results are consistent with previous studies in that the magnitude of transaction costs has a negative impact on optimal hedge ratios. Findings also suggest that \"cost risk\" has little influence on hedge ratios.

A confiabilidade de estruturas apresenta sólidos desenvolvimentos teóricos e crescentes aplicações práticas. Durante os últimos anos, avanços significativos foram obtidos em termos dos métodos de transformação (FORM, SORM), bem como em termos das técnicas de simulação de Monte Carlo. Métodos de transformação se mostraram eficientes para problemas de dimensões e não-linearidades moderadas. Já técnicas de simulação sempre permitiram a solução de problemas de grandes dimensões e fortemente não lineares, embora o custo computacional possa ser uma séria limitação. Com o avanço da capacidade de processamento dos computadores e com o desenvolvimento de técnicas de amostragem inteligente, a simulação de Monte Carlo passa a ser cada vez mais viável. Este trabalho tem por objetivo estudar e programar em computador técnicas de amostragem inteligente em simulação de Monte Carlo. O StRAnD é um programa de computador que já possui implementadas as técnicas de simulação de Monte Carlo Bruto e com Amostragem por Importância, ambas utilizando a Amostragem Simples na geração das variáveis básicas. Assim, são adicionadas, ao StRAnD, as técnicas de Amostragem Assintótica, Amostragem Melhorada e Simulação de Subconjuntos. Além disso, são programadas as técnicas de Amostragem por Hipercubo Latino e Amostragem por Variáveis Antitéticas. Nesta dissertação, são analisados seis problemas distintos, de forma que as vantagens e desvantagens de cada técnica sejam avaliadas, em termos da probabilidade de falha, do coeficiente de variação da probabilidade de falha, do erro relativo da probabilidade de falha e do tempo de processamento. / The structural reliability presents solid theoretical developments and increasing practical applications. During the past few years, significant advances were achieved in terms of transformation methods (FORM and SORM), as well as, in terms of Monte Carlo Simulation. Transformation methods are effective in problems with moderate dimensions and moderate nonlinearities. On the other hand, simulation techniques can be used to solve high-dimensional problems and highly nonlinear problems, although the computational cost could be a serious limitation. With the progress of computer processing capacity and with the development of intelligent sampling techniques, the Monte Carlo Simulation becomes increasingly feasible. This work aims to study and program intelligent sampling techniques in Monte Carlo simulation. The StRAnD (Structural Reliability Analysis and Design) software already has Crude Monte Carlo and Importance Sampling Monte Carlo, both using Simple Sampling as basic samples generator. Thus, the Asymptotic Sampling technique, the Enhanced Sampling technique and the Subset Simulation were added to the software. Moreover, the Latin Hypercube Sampling technique and the Antithetic Variates techniques were also added to the software. Six problems were evaluated in order to evaluate the advantages and disadvantages of each technique, in terms of probability of failure, coefficient of variation of the probability of failure, relative error and processing time.