The deformation behavior of unsaturated soil has been the subject of numerous experimental and theoretical investigations. However, this phenomenon is not fully understood. Problems, such as the adoption of the proper stress variables, reduction of suction inducing collapse, suction effect on soil stiffness, rate dependency and air trapped within the soil under rainfall infiltration still need additional studies. In the present studies, an elasto-viscoplastic model for unsaturated soil is used based on two stress variables: 1) the skeleton stress is adopted as the stress variable; 2) suction is incorporated into the constitutive model to describe the collapse behavior. In addition, to investigate the multiphase behavior of unsaturated soil, a three-phase coupled model has been proposed based on the Theory of Porous Media (TPM) and finite deformation theory. Van Genuchten type of equation is employed as a constitutive equation between the saturation and the suction. Three-dimensional multiphase simulations are carried out to reproduce the behavior of unsaturated soil during monotonic loading triaxial tests under drained and undrained conditions for water and air. Compared with experimental results and the simulated results, it is seen that the proposed formulation is very suitable to describe the mechanical behaviors of unsaturated soil. Cyclic behavior of unsaturated soil has attracted much attention during the past few years. An elasto-viscoplastic cyclic model for saturated soil is extended for modeling of unsaturated soil. Based on finite deformation theory, three-dimensional multiphase analyses for unsaturated soil under cyclic loading are presented. The simulations are verified with cyclic triaxial tests on unsaturated silty clay under undrained for water and air conditions. It shows that the proposed multiphase formulation can be used to simulate the behaviors of unsaturated soil under cyclic loading. The high expansiveness of bentonite is another significant problem in unsaturated soil mechanics. In this research, an elasto-viscoplastic model for unsaturated expansive soil has been developed. An evolutional equation is adopted for describing the absorption of water into interlayer of clay platelets. In addition, the internal compaction effect caused by swelling of clay unit is expressed with the expansion of overconsolidation boundary surface and static yield surface. Based on the model, one-dimensional finite element analysis is conducted to study the development of swelling pressure. Compared with experimental results and simulated results, it is found that the proposed model can reproduce the effects of dry density and initial water content on swelling behavior. Using the proposed swelling model, two-dimensional swelling behaviors of the waste barrier are simulated. / Kyoto University (京都大学) / 0048 / 新制・課程博士 / 博士(工学) / 甲第13775号 / 工博第2879号 / 新制||工||1425(附属図書館) / 25991 / UT51-2008-C691 / 京都大学大学院工学研究科社会基盤工学専攻 / (主査)教授 岡 二三生, 教授 松岡 俊文, 准教授 木元 小百合 / 学位規則第4条第1項該当

学位授与大学：京都大学 ; 取得学位: 博士(工学) ; 学位授与年月日: 2007-09-25 ; 学位の種類: 新制・課程博士 ; 学位記番号: 工博第2841号 ; 請求記号: 新制/工/1418 ; 整理番号: 25526 / The consolidation of peat is complex due to the resultant large strain associated with the highly compressible nature of natural peat deposits and to the rapid changes in soil properties during the consolidation process. In addition, the consolidation process is further complicated by the occurrence of secondary compression which significantly contributes to the overall settlement of peaty soil. Therefore, it is necessary to take these properties into account in order to obtain better predictions from peat consolidation analyses. In the present study, the consolidation behavior of peaty clay found in Sri Lanka is extensively studied using a model based on the elasto-viscoplastic theory. The model can describe the prominent creep behavior of peaty soil as a continuous process. In addition, the model can accommodate the effect of structural degradation on the consolidation process. The analysis takes into account all the main features involved in the peat consolidation process, namely, finite strain, variable permeability, and the effect of secondary compression. Also, it considers the variable compressibility for stage-constructed embankments which exert high levels of pressure on the peaty subsoil. The constitutive equations used in the model and the procedure adapted to account for the above-mentioned features of the analysis are described. The constitutive model is based on Perzyna’s type viscoplastic theory and the Cambridge elasto-plastic theory combined with empirical evidence. In the finite element formulations, which are based on the finite deformation theory, an updated Lagrangian method is adopted. A description of the material parameters used in the model and the procedures applied to evaluate them, with standard laboratory and field tests, are explained. In addition, a performance of the model incorporating the original and the modified Cam-clay theory is evaluated by simulating triaxial test results. A comparison shows that with the present definition of the parameters, the original model yields more representative results than the model based on the modified Cam-clay theory. Initially, the capability of the constitutive model to capture the consolidation behavior is verified using the consolidation model test data on peaty clay found in Sri Lanka. It is confirmed that the constitutive model is able to predict the observed creep characteristics and the effect of sample thickness on settlement predictions for the material under consideration. The performance of the model in predicting the consolidation behavior under field conditions is studied using field data on instrumented earth fill constructed on peaty clay. One-dimensional compression is assumed for the peaty clay due to the large plane area of the fill. Separate analyses are carried out by the model considering the infinitesimal strain theory, the finite strain theory, and the finite strain theory together with the effect of structural degradation in order to explore how these features describe the observed field behavior. Analyses reveal that it is necessary to consider finite deformation together with the effect of structural degradation in order to successfully simulate the resultant large strain and the stagnated pore water pressure observed in the field. The construction of road embankments over peat deposits is quite problematic, and thus, it is often done after first improving the properties of the peaty soil through the utilization of appropriate ground-improvement techniques. Understanding the field response of peaty clay during this improvement process is naturally of great importance. A constitutive model is applied to predict the field performance of embankments constructed on peaty clay using different ground-improvement techniques. The back analysis of embankments constructed with the preloading method indicates that the model can be successfully applied to predict both the deformation and the stability of structures constructed on peaty clays. The stability of the embankment during and after construction is verified by investigating the stress-strain characteristics of the subsoil. The model applications used to predict the consolidation behavior of embankments constructed by the preloading method, combined with other ground-improvement techniques, are then discussed. Embankments constructed with prefabricated vertical drains (PVDs) and sand compaction piles (SCPs) are considered, and finite element analyses are carried out in all cases by converting the actual three-dimensional conditions that exist around the drains into simplified two-dimensional plane strain conditions. The field behavior when PVDs are installed in the peaty clay is simulated using the equivalent vertical permeability for the PVD-improved subsoil. In the case of SCPs, a conversion scheme is used to transform the axisymmetric nature of sand columns into equivalent plane strain conditions. A comparison of the predicted results with the field observations shows a reasonable agreement. An analysis of the PVD-improved foundation indicates that the installation of PVDs not only accelerates the rate of consolidation, but influences the deformation pattern of the subsoil due to embankment loading. The analysis also shows that the use of PVDs can significantly increase embankment stability. The model prediction for the SCP-improved foundation reveals that the stiffness and the area replacement ratio used in the conversion scheme play vital roles in predicting the behavior of SCP-improved soft grounds. The observed improvements in the bearing capacity of the subsoil and in the stability of the embankment, brought about by the installation of SCPs, can be simulated by the model. / Kyoto University (京都大学) / 0048 / 新制・課程博士 / 博士(工学) / 甲第13370号 / 工博第2841号 / 新制||工||1418(附属図書館) / 25526 / UT51-2007-Q771 / 京都大学大学院工学研究科社会基盤工学専攻 / (主査)教授 岡 二三生, 教授 田村 武, 准教授 木元 小百合 / 学位規則第4条第1項該当

Cette thèse développe une méthodologie originale et non intrusive de construction de modèles de substitution applicable à des modèles physiques multiparamétriques.La méthodologie proposée permet d’approcher en temps réel, sur l’ensemble du domaine paramétrique, de multiples quantités d’intérêt hétérogènes issues de modèles physiques.Les modèles de substitution sont basés sur des représentations en train de tenseurs obtenues lors d'une phase hors ligne de calculs intensifs.L'idée essentielle de la phase d'apprentissage est de construire simultanément les approximations en se basant sur un nombre limité de résolutions du modèle physique lancées à la volée.L'exploration parcimonieuse du domaine paramétrique couplée au format compact de train de tenseurs permet de surmonter le fléau de la dimension.L'approche est particulièrement adaptée pour traiter des modèles présentant un nombre élevé de paramètres définis sur des domaines étendus.Les résultats numériques sur des lois élasto-viscoplastiques non linéaires montrent que des modèles de substitution compacts en mémoire qui approchent précisément les différentes variables mécaniques dépendantes du temps peuvent être obtenus à des coûts modérés.L'utilisation de tels modèles exploitables en temps réel permet la conception d'outils d'aide à la décision destinés aux experts métiers dans le cadre d'études paramétriques et visent à améliorer la procédure de calibration des lois matériaux. / This thesis presents a novel non-intrusive methodology to construct surrogate models of parametric physical models.The proposed methodology enables to approximate in real-time, over the entire parameter space, multiple heterogeneous quantities of interest derived from physical models.The surrogate models are based on tensor train representations built during an intensive offline computational stage.The fundamental idea of the learning stage is to construct simultaneously all tensor approximations based on a reduced number of solutions of the physical model obtained on the fly.The parsimonious exploration of the parameter space coupled with the compact tensor train representation allows to alleviate the curse of dimensionality.The approach accommodates particularly well to models involving many parameters defined over large domains.The numerical results on nonlinear elasto-viscoplastic laws show that compact surrogate models in terms of memory storage that accurately predict multiple time dependent mechanical variables can be obtained at a low computational cost.The real-time response provided by the surrogate model for any parameter value allows the implementation of decision-making tools that are particularly interesting for experts in the context of parametric studies and aim at improving the procedure of calibration of material laws.

Esta dissertação investiga numericamente a influência dos efeitos viscosos e elásticos em escoamentos de materiais viscoplásticos no interior de uma cavidade dirigida. O modelo mecânico empregado é constituído pelas equações de conservação de massa e pelo princípio da quantidade de movimento linear, para fluidos incompressíveis, acoplado à equação constitutiva. Esta equação modifica o modelo viscoelástico de Oldroyd-B de modo a acomodar que os tempos de relaxação e retardo do material, bem como sua viscosidade viscoplástica, dependam das mudanças de sua microestrutura. A aproximação numérica do modelo emprega o método multi-campos de Galerkin mínimos-quadrados em termos do tensor de tensão extra, do vetor velocidade e do campo de pressão. Os resultados objetivam a determinção do tamanho e localização das regiões aparentemente não-escoadas do material, bem como sua deformação elástica, intensidade de tensão, e a sua vorticidade no interior da cavidade. Os resultados claramente indicam que o padrão do escoamento é fortemente influenciado pela variação dos efeitos elásticos (variação do tempo de relaxação adimensional, θ0 * ), viscosos (variação do índice de power-law, n) e cinemáticos (variação da velocidade adimensional, U* , do escoamento) no interior da cavidade. / This dissertation investigated numerically the influence of viscous and elastic effects on flows of viscoplastic materials within a lid-driven cavity. The mechanical model used is made up of mass and momentum balance equations, coupled with the constitutive equation. This equation modifies the viscoelastic Oldroyd-B model to accommodate both relaxation and retardation times, and viscosity function, dependent on the microstructure changes. Numerical approximations of the model make use a three-field Galerkin least squares method in terms of the extra stress tensor, velocity vector and pressure field. Computations focus on the determination of the size and position of apparently unyielded regions as well as the elastic deformation, stress intensity, and the vorticity within of the cavity. Results clearly indicate that the flow pattern is strongly influenced by the elastic (variation of the dimensionless relaxation time, θ0 * ), viscous (variation of the power-law index, n) and kinematic (variation of the dimensionless flow velocity, U* ) effects within the cavity.

Esta dissertação investiga numericamente a influência dos efeitos viscosos e elásticos em escoamentos de materiais viscoplásticos no interior de uma cavidade dirigida. O modelo mecânico empregado é constituído pelas equações de conservação de massa e pelo princípio da quantidade de movimento linear, para fluidos incompressíveis, acoplado à equação constitutiva. Esta equação modifica o modelo viscoelástico de Oldroyd-B de modo a acomodar que os tempos de relaxação e retardo do material, bem como sua viscosidade viscoplástica, dependam das mudanças de sua microestrutura. A aproximação numérica do modelo emprega o método multi-campos de Galerkin mínimos-quadrados em termos do tensor de tensão extra, do vetor velocidade e do campo de pressão. Os resultados objetivam a determinção do tamanho e localização das regiões aparentemente não-escoadas do material, bem como sua deformação elástica, intensidade de tensão, e a sua vorticidade no interior da cavidade. Os resultados claramente indicam que o padrão do escoamento é fortemente influenciado pela variação dos efeitos elásticos (variação do tempo de relaxação adimensional, θ0 * ), viscosos (variação do índice de power-law, n) e cinemáticos (variação da velocidade adimensional, U* , do escoamento) no interior da cavidade. / This dissertation investigated numerically the influence of viscous and elastic effects on flows of viscoplastic materials within a lid-driven cavity. The mechanical model used is made up of mass and momentum balance equations, coupled with the constitutive equation. This equation modifies the viscoelastic Oldroyd-B model to accommodate both relaxation and retardation times, and viscosity function, dependent on the microstructure changes. Numerical approximations of the model make use a three-field Galerkin least squares method in terms of the extra stress tensor, velocity vector and pressure field. Computations focus on the determination of the size and position of apparently unyielded regions as well as the elastic deformation, stress intensity, and the vorticity within of the cavity. Results clearly indicate that the flow pattern is strongly influenced by the elastic (variation of the dimensionless relaxation time, θ0 * ), viscous (variation of the power-law index, n) and kinematic (variation of the dimensionless flow velocity, U* ) effects within the cavity.

De nombreuses applications industrielles mettent en jeu des fluides complexes qui possèdent souvent un seuil d'écoulement leur permettant de résister à des efforts finis sans s'écouler. Par ailleurs, ces fluides peuvent glisser aux parois lorsque les conditions interfaciales sont favorables. Toutes ces propriétés influencent leurs écoulements autour d'obstacles. Cette thèse se propose de comprendre ces écoulements dans le domaine où les vitesses d'écoulement sont telles que les effets inertiels peuvent être négligés devant les effets visqueux eux-mêmes faibles par rapport aux effets plastiques. Elle analyse l'influence de la vitesse et du glissement sur la force de traînée et les champs cinématiques générés par l'écoulement très lent et en régime permanent d'un fluide à seuil autour d'obstacles aux surfaces adhérentes ou glissantes. Les géométries considérées sont le disque, la sphère, le cône et la plaque plane. Le fluide utilisé a un comportement élasto-viscoplastique pouvant être décrit par les modèles de Herschel-Bulkley et de Hooke. Ce comportement a été caractérisé en volume et en présence de glissement par des mesures rhéométriques. Le nombre adimensionnel clé de l'étude est le nombre d'Oldroyd, ratio entre les effets plastiques et les effets visqueux, compris ici entre 10 et 200. Les mesures de forces de traînée ont montré qu'indépendamment de l'obstacle et des conditions interfaciales, le coefficient de traînée diminue avec le nombre d'Oldroyd et tend vers une valeur asymptotique. Cette valeur montre qu'au-delà d'un certain nombre d'Oldroyd, ce coefficient n'est plus gouverné par la vitesse mais dépend uniquement du seuil et de la surface caractéristique de l'obstacle. Elle permet de calculer un critère de stabilité pour lequel l'objet est maintenu en suspension. Les champs cinématiques déterminés par PIV ont permis de caractériser la forme et l'étendue des zones rigides et cisaillées. Les mesures de forces de traînée et de champs cinématiques ont permis de quantifier la contribution des contraintes normales et tangentielles dans la force de traînée totale. La présence de glissement aux parois de l'obstacle diminue significativement le coefficient de traînée et modifie la morphologie de l'écoulement en réduisant l'étendue des zones cisaillées. Une simulation numérique a été menée dans le cas de la plaque plane avec un modèle élasto-viscoplastique et un code à éléments finis avec points d'intégration Lagrangiens. / Many industrial processes include numerous complex fluids often presenting a yield stress. Those fluids can also slip when interfacial conditions are favorable. All these properties affect their flows around obstacles. This thesis aims to understand such flows in a domain where the flow velocities are so low that inertia effects can be neglected compared to viscous effects which are substantially low compared to plastic effects. It analyzes the influence of the velocity and the slip on the drag force and the kinematic fields of the creeping flow of a yield stress fluid around obstacles either with adhesive or slippery wall. The flow is analyzed in steady regime. The considered geometries are the disc, the sphere, the cone and the flat plate. The fluid used has an elasto-viscoplastic behavior which is modelled by the Herschel-Bulkley and Hooke models. This behavior has been characterized by rheometrical tests performed with adherence and slip conditions. The main non-dimensional number is the Oldroyd number, i.e. the ratio between plastic and viscous effects, which ranges from 10 to 200. The drag forces measurements have shown that regardless of the obstacle and the interfacial conditions, the drag coefficient decreases with the Oldroyd number before tending towards to an asymptotical value. This asymptotical value highlights that for high Oldroyd numbers the drag coefficient is no longer governed by the velocity but depends only on the yield stress and the characteristic section of the obstacle. A stability criterion for which the obstacle is held in suspension has been calculated from it. The kinematic fields determined by PIV have enabled to characterize the shape and the extent of the sheared and static regions. The drag forces and the kinematic fields measurements have enabled to quantify the contribution of the normal and tangential stresses in the total drag force. The wall slip significantly reduces the drag coefficient and also reduces the extent of the sheared zones. A numerical simulation has been performed with an elasto-viscoplastic model by means of a code using finite elements method with Lagrangian integration points in the case of an adhesive flat plane.

Esta dissertação investiga numericamente a influência dos efeitos viscosos e elásticos em escoamentos de materiais viscoplásticos no interior de uma cavidade dirigida. O modelo mecânico empregado é constituído pelas equações de conservação de massa e pelo princípio da quantidade de movimento linear, para fluidos incompressíveis, acoplado à equação constitutiva. Esta equação modifica o modelo viscoelástico de Oldroyd-B de modo a acomodar que os tempos de relaxação e retardo do material, bem como sua viscosidade viscoplástica, dependam das mudanças de sua microestrutura. A aproximação numérica do modelo emprega o método multi-campos de Galerkin mínimos-quadrados em termos do tensor de tensão extra, do vetor velocidade e do campo de pressão. Os resultados objetivam a determinção do tamanho e localização das regiões aparentemente não-escoadas do material, bem como sua deformação elástica, intensidade de tensão, e a sua vorticidade no interior da cavidade. Os resultados claramente indicam que o padrão do escoamento é fortemente influenciado pela variação dos efeitos elásticos (variação do tempo de relaxação adimensional, θ0 * ), viscosos (variação do índice de power-law, n) e cinemáticos (variação da velocidade adimensional, U* , do escoamento) no interior da cavidade. / This dissertation investigated numerically the influence of viscous and elastic effects on flows of viscoplastic materials within a lid-driven cavity. The mechanical model used is made up of mass and momentum balance equations, coupled with the constitutive equation. This equation modifies the viscoelastic Oldroyd-B model to accommodate both relaxation and retardation times, and viscosity function, dependent on the microstructure changes. Numerical approximations of the model make use a three-field Galerkin least squares method in terms of the extra stress tensor, velocity vector and pressure field. Computations focus on the determination of the size and position of apparently unyielded regions as well as the elastic deformation, stress intensity, and the vorticity within of the cavity. Results clearly indicate that the flow pattern is strongly influenced by the elastic (variation of the dimensionless relaxation time, θ0 * ), viscous (variation of the power-law index, n) and kinematic (variation of the dimensionless flow velocity, U* ) effects within the cavity.

No description available.

L'objectif de cette thèse est d'étudier expérimentalement les phénomènes de vieillissement dus à la diffusion des atomes en solution dans les alliages d’aluminium et les instabilités qui leur sont associées comme le phénomène Portevin-Le Châtelier ou les bandes de Lüders et de proposer une modélisation de ces phénomènes dans le cadre de la thermodynamique des processus irréversibles. Une étude expérimentale détaillée est alors entreprise sur les alliages d’aluminium AA5083-H116 et AA5182-O. Le comportement du premier présente l’effet PLC de façon prononcée et les deux types d’instabilités sont observées simultanément pour le second. La corrélation d'images numériques et la thermographie infrarouge sont essentiellement employées pour détecter et caractériser les aspects spatiotemporels des instabilités observées. La déformation non homogène due à l’apparition et la propagation de bandes de localisation est mise en évidence. Ces bandes de déformation sont visualisées, permettant à leurs diverses caractéristiques (vitesse, orientation, largeur, vitesse de déformation à l'intérieur des bandes, augmentation de la température à l'intérieur des bandes) d'être mesurées sur des éprouvettes lisses (plates, cylindriques ou prismatiques). Dans le cas des éprouvettes plates, l’effet des épaisseurs de l’éprouvette a aussi été examiné. Les caractéristiques des bandes sont aussi analysées sur d’autres géométries d’éprouvette (entaillées avec divers types d’entailles) à des vitesses de déformation différentes pour exhiber leur morphologie en présence de chargements multiaxiaux. Des simulations non linéaires et tridimensionnelles ont été effectuées en utilisant le modèle de McCormick pour montrer comment la prise en compte des phénomènes de vieillissement, même partielle, permet de décrire les hétérogénéités et le mouvement des bandes de déformation ainsi que de prévoir leurs différentes caractéristiques. Enfin, en se basant sur les mécanismes physiques à la base des phénomènes de vieillissement et en soulignant les limites du modèle de McCormick, un modèle élasto-viscoplastique prenant en considération les phénomènes de vieillissement est proposé dans le cadre de la thermodynamique des processus irréversibles. / The objective of this thesis is to study experimentally the ageing phenomena due to the atoms diffusion in solid solutions in aluminium alloys and associated instabilities as the Portevin-Le Châtelier phenomenon or the Lüders bands and to propose a model of these phenomena within the framework of thermodynamics of irreversible processes. Digital Image Correlation (DIC) and Digital Infrared Thermography (DIT) are essentially employed to capture and characterize the spatio-temporal aspects of these instabilities. Inhomogeneous deformation, due to the appearance and the propagation of various localization bands are observed. These deformation bands are visualized, allowing their various characteristics (velocity, orientation, width, strain rate inside the bands, temperature increase inside the bands) to be measured on smooth specimens (flat, cylindrical or prismatic). In the case of flat test, the effect of the specimen thickness was also examined. The band characteristics are also analyzed on other specimen geometries (notched with various kinds of notches) at different strain rates to exhibit their morphology in the presence of multiaxial loadings. Some nonlinear and three-dimensional simulations were carried out with McCormick model to show how the inclusions of ageing phenomena, even partial, makes possible to describe heterogeneities and the deformation bands but also reproducing qualitatively their various characteristics. Finally, based on the physical mechanisms of aging phenomena and underlying the limitations of the McCormick model, an elasto-viscoplastic model taking into account the aging phenomena is proposed in the framework of thermodynamics of irreversible processes.