No description available.

Este trabalho insere-se no âmbito das formulações não convencionais em elementos finitos. Particularmente, introduzem-se alguns aspectos do método dos elementos finitos generalizados (MEFG) e do clássico refino-p na consagrada formulação híbrida-Trefftz de tensão para a elasticidade bidimensional. A formulação apresentada aproxima diretamente dois campos independentes: o de tensões no domínio dos elementos e o de deslocamentos nas fronteiras dos elementos. Baseado na estrutura de enriquecimento centrada em nuvens, proposta pelo MEFG, podem ser selecionadas oportunamente regiões, formadas por um conjunto de elementos e fronteiras de elementos, onde o espaço da aproximação é adequadamente enriquecido mediante o refino-p. Neste contexto campos auto-equilibrados de tensões, derivados da solução da equação de Navier, são utilizados para compor a aproximação no domínio dos elementos, enquanto nas fronteiras dos elementos o campo de deslocamentos é construído a partir de bases específicas de aproximação; seja a base inicial, formada por funções de forma lineares, ou bases enriquecidas com polinômios hierárquicos, não hierárquicos e funções trigonométricas. Aborda-se também, ainda que preliminarmente, um estudo de painéis com múltiplas fissuras pelo método da partição em formulação híbrida-Trefftz com enriquecimento seletivo. As análises numéricas realizadas revelaram, em geral, uma formulação de ótimo desempenho, caracterizada por uma notável capacidade de aproximação dos campos de tensões e deslocamentos, elevada robustez numérica e reduzido dispêndio computacional. / This work is inserted in the context of unconventional formulations in the finite elements method. Particularly, some aspects of the generalized finite elements method (GFEM) and the classic p-refinement are introduced in the well known hybrid-Trefftz stress formulation for the two dimensional elasticity. The presented formulation approximates two independent fields: the one of stresses in the elements domain and the one of displacements in the boundaries of the elements. Based on the enrichment structure centered in clouds, proposed by the GFEM, some regions, formed by a group of elements and boundaries of elements where the approximation space is adequately enriched by the p-refinement, can be opportunely selected. In this context, self-equilibrated stress fields, derived from the solution of the Navier equation, are used to compose the approximation in the elements domain, whereas the displacements field in the borders of the elements is built from specific approximation bases, that is, the initial base formed by linear shape functions, or, bases enriched with hierarchical polynomials, nonhierarchical ones and trigonometric functions. Also, although preliminarily, a study of the multiple-cracked panels is done using the Splitting Method with a hybrid-Trefftz formulation and a selective enrichment. The numeric analyses done revealed, in general, a high performance formulation characterized by a great capacity of approximation the stress fields and displacements, high numeric robustness and reduced computer expenditure.

No description available.

Este trabalho insere-se no âmbito das formulações não convencionais em elementos finitos. Particularmente, introduzem-se alguns aspectos do método dos elementos finitos generalizados (MEFG) e do clássico refino-p na consagrada formulação híbrida-Trefftz de tensão para a elasticidade bidimensional. A formulação apresentada aproxima diretamente dois campos independentes: o de tensões no domínio dos elementos e o de deslocamentos nas fronteiras dos elementos. Baseado na estrutura de enriquecimento centrada em nuvens, proposta pelo MEFG, podem ser selecionadas oportunamente regiões, formadas por um conjunto de elementos e fronteiras de elementos, onde o espaço da aproximação é adequadamente enriquecido mediante o refino-p. Neste contexto campos auto-equilibrados de tensões, derivados da solução da equação de Navier, são utilizados para compor a aproximação no domínio dos elementos, enquanto nas fronteiras dos elementos o campo de deslocamentos é construído a partir de bases específicas de aproximação; seja a base inicial, formada por funções de forma lineares, ou bases enriquecidas com polinômios hierárquicos, não hierárquicos e funções trigonométricas. Aborda-se também, ainda que preliminarmente, um estudo de painéis com múltiplas fissuras pelo método da partição em formulação híbrida-Trefftz com enriquecimento seletivo. As análises numéricas realizadas revelaram, em geral, uma formulação de ótimo desempenho, caracterizada por uma notável capacidade de aproximação dos campos de tensões e deslocamentos, elevada robustez numérica e reduzido dispêndio computacional. / This work is inserted in the context of unconventional formulations in the finite elements method. Particularly, some aspects of the generalized finite elements method (GFEM) and the classic p-refinement are introduced in the well known hybrid-Trefftz stress formulation for the two dimensional elasticity. The presented formulation approximates two independent fields: the one of stresses in the elements domain and the one of displacements in the boundaries of the elements. Based on the enrichment structure centered in clouds, proposed by the GFEM, some regions, formed by a group of elements and boundaries of elements where the approximation space is adequately enriched by the p-refinement, can be opportunely selected. In this context, self-equilibrated stress fields, derived from the solution of the Navier equation, are used to compose the approximation in the elements domain, whereas the displacements field in the borders of the elements is built from specific approximation bases, that is, the initial base formed by linear shape functions, or, bases enriched with hierarchical polynomials, nonhierarchical ones and trigonometric functions. Also, although preliminarily, a study of the multiple-cracked panels is done using the Splitting Method with a hybrid-Trefftz formulation and a selective enrichment. The numeric analyses done revealed, in general, a high performance formulation characterized by a great capacity of approximation the stress fields and displacements, high numeric robustness and reduced computer expenditure.

We develop and analyse explicit in time schemes for the computation of compressible flows, based on staggered in space. Upwinding is performed equation by equation only with respect to the velocity. The pressure gradient is built as the transpose of the natural divergence. For the barotropic Euler equations, the velocity convection is built to obtain a discrete kinetic energy balance, with residual terms which are non-negative under a CFL condition. We then show that, in 1D, if a sequence of discrete solutions converges to some limit, then this limit is the weak entropy solution. For the full Euler equations, we choose to solve the internal energy balance since a discretization of the total energy is rather unnatural on staggered meshes. Under CFL-like conditions, the density and internal energy are kept positive, and the total energy cannot grow. To obtain correct weak solutions with shocks satisfying the Rankine-Hugoniot conditions, we establish a kinetic energy identity at the discrete level, then choose the source term of the internal energy equation to recover the total energy balance at the limit. More precisely speaking, we prove that in 1D, if we assume the L∞ and BV-stability and the convergence of the scheme, passing to the limit in the discrete kinetic and internal energy equations, we show that the limit of the sequence of solutions is a weak solution. Finally, we consider the computation of radial flows, governed by Euler equations in axisymetrical (2D) or spherical (3D) coordinates, and obtain similar results to the previous sections. In all chapters, we show numerical tests to illustrate for theoretical results. / We develop and analyse explicit in time schemes for the computation of compressible flows, based on staggered in space. Upwinding is performed equation by equation only with respect to the velocity. The pressure gradient is built as the transpose of the natural divergence. For the barotropic Euler equations, the velocity convection is built to obtain a discrete kinetic energy balance, with residual terms which are non-negative under a CFL condition. We then show that, in 1D, if a sequence of discrete solutions converges to some limit, then this limit is the weak entropy solution. For the full Euler equations, we choose to solve the internal energy balance since a discretization of the total energy is rather unnatural on staggered meshes. Under CFL-like conditions, the density and internal energy are kept positive, and the total energy cannot grow. To obtain correct weak solutions with shocks satisfying the Rankine-Hugoniot conditions, we establish a kinetic energy identity at the discrete level, then choose the source term of the internal energy equation to recover the total energy balance at the limit. More precisely speaking, we prove that in 1D, if we assume the L∞ and BV-stability and the convergence of the scheme, passing to the limit in the discrete kinetic and internal energy equations, we show that the limit of the sequence of solutions is a weak solution. Finally, we consider the computation of radial flows, governed by Euler equations in axisymetrical (2D) or spherical (3D) coordinates, and obtain similar results to the previous sections. In all chapters, we show numerical tests to illustrate for theoretical results.

No description available.

A novel system identification (SI)-based structural health assessment (SHA) procedure has been developed integrating several theoretical and implementation aspects. The procedure assesses health of structures using limited noise-contaminated dynamic responses and without using input excitation information. Since most practical structures are three dimensional (3D), the procedure has been developed for general 3D structures, represented by finite elements (FEs). The procedure identifies defects by tracking the changes in the stiffness of the elements in the FE representation. Once a defective element is identified, defect spot can be identified accurately within the defective element. The procedure is denoted as 3D Generalized Iterative Least-Squares Extended Kalman Filter with Unknown Input (3D GILS-EKF-UI) and implemented in two stages. In Stage 1, based on the available responses, substructure(s) are selected and the 3D GILS-UI procedure is used to generate the unknown input excitation, stiffness parameters of the elements in the substructure, and two Rayleigh damping coefficients. Using information from Stage 1, stiffness parameters for the whole structure are identified using EKF with Weighted Global Iteration (EKF-WGI) in Stage 2. The procedure accurately identified defect-free and defective states of various 3D structures using only analytically generated limited responses. To increase the robustness, 3D GILS-EKF-UI has been extended to develop an integrated structural health assessment strategy, denoted as Iterative Least-Squares Extended Kalman Filter with Unknown Input and Advanced Digital Integration Technique (ILS-EKF-UI-ADIT). The procedure has been implemented in three stages. In Stage 1, an advanced digital integration technique (ADIT) is implemented for post-processing of noise-contaminated acceleration time-histories, addressing all major challenges of digital integration. It also overcomes non-convergence issue in Stage 2 that arises due to phase-shift and amplitude errors. In Stage 2, substructure(s) are identified using the least-squares procedure. In Stage 3, stiffness parameters for the whole structure are identified using the EKF-WGI procedure. ILS-EKF-UI-ADIT has been verified in presence of relatively large noise in the acceleration time-histories, measured at small part(s) of defect-free and defective structures, without using excitation information. The SHA procedure is robust and has the potential to be applied for the health assessment, maintenance, retrofitting, and life extension of existing structural systems.

Interpolation error estimates in terms of geometric quality measures are established for harmonic coordinates on polytopes in two and three dimensions. First we derive interpolation error estimates over convex polygons that depend on the geometric quality of the triangles in the constrained Delaunay triangulation of the polygon. This characterization is sharp in the sense that families of polygons with poor quality triangles in their constrained Delaunay triangulations are shown to produce large error when interpolating a basic quadratic function. Non-convex polygons exhibit a similar limitation: large constrained Delaunay triangles caused by vertices approaching a non-adjacent edge also lead to large interpolation error. While this relationship is generalized to convex polyhedra in three dimensions, the possibility of sliver tetrahedra in the constrained Delaunay triangulation prevent the analogous estimate from sharply reflecting the actual interpolation error. Non-convex polyhedra are shown to be fundamentally different through an example of a family of polyhedra containing vertices which are arbitrarily close to non-adjacent faces yet the interpolation error remains bounded.

L'ostéoporose est un problème majeur de santé publique (150000 fractures par an en France), qui se traduit par une diminution de la masse et de la « qualité » osseuse entraînant une augmentation du risque fracturaire. Basés uniquement sur une mesure de densité minérale osseuse, les moyens cliniques actuels de diagnostic de cette pathologie apparaissent insuffisants pour quantifier précisément la résistance mécanique de l'os. Dans cette optique, il est primordial de pouvoir évaluer non seulement l'évolution de masse osseuse au cours du temps, mais également ses propriétés mécaniques tant macro que micro. Cette thèse est consacrée au développement d'un modèle numérique de remodelage osseux intégrant des paramètres mécano-biologiques permettant de simuler à la fois la dégradation d'une microstructure d'os trabéculaire au cours du temps ainsi que son adaptation mécanique. Une modification de ces paramètres permet de simuler des cas pathologiques avec ou sans traitements anti-ostéoporotiques. Le remodelage est ensuite couplé au processus de mécanotransduction, en intégrant l'influence du chargement mécanique sur la réponse cellulaire. Différents cas de chargement ont été ainsi modélisés, en particulier faible chargement tel que l'alitement, ou surchargement tel que la pratique d'activités sportives. L'hyperminéralisation ainsi que l'hétérogénéité du tissu osseux ont été intégrés au modèle et leur influence sur la localisation des sites de remodelage a été également analysée. Les résultats des simulations obtenus à partir de microarchitectures d'os humain sur plusieurs années sont encourageants car ils concordent avec les observations cliniques. / Osteoporosis is characterized by a low bone mass density but also an alteration of mechanical properties. The clinical diagnostic is made from the measure of the bone mineral density (BMD) but this examen seems insufficient to quantify bone resistance. In this work, a numerical model of cancellous bone degradation, aging and mechanical adaptation is proposed. Based on hypermineralization, this model simulates the cancellous bone remodeling process over many years. This model allows to predict the behavior of cancellous adaptation in a mechanical low loading case for instance. Results are similar to clinicial tendancy.

Comme son nom l'indique, le tolérancement flexible a pour objectif de tenir compte de la souplesse des pièces dans un processus de tolérancement. Il permet d'évaluer les défauts géométriques admissibles par des critères aussi bien géométriques que mécaniques. Ces travaux abordent en premier lieu l'élaboration de modèles adéquats permettant de prédire le comportement mécanique d'un assemblage de grandes pièces relativement souples lorsqu'elles sont sujettes à des défauts géométriques issues du procédé de fabrication. Une méthode a alors été proposée pour y implémenter des variations géométriques aléatoires réalistes vis-à-vis de ces hypothétiques défauts géométriques. Pour simuler les opérations d'assemblage, le phénomène de contact unilatéral et les variations de rigidité dues aux variabilités géométriques ont été prises en compte. En fonction de ces hypothèses, les stratégies d'analyse de tolérance avec Monte Carlo ou la méthode des coefficients d'influence ont été comparées afin de choisir celle minimisant les coûts de calcul tout en conservant la justesse des résultats. Finalement, ces travaux s'achèvent sur une étude de cas industriel : un treillis supportant des équipements sous le plancher du fuselage d'un Airbus A350. / As indicated by its name, the purpose of flexible tolerancing is to take into account the flexibility of the parts in a tolerancing process. It allows to evaluate the permissibles geometrical defects by both geometrical and mechanical criteria. These works deal first with the elaboration of a model able to predict the mechanical behaviour of an assembly of larges and flexibles parts which are suject to geometrical defects induced by their manufacturing process. A method has thus been suggested in order to implement geometrical variations which are realistic according to these hypothetical geometrical defects. To simulate the assembling process, the unilateral contact phenomenon and the stiffness variations induced by the geometrical variability have been taken into account. Depending on these assumptions, some strategies of tolerance analysis with Monte Carlo or the method of influence coefficients have been compared in order to choose the one minimizing the computational costs while maintaining the accuracy of the results. Finally, these works are completed with an industrial study case: a truss supporting equipments and hanging under the fuselage floor of an Airbus A350.