Hip power analysis in individuals with transfemoral amputation: A different strategy different from stabilisation during gait stance

Introduction and Objectives Joint moments and joint powers during gait are widely used to determine the effects of rehabilitation programs as well as prosthetic fitting. Following the definition of power (dot product of joint moment and joint angular velocity) it has been previously proposed to analyse the 3D angle between both vectors, αMw. Basically, joint power is maximised when both vectors are parallel and cancelled when both vectors are orthogonal. In other words, αMw < 60° reveals a propulsion configuration (more than 50% of the moment contribute to positive power) while αMw > 120° reveals a resistance configuration (more than 50% of the moment contribute to negative power). A stabilisation configuration (less than 50% of the moment contribute to power) corresponds to 60° < αMw < 120°. Previous studies demonstrated that hip joints of able-bodied adults (AB) are mainly in a stabilisation configuration (αMw about 90°) during the stance phase of gait. [1, 2] Individuals with transfemoral amputation (TFA) need to maximise joint power at the hip while controlling the prosthetic knee during stance. Therefore, we tested the hypothesis that TFAs should adopt a strategy that is different from a continuous stabilisation. The objective of this study was to compute joint power and αMw for TFA and to compare them with AB. Methods Three trials of walking at self-selected speed were analysed for 8 TFAs (7 males and 1 female, 46±10 years old, 1.78±0.08 m 82±13 kg) and 8 ABs (males, 25±3 years old, 1.75±0.04, m 67±6 kg). The joint moments are computed from a motion analysis system (Qualisys, Goteborg, Sweden) and a multi-axial transducer (JR3, Woodland, USA) mounted above the prosthetic knee for TFAs and from a motion analysis system (Motion Analysis, Santa Rosa, USA) and force plates (Bertec, Columbus, USA) for ABs. The TFAs were fitted with an OPRA (Integrum, AB, Gothengurg, Sweden) osseointegrated implant system and their prosthetic designs include pneumatic, hydraulic and microprocessor knees. Previous studies showed that the inverse dynamics computed from the multi-axial transducer is the proper method considering the absorption at the foot and resistance at the knee. Results The peak of positive power at loading response (H1) was earlier and lower for TFA compared to AB. Although the joint power is lower, the 3D angle between joint moment and joint angular velocity, αMw, reveals an obvious propulsion configuration (mean αMw about 20°) for TFA compared to a stabilisation configuration (mean αMw about 70°) for AB. The peaks of negative power at midstance (H2) and of positive power at preswing / initial swing (H3) occurred later, lower and longer for TFA compared to AB. Again, the joint powers are lower for TFA but, in this case, αMw is almost comparable (with a time lag), demonstrating a stabilisation (almost a resistance for TFA, mean αMw about 120°) and a propulsion configuration, respectively. The swing phase is not analysed in the present study. Conclusion The analysis of hip joint power may indicate that TFAs demonstrated less propulsion and resistance than ABs during the stance phase of gait. This is true from a quantitative point of view. On the contrary, the 3D angle between joint moment and joint angular velocity, αMw, reveals that TFAs have a remarkable propulsion strategy at loading response and almost a resistance strategy at midstance while ABs adopted a stabilisation strategy. The propulsion configuration, with αMw close to 0°, seems to aim at maximising the positive joint power. The configuration close to resistance, with αMw far from 180°, might aim at unlocking the prosthetic knee before swing while minimising the negative power. This analysis of both joint power and 3D angle between the joint moment and the joint angular velocity provides complementary insights into the gait strategies of TFA that can be used to support evidence-based rehabilitation and fitting of prosthetic components.

Some mechanisms of excitation of a railway wheel

Railway wheel vibrations are caused by a number of mechanisms. Two of these are considered: (a) gravitational load reaction acting on different points of the wheel rim, as the wheel rolls on, and (b) random fluctuating forces generated at the contact patch by roughness on the mating surfaces of the wheel and rail. The wheel is idealized as a thin ring, and the analysis is limited to a single wheel rolling on a rail. It is shown that the first mechanism results in a stationary pattern of vibration, which would not radiate any sound. The acceleration caused by roughness-excited forces is much higher at higher frequencies, but is of the same order as that caused by load reaction at lower frequencies. The computed acceleration level (and hence the radiated SPL) caused by roughness is comparable with the observed values, and is seen to increase by about 10 dB for a doubling of the wagon speed. The driving point impedance of the periodic rail-sleeper system at the contact patch, which is used in the analysis, is derived in a companion paper.

A new formalism for representation of binary thermodynamic data

The applicability of a formalism involving an exponential function of composition x1 in interpreting the thermodynamic properties of alloys has been studied. The excess integral and partial molar free energies of mixing are expressed as: $$\begin{gathered} \Delta F^{xs} = a_o x_1 (1 - x_1 )e^{bx_1 } \hfill \\ RTln\gamma _1 = a_o (1 - x_1 )^2 (1 + bx_1 )e^{bx_1 } \hfill \\ RTln\gamma _2 = a_o x_1^2 (1 - b + bx_1 )e^{bx_1 } \hfill \\ \end{gathered} $$ The equations are used in interpreting experimental data for several relatively weakly interacting binary systems. For the purpose of comparison, activity coefficients obtained by the subregular model and Krupkowski’s formalism have also been computed. The present equations may be considered to be convenient in describing the thermodynamic behavior of metallic solutions.

Electron-electron interactions in polyacetylene

Experimental evidence for strong electron-electron interactions in polyacetylene is presented. These include (i) observation of a dipole forbidden state below the optical gap, (ii) observation of negative spin densities at sites at which noninteracting models predict zero spin density (iii) vanishing optical gap, in the infinite chain limit, in the closely related symmetrical linear cyanine dyes. To correctly explain these features it is necessary to solve correlated model Hamiltonians. Using diagrammatic valence bond method model exact solutions of correlated models of finite-size systems can be obtained and various physical properties of the low-lying states can be computed. These properties, when extrapolated to the infinite chain limit explain many of the experimental features observed in polyacetylene.

Fracture Behavior of Ferrocement Beams

An attempt is made to study the fracture behavior of ferrocement beams using J-integral and critical crack opening displacement approaches. Ferrocement beams with three different relative notch depths and different percentages of mesh reinforcement were tested in four-point bending (third-point loading). The experimental results were used to evaluate the apparent J-integral and CODc values. Results show that the apparent J-integral does not seem to follow any particular trend in variation with notch depth, but is sensitive to the increase of mesh reinforcement. Hence, the apparent J-integral appears to be a useful fracture criterion for ferrocement. The computed values of CODt are found to be dependent on the depth of notch and, thus, cannot possibly be considered as a suitable fracture criterion for ferrocement.

Complexity analysis of EEG in patients with schizophrenia using fractal dimension

We computed Higuchi's fractal dimension (FD) of resting, eyes closed EEG recorded from 30 scalp locations in 18 male neuroleptic-naive, recent-onset schizophrenia (NRS) subjects and 15 male healthy control (HC) subjects, who were group-matched for age. Schizophrenia patients showed a diffuse reduction of FD except in the bilateral temporal and occipital regions, with the reduction being most prominent bifrontally. The positive symptom (PS) schizophrenia subjects showed FD values similar to or even higher than HC in the bilateral temporo-occipital regions, along with a co-existent bifrontal FD reduction as noted in the overall sample of NRS. In contrast, this increase in FD values in the bilateral temporo-occipital region was absent in the negative symptom (NS) subgroup. The regional differences in complexity suggested by these findings may reflect the aberrant brain dynamics underlying the pathophysiology of schizophrenia and its symptom dimensions. Higuchi's method of measuring FD directly in the time domain provides an alternative for the more computationally intensive nonlinear methods of estimating EEG complexity.

A conditionally linearized Monte Carlo filter in non-linear structural dynamics

State and parameter estimations of non-linear dynamical systems, based on incomplete and noisy measurements, are considered using Monte Carlo simulations. Given the measurements. the proposed method obtains the marginalized posterior distribution of an appropriately chosen (ideally small) subset of the state vector using a particle filter. Samples (particles) of the marginalized states are then used to construct a family of conditionally linearized system of equations and thus obtain the posterior distribution of the states using a bank of Kalman filters. Discrete process equations for the marginalized states are derived through truncated Ito-Taylor expansions. Increased analyticity and reduced dispersion of weights computed over a smaller sample space of marginalized states are the key features of the filter that help achieve smaller sample variance of the estimates. Numerical illustrations are provided for state/parameter estimations of a Duffing oscillator and a 3-DOF non-linear oscillator. Performance of the filter in parameter estimation is also assessed using measurements obtained through experiments on simple models in the laboratory. Despite an added computational cost, the results verify that the proposed filter generally produces estimates with lower sample variance over the standard sequential importance sampling (SIS) filter.

Determination of baryon magnetic moments from QCD sum rules

The magnetic moment μB of a baryon B with quark content (aab) is written as μB=4ea(1+δB)eħ/2cMB, where ea is the charge of the quark of flavor type a. The experimental values of δB have a simple pattern and have a natural explanation within QCD. Using the ratio method, the QCD sum rules are analyzed and the values of δB are computed. We find good agreement with data (≊10%) for the nucleons and the Σ multiplet while for the cascade the agreement is not as good. In our analysis we have incorporated additional terms in the operator-product expansion as compared to previous authors. We also clarify some points of disagreement between the previous authors. External-field-induced correlations describing the magnetic properties of the vacuum are estimated from the baryon magnetic-moment sum rules themselves as well as by independent spectral representations and the results are contrasted.

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.

The harmonic torsional oscillations of a thin disk submerged in a fluid with a surfactant surface layer

A formulation in terms of a Fredholm integral equation of the first kind is given for the axisymmetric problem of a disk oscillating harmonically in a viscous fluid whose surface is contaminated with a surfactant film. The equation of the first kind is converted to a pair of coupled integral equations of the second kind, which are solved numerically. The resistive torque on the disk is evaluated and surface velocity profiles are computed for varying values of the ratio of the coefficient of surface shear viscosity to the coefficient of viscosity of the substrate fluid, and the depth of the disk below the surface.

Optimal simultaneous maximum a posterioriestimation of states, noise statistics and parameters I. Algorithm

The simultaneous state and parameter estimation problem for a linear discrete-time system with unknown noise statistics is treated as a large-scale optimization problem. The a posterioriprobability density function is maximized directly with respect to the states and parameters subject to the constraint of the system dynamics. The resulting optimization problem is too large for any of the standard non-linear programming techniques and hence an hierarchical optimization approach is proposed. It turns out that the states can be computed at the first levelfor given noise and system parameters. These, in turn, are to be modified at the second level.The states are to be computed from a large system of linear equations and two solution methods are considered for solving these equations, limiting the horizon to a suitable length. The resulting algorithm is a filter-smoother, suitable for off-line as well as on-line state estimation for given noise and system parameters. The second level problem is split up into two, one for modifying the noise statistics and the other for modifying the system parameters. An adaptive relaxation technique is proposed for modifying the noise statistics and a modified Gauss-Newton technique is used to adjust the system parameters.

Gravity-induced coronal plane joint moments in adolescent idiopathic scoliosis

Background Adolescent Idiopathic Scoliosis is the most common type of spinal deformity, and whilst the risk of progression appears to be biomechanically mediated (larger deformities are more likely to progress), the detailed biomechanical mechanisms driving progression are not well understood. Gravitational forces in the upright position are the primary sustained loads experienced by the spine. In scoliosis they are asymmetrical, generating moments about the spinal joints which may promote asymmetrical growth and deformity progression. Using 3D imaging modalities to estimate segmental torso masses allows the gravitational loading on the scoliotic spine to be determined. The resulting distribution of joint moments aids understanding of the mechanics of scoliosis progression. Methods Existing low-dose CT scans were used to estimate torso segment masses and joint moments for 20 female scoliosis patients. Intervertebral joint moments at each vertebral level were found by summing the moments of each of the torso segment masses above the required joint. Results The patients’ mean age was 15.3 years (SD 2.3; range 11.9 – 22.3 years); mean thoracic major Cobb angle 52° (SD 5.9°; range 42°-63°) and mean weight 57.5 kg (SD 11.5 kg; range 41 – 84.7 kg). Joint moments of up to 7 Nm were estimated at the apical level. No significant correlation was found between the patients’ major Cobb angles and apical joint moments. Conclusions Patients with larger Cobb angles do not necessarily have higher joint moments, and curve shape is an important determinant of joint moment distribution. These findings may help to explain the variations in progression between individual patients. This study suggests that substantial corrective forces are required of either internal instrumentation or orthoses to effectively counter the gravity-induced moments acting to deform the spinal joints of idiopathic scoliosis patients.

Exploring metrics to express energy expenditure of physical activity in youth

Background Several approaches have been used to express energy expenditure in youth, but no consensus exists as to which best normalizes data for the wide range of ages and body sizes across a range of physical activities. This study examined several common metrics for expressing energy expenditure to determine whether one metric can be used for all healthy children. Such a metric could improve our ability to further advance the Compendium of Physical Activities for Youth. Methods A secondary analysis of oxygen uptake (VO2) data obtained from five sites was completed, that included 947 children ages 5 to 18 years, who engaged in 14 different activities. Resting metabolic rate (RMR) was computed based on Schofield Equations [Hum Nutr Clin Nut. 39(Suppl 1), 1985]. Absolute oxygen uptake (ml.min-1), oxygen uptake per kilogram body mass (VO2 in ml.kg-1.min-1), net oxygen uptake (VO2 – resting metabolic rate), allometric scaled oxygen uptake (VO2 in ml.kg-0.75.min-1) and YOUTH-MET (VO2.[resting VO2] -1) were calculated. These metrics were regressed with age, sex, height, and body mass. Results Net and allometric-scaled VO2, and YOUTH-MET were least associated with age, sex and physical characteristics. For moderate-to-vigorous intensity activities, allometric scaling was least related to age and sex. For sedentary and low-intensity activities, YOUTH-MET was least related to age and sex. Conclusions No energy expenditure metric completely eliminated the influence of age, physical characteristics, and sex. The Adult MET consistently overestimated EE. YOUTH-MET was better for expressing energy expenditure for sedentary and light activities, whereas allometric scaling was better for moderate and vigorous intensity activities. From a practical perspective, The YOUTH-MET may be the more feasible metric for improving of the Compendium of Physical Activities for Youth.

The nonequilibrium region of a mixing layer

Experiments were conducted in the nonequilibrium region of a free mixing layer with unequal freestream velocities. Four velocity ratios U(1)/U(2) of 0.32, 0.46, 0.74, and 0.96 were used in this investigation. The growth of the shear layer as well as the velocity adjustment in the near wake were examined. There was reasonable agreement between the measured mean velocity profiles and those computed using the K-epsilon turbulence model. Some periodic turbulence velocity fluctuations were observed in the mixing layer, but their frequency remained the same along the flow.

Design of laterally loaded piles in clays based on cone penetration test data: a reliability-based approach

The behaviour of laterally loaded piles is considerably influenced by the uncertainties in soil properties. Hence probabilistic models for assessment of allowable lateral load are necessary. Cone penetration test (CPT) data are often used to determine soil strength parameters, whereby the allowable lateral load of the pile is computed. In the present study, the maximum lateral displacement and moment of the pile are obtained based on the coefficient of subgrade reaction approach, considering the nonlinear soil behaviour in undrained clay. The coefficient of subgrade reaction is related to the undrained shear strength of soil, which can be obtained from CPT data. The soil medium is modelled as a one-dimensional random field along the depth, and it is described by the standard deviation and scale of fluctuation of the undrained shear strength of soil. Inherent soil variability, measurement uncertainty and transformation uncertainty are taken into consideration. The statistics of maximum lateral deflection and moment are obtained using the first-order, second-moment technique. Hasofer-Lind reliability indices for component and system failure criteria, based on the allowable lateral displacement and moment capacity of the pile section, are evaluated. The geotechnical database from the Konaseema site in India is used as a case example. It is shown that the reliability-based design approach for pile foundations, considering the spatial variability of soil, permits a rational choice of allowable lateral loads.