Les réactions de catalyse hétérogène impliquant un gaz et un liquide sont mises en oeuvre dans des réacteurs à lit fixes. Ces réacteurs peuvent être assimilés à un milieu poreux. La nature complexe de ce milieu rend la compréhension des interactions entre phases difficile, et nécessite une étude exhaustive à l’échelle globale et locale afin d’identifier les paramètres clés de l’hydrodynamique, des transferts de chaleur et de matière. Nous avons donc développé une cellule miniaturisée bidimensionnelle transparente, qui permet l’observation directe des écoulements avec une très bonne résolution spatiale et temporelle. En faisant varier le débit total et le rapport des débits gaz/liquide, nous avons mis en évidence l’apparition des régimes ruisselant et pulsé, observés par ailleurs dans des systèmes tridimensionnels. Grâce à une méthode d’analyse d’image, nous sommes capables de quantifier et cartographier la saturation liquide locale apparente et la morphologie des phases. L’analyse des variances a permis d’étudier les transitions de régimes pour différentes propriétés de la phase liquide. Cette approche a permis de comparer avec la littérature, mais aussi de s’intéresser aux mécanismes de génération et propagation des instabilités lors des transitions. Il a été mis en évidence que l’apparition des instabilités responsables de la déstabilisation du régime ruisselant intervient pour un nombre deWeber liquide fixé, indiquant que le régime pulsé apparait suite à la déstabilisation des interfaces gaz /liquides par les forces inertielles. Enfin, une étude préliminaire des transferts thermiques dans la cellule a été réalisée. La cellule est utilisée pour réaliser la réaction exothermique d’hydrogénation de l’alpha-methylstyrène. Un modèle simple de transferts thermique a été utilisé pour caractériser l’augmentation de température dans la cellule. Bien que cette approche ne permette pas des mesures quantitatives, elle ouvre à la perspective de suivi de transferts thermiques par caméra infra-rouge. / Heterogeneously catalyzed reactions involving a gas and a liquid phase are frequently achieved in fixed bed reactors. These reactors can be described as a porous medium. The complex nature of this medium makes the understanding of the interplay between phases difficult, and requires a thorough study at the global andlocal scale to identify the key parameters of hydrodynamics, mass and heat transfers. Therefore, we have developed a miniaturized two-dimensional system that is transparent to allow the direct observation of the flow with very high spatial and temporal resolution. While varying the total flow rate and the gas/liquid flow rate ratio, we observe the appearance of the trickle and the pulse flow regime, which can be observed in threedimensional beds. Thanks to some image analysis techniques, we are able to quantify and to map the local apparent liquid saturation and the morphology of the phases. Variances analysis allowed the study of the transition for different liquid properties. This approach allowed the comparison with the existing state of art, but also the study of the onset and propagation mechanisms of the instabilities during the transition. We report that the onset of instabilities responsible for the destabilization of the trickle flow regime occurs at a fixed Weber number. This indicates that the pulsed regime is due to the destabilization of the gas/liquid interface by inertial forces. Finally, a preliminary study of thermal transfers in the device was realized. The device was used to perform the exothermic hydrogenation of alpha-methylstyrene. A simple model is used to characterize the temperature increase in the device. Although this approach does not allow quantitative measurements, it opens up the perspective of monitoring thermal transfers with an infra-red camera.

The structural optimization of a cantilever aircraft wing with curvilinear spars and ribs and stiffeners is described. The design concept of reinforcing the wing structure using curvilinear stiffening members has been explored due to the development of novel manufacturing technologies like electron-beam-free-form-fabrication (EBF3). For the optimization of a complex wing, a common strategy is to divide the optimization procedure into two subsystems: the global wing optimization which optimizes the geometry of spars, ribs and wing skins; and the local panel optimization which optimizes the design variables of local panels bordered by spars and ribs. The stiffeners are placed on the local panels to increase the stiffness and buckling resistance. The panel thickness, size and shape of stiffeners are optimized to minimize the structural weight. The geometry of spars and ribs greatly influences the design of stiffened panels. During the local panel optimization, the stress information is taken from the global model as a displacement boundary condition on the panel edges using the so-called "Global-Local Approach". The aircraft design is characterized by multiple disciplines: structures, aeroelasticity and buckling. Particle swarm optimization is used in the integration of global/local optimization to optimize the SpaRibs. The interaction between the global wing optimization and the local panel optimization is usually computationally expensive. A parallel computing technology has been developed in Python programming to reduce the CPU time. The license cycle-check method and memory self-adjustment method are two approaches that have been applied in the parallel framework in order to optimize the use of the resources by reducing the license and memory limitations and making the code robust. The integrated global-local optimization approach has been applied to subsonic NASA common research model (CRM) wing, which proves the methodology's application scaling with medium fidelity FEM analysis. Both the global wing design variables and local panel design variables are optimized to minimize the wing weight at an acceptable computational cost. The structural weight of the wing has been, therefore, reduced by 40% and the parallel implementation allowed a reduction in the CPU time by 89%. The aforementioned Global-Local Approach is investigated and applied to a composite panel with crack at its center. Because of composite laminates' heterogeneity, an accurate analysis of these requires very high time and storage space. In the presence of structural discontinuities like cracks, delaminations, cutouts etc., the computational complexity increases significantly. A possible alternative to reduce the computational complexity is the global-local analysis which involves an approximate analysis of the whole structure followed by a detailed analysis of a significantly smaller region of interest. We investigate here the performance of the global-local scheme based on the finite element method by comparing it to the traditional finite element method. To do so, we conduct a 2D structural analysis of a composite square plate, with a thin rectangular notch at its center, subjected to a uniform transverse pressure, using the commercial software ABAQUS. We show that the presence of the thin notch affects only the local response of the structure and that the size of the affected area depends on the notch length. We investigate also the effect of the notch shape on the response of the structure. Stiffeners attached to composite panels may significantly increase the overall buckling load of the resultant stiffened structure. Buckling analysis of a composite panel with attached longitudinal stiffeners under compressive loads is performed using Ritz method with trigonometric functions. Results are then compared to those from ABAQUS FEA for different shell elements. The case of composite panel with one, two, and three stiffeners is investigated. The effect of the distance between the stiffeners on the buckling load is also studied. The variation of the buckling load and buckling modes with the stiffeners' height is investigated. It is shown that there is an optimum value of stiffeners' height beyond which the structural response of the stiffened panel is not improved and the buckling load does not increase. Furthermore, there exist different critical values of stiffener's height at which the buckling mode of the structure changes. Next, buckling analysis of a composite panel with two straight stiffeners and a crack at the center is performed. Finally, buckling analysis of a composite panel with curvilinear stiffeners and a crack at the center is also conducted. ABAQUS is used for these two examples and results show that panels with a larger crack have a reduced buckling load. It is shown also that the buckling load decreases slightly when using higher order 2D shell FEM elements. A damage tolerance framework, EBF3PanelOpt, has been developed to design and analyze curvilinearly stiffened panels. The framework is written with the scripting language PYTHON and it interacts with the commercial software MSC. Patran (for geometry and mesh creation), MSC. Nastran (for finite element analysis), and MSC. Marc (for damage tolerance analysis). The crack location is set to the location of the maximum value of the major principal stress while its orientation is set normal to the major principal axis direction. The effective stress intensity factor is calculated using the Virtual Crack Closure Technique and compared to the fracture toughness of the material in order to decide whether the crack will expand or not. The ratio of these two quantities is used as a constraint, along with the buckling factor, Kreisselmeier and Steinhauser criteria, and crippling factor. The EBF3PanelOpt framework is integrated within a two-step Particle Swarm Optimization in order to minimize the weight of the panel while satisfying the aforementioned constraints and using all the shape and thickness parameters as design variables. The result of the PSO is used then as an initial guess for the Gradient Based Optimization using only the thickness parameters as design variables. The GBO is applied using the commercial software VisualDOC. / Ph. D.

Ce travail porte sur l’étude du risque d’amorçage de fissure en fretting/fatigue du contact aubes/disque de soufflante en alliage de titane Ti-6Al-4V d’un turboréacteur civil. De part les vibrations du moteur et le passage du flux d’air, le contact aube/disque est sujet à de microdéplacements: le fretting. Ce problème complexe concerne des zones en contact de faibles dimensions subissant des pressions de contact importantes et des micro-déplacements associées à des structures soumises à des chargements eux-mêmes complexes. Des essais de fretting simple, fretting précontraint et fretting fatigue sont réalisés en vue d’identifier les conditions d’amorçage de fissure pour plusieurs géométries présentant des gradients de contraintes différents. Les résultats expérimentaux mettent en évidence l’effet du gradient de contrainte sur l’amorçage. Ainsi, via des essais présentant une gamme étendue de gradient, il a été confirmé que les gradients de contraintes retardent l’initiation de fissure. Une approche non-locale couplée à un critère de fatigue multiaxiale et identifiée sur une unique condition expérimentale de fretting simple a permis la prédiction de l’ensemble des résultats expérimentaux. Ainsi, une méthodologie prédictive des durées de vie à amorçage a été proposée. Cette approche numérique intègre la multiaxialité du chargement, l’effet du gradient des contraintes, la réponse élastoplastique de l’interface et les limites de fatigue en traction/compression et torsion alternée du matériau. L’interface aube/disque grenaillée du moteur subit une alternance de sollicitations de faibles amplitudes générant l’amorçage de fissure, et de grandes amplitudes favorisant l’usure des zones en contact. Le phénomène de compétition usure/fissure a alors été étudié via les courbes en cloches (durée de vie en fonction de l’amplitude de déplacement) ainsi que l’influence des contraintes résiduelles de grenaillage. Une première conclusion sur l’effet bénéfique de l’usure a pu être apportée dans le cas d’un matériau sans traitement de surface. Si les contraintes résiduelles de grenaillage permettent d’apporter un gain notable sur les durées de vie pour des sollicitations de glissement partiel, à l’inverse, il a été observé que le matériau avec le traitement de surface de grenaillage présente des durées de vie plus faibles pour des sollicitations favorisant l’usure. À travers d’essais inédits avec un dispositif expérimental à trois vérins, le chargement complexe de l’interface aube/disque du moteur a pu être reproduit à l’échelle du laboratoire. Ainsi, la durabilité globale du contact représentatif de l’interface aube/disque de soufflante a été étudiée. Une synergie des sollicitations oligocycliques et polycycliques à l’interface a été constatée provoquant des phénomènes de réorganisation de contraintes résiduelles et l’accélération des cinétiques d’usure. Cette étude a montré l’importance de prendre en compte les sollicitations polycycliques dans le dimensionnement des structures aubes/disque. Ces sollicitations couplent les aspects de fatigue « fissuration » aux phénomènes de cumul de dommage et de cinétique d’usure des interfaces. / This work aims at studying the crack initiation risk of a blade/disk contact under fretting/fatigue loading. The fan stage of a civil engine is studied and the material used is a titanium alloy: Ti-6Al-4V. This complex issue concerns small contacts which are subjected to high contact pressures and micro-displacements. Thus, a complex multiaxial loading occurs on the structure. Fretting, static stressed fretting and fretting fatigue testing are performed in order to quantify the crack nucleation thresholds. In addition, several geometries which present different stress gradient values are used. The experimental results show a stress gradient effect on the crack initiation. Thus, using a wide range of stress gradient, it has been confirmed that stress gradient delays crack nucleation. A multiaxial fatigue criterion is used to determine the equivalent stress field under the contact. Then, a non-local approach, identified thanks to one fretting experimental condition, is employed in order to consider the stress gradient effect. This approach provides the prediction of the whole experimental results. In this way, a predictive method of the initiation fatigue life has been introduced. This numerical approach takes into account the multiaxial loading, the stress gradient effect, the plastic-elastic behavior of the interface and the fatigue strength limits of the material. The blade/disk contacts of the engine are shot-peened. During the flight, interfaces are subjected to low displacement amplitudes leading to crack initiation. On the other hand, during landing and take-off, contacts are submitted to high displacement amplitudes leading to the interface wear. As a consequence, competition between wear kinetic and nucleation one is studied and “bell curves” are plotted (i.e. fatigue life time versus displacement amplitude).Moreover, the effect of the shot-peening residual stresses on fatigue life time is observed. In case of un-treated material, a beneficial impact on the fatigue life time is observed due to wear process. Compressive residual stresses of shot-peening increase the fatigue life when fretting/fatigue on partial slip regime occurs. Nevertheless, a decrease of the fatigue life time is observed when wear process is activated by the gross slip condition. Finally, the complex loading of the engine blade/disk contact has been reproduced in the laboratory. To achieve this, a new testing, using an experimental machine composed of three hydraulic actuators is developed. Thus, the global life time of the representative interface of the fan stage is studied. Oligocyclic (low frequency) and polycyclic (high frequency) solicitations interact themselves and lead, firstly, to a rearrangement of the residual stresses, and secondly, to an increase of the wear kinetics. This research work highlights the interest to consider the polycyclic loadings to design the blade/disk structure. These solicitations reproduce the “cracking” fatigue phenomenon, the cumulative damages and the wear kinetics of the interface.

[pt] Esta dissertação revisita quatro modelos multiaxiais de dano por fadiga (SWTD, SWT, Findley e Crossland) e discute a sua aplicação na previsão do risco de nucleação de trinca sob solicitações de fadiga por fretting segundo as abordagens local e não-local. As previsões do risco de nucleação de trinca obtidas através dos modelos estudados foram confrontadas com dados experimentais disponíveis na literatura para o contato cilindro/plano em aços AISI 52100 e AISI 1034 respectivamente. As melhores previsões foram obtidas através da versão desviadora SWTD. Segundo a abordagem local, o modelo SWTD fez previsões menos conservativas que os demais, sendo estas com erros conservativos inferiores a 15 porcento. Segundo a abordagem não-local, os quatro modelos geram previsões com erros inferiores a 15 porcento se os mesmos forem calibrados de forma independente. É importante mencionar que as melhores previsões obtidas através da abordagem não-local (que requer a calibração de um parâmetro adicional) são similares àquelas obtidas por SWTD local. Especificamente para os testes de fadiga por fretting com carga remota, o modelo SWTD não-local gerou previsões com 100 porcento de acertos (13 testes). Segundo a abordagem não-local, foi verificado que, para o conjunto de testes estudado, a precisão nas previsões é mais influenciada pela dimensão característica (que deve ser calibrada) que pelo modelo de fadiga multiaxial ou método de integração considerado. Por fim, verificou-se que previsões não conservativas podem ser obtidas ao se considerar uma dimensão característica assumida como parâmetro material (sem calibração) tal qual a distância crítica proposta por Taylor, definida em termos do método do ponto. / [en] This work presents and discusses four multiaxial fatigue damage models (SWTD, SWT, Findley and Crossland) applied to crack nucleation prediction under fretting fatigue loadings according to the local and the non-local approaches. The predictions have been compared with experimental data available in the literature for the cylinder on flat configuration made of 52100 and 1034 AISI steel s respectively. In particular, for the data set analyzed, SWTD parameter (deviatoric version of SWT model) predictions were less conservative than the ones made by the other three models and have better correlated the experimental results. Considering the local approach, the error observed on the SWTD predictions was lower than 15 percent. On the other hand, for the non-local approach, it has been demonstrated that the four multiaxial fatigue models can provide crack nucleation predictions with similar level of accuracy (error lower than 15 percent) if the characteristic dimension is calibrated model-by-model. Note that this same level of accuracy was obtained with SWTD parameter defined in terms of the local approach, which requires one less parameter to be calibrated (characteristic dimension). In particular, considering only the fretting fatigue tests with bulk stresses, SWTD non-local (area method) prediction accuracy was 100 percent (13 tests). For the material studied, it has also been concluded that the non-local predictions accuracy depends more on the calibrated characteristic dimension than on the multiaxial model or integration method considered. Moreover, it has been verified that non-conservative predictions can be made if the characteristic dimension is assumed as a material property as proposed by Taylor on the critical distances theory (point method).

Cette thèse est une contribution à la construction de l’Approche Locale de la rupture à l’échelle microscopique à l’aide de la modélisation d’agrégats polycristallins. Elle consiste à prendre en compte la variabilité spatiale de la microstructure du matériau. Pour ce faire, la modélisation micromécanique du matériau est réalisée par la simulation d’agrégats polycristallins par éléments finis. Les champs aléatoires de contrainte (principale maximale et de clivage) dans le matériau qui représentent la variabilité spatiale de la microstructure sont ensuite modélisés par un champ aléatoire gaussien stationnaire ergodique. Les propriétés de variabilité spatiale de ces champs sont identifiés par une méthode d’identification, e.g. méthode du périodogramme, méthode du variogramme, méthode du maximum de vraisemblance. Des réalisations synthétiques des champs de contraintes sont ensuite simulées par une méthode de simulation, e.g. méthode Karhunen-Loève discrète, méthode “Circulant Embedding”, méthode spectrale, sans nouveau calcul aux éléments finis. Enfin, le modèle d’Approche Locale de la rupture par simulation de champ de contrainte de clivage permettant d’y intégrer les réalisations simulées du champ est construit pour estimer la probabilité de rupture du matériau. / This thesis is a contribution to the construction of the Local Approach to fracture at the microscopic scale using polycrystalline aggregate modeling. It consists in taking into account the spatial variability of the microstructure of the material. To do this, the micromechanical modeling is carried out by finite element analysis of polycrystalline aggregates. The random stress fields (maximum principal et cleavage stress) in the material representing the spatial variability of the microstructure are then modeled by a stationary ergodic Gaussian random field. The properties of the spatial variability of these fields are identified by an identification method, e.g. periodogram method, variogram method, maximum likelihood method. The synthetic realizations of the stress fields are then simulated by a simulation method, e.g. discrete Karhunen-Loève method, circulant embedding method, spectral method, without additional finite element calculations. Finally, a Local Approach to fracture by simulation of the cleavage stress field using the simulated realizations is constructed to estimate the rupture probability of the material.