Determination of strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes using atomic force microscopy and inverse finite element analysis

The aim of this paper is to determine the strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes, in vitro. Firstly, Atomic Force Microscopy (AFM) was used to obtain the force-indentation curves of these single cells at four different strain-rates. These results were then employed in inverse finite element analysis (FEA) using Modified Standard neo-Hookean Solid (MSnHS) idealization of these cells to determine their mechanical properties. In addition, a FEA model with a newly developed spring element was employed to accurately simulate AFM evaluation in this study. We report that both cytoskeleton (CSK) and intracellular fluid govern the strain-rate-dependent mechanical property of living cells whereas intracellular fluid plays a predominant role on fixed cells’ behavior. In addition, through the comparisons, it can be concluded that osteocytes are stiffer than chondrocytes at all strain-rates tested indicating that the cells could be the biomarker of their tissue origin. Finally, we report that MSnHS is able to capture the strain-rate-dependent mechanical behavior of osteocyte and chondrocyte for both living and fixed cells. Therefore, we concluded that the MSnHS is a good model for exploration of mechanical deformation responses of single osteocytes and chondrocytes. This study could open a new avenue for analysis of mechanical behavior of osteocytes and chondrocytes as well as other similar types of cells.

Examining the cost efficiency of Chinese hydroelectric companies using a finite mixture model

This paper evaluates the operational activities of Chinese hydroelectric power companies over the period 2000-2010 using a finite mixture model that controls for unobserved heterogeneity. In so doing, a stochastic frontier latent class model, which allows for the existence of different technologies, is adopted to estimate cost frontiers. This procedure not only enables us to identify different groups among the hydro-power companies analysed, but also permits the analysis of their cost efficiency. The main result is that three groups are identified in the sample, each equipped with different technologies, suggesting that distinct business strategies need to be adapted to the characteristics of China's hydro-power companies. Some managerial implications are developed. © 2012 Elsevier B.V.

Porohyperelastic finite element model for the kangaroo humeral head cartilage based on experimental study and the consolidation theory

Solid-extracellular fluid interaction is believed to play an important role in the strain-rate dependent mechanical behaviors of shoulder articular cartilages. It is believed that the kangaroo shoulder joint is anatomically and biomechanically similar to human shoulder joint and it is easy to get in Australia. Therefore, the kangaroo humeral head cartilage was used as the suitable tissue for the study in this paper. Indentation tests from quasi-static (10-4/sec) to moderately high strain-rate (10-2/sec) on kangaroo humeral head cartilage tissues were conduced to investigate the strain-rate dependent behaviors. A finite element (FE) model was then developed, in which cartilage was conceptualized as a porous solid matrix filled with incompressible fluids. In this model, the solid matrix was modeled as an isotropic hyperelastic material and the percolating fluid follows Darcy’s law. Using inverse FE procedure, the constitutive parameters related to stiffness, compressibility of the solid matrix and permeability were obtained from the experimental results. The effect of solid-extracellular fluid interaction and drag force (the resistance to fluid movement) on strain-rate dependent behavior was investigated by comparing the influence of constant, strain dependent and strain-rate dependent permeability on FE model prediction. The newly developed porohyperelastic cartilage model with the inclusion of strain-rate dependent permeability was found to be able to predict the strain-rate dependent behaviors of cartilages.

Soft tissue human thigh and buttock finite element model to simulate vehicle seat cushion indentation

Accurate modelling of automotive occupant posture is strongly related to the mechanical interaction between human body soft tissue and flexible seat components. This paper presents a finite-element study simulating the deflection of seat cushion foam and supportive seat structures, as well as human buttock and thigh soft tissue when seated. The thigh-buttock surface shell model was based on 95th percentile male subject scan data and made of two layers, covering thin to moderate thigh and buttock proportions. To replicate the effects of skin and fat, the neoprene rubber layer was modelled as a hyperelastic material with viscoelastic behaviour. The analytical seat model is based on a Ford production seat. The result of the finite-element indentation simulation is compared to a previous simulation of an indentation with a hard shell human model of equal geometry, and to the physical indentation result. We conclude that SAE composite buttock form and human-seat indentation of a suspended seat cushion can be validly simulated.

Evaluation of a patient-specific finite-element model to simulate conservative treatment in adolescent idiopathic scoliosis

Study design Retrospective validation study. Objectives To propose a method to evaluate, from a clinical standpoint, the ability of a finite-element model (FEM) of the trunk to simulate orthotic correction of spinal deformity and to apply it to validate a previously described FEM. Summary of background data Several FEMs of the scoliotic spine have been described in the literature. These models can prove useful in understanding the mechanisms of scoliosis progression and in optimizing its treatment, but their validation has often been lacking or incomplete. Methods Three-dimensional (3D) geometries of 10 patients before and during conservative treatment were reconstructed from biplanar radiographs. The effect of bracing was simulated by modeling displacements induced by the brace pads. Simulated clinical indices (Cobb angle, T1–T12 and T4–T12 kyphosis, L1–L5 lordosis, apical vertebral rotation, torsion, rib hump) and vertebral orientations and positions were compared to those measured in the patients' 3D geometries. Results Errors in clinical indices were of the same order of magnitude as the uncertainties due to 3D reconstruction; for instance, Cobb angle was simulated with a root mean square error of 5.7°, and rib hump error was 5.6°. Vertebral orientation was simulated with a root mean square error of 4.8° and vertebral position with an error of 2.5 mm. Conclusions The methodology proposed here allowed in-depth evaluation of subject-specific simulations, confirming that FEMs of the trunk have the potential to accurately simulate brace action. These promising results provide a basis for ongoing 3D model development, toward the design of more efficient orthoses.

A model for the interaction of trifluoroethanol with peptides and proteins

The structural stabilizing property of 2,2,2-trifluoroethanol (TFE) in peptides has been widely demonstrated, More recently, TFE has been shown to enhance secondary structure content in globular proteins, and to influence quaternary interactions in protein multimers. The molecular mechanisms by which TFE exerts its Influence on peptide and protein structures remain poorly understood. The present analysis integrates the known physical properties of TFE with a variety of experimental observations on the interaction of TFE with peptides and proteins and on the properties of fluorocarbons. Two features of TFE, namely the hydrophobicity of the trifluoromethyl group and the hydrogen bonding character (strong donor and poor acceptor), emerge as the most important factors for rationalising the observed effects of TFE. A model is proposed for TFE interaction with peptides which involves an initial replacement of the hydration shell by fluoroalcohol molecules, a process driven by apolar interactions and favourable entropy of dehydration. Subsequent bifurcated hydrogen-bond formation with peptide carbonyl groups, which leave intramolecular interactions unaffected, promotes secondary structure formation.

Strong-coupling expansion for the momentum distribution of the Bose-Hubbard model with benchmarking against exact numerical results

A strong-coupling expansion for the Green's functions, self-energies, and correlation functions of the Bose-Hubbard model is developed. We illustrate the general formalism, which includes all possible (normal-phase) inhomogeneous effects in the formalism, such as disorder or a trap potential, as well as effects of thermal excitations. The expansion is then employed to calculate the momentum distribution of the bosons in the Mott phase for an infinite homogeneous periodic system at zero temperature through third order in the hopping. By using scaling theory for the critical behavior at zero momentum and at the critical value of the hopping for the Mott insulator–to–superfluid transition along with a generalization of the random-phase-approximation-like form for the momentum distribution, we are able to extrapolate the series to infinite order and produce very accurate quantitative results for the momentum distribution in a simple functional form for one, two, and three dimensions. The accuracy is better in higher dimensions and is on the order of a few percent relative error everywhere except close to the critical value of the hopping divided by the on-site repulsion. In addition, we find simple phenomenological expressions for the Mott-phase lobes in two and three dimensions which are much more accurate than the truncated strong-coupling expansions and any other analytic approximation we are aware of. The strong-coupling expansions and scaling-theory results are benchmarked against numerically exact quantum Monte Carlo simulations in two and three dimensions and against density-matrix renormalization-group calculations in one dimension. These analytic expressions will be useful for quick comparison of experimental results to theory and in many cases can bypass the need for expensive numerical simulations.

On a finite element model for the analysis of through cracks in laminated anisotropic cylindrical shells

A finite element model for the analysis of laminated composite cylindrical shells with through cracks is presented. The analysis takes into account anisotropic elastic behaviour, bending-extensional coupling and transverse shear deformation effects. The proposed finite element model is based on the approach of dividing a cracked configuration into triangular shaped singular elements around the crack tip with adjoining quadrilateral shaped regular elements. The parabolic isoparametric cylindrical shell elements (both singular and regular) used in this model employ independent displacement and rotation interpolation in the shell middle surface. The numerical comparisons show the evidence to the conclusion that the proposed model will yield accurate stress intensity factors from a relatively coarse mesh. Through the analysis of a pressurised fibre composite cylindrical shell with an axial crack, the effect of material orthotropy on the crack tip stress intensity factors is shown to be quite significant.

Finitary Abstract Elementary Classes

The research in model theory has extended from the study of elementary classes to non-elementary classes, i.e. to classes which are not completely axiomatizable in elementary logic. The main theme has been the attempt to generalize tools from elementary stability theory to cover more applications arising in other branches of mathematics. In this doctoral thesis we introduce finitary abstract elementary classes, a non-elementary framework of model theory. These classes are a special case of abstract elementary classes (AEC), introduced by Saharon Shelah in the 1980's. We have collected a set of properties for classes of structures, which enable us to develop a 'geometric' approach to stability theory, including an independence calculus, in a very general framework. The thesis studies AEC's with amalgamation, joint embedding, arbitrarily large models, countable Löwenheim-Skolem number and finite character. The novel idea is the property of finite character, which enables the use of a notion of a weak type instead of the usual Galois type. Notions of simplicity, superstability, Lascar strong type, primary model and U-rank are inroduced for finitary classes. A categoricity transfer result is proved for simple, tame finitary classes: categoricity in any uncountable cardinal transfers upwards and to all cardinals above the Hanf number. Unlike the previous categoricity transfer results of equal generality the theorem does not assume the categoricity cardinal being a successor. The thesis consists of three independent papers. All three papers are joint work with Tapani Hyttinen.

Studies on the Effects of Property Taxation, Rent Control and Housing Allowances

Tämä väitöskirja koostuu asuntomarkkinoiden taloustieteellistä analyysia esittelevästä johdantoluvusta ja kolmesta tutkimuksesta, joissa analysoidaan asuntomarkkinoihin vaikuttavia politiikkatoimenpiteitä. Luvussa 2 tutkitaan Suomen kiinteistöverojärjestelmän vaikutusta asuntorakentamiseen. Vuonna 2001 tehtiin uudistus, jonka myötä kunnat voivat verottaa rakentamatonta asuintonttia korkeammalla veroasteella kuin rakennettua tonttia. Maanomistajan rakentamispäätöksen teoreettisen mallin mukaan rakentamattoman tontin korotettu kiinteistöveron pitäisi nopeuttaa rakentamista, mutta toisaalta myös rakentamiseen investoitu rahamäärä saattaa muuttua. Asuintonttien kiinteistöverojen yleinen taso ei vaikuta maanomistajan käyttäytymiseen, sillä tontin verotusarvo ei riipu rakentamispäätöksestä. Vain rakentamattoman ja rakennetun tontin veroasteiden erolla on merkitystä. Empiiriset tulokset ovat sopusoinnussa teorian kanssa. Tulosten mukaan prosenttiyksikön nousu rakentamattoman ja rakennetun tontin veroasteiden erossa lisää omakotialoitusten määrää viidellä prosentilla lyhyellä aikavälillä. Luvussa 3 analysoidaan vuokrasääntelystä vuokralaisille aiheutuvia hyötyjä ja haittoja. Vuokrasäännellyissä asunnoissa asuvat kotitaloudet hyötyvät vuokrasääntelystä alhaisen vuokran muodossa. Heille saattaa kuitenkin koitua myös haittaa siitä, että toiveita vastaavan asunnon löytäminen on vuokrasääntelytilanteessa vaikeaa, koska vapaille asunnoille on suuri määrä ottajia. Vapaarahoitteisen vuokra-asuntokannan vuokrien sääntely purettiin Suomessa asteittain vuosina 1992–1995. Tutkimuksen empiiriset tulokset viittaavat siihen, että vuokrasääntelyn aiheuttamista suurista eroista halutun ja todellisen asuntokulutuksen välillä koituvat hyvinvointitappiot kumosivat merkittävän osan matalien vuokrien hyödyistä vuokralaisille. Luvussa 4 tutkitaan Suomen asumistukijärjestelmän kannustinvaikutuksia. Asumistuen määrää rajoittavat asunnon pinta-alalle ja neliövuokralle asetetut ylärajat. Neliövuokrarajoite voidaan tulkita asumisen laatua rajoittavana tekijänä. Tutkimuksen teoreettisessa osassa osoitetaan, että asumistukijärjestelmä luo vahvat kannustimet muuttaa asuntoihin, joissa pinta-ala- ja laaturajoitteet purevat. Empiiristen tulosten mukaan asumistukeen oikeutetut kotitaloudet eivät näytä reagoivan kannusteisiin. Tukeen oikeutettujen kotitalouksien asumisvalinnat suhteessa pinta-ala ja laaturajoitteisiin vastaavat muiden kotitalouksien valintoja ja asunnonvaihdon mahdollistama potentiaalinen asumistuen lisäys ei nosta muuttotodennäköisyyttä. Muuttamiseen liittyvät kustannukset ja vajavaiset tiedot tukijärjestelmästä saattavat selittää heikkoa reagointia asumistuen luomiin kannustimiin. Toinen mahdollinen selitys on asumistuen vajaakäyttö. Tutkimuksen mukaan vain 70–80 prosenttia asumistukeen oikeutetuista kotitalouksista nostaa tukea. Asumistuen hyödyntämisen todennäköisyys riippuu koulutustasosta, tuen määrästä ja tulo-odotuksista.

The Line Spectral Frequency Model of a Finite-Length Sequence

The line spectral frequency (LSF) of a causal finite length sequence is a frequency at which the spectrum of the sequence annihilates or the magnitude spectrum has a spectral null. A causal finite-length sequencewith (L + 1) samples having exactly L-LSFs, is referred as an Annihilating (AH) sequence. Using some spectral properties of finite-length sequences, and some model parameters, we develop spectral decomposition structures, which are used to translate any finite-length sequence to an equivalent set of AH-sequences defined by LSFs and some complex constants. This alternate representation format of any finite-length sequence is referred as its LSF-Model. For a finite-length sequence, one can obtain multiple LSF-Models by varying the model parameters. The LSF-Model, in time domain can be used to synthesize any arbitrary causal finite-length sequence in terms of its characteristic AH-sequences. In the frequency domain, the LSF-Model can be used to obtain the spectral samples of the sequence as a linear combination of spectra of its characteristic AH-sequences. We also summarize the utility of the LSF-Model in practical discrete signal processing systems.

A kinetic Ising model study of dynamical correlations in confined fluids: Emergence of both fast and slow time scales

Experiments and computer simulation studies have revealed existence of rich dynamics in the orientational relaxation of molecules in confined systems such as water in reverse micelles, cyclodextrin cavities, and nanotubes. Here we introduce a novel finite length one dimensional Ising model to investigate the propagation and the annihilation of dynamical correlations in finite systems and to understand the intriguing shortening of the orientational relaxation time that has been reported for small sized reverse micelles. In our finite sized model, the two spins at the two end cells are oriented in the opposite directions to mimic the effects of surface that in real system fixes water orientation in the opposite directions. This produces opposite polarizations to propagate inside from the surface and to produce bulklike condition at the center. This model can be solved analytically for short chains. For long chains, we solve the model numerically with Glauber spin flip dynamics (and also with Metropolis single-spin flip Monte Carlo algorithm). We show that model nicely reproduces many of the features observed in experiments. Due to the destructive interference among correlations that propagate from the surface to the core, one of the rotational relaxation time components decays faster than the bulk. In general, the relaxation of spins is nonexponential due to the interplay between various interactions. In the limit of strong coupling between the spins or in the limit of low temperature, the nature of relaxation of the spins undergoes a qualitative change with the emergence of a homogeneous dynamics where decay is predominantly exponential, again in agreement with experiments. (C) 2010 American Institute of Physics. doi: 10.1063/1.3474948]

Aperiodic Variability in the Cane-Zebiak Model: A Diagnostic Study

The finite predictability of the coupled ocean-atmosphere system is determined by its aperiodic variability. To gain insight regarding the predictability of such a system, a series of diagnostic studies has been carried out to investigate the role of convergence feedback in producing the aperiodic behavior of the standard version of the Cane-Zebiak model. In this model, an increase in sea surface temperature (SST) increases atmospheric heating by enhancing local evaporation (SST anomaly feedback) and low-level convergence (convergence feedback). The convergence feedback is a nonlinear function of the background mean convergence field. For the set of standard parameters used in the model, it is shown that the convergence feedback contributes importantly to the aperiodic behaviour of the model. As the strength of the convergence feedback is increased from zero to its standard value, the model variability goes from a periodic regime to an aperiodic regime through a broadening of the frequency spectrum around the basic periodicity of about 4 years. Examination of the forcing associated with the convergence feedback reveals that it is intermittent, with relatively large amplitude only during 2 or 3 months in the early part of the calendar year. This seasonality in the efficiency of the convergence feedback is related to the strong seasonality of the mean convergence over the eastern Pacific. It is shown that if the mean convergence field is fixed at its March value, aperiodic behavior is produced even in the absence of annual cycles in the other mean fields. On the, other hand, if the mean convergence field is fixed at its September value, the coupled model evolution remains close to periodic, even in the presence of the annual cycle in the other fields. The role of convergence feedback on the aperiodic variability of the model for other parameter regimes is also examined. It is shown that a range exists in the strength of the SST anomaly feedback for which the model variability is aperiodic even without the convergence feedback. It appears that in the absence of convergence feedback, enhancement of the strength of the air-sea coupling in the model through other physical processes also results in aperiodicity in the model.

Reconstructing Solid Model from 2D Scanned Images of Biological Organs for Finite Element Simulation

This work presents a methodology to reconstruct 3D biological organs from image sequences or other scan data using readily available free softwares with the final goal of using the organs (3D solids) for finite element analysis. The methodology deals with issues such as segmentation, conversion to polygonal surface meshes, and finally conversion of these meshes to 3D solids. The user is able to control the detail or the level of complexity of the solid constructed. The methodology is illustrated using 3D reconstruction of a porcine liver as an example. Finally, the reconstructed liver is imported into the commercial software ANSYS, and together with a cyst inside the liver, a nonlinear analysis performed. The results confirm that the methodology can be used for obtaining 3D geometry of biological organs. The results also demonstrate that the geometry obtained by following this methodology can be used for the nonlinear finite element analysis of organs. The methodology (or the procedure) would be of use in surgery planning and surgery simulation since both of these extensively use finite elements for numerical simulations and it is better if these simulations are carried out on patient specific organ geometries. Instead of following the present methodology, it would cost a lot to buy a commercial software which can reconstruct 3D biological organs from scanned image sequences.