Feasibility of agent-based modelling of articular cartilage including a conceptual representation of its structure

Articular cartilage is a complex structure with an architecture in which fluid-swollen proteoglycans constrained within a 3D network of collagen fibrils. Because of the complexity of the cartilage structure, the relationship between its mechanical behaviours at the macroscale level and its components at the micro-scale level are not completely understood. The research objective in this thesis is to create a new model of articular cartilage that can be used to simulate and obtain insight into the micro-macro-interaction and mechanisms underlying its mechanical responses during physiological function. The new model of articular cartilage has two characteristics, namely: i) not use fibre-reinforced composite material idealization ii) Provide a framework for that it does probing the micro mechanism of the fluid-solid interaction underlying the deformation of articular cartilage using simple rules of repartition instead of constitutive / physical laws and intuitive curve-fitting. Even though there are various microstructural and mechanical behaviours that can be studied, the scope of this thesis is limited to osmotic pressure formation and distribution and their influence on cartilage fluid diffusion and percolation, which in turn governs the deformation of the compression-loaded tissue. The study can be divided into two stages. In the first stage, the distributions and concentrations of proteoglycans, collagen and water were investigated using histological protocols. Based on this, the structure of cartilage was conceptualised as microscopic osmotic units that consist of these constituents that were distributed according to histological results. These units were repeated three-dimensionally to form the structural model of articular cartilage. In the second stage, cellular automata were incorporated into the resulting matrix (lattice) to simulate the osmotic pressure of the fluid and the movement of water within and out of the matrix; following the osmotic pressure gradient in accordance with the chosen rule of repartition of the pressure. The outcome of this study is the new model of articular cartilage that can be used to simulate and study the micromechanical behaviours of cartilage under different conditions of health and loading. These behaviours are illuminated at the microscale level using the socalled neighbourhood rules developed in the thesis in accordance with the typical requirements of cellular automata modelling. Using these rules and relevant Boundary Conditions to simulate pressure distribution and related fluid motion produced significant results that provided the following insight into the relationships between osmotic pressure gradient and associated fluid micromovement, and the deformation of the matrix. For example, it could be concluded that: 1. It is possible to model articular cartilage with the agent-based model of cellular automata and the Margolus neighbourhood rule. 2. The concept of 3D inter connected osmotic units is a viable structural model for the extracellular matrix of articular cartilage. 3. Different rules of osmotic pressure advection lead to different patterns of deformation in the cartilage matrix, enabling an insight into how this micromechanism influences macromechanical deformation. 4. When features such as transition coefficient were changed, permeability (representing change) is altered due to the change in concentrations of collagen, proteoglycans (i.e. degenerative conditions), the deformation process is impacted. 5. The boundary conditions also influence the relationship between osmotic pressure gradient and fluid movement at the micro-scale level. The outcomes are important to cartilage research since we can use these to study the microscale damage in the cartilage matrix. From this, we are able to monitor related diseases and their progression leading to potential insight into drug-cartilage interaction for treatment. This innovative model is an incremental progress on attempts at creating further computational modelling approaches to cartilage research and other fluid-saturated tissues and material systems.

Performance of buried tunnels subjected to surface blast incorporating fluid structure interaction

This paper uses finite element techniques to investigate the performance of buried tunnels subjected to surface blasts incorporating fully coupled Fluid Structure Interaction and appropriate material models which simulate strain rate effects. Modelling techniques are first validated against existing experimental results and then used to treat the blast induced shock wave propagation and tunnel response in dry and saturated sands. Results show that the tunnel buried in saturated sand responds earlier than that in dry sand. Tunnel deformations decrease with distance from explosive in both sands, as expected. In the vicinity of the explosive, the tunnel buried in saturated sand suffered permanent deformation in both axial and circumferential directions, whereas the tunnel buried in dry sand recovered from most of the axial deformation. Overall, response of the tunnel in saturated sand is more severe for a given blast event and shows the detrimental effect of pore water on the blast response of buried tunnels. The validated modelling techniques developed in this paper can be used to investigate the blast response of tunnels buried in dry and saturated sands.

Deformation of a three-dimensional red blood cell in a stenosed micro-capillary

Red blood cells (RBCs) exhibit different types of motions and deformations when the blood flows through capillaries. Interestingly, due to the complex three-dimensional structure of the RBC membrane, RBCs show three-dimensional motions and deformations in the blood flow. These motions and deformations of the RBCs highly depend on the stiffness of the RBC membrane and on the geometrical parameters of the capillary through which blood flows. However, capillaries always do not have uniform cross sections and some capillaries have stenosed segments, where cross sectional area suddenly reduces. Further, some diseases can alter the stiffness of the RBC membrane drastically. In this study, the deformation behaviour of a single three-dimensional RBC is examined, when it moves through a stenosed capillary. A three-dimensional spring network is used to model the RBC membrane. The RBC’s inside and outside fluids are discretized into a finite number of mass points and treated by smoothed particle hydrodynamics (SPH) method. The capillary is considered as a rigid tube with a stenosed section. The deformation index, mean velocity and total energy of the RBC are analysed when it flows through the stenosed capillary. Further, motion and deformation of the RBCs with different membrane stiffness (KB) are compared when they flow through the stenosed segment of the capillary. The simulation results demonstrate the RBCs are subjected to a larger deformation when they move through the stenosed part of the capillary and the RBCs with lower KBvalues easily pass through the stenosed segment of the capillary. Further, RBCs having higher KBvalues have a lower mean velocity and it leads to slow down the overall blood flow rate

A tensor-based morphometry study of genetic influences on brain structure using a new fluid registration method

We incorporated a new Riemannian fluid registration algorithm into a general MRI analysis method called tensor-based morphometry to map the heritability of brain morphology in MR images from 23 monozygotic and 23 dizygotic twin pairs. All 92 3D scans were fluidly registered to a common template. Voxelwise Jacobian determinants were computed from the deformation fields to assess local volumetric differences across subjects. Heritability maps were computed from the intraclass correlations and their significance was assessed using voxelwise permutation tests. Lobar volume heritability was also studied using the ACE genetic model. The performance of this Riemannian algorithm was compared to a more standard fluid registration algorithm: 3D maps from both registration techniques displayed similar heritability patterns throughout the brain. Power improvements were quantified by comparing the cumulative distribution functions of the p-values generated from both competing methods. The Riemannian algorithm outperformed the standard fluid registration.

Best individual template selection from deformation tensor minimization

We study the influence of the choice of template in tensor-based morphometry. Using 3D brain MR images from 10 monozygotic twin pairs, we defined a tensor-based distance in the log-Euclidean framework [1] between each image pair in the study. Relative to this metric, twin pairs were found to be closer to each other on average than random pairings, consistent with evidence that brain structure is under strong genetic control. We also computed the intraclass correlation and associated permutation p-value at each voxel for the determinant of the Jacobian matrix of the transformation. The cumulative distribution function (cdf) of the p-values was found at each voxel for each of the templates and compared to the null distribution. Surprisingly, there was very little difference between CDFs of statistics computed from analyses using different templates. As the brain with least log-Euclidean deformation cost, the mean template defined here avoids the blurring caused by creating a synthetic image from a population, and when selected from a large population, avoids bias by being geometrically centered, in a metric that is sensitive enough to anatomical similarity that it can even detect genetic affinity among anatomies.

Analysis of effects of axial shortening of steel columns in frame structure

Steel columns in frame structure always carry heavy upcoming compressive forces. As a consequence, axial shortening becomes a common phenomenon in a multistoried steel structure. A 100 storied steel structure is analyzed in SAP2000 to study the magnitude overall effects of column shortening. It was found from the study that the maximum axial shortening occurs at the columns of top storey of the steel structure and at the columns of bottom storey, the axial deformation is negligible. The increasing rate of axial shortening is significant at the initial levels. However, at the upper levels, the amount of axial shortening in steel columns differs insignificantly. In the selected rigid frame structure, the axial shortening of adjacent steel columns is found to influence significantly the differential shortening of the structure. The consequent effect of differential shortening leads to develop excessive stress in the corner joints which ultimately hamper the normal behavior of the structural systems. The results are discussed elaborately in the paper.

The Crystal and Molecular Structure of 4-Aza-5-oxo-tricyclo[4.4.0.03,8]decane (Aza-twistanone)

Die kristalline Struktur von Aza-twistanon wurde durch eine Röntgenstruktur-analyse untersucht. Die Kristalle gehören zur monoklinen Raumgruppe P21/n mit den Zelldimensionen a = 6,662(6), b = 13,36(2), c = 8,606(9) Å, = 98,97(2)°, V = 757 Å3, Z = 4. Die Struktur wurde mit Direktmethoden gelöst und bis zu R = 0,035 verfeinert (mittlere (c) = 0,003 Å3).Die cis-Amidgruppe ist relativ stark deformiert und hat einen Torsionswinkel C -C -N-C von 14,5(4)° (Deformation aus der Ebene c = 5,0(5)° und N = 13,5(4,0)°). Die gegenüberliegende äthylenbrücke weist einen Torsionswinkel von 25,1(5)° auf. Die entsprechenden Winkel in Twistan betragen je 20°. Das tricyclische Gerüst von Aza-twistanon hat approximative.

Structure, function and intracellular dynamics of alphavirus replication complexes

Intracellular membrane alterations are hallmarks of positive-sense RNA (+RNA) virus replication. Strong evidence indicates that within these exotic compartments, viral replicase proteins engage in RNA genome replication and transcription. To date, fundamental questions such as the origin of altered membranes, mechanisms of membrane deformation and topological distribution and function of viral components, are still waiting for comprehensive answers. This study addressed some of the above mentioned questions for the membrane alterations induced during Semliki Forest virus (SFV) infection of mammalian cells. With the aid of electron and fluorescence microscopy coupled with radioactive labelling and immuno-cytochemistry techniques, our group and others showed that few hours after infection the four non structural proteins (nsP1-4) and newly synthesized RNAs of SFV colocalized in close proximity of small membrane invaginations, designated as spherules . These 50-70 nm structures were mainly detected in the perinuclear area, at the limiting membrane of modified endosomes and lysosomes, named CPV-I (cytopathic vacuoles type I). More rarely, spherules were also found at the plasma membrane (PM). In the first part of this study I present the first three-dimensional reconstruction of the CPV-I and the spherules, obtained by electron tomography after chemical or cryo-fixation. Different approaches for imaging these macromolecular assemblies to obtain better structure preservation and higher resolution are presented as unpublished data. This study provides insights into spherule organization and distribution of viral components. The results of this and other experiments presented in this thesis will challenge currently accepted models for virus replication complex formation and function. In a revisitation of our previous models, the second part of this work provides the first complete description of the biogenesis of the CPV-I. The results demonstrate that these virus-induced vacuoles, where hundreds of spherules accumulate at late stages during infection, represent the final phase of a journey initiated at the PM, which apparently serves as a platform for spherule formation. From the PM spherules were internalized by an endocytic event that required the activity of the class I PI3K, caveolin-1, cellular cholesterol and functional actin-myosin network. The resulting neutral endocytic carrier vesicle delivered the spherules to the membrane of pre-existing acidic endosomes via multiple fusion events. Microtubule based transport supported the vectorial transfer of these intermediates to the pericentriolar area where further fusions generated the CPV-I. A signal for spherule internalization was identified in one of the replicase proteins, nsP3. Infections of cells with viruses harbouring a deletion in a highly phosphorylated region of nsP3 did not result in the formation of CPV-Is. Instead, thousands of spherules remained at the PM throughout the infection cycle. Finally, the role of the replicase protein nsP2 during viral RNA replication and transcription was investigated. Three enzymatic activities, protease, NTPase and RNA-triphosphatase were studied with the aid of temperature sensitive mutants in vitro and, when possible, in vivo. The results highlighted the interplay of the different nsP2 functions during different steps of RNA replication and sub-genomic promoter regulation, and suggest that the protein could have different activities when participating in the replication complex or as a free enzyme.

Studies on the Cell Wall Structure and on the Mechanical Properties of Norway Spruce

The structure and the mechanical properties of wood of Norway spruce (Picea abies [L.] Karst.) were studied using small samples from Finland and Sweden. X-ray diffraction (XRD) was used to determine the orientation of cellulose microfibrils (microfibril angle, MFA), the dimensions of cellulose crystallites and the average shape of the cell cross-section. X-ray attenuation and x-ray fluorescence measurements were used to study the chemical composition and the trace element content. Tensile testing with in situ XRD was used to characterise the mechanical properties of wood and the deformation of crystalline cellulose within the wood cell walls. Cellulose crystallites were found to be 192 284 Å long and 28.9 33.4 Å wide in chemically untreated wood and they were longer and wider in mature wood than in juvenile wood. The MFA distribution of individual Norway spruce tracheids and larger samples was asymmetric. In individual cell walls, the mean MFA was 19 30 degrees, while the mode of the MFA distribution was 7 21 degrees. Both the mean MFA and the mode of the MFA distribution decreased as a function of the annual ring. Tangential cell walls exhibited smaller mean MFA and mode of the MFA distribution than radial cell walls. Maceration of wood material caused narrowing of the MFA distribution and removed contributions observed at around 90 degrees. In wood of both untreated and fertilised trees, the average shape of the cell cross-section changed from circular via ambiguous to rectangular as the cambial age increased. The average shape of the cell cross-section and the MFA distribution did not change as a result of fertilisation. The mass absorption coefficient for x-rays was higher in wood of fertilised trees than in that of untreated trees and wood of fertilised trees contained more of the elements S, Cl, and K, but a smaller amount of Mn. Cellulose crystallites were longer in wood of fertilised trees than in that of untreated trees. Kraft cooking caused widening and shortening of the cellulose crystallites. Tensile tests parallel to the cells showed that if the mean MFA is initially around 10 degrees or smaller, no systematic changes occur in the MFA distribution due to strain. The role of mean MFA in defining the tensile strength or the modulus of elasticity of wood was not as dominant as that reported earlier. Crystalline cellulose elongated much less than the entire samples. The Poisson ratio νca of crystalline cellulose in Norway spruce wood was shown to be largely dependent on the surroundings of crystalline cellulose in the cell wall, varying between -1.2 and 0.8. The Poisson ratio was negative in kraft cooked wood and positive in chemically untreated wood. In chemically untreated wood, νca was larger in mature wood and in latewood compared to juvenile wood and earlywood.

Structure of the hydrogen-bonded water molecule in crystals

A correlation of the structural data on IS hydrates obtained by x-ray diffraction, neutron diffraction, and proton magnetic resonance reveals that when a water molecule is hydrogen bonded into a crystal structure and the angle subtended at the donor water oxygen by the acceptor atoms deviates from the vapor H-O-H angle, bent hydrogen bonds are formed in preference to distortion of the H-O-H angle. Theoretical justification for this result is obtained from energy considerations by calculating the energy of formation of bent hydrogen bonds on the basis of the Lippincott-Schroeder potential function model for the hydrogen bond and the energy of deformation of the H-O-H angle from spectroscopic force constants.

The effect of cooling rate on the structure and properties of closed-cell aluminium foams

The application of different cooling rates as a strategy to enhance the structure of aluminium foams is studied. The potential to influence the level of morphological defects and cell size non-uniformities is investigated. AlSi6Cu4 alloy was foamed through the powder compact route and then solidified, applying three different cooling rates. Foam development was monitored in situ by means of X-ray radioscopy while foaming inside a closed mould. The macro-structure of the foams was analysed in terms of cell size distribution as determined by X-ray tomography. Compression tests were conducted to assess the mechanical performance of the foams and measured properties were correlated with structural features of the foams. Moreover, possible changes in the ductile brittle nature of deformation with cooling rate were analysed by studying the initial stages of deformation. We observed improvements in the cell size distributions, reduction in microporosity and grain size at higher cooling rates, which in turn led to a notable enhancement in compressive strength. (C) 2010 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Deformation and structural densification in Al2O3-Y2O3 glass

In this article, we examine the unusual plastic deformation under uniaxial compression of an Al2O3-15 mol % Y2O3 (A15Y) glass synthesized by a wet chemical route At a low temperature of 650-725 degrees C, plastic deformation of this glass is largely non-viscous through shear instabilities In contrast deformation near the crystallization temperature (850 degrees C) occurs homogeneously with work hardening and with a monotonic increase in the true density of the glass by 10-12% accompanied by an increase in hardness (H) and elastic modulus (E) of up to 100% We hypothesize a phenomenon of molecular densification of the amorphous structure through a hierarchy of multiple phases, analogous to density- or entropy-driven amorphous to-amorphous phase transitions (polyamorphism) These results suggest that the present method of preparation and the unusual behavior can trigger a search for many more systems that display such behavior (C) 2010 Acta Materialia Inc Published by Elsevier Ltd All rights reserved

A Non-contact measurement technique to measure micro surface stress and obtain deformation profiles of the order of 1nm in Micro-cantilever based structures by single image Optical Diffraction method

A new method based on analysis of a single diffraction pattern is proposed to measure deflections in micro-cantilever (MC) based sensor probes, achieving typical deflection resolutions of 1nm and surface stress changes of 50 mu N/m. The proposed method employs a double MC structure where the deflection of one of the micro-cantilevers relative to the other due to surface stress changes results in a linear shift of intensity maxima of the Fraunhofer diffraction pattern of the transilluminated MC. Measurement of such shifts in the intensity maxima of a particular order along the length of the structure can be done to an accuracy of 0.01mm leading to the proposed sensitivity of deflection measurement in a typical microcantilever. This method can overcome the fundamental measurement sensitivity limit set by diffraction and pointing stability of laser beam in the widely used Optical Beam Deflection method (OBDM).

Hot deformation characteristics of INCONEL alloy MA 754 and development of a processing map

The characteristics of hot deformation of INCONEL alloy MA 754 have been studied processing maps obtained on the basis of flow stress data generated in compression in the temperature range 700-degrees-C to 1150-degrees-C and strain rate range 0.001 to 100 s-1. The map exhibited three domains. (1) A domain of dynamic recovery occurs in the temperature range 800-degrees-C to 1075-degrees-C and strain rate range 0.02 to 2 s-1, with a peak efficiency of 18 pct occurring at 950-degrees-C and 0.1 s-1. Transmission electron microscope (TEM) micrographs revealed stable subgrain structure in this domain with the subgrain size increasing exponentially with an increase in temperature. (2) A domain exhibiting grain boundary cracking occurs at temperatures lower than 800-degrees-C and strain rates lower than 0.01 s-1. (3) A domain exhibiting intense grain boundary cavitation occurs at temperatures higher than 1075-degrees-C. The material did not exhibit a dynamic recrystallization (DRX) domain, unlike other superalloys. At strain rates higher than about 1 s-1, the material exhibits flow instabilities manifesting as kinking of the elongated grains and adiabatic shear bands. The material may be safely worked in the domain of dynamic recovery but can only be statically recrystallized.

An assessment of creep deformation and fracture behavior of 2.25Cr-1Mo similar and dissimilar weld joints

The evaluation of the creep deformation and fracture behavior of a 2.25Cr-1Mo steel base metal, a 2.25Cr-1Mo/2.25Cr-1Mo similar weld joint, and a 2.25Cr-1Mo/Alloy 800 dissimilar weld joint at 823 K over a stress range of 90 to WO MPa has been carried out. The specimens for creep testing were taken from single-V weld pads fabricated by a shielded metal arc-welding process using 2.25Cr-1Mo steel (for similar-joint) and INCONEL 182 (for dissimilar-joint) electrodes. The weld pads were subsequently given a postweld hear treatment (PWHT) of 973 K for I hour. The microstructure and microhardness of the weld joints were evaluated in the as-welded, postweld heat-treated, and creep-tested conditions. The heat-affected zone (HAZ) of similar weld joint consisted of bainite in the coarse-prior-austenitic-grain (CPAG) region near the fusion line, followed by bainite in the fine-prior-austenitic-grain (FPAG) and intercritical regions merging with the unaffected base metal. In addition to the HAZ structures in the 2.25Cr-1Mo steel, the dissimilar weld joint displayed a definite INCONEL/2.25Cr-1Mo weld interface structure present either as a sharp line or as a diffuse region. A hardness trough was observed in the intercritical region of the HAZ in both weld joints, while a maxima in hardness was seen at the weld interface of the dissimilar weld joint. Both weld joints exhibited significantly lower rupture lives compared to the 2.25Cr-1Mo base metal. The dissimilar weld joint exhibited poor rupture life compared to the similar weld joint, at applied stresses lower than 130 MPa. In both weld joints, the strain distribution across the specimen gage length during creep testing varied significantly. During creep testing, localization of deformation occurred in the intercritical HAZ. In the similar weld joint, at all stress levels investigated, and in the dissimilar weld joint, at stresses greater than or equal to 150 MPa, the creep failure occulted in the intercritical HAZ. The fracture occurred by transgranular mode with a large number of dimples. At stresses below 150 MPa, the failure in the dissimilar weld joint occurred in the CPAG HAZ near to the weld interface. The failure occurred by extensive intergranular creep cavity formation.