Statistical finite element model for bone shape and biomechanical properties

We present a framework for statistical finite element analysis combining shape and material properties, and allowing performing statistical statements of biomechanical performance across a given population. In this paper, we focus on the design of orthopaedic implants that fit a maximum percentage of the target population, both in terms of geometry and biomechanical stability. CT scans of the bone under consideration are registered non-rigidly to obtain correspondences in position and intensity between them. A statistical model of shape and intensity (bone density) is computed by means of principal component analysis. Afterwards, finite element analysis (FEA) is performed to analyse the biomechanical performance of the bones. Realistic forces are applied on the bones and the resulting displacement and bone stress distribution are calculated. The mechanical behaviour of different PCA bone instances is compared.

Effects of glenoid component version on humeral head displacement and joint reaction forces: an experimental study

The purpose of this study was to determine whether changes in glenoid version are associated with humeral head displacement and changes in the joint reaction forces, as these might contribute to instability or loosening in total shoulder replacement. A total shoulder prosthesis was implanted in neutral version in 6 cadaveric shoulders. Glenoid version was then changed in steps of 4 degrees toward more anteversion and retroversion. An increase in anteversion resulted in anterior translation of the humeral head and in eccentric loading of the anterior part of the glenoid. Retroversion was associated with posterior displacement and posterior loading of the glenoid. A change in rotation of the humeral component did not compensate for altered version of the glenoid component. These results suggest that both instability and glenoid component loosening may be related to the version of the glenoid component. Therefore, assessment of loosening and instability justifies precise assessment of glenoid component version.

The correlation of urethral resistance pressure with maximum urethral closure pressure and stress incontinence

OBJECTIVE: Aim of the study was to correlate urethral retro resistance pressure with the maximum urethral closure pressure (MUCP) and functional urethral length (FUL) in patients with urinary incontinence and healthy individuals. STUDY DESIGN: Two hundred and twenty patients with the complaint of urinary incontinence had a urodynamic examination including urethral pressure profiles and URP. Additionally, 15 healthy individuals without the complaint of any incontinence had their URP and urethral pressure profiles measured. The correlation of MUCP, FUL and URP were calculated using Graph Pad Instat 4.0 for windows. RESULTS: URP correlates well with the diagnosis of urodynamic stress incontinence. Correlation coefficient between URP and MUCP is 0.9262. Healthy individuals have significantly higher values for URP and MUCP. CONCLUSION: URP is a valuable less invasive test than conventional urethral function tests for the diagnosis of urodynamic incontinence with an excellent correlation of MUCP and URP.

Gunshot-related displacement of skin particles and bacteria from the exit region back into the bullet path

In previous studies, it was shown that there is a gunshot-related transport of skin particles and microorganisms from the entrance region into the depth of the bullet path. The present study deals with the question of whether gunshots may also cause a retrograde transport of skin particles and microorganisms from the bullet exit region back into the bullet path. For this purpose, we used a composite model consisting of rectangular gelatin blocks and pig skin. The skin pieces were firmly attached to the gelatin blocks on the side where the bullet was to exit. Prior to the test shots, the outer surface of the pig skin was contaminated with a thin layer of a defined bacterial suspension. After drying the skin, test shots were fired from a distance of 10 m using cartridges calibre .38 spec. with different bullet types. Subsequent analyses showed that in all shots with full penetration of the composite model, the bullet path contained displaced skin particles and microorganisms from the skin surface at the exit site. These could be regularly detected in the distal 6-8 cm of the track, occasionally up to a distance of 18 cm from the exit hole. The distribution of skin particles and microorganisms is presented and the possible mechanism of this retrograde transport is discussed.

Isometric rate of force development, maximum voluntary contraction, and balance in women with and without joint hypermobility

OBJECTIVE: To determine differences between hypermobile subjects and controls in terms of maximum strength, rate of force development, and balance. METHODS: We recruited 13 subjects with hypermobility and 18 controls. Rate of force development and maximal voluntary contraction (MVC) during single leg knee extension of the right knee were measured isometrically for each subject. Balance was tested twice on a force plate with 15-second single-leg stands on the right leg. Rate of force development (N/second) and MVC (N) were extracted from the force-time curve as maximal rate of force development (= limit Deltaforce/Deltatime) and the absolute maximal value, respectively. RESULTS: The hypermobile subjects showed a significantly higher value for rate of force development (15.2% higher; P = 0.038, P = 0.453, epsilon = 0.693) and rate of force development related to body weight (16.4% higher; P = 0.018, P = 0.601, epsilon = 0.834) than the controls. The groups did not differ significantly in MVC (P = 0.767, P = 0.136, epsilon = 0.065), and MVC related to body weight varied randomly between the groups (P = 0.921, P = 0.050, epsilon = 0.000). In balance testing, the mediolateral sway of the hypermobile subjects showed significantly higher values (11.6% higher; P = 0.034, P = 0.050, epsilon = 0.000) than that of controls, but there was no significant difference (4.9% difference; P = 0.953, P = 0.050, epsilon = 0.000) in anteroposterior sway between the 2 groups. CONCLUSION: Hypermobile women without acute symptoms or limitations in activities of daily life have a higher rate of force development in the knee extensors and a higher mediolateral sway than controls with normal joint mobility.

Comparative study of anchorage strengths of epoxy coated hooked bars

Hooked reinforcing bars (rebar) are used frequently to carry the tension forces developed in beams and transferred to columns. Research into epoxy coated hooked bars has only been minimally performed and no research has been carried out incorporating the coating process found in ASTM A934. This research program compares hooked rebar that are uncoated, coated by ASTM A775, and coated by ASTM A934. In total, forty-two full size beam-column specimens were created, instrumented and tested to failure. The program was carried out in three phases. The first phase was used to refine the test setup and procedures. Phase two explored the spacing of column ties within the joint region. Phase three explored the three coating types found above. Each specimen included two hooked rebar which were loaded and measured independently for relative rebar slip. The load and displacement of the hooked rebar were analyzed, focusing on behavior at the levels of 30 ksi, 42 ksi and 60 ksi of rebar stress. Statistical and general comparisons were made using the coating types, tie spacing, and rebar stress level. Many of the parameters composing the rebar and concrete were also tested to characterize the components and specimens. All rebar tested met ASTM standards for tensile strength, but the newer ASTM A934 method seemed to produce slightly lower yield strengths. The A934 method also produced coating thicknesses that were very inconsistent and were higher than ASTM maximum limits in many locations. Continuity of coating surfaces was found to be less than 100% for both A775 and A934 rebar, but for different reasons. The many comparisons made did not always produce clear conclusions. The data suggests that the ACI Code (318-05) parameter of 1.2 for including epoxy coating on hooked rebar may need to be raised, possibly to 2.5, but more testing needs to be performed before such a large value change is set forth. This is particularly important as variables were identified which may have a larger influence on rebar capacity than the Development Length, of which the current 1.2 factor modifies. Many suggestions for future work are included throughout the thesis to help guide other researchers in carrying out successful and productive programs which will further the highly understudied topic of hooked rebar.

Mechanisms contributing to the deep chlorophyll maximum in Lake Superior

The seasonal appearance of a deep chlorophyll maximum (DCM) in Lake Superior is a striking phenomenon that is widely observed; however its mechanisms of formation and maintenance are not well understood. As this phenomenon may be the reflection of an ecological driver, or a driver itself, a lack of understanding its driving forces limits the ability to accurately predict and manage changes in this ecosystem. Key mechanisms generally associated with DCM dynamics (i.e. ecological, physiological and physical phenomena) are examined individually and in concert to establish their role. First the prevailing paradigm, “the DCM is a great place to live”, is analyzed through an integration of the results of laboratory experiments and field measurements. The analysis indicates that growth at this depth is severely restricted and thus not able to explain the full magnitude of this phenomenon. Additional contributing mechanisms like photoadaptation, settling and grazing are reviewed with a one-dimensional mathematical model of chlorophyll and particulate organic carbon. Settling has the strongest impact on the formation and maintenance of the DCM, transporting biomass to the metalimnion and resulting in the accumulation of algae, i.e. a peak in the particulate organic carbon profile. Subsequently, shade adaptation becomes manifest as a chlorophyll maximum deeper in the water column where light conditions particularly favor the process. Shade adaptation mediates the magnitude, shape and vertical position of the chlorophyll peak. Growth at DCM depth shows only a marginal contribution, while grazing has an adverse effect on the extent of the DCM. The observed separation of the carbon biomass and chlorophyll maximum should caution scientists to equate the DCM with a large nutrient pool that is available to higher trophic levels. The ecological significance of the DCM should not be separated from the underlying carbon dynamics. When evaluated in its entirety, the DCM becomes the projected image of a structure that remains elusive to measure but represents the foundation of all higher trophic levels. These results also offer guidance in examine ecosystem perturbations such as climate change. For example, warming would be expected to prolong the period of thermal stratification, extending the late summer period of suboptimal (phosphorus-limited) growth and attendant transport of phytoplankton to the metalimnion. This reduction in epilimnetic algal production would decrease the supply of algae to the metalimnion, possibly reducing the supply of prey to the grazer community. This work demonstrates the value of modeling to challenge and advance our understanding of ecosystem dynamics, steps vital to reliable testing of management alternatives.

ANKLE IMPEDANCE AND ANKLE ANGLES DURING STEP TURN AND STRAIGHT WALK: IMPLICATIONS FOR THE DESIGN OF A STEERABLE ANKLE-FOOT PROSTHETIC ROBOT

During locomotion, turning is a common and recurring event which is largely neglected in the current state-of-the-art ankle-foot prostheses, forcing amputees to use different steering mechanisms for turning, compared to non-amputees. A better understanding of the complexities surrounding lower limb prostheses will lead to increased health and well-being of amputees. The aim of this research is to develop a steerable ankle-foot prosthesis that mimics the human ankle mechanical properties. Experiments were developed to estimate the mechanical impedance of the ankle and the ankles angles during straight walk and step turn. Next, this information was used in the design of a prototype, powered steerable ankle-foot prosthesis with two controllable degrees of freedom. One of the possible approaches in design of the prosthetic robots is to use the human joints’ parameters, especially their impedance. A series of experiments were conducted to estimate the stochastic mechanical impedance of the human ankle when muscles were fully relaxed and co-contracting antagonistically. A rehabilitation robot for the ankle, Anklebot, was employed to provide torque perturbations to the ankle. The experiments were performed in two different configurations, one with relaxed muscles, and one with 10% of maximum voluntary contraction (MVC). Surface electromyography (sEMG) was used to monitor muscle activation levels and these sEMG signals were displayed to subjects who attempted to maintain them constant. Time histories of ankle torques and angles in the lateral/medial (LM) directions, inversion-eversion (IE), and dorsiflexionplantarflexion (DP) were recorded. Linear time-invariant transfer functions between the measured torques and angles were estimated providing an estimate of ankle mechanical impedance. High coherence was observed over a frequency range up to 30 Hz. The main effect of muscle activation was to increase the magnitude of ankle mechanical impedance in all degrees of freedom of the ankle. Another experiment compared the three-dimensional angles of the ankle during step turn and straight walking. These angles were measured to be used for developing the control strategy of the ankle-foot prosthesis. An infrared camera system was used to track the trajectories and angles of the foot and leg. The combined phases of heel strike and loading response, mid stance, and terminal stance and pre-swing were determined and used to measure the average angles at each combined phase. The Range of motion (ROM) in IE increased during turning while ML rotation decreased and DP changed the least. During the turning step, ankle displacement in DP started with similar angles to straight walk and progressively showed less plantarflexion. In IE, the ankle showed increased inversion leaning the body toward the inside of the turn. ML rotation initiated with an increased medial rotation during the step turn relative to the straight walk transitioning to increased lateral rotation at the toe off. A prototype ankle-foot prosthesis capable of controlling both DP and IE using a cable driven mechanism was developed and assessed as part of a feasibility study. The design is capable of reproducing the angles required for straight walk and step turn; generates 712N of lifting force in plantarflexion, and shows passive stiffness comparable to a nonload bearing ankle impedance. To evaluate the performance of the ankle-foot prosthesis, a circular treadmill was developed to mimic human gait during steering. Preliminary results show that the device can appropriately simulate human gait with loading and unloading the ankle joint during the gait in circular paths.

Maximum principle preserving high order schemes for convection-dominated diffusion equations

The maximum principle is an important property of solutions to PDE. Correspondingly, it's of great interest for people to design a high order numerical scheme solving PDE with this property maintained. In this thesis, our particular interest is solving convection-dominated diffusion equation. We first review a nonconventional maximum principle preserving(MPP) high order finite volume(FV) WENO scheme, and then propose a new parametrized MPP high order finite difference(FD) WENO framework, which is generalized from the one solving hyperbolic conservation laws. A formal analysis is presented to show that a third order finite difference scheme with this parametrized MPP flux limiters maintains the third order accuracy without extra CFL constraint when the low order monotone flux is chosen appropriately. Numerical tests in both one and two dimensional cases are performed on the simulation of the incompressible Navier-Stokes equations in vorticity stream-function formulation and several other problems to show the effectiveness of the proposed method.