At Saab Aerostructures they are manufacturing a lot of parts for Airbus and Boeing. When these components are investigated using finite element analysis four-node Kirchhoff shell elements and a very fine mesh is often used. In order to make the pre-processing easier Saab would like to offset the shell mid surface from the nodal plane (the modeling surface) rather than to extract mid surfaces for the entire component. This would also make it easier to model a component which needs a thickness change later on, this since the original modeling surface could be used but with an offset of the elements in order to represent the new geometry. When offset is used in Nastran multi point constraints are created between the nodes and the shell mid surface points. All loads, which are applied in the nodal plane, are then transformed to the mid surface where the stiffness matrices, displacements and stresses are calculated. In order to be able to use this method more knowledge about its effects are needed, which is the reason for this thesis work. The offset is studied for two simpler cases, thickness variation and a 90°corner, as well as fora more complicated component called a C-bar. This is a hinge connecting the flaps to the wings of an airplane. The simpler cases are modeled using both mid surface and offset models subject to either a transverse load, an in-plane load or a bending moment. These are compared to a solid model in order to determine which is the most accurate. When mid surface modeling is used fort he thickness variation the surfaces are connected using rigid links. The conclusion made from these simulations is that using offset may give different results if the load is an in-plane load. This kind of load leads to the creation of a bending moment, which is linearly dependent on the amount of offset. The severity of this depends on the overall geometry and how this load is applied.

No description available.

No description available.

Thin and slender structures are widely occurring both in nature and in human creations. Clever geometries of thin structures can produce strong constructions while using a minimal amount of material. Computer modeling and analysis of thin and slender structures has its own set of problems stemming from assumptions made when deriving the equations modeling their behavior from the theory of continuum mechanics. In this thesis we consider two kinds of thin elastic structures; threads and plates. Real-time simulation of threads are of interest in various types of virtual simulations such as surgery simulation for instance. In the first paper of this thesis we develop a thread model for use in interactive applications. By viewing the thread as a continuum rather than a truly one dimensional object existing in three dimensional space we derive a thread model that naturally handles both bending, torsion and inertial effects. We apply a corotational framework to simulate large deformation in real-time. On the fly adaptive resolution is used to minimize corotational artifacts. Plates are flat elastic structures only allowing deflection in the normal direction. In the second paper in this thesis we propose a family of finite elements for approximating solutions to the Kirchhoff-Love plate equation using a continuous piecewise linear deflection field. We reconstruct a discontinuous piecewise quadratic deflection field which is applied in a discontinuous Galerkin method. Given a criterion on the reconstruction operator we prove a priori estimates in energy and L2 norms. Numerical results for the method using three possible reconstructions are presented.

Orientador: Rubén Augusto Romero Lázaro / Resumo: A análise de sistemas altamente complexos quando e analizado o problema de planejamento de expansão de redes de transmissão de longo prazo, é o foco principal deste trabalho. Os modelos e metodos propostos são aplicados ao problema de planejamento estático tradicional, que é um problema de otimização matemática classificado como NP-completo, não-linear inteiro misto. O qual envolve no investimento, variáveis operacionais contínuas e variáveis inteiras. O comportamento normal de cada sistema pode conter informação essencial para a criação de novos métodos, como os planos de corte baseados em cortes de diferença de ângulos para problemas de grande escala, o que é a base é o ponto de partida deste trabalho, derivando em desigualdades válidas é ciclos críticos. Os cortes angulares básicos reduzem o espaço de busca do problema e o tempo total de cálculo deste problema, enquanto ao método de inequações válidas que pode ser usado para fornecer limites inferiores sólidos no investimento ótimo do planejamento de transmissão, já que a diferença entre o modelo DC (modelo exato) e o modelo de transporte (modelo mais relaxado) são as restrições angulares. Os ciclos críticos têm sido desenvolvidos para melhoraralguns dos modelos tradicionais do problemas de planejamento da expansão da rede de transmissão de longo prazo. A razão por trás disso é a ausência da segunda lei de Kirchhoff, que completa a representação do sistema, mas aumenta a complexidade. Para resolver os problemas resultantes... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: The analysis of highly complex systems when solving the long-term transmission network expansion planning problem is the main focus of this work. The proposed improved models and methodology are applied to the traditionalstatic planning problem, which is a mathematical optimization problem classified as NP-complete and mixed-integer nonlinear problem. It involves continuousoperating variables and integer investment variables. The normal behavior of each system can be shown essential information to the creation of new methods, as the cutting-planes based in bus-angle difference cuts for large-scale problems which were the starting point of this work, deriving in valid inequalities and critic cycles. The angular cuts aim to reduce the search space of the problem and the total computation time of this NP-hard problem as for the valid inequalities methodthat can be used to provide strong lower bounds on the optimal investment of the transmissionplanning, since the difference between the DC model (exact model) and the transport model (more relaxed model) are the angular constraints. Critic cycles has been develop in order to improve some of the traditional long-term transmission network expansion planning problem models. The reason behind it is the absence of second Kirchhoff’s law which completes the representationof the system, but increase the complexity. In order to solve the resulting problems, this work uses the modeling language AMPL with the solver CPLEX. In test systems w... (Complete abstract click electronic access below) / Doutor

Ce travail présente et valide une théorie nonlocale nouvelle et généralisée, de la propagation acoustique dans les milieux poreux à structure rigide, saturés par un fluide viscothermique. Cette théorie linéaire permet de dépasser les limites de la théorie classique basée sur la théorie de l'homogénéisation. Elle prend en compte non seulement les phénomènes de dispersion temporelle, mais aussi ceux de dispersion spatiale. Dans le cadre de la nouvelle approche, une nouvelle procédure d'homogénéisation est proposée, qui permet de trouver les propriétés acoustiques à l'échelle macroscopique, en résolvant deux problèmes d'action-réponse indépendants, posés à l'échelle microscopique de Navier-Stokes-Fourier. Contrairement à la méthode classique d'homogénéisation, aucune contrainte de séparation d'échelle n'est introduite. En l'absence de structure solide, la procédure redonne l'équation de dispersion de Kirchhoff-Langevin, qui décrit la propagation des ondes longitudinales dans les fluides viscothermiques. La nouvelle théorie et procédure d'homogénéisation nonlocale sont validées dans trois cas, portant sur des microgéométries significativement différentes. Dans le cas simple d'un tube circulaire rempli par un fluide viscothermique, on montre que les nombres d'ondes et les impédances prédits par la théorie nonlocale, coïncident avec ceux de la solution exacte de Kirchhoff, connue depuis longtemps. Au contraire, les résultats issus de la théorie locale (celle de Zwikker et Kosten, découlant de la théorie classique d'homogénéisation) ne donnent que le mode le plus attenué, et encore, seulement avec le petit désaccord existant entre la solution simplifiée de Zwikker et Kosten et celle exacte de Kirchhoff. Dans le cas où le milieu poreux est constitué d'un réseau carré de cylindres rigides parallèles, plongés dans le fluide, la propagation étant regardée dans une direction transverse, la vitesse de phase du mode le plus atténué peut être calculée en fonction de la fréquence en suivant les approches locale et nonlocale, résolues au moyen de simulations numériques par la méthode des Eléments Finis. Elle peut être calculée d'autre part par une méthode complètement différente et quasi-exacte, de diffusion multiple prenant en compte les effets viscothermiques. Ce dernier résultat quasi-exact montre un accord remarquable avec celui obtenu par la théorie nonlocale, sans restriction de longueur d'onde. Avec celui de la théorie locale, l'accord ne se produit que tant que la longueur d'onde reste assez grande.

Submitted by Johnny Rodrigues (johnnyrodrigues@ufcg.edu.br) on 2018-07-24T14:12:47Z No. of bitstreams: 1 NATAN DE ASSIS LIMA - DISSERTAÇÃO PPGMAT 2010..pdf: 614405 bytes, checksum: d19b00bf4d0fb78e21179e363cfc96f8 (MD5) / Made available in DSpace on 2018-07-24T14:12:47Z (GMT). No. of bitstreams: 1 NATAN DE ASSIS LIMA - DISSERTAÇÃO PPGMAT 2010..pdf: 614405 bytes, checksum: d19b00bf4d0fb78e21179e363cfc96f8 (MD5) Previous issue date: 2010-03 / CNPq / Neste trabalho usaremos algumas técnicas da Análise Funcional Não-Linear para estudar a existência de solução para os chamados Problemas Elípticos Não-Locais, entre os quais destacamos aqueles que incluem o operador de Kirchhoff [...]. * Para visualizar o resumo recomendamos do download do arquivo uma vez que o mesmo utiliza formulas ou equações matemáticas que não puderam ser transcritas neste espaço. / In this work we will use same techniques of Nonlinear Analysis Functional to study the existence of solutions for the some Nonlocal Elliptic Problems, among then those which include Kirchhoff operator [...]. * To preview the summary we recommend downloading the file since it uses mathematical formulas or equations that could not be transcribed in this space.

À l’heure actuelle, les modèles photochimiques utilisés dans la correction des effets de proximitéoptique (OPC) en photolithographie sont devenus complexes et moins physiques afin de permettre decapturer rapidement le maximum d’effets optiques et chimiques. La question de la stabilité de tels modèlespurement empiriques est devenue d’actualité. Dans ce mémoire, nous avons étudié la stabilité desmodèles photochimiques actuels en examinant les différentes causes d’instabilité vis-à-vis des paramètresdu procédé. Dans la suite, nous avons développé une méthode perturbative permettant d’évaluer le critèrede la stabilité. L’obtention de modèles simples et stables nous conduit à séparer les effets optiques desautres effets chimiques. De plus, les approximations utilisées dans la modélisation des systèmes optiquesopérant à grande ouverture numérique entraînent des erreurs résiduelles pouvant dégrader la précisionet la stabilité des modèles OPC. Ainsi, nous nous sommes intéressés à étudier les limites de validitéde l’approximation de Kirchhoff, méthode qui, jusqu’à présent, est la plus utilisée dans la modélisationdu champ proche d’un masque. D’autres méthodes semi-rigoureuses, permettant de modéliser les effetstopographiques, ont été également évaluées. Ces méthodes approchées permettent de gagner en précisionmais dégradent le temps de calcul. Nous avons ainsi proposé différentes façons de corriger les effetstopographiques du masque, tout en gardant l’approximation de Kirchhoff dans la modélisation de la partieoptique. Parmi les méthodes proposées, nous exploitons celle permettant de réduire les erreurs liéesaux effets topographiques du masque par l’intermédiaire d’un second modèle empirique. Nous montronsque pour garantir une précision adéquate, il est nécessaire d’augmenter la complexité du modèle en rajoutantdes termes additionnels. Enfin, pour garantir la stabilité numérique du modèle empirique, nousintroduirons une nouvelle méthode approchée hybride rapide et précise, la méthode des multi-niveaux,permettant d’inclure les effets topographiques par décomposition multi-niveaux du masque fin et discuteronsses avantages et ses limites. / At present, common resist models utilized in photolithography to correct for optical proximity effects(OPC) became complex and less physical in order to capture the maximum of optical and chemical effectsin shorter times. The question on the stability of such models, purely empirical, become topical. In thisthesis, we study the stability of existing OPC models by examining the origins of model instability towardsprocess parameters. Thus, we have developed a perturbative method in order to evaluate the stabilitycriterion. However, achieving stable and simple models needs a separation between optical and otherchemical effects. Besides, multiple approximations, widely utilized in the modeling of optical systemsoperating at high numerical aperture, lead to residual errors which can degrade OPC model accuracyand stability. Thus, we were interested to study the limits of validity of the Kirchhoff approximation,a method which, so far, is the most commonly used in mask near-field modeling. Other semi-rigorousmethods for mask topography effect modeling were also evaluated. These approximate methods canimprove the accuracy but degrades the run time. We then suggested different techniques to correct formask topography effects, while keeping the Kirchhoff approximation in the modeling of the optical part.Among them, we showed that errors due to mask topography effects can be partially captured by asecond empirical model. However, in order to ensure a good accuracy, it is necessary to increase themodel complexity by using more additional empirical terms. Finally, in order to achieve a numericalstability of the empirical model, we introduced a new hybrid fast and accurate method, the multi-levelmethod, which allows us to correct for mask topography effects through a multi-level decomposition ofthe thin mask and discussed its advantages and drawbacks.

Les structures en forme de tige sont omniprésentes dans le monde aujourd'hui. Désormais, prédire avec précision leur comportement pour l'ingénierie et les environnements virtuels est indispensable pour de nombreuses industries, notamment l'infographie, l'animation par ordinateur et la conception informatique. Dans ce mémoire, nous explorons un nouveau modèle de calcul pour les tiges élastiques qui exploite les données de simulation pour reproduire les effets de coque et de solide présents dans les tiges qui brisent les hypothèses de la théorie classique de la tige de Kirchhoff, présentant ainsi une voie d'amélioration possible pour de nombreux états de l'art techniques. Notre approche consiste à prendre un ensemble de données de simulations à partir de solides volumétriques ou de coques pour former un nouveau modèle d'énergie définie positive polynomiale d'ordre élevé pour une tige élastique. Cette nouvelle énergie élargit la gamme des comportements des matériaux qui peuvent être modélisés pour la tige, permettant ainsi de capturer une plus large gamme de phénomènes. Afin de proposer et tester ce modèle, nous concevons un pipeline expérimental pour tester les limites de la théorie linéaire des tiges et étudier les géométries d'interface entre les cas coque à tige et volume à coque pour observer les effets d'un modèle de matériau non linéaire et une section transversale non elliptique dans la déformation de la tige. Nous étudions également la relation entre la courbure de la tige et la déformation de la section transversale et la courbure pour introduire une modification sur le terme de flexion de l'énergie. Cela nous permet de reproduire à la fois le comportement de flexion asymétrique présent dans les poutres volumétriques minces et les poutres à coque avec des sections transversales non convexes. Des suggestions pour de nouvelles améliorations des modèles et des techniques expérimentales sont également données. / Rod-like structures are ubiquitous in the world today. Henceforth accurately predicting their behavior for engineering and virtual environments are indispensable for many industries including computer graphics, computer animation, and computational design. In this thesis we explore a new computational model for elastic rods that leverages simulation data to reproduce shell and solid-like effects present in rods that break the assumptions of the classical Kirchhoff rod theory, thus presenting a possible improvement avenue to many states-of-the-art techniques. Our approach consists of taking a data set of simulations from both volumetric solids or shells to train a novel high-order polynomial positive-definite energy model for an elastic rod. This new energy increases the range of material behaviors that can be modeled for the rod, thus allowing for a larger range of phenomena to be captured. In order to propose and test this model, we design an experimental pipeline to test the limits of the linear theory of rods and investigate the interface geometries between the Shell-Rod and Volume-Shell cases to observe the effects of a nonlinear material model and a non-elliptical cross-section in the rod deformation. We also investigate the relation between rod curvature and deformation of the cross-section and curvature to introduce a modification on the bending term of the energy. This allows us to reproduce both the asymmetric bending behavior present in thin volumetric solid and shell beams with non-convex cross-sections. Suggestions for further improvements in models and experimental techniques are also given.

The prediction of fracture in thin-walled structures is decisive for a wide range of applications. Modeling methods such as the phase-field method usually consider cracks to be constant over the thickness which, especially in load cases involving bending, is an imperfect approximation. In this contribution, fracture phenomena along the thickness direction of structural elements (plates or shells) are addressed with a phase-field modeling approach. For this purpose, a new, so called “mixed-dimensional” model is introduced, which combines structural elements representing the displacement field in the two-dimensional shell midsurface with continuum elements describing a crack phase-field in the three-dimensional solid space. The proposed model uses two separate finite element discretizations, where the transfer of variables between the coupled twoand three-dimensional fields is performed at the integration points which in turn need to have corresponding geometric locations. The governing equations of the proposed mixed-dimensional model are deduced in a consistent manner from a total energy functional with them also being compared to existing standard models. The resulting model has the advantage of a reduced computational effort due to the structural elements while still being able to accurately model arbitrary through-thickness crack evolutions as well as partly along the thickness broken shells due to the continuum elements. Amongst others, the higher accuracy aswell as the numerical efficiency of the proposed model are tested and validated by comparing simulation results of the new model to those obtained by standard models using numerous representative examples.