Dans l’industrie pharmaceutique, la granulation sèche par compactage à rouleaux est un procédé d’agglomération de poudres en granulés pour améliorer les propriétés d’écoulement nécessaire pour le procédé de compression en matrice. Comprendre le procédé de compactage à rouleaux et optimiser l’efficacité de production est limitée par l’utilisation de l’approche expérimentale à cause du coût élevé des poudres, le temps des essais et la complexité du procédé. Dans ce travail, une méthode d’éléments finis en 3D, est développée dans le but d’identifier les paramètres critiques du matériau et du procédé pour le contrôle de la qualité de la production. Le modèle de comportement de Drucker-Prager Cap est utilisé pour décrire le comportement en compression de poudres et sa calibration est déterminée à partir des essais standard. Pour surmonter la complexité liée à l’existence de deux mécanismes différents, l’alimentation en poudre par une vis sans fin et le compactage entre les rouleaux, une nouvelle méthode d’interfaçage entre la méthode des éléments discrets (DEM) employée pour décrire l’écoulement dans l’alimentation et la méthode des éléments finis (FEM) utilisée pour le compactage entre les rouleaux est développée. Enfin, pour une modélisation de compactage de rouleaux plus réaliste, prenant en compte la variation de l’entrefer entre les rouleaux, une nouvelle approche de couplage Euler-Lagrange est proposée. Les résultats de simulations par éléments finis montrent clairement l’effet des différents paramètres du procédé sur les distributions de pression et de densité dans la zone de compactage. En outre, les résultats montrent que l'utilisation de plaques de confinement de la poudre entre les rouleaux, développe une distribution de pression et de densité non homogène dans le compact, avec une densité plus élevée au centre et plus faible aux bords. D'autre part, l’utilisation de rouleaux dont l’un est surmonté d’une jante de confinement, a montré une distribution de propriétés globalement plus uniforme sur la largeur du compact avec des valeurs légèrement plus élevées aux bords qu’au centre. La méthodologie combinant les méthodes DEM & FEM montre clairement une corrélation directe entre la vitesse des particules entraînées par la vis dans la zone d’alimentation et la pression du rouleau. Tous les deux oscillent avec la même période. Cela se traduit par un compact anisotrope avec un profile de densité variant de manière sinusoïdale le long de sa largeur. Afin d'étudier la capacité du modèle à prédire les propriétés des compacts produits par compactage à rouleaux, les prédictions par simulations numériques sont comparées aux données de la littérature et validées par des mesures spécifiques. / In the pharmaceutical industry, dry granulation by roll compaction is a process of size enlargement of powder into granules with good flowability for subsequent die compaction process. Understanding the roll compaction process and optimizing manufacturing efficiency is limited using the experimental approach due to the high cost of powder, time-consuming and the complexity of the process. In this work, a 3D Finite Element Method (FEM) model was developed to identify the critical material properties, roll press designs and process parameters controlling the quality of the product. The Drucker-Prager Cap (DPC) model was used to describe the powder compaction behavior and was determined based on standard calibration method. To overcome the complexity involving two different mechanisms of powder feeding by the screw and powder compaction between rolls, a novel combined approach of Discrete Element Method (DEM), used to predict the granular material flow in the feed zone and the Finite Elements Method (FEM) employed for roll compaction, was developed. Lastly, for a more realistic roll compaction modelling, allowing the fluctuation of the gap between rolls, a Coupled-Eulerian Lagrangian (CEL) approach was developed. FEM simulation results clearly show the effect of different process parameters on roll pressure and density distribution in the compaction zone of powder between the rolls. Moreover, results show that using a cheek-plates sealing system causes a nonuniform roll pressure and density distribution with the highest values in the middle and the lowest at the edges. On the other hand, the resultant pressure and density distributions with the rimmed-roll obtained higher values in the edges than in the middle and overall a more uniform distribution. The combined DEM-FEM methodology clearly shows a direct correlation between the particle velocity driven by the screw conveyor to the feed zone and the roll pressure, both oscillating in the same period. This translates into an anisotropic ribbon with a density profile varying sinusoidally along its length. To validate the results, the simulations are compared with literature and experimentally measured values in order to assess the ability of the model to predict the properties of the produced ribbons.

This work is dedicated to faulty induction motors. These motors are often used in industrial applications thanks to their usability and their robustness. However, nowadays optimisation of production becomes so critical that the conceptual reliability of the motor is not sufficient anymore. Motor condition monitoring is expanding to serve maintenance planning and uptime maximisation. Moreover, the use of drive control sensors (namely stator current and voltage) can avoid the installation and maintenance of dedicated sensors for condition monitoring.<p><p>Many authors are working in this field but few approach the diagnosis from a detailed and clear physical understanding of the localised phenomena linked to the faults. Broken bars are known to modulate stator currents but it is shown in this work that it also changes machine saturation level in the neighbourhood of the bar. Furthermore, depending on the voltage level, this change in local saturation affects the amplitude and the phase of the modulation. This is of major importance as most diagnosis techniques use this feature to detect and quantify broken bars. For stator short-circuits, a high current is flowing in the short-circuited coil due to mutual coupling with the other windings and current spikes are flowing in the rotor bars as they pass in front of the short-circuited conductors. In the case of rotor eccentricities, the number of pole-pairs and the connection of these pole-pairs greatly affect the airgap flux density distribution as well as the repartition of the line currents in the different pole-pairs.<p><p>These conclusions are obtained through the use of time-stepping finite element models of the faulty motors. Moreover, circuit models of faulty machines are built based on the conclusions of previously explained fault analysis and on classical Park models. A common mathematical description is used which allows objective comparison of the models for representation of the machine behaviour and computing time.<p><p>The identifiability of the parameters of the models as well as methods for their identification are studied. Focus is set on the representation of the machine behaviour using these parameters more than the precise identification of the parameters. It is shown that some classical parameters can not be uniquely identified using only stator measurements.<p><p>Fault detection and identification using computationally cheap models are compared to advanced detection through motor stator current spectral analysis. This last approach allows faster detection and identification of the fault but leads to incorrect conclusions in low load conditions, in transient situations or in perturbed environments (i.e. fluctuating load torque and unideal supply). Efficient quantification of the fault can be obtained using detection techniques based on the comparison of the process to a model.<p><p>Finally, the work provides guidelines for motor supervision strategies depending on the context of motor utilisation. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Excavations in soft clays come with a number of specific challenges. The use of lime-cement columns in the passive zone of such works represents a promising solution to some of these problems, such as basal heave stability or excessive deformations. Nevertheless, lime-cement columns in the passive zone of excavations in soft clays are not yet widely used and studied. Knowledge about this improvement method is still limited, particularly when it comes to numerical modelling. They have mostly been studied numerically using simple geometries or constitutive models that do not allow to represent the key features of the material, that could influence the behaviour of the structure. Thus, there is a need for more investigations regarding precise modelling tools for this type of problems. In this study, an advanced constitutive model named the Concrete model was employed to capture the mechanical behaviour of the lime-cement improved soil material. A boundary value problem representing an excavation in soft clay supported by sheet pile walls and lime-cement columns was studied using a three dimensional finite element model that included the Concrete model for lime-cement column panels. The results show that the Concrete model is able to capture the behaviour of lime-cement improved clay very well. Key features such as strain hardening and strain softening are well represented. The results of the boundary value problem were compared to field measurements and it was showed that the Concrete model employed in a three dimensional finite element representation of the problem is able to give realistic results. / Schakter i lös lera medför ett antal specifika utmaningar. Användningen av kalk-cementpelare i den passiva zonen i sådana fallen är en lovande lösning på några av dessa problem, t.ex. jordhävning eller stora deformationer. Kalk-cement pelare i den passiva zonen av schakter i lös lera har dock ännu inte använts och studerats i någon större utsträckning i Sverige. Kunskapen om denna metod är fortfarande begränsad, särskilt när det gäller numerisk modellering. Denna förstärkningsmetod har studerats numeriskt med hjälp av enkla geometrier eller konstitutiva modeller som inte gör det möjligt att representera materialets viktigaste egenskaper som kan påverka konstruktionens beteende. Det finns därför ett behov av mer forskning för en precis modellering för denna typ av problem. I den här studien användes en avancerad konstitutiv modell (Concrete model) för att simulera det mekaniska beteendet hos det kalkcementstabiliserade jordmaterialet. Ett gränsvärdesproblem som representerar en schakt i lös lera som stöds av spontväggar och kalkcementpelare studerades med hjälp av en tredimensionell finita elementmodell. Resultaten visar att den advancerade konstitutiva modellen kan simulera beteendet hos kalkcementstabiliserad lera mycket väl. Viktiga egenskaper som t.ex. töjningshärdning och töjningsuppmjukning är väl representerade. Resultaten av gränsvärdesproblemet modellering jämfördes med fältmätningar och det visades att konstitutiva modellen som används kan ge realistiska resultat. / Les excavations dans les argiles molles présentent de nombreuses difficultés. L’utilisation de colonnes en chaux-ciment dans la zone passive de ces structures représente une solution prometteuse à certains de ces défis, tels que le soulèvement de la base ou les déformations excessives. Néanmoins, l’utilisation et l’étude des colonnes en chaux-ciment dans la zone passive des excavations dans les argiles molles ne sont pas généralisées et les connaissances à ce sujet sont encore limitées, notamment en ce qui concerne la modélisation numérique. La plupart des études numériques se basent sur des géométries simples ou des modèles constitutifs ne permettant pas de représenter certaines particularités du comportement du matériau, qui pourraient influencer le comportement de la structure. Il est donc nécessaire d’approfondir les recherches sur les outils de modélisation pour ce type de problème. Dans cette étude, un modèle constitutif avancé (Concrete model) a été utilisé pour reproduire le comportement mécanique de l’argile stabilisée à la chaux et au ciment. Un problème de valeurs aux limites représentant une excavation dans de l’argile molle soutenue par des murs de palplanches et des colonnes en chaux-ciment a été étudié à l’aide d’un modèle d’éléments finis tridimensionnel. Les résultats montrent que le modèle employé permet de reproduire le comportement de l’argile stabilisée à la chaux-ciment de manière très satisfaisante. Les particularités du comportement du matériau, telles que l’écrouissage et l’adoucissement, sont bien représentées. Les résultats du problème de valeurs aux limites ont été comparés aux mesures sur le terrain et il semble que ce modèle constitutif, employé dans une représentation tridimensionnelle par éléments finis du problème, soit capable de donner des résultats réalistes.

L'objet de cette thèse est de présenter une étude du comportement thermomécanique des disques de frein automobiles pour la prédiction de leur tenue en fatigue. La stratégie de calcul numérique est repose sur le code de calcul Ansys v.11. Ce dernier qui est basé sur la méthode des éléments finis et qui possède des algorithmes de gestion du contact avec frottement est utilisé pour simuler dans l'application du freinage le comportement du mécanisme malgré son complexité. Dans un premier temps, est présentée une analyse des phénomènes thermiques opérant dans un disque de frein en service (flux de chaleur généré par frottement, gradients thermiques élevés, élévation de température). Cette modélisation est effectuée en tenant compte l'influence d'un certains nombre de paramètre tel que le type de freinage, le mode de refroidissement, les matériaux de conceptions. Ensuite, une étude purement mécanique du contact sec entre le disque et plaquettes est développée avec une bonne prédiction devient un enjeu majeur pour les industriels tout en modélisant le chargement et les conditions aux limites autours du disque. Nous avons utilisé le même code de calcul pour visualiser les déplacements, les déformations globales dans le disque, les contraintes de cisaillement, les contraintes de Von Mises et les outils de contact des plaquettes tout en effectuant une étude paramétrique telle que ( le module de Young des plaquettes, le coefficient de frottement , le type de chargement , la vitesse de rotation du disque ,..) pour voir sa sensibilité sur les résultats de calcul. Ainsi, les analyses faites sur le comportement thermique et mécanique que de ces prototypes montrent que ces types de solutions technologiques représentent de réelles pistes d'amélioration qui répond au besoin de l'ingénieur en charge de la conception des disques de frein.

L’objectif principal de cette thèse consiste à étudier plusieurs thématiques : l’étude de l’observation et la commande d’un système de structure flexible et l’étude de la stabilité asymptotique d’un système d’échangeurs thermiques. Ce travail s’inscrit dans le domaine du contrôle des systèmes décrits par des équations aux dérivées partielles (EDP). On s’intéresse au système du corps-poutre en rotation dont la dynamique est physiquement non mesurable. On présente un observateur du type Luenberger de dimension infinie exponentiellement convergent afin d’estimer les variables d’état. L’observateur est valable pour une vitesse angulaire en temps variant autour d’une constante. La vitesse de convergence de l’observateur peut être accélérée en tenant compte d’une seconde étape de conception. La contribution principale de ce travail consiste à construire un simulateur fiable basé sur la méthode des éléments finis. Une étude numérique est effectuée pour le système avec la vitesse angulaire constante ou variante en fonction du temps. L’influence du choix de gain est examinée sur la vitesse de convergence de l’observateur. La robustesse de l’observateur est testée face à la mesure corrompue par du bruit. En mettant en cascade notre observateur et une loi de commande stabilisante par retour d’état, on souhaite obtenir une stabilisation globale du système. Des résultats numériques pertinents permettent de conjecturer la stabilité asymptotique du système en boucle fermée. Dans la seconde partie, l’étude est effectuée sur la stabilité exponentielle des systèmes d’échangeurs thermiques avec diffusion et sans diffusion. On établit la stabilité exponentielle du modèle avec diffusion dans un espace de Banach. Le taux de décroissance optimal du système est calculé pour le modèle avec diffusion. On prouve la stabilité exponentielle dans l’espace Lp pour le modèle sans diffusion. Le taux de décroissance n’est pas encore explicité dans ce dernier cas. / The main objective of this thesis consists to investigate the following themes : observation and control of a flexible structure system and asymptotic stability of a heat exchangers system. This work is placed in the field of the control of systems described by partial differential equations (PDEs). We consider a rotating body-beam system whose dynamics are not physically measurable. An infinite-dimensional exponentially convergent Luenberger-like observer is presented in order to estimate the state variables. The observer is also valid for a time-varying angular velocity around some constant. We can accelerate the decay rate of the observer by a second step design. The main contribution of this work consists in building a numerical simulator based on the finite element method (FEM). A numerical investigation is carried out for the system with constant or time-varying angular velocity. We examine the influence of the gain choice on the decay rate of the observer. The robustness of the observer is tested with the measurement corrupted by noise. By cascading our observer and a feedback control law, we wish to obtain a global stabilization of the rotating bodybeam system. The relevant numerical results make it possible for us to conjecture that the closed-loop system is locally asymptotically stable. We investigate the exponential stability of the heat exchangers systems with diffusion or without diffusion. We establish the exponential stability of the model with diffusion in a Banach space. Moreover, the optimal decay rate of the system is computed for the model with diffusion. We prove exponential stability in (C[0, 1])4 space for the model without diffusion. The optimal decay rate in the latter case is not yet found.

Gliomas are the most common primary brain tumours and are associated with poor prognosis. Among them, diffuse gliomas – which include their most aggressive form glioblastoma (GBM) – are known to be highly infiltrative. The diagnosis and follow-up of gliomas rely on positron emission tomography (PET) and magnetic resonance imaging (MRI). However, these imaging techniques do not currently allow to assess the whole extent of such infiltrative tumours nor to anticipate their preferred invasion patterns, leading to sub-optimal treatment planning. Mathematical tumour growth modelling has been proposed to address this problem. Reaction-diffusion tumour growth models, which are probably the most commonly used for diffuse gliomas growth modelling, propose to capture the proliferation and migration of glioma cells by means of a partial differential equation. Although the potential of such models has been shown in many works for patient follow-up and therapy planning, only few limited clinical applications have seemed to emerge from these works. This thesis aims at revisiting reaction-diffusion tumour growth models using state-of-the-art medical imaging and data processing technologies, with the objective of integrating multi-parametric PET/MRI data to further personalise the model. Brain tissue segmentation on MR images is first addressed with the aim of defining a patient-specific domain to solve the model. A previously proposed method to derive a tumour cell diffusion tensor from the water diffusion tensor assessed by diffusion-tensor imaging (DTI) is then implemented to guide the anisotropic migration of tumour cells along white matter tracts. The use of dynamic [S-methyl-11C]methionine ([11C]MET) PET is also investigated to derive patient-specific proliferation potential maps for the model. These investigations lead to the development of a microscopic compartmental model for amino acid PET tracer transport in gliomas. Based on the compartmental model results, a novel methodology is proposed to extract parametric maps from dynamic [11C]MET PET data using principal component analysis (PCA). The problem of estimating the initial conditions of the model from MR images is then addressed by means of a translational MRI/histology study in a case of non-operated GBM. Numerical solving strategies based on the widely used finite difference and finite element methods are finally implemented and compared. All these developments are embedded within a common framework allowing to study glioma growth in silico and providing a solid basis for further research in this field. However, commonly accepted hypothesis relating the outlines of abnormalities visible on MRI to tumour cell density iso-contours have been invalidated by the translational study carried out, leaving opened the questions of the initialisation and the validation of the model. Furthermore, the analysis of the temporal evolution of real multi-treated glioma patients demonstrates the limitations of the formulated model. These latter statements highlight current obstacles to the clinical application of reaction-diffusion tumour growth models and pave the way to further improvements. / Les gliomes sont les tumeurs cérébrales primitives les plus communes et sont associés à un mauvais pronostic. Parmi ces derniers, les gliomes diffus – qui incluent la forme la plus agressive, le glioblastome (GBM) – sont connus pour être hautement infiltrants. Le diagnostic et le suivi des gliomes s'appuient sur la tomographie par émission de positons (TEP) ainsi que l'imagerie par résonance magnétique (IRM). Cependant, ces techniques d'imagerie ne permettent actuellement pas d'évaluer l'étendue totale de tumeurs aussi infiltrantes ni d'anticiper leurs schémas d'invasion préférentiels, conduisant à une planification sous-optimale du traitement. La modélisation mathématique de la croissance tumorale a été proposée pour répondre à ce problème. Les modèles de croissance tumorale de type réaction-diffusion, qui sont probablement les plus communément utilisés pour la modélisation de la croissance des gliomes diffus, proposent de capturer la prolifération et la migration des cellules tumorales au moyen d'une équation aux dérivées partielles. Bien que le potentiel de tels modèles ait été démontré dans de nombreux travaux pour le suivi des patients et la planification de thérapies, seules quelques applications cliniques restreintes semblent avoir émergé de ces derniers. Ce travail de thèse a pour but de revisiter les modèles de croissance tumorale de type réaction-diffusion en utilisant des technologies de pointe en imagerie médicale et traitement de données, avec pour objectif d'y intégrer des données TEP/IRM multi-paramétriques pour personnaliser davantage le modèle. Le problème de la segmentation des tissus cérébraux dans les images IRM est d'abord adressé, avec pour but de définir un domaine propre au patient pour la résolution du modèle. Une méthode proposée précédemment permettant de dériver un tenseur de diffusion tumoral à partir du tenseur de diffusion de l'eau évalué par imagerie DTI a ensuite été implémentée afin de guider la migration anisotrope des cellules tumorales le long des fibres de matière blanche. L'utilisation de l'imagerie TEP dynamique à la [S-méthyl-11C]méthionine ([11C]MET) est également investiguée pour la génération de cartes de potentiel prolifératif propre au patient afin de nourrir le modèle. Ces investigations ont mené au développement d'un modèle compartimental pour le transport des traceurs TEP dérivés des acides aminés dans les gliomes. Sur base des résultats du modèle compartimental, une nouvelle méthodologie est proposée utilisant l'analyse en composantes principales pour extraire des cartes paramétriques à partir de données TEP dynamiques à la [11C]MET. Le problème de l'estimation des conditions initiales du modèle à partir d'images IRM est ensuite adressé par le biais d'une étude translationelle combinant IRM et histologie menée sur un cas de GBM non-opéré. Différentes stratégies de résolution numérique basées sur les méthodes des différences et éléments finis sont finalement implémentées et comparées. Tous ces développements sont embarqués dans un framework commun permettant d'étudier in silico la croissance des gliomes et fournissant une base solide pour de futures recherches dans le domaine. Cependant, certaines hypothèses communément admises reliant les délimitations des anormalités visibles en IRM à des iso-contours de densité de cellules tumorales ont été invalidée par l'étude translationelle menée, laissant ouverte les questions de l'initialisation et de la validation du modèle. Par ailleurs, l'analyse de l'évolution temporelle de cas réels de gliomes multi-traités démontre les limitations du modèle. Ces dernières affirmations mettent en évidence les obstacles actuels à l'application clinique de tels modèles et ouvrent la voie à de nouvelles possibilités d'amélioration. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished