Mechanical analysis of foot plantar fascia during normal walking condition

Foot plantar fascia is an important foot tissue in stabilizing the longitudinal arch of human foot. Direct measurement to monitor the mechanical situation of plantar fascia at human locomotion is difficult. The purpose of this study was to construct a three-dimensional finite element model of the foot to calculate the internal stress/strain value of plantar fascia during different stage of gait. The simulated stress distribution of plantar fascia was the lowest at heel-strike, which concentrated on the medial side of calcaneal tubercle. The peak stress of plantar fascia was appeared at push-off, and the value is more than 5 times of the heel-strike position. Current FE model was able to explore the plantar fascia tension trend at the main sub-phases of foot. More detailed fascia model and intrinsic muscle forces could be developed in the further study.

MRI-based fatigue analysis predicts plaque vulnerability: A study comparing symptomatic and asymptomatic individuals

BACKGROUND: Rupture of atheromatous plaque in the carotid artery often leads to thrombosis and subsequent stroke. The mechanism of plaque rupture is not entirely clear but is thought to be a multi-factorial process involving thinning and weakening of the fibrous cap and biomechanical stress as the trigger leading to plaque rupture. As the cardiovascular system is a classic fatigue environment, the weakening of plaque leading to rupture may be a fatigue process, which is a symptomatically quiescent but potentially progressive failure process. In this study, we used a fatigue analysis based on in vivo magnetic resonance imaging (MRI) to investigate the rupture initiation location, crack propagation path and fatigue life within plaques of asymptomatic and symptomatic individuals. METHODS: Forty non-consecutive subjects (20 symptomatic and 20 asymptomatic) underwent high-resolution multi-sequence in vivo MRI of the carotid bifurcation. Fatigue analysis was performed based on the plaque geometry derived from in vivo MRI of the carotid artery at the point of maximum stenosis. Paris’ Law in fracture mechanics is adopted to determine the fatigue crack growth rate. Incremental crack propagation was dynamically simulated based on stress distributions. Plaque initiation location, crack propagation path and fatigue cycle of symptomatic and asymptomatic individuals were compared. RESULTS: Cracks were often found to begin at the lumen wall at areas of stress concentration. The preferred rupture direction was radial from the lumen center. The crack initially advanced slowly but accelerated as it developed, depending on plaque morphology. The fatigue cycles of symptomatic plaques were significantly less than those in the asymptomatic group (2.3 ± 0.9 vs 3.1 ± 0.7 (x106); p = 0.003). CONCLUSIONS: The number of cycles to rupture in symptomatic patients was higher than those predicted in asymptomatic patients by fatigue analysis, suggesting the possibility that plaques with a less fatigue life may be more prone to be symptomatic and rupture. If further validated by large-scale longitudinal studies, fatigue analysis based on high resolution in vivo MRI could potentially act as a useful tool for risk assessment of carotid atheroma.

Finite element analysis of deep transverse metatarsal ligaments mechanical response during landing

The deep transverse metatarsal ligaments (DTML) play an important role in stabilizing the metatarsal bones and manipulating foot transverse arch deformation. However, the biomechanical research about DTML in the foot maneuver is quite few. Due to the difficulties and lack of better measurement technology for these ligaments experimental monitor, the load transfer mechanism and internal stress state also hadn't been well addressed. The purpose of this study was to develop a detailing foot finite element model including DTML tissues, to investigate the mechanical response of DTML during the landing condition. The DTML was considered as hyperelastic material model was used to represent the nonlinear and nearly incompressible nature of the ligament tissue. From the simulation results, it is clearly to find that the peak maiximal principal stress of DTML was between the third and fourth metatarsals. Meanwhile, it seems the DTML in the middle position experienced higher tension than the sides DTML.

Association between aneurysm shoulder stress and abdominal aortic aneurysm expansion: A longitudinal follow-up study

Background Aneurysm expansion rate is an important indicator of the potential risk of abdominal aortic aneurysm (AAA) rupture. Stress within the AAA wall is also thought to be a trigger for its rupture. However, the association between aneurysm wall stresses and expansion of AAA is unclear. Methods and Results Forty-four patients with AAAs were included in this longitudinal follow-up study. They were assessed by serial abdominal ultrasonography and computed tomography scans if a critical size was reached or a rapid expansion occurred. Patient-specific 3-dimensional AAA geometries were reconstructed from the follow-up computed tomography images. Structural analysis was performed to calculate the wall stresses of the AAA models at both baseline and final visit. A nonlinear large-strain finite element method was used to compute the wall-stress distribution. The relationship between wall stresses and expansion rate was investigated. Slowly and rapidly expanding aneurysms had comparable baseline maximum diameters (median, 4.35 cm [interquartile range, 4.12 to 5.0 cm] versus 4.6 cm [interquartile range, 4.2 to 5.0 cm]; P=0.32). Rapidly expanding AAAs had significantly higher shoulder stresses than slowly expanding AAAs (median, 300 kPa [interquartile range, 280 to 320 kPa] versus 225 kPa [interquartile range, 211 to 249 kPa]; P=0.0001). A good correlation between shoulder stress at baseline and expansion rate was found (r=0.71; P=0.0001). Conclusion A higher shoulder stress was found to have an association with a rapidly expanding AAA. Therefore, it may be useful for estimating the expansion of AAAs and improve risk stratification of patients with AAAs.

Computed wall stress may predict the growth of abdominal aortic aneurysm

Growth rate of abdominal aortic aneurysm (AAA) is thought to be an important indicator of the potential risk of rupture. Wall stress is also thought to be a trigger for its rupture. However, stress change during the expansion of an AAA is unclear. Forty-four patients with AAAs were included in this longitudinal follow-up study. They were assessed by serial abdominal ultrasonography and computerized tomography (CT) scans if a critical size was reached or a rapid expansion occurred. Patient-specific 3-dimensional AAA geometries were reconstructed from the follow-up CT images. Structural analysis was performed to calculate the wall stresses of the AAA models at both baseline and final visit. A non-linear large-strain finite element method was used to compute the wall stress distribution. The average growth rate was 0.66cm/year (range 0-1.32 cm/year). A significantly positive correlation between shoulder tress at baseline and growth rate was found (r=0.342; p=0.02). A higher shoulder stress is associated with a rapidly expanding AAA. Therefore, it may be useful for estimating the growth expansion of AAAs and further risk stratification of patients with AAAs.

Energy-momentum conserving algorithm for nonlinear transient analysis within the framework of hybrid elements

This work deals with the formulation and implementation of an energy-momentum conserving algorithm for conducting the nonlinear transient analysis of structures, within the framework of stress-based hybrid elements. Hybrid elements, which are based on a two-field variational formulation, are much less susceptible to locking than conventional displacement-based elements within the static framework. We show that this advantage carries over to the transient case, so that not only are the solutions obtained more accurate, but they are obtained in fewer iterations. We demonstrate the efficacy of the algorithm on a wide range of problems such as ones involving dynamic buckling, complicated three-dimensional motions, et cetera.

The impact of wall shear stress and pressure drop on the stability of the atherosclerotic plaque

Rupture of vulnerable atheromatous plaque in the carotid and coronary arteries often leads to stroke and heart attack respectively. The mechanism of blood flow and plaque rupture in stenotic arteries is still not fully understood. A three dimensional rigid wall model was solved under steady state conditions and unsteady conditions by assuming a time-varying inlet velocity profile to investigate the relative importance of axial forces and pressure drops in arteries with asymmetric stenosis. Flow-structure interactions were investigated for the same geometry and the results were compared with those retrieved with the corresponding 2D cross-section structural models. The Navier-Stokes equations were used as the governing equations for the fluid. The tube wall was assumed hyperelastic, homogeneous, isotropic and incompressible. The analysis showed that the three dimensional behavior of velocity, pressure and wall shear stress is in general very different from that predicted by cross-section models. Pressure drop across the stenosis was found to be much higher than shear stress. Therefore, pressure may be the more important mechanical trigger for plaque rupture other than shear stress, although shear stress is closely related to plaque formation and progression.

Identifying vulnerable carotid plaques in vivo using high resolution magnetic resonance imaging-based finite element analysis

Object. Individuals with carotid atherosclerosis develop symptoms following rupture of vulnerable plaques. Biomechanical stresses within this plaque may increase vulnerability to rupture. In this report the authors describe the use of in vivo carotid plaque imaging and computational mechanics to document the magnitude and distribution of intrinsic plaque stresses. Methods. Ten (five symptomatic and five asymptomatic) individuals underwent plaque characterization magnetic resonance (MR) imaging. Plaque geometry and composition were determined by multisequence review. Intrinsic plaque stress profiles were generated from 3D meshes by using finite element computational analysis. Differences in principal (shear) stress between normal and diseased sections of the carotid artery and between symptomatic and asymptomatic plaques were noted. Results. There was a significant difference in peak principal stress between diseased and nondiseased segments of the artery (mean difference 537.65 kPa, p < 0.05). Symptomatic plaques had higher mean stresses than asymptomatic plaques (627.6 kPa compared with 370.2 kPa, p = 0.05), which were independent of luminal stenosis and plaque composition. Conclusions. Significant differences in plaque stress exist between plaques from symptomatic individuals and those from asymptomatic individuals. The MR imaging-based computational analysis may therefore be a useful aid to identification of vulnerable plaques in vivo.

Correlation of macrophage location and plaque stress distribution using USPIO-enhanced MRI in a patient with symptomatic severe carotid stenosis: A new insight into risk stratification

High resolution, USPIO-enhanced MR imaging can be used to identify inflamed atherosclerotic plaque. We report a case of a 79-year-old man with a symptomatic carotid stenosis of 82%. The plaque was retrieved for histology and finite element analysis (FEA) based on the preoperative MR imaging was used to predict maximal Von Mises stress on the plaque. Macrophage location correlated with maximal predicted stresses on the plaque. This supports the hypothesis that macrophages thin the fibrous cap at points of highest stress, leading to an increased risk of plaque rupture and subsequent stroke.

Atheroma: Is calcium important or not? A modelling study of stress within the atheromatous fibrous cap in relation to position and size of calcium deposits

Atheromatous plaque rupture h the cause of the majority of strokes and heart attacks in the developed world. The role of calcium deposits and their contribution to plaque vulnerability are controversial. Some studies have suggested that calcified plaque tends to be more stable whereas others have suggested the opposite. This study uses a finite element model to evaluate the effect of calcium deposits on the stress within the fibrous cap by varying their location and size. Plaque fibrous cap, lipid pool and calcification were modeled as hyperelastic, Isotropic, (nearly) incompressible materials with different properties for large deformation analysis by assigning time-dependent pressure loading on the lumen wall. The stress and strain contours were illustrated for each condition for comparison. Von Mises stress only increases up to 1.5% when varying the location of calcification in the lipid pool distant to the fibrous cap. Calcification in the fibrous cap leads to a 43% increase of Von Mises stress when compared with that in the lipid pool. An increase of 100% of calcification area leads to a 15% stress increase in the fibrous cap. Calcification in the lipid pool does not increase fibrous cap stress when it is distant to the fibrous cap, whilst large areas of calcification close to or in the fibrous cap may lead to a high stress concentration within the fibrous cap, which may cause plaque rupture. This study highlights the application of a computational model on a simulation of clinical problems, and it may provide insights into the mechanism of plaque rupture.

Analysis of a finite composite plate with smooth rigid pin

A continuum method of analysis is presented in this paper for the problem of a smooth rigid pin in a finite composite plate subjected to uniaxial loading. The pin could be of interference, push or clearance fit. The plate is idealized to an orthotropic sheet. As the load on the plate is progressively increased, the contact along the pin-hole interface is partial above certain load levels in all three types of fit. In misfit pins (interference or clearance), such situations result in mixed boundary value problems with moving boundaries and in all of them the arc of contact and the stress and displacement fields vary nonlinearly with the applied load. In infinite domains similar problems were analysed earlier by ‘inverse formulation’ and, now, the same approach is selected for finite plates. Finite outer domains introduce analytical complexities in the satisfaction of boundary conditions. These problems are circumvented by adopting a method in which the successive integrals of boundary error functions are equated to zero. Numerical results are presented which bring out the effects of the rectangular geometry and the orthotropic property of the plate. The present solutions are the first step towards the development of special finite elements for fastener joints.

Predicting survival of western rock lobsters Panulirus cygnus using discriminant analysis of hemolymph parameters taken immediately following simulated handling treatments

Instances of morbidity amongst rock lobsters (Panulirus cygnus) arriving at factories in Western Australia (WA) have been attributed to stress during post-harvest handling. This study used discriminant analysis to determine whether physiological correlates of stress following a period of simulated post-harvest handling had any validity as predictors of future rejection or morbidity of western rock lobsters. Groups of 230 western rock lobsters were stored for 6 h in five environments (submerged/flowing sea water, submerged/re-circulating sea water, humid air, flowing sea water spray, and re-circulated sea water spray). The experiment was conducted in late spring (ambient sea water 22°C), and repeated again in early autumn (ambient sea water 26°C). After 6 h treatment, each lobster was graded for acceptability for live export, numbered, and its hemolymph was sampled. The samples were analysed for a number of physiological and health status parameters. The lobsters were then stored for a week in tanks in the live lobster factory to record mortality. The mortality of lobsters in the factory was associated with earlier deviations in hemolymph parameters as they emerged from the storage treatments. Discriminant analysis (DA) of the hemolymph assays enabled the fate of 80-90% of the lobsters to be correctly categorised within each experiment. However, functions derived from one experiment were less accurate at predicting mortality when applied to the other experiments. One of the reasons for this was the higher mortality and the more severe patho-physiological changes observed in lobsters stored in humid air or sprays at the higher temperature. The analysis identified lactate accumulation during emersion and associated physiological and hemocyte-related effects as a major correlate of mortality. Reducing these deviations, for example by submerged transport, is expected to ensure high levels of survival. None of the indicators tested predicted mortality with total accuracy. The simplest and most accurate means of comparing emersed treatments was to count the mortality afterwards.

On an analytical-numerical procedure for the analysis of cylindrical shells with arbitrarily oriented cracks

An analytical-numerical procedure for obtaining stress intensity factor solutions for an arbitrarily oriented crack in a long, thin circular cylindrical shell is presented. The method of analysis involves obtaining a series solution to the governing shell equation in terms of Mathieu and modified Mathieu functions by the method of separation of variables and satisfying the crack surface boundary conditions numerically using collocation. The solution is then transformed from elliptic coordinates to polar coordinates with crack tip as the origin through a Taylor series expansion and membrane and bending stress intensity factors are computed. Numerical results are presented and discussed for the pressure loading case.

A sequential analysis of "nú" and "núna" in Icelandic conversation

This thesis is an empirical study of how two words in Icelandic, "nú" and "núna", are used in contemporary Icelandic conversation. My aims in this study are, first, to explain the differences between the temporal functions of "nú" and "núna", and, second, to describe the non-temporal functions of "nú". In the analysis, a focus is placed on comparing the sequential placement of the two words, on their syntactical distribution, and on their prosodic realization. The empirical data comprise 14 hours and 11 minutes of naturally occurring conversation recorded between 1996 and 2003. The selected conversations represent a wide range of interactional contexts including informal dinner parties, institutional and non-institutional telephone conversations, radio programs for teenagers, phone-in programs, and, finally, a political debate on television. The theoretical and methodological framework is interactional linguistics, which can be described as linguistically oriented conversation analysis (CA). A comparison of "nú" and "núna" shows that the two words have different syntactic distributions. "Nú" has a clear tendency to occur in the front field, before the finite verb, while "núna" typically occurs in the end field, after the object. It is argued that this syntactic difference reflects a functional difference between "nú" and "núna". A sequential analysis of "núna" shows that the word refers to an unspecified period of time which includes the utterance time as well as some time in the past and in the future. This temporal relation is referred to as reference time. "Nú", by contrast, is mainly used in three different environments: a) in temporal comparisons, 2) in transitions, and 3) when the speaker is taking an affective stance. The non-temporal functions of "nú" are divided into three categories: a) "nú" as a tone particle, 2) "nú" as an utterance particle, and 3) "nú" as a dialogue particle. "Nú" as a tone particle is syntactically integrated and can occur in two syntactic positions: pre-verbally and post-verbally. I argue that these instances are employed in utterances in which a speaker is foregrounding information or marking it as particularly important. The study shows that, although these instances are typically prosodically non-prominent and unstressed, they are in some cases delivered with stress and with a higher pitch than the surrounding talk. "Nú" as an utterance particle occurs turn-initially and is syntactically non-integrated. By using "nú", speakers show continuity between turns and link new turns to prior ones. These instances initiate either continuations by the same speaker or new turns after speaker shifts. "Nú" as a dialogue particle occurs as a turn of its own. The study shows that these instances register informings in prior turns as unexpected or as a departure from the normal state of affairs. "Nú" as a dialogue particle is often delivered with a prolonged vowel and a recognizable intonation contour. A comparative sequential and prosodic analysis shows that in these cases there is a correlation between the function of "nú" and the intonation contour by which it is delivered. Finally, I argue that despite the many functions of "nú", all the instances can be said to have a common denominator, which is to display attention towards the present moment and the utterances which are produced prior or after the production of "nú". Instead of anchoring the utterances in external time or reference time, these instances position the utterance in discourse internal time, or discourse time.

Remote sensing of nitrogen and water stress in wheat

Nitrogen (N) is the largest agricultural input in many Australian cropping systems and applying the right amount of N in the right place at the right physiological stage is a significant challenge for wheat growers. Optimizing N uptake could reduce input costs and minimize potential off-site movement. Since N uptake is dependent on soil and plant water status, ideally, N should be applied only to areas within paddocks with sufficient plant available water. To quantify N and water stress, spectral and thermal crop stress detection methods were explored using hyperspectral, multispectral and thermal remote sensing data collected at a research field site in Victoria, Australia. Wheat was grown over two seasons with two levels of water inputs (rainfall/irrigation) and either four levels (in 2004; 0, 17, 39 and 163 kg/ha) or two levels (in 2005; 0 and 39 kg/ha N) of nitrogen. The Canopy Chlorophyll Content Index (CCCI) and modified Spectral Ratio planar index (mSRpi), two indices designed to measure canopy-level N, were calculated from canopy-level hyperspectral data in 2005. They accounted for 76% and 74% of the variability of crop N status, respectively, just prior to stem elongation (Zadoks 24). The Normalised Difference Red Edge (NDRE) index and CCCI, calculated from airborne multispectral imagery, accounted for 41% and 37% of variability in crop N status, respectively. Greater scatter in the airborne data was attributable to the difference in scale of the ground and aerial measurements (i.e., small area plant samples against whole-plot means from imagery). Nevertheless, the analysis demonstrated that canopy-level theory can be transferred to airborne data, which could ultimately be of more use to growers. Thermal imagery showed that mean plot temperatures of rainfed treatments were 2.7 °C warmer than irrigated treatments (P < 0.001) at full cover. For partially vegetated fields, the two-Dimensional Crop Water Stress Index (2D CWSI) was calculated using the Vegetation Index-Temperature (VIT) trapezoid method to reduce the contribution of soil background to image temperature. Results showed rainfed plots were consistently more stressed than irrigated plots. Future work is needed to improve the ability of the CCCI and VIT methods to detect N and water stress and apply both indices simultaneously at the paddock scale to test whether N can be targeted based on water status. Use of these technologies has significant potential for maximising the spatial and temporal efficiency of N applications for wheat growers. ‘Ground–breaking Stuff’- Proceedings of the 13th Australian Society of Agronomy Conference, 10-14 September 2006, Perth, Western Australia.