A study is made of the problem of metamaterial homogenization, which is the attempt to represent an artificially fabricated inhomogeneous periodic structure as a homogeneous medium with an electromagnetic response described by a number of constitutive parameters (permittivity, permability, etc.) In particular, the importance of spatial dispersion in metamaterials and the need to characterize metamaterials with wavevector dependent constitutive parameters is explained an examined. A brief survey of important previous attempts at metamaterial homogenization is presented. This is followed by a discussion of spatial dispersion in metamaterial crystals. The importance of spatial dispersion in metamaterials is justified and some manifestations of spatial dispersion described. In particular the little known phenomenon of bianisotropy in centrosymmetric crystals due to spatial dispersion is explained. Also, the effects of spatial dispersion on physical quantities such as energy flux and dissipation are identified. We then describe a new method for solving for the free eigenmodes of a metamaterial crystal with a complex wavevector eigenvalue simulation. Next, two different theoretical attempts by the author at metamaterial homogenization are described, both accompanied by tests of the calculated constitutive parameters and critical examination of the strengths and weaknesses of each approach. Finally, strong evidence of the presence and importance of spatial dispersion in metamaterials is presented. / text

This work presents an effective methodology for measuring the depth-resolved 3D full-field deformation of semitransparent, light scattering soft tissues such as vertebrate eye cornea. This was obtained by performing digital volume correlation on optical coherence tomography volume reconstructions of silicone rubber phantoms and porcine cornea samples. Both the strip tensile tests and the posterior inflation tests have been studied. Prior to these tests, noise effect and strain induced speckle decorrelation were first studied using experimental and simulation methods. The interpolation bias in the strain results has also been analyzed. Two effective approaches have been introduced to reduce the interpolation bias. To extract material constitutive parameters from the 3D full-field deformation measurements, the virtual fields method has been extended into 3D. Both manually defined virtual fields and the optimized piecewise virtual fields have been developed and compared with each other. Efforts have also been made in developing a method to correct the refraction induced distortions in the optical coherence tomography reconstructions. Tilt tests of different silicone rubber phantoms have been implemented to evaluate the performance of the refraction correction method in correcting the distorted reconstructions.

Análises numéricas confiáveis do comportamento mecânico de materiais como borrachas, dependem muito de uma calibração precisa do modelo constitutivo hiperelástico utilizado. Estes modelos são calibrados ajustando as curvas teóricas aos dados experimentais, obtidos por meio de ensaios usuais. Em muitos casos as amostras de matéria prima desses elastômeros já não se encontram disponíveis ou é impossível fabricar os corpos de prova requeridos. O objetivo deste trabalho é verificar a possibilidade de encontrar constantes constitutivas testando o próprio componente, ao invés dos usuais ensaios de tração, compressão e cisalhamento. A abordagem proposta consiste em criar uma rotina de programação associada a uma função custo onde, a partir de uma estimativa inicial de constantes constitutivas, sejam realizados processos iterativos de otimização buscando aproximar as curvas de força × deslocamento teórica e experimental. Um componente automotivo será utilizado nos estudos e dois modelos constitutivos hiperelásticos serão testados. As equações de tensões nominais dos modelos hiperelásticos serão utilizadas para predizer o comportamento teórico dos ensaios usuais, de forma a verificar a qualidade das constantes obtidas. Conclui-se que é possível utilizar o ensaio da própria peça para caracterizar o material hiperelástico, com resultados comparáveis aos que seriam obtidos com os ensaios típicos para esta aplicação. / The reliable numerical analysis of the mechanical behavior of rubber-like materials depends strongly on accurately calibrated hyperelastic constitutive models. Such models are calibrated by fitting theoretical curves against experimental data obtained in well known tests. In many cases samples of the original elastomer are no longer available or it is impossible to manufacture the specimens required by the standard tests. The aim of this work is verify the possibility of finding the constitutive constants by testing the actual component instead of the usual tensile, compression and shear tests. The proposed approach consists in creating a programming routine with a cost function that, starting from an initial estimate of the constitutive constants, iterate through an optimization algorithm in order to fit the theoretical force × displacement curves to the experimental ones. An automotive component will be used during the studies and two hyperelastic constitutive models will be tested. The nominal stress equations for the hyperelastic models are used to predict the standard tests behavior, to assess the quality of the constants obtained. The results shown that is possible to characterize an hyperelastic material by testing the actual component, with results comparable to those which would be obtained with standard tests.

Análises numéricas confiáveis do comportamento mecânico de materiais como borrachas, dependem muito de uma calibração precisa do modelo constitutivo hiperelástico utilizado. Estes modelos são calibrados ajustando as curvas teóricas aos dados experimentais, obtidos por meio de ensaios usuais. Em muitos casos as amostras de matéria prima desses elastômeros já não se encontram disponíveis ou é impossível fabricar os corpos de prova requeridos. O objetivo deste trabalho é verificar a possibilidade de encontrar constantes constitutivas testando o próprio componente, ao invés dos usuais ensaios de tração, compressão e cisalhamento. A abordagem proposta consiste em criar uma rotina de programação associada a uma função custo onde, a partir de uma estimativa inicial de constantes constitutivas, sejam realizados processos iterativos de otimização buscando aproximar as curvas de força × deslocamento teórica e experimental. Um componente automotivo será utilizado nos estudos e dois modelos constitutivos hiperelásticos serão testados. As equações de tensões nominais dos modelos hiperelásticos serão utilizadas para predizer o comportamento teórico dos ensaios usuais, de forma a verificar a qualidade das constantes obtidas. Conclui-se que é possível utilizar o ensaio da própria peça para caracterizar o material hiperelástico, com resultados comparáveis aos que seriam obtidos com os ensaios típicos para esta aplicação. / The reliable numerical analysis of the mechanical behavior of rubber-like materials depends strongly on accurately calibrated hyperelastic constitutive models. Such models are calibrated by fitting theoretical curves against experimental data obtained in well known tests. In many cases samples of the original elastomer are no longer available or it is impossible to manufacture the specimens required by the standard tests. The aim of this work is verify the possibility of finding the constitutive constants by testing the actual component instead of the usual tensile, compression and shear tests. The proposed approach consists in creating a programming routine with a cost function that, starting from an initial estimate of the constitutive constants, iterate through an optimization algorithm in order to fit the theoretical force × displacement curves to the experimental ones. An automotive component will be used during the studies and two hyperelastic constitutive models will be tested. The nominal stress equations for the hyperelastic models are used to predict the standard tests behavior, to assess the quality of the constants obtained. The results shown that is possible to characterize an hyperelastic material by testing the actual component, with results comparable to those which would be obtained with standard tests.

Análises numéricas confiáveis do comportamento mecânico de materiais como borrachas, dependem muito de uma calibração precisa do modelo constitutivo hiperelástico utilizado. Estes modelos são calibrados ajustando as curvas teóricas aos dados experimentais, obtidos por meio de ensaios usuais. Em muitos casos as amostras de matéria prima desses elastômeros já não se encontram disponíveis ou é impossível fabricar os corpos de prova requeridos. O objetivo deste trabalho é verificar a possibilidade de encontrar constantes constitutivas testando o próprio componente, ao invés dos usuais ensaios de tração, compressão e cisalhamento. A abordagem proposta consiste em criar uma rotina de programação associada a uma função custo onde, a partir de uma estimativa inicial de constantes constitutivas, sejam realizados processos iterativos de otimização buscando aproximar as curvas de força × deslocamento teórica e experimental. Um componente automotivo será utilizado nos estudos e dois modelos constitutivos hiperelásticos serão testados. As equações de tensões nominais dos modelos hiperelásticos serão utilizadas para predizer o comportamento teórico dos ensaios usuais, de forma a verificar a qualidade das constantes obtidas. Conclui-se que é possível utilizar o ensaio da própria peça para caracterizar o material hiperelástico, com resultados comparáveis aos que seriam obtidos com os ensaios típicos para esta aplicação. / The reliable numerical analysis of the mechanical behavior of rubber-like materials depends strongly on accurately calibrated hyperelastic constitutive models. Such models are calibrated by fitting theoretical curves against experimental data obtained in well known tests. In many cases samples of the original elastomer are no longer available or it is impossible to manufacture the specimens required by the standard tests. The aim of this work is verify the possibility of finding the constitutive constants by testing the actual component instead of the usual tensile, compression and shear tests. The proposed approach consists in creating a programming routine with a cost function that, starting from an initial estimate of the constitutive constants, iterate through an optimization algorithm in order to fit the theoretical force × displacement curves to the experimental ones. An automotive component will be used during the studies and two hyperelastic constitutive models will be tested. The nominal stress equations for the hyperelastic models are used to predict the standard tests behavior, to assess the quality of the constants obtained. The results shown that is possible to characterize an hyperelastic material by testing the actual component, with results comparable to those which would be obtained with standard tests.

Cette thèse traite de l'analyse limite des structures de châssis en acier, qui s'utilise souvent comme la structure principale de support des bâtiments. La structure du cadre en acier est caractérisée par une réponse très ductile et un grand potentiel pour dissiper l'énergie, ce qui est crucial pour la résistance par rapport aux tremblements de terre. La ductilité dans la réponse de la structure est la cause du comportement du matériau lui-même et du comportement des connexions entre les éléments de la structure. Les connexions entre les poutres et les poteaux peuvent influencer de manière significative la réponse de la structure du cadre en acier, parfois jusqu'à 30%. L'idée est de intégrer le comportement des connexions par les éléments de poutres qui seront situés dans les coins du cadre et la modélisation du reste serra fait avec des éléments de poutres non-linéaires qui décrirons le comportement des poutres en acier. Cette recherche est composée de deux parties. La première partie est consacrée au comportement des connexions structurelles, la deuxième partie présente le développement de l'élément fini du faisceau non linéaire capable de représenter le comportement ductile d'un élément de la structure en acier. Dans la première partie de la thèse, nous définissons la procédure d'identification des paramètres constitutifs pour le modèle couplé de plasticité-dégâts avec dix-huit inconnus. Ce modèle constitutif est très robuste et capable de représenter une large gamme de problèmes. La procédure définie a été utilisée dans la préparation de tests expérimentaux pour trois types de connexions en acier structuré. Les tests expérimentaux ont été effectués pour deux cas de charge. Pour la première, la charge a été appliquée dans un sens avec les cycles de chargement et de déchargement. À partir des mesures expérimentales, nous avons conclu que le modèle de plasticité peut bien représentée le comportement de la connexion structurale. Paramètres constitutifs ont été déterminés à partir des résultats de l'expérimentation, on a utilisé une poutre géométriquement exacte avec la loi bilinéaires renforcement du matériel et la loi linéaire pour le ramollissement. Également, on a effectué des essais expérimentaux de deux types de raccords en acier en cas de chargement cyclique. Les données mesurées montrent que le modèle de la plasticité n'est pas assez bon pour décrire le comportement de connexion pour ce type de charge. A savoir, en raison de changements du sens de l'application du chargement, les connexions montrent moins de rigidité, qui peut être décrite avec un modèle constitutif de dommages. Pour cette raison, nous avons développé un nouveau modèle plasticité-dommages qui est capable d'inclure le phénomène mentionné ci-dessus. A la fin de cette section est faite l'identification des paramètres constitutifs. La deuxième partie de la thèse de doctorat est composé de formulations théoriques et la mise en œuvre numérique des faisceaux géométriquement exacte. La réponse de durcissement de la poutre comprend l'interaction entre les forces de la section résultant du stress (N, T et M), et la réponse de ramollissement est définit par la loi non linéaire. Ce type d'élément fini de poutre est capable de décrire le comportement ductile des structures en acier et inclure les effets du second ordre, qui sont très importantes pour l'analyse ultime des structures de cadre en acier. L'élément fini développé de poutre géométriquement exacte et les lois définies de liaison de comportement dans la construction en acier, offrant la possibilité d'une analyse de haute qualité des structures en acier. En utilisant les modèles de poutre proposé et la méthodologie de modélisation des structures de châssis en acier, il est possible de déterminer une distribution réaliste des forces de section transversale , y compris la redistribution due à la formation de rotules plastiques. / This thesis deals with the ultimate load limit analysis of steel frame structures. The steel frame structure has a very ductile response and a large potential to dissipate energy, which is crucial in the case of earthquakes. The ductility in the response of the structure comes from the behavior of the material itself and the behavior of the semi-rigid structural connections. The semi-rigid connections between beams and columns can significantly influence the response of the structure, sometimes up to 30%. In this thesis, we propose a methodology for modeling steel frame structures with included connection behavior. The idea is to model the behavior of the structural connections by the beam elements positioned in the corners of the steel frame structure. Other members of the steel frame structure, steel beams, and columns, will be modeled with nonlinear beam elements. This research consists of two parts. The first part deals with the behavior of the structural steel connections. In the second part, we present the development of the nonlinear beam element capable of representing the ductile behavior of steel structural elements, beams and columns. In the first part of the thesis, we define constitutive parameters identification procedure for the coupled plasticity-damage model with eighteen unknowns. This constitutive model is very robust and capable of representing a wide range of problems. The identification procedure was used in the preparation of experimental tests for three different types of structural steel connections. The experimental tests have been performed for two load cases. In the first, the load was applied in one direction with both the loading and unloading cycles. From the experimental measurements, we have concluded that the response of the experimental structure can be represented by the plasticity model only because no significant change in the elastic response throughout the loading program was observed. Therefore, we have chosen an elastoplastic geometrically exact beam to describe connection behavior. The hardening response of the beam is governed by bilinear law, and the softening response is governed by nonlinear exponential law. The identification of the parameters has been successfully done with fifteen unknown parameters identified. The two types of the experimental structures were also exposed to the cyclic loading. Measured experimental data shows complex connection behavior that cannot be described by the plasticity model alone. Namely, after changing load direction stiffness of the connection decreases. This suggests that the damage model should be incorporated in the constitutive law for the connections behavior as well. Therefore, we propose a new coupled plasticity-damage model capable of representing the loss in the stiffness of the connection with the changing of the load direction. At the end of this part, we also give the constitutive parameters identification for the proposed model. The second part of the thesis deals with the theoretical formulation and numerical implementation of the elastoplastic geometrically exact beam. The hardening response of the beam includes interaction between stress resultant section forces (N, T and M), and the softening response of the beam, which is governed by the nonlinear law. This type of the beam element is capable of representing the ductile behavior of a steel frame structure, and it takes into account second order theory effects. Performed numerical simulations show that the proposed geometrically nonlinear beam element is very robust and is able to provide a more precise limit load analysis of steel frame structures. By using proposed methodology for modeling steel structures, we are able to obtain the real distribution of section forces, including their redistribution caused by forming of the hinges and the connections behavior.

Ce travail présente une méthodologie efficace pour mesurer le champ de déformation total en 3-D des tissus mous semi-transparents et diffusant la lumière, tels que la cornée de l'œil vertébré. Cela a été obtenu en utilisant la tomographie par cohérence optique couplée avec la corrélation volumique numérique sur des fantômes en caoutchouc de silicone et de la cornée porcine. Des tests de traction et des tests d'inflation postérieure ont été étudiés. Avant ces tests, l'effet du bruit et la décorrélation des speckles due à la déformation sont d'abord étudiés en utilisant des méthodes expérimentales et numériques. Le biais d'interpolation dans les résultats de déformation a été analysé. Deux approches efficaces ont été prises pour réduire le biais d'interpolation. Pour extraire les paramètres constitutifs des mesures de déformation 3-D, la méthode des champs virtuels a été étendue en 3-D. Les champs virtuels définis manuellement et les champs virtuels optimisés par morceaux ont été développés et comparés entre eux. Des efforts ont également été déployés pour corriger les distorsions induites par réfraction dans les reconstructions de la tomographie par cohérence optique. Des tests d'inclinaison des différents fantômes de silicone ont été introduits afin d’évaluer la performance de la méthode pour corriger les reconstructions distordues / This work presents an effective methodology for measuring the depth-resolved 3-D full-field deformation of semitransparent, light scattering soft tissues such as vertebrate eye cornea. This was obtained by performing digital volume correlation on optical coherence tomography volume reconstructions of silicone rubber phantoms and porcine cornea samples. Both the strip tensile tests and the posterior inflation tests have been studied. Prior to these tests, noise effect and strain induced speckle decorrelation were first studied using experimental and simulation methods. The interpolation bias in the strain results has also been analyzed. Two effective approaches have been introduced to reduce the interpolation bias. To extract material constitutive parameters from the 3-D full-field deformation measurements, the virtual fields method has been extended into 3-D. Both manually defined virtual fields and the optimized piecewise virtual fields have been developed and compared with each other. Efforts have also been made in developing a method to correct the refraction induced distortions in the optical coherence tomography reconstructions. Tilt tests of different silicone rubber phantoms have been implemented to evaluate the performance of the refraction correction method in correcting the distorted reconstructions