A constitutive model for cavitation and cavity growth in rubber-like materials under arbitrary tri-axial loading

In this paper, we attempted to construct a constitutive model to deal with the phenomenon of cavitation and cavity growth in a rubber-like material subjected to an arbitrary tri-axial loading. To this end, we considered a spherical elementary representative volume in a general Rivlin's incompressible material containing a central spherical cavity. The kinematics proposed by [Hou, H.S., Abeyaratne, R., 1992. Cavitation in elastic and elastic-plastic solids. J. Mech. Phys. Solids 40, 571-722] was adopted in order to construct an approximate but optimal field. In order to establish a suitable constitutive law for this class of materials, we utilized the homogenisation technique that permits us to calculate the average strain energy density of the volume. The cavity growth was considered through a physically realistic failure criterion. Combination of the constitutive law and the failure criterion enables us to describe correctly the global behaviour and the damage evolution of the material under tri-axial loading. It was shown that the present models can efficiently reproduce different stress states, varying from uniaxial to tri-axial tensions, observed in experimentations. Comparison between predicted results and experimental data proves that the proposed model is accurate and physically reasonable. Another advantage is that the proposed model does not need special identification work, the initial Rivlin's law for the corresponding incompressible material is sufficient to form the new law for the compressible material resulted from cavitation procedure. (C) 2007 Elsevier Ltd. All rights reserved.

Catastrophic Rupture Induced Damage Coalescence in Heterogeneous Brittle Media

In heterogeneous brittle media, the evolution of damage is strongly influenced by the multiscale coupling effect. To better understand this effect, we perform a detailed investigation of the damage evolution, with particular attention focused on the catastrophe transition. We use an adaptive multiscale finite-element model (MFEM) to simulate the damage evolution and the catastrophic failure of heterogeneous brittle media. Both plane stress and plane strain cases are investigated for a heterogeneous medium whose initial shear strength follows the Weibull distribution. Damage is induced through the application of the Coulomb failure criterion to each element, and the element mesh is refined where the failure criterion is met. We found that as damage accumulates, there is a stronger and stronger nonlinear increase in stress and the stress redistribution distance. The coupling of the dynamic stress redistribution and the heterogeneity at different scales result in an inverse cascade of damage cluster size, which represents rapid coalescence of damage at the catastrophe transition.

Study on the Capacity Degradation of Bucket Foundation in Liquefied Sand Layer under Cyclic Loads

The capacity degradation of bucket foundation in liquefied sand layer under cyclic loads such as equivalent dynamic ice-induced loads is studied. A simplified numerical model of liquefied sand layer has been presented based on the dynamic centrifuge experiment results. The ice-induced dynamic loads are modeled as equivalent sine cyclic loads, the liquefaction degree in different position of sand layer and effects of main factors are investigated. Subsequently, the sand resistance is represented by uncoupled, non-linear sand springs which describe the sub-failure behavior of the local sand resistance as well as the peak capacity of bucket foundation under some failure criterion. The capacity of bucket foundation is determined in liquefied sand layer and the rule of capacity degradation is analyzed. The capacity degradation in liquefied sand layer is analyzed comparing with that in non-liquefied sand layer. The results show that the liquefaction degree is 0.9 at the top and is only 0.06 at the bottom of liquefied sand layer. The numerical results are agreement well with the centrifugal experimental results. The value of the degradation of bucket capacity is 12% in numerical simulating whereas it is 17% in centrifugal experiments.

Redundancy quantification in the safety assessment of existing concrete beam-and-slab bridges

In current practice the strength evaluation of a bridge system is typically based on firstly using elastic analysis to determine the distribution of load effects in the elements and then checking the ultimate section capacity of those elements. Ductility of the components in most bridge structures permits local yield and subsequent redistribution of the applied loads from the most heavily loaded elements. As a result a bridge can continue to carry additional loading even after one member has yielded, which has conventionally been adopted as the "failure criterion" in bridge strength evaluation. This means that a bridge with inherent redundancy has additional reserves of strength such that the failure of one element does not result in the failure of the complete system. For these bridges warning signs will show up and measures can be undertaken before the ultimate collapse is happening. This paper proposes a rational methodology for calculating the ultimate system strength and including in bridge evaluation the warning level due to redundancy. © 2004 Taylor & Francis Group, London.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

Fracturing pressure in clay

Hydraulic fracturing in clayey soils can be triggered by either tensile or shear failure. In this paper, the physical meanings of various equations to predict fracture initiation pressure proposed in the past are discussed using the cavity expansion theory. In particular, when fracturing pressure is plotted against initial confining pressure, published laboratory test results as well as analytical models show a linear relationship. When the slope is close to 2, fracture is initiated by tensile failure of the clay, whereas when the slope is close to 1, it is initiated by shear failure of the clay. In this study, the analytical models, validated only on laboratory test data to date, were applied to unique data from field grouting work in which extensive soil fracturing was carried out to improve the mechanical characteristics of the soft silty clay underlying a bell tower in Venice, Italy. By a careful assessment of initial confining pressure in the field, the variation in recorded injection pressures with confining pressure was examined. Results suggest that the fractures at this site were likely to be initiated by shear failure of the clay, and the values were similar to what was predicted by the model with the shear failure criterion. © 2013 American Society of Civil Engineers.

典型人类活动对边坡变形及稳定性的影响研究

China is a mountainous country in which geological hazards occurred frequently, especially in the east of China. Except the geology, topography and extreme climate, the large scale human activities have become a major factor to landslides. Typical human activities which induced landslides are fill, cut and underground mining. On the topic of the deformation mechanism and slope stability, taking three different man-made slopes as examples, deformation mechanism and slope stability were studied by several methods, such as field work, numerical modeling and monitor. The details are as following: (1) The numerical modeling approach advantages over other conventional methods such as limit methods, so the numerical modeling is the major tool in this thesis. So far, there is no uniform failure criterion for numerical simulation. The failure criterion were summarized and analyzed firstly, subsequently the appropriate criterion was determinated. (2) Taking 220kV Yanjin transformation substation fill slope as example, the deformable characteristic, unstable mode and laboratory tests were studied systematically. The results show: the slope deformation was probably caused by a combination effect of unfavorable topographic, geological and hydro geological conditions, and external loading due to filling. It was concluded that the creep deformation of the slope was triggered by external loading applied at the back of the slope. In order to define the calculating parameters, a set of consolidated drained (CD) tests, consolidated undrained (CU) tests, repeated direct shear tests and UCS tests were carried out. The stability of the slope before and after reinforcement was assessed using 3D numerical modeling and shear strength reduction technique. The numerical modeling results showed: the factor of safety (FOS) of the slope was 1.10 in the natural state, and reduced to 1.03 after fill, which was close to the critical state and it caused creeping slip or deformation under rainfall. The failure surface in the slope is in active shear failure, whereas tensile failure occurs at the slope crest. After the site was reinforced with piles, the FOS was 1.27. Therefore, the slope is stable after reinforcement measures were taken. (3) The cut slope stability is a complex problem. Taking the left cut slope of Xiangjiaba as example in this thesis, the deformation and slope stability were studied systematically by numerical modeling and monitor methods. The numerical results show: the displacement is gradually increasing along with the cutting, and the largest displacement is 27.5mm which located at the bench between the elevation 340 and 380. Some failure state units distribute near the undermining part and there is no linked failure state occurred from crest to bottom during cutting. After cutting, some failure units appeared at the ground surface between elevation 340 and 360. The increasing tense stress made the disturbed rock failed. The slope is stable after cutting by the monitor method, such as surface monitor, multipoint displacement meter, inclinometer and anchor cable tensometer. (4) The interaction between underground mining and slope stability is a common situation in mountainous. The slope deformation mechanism induced by underground mining may contributed significantly to slope destabilization. The Mabukan slope in xiangjiaba was analyzed to illustrate this. Failure mechanism and the slope stability were presented by numerical modeling and residual deformation monitor. The results show: the roof deformed to the free face and the floor uplift lightly to the free face. The subsidence basin is formed, but the subsidence and the horizontal movement is small, and there is no failure zone occurred. When the underground mining is going on, the roof deformation, subsidence and the horizontal movements begin increasing. The rock deformation near the free face is larger than the ground surface, and the interaction between these coal seams appeared. There are some tensile failures and shear failures occurred on the roof and floor, and a majority of failure is tensile failure. The roof deformation, subsidence and the horizontal movements increased obviously along with the underground mining. The failure characteristic is shear failure which means the tensile stress transformed to the compressive stress. So the underground mining will induced tensile stress first which lead to structure crack, subsequently the compressive stress appeared which result in slippage. The crest was subjected to horizontal tension which made the rock crack along with the joint. The long term residual deformation monitor demonstrates that the slope is stable after the underground mining stopped.

Predicting the ultimate load of a CFRP wingbox

The finite element method in conjunction with the Soutis-Fleck model is used to predict the residual strength after impact of a carbon-fibre reinforced plastic wingbox subjected to a cantilever type loading. The maximum stress failure criterion further validates the Soutis-Fleck model predictions. The Soutis-Fleck model predicts that the wingbox fails at a tip load of 99.2 kN, approximately 5.5% less than the experimental observation

Multiple cutout optimization in composite plates using evolutionary structural optimization

The optimization of cutouts in composite plates was investigated by implementing a procedure known as Evolutionary Structural Optimization. Perforations were introduced into a finite element mesh of the plate from which one or more cutouts of a predetermined size were evolved. In the examples presented, plates were rejected from around each evolving cutout based on a predefined rejection criterion. The Limiting ply within each plate element around the cutout was determined based on the Tsai-Hill failure criterion. Finite element plates with values below the product of the average Tsai-Hill number and a rejection criterion were subsequently removed. This process was iterated until a steady state was reached and the rejection criterion was then incremented by an evolutionary rate and the above steps repeated until the desired cutout area was achieved. Various plates with differing lay-up and loading parameters were investigated to demonstrate the generality and robustness of this optimization procedure.

Mechanical properties of structures 3D printed with cementitious powders

The three dimensional (3D) printing technology has undergone rapid development in the last few years and it is now possible to print engineering structures. This paper presents a study of the mechanical behavior of 3D printed structures using cementitious powder. Microscopic observation reveals that the 3D printed products have a layered orthotropic microstructure, in which each layer consists of parallel strips. Compression and flexural tests were conducted to determine the mechanical properties and failure characteristics of such materials. The test results confirmed that the 3D printed structures are laminated with apparent orthotropy. Based on the experimental results, a stress-strain relationship and a failure criterion based on the maximum stress criterion for orthotropic materials are proposed for the structures of 3D printed material. Finally, a finite element analysis was conducted for a 3D printed shell structure, which shows that the printing direction has a significant influence on the load bearing capacity of the structure.

Modélisation radiobiologique pour la planification des traitements en radiothérapie à partir de données d’imagerie spécifiques aux patients

Un modèle de croissance et de réponse à la radiothérapie pour le glioblastome multiforme (GBM) basé le formalisme du modèle de prolifération-invasion (PI) et du modèle linéaire-quadratique a été développé et implémenté. La géométrie spécifique au patient est considérée en modélisant, d'une part, les voies d'invasion possibles des GBM avec l'imagerie du tenseur de diffusion (DTI) et, d'autre part, les barrières à la propagation à partir des images anatomiques disponibles. La distribution de dose réelle reçue par un patient donné est appliquée telle quelle dans les simulations, en respectant l'horaire de traitement. Les paramètres libres du modèle (taux de prolifération, coefficient de diffusion, paramètres radiobiologiques) sont choisis aléatoirement à partir de distributions de valeurs plausibles. Un total de 400 ensembles de valeurs pour les paramètres libres sont ainsi choisis pour tous les patients, et une simulation de la croissance et de la réponse au traitement est effectuée pour chaque patient et chaque ensemble de paramètres. Un critère de récidive est appliqué sur les résultats de chaque simulation pour identifier un lieu probable de récidive (SPR). La superposition de tous les SPR obtenus pour un patient donné permet de définir la probabilité d'occurrence (OP). Il est démontré qu'il existe des valeurs de OP élevées pour tous les patients, impliquant que les résultats du modèle PI ne sont pas très sensibles aux valeurs des paramètres utilisés. Il est également démontré comment le formalisme développé dans cet ouvrage pourrait permettre de définir un volume cible personnalisé pour les traitements de radiothérapie du GBM.

Modeling, Estimation, and Control of Robot-Soil Interactions

This thesis presents the development of hardware, theory, and experimental methods to enable a robotic manipulator arm to interact with soils and estimate soil properties from interaction forces. Unlike the majority of robotic systems interacting with soil, our objective is parameter estimation, not excavation. To this end, we design our manipulator with a flat plate for easy modeling of interactions. By using a flat plate, we take advantage of the wealth of research on the similar problem of earth pressure on retaining walls. There are a number of existing earth pressure models. These models typically provide estimates of force which are in uncertain relation to the true force. A recent technique, known as numerical limit analysis, provides upper and lower bounds on the true force. Predictions from the numerical limit analysis technique are shown to be in good agreement with other accepted models. Experimental methods for plate insertion, soil-tool interface friction estimation, and control of applied forces on the soil are presented. In addition, a novel graphical technique for inverting the soil models is developed, which is an improvement over standard nonlinear optimization. This graphical technique utilizes the uncertainties associated with each set of force measurements to obtain all possible parameters which could have produced the measured forces. The system is tested on three cohesionless soils, two in a loose state and one in a loose and dense state. The results are compared with friction angles obtained from direct shear tests. The results highlight a number of key points. Common assumptions are made in soil modeling. Most notably, the Mohr-Coulomb failure law and perfectly plastic behavior. In the direct shear tests, a marked dependence of friction angle on the normal stress at low stresses is found. This has ramifications for any study of friction done at low stresses. In addition, gradual failures are often observed for vertical tools and tools inclined away from the direction of motion. After accounting for the change in friction angle at low stresses, the results show good agreement with the direct shear values.

Extensional rheology and stability of gas cell walls in bread doughs at elevated temperatures in relation to breadmaking performance

The rheological properties of gas cell walls in bread doughs are considered to be important in relation to their stability and gas retention during proof and baking. Large deformation rheological properties of gas cell walls were measured using biaxial extension for a number of doughs of varying breadmaking quality at constant strain rate and elevated temperatures of 25-60degreesC. Strain hardening and failure strain of cell walls both decreased with temperature, with cell walls in good breadmaking doughs remaining stable and retaining their strain hardening properties at higher temperatures (60degreesC), while the cell walls of poor breadmaking doughs became unstable at lower temperatures (45-50degreesC) and had lower strain hardening. Strain hardening measured at 50degreesC gave good correlations with baking volume, with the best correlations achieved between rheological measurements and baking tests that used similar mixing conditions. As predicted by the considered failure criterion, a strain hardening value of I defines a region below which gas cell walls become unstable, and discriminates well between the baking quality of a range of commercial flour blends of varying quality. This indicates that the stability of gas cell walls during baking is strongly related to strain hardening properties, and that extensional rheological measurements can be used as indicators of baking quality.

Rheology and the breadmaking process

The applications of rheology to the main processes encountered during breadmaking (mixing, sheeting, fermentation and baking) are reviewed. The most commonly used rheological test methods and their relationships to product functionality are reviewed. It is shown that the most commonly used method for rheological testing of doughs, shear oscillation dynamic rheology, is generally used under deformation conditions inappropriate for breadmaking and shows little relationship with end-use performance. The frequency range used in conventional shear oscillation tests is limited to the plateau region, which is insensitive to changes in the HMW glutenin polymers thought to be responsible for variations in baking quality. The appropriate deformation conditions can be accessed either by long-time creep or relaxation measurements, or by large deformation extensional measurements at low strain rates and elevated temperatures. Molecular size and structure of the gluten polymers that make up the major structural components of wheat are related to their rheological properties via modern polymer rheology concepts. Interactions between polymer chain entanglements and branching are seen to be the key mechanisms determining the rheology of HMW polymers. Recent work confirms the observation that the dynamic shear plateau modulus is essentially independent of variations in MW of glutens amongst wheat varieties of varying baking performance and also that it is not the size of the soluble glutenin polymers, but the secondary structural and rheological properties of the insoluble polymer fraction that are mainly responsible for variations in baking performance. Extensional strain hardening has been shown to be a sensitive indicator of entanglements and long-chain branching in HMW polymers, and is well related to baking performance of bread doughs. The Considere failure criterion for instability in extension of polymers defines a region below which bubble walls become unstable, and predicts that when strain hardening falls below a value of around 1, bubble walls are no longer stable and coalesce rapidly, resulting in loss of gas retention and lower volume and texture. Strain hardening in doughs has been shown to reach this value at increasingly higher temperatures for better breadmaking varieties and is directly related to bubble stability and baking performance. (C) 2003 Elsevier Ltd. All rights reserved.

The physics of baking: rheological and polymer molecular structure-function relationships in breadmaking

Molecular size and structure of the gluten polymers that make up the major structural components of wheat are related to their rheological properties via modem polymer rheology concepts. Interactions between polymer chain entanglements and branching are seen to be the key mechanisms determining the rheology of HMW polymers. Recent work confirms the observation that dynamic shear plateau modulus is essentially independent of variations in MW amongst wheat varieties of varying baking performance and is not related to variations in baking performance, and that it is not the size of the soluble glutenin polymers, but the structural and rheological properties of the insoluble polymer fraction that are mainly responsible for variations in baking performance. The rheological properties of gas cell walls in bread doughs are considered to be important in relation to their stability and gas retention during proof and baking, in particular their extensional strain hardening properties. Large deformation rheological properties of gas cell walls were measured using biaxial extension for a number of doughs of varying breadmaking quality at constant strain rate and elevated temperatures in the range 25-60 degrees C. Strain hardening and failure strain of cell walls were both seen to decrease with temperature, with cell walls in good breadmaking doughs remaining stable and retaining their strain hardening properties to higher temperatures (60 degrees C), whilst the cell walls of poor breadmaking doughs became unstable at lower temperatures (45-50 degrees C) and had lower strain hardening. Strain hardening measured at 50 degrees C gave good correlations with baking volume, with the best correlations achieved between those rheological measurements and baking tests which used similar mixing conditions. As predicted by the Considere failure criterion, a strain hardening value of I defines a region below which gas cell walls become unstable, and discriminates well between the baking quality of a range of commercial flour blends of varying quality. This indicates that the stability of gas cell walls during baking is strongly related to their strain hardening properties, and that extensional rheological measurements can be used as predictors of baking quality. (C) 2004 Elsevier B.V. All rights reserved.