Submitted by Claudinei Pracidelli (cpracide@ipen.br) on 2016-06-22T12:55:19Z No. of bitstreams: 0 / Made available in DSpace on 2016-06-22T12:55:19Z (GMT). No. of bitstreams: 0 / Dissertação (Mestrado em Tecnologia Nuclear) / IPEN/D / Instituto de Pesquisas Energeticas e Nucleares - IPEN-CNEN/SP

Submitted by Philipp Herbst Minarelli null (hminareli@hotmail.com) on 2016-11-09T21:28:11Z No. of bitstreams: 1 Philipp_Defesa do mestrado_Definitivo.vf2.pdf: 1053422 bytes, checksum: c6a949c72f961e34200d0497370c1ee3 (MD5) / Rejected by LUIZA DE MENEZES ROMANETTO null (luizaromanetto@hotmail.com), reason: Recusado a pedido do autor. Solicitamos que realize uma nova submissão seguindo as orientações abaixo: A data (mês) que consta na capa e folha de rosto do trabalho deve ser a mesma que consta na folha de aprovação. Incluir o processo de financiamento nos agradecimentos da dissertação/tese. Corrija estas informações e realize uma nova submissão com o arquivo correto. Agradecemos a compreensão. on 2016-11-16T11:35:33Z (GMT) / Submitted by Philipp Herbst Minarelli null (hminareli@hotmail.com) on 2016-11-16T14:56:55Z No. of bitstreams: 1 Philipp_Defesa do mestrado_Definitivo.vf3.pdf: 1054510 bytes, checksum: 18a84d22cc65a0da2332b20b03843a00 (MD5) / Approved for entry into archive by Felipe Augusto Arakaki (arakaki@reitoria.unesp.br) on 2016-11-23T12:12:39Z (GMT) No. of bitstreams: 1 minarelli_ph_me_bot.pdf: 1054510 bytes, checksum: 18a84d22cc65a0da2332b20b03843a00 (MD5) / Made available in DSpace on 2016-11-23T12:12:39Z (GMT). No. of bitstreams: 1 minarelli_ph_me_bot.pdf: 1054510 bytes, checksum: 18a84d22cc65a0da2332b20b03843a00 (MD5) Previous issue date: 2016-09-09 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / A inclusão das incertezas em um modelo de análise apresenta informações relevantes aos fornecedores de cana de açúcar para a região de Jaú/SP e ganha importância em estudos onde são comparados vários projetos de investimento, pois permite o uso de mais uma ferramenta para a tomada de decisão. A análise de indicadores econômicos para projetos de investimento no setor visa identificar a competitividade dos custos de produção de cana de açúcar. Tendo em vista as grandes incertezas que permeia o setor sucroenergético, como política interna, produtividade, qualidade da matéria prima. O trabalho teve como objetivo determinar, mediante o cálculo do valor presente líquido e da taxa interna de retorno, a viabilidade da produção de cana de açúcar na região de Jaú/SP, na condição de fornecedor, bem como identificar, usando o método de Monte Carlo, o risco dessa atividade. Os indicadores econômicos avaliados oscilaram em função de quatro variáveis críticas: (1) produtividade da cana de açúcar em diferentes estágios de cortes; (2) qualidade da matéria prima, medida em ATR; (3) preço pago ao fornecedor; (4) taxa de juros. As funções de distribuição de probabilidade dessas variáveis críticas foram estimadas e incluídas nas equações de VPL e TIR e, utilizando o método de simulação de Monte Carlo, foi estimado o seguinte conjunto de indicadores: análise de sensibilidade dos indicadores em relação às variáveis críticas, variabilidade dos valores de VPL e TIR e, por fim, mapeamento de risco para os investidores em fornecer cana de açúcar segundo o modelo proposto nesse trabalho para o EDR de Jáu. Esses indicadores permitiram desenhar cenários econômicos associados à sua probabilidade de ocorrência. Os resultados mostram que o fornecimento de cana de açúcar na região de Jaú/SP é, para o modelo adotado, uma atividade de alto risco sendo o preço pago ao fornecedor, a taxa de juros (SELIC) e a qualidade da matéria prima os maiores ofensores para que se atinjam níveis de rentabilidades satisfatórias, sendo que as rentabilidades positivas foram encontradas para probabilidades menores do que 40%. A variável com maior efeito sobre a rentabilidade foi o preço de venda do produto, e também importantes o custo de corte, carregamento e transporte. / The inclusion of uncertainty in an analysis model presents information relevant to sugar cane suppliers to the region of Jaú(São Paulo) and gains importance in studies where they compared several investment projects, as it allows the use of another tool for making decision. The analysis of economic indicators for investment projects in the sector aims to identify the competitiveness of sugar cane production costs. Given the great uncertainty that pervades the sugarcane industry, as domestic policy, productivity, quality of raw material. The study aimed to determine, by calculating the net present value and internal rate of return, the viability of sugarcane production in region of Jaú of grower status under risk conditions, using Monte Carlo simulation. The assessed economic indicators fluctuated against four critical variables: (1) productivity of sugar cane at different ages; (2) quality of raw material (TRS - Total Recoverable Sugar); (3) price paid to the supplier; (4) interest rate. The probability of these critical variables distribution functions were estimated and included in the Net Present Value (NPV) and Internal Return Rate (IRR) equations and using the Monte Carlo simulation, the following set of indicators were estimated: the indicators sensitivity analysis in relation to the critical variables, variability NPV and the IRR values and finally risk mapping for investors to provide cane sugar based on the model in this work for the region of Jaú. These indicators allowed to draw economic scenarios associated with their probability of occurrence. The results show that the grower of sugar cane in the region of Jaú/SP is, for the adopted model, a high-risk activity and the price paid to the supplier, the interest rate (SELIC) and the quality of the raw material the largest offenders order to achieve satisfactory profitability levels, wherein the positive benefits were found to lower probability than 40%. The variable with the greatest effect on profitability was the selling price of the product, and also the important cost cutting, loading and transportation.

Neste trabalho o transporte de nêutrons é simulado em materiais presentes no núcleo de reatores. O espectro de nêutrons é decomposto como uma soma de três distribuições de probabilidade. Duas das distribuições preservam sua forma com o tempo, mas não necessariamente sua integral. Uma das duas distribuições é devido ao espectro de fissão, isto é, altas energias de nêutrons, a outra é uma distribuição de Maxwell-Boltzmann para nêutrons de baixas energias (térmicos). A terceira distribuição tem uma forma a priori desconhecida e que pode variar com o tempo, sendo determinada a partir de uma simulação Monte Carlo com acompanhamento dos nêutrons e suas interações com dependência contínua de energia. Isto é obtido pela parametrização das seções de choque dos materiais do reator com funções contínuas, incluindo as regiões de ressonâncias resolvidas e não resolvidas. O objetivo deste trabalho é implementar efeitos de up-scattering através do tratamento estat ístico da população de nêutrons na distribuição térmica. O programa de simulação calcula apenas down-scattering, pois o cálculo do up-scattering microscópico aumenta signi_cativamente tempo de processamento computacional. Além de contornar esse problema, pode-se reconhecer que up-scattering é dominante na região de energia mais baixa do espectro, onde assume-se que as condições de equilíbrio térmico para nêutrons imersos em seu ambiente são válidas. A otimização pode, assim, ser atingida pela manutenção do espectro de Maxwell- Boltzmann, isto é, up-scattering é simulado por um tratamento estatístico da população de nêutrons. Esta simulação é realizada utilizando-se dependência energética contínua, e, como um primeiro caso a ser estudado assume-se um regime recorrente. As três distribuições calculadas são então utilizadas no código Monte Carlo para calcular os passos Monte Carlo subsequentes. / In this work the neutron transport is simulated in reactor core materials. The neutron spectrum is decomposed as a sum of three probability distributions. Two of the distributions preserve shape with time but not necessarily the integral. One of the two distributions is due to prompt ssion, i.e. high neutron energies and the second a Maxwell-Boltzmann distribution for low (thermal) neutron energies. The third distribution has an a priori unknown and possibly variable shape with time and is determined from a Monte Carlo simulation with tracking and interaction with continuous energy dependence. This is done by the parametrization of the material cross sections with continuous functions, including the resolved and unresolved resonances region. The objective of this work is to implement up-scattering e ects through the treatment of the neutron population in the thermal distribution. The simulation program only computes down-scattering, for the calculation of microscopic upscattering increases signi cantly computational processing time. In order to circumvent this problem, one may recognize that up-scattering is dominant towards the lower energy end of the spectrum, where we assume that thermal equilibrium conditions for neutrons immersed in their environment holds. The optimization may thus be achieved by the maintenance of the Maxwell-Boltzmann spectrum, i.e. up-scattering is simulated by a statistical treatment of the neutron population. This simulation is performed using continuous energy dependence, and as a rst case to be studied we assume a recurrent regime. The three calculated distributions are then used in the Monte Carlo code to compute the Monte Carlo steps with subsequent updates.

Au cours des dernières décennies, la recherche de nouveaux matériaux multiferroïques nanostructurés avec des propriétés optimisées a conduit à l'élaboration d'une grande variété de modèles théoriques et des approches de simulation. Allant des modèles ab initio capables de décrire les propriétés à la température nulle des composés artificiels monocristallins à des approximations phénoménologiques pour la description des composites à la mésoscopique, ces recherches ont soulevé la question fondamentale de la relation entre la géométrie de la structure des systèmes hétérogènes et les propriétés des leurs transitions de phase. Cependant, malgré des progrès significatifs en la matière,cette question n'a pas encore été élucidée et les relations entre les modèles à différentes échelles ne sont pas entièrement distingués. La présente étude est consacrée à lier l’ensemble des modèles décrivant les matériaux nanocomposites multiferroïques à différentes échelles. Tout d'abord, nous présentons un développement méthodologique de l'approche Hamiltonien effectif couramment utilisé pour étudier les transitions de phase structurales. Les modifications introduites permettent d'étendre cette méthode pour prédire les propriétés à la température finie des systèmes hétérogènes. Le modèle construit est ensuite utilisé pour étudier les propriétés des nanostructures et solutions solides (BiFeO3)(BaTiO3). Recourant à des simulations Monte-Carlo, nous montrons que notre modèle fournit des résultats qui sont en ligne avec les observations expérimentales récentes et qu’il permet de prédire théoriquement les propriétés d'une large gamme de systèmes avec différentes géométries composites. La deuxième partie de l'étude consiste en l'application de la théorie de Ginzburg-Landau des transitions de phase à l’étude des propriétés des multicouches ferroélectriques et ferromagnétiques avec des interfaces épitaxiales. Plus précisément, nous décrivons théoriquement l’effet magnétoélectrique exhibé par les hétérostructures autonomes Pb(Zr0.5 Ti0.5) O3-FeGaB et BaTiO3-FeGaB. Enfin, nous montrons que la géométrie multicouche d'un nanocomposite ferroélectrique et ferromagnétique ouvre la voie à une amélioration radicale du signal de charge de sortie. / During past decades, the search for new nanostructured multiferroic materials with optimized properties has lead to the development of a vast variety of theoretical models and simulation approaches. Spreading from first principles based models able to describe zero-temperature properties of artificial single crystal compounds to phenomenological approximations for composites with mesoscale morphology, these investigations have raised the fundamental question of how the geometry of the structure affects the properties of phase transitions exhibited by heterogeneous systems. However, despite significant progress, the answer to this question still lacks clarity and the bridge connecting models at different scales is not fully constructed. The current study is devoted to linking together models of multiferroic nanocomposite materials applicable at different scales. First, we present a methodological development of effective Hamiltonian approach commonly used to study structural phase transitions. The introduced modifications allow to extend this widely used method to predict finite-temperature properties of compositionally heterogeneous systems. The constructed model is then used to study properties of (BiFeO3)(BaTiO3) nanostructures and solid-solutions. Resorting to Monte-Carlo simulations, we show that our model provides results that are in-line with recent experimental observations and allows to theoretically predict properties of a wide range of systems with different composite geometries. The second part of the study consists inapplication of Landau theory of phase transitions to investigate the properties of ferroelectric-ferromagnetic multilayerswith epitaxial interfaces. Specifically, we theoretically describe the strain-mediated direct ME effect exhibited byfree-standing Pb(Zr0.5 Ti0.5 )O3 -FeGaB and BaTiO3 -FeGaB heterostructures. Finally, we show that the multilayer geometry of a ferroelectric-ferromagnetic nanocomposite opens the way for a drastic enhancement of the output charge signal.

De nombreuses études, aussi bien expérimentales que théoriques, ont été menées pour l'industrie électrique afin de comprendre et d'être capable de prédire les mécanismes intervenant lors d'une coupure électrique par arc. Ces études ont montré la grande diversité et complexité des phénomènes physiques et chimiques mis en œuvre. L'arc électrique créé juste après la séparation des contacts est poussé par les forces électromagnétiques vers une zone d'extinction. Durant son développement et sa propagation, l'énergie lui est fournie par effet Joule et est dissipée par différents modes de transfert thermique. Son expansion rapide induit des effets de compressibilité avec la propagation d'ondes de pression. Le rayonnement intense du plasma créé, ainsi que les phénomènes aux pieds des électrodes, induisent une ablation des parois qui change la composition chimique du milieu et rend plus complexe la modélisation de l'ensemble des phénomènes couplés par hydrodynamique, électromagnétisme, transferts thermiques et diffusion d'espèces chimiques en régime fortement instationnaire. Le rôle du transfert radiatif est primordial dans la mesure où, d'une part, il conditionne le champ de température et donc les propriétés de transport, notamment électriques, et, d'autre part, participe pour une grande part à la thermo-dégradation des parois.La modélisation de l'ensemble des phénomènes physico-chimiques a beaucoup progressé durant les vingt dernières années mais le calcul du transfert radiatif demeure un point bloquant dans l'avancement des méthodes de modélisation.En effet, le champ de rayonnement est caractérisé par une luminance qui dépend de la longueur d'onde, de la position spatiale, de la direction de propagation, ainsi que du temps au travers ici des variations des champs de température et de la composition chimique.La prise en compte rigoureuse de toutes ces dépendances demeure inaccessible à l'heure actuelle pour des simulations complètes d'extinction d'arc, en particulier à cause de la complexité des spectres d'émission et d'absorption des milieux plasmas. Le recours à des modèles approchés pour le traitement spectral et/ou pour les dépendances géométriques et directionnelles est nécessaire.Ce travail a été mené en collaboration entre le Laboratoire EM2C et la société Schneider Electric. Il fait suite à des travaux de collaboration antérieurs dont l'objectif était la détermination des propriétés radiatives fondamentales des plasmas d'arc (Thèse S. Chauveau).Le but principal du présent travail est de développer des modèles approchés mais précis, ainsi que des outils de simulation numérique avec différents degrés de finesse, pour le calcul du champ de puissance radiative et des flux pariétaux dans des chambres de coupure électrique basse tension. Les modèles et outils développés doivent être implémentés dans des codes de simulation hydrodynamique dédiés soit à des géométries bidimensionnelles simplifiées, soit à des géométries industrielles complexes et tridimensionnelles. Ces outils doivent Les études expérimentales menées sur des maquettes représentatives des dispositifs industriels montrent la présence de vapeurs plastiques et métalliques dues à l'ablation et à l'érosion des parois. La composition chimique change par ailleurs fortement avec le temps entre la naissance et l'extinction de l'arc électrique. Bien que les modèles développés ici puissent être adaptés à cette composition chimique complexe, nous supposerons tout le long de ce travail que le milieu gazeux est un plasma d'air à l'équilibre thermique et chimique local. / Non fourni.

Apresentamos dois autômatos celulares com regras de interação locais que permitem descrever a dinâmica de população de um sistema predador-presa. Os modelos são definidos sobre uma rede regular quadrada e se diferenciam pelo caráter isotrópico ou anisotrópico da interação entre os sítios. A cada sítio é associada uma variável estocástica, que pode assumir três estados - vazio, presa ou predador. A dinâmica de competição entre espécies animais que nos interessa é a mesma descrita pelo modelo de Lotka-Volterra no qual as populações de presas e predadores oscilam temporalmente. Nosso objetivo é a análise dessas oscilações, como se comportam com o aumento da rede e se permanecem estáveis. Para a obtenção das séries temporais realizamos simulações de Monte Carlo. Para o autômato definido sobre o espaço isotrópico, também realizamos análise de campo médio dinâmico. Os resultados indicam que a oscilação é um efeito local (não sobrevive em sistemas infinitos), e é mais significativo devido à migração das espécies pelos subsistemas. O estudo da anisotropia revela alguns padrões espaciais organizados e que as oscilações são menos intensas do que no caso isotrópico e como consequência a fase ativa é mais abrangente. / We present two cellular automata with local interaction rules which allow us to describe the dynamical population of a predator-prey system. The models are defined on a regular square lattice and are distinguished by the isotropic or anisotropic character of the interaction between sites. To each site a stochastic variable is associated, which can assume three states- void, prey or predator. The competition dynamics between animal species which interest us is the same described by the Lotka-Volterra model in which the populations of preys and predators oscillate in time. Our aim is the analysis of these oscillations, how they behave with an increasing lattice and if they remain stable. In order to obtain temporal series we perform Monte Carlo simulations. For the automaton defined on isotropic space, dynamical mean field analysis was also performed. Results indicate that the oscillation is a local effect ( vanishing in infinite systems), and is more significant due to migration of species through the subsystems. The study of anisotropy reveals some organized spatial patterns and that oscillations are less intense than in the isotropic case and as a consequence the active phase is more comprehensive.

O uso de nanomateriais em aplicações como catálise e medicina, despertou nos últimos anos o interesse no estudo das propriedades de nanoclusters. O estudo das propriedades termodinâmicas desses sistemas é essencial, pois mudanças estruturais originadas de mudanças de fase podem alterar propriedades como atividade catalítica, momento magnético e propriedades óticas. A dinâmica molecular vêm sendo utilizada para o estudo computacional das propriedades termodinâmicas de diversos nanomateriais, enquanto o uso de métodos de Monte Carlo (MC), nesse contexto, tem se restringido ao estudo de nanoclusters de Lennard-Jones (LJ). Para avaliar a viabilidade do uso de métodos de MC no estudo de propriedades de sistemas reais, uma implementação do método Parallel Tempering Monte Carlo (PTMC) utilizando algoritmos do estado da arte para realização de trocas, determinação de temperaturas e ajuste de deslocamentos foi construída. Através de testes, é mostrado que alguns dos algoritmos implementados podem não ser adequados ao estudo do problema em questão. A implementação foi validada com o estudo das propriedades termodinâmicas de nanoclusters de LJ com 38, 55 e 147 átomos, que possuem resultados conhecidos na literatura. Além disso, resultados para propriedades do nanocluster LJ98 são apresentados, e devido à características estruturais desse sistema, é observada uma transição sólido-sólido entre as estruturas tetraédricas e icosaédricas em temperatura abaixo da de fusão. A hipótese do uso do algoritmo PTMC para o estudo de propriedades de materiais reais, foi testada nas nanoligas (PtCo)55 e (PtNi)55, descritas pelo potencial de Gupta. Através da comparação das estruturas de mais baixa energia com resultados de teoria do funcional da densidade (DFT, do inglês), é mostrado que o uso do potencial de Gupta pode ser justificado, dados os baixos desvios no comprimento de ligação (menores que 2.4%) e a semelhança de outras características estruturais. Os resultados obtidos indicam que o método PTMC é capaz de identificar as mudanças de fase das nanoligas estudadas. Essas mudanças são ilustradas e analisadas com o uso de um algoritmo para a comparação da similaridade de estruturas, com o qual foi possível analisar a fusão dos nanoclusters Co55, Ni55, Pt30Co25 e Pt40Ni15 (obtida em temperaturas entre 900 e 1100 K); além da fusão, a 727 K, e transição sólido-sólido, a 300 K, para a Pt55. Com as estruturas mais frequentes, encontradas através da análise de similaridade, e com a realização de cálculos DFT, foi possível observar um deslocamento do centro da banda d em direção ao HOMO (Highest Occupied Molecular Orbital) causado pelo aumento da temperatura. Esse deslocamento, segundo o modelo da banda d válido para superfícies, pode indicar uma maior reatividade dos nanoclusters nesses casos. / The use of nanomaterials in applications such as catalysis and medicine, aroused in the last years interest in studying properties of nanoclusters. The study of thermodynamic properties of these systems is essential, since structural changes originated from phase changes can alter properties such as catalytic activity, magnetic moment and optical properties. Molecular dynamics have been used for the computational study of thermodynamic properties of various nanomaterials, while the use of Monte Carlo methods (MC), in this context, has been restricted to the study of Lennard-Jones (LJ) nanoclusters. To evaluate the feasibility of using MC methods to study properties of real systems, an implementation of the Parallel Tempering Monte Carlo (PTMC) method using state of the art algorithms to perform exchanges, determine the temperature set and adjust the maximum displacement, was built. Through testing, it is shown that some of the implemented algorithms may not be suitable for the study of the problem in question. The implementation was validated by studying the thermodynamic properties of LJ nanoclusters with 38, 55 and 147 atoms, which have results known in the literature. In addition, results for the properties of the LJ98 nanocluster are reported, and due to the structural features of this system, a solid-solid transition between the tetrahedral and icosahedral structures in a temperature below melting is observed. The possibility of using the PTMC algorithm in the study of properties of real materials, is tested in the (PtCo)55 and (PtNi)55 nanoalloys, described by the Gupta potential. By comparing the lowest energy structures with density functional theory (DFT) results, it is shown that the use of the Gupta potential can be justified, given the small deviation in the bond lenght (less than 2.4%) and the similarity of other structural features. The results indicate that the PTMC method is able to identify the phase changes in the studied nanoalloys. These changes are illustrated and analyzed with the use of an algorithm for comparing the structure similarity, which made possible the analysis of the melting of the Co55, Ni55, Pt30Co25 e Pt40Ni15 nanoclusters (obtained at temperatures between 900 e 1100 K); and the melting at 727 K, and solid-solid transition at 300 K, for Pt55. With the most frequent structures, obtained by the similarity analysis, and through DFT calculations, it was possible to observe a shift in the d band center to the HOMO (Highest Occupied Molecular Orbital) caused by the temperature increase. This shift, following the d band model valid for surfaces, may indicate a higher reactivity of the nanocluster in these cases.

Censura intervalar resulta quando os tempos de sobrevivência não são exatamente conhecidos, sabe-se apenas que eles ocorreram dentro de um intervalo. Dados de sobrevivência agrupados são casos particulares de censura intervalar quando os indivíduos são avaliados nos mesmos intervalos de tempo, ocasionando um grande número de empates. Um procedimento comum para a análise desse tipo de dados é ignorar a natureza de censura intervalar dos dados, ou seja, tratar a variável aleatória tempo como contínua e assumir que o evento ocorreu no início, no ponto médio ou no final do intervalo e, então, usar um método padrão de análise de sobrevivência. Neste estudo, simulações de Monte Carlo, com o modelo de Weibull, foram realizadas para comparar esses três procedimentos e um método novo proposto que é uma combinação desses três métodos e é orientado pela observação do histograma do tempo versus a freqüência de cada intervalo para a decisão de qual valor a ser usado. Considera-se também a análise dos dados como censura intervalar. Os resultados mostram que analisar os dados exatamente como censura intervalar é a forma correta. Entretanto, quando a taxa de falha aumenta o ponto médio poderia ser usado. A natureza discreta dos tempos de falha deve ser reconhecida quando existe um grande número de empates. Métodos de regressão para tratar dados agrupados são apresentados por Lawless (2003) e Collett (2003), cuja estrutura é especificada em termos da probabilidade de um indivíduo falhar em um intervalo, condicionada à sua sobrevivência ao intervalo anterior. Os modelos considerados na literatura são o de riscos proporcionais de Cox ou o logístico. O modelo de Weibull é proposto, neste trabalho, como uma alternativa ao modelo de Cox para ajustar dados de sobrevivência com censura intervalar no contexto de modelos discretos. Através de simulações foram construídas as estatísticas da razão de verossimilhança e do teste escore para a discriminação entre esses dois modelos. Para ilustrar as simulações duas aplicações em dados agronômicos foram utilizadas. / Interval-censored results when survival times are not exactly known, knowing only that they occur in an interval. Grouped survival data are particular cases of intervalcensored when individuals are evaluated in the same time-intervals, causing a great number of ties. A common procedure for the analysis this type of data is to ignore the nature of interval-censored data, or rather, treat the random variable time as continuous, and assume that the event occurred in the beginning, midpoint or interval end, and then use a standard method of survival analysis. In this study, Monte Carlo simulations according to Weibull model, were performed in order to compare these three procedures and a new method proposed which is a combination of the three, and is directed by the observation of time histogram versus each interval frequency in order to decide which value be used. Interval-censored data is also considered. The results show that to analyse the data exactly as interval-censored is the correct form. However, when the failure rate increase the midpoint could be used. The discrete nature of failure time must be recognized when there are a great number of ties. Regression methods to treat grouped data are presented by Lawless (2003) and Collett (2003), whose structure is specified in terms of the probability of an individual failing in an interval, conditioned to his survival to previous interval. The models considered in literature are either those of Cox proportional hazards or the logistic one. Weibull model is proposed in this study as an alternative to Cox model in order to adjust survival data with interval-censored in the context of discrete models. Through simulations were built the statistics of ratio likelihood and score test to distinguish between these two models. To illustrate the simulations two applications in agronomy data were used.

O conceito de tabuleiros laminados de madeira protendidos transversalmente foi usado inicialmente no Canadá na década de 70. Desde então, foi largamente utilizado em um número crescente de países. No Brasil, esse sistema foi utilizado pela primeira vez com a construção da ponte sobre o rio Monjolinho, na região metropolitana de São Carlos, estado de São Paulo. A importância deste sistema estrutural requer um maior conhecimento de sua segurança estrutural. Assim, este trabalho teve como objetivo estudar a confiabilidade estrutural de pontes protendidas de madeira, com enfoque na resistência da estrutura à flexão transversal. As pontes foram dimensionadas através de três métodos, o de Ritter, Eurocode e OTB, sendo os dois primeiros conhecidos na literatura e o último baseado na solicitação encontrada através de um software de análise de placas ortotrópicas, OTB. Para obter índices de confiabilidade foi utilizado o método FORM, e o método de Monte Carlo para simular a utilização das fórmulas de obtenção de resistências características sugeridas na NBR 7190:1997 e DIN 68364. Os resultados indicaram que o sistema estrutural é confiável dentro do modo de falha estudado. / Stress laminated timber bridges were first built in Canadá in the seventies. Since then, this kind of structure has been increasingly used in many countries. In Brazil, this structure was first designed in the bridge over Monjolinho stream, in São Carlos, São Paulo state. The importance of this system requires a better knowledge about its structural safety. Therefore, the aim of this work is to research the structural reliability of stress laminated timber bridges, with special focus on transversal bending strength. Bridges were designed by three methods: Ritter, Eurocode and OTB. Ritter and Eurocode are well known design methods. OTB is based on an ortotropic timber bridges analysis software. In order to obtain reliability results, FORM method was used. Monte Carlo method was also considered to simulate characteristic values outputs by brazilian code, NBR 7190:1997, and german code, DIN 68364. Results indicate that the structure os reliable for the limit state studied.

The dose limits are set in terms of protection quantities. However, these quantities are not directly measurable. Thus, it is necessary to calculate the dose conversion coefficients (CCs), which relate these protection quantities with measurable physical quantities. In this work, were calculated CCs for equivalent dose and effective dose in terms of air kerma in a male voxel simulator at standing and sitting postures in air and vacuum environments. The Visual Monte Carlo (VMC) code was used to develop irradiation scenarios with a plane source of monoenergetic photons with energy ranging from 0.01 MeV to 2.0 MeV for irradiation geometries anteroposterior (AP), postero-anterior (PA), left lateral (LLAT), right lateral (rlat) and rotational (ROT). In order to observe possible variations between the CCs calculated, these values were compared for each environment and each posture. Were also compared the values of CCs of the simulator in standing posture with the CCs presented in publications 74 and 116 of the ICRP. The comparison between the values of CCs obtained with the simulator in the standing posture used in this work and the values presented in publications 74 and 116 of the ICRP shows variations in practically all organs, due mainly to differences between the simulators used such as organs mass and anatomy. By comparing the CCs values in air and vacuum environments, variations were observed in the CCs up to 83% at low energies (energies below 0.05 MeV). These variations occur mainly in the AP irradiation geometry. There were also significant differences primarily in AP, PA and ROT irradiation geometries, with variations of up to 100% in the CCs between the two postures in 7 of the 18 organs studied. The comparison between the CCs for effective dose between the sitting and standing posture showed differences in irradiation geometries AP, PA and ROT, with variations of up to 53 %. The divergences found in CCs between sitting and standing posture are attributed mainly to the difference in distribution and geometry of the organs when the posture of the simulator is modified. It is verified thus the importance of obtaining dose conversion coefficients using radiation scenarios in which the simulator are in different environments and postures than those found in the literature / Os limites de dose são estabelecidos em termos de grandezas de proteção. Entretanto, essas grandezas não são diretamente mensuráveis. Dessa forma, é necessário o cálculo de coeficientes de conversão de dose (CCs), que relacionam essas grandezas de proteção com grandezas físicas mensuráveis. Neste trabalho foram calculados CCs de dose equivalente e dose efetiva em termos de kerma no ar em um simulador voxel masculino na postura em pé e sentada em ambientes com ar e vácuo. Foi utilizado o código Visual Monte Carlo (VMC) para elaborar os cenários de irradiação com uma fonte plana de fótons monoenergéticos variando sua energia de 0,01 MeV até 2,0 MeV para as geometrias de irradiação antero-posterior (AP), postero-anterior (PA), lateral esquerdo (LLAT), lateral direito (RLAT) e rotacional (ROT). A fim de observar possíveis variações entre os CCs calculados, os valores para cada ambiente e cada postura foram comparados. Também foram comparados os valores dos CCs do simulador na postura em pé com os CCs apresentados nas publicações nº 74 e nº 116 da ICRP. A comparação entre os valores dos CCs obtidos com o simulador na postura em pé utilizado neste trabalho e os valores apresentados nas publicações nº 74 e nº 116 da ICRP apresenta divergências em praticamente todos os órgãos, em virtude, principalmente, de diferenças entre os simuladores utilizados, como massa e anatomia dos órgãos. Quando comparados os valores dos CCs em ambientes com vácuo e ar, foram observadas variações de até 83 % em baixas energias, energias abaixo de 0,05 MeV. Estas variações ocorrem principalmente na geometria de irradiação AP. Observaram-se, também, diferenças significativas principalmente nas geometrias de irradiação AP, PA e ROT, com variações de até 100 %, nos CCs entre as duas posturas em 7 órgãos dos 18 estudados. A comparação entre os CCs de dose efetiva entre a postura em pé e sentada apresentaram diferenças significativas nas geometrias de irradiação AP, PA e ROT, com variações de até 53 %. As divergências encontradas nos CCs entre a postura em pé e sentada são atribuídas principalmente à diferença na distribuição e geometria dos órgãos quando a postura do simulador é modificada. Verifica-se, dessa forma, a importância da obtenção de coeficientes de conversão de doses utilizando cenários de radiação nos quais os simuladores estão em ambientes e posturas diferentes daquelas encontradas na literatura.