Modelling of some biological materials using continuum mechanics

Continuum mechanics provides a mathematical framework for modelling the physical stresses experienced by a material. Recent studies show that physical stresses play an important role in a wide variety of biological processes, including dermal wound healing, soft tissue growth and morphogenesis. Thus, continuum mechanics is a useful mathematical tool for modelling a range of biological phenomena. Unfortunately, classical continuum mechanics is of limited use in biomechanical problems. As cells refashion the �bres that make up a soft tissue, they sometimes alter the tissue's fundamental mechanical structure. Advanced mathematical techniques are needed in order to accurately describe this sort of biological `plasticity'. A number of such techniques have been proposed by previous researchers. However, models that incorporate biological plasticity tend to be very complicated. Furthermore, these models are often di�cult to apply and/or interpret, making them of limited practical use. One alternative approach is to ignore biological plasticity and use classical continuum mechanics. For example, most mechanochemical models of dermal wound healing assume that the skin behaves as a linear viscoelastic solid. Our analysis indicates that this assumption leads to physically unrealistic results. In this thesis we present a novel and practical approach to modelling biological plasticity. Our principal aim is to combine the simplicity of classical linear models with the sophistication of plasticity theory. To achieve this, we perform a careful mathematical analysis of the concept of a `zero stress state'. This leads us to a formal de�nition of strain that is appropriate for materials that undergo internal remodelling. Next, we consider the evolution of the zero stress state over time. We develop a novel theory of `morphoelasticity' that can be used to describe how the zero stress state changes in response to growth and remodelling. Importantly, our work yields an intuitive and internally consistent way of modelling anisotropic growth. Furthermore, we are able to use our theory of morphoelasticity to develop evolution equations for elastic strain. We also present some applications of our theory. For example, we show that morphoelasticity can be used to obtain a constitutive law for a Maxwell viscoelastic uid that is valid at large deformation gradients. Similarly, we analyse a morphoelastic model of the stress-dependent growth of a tumour spheroid. This work leads to the prediction that a tumour spheroid will always be in a state of radial compression and circumferential tension. Finally, we conclude by presenting a novel mechanochemical model of dermal wound healing that takes into account the plasticity of the healing skin.

A neutron diffraction investigation of residual stresses in rail ends after severe deformation of rail surfaces

Insulated rail joints (IRJs) are a primary component of the rail track safety and signalling systems. Rails are supported by two fishplates which are fastened by bolts and nuts and, with the support of sleepers and track ballast, form an integrated assembly. IRJ failure can result from progressive defects, the propagation of which is influenced by residual stresses in the rail. Residual stresses change significantly during service due to the complex deformation and damage effects associated with wheel rolling, sliding and impact. IRJ failures can occur when metal flows over the insulated rail gap (typically 6-8 mm width), breaks the electrically isolated section of track and results in malfunction of the track signalling system. In this investigation, residual stress measurements were obtained from rail-ends which had undergone controlled amounts of surface plastic deformation using a full scale wheel-on-track simulation test rig. Results were compared with those obtained from similar investigations performed on rail ends associated with ex-service IRJs. Residual stresses were measured by neutron diffraction at the Australian Nuclear Science and Technology Organisation (ANSTO). Measurements with constant gauge volume 3x3x3 mm3 were carried in the central vertical plane on 5mm thick sliced rail samples cut by an electric discharge machine (EDM). Stress evolution at the rail ends was found to exhibit characteristics similar to those of the ex-service rails, with a compressive zone of 5mm deep that is counterbalanced by a tension zone beneath, extending to a depth of around 15mm. However, in contrast to the ex-service rails, the type of stress distribution in the test-rig deformed samples was apparently different due to the localization of load under the particular test conditions. In the latter, in contrast with clear stress evolution, there was no obvious evolution of d0. Since d0 reflects rather long-term accumulation of crystal lattice damage and microstructural changes due to service load, the loading history of the test rig samples has not reached the same level as the ex-service rails. It is concluded that the wheel-on-rail simulation rig provides the potential capability for testing the wheel-rail rolling contact conditions in rails, rail ends and insulated rail joints.

Residual stresses in rail-ends from the in-service insulated rail joints using neutron diffraction

Insulated rail joints (IRJs) are an integral part of the rail track signaling system and pose significant maintenance and replacement costs due to their low and fluctuating service lives. Failure occurs mainly in rail head region, bolt- holes of fishplates and web-holes of the rails. Propagation of cracks is influenced by the evolution of internal residual stresses in rails during rail manufacturing (hot-rolling, roller-straightening, and head-hardening process), and during service, particularly in heavy rail haul freight systems where loads are high. In this investigation, rail head accumulated residual stresses were analysed using neutron diffraction at the Australian Nuclear Science and Technology Organisation (ANSTO). Two ex-service two head-hardened rail joints damaged under different loading were examined and results were compared with those obtained from an unused rail joint reference sample in order to differentiate the stresses developed during rail manufacturing and stresses accumulated during rail service. Neutron diffraction analyses were carried out on the samples in longitudinal, transverse and vertical directions, and on 5mm thick sliceed samples cut by Electric Discharge Machining (EDM). For the rail joints from the service line, irrespective of loading conditions and in-service times, results revealed similar depth profiles of stress distribution. Evolution of residual stress fields in rails due to service was also accompanied by evidence of larger material flow based on reflected light and scanning electron microscopy studies. Stress evolution in the vicinity of rail ends was characterised by a compressive layer, approximately 5 mm deep, and a tension zone located approximately 5- 15mm below the surfaces. A significant variation of d0 with depth near the top surface was detected and was attributed to decarburization in the top layer induced by cold work. Stress distributions observed in longitudinal slices of the two different deformed rail samples were found to be similar. For the undeformed rail, the stress distributions obtained could be attributed to variations associated with thermo-mechanical history of the rail.

Influence of surface stress on elastic constants of nanohoneycombs

Surface effect on the four independent elastic constants of nanohoneycombs is investigated in this paper. The axial deformation of the horizontal cell wall is included, comparing to the Gibson's method, and the contributions of the two components of surface stress (i.e. surface residual stress and surface elasticity) are discussed. The result shows that the regular hexagonal honeycomb is not isotropic but orthotropic. An increase in the cell-wall thickness t leads to an increase in the discrepancy of the Young's moduli in both directions. Furthermore, the surface residual stress dominates the surface effect on the elastic constants when t < 15 nm (or the relative density <0.17), which is in contrast to that the surface elasticity does when t > 15 nm (or the relative density > 0.17) for metal Al. The present structure and theory may be useful in the design of future nanodevices.

Effect of carbon equivalent and inoculation on residual stresses in grey iron castings

It is virtually impossible to produce castings free from internal stresses using conventional methods of founding. Castings with appreciable stresses distort during storage, transportation, machining and service. Though composition and melt treatment are known to affect the magnitude of residual stress in castings, the data on the effect of carbon equivalent and inoculation on the magnitude of residual stress in castings are limited. In the present investigation, an attempt is made to study (i) the effect of carbon equivalent on residual stress in cast iron castings, and (ii) the effect of inoculants such as calcium silicide and ferrosilicon on residual stress in iron castings in the carbon equivalent range 3.0–4.0%. The results of the investigation indicate the following: (i) the residual strains decrease linearly with increase in carbon equivalent in the uninoculated and inoculated irons; (ii) the tensile residual stresses decrease linearly with increase in carbon equivalent value of the uninoculated, calcium silicide-inoculated and ferrosilicon-inoculated cast iron castings; (iii) the ratio of UTS to residual stress increased on inoculating the grid castings. This increase is higher for calcium silicide-inoculated grids than for ferrosilicon-inoculated grid castings. This implies that from the residual stress point of view, inoculation of the iron with calcium silicide is beneficial.

Effect of residual stresses and metallographic stability on the over all performance of integral diaphragm material

The integral diaphragm pressure transducer consists of a diaphragm machined from precipitation hardened martensitic (APX4) steel. Its performance is quite significant as it depends upon various factors such as mechanical properties including induced residual stress levels, metallurgical and physical parameters due to different stages of processing involved. Hence, the measurement and analysis of residual stress becomes very important from the point of in-service assessment of a component. In the present work, the stress measurements have been done using the X-ray diffraction (XRD) technique, which is a non-destructive test (NDT). This method is more reliable and widely used compared to the other NDT techniques. The metallurgical aspects have been studied by adopting the conventional metallographic practices including examination of microstructure using light microscope. The dimensional measurements have been carried out using dimensional gauge. The results of the present investigation reveals that the diaphragm material after undergoing series of realization processes has yielded good amount of retained austenite in it. Also, the presence of higher compressive stresses induced in the transducer results in non-linearity, zero shift and dimensional instability. The problem of higher retained austenite content and higher compressive stress have been overcome by adopting a new realization process involving machining and cold and hot stabilization soak which has brought down the retained austenite content to about 5–6% and acceptable level of compressive stress in the range −100 to −150 MPa with fine tempered martensitic phase structure and good dimensional stability. The new realization process seems to be quite effective in terms of controlling retained austenite content, residual stress, metallurgical phase as well as dimensional stability and this may result in minimum zero shift of the diaphragm system.

Stress Evolution With Temperature in Titanium-Rich NiTi Shape Memory Alloy Films

The effect of deposition temperature on residual stress evolution with temperature in Ti-rich NiTi films deposited on silicon substrates was studied. Ti-rich NiTi films were deposited on 3? Si (100) substrates by DC magnetron sputtering at three deposition temperatures (300, 350 and 400 degrees C) with subsequent annealing in vacuum at their respective deposition temperatures for 4 h. The initial value of residual stress was found to be the highest for the film deposited and annealed at 400 degrees C and the lowest for the film deposited and annealed at 300 degrees C. All the three films were found to be amorphous in the as-deposited and annealed conditions. The nature of the stress response with temperature on heating in the first cycle (room temperature to 450 degrees C) was similar for all three films although the spike in tensile stress, which occurs at similar to 330 degrees C, was significantly higher in the film deposited and annealed at 300 degrees C. All the films were also found to undergo partial crystallisation on heating up to 450 degrees C and this resulted in decrease in the stress values around 5560 degrees C in the cooling cycle. The stress response with temperature in the second thermal cycle (room temperature to 450 degrees C and back), which is reflective of the intrinsic film behaviour, was found to be similar in all cases and the elastic modulus determined from the stress response was also more or less identical. The three deposition temperatures were also not found to have a significant effect on the transformation characteristics of these films such as transformation start and finish temperatures, recovery stress and hysteresis.

Effect of Stress and Interface defects on Photo Luminescence of Si nano-crystals embedded in SiO2

Effect of stress and interface defects on photo luminescence property of a silicon nano-crystal (Si-nc) embedded in amorphous silicon dioxide (a-SiO2) are studied in this paper using a self-consistent quantum-continuum based modeling framework. Si-ncs or quantum dots show photoluminescence at room temperature. Whether its origin is due to Si-nc/a-SiO2 interface defects or quantum confinement of carriers in Si-nc is still an outstanding question. Earlier reports have shown that stresses greater than 12 GPa change the indirect energy band gap structure of bulk Si to a direct energy band gap structure. Such stresses are observed very often in nanostructures and these stresses influence the carrier confinement energy significantly. Hence, it is important to determine the effect of stress in addition to the structure of interface defects on photoluminescence property of Si-nc. In the present work, first a Si-nc embedded in a-SiO2 is constructed using molecular dynamics simulation framework considering the actual conditions they are grown so that the interface and residual stress in the structure evolves naturally during formation. We observe that the structure thus created has an interface of about 1 nm thick consisting of 41.95% of defective states mostly Sin+ (n = 0 to 3) coordination states. Further, both the Si-nc core and the embedding matrix are observed to be under a compressive strain. This residual strain field is applied in an effective mass k.p Hamiltonian formulation to determine the energy states of the carriers. The photo luminescence property computed based on the carrier confinement energy and interface energy states associated with defects will be analysed in details in the paper.

Stress states in individual Si particles of a cast Al-Si alloy: Micro-Raman analysis and microstructure based modeling

The stress states in Si particles of cast Al-Si based alloys depend on its morphology and the heat treatment given to the alloy. The Si particles fracture less on modification and fracture more in the heat treated condition. An attempt has been made in this work to study the effect of heat treatment and Si modification on the stress states of the particles. Such understanding will be valuable for predicting the ductility of the alloy. The stress states of Si particles are estimated by Raman technique and compared with the microstructure-based FEM simulations. Combination of Electron Back-Scattered Diffraction (EBSD) and frequency shift, polarized micro-Raman technique is applied to determine the stress states in Si particles with (111) orientations. Stress states are measured in the as-received state and under uniaxial compression. The residual stress, the stress in the elastic-plastic regime and the stress which causes fracture of the particles is estimated by Raman technique. FEM study demonstrates that the stress distribution is uniform in modified Si, whereas the unmodified Si shows higher and more complex stress states. The onset of plastic flow is observed at sharp corners of the particles and is followed by localization of strain between particles. Clustering of particles generates more inhomogeneous plastic strain in the matrix. Particle stress estimated by Raman technique is in agreement with FEM calculations. (C) 2014 Elsevier B.V. All rights reserved.

Thermal Residual Stresses in Particulate Composites and Its Toughening Effect

In this paper, an accurate formula for calculating the thermal residual stress field in a particle-reinforced composite are presented. Numerical examples are given to show r-variations of the thermal residual stresses. The increase in fracture toughness of matrix predicted by the thermal residual stress field is compared well with the experimentally measured increase.

Deflections and Curvatures of a Film–Substrate Structure with the Presence of Gradient Stress in MEMS Applications

The close form solutions of deflections and curvatures for a film–substrate composite structure with the presence of gradient stress are derived. With the definition of more precise kinematic assumption, the effect of axial loading due to residual gradient stress is incorporated in the governing equation. The curvature of film–substrate with the presence of gradient stress is shown to be nonuniform when the axial loading is nonzero. When the axial loading is zero, the curvature expressions of some structures derived in this paper recover the previous ones which assume the uniform curvature. Because residual gradient stress results in both moment and axial loading inside the film–substrate composite structure, measuring both the deflection and curvature is proposed as a safe way to uniquely determine the residual stress state inside a film–substrate composite structure with the presence of gradient stress.

Influence of thermal residual stresses on the composite macroscopic behavior

The influence of the thermal residual stress on the deformation behavior of a composite has been analyzed with a new micromechanical method. The method is based on secant moduli approximation and a new homogenized effective stress to characterize the plastic state of the matrix. It is found that the generated thermal residual stresses after cooling and their influence on the subsequent deformation behavior depends significantly on the aspect ratio of the inclusions. With prolate inclusions, the presence of thermal residual stresses generate a higher compressive hardening curves of the composite, but it is reversed with oblate inclusions. For particle reinforced composite, thermal residual stresses induce a tensile hardening curve higher than the compressive one and this is in agreement with experimental observations. (C) 1998 Elsevier Science Ltd.

Deflections and curvatures of a film-substrate structure with the presence of gradient stress in MEMS applications

The close form solutions of deflections and curvatures for a film-substrate composite structure with the presence of gradient stress are derived. With the definition of more precise kinematic assumption, the effect of axial loading due to residual gradient stress is incorporated in the governing equation. The curvature of film-substrate with the presence of gradient stress is shown to be nonuniform when the axial loading is nonzero. When the axial loading is zero, the curvature expressions of some structures derived in this paper recover the previous ones which assume the uniform curvature. Because residual gradient stress results in both moment and axial loading inside the film-substrate composite structure, measuring both the deflection and curvature is proposed as a safe way to uniquely determine the residual stress state inside a film-substrate composite structure with the presence of gradient stress.

Stress analysis of ZrO2/SiO2 multilayers deposited on different substrates with different thickness periods

The effects of repeating thickness periods on stress are studied in ZrO2/SiO2 multilayers deposited by electron-beam evaporation on BK7 glass and fused-silica substrates. The results show that the residual stress is compressive and decreases with an increase of the periods of repeating thickness in the ZrO2/SiO2 multilayers. At the same time, the residual stress in multilayers deposited on BK7 glass is less than that of samples deposited on fused silica. The variation of the microstructure examined by x-ray diffraction shows that microscopic deformation does not correspond to macroscopic stress, which may be due to variation of the interface stress. (c) 2005 Optical Society of America.