Em geral, as formulações disponíveis na literatura para a análise de cascas laminadas reforçadas com fibras substituem o meio original heterogêneo por um homogêneo equivalente, que dificulta a identificação das tensões fibra-matriz, ou requerem que a malha de elementos finitos seja disposta de modo que os nós dos elementos finitos de fibra coincidam com os nós dos elementos finitos de casca, que é uma exigência bastante restritiva e que aumenta o número de graus de liberdade do sistema de equações resultante. Neste sentido, o objetivo geral desta tese consiste em desenvolver uma formulação para a inclusão de fibras longas e curtas aleatórias nas diversas lâminas de cascas laminadas anisotrópicas com não linearidade geométrica utilizando o método dos elementos finitos sem aumentar o número de graus de liberdade do sistema de equações resultante e sem a necessidade de coincidência de nós na discretização das fibras e da matriz. Nesta formulação, o elemento finito triangular de casca laminada utilizado para discretizar a matriz possui dez nós e sete graus de liberdade por nó, sendo três translações, três componentes do vetor generalizado e a taxa de variação linear da deformação ao longo da espessura. As fibras curvas, curtas aleatórias ou longas, são introduzidas, em qualquer camada do laminado, por meio de relações cinemáticas que garantem sua aderência à matriz sem a introdução de novos graus de liberdade no sistema de equações resultante. Para discretizá-las são utilizados elementos finitos unidimensionais de ordem qualquer com três graus de liberdade por nó e que consideram consistentemente a não linearidade geométrica. Todas as grandezas envolvidas são escritas em relação à configuração inicial do corpo, caracterizando a descrição Lagrangeana total ou material do movimento. Para modelar o comportamento do material adota-se a Lei Constitutiva de Saint-Venant-Kirchhoff que relaciona de forma linear o tensor de tensões de Piolla-Kirchhoff de segunda espécie e o tensor de deformações de Green-Lagrange. O equilíbrio é encontrado a partir do Princípio da Mínima Energia Potencial Total e o sistema não linear de equações resultante é resolvido utilizando-se o procedimento iterativo de Newton-Raphson. As ações externas podem ser introduzidas ao sistema de forma total ou incremental e a contribuição das fibras para a energia do sistema é adicionada na matriz global do problema. Os exemplos numéricos testados validam e demonstram as potencialidades da formulação proposta. / In general, the Finite Element (FE) formulations available in the literature for the analysis of fibre reinforced laminated shells replace the original heterogeneous medium by an equivalent homogeneous one, which makes difficult the identification of fiber-matrix stress distribution, or require that the finite element mesh is arranged in a way that the fibre finite element nodes coincide with the shell finite element ones, which is a very restrictive requirement and increases the number of degrees of freedom of the resulting system of equations. In this sense, the objective of this thesis is to develop a formulation for the inclusion of long and random short fibres in any layer of FE laminated anisotropic shells developing large displacement and rotations without increasing the number of degrees of freedom and the necessity of matching nodes in the discretization of the fibre and the matrix. In this formulation, the triangular laminated shell finite element used to discretize the matrix has ten nodes and seven degrees of freedom per node, that are, three translations, three components of a generalized vector and the linear rate of strain variation along the thickness. The curved fibres, long or random short, are introduced in any layer of the laminate shell by means of kinematic relation to ensure its adherence to the matrix without introducing new degrees of freedom in the resulting system of equations. To discretize them, any order one-dimensional finite elements with three degrees of freedom per node are used. These fibres elements are consistently considered by Geometric nonlinearity. All involved variables are written with respect to the initial configuration of the body, characterizing the Total Lagrangian description. To model the behavior of the material we use the Saint-VenantKirchhoff Constitutive Law that relates linearly the second Piolla-Kirchhoff stress tensor and Green-Lagrange strain tensor. The equilibrium is achieved from the Principle of Minimum Potential Energy and the non-linear system of equations is solved by the Newton-Raphson iterative procedure. External loads may be introduced to the system by one or various steps and the contribution of fibres to the energy of the system is added to the global matrix of the problem. The numerical examples validate and demonstrate the potential of the proposed formulation.

Em geral, as formulações disponíveis na literatura para a análise de cascas laminadas reforçadas com fibras substituem o meio original heterogêneo por um homogêneo equivalente, que dificulta a identificação das tensões fibra-matriz, ou requerem que a malha de elementos finitos seja disposta de modo que os nós dos elementos finitos de fibra coincidam com os nós dos elementos finitos de casca, que é uma exigência bastante restritiva e que aumenta o número de graus de liberdade do sistema de equações resultante. Neste sentido, o objetivo geral desta tese consiste em desenvolver uma formulação para a inclusão de fibras longas e curtas aleatórias nas diversas lâminas de cascas laminadas anisotrópicas com não linearidade geométrica utilizando o método dos elementos finitos sem aumentar o número de graus de liberdade do sistema de equações resultante e sem a necessidade de coincidência de nós na discretização das fibras e da matriz. Nesta formulação, o elemento finito triangular de casca laminada utilizado para discretizar a matriz possui dez nós e sete graus de liberdade por nó, sendo três translações, três componentes do vetor generalizado e a taxa de variação linear da deformação ao longo da espessura. As fibras curvas, curtas aleatórias ou longas, são introduzidas, em qualquer camada do laminado, por meio de relações cinemáticas que garantem sua aderência à matriz sem a introdução de novos graus de liberdade no sistema de equações resultante. Para discretizá-las são utilizados elementos finitos unidimensionais de ordem qualquer com três graus de liberdade por nó e que consideram consistentemente a não linearidade geométrica. Todas as grandezas envolvidas são escritas em relação à configuração inicial do corpo, caracterizando a descrição Lagrangeana total ou material do movimento. Para modelar o comportamento do material adota-se a Lei Constitutiva de Saint-Venant-Kirchhoff que relaciona de forma linear o tensor de tensões de Piolla-Kirchhoff de segunda espécie e o tensor de deformações de Green-Lagrange. O equilíbrio é encontrado a partir do Princípio da Mínima Energia Potencial Total e o sistema não linear de equações resultante é resolvido utilizando-se o procedimento iterativo de Newton-Raphson. As ações externas podem ser introduzidas ao sistema de forma total ou incremental e a contribuição das fibras para a energia do sistema é adicionada na matriz global do problema. Os exemplos numéricos testados validam e demonstram as potencialidades da formulação proposta. / In general, the Finite Element (FE) formulations available in the literature for the analysis of fibre reinforced laminated shells replace the original heterogeneous medium by an equivalent homogeneous one, which makes difficult the identification of fiber-matrix stress distribution, or require that the finite element mesh is arranged in a way that the fibre finite element nodes coincide with the shell finite element ones, which is a very restrictive requirement and increases the number of degrees of freedom of the resulting system of equations. In this sense, the objective of this thesis is to develop a formulation for the inclusion of long and random short fibres in any layer of FE laminated anisotropic shells developing large displacement and rotations without increasing the number of degrees of freedom and the necessity of matching nodes in the discretization of the fibre and the matrix. In this formulation, the triangular laminated shell finite element used to discretize the matrix has ten nodes and seven degrees of freedom per node, that are, three translations, three components of a generalized vector and the linear rate of strain variation along the thickness. The curved fibres, long or random short, are introduced in any layer of the laminate shell by means of kinematic relation to ensure its adherence to the matrix without introducing new degrees of freedom in the resulting system of equations. To discretize them, any order one-dimensional finite elements with three degrees of freedom per node are used. These fibres elements are consistently considered by Geometric nonlinearity. All involved variables are written with respect to the initial configuration of the body, characterizing the Total Lagrangian description. To model the behavior of the material we use the Saint-VenantKirchhoff Constitutive Law that relates linearly the second Piolla-Kirchhoff stress tensor and Green-Lagrange strain tensor. The equilibrium is achieved from the Principle of Minimum Potential Energy and the non-linear system of equations is solved by the Newton-Raphson iterative procedure. External loads may be introduced to the system by one or various steps and the contribution of fibres to the energy of the system is added to the global matrix of the problem. The numerical examples validate and demonstrate the potential of the proposed formulation.

[pt] Este trabalho apresenta aplicações da técnica de otimização topológica para problemas de escoamento com fluidos não Newtonianos, usando o método dos elementos virtuais (VEM) em domínios bidimensionais arbitrários. O objetivo é projetar a trajetória ótima, a partir da minimização da energia dissipativa, de um escoamento governado pelas equações de Navier-Stokes-Brinkman e do modelo não Newtoniano de Carreau-Yasuda. A abordagem de porosidade proposta por (Borrvall e Petersson, 2003) [1] é usada na formulação do problema de otimização topológica. Para resolver este problema numericamente é usado o método VEM, recentemente proposto. A principal característica que diferencia o VEM do método dos elementos finitos (FEM) é que as funções de interpolação no interior dos elementos não precisam ser computadas explicitamente. Isso ocorre porque a integração é feita em funções polinomiais e bases de ordem inferior, permitindo assim uma grande flexibilidade no que diz respeito ao uso de elementos não convexos. Portanto, o cálculo das matrizes e vetores elementares se reduz à avaliação de grandezas geométricas nos contornos desses elementos. Finalmente, são apresentados exemplos numéricos representativos para demonstrar a eficiência do VEM em comparação com o FEM e a aplicabilidade da otimização topológica para esta classe de problemas de escoamento. / [en] This work presents selected applications of topology optimization for non-Newtonian fluid flow problems using the virtual element method (VEM) in arbitrary two-dimensional domains. The objective is to design an optimal layout into a fluid flow domain to minimize dissipative energy governed by the Navier-Stokes-Brinkman and non-Newtonian Carreau-Yasuda model equations. The porosity approach proposed by (Borrvall and Petersson, 2003) [1] is used in the topology optimization formulation. To solve this problem numerically, the recently proposed VEM method is used. The key feature that distinguishes VEM from the standard finite element method (FEM) is that the interpolation functions in the interior of the elements do not need to be computed explicitly. This is because the integration is on lower-order polynomial and basis functions, and there is great flexibility by using a non-convex element. Therefore, the computation of the main element matrices and vectors are reduced to the evaluation of geometric quantities on the boundary of the elements. Finally, several numerical examples are provided to demonstrate the efficiency of the VEM compared to FEM and the applicability of the topology optimization to fluid flow problems.

L'objectif de cette thèse est de développer un nouveau code pour prédire l'équilibre final d'un processus chimique complexe impliquant beaucoup de produits, plusieurs phases et plusieurs processus chimiques. Des méthodes numériques ont été développées au cours des dernières décennies pour prédire les équilibres chimiques finaux en utilisant le principe de minimisation de l'enthalpie libre du système. La plupart des méthodes utilisent la méthode des « multiplicateurs de Lagrange » et résolvent les équations en employant une approximation du problème de Lagrange et en utilisant un algorithme de convergence pas à pas de type Newton-Raphson. Les équations mathématiques correspondantes restent cependant fortement non linéaires, de sorte que la résolution, notamment de systèmes multiphasiques, peut être très aléatoire. Une méthode alternative de recherche du minimum de l’énergie de Gibbs (MCGE) est développée dans ce travail, basée sur une technique de Monte-Carlo associée à une technique de Pivot de Gauss pour sélectionner des vecteurs composition satisfaisant la conservation des atomes. L'enthalpie libre est calculée pour chaque vecteur et le minimum est recherché de manière très simple. Cette méthode ne présente a priori pas de limite d’application (y compris pour las mélanges multiphasiques) et l’équation permettant de calculer l’énergie de Gibbs n’a pas à être discrétisée. Il est en outre montré que la précision des prédictions dépend assez significativement des valeurs thermodynamiques d’entrée telles l'énergie de formation des produits et les paramètres d'interaction moléculaire. La valeur absolue de ces paramètres n'a pas autant d’importance que la précision de leur évolution en fonction des paramètres du process (pression, température, ...). Ainsi, une méthode d'estimation cohérente est requise. Pour cela, la théorie de la « contribution de groupe » est utilisée (ceux de UNIFAC) et a été étendue en dehors du domaine d'interaction moléculaire traditionnel, par exemple pour prédire l'énergie de formation d’enthalpie libre, la chaleur spécifique... Enfin, l'influence du choix de la liste finale des produits est discutée. On montre que la prédictibilité dépend du choix initial de la liste de produits et notamment de son exhaustivité. Une technique basée sur le travail de Brignole et Gani est proposée pour engendrer automatiquement la liste des produits stable possibles. Ces techniques ont été programmées dans un nouveau code : CIRCE. Les travaux de Brignole et de Gani sont mis en œuvre sur la base de la composition atomique des réactifs pour prédire toutes les molécules « réalisables ». La théorie de la « contribution du groupe » est mise en œuvre pour le calcul des propriétés de paramètres thermodynamiques. La méthode MCGE est enfin utilisée pour trouver le minimum absolu de la fonction d'enthalpie libre. Le code semble plus polyvalent que les codes traditionnels (CEA, ASPEN, ...) mais il est plus coûteux en termes de temps de calcul. Il peut aussi être plus prédictif. Des exemples de génie des procédés illustrent l'étendue des applications potentielles en génie chimique. / The objective of this work is to develop a new code to predict the final equilibrium of a complex chemical process with many species/reactions and several phases. Numerical methods were developed in the last decades to predict final chemical equilibria using the principle of minimizing the Gibbs free energy of the system. Most of them use the “Lagrange Multipliers” method and solve the resulting system of equations under the form of an approximate step by step convergence technique. Notwithstanding the potential complexity of the thermodynamic formulation of the “Gibbs problem,” the resulting mathematical formulation is always strongly non-linear so that solving multiphase systems may be very tricky and having the difficult to reach the absolute minimum. An alternative resolution method (MCGE) is developed in this work based on a Monte Carlo technique associated to a Gaussian elimination method to map the composition domain while satisfying the atom balance. The Gibbs energy is calculated at each point of the composition domain and the absolute minimum can be deduced very simply. In theory, the technique is not limited, the Gibbs function needs not be discretised and multiphase problem can be handled easily. It is further shown that the accuracy of the predictions depends to a significant extent on the “coherence” of the input thermodynamic data such the formation Gibbs energy of the species and molecular interaction parameters. The absolute value of such parameters does not matter as much as their evolution as function of the process parameters (pressure, temperature, …). So, a self-consistent estimation method is required. To achieve this, the group contribution theory is used (UNIFAC descriptors) and extended somewhat outside the traditional molecular interaction domain, for instance to predict the Gibbs energy of formation of the species, the specific heat capacity… Lastly the influence of the choice of the final list of products is discussed. It is shown that the relevancy of the prediction depends to a large extent on this initial choice. A first technique is proposed, based on Brignole and Gani‘s work, to avoid omitting species and another one to select, in this list, the products likely to appear given the process conditions. These techniques were programmed in a new code name CIRCE. Brignole and Gani-‘s method is implemented on the basis of the atomic composition of the reactants to predict all “realisable” molecules. The extended group contribution theory is implemented to calculate the thermodynamic parameters. The MCGE method is used to find the absolute minimum of the Gibbs energy function. The code seems to be more versatile than the traditional ones (CEA, ASPEN…) but more expensive in calculation costs. It can also be more predictive. Examples are shown illustrating the breadth of potential applications in chemical engineering.

For planning issues, it is useful to know the upper limit for photovoltaics (PV) in the electrical grid with current design and operation (defined as hosting capacity) and how this limit can be increased. Future costs for grid reinforcement can be avoided if measures are taken to implement smart grid technology in the distribution grid. The aim of this project is to identify challenges in Swedish electricity distribution grids with a high penetration of local generation of electricity from PV. The aim is also to help Swedish Distribution System Operators (DSOs) to better understand hosting capacity issues, and to see which room for PV integration there is before there is need for actions to maintain power quality. Three distribution grids are modelled and simulated in Matlab: Rural area, Residential area and City (Stockholm Royal Seaport). Since the project is a cooperation between Uppsala University and Fortum, three different representative grids from Fortum’s grid software ”Power Grid” have been used as input to a flexible simulation program developed at Uppsala University. The simulation includes Newton-Raphson power-flow computing but has also been improved with a model of the temperature dependency of the resistance. The results show that there is room for a lot of PV systems in the Swedish grids. When using voltage rise above 1.1 p.u. voltage as limitation, the hosting capacity 60% PV electricity generation as a fraction of the yearly load were determined for the rural grid and the suburban grid. For the city grid, which is very robust, the hosting capacity 325% was determined. When using overload as limitation, the hosting capacities 70%, 20% and 25%, were determined for the same grids.

The aim of the thesis is to create an analytical 2D model of a multi-storey building and its consequent loading until the point of collapse which occures due to formation of plastic hinges. The first part is going to present a linear analysis of the problem; it focuses on location and level of load when the plastic hinges are formed. Then a nonlinaer analysis is going to be performed by RFEM programme using postcritical analysis and dynamic relaxation. Differences between the results of mentioned types and methods of analysis are going to be compared and an impact of variables is going to be evaluated .

La thèse CIFRE s'est déroulée dans un contexte de métallurgie industrielle et de simulation numérique. La modélisation de la solidification, via l'équation de la chaleur et avec des méthodes de réduction de modèle, était un objectif majeur. L'entreprise Montupet, spécialisée dans la fonderie d'aluminium, est le porteur du projet et financeur de la thèse. L'université de Technologie de Compiègne (UTC) a réalisé l'accompagnement académique. La méthode PGD “Proper General Decomposition”, basée sur une séparation de variables, est l'objet de nombreuses recherches. Nous avons proposé, concernant des propriétés matériaux non-linéaires, une discrétisation spatio-temporelle des matrices matériaux. Avec une formulation en température, sans chaleur latente, les gains sont élevés. L'introduction de la chaleur latente réduit fortement les gains. Nous formulons l'hypothèse que la difficulté de convergence de la méthode PGD dans le cas de la solidification repose sur une formulation en température inadaptée. Nous décidons d'opter une formulation en enthalpie. Il s'avère que cette formulation offre des perspectives encourageantes, mais nécessite encore beaucoup de développements. En parallèle de ces développements, un séjour de 5 mois aux États-Unis a été réalisé afin d'obtenir une meilleure caractérisation de la fraction solide. La recherche s'est portée sur l'évolution de la fraction solide en fonction de la vitesse de refroidissement. Sous réserve de mesures complémentaires, les essais ont mis en évidence une modification de la courbe de fraction solide en fonction de la vitesse de refroidissement, notamment un agrandissement de l'intervalle de solidification. / The PhD Thesis was carried out in a metallurgy and numerical simulation environment. The main topic was to model solidification, thought heat equation formulation and reduced order model PGD resolution. Montupet, specialized in aluminium alloys foundry hold and financed the project, the Université de Technologie de Compiègne did the acadernic supervising. The PGD method "Proper General Decomposition" is a hot topic based on variable separation. We proposed, regarding the non-linear materials, a space-time discretization of material matrix. With a temperature formulation, without latent heat, gains are high. With latent heat, gains fall drastically. We proposed the hypothesis that temperature could be an inadapted formulation. We decided to use the enthalpy formulation. This formulation offers good perspectives but needs more developments. During the thesis, five months were spent in the USA to get a better caracterisation of the solid fraction, particularly its variation in function of the cooling rate. Under reservation, the samples show a modification of solid fraction curves and particularly a change of solidification interval in function of cooling rate.

The aim of the thesis is to compare the iterative methods which program RFEM 5 uses the non-linear calculations of structures, namely the analysis of large deformations and post critical analysis. Comparison should serve as a basis for which calculation method is the most accurate, fastest and most reliable in terms of getting results. Time-consuming will be judged according to the calculation of the solution and the time needed to compute one iterativ. Robustness we will compare the reliability of methods in in normal use. Accuracy of the calculation will be determined by comparing the maximum deformation structures. Comparison will be made with examples from practice.

Este trabajo de investigación busca mejorar la estimación de la fragmentación de roca volada en minería a cielo abierto, específicamente en la empresa Minera Aruntani S.A.C. Unidad Arasi. La investigación se lleva a cabo mediante una reestructuración a la ecuación de predicción del tamaño de rocas, Kuz-Ram, modelo tradicional, cambiando el coeficiente de tierra, permitiendo disminuir los márgenes de error tomando en consideración aspectos relevantes del diseño como los parámetros geológicos, así como parámetros de perforación y voladura. Para efectuar dicho logro, se planteó una función objetivo y su posterior solución mediante el método iterativo de Newton-Raphson desarrollado en un entorno del software Matlab. El pertinente planteamiento y desarrollo permitió obtener un nuevo valor para el coeficiente de tierra con el cual se logra disminuir los porcentajes de error en la estimación de fragmentación de 40%, con el método tradicional, a un 8% en promedio logrando así una mejora considerable. / This research work seeks to improve the prediction of the fragmentation of blasted rock in open pit mining, specifically in the company Minera Aruntani S.A.C. Arasi unit. The research is carried out through a modification to the Kuz-Ram fragmentation estimation model, a traditional model, changing the land coefficient, allowing to reduce the margins of error taking into consideration relevant aspects of the design such as geological parameters, as well as parameters of drilling and blasting. To carry out this achievement, an objective function was proposed and its subsequent solution by means of the iterative Newton-Raphson method developed in a Matlab software environment. The pertinent approach and development allowed obtaining a new value for the land coefficient with which it is possible to reduce the error percentages in the fragmentation estimation from 40%, with the traditional method, to 8% on average, achieving a considerable improvement. / Tesis

No description available.