Steady-state Force Research Articles

Differential effects of myosin activators on myocardial contractile function in non-failing and failing human hearts.

The second-generation myosin activator danicamtiv (DN) has shown improved function compared to the first generation myosin activator omecamtiv mecarbil (OM) in non-failing myocardium by enhancing cardiac force generation but attenuating slowed relaxation. However, whether the functional improvement with DN compared to OM persists in remodeled failing myocardium remain unknown. Therefore, this study aimed to investigate the differential contractile response to myosin activators in non-failing and failing myocardium. Mechanical measurements were performed in detergent-skinned myocardium isolated from donor and failing human hearts. Steady-state force, stretch activation responses, and loaded shortening velocity were analyzed at submaximal [Ca2+] in the absence or presence of 0.5 µmol/L OM or 2 µmol/L DN. The effects of DN and OM on Ca2+-sensitivity of force generation were determined by incubating myocardial preparations at various [Ca2+]. The inherent impairment in force generation and cross-bridge behavior sensitized failing myocardium to the effects of myosin activators. Specifically, increased Ca2+-sensitivity of force generation, slowed rates of cross-bridge recruitment and detachment following acute stretch, slowed loaded shortening velocity, and diminished power output were more prominent following treatment with OM or DN in failing myocardium compared to donor myocardium. Although these effects were less pronounced with DN compared to OM in failing myocardium, DN impaired contractile properties in failing myocardium that were not affected in donor myocardium. Our results indicate that similar to first-generation myosin activators, the DN-induced slowing of cross-bridge kinetics may result in a prolongation of systolic ejection and delayed diastolic relaxation in the heart failure setting.

Simulation analysis and optimization of pilot type electro-hydraulic proportional directional valve based on multi field coupling

The pilot-operated electro-hydraulic proportional directional valve is widely used in complex and precision engineering fields. Its operational stability and reliability determine the overall performance of the hydraulic system. During the working process of the pilot-operated electro-hydraulic proportional directional valve, the spool is deformed due to the influence of pressure and temperature, resulting in problems such as clamping or stagnation. In order to solve this problem, this study uses the fluid-solid-thermal coupling simulation method to analyze the flow field characteristics and thermal deformation of the spool with different structural throttling grooves, and optimizes the U-shaped throttling groove structure. The results indicate that as the valve opening increases, the front section of the throttling groove becomes the primary region where the maximum radial deformation of the spool occurs. In this paper, the steady-state flow force of the valve core is taken as the measurement standard, the constraint range of the structural parameters is determined, and the radial thermal deformation is included in the optimization analysis of the structural size. It’s verified that the optimized flow characteristics are improved by about 32 % compared with that before optimization, and the maximum steady-state flow force after considering thermal deformation is reduced by about 4N compared with that before optimization. This provides an important design basis and reference for the structural design and optimization of the pilot electro-hydraulic proportional directional valve.

Change in EEG-EMG synchronization reflecting abnormal functional corticomuscular coupling following stroke: A pilot study

Stroke is a common clinical cardiovascular and cerebrovascular disease, which can cause severe motor dysfunction. Therefore, it is an important topic to investigate the abnormality mechanism of the cerebral cortex and the corresponding muscles after stroke. In this study, we investigated the functional corticomuscular coupling (FCMC) at specific frequencies by analyzing differences between stroke patients and healthy controls in hand movements. The transfer spectral entropy (TSE) method was used to analyze simultaneous electroencephalography (EEG) and electromyography (EMG) in the right-hand steady state force task. The results illustrated that healthy subjects had the highest TSE values at the beta band in the EMG→EEG and EEG→EMG directions, and the TSE value in the EEG→EMG direction was higher than that in the EMG→EEG direction. In contrast, for stroke patients, beta band coupling was weakened, and there was a notably higher enhancement of alpha and gamma bands in the EMG→EEG direction relative to the EEG→EMG direction. Further analysis found significant correlations between TSE area values at beta2 and gamma2 bands and clinical rating scales. This study demonstrates the frequency specificity properties of FCMC estimated by TSE can assess the rehabilitation status of stroke patients and contribute to our comprehension of the potential mechanism of motor control systems.

Phenomenological Muscle Constitutive Model With Actin-Titin Binding for Simulating Active Stretching.

The force produced by a muscle depends on its contractile history, yet human movement simulations typically employ muscle models that define the force-length relationship from measurements of fiber force during isometric contractions. In these muscle models, the total force-length curve can have a negative slope at fiber lengths greater than the fiber length at which peak isometric force is produced. This region of negative stiffness can cause numerical instability in simulations. Experiments have found that the steady-state force in a muscle fiber following active stretching is greater than the force produced during a purely isometric contraction. This behavior is called residual force enhancement. We present a constitutive model that exhibits force enhancement, implemented as a hyperelastic material in the febio finite element software. There is no consensus on the mechanisms responsible for force enhancement; we adopt the assumption that the passive fiber force depends on the sarcomere length at the instant that the muscle is activated above a threshold. We demonstrate the numerical stability of our model using an eigenvalue analysis and by simulating a muscle whose fibers are of different lengths. We then use a three-dimensional muscle geometry to verify the effect of force enhancement on the development of stress and the distribution of fiber lengths. Our proposed muscle material model is one of the few models available that exhibits force enhancement and is suitable for simulations of active lengthening. We provide our implementation in febio so that others can reproduce and extend our results.

90-degree peeling of elastic thin films from elastic soft substrates: Theoretical solutions and experimental verification

Peeling of thin films has been widely used in adhesion measurement, film transfer and bio-inspired design. Most previous studies focused on the peeling of thin films from rigid substrates, but soft substrates are common in practical applications. Herein, we propose a two-dimensional model based on the bilinear cohesive law to characterize the 90-degree peeling of elastic thin films from elastic soft substrates, and obtain theoretical solutions expressed in terms of the Chebyshev series. The theoretical solutions match well with the finite element method results, including the load-displacement curves and the bulging deformation of soft substrates. We find that with decreasing substrate modulus, the maximum peeling force (Pmax) decreases but the steady-state peeling force remains unchanged. With the present solutions, the interfacial strength and fracture energy can be extracted simultaneously from the 90-degree peeling experiments of thin film/soft substrate systems, and then the experimentally measured Pmax for different film thicknesses can be well predicted. Furthermore, we obtain a new power scaling law of Pmax, where the scaling exponent depends on substrate elasticity. These results can help us measure the interfacial properties of thin film/soft substrate systems via peel tests, and regulate their peeling behaviors by interface design.

Impairment of Multi-Finger Force Control in Adult Patients with Dyskinetic Cerebral Palsy

PURPOSE: This study aimed to investigate the impact of cerebral palsy (CP) on motor control strategies, with a focus on understanding the altered capability of multi-finger force control in adult patients with dyskinetic CP.METHODS: Eight adults with dyskinetic CP and ten age- and sex-matched healthy controls participated in this study. The participants performed three force production tasks using four fingers of the dominant hand under isometric conditions: maximum voluntary contraction (MVC) to assess maximum force production, a ramp task to evaluate finger interdependency (enslaving), and a pulse task to examine the ability to rapidly and precisely adjust force. The co-contraction index (CCI) of finger flexor and extensor was also determined during steady state and pulse force production.RESULTS: Patients with dyskinetic CP demonstrated significant impairments in finger force control compared to healthy controls. Specifically, the accuracy and consistency of force production were reduced in the CP group, and they exhibited a longer time to peak pulse force. The CP group also showed increased finger interdependency and elevated CCI during pulse force production, suggesting a greater reliance on co-contraction and less efficient motor control strategies.CONCLUSIONS: This study highlights the distinct motor control deficits in individuals with dyskinetic CP, particularly in tasks requiring rapid and precise finger force adjustment. These results have important implications for the functional rehabilitation of patients with CP. Therapies aimed at reducing excessive co-contraction and decreasing finger interdependency may be beneficial for improving muscle coordination and overall motor function.

Multiscale modeling shows how 2’-deoxy-ATP rescues ventricular function in heart failure

2'-deoxy-ATP (dATP) improves cardiac function by increasing the rate of crossbridge cycling and Ca[Formula: see text] transient decay. However, the mechanisms of these effects and how therapeutic responses to dATP are achieved when dATP is only a small fraction of the total ATP pool remain poorly understood. Here, we used a multiscale computational modeling approach to analyze the mechanisms by which dATP improves ventricular function. We integrated atomistic simulations of prepowerstroke myosin and actomyosin association, filament-scale Markov state modeling of sarcomere mechanics, cell-scale analysis of myocyte Ca[Formula: see text] dynamics and contraction, organ-scale modeling of biventricular mechanoenergetics, and systems level modeling of circulatory dynamics. Molecular and Brownian dynamics simulations showed that dATP increases the actomyosin association rate by 1.9 fold. Markov state models predicted that dATP increases the pool of myosin heads available for crossbridge cycling, increasing steady-state force development at low dATP fractions by 1.3 fold due to mechanosensing and nearest-neighbor cooperativity. This was found to be the dominant mechanism by which small amounts of dATP can improve contractile function at myofilament to organ scales. Together with faster myocyte Ca[Formula: see text] handling, this led to improved ventricular contractility, especially in a failing heart model in which dATP increased ejection fraction by 16% and the energy efficiency of cardiac contraction by 1%. This work represents a complete multiscale model analysis of a small molecule myosin modulator from single molecule to organ system biophysics and elucidates how the molecular mechanisms of dATP may improve cardiovascular function in heart failure with reduced ejection fraction.

AcousTac: Tactile Sensing with Acoustic Resonance for Electronics-Free Soft Skin.

Sound is a rich information medium that transmits through air; people communicate through speech and can even discern material through tapping and listening. To capture frequencies in the human hearing range, commercial microphones typically have a sampling rate of over 40 kHz. These accessible acoustic technologies are not yet widely adopted for the explicit purpose of giving robots a sense of touch. Some researchers have used sound to sense tactile information, both monitoring ambient soundscape and with embedded speakers and microphones to measure sounds within structures. However, these options commonly do not provide a direct measure of steady state force or require electronics integrated somewhere near the contact location. In this work, we present AcousTac, an acoustic tactile sensor for electronics-free, force-sensitive soft skin. Compliant silicone caps and plastic tubes compose the resonant chambers that emit pneumatic-driven sound measurable with a conventional off-board microphone. The resulting frequency changes depend on the external loads on the compliant endcaps. The compliant cap vibrates with the resonant pressure waves and is a nonidealized boundary condition, initially producing a nonmonotonic force response. We characterize two solutions-adding a distal hole and mass to the cap-resulting in monotonic and nonhysteretic force readings with this technology. We can tune each AcousTac taxel to specific force and frequency ranges, based on geometric parameters including tube length, and thus uniquely sense each taxel simultaneously in an array. We demonstrate AcousTac's functionality on two robotic systems: a 4-taxel array and a 3-taxel astrictive gripper. Simple to implement with off-the-shelf parts, AcousTac is a promising concept for force sensing on soft robotic surfaces, especially in situations where electronics near the contact are not suitable. Equipping robots with tactile sensing and soft skin provides them with a sense of touch and the ability to safely interact with their surroundings.

Modeling Torsional Stiffness and Damping of Liquid Slosh in Spacecraft Diaphragm Tanks

This study presents a pendulum-based empirical model of liquid slosh in spacecraft diaphragm tanks. Diaphragm torsional stiffness and effective damping are predicted within a range of lateral excitation frequency typical of launch or on-orbit operations. A pendulum-analog model with two primary slosh modes is extracted from diaphragm tank slosh data corresponding to steady-state triaxial force measurements at the support points of the tank under sinusoidal excitation. The approach leads to repeatable and accurate determination of pendulum model parameters. A gradient-based algorithm and a procedure for setting initial conditions for optimal parameter estimation are proposed to predict the liquid slosh forces and torques acting on the tank as frequency response functions of apparent mass and apparent mass moment. The pendulum models were assessed by comparing model-predicted force/torques with measurements, showing good agreement (1–5% prediction error in the time domain) under most test conditions (fill level, excitation amplitude, and frequency). An empirical model of the stiffness and damping is also proposed for the tested range of operational conditions. Because of the nonlinear nature of liquid slosh behavior in a diaphragm tank, both the stiffness and damping are local values for a given set of test conditions.

Determination traction-separation parameters and its sensitivity analysis of bilinear cohesive zone model through experimental and numerical peeling analysis of electrodes

The optimization of traction-separation parameters within the cohesive zone model (CZM) is crucial for accurately simulating electrode peeling failure in lithium-ion batteries. This study performed 180° peeling tests on NCM electrodes and used Abaqus software alongside the sim-flow optimization tool for simulation analysis, aiming to accurately replicate the failure behavior at the lithium-ion battery electrode interface. Experimental results were used to determine the fracture energy of the electrode interface, while numerical simulations analyzed the peeling performance of the electrode. By fitting the force-displacement curves from both simulation and experiment, the traction-separation parameters within the bilinear CZM were optimized. Additionally, the study analyzed the effects of traction-separation parameters, as well as the modulus and thickness of the current collector, on peeling behavior. The results showed that a significant increase in the interface stiffness and strength increased the steady-state peeling force. Among these, the fracture energy had the most significant influence on the steady-state peeling force. Furthermore, while the modulus and thickness of the current collector affected the maximum peeling force, their effect on the steady-state peeling force was negligible. These findings enhance the accuracy of simulating the failure behavior of the electrode interface.

Dynamic biaxial loading of vascular smooth muscle cell seeded tissue equivalents

An intricate reciprocal relationship exists between adherent synthetic cells and their extracellular matrix (ECM). These cells deposit, organize, and degrade the ECM, which in turn influences cell phenotype via responses that include sensitivity to changes in the mechanical state that arises from changes in external loading. Collagen-based tissue equivalents are commonly used as simple but revealing model systems to study cell–matrix interactions. Nevertheless, few quantitative studies report changes in the forces that the cells establish and maintain in such gels under dynamic loading. Moreover, most prior studies have been limited to uniaxial experiments despite many soft tissues, including arteries, experiencing multiaxial loading in vivo. To begin to close this gap, we use a custom biaxial bioreactor to subject collagen gels seeded with primary aortic smooth muscle cells to different biaxial loading conditions. These conditions include cyclic loading with different amplitudes as well as different mechanical constraints at the boundaries of a cruciform sample. Irrespective of loading amplitude and boundary condition, similar mean steady-state biaxial forces emerged across all tests. Additionally, stiffness-force relationships assessed via intermittent equibiaxial force–extension tests showed remarkable similarity for ranges of forces to which the cells adapted during periods of cyclic loading. Taken together, these findings are consistent with a load-mediated homeostatic response by vascular smooth muscle cells.

Danicamtiv affected isometric force and cross-bridge kinetics similarly in skinned myocardial strips from male and female rats.

Myotropes are pharmaceuticals that have recently been developed or are under investigation for the treatment of heart diseases. Myotropes have had varied success in clinical trials. Initial research into myotropes have widely focused on animal models of cardiac dysfunction in comparison with normal animal cardiac physiology-primarily using males. In this study we examined the effect of danicamtiv, which is one type of myotrope within the class of myosin activators, on contractile function in permeabilized (skinned) myocardial strips from male and female Sprague-Dawley rats. We found that danicamtiv increased steady-state isometric force production at sub-maximal calcium levels, leading to greater Ca2+-sensitivity of contraction for both sexes. Danicamtiv did not affect maximal Ca2+-activated force for either sex. Sinusoidal length-perturbation analysis was used to assess viscoelastic myocardial stiffness and cross-bridge cycling kinetics. Data from these measurements did not vary with sex, and the data suggest that danicamtiv slows cross-bridge cycling kinetics. These findings imply that danicamtiv increases force production via increasing cross-bridge contributions to activation of contraction, especially at sub-maximal Ca2+-activation. The inclusion of both sexes in animal models during the formative stages of drug development could be helpful for understanding the efficacy or limitation of a drug's therapeutic impact on cardiac function.

Analysis of Nonlinear Vibration Characteristics and Whirl Behavior of Dual-Rotor Systems with Inter-Shaft Rub Impact

Previous studies on rub-impact faults have mainly focused on the rub-impact between rotors and stators, with less research on inter-rotor rub impact. The impact of inter-rotor rub impact on rotor nonlinear vibration is particularly significant. This study investigates the effects of inter-shaft rub impact on the vibration characteristics and whirl behavior of dual-rotor systems. Initially, a dual-rotor model with inter-shaft bearings is established using the finite element method, and inter-shaft rub-impact forces are derived based on contact mechanics. Next, the system response is solved using the Newmark method. Vibration characteristics are analyzed through Campbell diagrams, 3D waterfall plots, time-frequency domain plots, and steady-state rub-impact force plots. Finally, the influence of inter-shaft rub impact on the whirl behavior of the dual-rotor system is studied based on the theory of full-spectrum analysis. The study concludes that inter-shaft rub-impact faults shift the system’s resonance points backward, increase harmonic and combination frequency components, and significantly affect the system response under dual-rotor co-rotation. Excessive friction can lead to self-excited vibrations and sudden amplitude increases, particularly in the LP rotor frequency. Additionally, inter-shaft rub impact primarily affects the whirl behavior of the LP-compressor disk1, showing multiple cycles of forward and backward whirl alternation during acceleration due to combined unbalanced and rub-impact excitations.

Addressing Disparities Between Physician Rated and Self-Reported Dyspnea Using Machine Learning

Nearly 50% of intensive care and 25% of tertiary care patients experience dyspnea, or the sensation of unsatisfied breathing. Dyspnea decreases quality of life, impacts clinical outcomes, and can have long standing psychological consequences for some patients. Current methods of identifying dyspnea severity (such as the respiratory distress observation scale or RDOS) lack suffcient quantifiable measurements of dyspnea to accurately predict its severity. In our testing, while the RDOS does have predictive power, it frequently produces underestimates compared to self-reported dyspnea severity (r = 0.32; p = 0.0081). Furthermore, dyspnea severity can change throughout the day, requiring healthcare professionals to conduct RDOS evaluations multiple times to attempt to meet patient needs. Therefore, there is an urgent need for continuous, accurate monitoring of dyspnea severity in the hospital setting. To address this need, we are using non-invasive data collected during induced bouts of dyspnea to train a machine learning (ML) algorithm in dyspnea prediction. We employ a steady-state end-tidal forcing breathing system and adjust end-tidal gas tensions manually while collecting PO­­­­2, PCO2, SpO2, and ventilation data. Preliminary findings show several physiological and demographic factors independently correlate to dyspnea severity, including BMI (r = 0.49; p = 0.007), end-tidal PCO­­2 (r = 0.33; p = 0.008), and inspired PO2 (r = -0.25; p = 0.043). We achieved an accuracy of 84%, 83.3%, and 82% in-dataset predictions using k-nearest neighbor, random forest, and logistic regression approaches respectively. Several physiological biomarkers and demographic factors were used as inputs and self-reported dyspnea was used as the predicted output. These findings show promise for the use of ML in dyspnea prediction. With further testing and use of a broader participant population, it may be possible to use ML models such as these to improve patient outcomes and overall quality of life. University of California, Riverside. This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Analysis of Wind Field Response Characteristics of Tethered Balloon Systems

Tethered balloon systems encounter various complex wind field environments during flight. To investigate the conditions under which the system can operate safely and smoothly, a longitudinal dynamic model for tethered balloon systems is established. The model incorporates a streamlined balloon shape with its aerodynamic center at the body’s center. Steady-state aerodynamic force coefficients are calculated through simulations and fitted to a function based on the angle of attack within a specified range. The complex cable model is simplified using the lumped mass method, considering the influence of branch cables on the main node position. Experimental results from windless oscillation tests on scaled tethered balloon systems are compared with numerical solutions obtained using the dynamic model under the same conditions, validating the feasibility of the model for simulating different wind field scenarios. Finally, the motion characteristics of tethered balloon systems in different wind fields are analyzed. The numerical simulation results show that in a horizontal step wind field, the cable tension and cable inclination angle increase with the wind speed, and the slower the wind field changes, the shorter the time required for system stabilization. Updrafts greatly increase the likelihood of balloon escape, while downdrafts greatly increase the likelihood of the system making contact with the ground. The findings of this study can provide a basis for selecting suitable wind field conditions and issuing risk warnings for tethered balloon systems.

Numerical Simulations of Ionic Wind Induced by Positive DC-Corona Discharges

This paper analyzes ionic wind production and propulsive force in various electrode configurations under atmospheric conditions. By considering the aerodynamic forces in addition to previously considered electric ones, new predictions for steady-state forces and ionic wind flow velocity are successfully compared with experimental measurements, providing convincing quantitative evidence of the predictive capabilities of drift-diffusion modeling associated with one-way Coulomb forcing of Navier–Stokes equations for ionic wind generation. Furthermore, various electrode configurations are analyzed, some of them streamlined, reducing wakes downstream collectors on the one hand and providing additional thrust on the other. The quantification of these additional thrusts is analyzed, physically discussed, and explored in various configurations.

Research on vibration characteristics of cone valve based on steady state hydrodynamics

The steady-state fluid dynamic force of the cone valve core reduces the stability of the cone valve, thereby affecting its vibration characteristics. This study first derived the formula for calculating the steady-state fluid dynamic force of the valve core using the momentum theorem. Then, the study used the Computational Fluid Dynamics (CFD) method to explore the variation of steady-state fluid dynamic force with the opening degree of the cone valve and the pressure difference at the inlet and outlet. The study also analyzed the steady-state fluid dynamic force of the improved structure cone valve. Combining the CFD calculation results with the established axial vibration and dynamic models of the valve core spring system, a system dynamics simulation model was developed. Innovatively, the study combined experimental and simulation methods to investigate the vibration characteristics of the cone valve and the impact of steady-state fluid dynamic force on these characteristics. The experimental results aligned well with the simulation results, showing that the influence of steady-state fluid dynamic force on vibration characteristics is similar to that of spring force. At a small opening degree, the compensation structure of steady-state fluid dynamic force can reduce axial vibration amplitude. Increasing the inlet pressure from 2.0 MPa to 2.5 MPa results in a decrease in axial amplitude and an increase in the numerical value of the vibration balance position. Increasing the outlet pressure from 0.5 MPa to 0.9 MPa initially decreases axial amplitude and axial waveform factor, followed by an increase. Increasing the spring preload from 9.4 mm to 13.4 mm leads to a decrease in the numerical value of the balance position, a reduction in amplitude, and a decrease in the waveform factor of the axial vibration of the valve core. The study reveals systematic effects of inlet pressure, outlet pressure, and spring preload on the vibration characteristics of the cone valve. The research results in this paper can provide theoretical support for the structural design of hydraulic valves, reduce the impact of fluid dynamics on vibration, and thereby improve the stability of hydraulic system operation.

Design and optimization of hydraulic slide valve spool structure based on steady state flow force

This article studies the structure of hydraulic spool valve and the influence of different pressure differences on steady-state flow force. An improved scheme of designing a jet guide groove on the valve core is proposed to compensate for steady-state flow force. This article uses the CFD method to study the influence of spool valve structural parameters on steady-state flow force. The response surface method is used to establish the Polynomial regression model of the steady-state flow force, and the optimal structure size of the jet guide slot is obtained. The results indicate that the flow force increases with the increase of pressure difference. Moreover, under the same operating conditions, the absolute flow force of the spool valve with jet guide groove is smaller than that of the spool valve without jet guide groove. Therefore, the jet guide groove can effectively compensate for steady-state flow force, improve the control accuracy of the valve, and improve the stability of the spool valve.

Orbital evolution of eccentric perturbers under dynamical friction: crossing the sound barrier

ABSTRACT In a gaseous medium, dynamical friction (DF) reaches a maximum when the orbital speed of a (point-like) perturber moving on a circular orbit is close to the sound speed. Therefore, in a quasi-steady state, eccentric orbits of perturbers approaching the sound barrier (from below) should rapidly circularize as they experience the strongest drag at pericentre passage. To investigate this effect, we extend the solution for circular DF in a uniform gaseous medium to eccentric Keplerian orbits. We derive an approximation to the steady-state DF force, which is valid for eccentricities as high as e = 0.9 in a limited range of Mach number around the transition to supersonic regime. We validate our analytical result with 3D simulations of the gas density response. Although gaseous DF generally dissipates orbital energy, we find that it can be directed along the motion of the perturber near pericentre passage when the eccentricity is e ≳ 0.9. We apply our results to compute the long-time evolution of the orbital parameters. Most trajectories tend to circularize as the perturber moves into the supersonic regime. However, orbits with eccentricities e ≳ 0.8 below the sound barrier experience a slight increase in eccentricity as they loose orbital energy. Possible extensions to our analytical approach are also discussed.

Residual force depression is not related to positive muscle fascicle work during submaximal voluntary dorsiflexion contractions in humans.

Residual force depression (rFD) following active muscle shortening is assumed to correlate most strongly with muscle work, but this has not been tested during voluntary contractions in humans. Using dynamometry, we compared steady-state ankle joint torques (N=16) following tibialis anterior (TA) muscle-tendon unit (MTU) lengthening and shortening to the time-matched torque during submaximal voluntary fixed-end dorsiflexion reference contractions (REF) at a matched MTU length and EMG amplitude. Ultrasound revealed significantly reduced (P<0.001) TA fascicle shortening amplitudes during MTU lengthening without a preload over small and medium amplitudes, respectively, relative to REF. MTU lengthening with a preload over a large amplitude significantly (P<0.001) increased fascicle shortening relative to REF, as well as stretch amplitudes relative to MTU lengthening without a preload (P=0.001). Significant (P=0.028) steady-state fascicle force enhancement relative to REF was observed following MTU lengthening, and was similar among MTU lengthening-hold conditions (3-5%). MTU shortening with and without a preload over small and large amplitudes significantly (P<0.001) increased positive fascicle and MTU work relative to REF, but significant (P=0.006) rFD was observed following MTU shortening with a preload (7-10%) only. rFD was linearly related to positive MTU work [rrm (47)=0.48, P<0.001], but not positive fascicle work [rrm (47)=0.16, P=0.277]. Our findings indicate that MTU lengthening without substantial fascicle stretch enhances steady-state force output, which might arise from less shortening-induced rFD. Our findings also indicate similar rFD following different amounts of positive fascicle/MTU work, which cautions against using work to predict rFD during submaximal voluntary contractions. KEY POINTS: Accurately predicting muscle force is challenging because active muscle shortening depresses force output. The residual force depression (rFD) that exists following active muscle shortening is commonly assumed to correlate strongly and positively with muscle work. We found that tibialis anterior muscle fascicle work and muscle-tendon unit work did not accurately predict rFD during submaximal voluntary dorsiflexion contractions. Fascicle shortening during fixed-end reference contractions also potentially induced rFD of 3-5%, which was similar to the rFD following muscle-tendon unit shortening without a preload. A higher number of active muscle fibres during shortening probably increased rFD, which suggests that motor unit recruitment during shortening might predict rFD.