Next Article in Journal
Monitoring the Early Strength Development of Cement Mortar with Piezoelectric Transducers Based on Eigenfrequency Analysis Method
Next Article in Special Issue
Assessing Electroencephalography as a Stress Indicator: A VR High-Altitude Scenario Monitored through EEG and ECG
Previous Article in Journal
Pendulum-Based River Current Energy Converter for Hydrometric Monitoring Systems
Previous Article in Special Issue
Automatic Cognitive Fatigue Detection Using Wearable fNIRS and Machine Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Respiratory Inductance Plethysmography to Assess Fatigability during Repetitive Work

1
Laboratório de Instrumentação, Engenharia Biomédica e Física da Radiação (LIBPhys-UNL), Departamento de Física, Faculdade de Ciências e Tecnologia, Universidade Nova de Lisboa, 2829-516 Caparica, Portugal
2
Associação Fraunhofer Portugal Research, Rua Alfredo Allen 455/461, 4200-135 Porto, Portugal
3
Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
4
MIT.nano Immersion Lab, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
5
Device Realization Laboratory, Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
6
Faculty of Medicine, Rīga Stradiņš University, 16 Dzirciema iela, LV-1007 Rīga, Latvia
7
Campus de Azurém, Minho University, 4800-058 Guimarães, Portugal
8
Brigham and Women’s Hospital, Cardiovascular Division, 75 Francis Street, Boston, MA 02115, USA
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(11), 4247; https://doi.org/10.3390/s22114247
Submission received: 1 May 2022 / Revised: 27 May 2022 / Accepted: 30 May 2022 / Published: 2 June 2022
(This article belongs to the Special Issue Wearable and Mobile Sensors and Data Processing)

Abstract

:
Cumulative fatigue during repetitive work is associated with occupational risk and productivity reduction. Usually, subjective measures or muscle activity are used for a cumulative evaluation; however, Industry 4.0 wearables allow overcoming the challenges observed in those methods. Thus, the aim of this study is to analyze alterations in respiratory inductance plethysmography (RIP) to measure the asynchrony between thorax and abdomen walls during repetitive work and its relationship with local fatigue. A total of 22 healthy participants (age: 27.0 ± 8.3 yrs; height: 1.72 ± 0.09 m; mass: 63.4 ± 12.9 kg) were recruited to perform a task that includes grabbing, moving, and placing a box in an upper and lower shelf. This task was repeated for 10 min in three trials with a fatigue protocol between them. Significant main effects were found from Baseline trial to the Fatigue trials (p < 0.001) for both RIP correlation and phase synchrony. Similar results were found for the activation amplitude of agonist muscle (p < 0.001), and to the muscle acting mainly as a joint stabilizer (p < 0.001). The latter showed a significant effect in predicting both RIP correlation and phase synchronization. Both RIP correlation and phase synchronization can be used for an overall fatigue assessment during repetitive work.

1. Introduction

Fatigue-inducing repetitive work is a widely known risk factor for occupational-related upper-limb musculoskeletal disorders [1,2,3], often requiring a consolidation of mental and physical demands [4]. Fatigue may interfere with the work itself by decreasing efficiency, and if it is excessive and extended, it may alter the worker’s mental and physical conditions, giving rise to illness or worsen a preexisting condition [5,6], even in tasks requiring high vigilance but low neuromuscular work [4]. This is a concern if we take into consideration that people worked at their head job on average 37.1 h a week in Europe [7] and 34.9 h in the United States [8]. Additionally, approximately one-third of the global workforce works more than 48 h per week [9]. Furthermore, work-related fatigue can bring consequences beyond the occupational risk and work-related illness, it fosters the occurrence of mistakes due to tiredness, which has severe implications in productivity and promotes absenteeism [10,11].
To avoid workers’ fatigue, ergonomic interventions such as a job rotation plan and breaks have been suggested, but the former do not demonstrate a significant effect on fatigue at lower exertion tasks [3]. On the other hand, the implementation of breaks requires fixed work/rest ratio or workload variation based on postures and movement.
The assessment of occupational risk can be performed with several methods classified according to three main approaches: self-reported, observational, and directly measured. Self-reported methods rely on questionnaires based on physical demand perception. However, they require interruptions of the worker’s activity and are prone to self-reporting bias arising from external perceived factors (e.g., behavioral and physiological factors). Observational methods are conducted by an external observer watching a worker performing his job and then filling appropriate standard instruments and ergonomic risk assessment tools (e.g., [12,13,14]). Direct methods use sensors that give physiological, movement and position, and/or environmental information (e.g., [15,16,17]). The latter has been gaining attention over the last years with the technological developments in wearables and data analytics surrounding the industry Operator 4.0 [18].
Regarding fatigue evaluation, empirical and theoretical models have been introduced, but they require individual calibration through fitting experimental data for a specific task [19]. Fatigue evaluation using wearables has recently attracted increasing attention as a promising approach to overcome the challenges associated with previous methods [20,21]. Wearables allow ubiquitous data collection in the workplace and combined with data analytics can provide useful information to enable healthy operators 4.0 [22,23]. In this sense, fatigue at low exertion load is usually evaluated by electromyography (EMG) [24]; nevertheless, also respiration biomarkers have been proposed based on the rib cage and abdominal asynchrony [25] and hyperventilation hypotheses [26]. These biomarkers can give information on overall tiredness rather than local, or even the effects of local fatigue on general fatigability and tiredness. We further provide an overview of such methods and identify their limitations.
EMG fatigue manifestation is observed by a decrease in signal median frequency and an increase in amplitude observed over time [24,27]. Compared to activities performed at high-level intensity, the quantification of fatigue due to low-level force is not so evident. However, fatigue can be used as a biomarker for cumulative exposure to repetitive work [28]. In amplitude analysis, standardized voluntary contractions have frequently been used, in which maximal voluntary contraction (MVC) tends to be the preferred reference choice (e.g., [29,30]). In the review of de Looze et al. [24], the MVC values that quantified a low load characterizing repetitive work ranged between 3 up to the cut-point of 20% of MVC in trapezius muscle, in which fatigue was found to be prevalent between 15 and 20% of MVC with a task duration varying between 1 h up to 9.5 h. Due to the low exertion during long periods associated with work unlike sports, amplitude probability distribution function [31] is commonly used to quantify local muscular workload [32,33]. The median frequency, when showing alterations, decreases ranging between 2 and 15% per hour [34]. However, in workplaces where the effort is almost at a lower level, this indicator could be insufficient or dependent on the posture fluctuations [35]. The EMG, regardless its relevance, presents several drawbacks, such as: (1) can only provide local muscle load, lacking on overall physical load [36]; and (2) it requires MVC procedure, which is difficult to adapt to the non-laboratory environment of industrial assembly lines settings.
To overcome overall physical load, the relationship between reserve oxygen consumption and heart rate (HR) has been a suitable option [37,38] since the former is the product between the arteriovenous difference and cardiac output (the product between HR and stroke volume). An indicator that works as an index of the intensity of physical work is the relative HR (e.g., [39,40]), i.e., the ratio of the difference between work and rest HR with the reserve HR (maximum HR–rest HR). From relative HR, the maximum acceptable work time (MAWT) balance between cardiovascular load and physical workload can be determined. Both relative HR and MAWT indicators were associated with the percentage of maximal oxygen consumption and work-related oxygen consumption [41].
Considering the previous paragraph, it is reasonable to wonder whether MAWT can be associated with other respiratory concepts such as asynchrony and hyperventilation theories. The former states that, during repetitive work, there is an asynchrony between the work rhythm and a worker’s internal physiological rhythms which is a source of both motor and cognitive stress [25]. The alterations in respiratory function can be either derived from intense or continued physical exertion, or from disease. Thoraco-abdominal (a)synchrony can be used to determine respiratory disorders and/or respiratory muscle dysfunctions associated with respiratory distress and increased breathing load [42]. The (a)synchrony mechanism depends on the heavy-handed activity of the respiratory muscles, such as inspiratory accessory muscles and the diaphragm, which bring forth changes in lung volume [43]. On the other hand, hyperventilation or over-breathing occurs when there is a drop in arterial carbon dioxide, exceeding metabolic demands of oxygen due to ventilation, which in turn causes a rise in blood pH, that is, respiratory alkalosis [26]. Additionally, hyperventilation is usually described by a shift from a diaphragmatic to a thoracic breathing pattern, which imposes biomechanical musculoskeletal stress on the neck/shoulder region, namely in the muscles that support thoracic breathing [44].
Several parameters have shown valuable use in discriminating healthy respiratory patterns from those with neuromuscular diseases, such as phase angle, percent of rib cage contribution to breathing, respiratory rate, and labored breathing index [45]. Therefore, it was demonstrated that the relationship between respiratory function and fatigability can be explored as an occupational risk indicator. Nevertheless, it remains unclear if internal respiratory synchrony is altered due to fatigue. Furthermore, there is a lack of literature on evaluating the feasibility of using respiration to monitor repetitive work from industrial operators. On the other hand, it has been shown that physiological and subjective fatigue outcomes are not always consistently correlated due to the large variability of requirements from different work settings [46]. Accordingly, this study aims to verify thoraco-abdominal (a)synchrony relationship with fatigability while performing simulated repetitive work. In addition, this work also analyzes the relationship between thoraco-abdominal (a)synchrony and local muscle load.
The rest of this work is organized as follows: Section 2 outlines the study methods; the results are presented in Section 3 and discussed in Section 4; finally, Section 5 summarizes the contributions of this work, along with some recommendations for future work.

2. Materials and Methods

To quantify the problem under study, i.e., the susceptibility of using the synchronization between chest and abdominal wall movements as an indicator of fatigability/fatigue, respiratory inductance plethysmography (RIP) was used. To supplement this assessment with local fatigue and alterations that can occur in the movement of the participants, EMG and accelerometry data were also collected. This section begins by presenting the characterization of the voluntaries that participated in this study. Then, the task that was used to simulate repetitive work is explained in detail. Next, the instrumentation of the participants, equipment used, and steps in processing the different signals are provided. Finally, the statistical analysis used in the study is described.

2.1. Participants

A total of 22 healthy participants (11 females and 11 males; age: 27.0 ± 8.3; height: 1.72 ± 0.09 m; mass: 63.4 ± 12.9 kg) were recruited for this study. The protocol was clearly explained, and all subjects gave their informed consent for inclusion before they participated in the study. The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of University of Porto.

2.2. Task

The protocol included three trials of 10 min duration each (Baseline, Fatigue 1, and Fatigue 2). Between trials, a fatigue protocol was performed to ensure a loss of muscle force during the second and third tasks. Data from two minutes of rest before and after the three 10-min trials were also acquired. The repetitive task performed during each 10-min trial was chosen to contain basic motions required to perform work, such as reaching, grasping, releasing, moving, and positioning an object. The task starts from the neutral pose in the upright position, grabbing a box from a lower shelf of a locker. Then, the participant moves the box, placing it on the upper shelf. Immediately after returning to the neutral pose, the participant grabbed the box on the upper shelf, moving it to the lower shelf (see Figure 1). This description corresponds to one cycle with an ascending phase (moving the box from the lower shelf to the upper one) and a descending phase (moving the box from the upper shelf to the lower one). The periodicity among cycles was self-paced. The shelves’ height was adjusted according to the subject’s morphology: the lower shelf was placed at the height of the medium finger while the participant was standing in the upright position; the upper shelf was adjusted so that the box was placed at the eyes level. The box was set at 5% of the participant’s mass.
The fatigue protocol comprised three trials of 10–15 repetitions each until achieving muscular failure. The resistance was set at 65–75% of maximal force, starting with an estimated load according to participants body mass. To ensure that muscle failure occurred in the interval of 10 to 15 repetitions, the following procedure was performed: when the participant reached muscle failure before being able to complete 10 repetitions, the load was decreased according to the number of repetitions performed, and this attempt had to be repeated; when the participant exceeded 15 repetitions, the load was increased so that muscle failure could occur in the mentioned interval, repeating the attempt. During the last repetition of each trial, the participants were encouraged to perform one last repetition with the researcher’s help.

2.3. Respiratory Inductance Plethysmography

RIP is a non-invasive method to assess lung volume using elastic bands around the thorax and abdomen [45,47]. Transducers (Plux, Lisbon, Portugal) were placed at the rib cage (at the level of the nipples in males, and over the chest passing underarms in females) and at the umbilicus to monitor changes in trunk volume during inspiration and expiration cycles [48].
The raw signal was sampled at a rate of 1000 Hz (this sample rate was used due to EMG synchronization) and Finite Impulse Response bandpass filtered between 0.15 and 0.45 Hz [49]. The signal was then downsampled to 40 Hz and submitted to Masked Sift Empirical Mode Decomposition [50,51,52] to extract the first intrinsic mode function. Finally, the signal was z-normalized and a smoothing procedure with a sliding average of 10 s was applied. Rolling window correlation and instantaneous phase synchrony between chest and abdominal respiration signals were computed.
The instantaneous phase synchrony measures the similarities regarding the phase (angles) between two signals at each point. Its calculation includes the following steps: (i) computation of discrete Hilbert transform of each signal (summarized in Equation (1)) [53], (ii) extraction of the phase of each signal (imaginary part of the transform, represented in Equation (2)), and (iii) calculation of the difference between the two phases (Equation (3)). In Figure 2, three different examples of respiration signals are presented: the first where the subject is in a resting state, the second during the first trial and the third and last during the last trial.
g ( n ) = { 2 π   t   o d d f ( t ) n t n   e v e n 2 π   t   e v e n f ( t ) n t n   o d d
= i m a g ( g ( n ) )
φ ( n ) = 1 sin ( a b s ( c h e s t ( n ) a b d o m i n a l ( n ) / 2 ) )

2.4. EMG

Surface EMG system (Plux, Lisbon, Portugal) with disposable adhesive Ag/AgCl electrodes were used to acquire bilaterally muscle activity from the upper trapezius, anterior deltoideus, biceps brachii, and triceps brachii. Data collection was performed using the acquisition software OpenSignals (Plux, Lisbon, Portugal) with a sample rate of 1000 Hz. Hair removal, abrasion, and alcohol cleaning were performed prior to electrode placement to improve skin-electrode conductivity. The participant was asked to contract the muscles of interest prior to fixation of the electrodes to help place the electrode over the muscle belly. The electrodes were placed at a 20 mm center-to-center distance parallel to the muscle fibers according to Seniam Recommendations [54]. Trapezius’ electrodes were placed at 50% (bulkier area of the muscle) on the line from the acromion to the spine on vertebra C7 as the reference; deltoideus anterior electrodes were set up at one finger width distal and anterior to the acromion with the orientation of the line between the acromion and the thumb; the electrodes on the biceps brachii were placed at 1/3 from the fossa cubit between this localization and the medial acromion. We only considered the dominant limb muscles. The reference electrode was placed on the olecranon surface.
The raw data were stored and digitally filtered with a 4th order bandpass Butterworth filter with cutoff frequencies of 10 and 500 Hz. Then, the root mean square was applied to the signal. The median frequency was used to quantify the fatigue in terms of frequency-domain with its power spectral density of the raw EMG, using Welch’s modified periodogram. Time-domain analysis comprised the normalization process based on MVC. The smoothed signals were normalized by the mean value of the 100 ms window around the maximum value acquired from MVC.
MVC procedure was performed in standing position for both muscles. Each muscle was submitted to three trials of 3–4 s with a rest interval of 20 s between them. Before starting to collect, participants were instructed about each MVC exercise position and how to produce force. The participants began producing force gradually against researcher’s resistance to avoid EMG bursts. For the biceps brachii MVC, the participants’ elbows were ipsilaterally placed at 90° with arms parallel and leaning against the trunk. The force production was exerted against the researcher’s hands that offered resistance in the region of the participant’s wrists. For trapezius’ MVC, the participant shrugged his shoulders as much as he could against the resistance of the investigator’s hands. During this exercise, the researcher was positioned posteriorly in relation to the participant. Verbal encouragement was used to motivate participants to attain maximum performance during each MVC exercise.

2.5. Accelerometer

A triaxial accelerometer (Plux, Lisbon, Portugal) was placed on the right wrist of each participant. Accelerometer XYZ axes were digitally filtered with a 4th order bandpass Butterworth filter with cutoff frequencies of 0.1 and 10 Hz. Then, each signal was smoothed with an average window of 0.2 s. The y-axis was used to establish the task cycles within each trial (baseline, fatigue 1, and fatigue 2). For this purpose, Hilbert transform was applied, and the normalized transition maximum was used to determine the start and the end of each cycle. Then, visual inspection of each trial was performed to ensure the adequate cycle segmentation. To quantify wrist motion, accelerometer axes were used to calculate the pitch and the roll through the following Euler equations:
Pitch = arctan ( A c c x A c c y 2 + A c c z 2 )
R o l l = a r c t a n ( A c c y A c c z )
where A c c x , A c c y , and A c c z correspond to the xyz axes acceleration.
Having as reference the moment when the participant is grabbing the box, the x-axis, y-axis, and the z-axis correspond to the medio-lateral, antero-posterior, and superior-inferior directions, respectively. Thus, the pitch can be essentially associated with the pronation/supination of the forearm and adduction/abduction of the arm. The roll can be related to flexion/extension of the forearm and arm.

2.6. Statistical Analysis

The statistical analysis comprises mixed effects regression modeling considering inter-participant variability as random effect. For this purpose, the Python package statsmodels was used [55]. The choice of this analysis was due to its greater flexibility in dealing with inter-subject variability over time and between trials in repeated measures designs. Two approaches were followed: (1) one of the dependent variables such as RIP correlation, RIP synchronization, EMG amplitude, EMG median frequency with the following independent variables: trials (Baseline, Fatigue 1, and Fatigue 2) and divisions of each trial (each trial was divided by deciles parts, corresponding about 1 min each); (2) as dependent variables the RIP correlation or RIP synchronization and as independent variables the EMG amplitude and EMG median frequency keeping the trials and the divisions in each trial as previously stated. The assumptions of independence, normality, and homoscedasticity of residuals were verified using the Wald test, Shapiro–Wilk test, and the White’s Lagrange Multiplier test, respectively. Product moment of Pearson was applied to quantify the correlation between significant muscle activity contribution (in the mixed effects model) and RIP correlation and/or synchrony. Normality correction was performed using Yeo-Johnson power transformation [56]. After this transformation, z-score was applied to avoid mean and standard deviation range bias among variables due to the previous power normalization.
Regarding accelerometer results, that is, to verify whether the number of cycles, pitch, and roll change over the trials, a repeated measures ANOVA was used. The sphericity was verified through the Mauchly’s test, and when it was not present the Greenhouse-Geisser correction was performed. Pairwise comparisons were performed with paired t-tests with Bonferroni correction. This analysis was done with the Python package pingouin [57]. The level of significance was set at 5%.

3. Results

This section begins by presenting the results of RIP correlation and phase synchronization and muscle activity separately. Then, the relationship of muscle activity with both RIP correlation and phase synchronization is explored. Finally, we present accelerometry results that account for changes in movement across the three trials. The figures presented for RIP correlation and phase synchronization, EMG amplitude, and median EMG frequency show real values. In contrast, the statistical results presented in the tables stem from the transformed values to ensure the applicability assumptions.

3.1. RIP

Figure 3 and Figure 4 illustrate the mean and standard deviation of rib cage and abdominal wall correlation and phase synchronization for each trial (Baseline, Fatigue 1, and Fatigue 2) across their deciles. Table 1 shows the mixed effect models results taking the rib cage and abdominal correlation and phase synchronization as dependent variables.
The mean value of correlation was 0.65 ± 0.20 for Baseline, 0.57 ± 0.21 for Fatigue 1, and 0.58 ± 0.20 for Fatigue 2. Phase synchrony showed coefficients of 0.65 ± 0.10, 0.61 ± 0.11, and 0.62 ± 0.09 for Baseline, Fatigue 1, and Fatigue 2, respectively.
Fatigue 1 and Fatigue 2 presented significant differences from baseline condition in both RIP (p < 0.001) correlation and phase synchronization (p < 0.001). The division only showed significant differences for RIP synchronization in decile 10 (p = 0.026) for the correlation.

3.2. EMG

Similar to the presentation for RIP results, Figure 5 and Figure 6 show the evolution of biceps and trapezius EMG amplitude using %MVC and median frequency, during each trial across their deciles, respectively. The biceps brachii was activated at 16.7 ± 5.9% MVC during Baseline, increasing to 19.0 ± 6.4% MVC in the Fatigue 1 trial, and then, achieving 18.2 ± 5.8% MVC in Fatigue 2, showing a significant main effect (p < 0.001). The median frequency was about 62.8 ± 8.4 HZ, 63.3 ± 8.3 HZ, and 62.2 ± 8.5 HZ, but the biceps only showed a significant effect in Fatigue 2 trial (p = 0.007) and was similar among trials (p > 0.050).
Trapezius muscle also increased significantly from trial to trial (p < 0.001), showing amplitude values of 12.7 ± 5.3% MVC, X ± Y 13.5% MVC, and 14.2 ± 6.3% MVC. Median frequency also presented a significant effect, but contrary to biceps brachii, it was present in both Fatigue trials (p < 0.010), presenting values of 78.5 ± 11.9 Hz, 80.8 ± 12.9, and 79.7 ± 13.1 Hz across the three trials. The statistical results can be seen in Table 2.

3.3. RIP vs. EMG

The following analysis was performed to verify whether the muscle activity has a significant effect on RIP correlation and phase synchronization despite having revealed significant differences across trials. As in previous models, there were significant main effects of the trials predicting RIP correlation and phase synchrony for both amplitude and frequency analysis p ∈ [<0.001, 0.003]); however, only trapezius activation showed a significant effect p ∈ [<0.001, 0.001]) in the prediction of these two parameters regarding amplitude of activation. On the contrary, only the biceps brachii had a significant association of the median frequency (p = 0.006), but just for phase synchrony prediction, not for correlation (p = 0.115). The statistical results can be seen in Table 3.
Figure 7 illustrates joint plots with the distribution and scatter plot of the muscle activity variables that significantly contributed to RIP correlation and phase synchronization in the previous mixed-effect models. Of these variables, a small and negative correlation was found between trapezius amplitude (average EMG) and both RIP correlation ( r = 0.282 ,   p < 0.001 ) and phase synchrony ( r = 0.295 ,   p < 0.001 ) and positive between and biceps median frequency and phase synchrony ( r = 0.241 ,   p < 0.001 ).

3.4. Accelerometer

Figure 8 illustrates the results for the number of cycles performed by the participants and both pitch and roll of the accelerometer placed on the wrist. No significant differences were found for the number of cycles ( F ( 1.436 ,   30.164 ) = 1.883 ; p = 0.165 ;   η p 2 = 0.082 ) despite the decreasing tendency.
Significant differences were found for the pitch ( F ( 2 ,   42 ) = 11.537 ; p < 0.001 ;   η p 2 = 0.355 ) between the baseline and the two fatigue trials (p = [0.001, 0.013]). The roll showed significant differences ( F ( 2 ,   42 ) = 4.403 ; p = 0.018 ;   η p 2 = 0.173 ) between the two fatigue trials ( p = 0.041 ), and it was tendencialy between the baseline and the Fatigue 2 trial ( p = 0.050 ).

4. Discussion

This study analyzes the evolution of RIP correlation and synchronization as fatigability increases during a simulated protocol of repetitive work. Accordingly, the RIP relationship with EMG (amplitude and median frequency) was also studied to analyze the association between respiration adjustments with local fatigue. To achieve these purposes, participants were asked to perform three trials of 10 min each while performing a repetitive task. Between each pair of trials, a fatigue protocol was performed to ensure that an agonist muscle will act after muscular failure. The participant performed the task after local fatigue which should influence the general body adjustments while, repeatedly, he/she is trying to accomplish the task fighting against tiredness. RIP correlation and phase synchronization represent a response to these adjustments. Our results showed that RIP correlation and phase synchronization change from baseline when compared to fatigue trials, showing tendentially a decrease in thoracic-abdominal synchronization across trials. These results were in agreement with the local fatigability of the agonist muscle (i.e., biceps brachii) and the shoulder fixator (i.e., upper trapezius) by means of EMG amplitude rather than median frequency. However, only the trapezius activation showed a significant association with both RIP correlation and phase synchrony. Although the median frequency of biceps brachii had a significant association with phase synchrony, the results were not consistent with RIP correlation. In addition, there was found inter-subject variability suggesting different motor strategies.

4.1. RIP

The thoracic-abdominal synchronization is dependent on the coordination of how respiratory muscles act as oscillators during respiratory cycles. The efficiency of this action is correlated with their inspiratory muscle strength [58]. Inspiratory muscle strength will be limited if the driving pressure is reduced or when vascular conductance reaches its maximum capacity regarding the requirements, resulting in some kind of alveolar hypoventilation and hypercapnia [59]. In case of fatigue or as the demand is increasing, the central nervous system through the respiratory centers in the brainstem stimulate the inspiratory neurons to increase ventilation in order to promote oxygen supply and remove carbon dioxide, maintaining the acid-basic balance. Thus, as the contraction requirement of non-respiratory muscles increases, so does the contribution of respiratory muscles. Although this RIP technique has limitations quantifying respiratory muscle fatigue [42], a strong relationship between physical demand and respiratory response was reported while performing construction tasks [21]. Thoracic dysfunctional breathing can be provoked by both physiological or psychological stress, in which biomechanical patterns may be respiratory and/or non-respiratory [60]. Thus, our article explores the relationship of thoracic-abdominal synchronization when non-respiratory muscles are susceptible to fatigability during repetitive tasks. Commonly, the variables that are used to quantify the thoracic-abdominal asynchrony are phase difference (φ), percentage of rib cage contribution, and labored breathing index (e.g., [61,62]).
Respiratory muscles act as mechanical effectors and can be divided into three main groups: (1) the inspiratory muscles, (2) the expiratory muscles, and (3) the accessory muscles of respiration [63]. The inspiratory muscle that has the greatest preponderance in inspiration is the diaphragm, whose contraction causes it to descend, increasing both abdominal pressure and thoracic volume. The two sheets of intercostal muscles, the internal and external intercostal muscles, contribute to expanding the rib cage during inspiration and to decreasing rib cage volume during expiration, respectively. The abdominal muscles mainly participate in expiration but also contribute to inspiration, having a relevant role in storing elastic recoil energy in the chest wall, assisting its spread for the next inspiration [63].
Respiratory muscle failure due to fatigue or disease results in a reduction of the capacity in generating force when these muscles’ energy supply is not enough during increased work demand, or when respiratory work requirements are exceeded [64,65]. In this scenario, the fine coordination between inspiratory and expiratory muscles is compromised in the role of avoiding rib cage distortion and lung pressure control, increasing thoraco-abdominal asynchrony [66]. In our study, the phase difference was the chosen metric complemented by the signal correlation. The highest phase asynchrony should occur when one compartment moves in the opposite direction to the other.
Our results showed synchrony alterations when fatigue is present, namely after action with higher local muscle requirements. Both RIP synchronization and correlation decrease when comparing baseline to fatigue trials. Regarding continuous working without stopping, only RIP correlation was sensitive to changes between rib cage and abdominal walls, showing a significant association during the 7th and 10th decile. This disrupted alteration across the trial could represent motor control adjustments due to proprioceptive afferent feedback.
Usually, two major types of fatigue are reported (e.g., [67,68]): central and peripheral. Central fatigue occurs from the synthesis and metabolism of central monoamines impairment, resulting in a reduction in motor unit recruitment, which, in turn, could negatively alter physical and mental efficiency [69]. Peripheral fatigue occurs either due to disabled impulse propagation across the neuromuscular junction and/or over the muscle surface membrane or due to the cross-bridges cycle impairment [65], which biomarkers depend on the type of effort and underlying muscle energy mechanism [70].
In the present study, two EMG methods were used to quantify muscle fatigue or fatigability: amplitude normalized from MVC and median frequency. It is expected that, as fatigue increases, the amplitude of the signal also increases, and the opposite occurs with the median frequency due to the shift toward lower frequencies in the power spectrum. In our study, the EMG amplitude showed higher sensitivity than median frequency regarding repetitive work. Although there is a decrease during each 10-min trial, the median frequency change was not significant among trials. One can point out that the shift of median frequency to lower frequencies depends on alterations of motor units’ potentials actions, showing an increase in their duration and a decrease in the conduction velocity [71]. However, we studied a dynamic rather than isometric that was performed at an intensity lower than 30% of MVC [72], compromising the visible shift of frequencies to lower parts of the frequency spectrum. The task of our study was designed to mimic (as much as possible) repetitive work that was reported to be lower than 20% of MVC [24].
The mean activity of the trapezius in our study was lower than 20% of MVC, but its amplitude showed significant effects on both RIP correlation and phase synchronization rather than biceps brachii. Interestingly, other studies have reported the fatigue of this muscle during repetitive work or long periods of standing [73,74]. This occurs due to the role of the shoulder allowing kinematic and muscular adjustments to promote performance sequentially during prolonged repetitive work [74]. On the other hand, the biceps brachii achieves higher activation in the second trial compared to the third, which suggests biomechanical adaptations due to localized fatigue that are compensated with other muscles or motor strategies, as previously demonstrated during a cyclic lifting task [75]. Indeed, the accelerometer placed on the wrist (Supplemental Material) showed from trial to trial a decrease in the anterior-posterior movement and an increase of medio-lateral. This could mean that as fatigability increases in the agonist muscle (the biceps brachii in this case) the main action of flexion/extension starts to become compromised leading to more abduction/adduction, and then higher activation of the shoulder complex muscles. This type of behavior is according to others that reported a repositioning of joints to maintain task performance during repetitive work [74,76].

4.2. RIP vs. Muscle Activity

To analyze the contribution of local fatigue with alterations in RIP synchronization and phase synchrony, mixed effect regression modeling was performed, adding muscle activity as independent variables. In EMG amplitude, only trapezius muscle showed a significant association with both RIP correlation (p < 0.001) and phase synchrony (p = 0.001). However, only the biceps brachii presented a significant contribution (p = 0.005) to the model of phase synchronization; but not for correlation (p = 0.126). These results highlight the role of trapezius in the shoulder to promote joint repositioning and compensatory action due to fatigue of agonist muscles; as stated in the previous paragraph. Accordingly, the muscles’ local fatigue produces changes in the general functional profile as a consequence of adaptations of motor control strategy [77]. Joint movement trade-offs to redistribute loads across the upper extremity and trunk during repetitive work were previously reported [1], allowing to continue task performance with the presence of fatigue.
The joint repositioning and compensatory alterations due to the local fatigue of the studied agonist muscle and the higher activation of trapezius are supported by the accelerometry results, although limited to the wrist information only. Both pitch and roll tend to increase in amplitude over the trials. In the case of pitch, the baseline differs from both fatigue trials, but not when fatigue is already installed (i.e., between fatigue 1 and fatigue 2 trials). In the case of roll, it is the last fatigue trial (i.e., fatigue 2) that differs from both the baseline and the previous fatigue state. This occurs without significant differences in the number of cycles over the three trials, which in turn reinforces that to maintain performance during accumulated fatigue, the participants performed adaptations of motor control strategy, as stated before [77]. On the one hand, the participants tend to increase pronation/supination of the forearm and/or adduction/abduction of the shoulder during fatigue conditions to compensate for the wear and tear of agonist muscle activity. On the other hand, to maintain performance when accumulated fatigue is present, paradoxically, it increases again the flexion/extension of the shoulder and elbow.
Additionally, due to non-ventilatory functions in tasks where arms are unsupported, some of the inspiratory muscles of the rib cage decrease their ventilatory contribution, shifting the dynamic work to the diaphragm and abdominal muscles of exhalation [78], which leads to RIP asynchrony. Moreover, the task used in this study was performed in the standing position requiring permanent re-establishment adaptation to keep balance while grabbing, placing, and moving a box concentrically and eccentrically.
Another issue that should be taken into consideration is that these motor strategies are not so easy to standardize because there is a high inter-subject variability [1]. We also encountered this issue in the present study, as shown in Figure 7. All the participants changed their behavior through the studied variables, but not in the same fashion.

4.3. Limitations

We highlight two limitations: (1) Statistical analysis was performed using data transformation to guarantee the assumptions of applicability due to inter-subject variability, instead of raw values; (2) a limitation that can compromise the reproducibility of studies is the determination of the respiratory cycle in dynamic actions, since the signal becomes contaminated with adaptations of the thoracic and abdominal grid wall due to the movement of the upper limbs. The expansion of the trunk volume is not merely dependent on the influence of the respiratory cycles, but also on the action of other body segments. For example, the action of flexion/extension and abduction/adduction movements of the arms will be incorporated into the respiratory signal of chest volume expansion/reduction. Unlike at rest, the disparity of different actions in labor context and the adaptation of the necessary processing may lead to divergent interpretations.

4.4. Future Work

The choice of the data collection setup with multiple sensors used in our study was due to the need to understand these indicators in a controlled environment such as a laboratory. Thus, the main recommendation for future work is to collect this type of data in a real environment setting, searching for changes as the working shift lasts. Likewise, comparing different types of workstations and workers’ morphological conditions is recommended. It should be relevant to study RIP correlation and synchronization in both between and within designs during different shifts and chronotypes since both sleep [79] and circadian cycles [80] were associated with work-related fatigue. Studies that comprise inter-subject variability also should be considered to better adjust intervention to individual motor strategies. A great challenge will be to establish the RIP guidelines in which dropping below a certain value is advisable for workers to take a break, namely considering the individual variability.

5. Conclusions

Respiratory inductance plethysmography correlation and synchronization during repetitive work change after muscle failure of tasks’ agonist muscle. Accordingly, both the agonist (biceps brachii in this case) and the shoulder stabilizer muscle (trapezius) increase their amplitude activity from Baseline to Fatigue status. The median frequency of the latter was more discriminative in Fatigue trials than biceps brachii, which in turn showed significant effects only during the second Fatigue task. Remarkably, the biceps brachii amplitude showed significant effects across the 10-min trials, namely from the middle to the end of the trial duration.
Despite the previous results for Bicep brachii amplitude, this muscle did not show a significant effect on RIP correlation and synchronization, unlike the trapezius muscle. This muscle showed a significant contribution to predicting both respiratory parameters. This result could be linked to compensatory activity of trapezius in repositioning the shoulder complex to maintain task performance during repetitive work.
This study supports that RIP can be an interesting quantitative approach to evaluate fatigability during repetitive work because it has been shown that this tool is sensitive to fatigue, is easy to use, and is a non-invasive solution. Thus, it might help decision-makers in ergonomic risk assessment studies and interventions as well as in tasks distribution throughout the day in the production line.

Author Contributions

L.S. and M.D.: conceptualization, methodology, formal analysis, data curation, writing—original draft; D.F. and M.N.: conceptualization, data curation, writing—review and editing, visualization; P.N. and B.A.: conceptualization, writing—review and editing, visualization; M.C.: project administration, writing—review and editing, and funding acquisition; D.C.: writing—review and editing, visualization; E.E. and H.G.: conceptualization, supervision, project administration, writing—review and editing, and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Project OPERATOR (NORTE01-0247-FEDER-045910), co-financed by the ERDF—European Regional Development Fund through the North Portugal Regional Operational Program and Lisbon Regional Operational Program and by the Portuguese Foundation for Science and Technology, under the MIT Portugal Program (2019 Open Call for Flagship projects).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and OPERATOR projected procedures were reviewed by the Ethics Committee of University of Porto, as partner of the project, under the approval number reference 2020/09-6.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Acknowledgments

The authors would like to thank Liliana Cunha and Ana Correia de Barros for their collaboration in the approval of the OPERATOR project procedures by the ethics committee of the University of Porto.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AcronymousDesignation
EMGElectromyography
HRHeart rate
MAWTMaximum acceptable work time
MVCMaximal voluntary contraction
RIPRespiratory inductance plethysmography
A c c x ,   A c c y ,   A c c z xyz axes acceleration

References

  1. Mulla, D.M.; McDonald, A.C.; Keir, P.J. Joint moment trade-offs across the upper extremity and trunk during repetitive work. Appl. Ergon. 2020, 88, 103142. [Google Scholar] [CrossRef] [PubMed]
  2. McDonald, A.C.; Keir, P.J. The Response of the Shoulder Complex to Repetitive Work: Implications for Workplace Design. Crit. Rev. Biomed. Eng. 2015, 43, 21–32. [Google Scholar] [CrossRef] [PubMed]
  3. Santos, J.; Baptista, J.S.; Monteiro, P.R.R.; Miguel, A.S.; Santos, R.; Vaz, M.A. The influence of task design on upper limb muscles fatigue during low-load repetitive work: A systematic review. Int. J. Ind. Ergon. 2016, 52, 78–91. [Google Scholar] [CrossRef] [Green Version]
  4. Tomei, G.; Cinti, M.E.; Cerratti, D.; Fioravanti, M. Attention, repetitive works, fatigue and stress. Ann. Ig. Med. Prev. Comunita 2006, 18, 417–429. [Google Scholar]
  5. Bepko, J.; Mansalis, K. Common occupational disorders: Asthma, COPD, dermatitis, and musculoskeletal disorders. Am. Fam. Physician. 2016, 93, 1000–1006. [Google Scholar]
  6. ILO (International Labour Organization). ILO List of Occupational Diseases. 2010. Available online: https://www.ilo.org/safework/info/publications/WCMS_125137/lang--en/index.htm (accessed on 24 January 2022).
  7. Eurostat 2020. Available online: https://ec.europa.eu/eurostat/statistics-explained/index.php?title=Hours_of_work_-_annual_statistics#In_2020_more_people_worked_fewer_hours_than_usual (accessed on 8 January 2022).
  8. Statista 2022. Available online: https://www.statista.com/statistics/215643/average-weekly-working-hours-of-all-employees-in-the-us-by-month/ (accessed on 15 April 2022).
  9. Messenger, J. Working Time and the Future of Work; ILO Future of Work Research Paper Series; International Labour Organization: Geneva, Switzerland, 2018; pp. 1–35. [Google Scholar]
  10. Leitão, J.; Pereira, D.; Gonçalves, A. Quality of work life and contribution to productivity: Assessing the moderator effects of burnout syndrome. Int. J. Environ. Res. Public Health 2021, 18, 2425. [Google Scholar] [CrossRef]
  11. Björklund, M.; Crenshaw, A.G.; Djupsjöbacka, M.; Johansson, H. Position sense acuity is diminished following repetitive low-intensity work to fatigue in a simulated occupational setting. Eur. J. Appl. Physiol. 2000, 81, 361–367. [Google Scholar] [CrossRef]
  12. Schaub, K.; Caragnano, G.; Britzke, B.; Bruder, R. The European assembly worksheet. Theor. Issues Ergon. Sci. 2013, 14, 616–639. [Google Scholar] [CrossRef]
  13. Hignett, S.; McAtamney, L. Rapid Entire Body Assessment (REBA). Appl. Ergon. 2000, 31, 201–205. [Google Scholar] [CrossRef]
  14. Joshi, M.; Deshpande, V. Investigative study and sensitivity analysis of Rapid Entire Body Assessment (REBA). Int. J. Ind. Ergon. 2020, 79, 103004. [Google Scholar] [CrossRef]
  15. Koshy, K.; Syed, H.; Luckiewicz, A.; Alsoof, D.; Koshy, G.; Harry, L. Interventions to improve ergonomics in the operating theatre: A systematic review of ergonomics training and intra-operative microbreaks. Ann. Med. Surg. 2020, 55, 135–142. [Google Scholar] [CrossRef]
  16. Radwan, A.; Barnes, L.; DeResh, R.; Englund, C.; Gribanoff, S. Effects of active microbreaks on the physical and mental well-being of office workers: A systematic review. Cogent Eng. 2022, 9, 2026206. [Google Scholar] [CrossRef]
  17. Ren, X.; Yu, B.; Lu, Y.; Zhang, B.; Hu, J.; Brombacher, A. LightSit: An unobtrusive health-promoting system for relaxation and fitness microbreaks at work. Sensors 2019, 19, 2162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Romero, D.; Stahre, J.; Taisch, M. The Operator 4.0: Towards socially sustainable factories of the future. Comput. Ind. Eng. 2020, 139, 106128. [Google Scholar] [CrossRef]
  19. El ahrache, K.; Imbeau, D.; Farbos, B. Percentile values for determining maximum endurance times for static muscular work. Int. J. Ind. Ergon. 2006, 36, 99–108. [Google Scholar] [CrossRef]
  20. Aryal, A.; Ghahramani, A.; Becerik-Gerber, B. Monitoring fatigue in construction workers using physiological measurements. Autom. Constr. 2017, 82, 154–165. [Google Scholar] [CrossRef]
  21. Sadat-Mohammadi, M.; Shakerian, S.; Liu, Y.; Asadi, S.; Jebelli, H. Non-invasive physical demand assessment using wearable respiration sensor and random forest classifier. J. Build. Eng. 2021, 44, 103279. [Google Scholar] [CrossRef]
  22. Sun, S.; Zheng, X.; Gong, B.; Garcia Paredes, J.; Ordieres-Meré, J. Healthy operator 4.0: A human cyber–physical system architecture for smart workplaces. Sensors 2020, 20, 2011. [Google Scholar] [CrossRef] [Green Version]
  23. Romero, D.; Stahre, J.; Wuest, T.; Noran, O.; Bernus, P.; Fast-Berglund, Å.; Gorecky, D. Towards an operator 4.0 typology: A human-centric perspective on the fourth industrial revolution technologies. In Proceedings of the International Conference on Computers and Industrial Engineering, Tianjin, China, 29–31 October 2016. [Google Scholar]
  24. Looze, M.D.; Bosch, T.; Dieën, J.V. Manifestations of shoulder fatigue in prolonged activities involving low-force contractions. Ergonomics 2009, 52, 428–437. [Google Scholar] [CrossRef]
  25. Henning, R.A.; Sauter, S.L.; Krieg, E.F. Work rhythm and physiological rhythms in repetitive computer work: Effects of synchronization on well-being. Int. J. Hum.-Comput. Interact. 1992, 4, 233–243. [Google Scholar] [CrossRef]
  26. Schleifer, L.M.; Ley, R.; Spalding, T.W. A hyperventilation theory of job stress and musculoskeletal disorders. Am. J. Ind. Med. 2002, 41, 420–432. [Google Scholar] [CrossRef] [PubMed]
  27. Toledo-Pérez, D.C.; Rodríguez-Reséndiz, J.; Gómez-Loenzo, R.A.; Jauregui-Correa, J.C. A study of computing zero-crossing methods and an improved proposal for EMG signals. Appl. Sci. 2019, 9, 4402. [Google Scholar] [CrossRef] [Green Version]
  28. Dennerlein, J.T.; Ciriello, V.M.; Kerin, K.J.; Johnson, P.W. Fatigue in the forearm resulting from low-level repetitive ulnar deviation. AIHA J. 2003, 64, 799–805. [Google Scholar] [CrossRef]
  29. Dahlqvist, C.; Enquist, H.; Löfqvist, L.; Nordander, C. The effect of two types of maximal voluntary contraction and two electrode positions in field recordings of forearm extensor muscle activity during hotel room cleaning. Int. J. Occup. Saf. 2020, 26, 595–602. [Google Scholar] [CrossRef]
  30. Bosch, T.; De Looze, M.P.; Kingma, I.; Visser, B.; Van Dieen, J.H. Electromyographical manifestations of muscle fatigue during different levels of simulated light manual assembly work. J. Electromyogr. Kinesiol. 2009, 19, e246–e256. [Google Scholar] [CrossRef]
  31. Hagberg, M. The amplitude distribution of surface EMG in static and intermittent static muscular performance. Eur. J. Appl. Physiol. Occup. Physiol. 1979, 40, 265–272. [Google Scholar] [CrossRef]
  32. Carnide, F.; Veloso, A.; Gamboa, H.; Caldeira, S.; Fragoso, I. Interaction of biomechanical and morphological factors on shoulder workload in industrial paint work. Clin. Biomech. 2006, 21, S33–S38. [Google Scholar] [CrossRef]
  33. Village, J.; Trask, C. Ergonomic analysis of postural and muscular loads to diagnostic sonographers. Int. J. Ind. Ergon. 2007, 37, 781–789. [Google Scholar] [CrossRef]
  34. Sundelin, G. Patterns of electromyographic shoulder muscle fatigue during MTM-paced repetitive arm work with and without pauses. Int. Arch. Occup. Environ. Health 1993, 64, 485–493. [Google Scholar] [CrossRef]
  35. Yu, N.; Hong, L.; Guo, J. Analysis of upper-limb muscle fatigue in the process of rotary handling. Int. J. Ind. Ergon. 2021, 83, 103109. [Google Scholar] [CrossRef]
  36. Tuček, D.; Dombeková, B. Local Muscular Load Measurement with the Help of a Datalogger. Acta Polytech. Hung. 2017, 14, 215–234. [Google Scholar]
  37. Lee, S.J.; Moon, M.K.; Choi, W.J.; Jang, T.W. Exercise and cardiovascular load in workers with high occupational physical activity. Arch. Environ. Occup. Health 2020, 75, 339–345. [Google Scholar] [CrossRef] [PubMed]
  38. Preisser, A.M.; Zhou, L.; Garrido, M.V.; Harth, V. Measured by the oxygen uptake in the field, the work of refuse collectors is particularly hard work: Are the limit values for physical endurance workload too low? Int. Arch. Occup. Environ. Health 2016, 89, 211–220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Hjelm, E.W.; Winkel, J.; Nygård, C.H.; Wiktorin, C.; Karlqvist, L. Can cardiovascular load in ergonomic epidemiology be estimated by self-report? Stockholm MUSIC 1 Study Group. J. Occup. Environ. Med. 1995, 37, 1210–1217. [Google Scholar] [CrossRef]
  40. Lunde, L.K.; Koch, M.; Veiersted, K.B.; Moen, G.H.; Wærsted, M.; Knardahl, S. Heavy physical work: Cardiovascular load in male construction workers. Int. J. Environ. Res. Public Health 2016, 13, 356. [Google Scholar] [CrossRef] [Green Version]
  41. Wu, H.C.; Wang, M.J.J. Relationship between maximum acceptable work time and physical workload. Ergonomics 2002, 45, 280–289. [Google Scholar] [CrossRef]
  42. Hammer, J.; Newth, C.J.L. Assessment of thoraco-abdominal asynchrony. Paediatr. Respir. Rev. 2009, 10, 75–80. [Google Scholar] [CrossRef]
  43. Epstein, S.K. An overview of respiratory muscle function. Clin. Chest Med. 1994, 15, 619–639. [Google Scholar] [CrossRef]
  44. Wixted, F.; O’Riordan, C.; O’Sullivan, L. Inhibiting the physiological stress effects of a sustained attention task on shoulder muscle activity. Int. J. Environ. Res. Public Health 2018, 15, 115. [Google Scholar] [CrossRef] [Green Version]
  45. Strang, A.; Ryan, L.; Rahman, T.; Balasubramanian, S.; Hossain, J.; Heinle, R.; Shaffer, T.H. Measures of respiratory inductance plethysmography (RIP) in children with neuromuscular disease. Pediatr. Pulmonol. 2018, 53, 1260–1268. [Google Scholar] [CrossRef]
  46. Mehta, R.K.; Peres, S.C.; Kannan, P.; Rhee, J.; Shortz, A.E.; Mannan, M.S. Comparison of objective and subjective operator fatigue assessment methods in offshore shiftwork. J. Loss Prev. Process Ind. 2017, 48, 376–381. [Google Scholar] [CrossRef]
  47. Atkins, J.H.; Mandel, J.E.; Weinstein, G.S.; Mirza, N. A pilot study of respiratory inductance plethysmography as a safe, noninvasive detector of jet ventilation under general anesthesia. Anesth. Analg. 2010, 111, 1168–1175. [Google Scholar] [CrossRef] [PubMed]
  48. Mannée, D.; Van Helvoort, H.; De Jongh, F. The feasibility of measuring lung hyperinflation with a smart shirt: An in vitro study. IEEE Sens. J. 2020, 20, 15154–15162. [Google Scholar] [CrossRef]
  49. Chen, Y.C.; Hsiao, T.C. Instantaneous phase difference analysis between thoracic and abdominal movement signals based on complementary ensemble empirical mode decomposition. Biomed. Eng. Online 2016, 15, 1–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Deering, R.; Kaiser, J.F. The use of a masking signal to improve empirical mode decomposition. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing, Philadelphia, PA, USA, 23 March 2005. [Google Scholar]
  51. Tsai, F.F.; Fan, S.Z.; Lin, Y.S.; Huang, N.E.; Yeh, J.R. Investigating power density and the degree of nonlinearity in intrinsic components of anesthesia EEG by the Hilbert-Huang transform: An example using ketamine and alfentanil. PLoS ONE 2016, 11, e0168108. [Google Scholar]
  52. Fabus, M.S.; Quinn, A.J.; Warnaby, C.E.; Woolrich, M.W. Automatic decomposition of electrophysiological data into distinct nonsinusoidal oscillatory modes. J. Neurophysiol. 2021, 126, 1670–1684. [Google Scholar] [CrossRef]
  53. Kak, S.C. The discrete Hilbert transform. Proc. IEEE 1970, 58, 585–586. [Google Scholar] [CrossRef]
  54. Hermens, H.J.; Freriks, B.; Disselhorst-Klug, C.; Rau, G. Development of recommendations for SEMG sensors and sensor placement procedures. J. Electromyogr. Kinesiol. 2000, 10, 361–374. [Google Scholar] [CrossRef]
  55. Seabold, S.; Perktold, J. Statsmodels: Econometric and statistical modeling with python. In Proceedings of the 9th Python in Science Conference, Austin, TX, USA, 28 June–3 July 2010. [Google Scholar]
  56. Yeo, I.-K.; Johnson, R.A. A New Family of Power Transformations to Improve Normality or Symmetry. Biometrika 2000, 87, 954–959. [Google Scholar] [CrossRef]
  57. Vallat, R. Pingouin: Statistics in Python. J. Open Source Softw. 2018, 3, 1026. [Google Scholar] [CrossRef]
  58. Weiner, P.; Suo, J.; Fernandez, E.; Cherniack, R.M. Efficiency of the Respiratory Muscles in Healthy Individuals. Am. Rev. Respir. Dis. 1989, 140, 2–396. [Google Scholar] [CrossRef] [PubMed]
  59. Roussos, C.; Koutsoukou, A. Respiratory failure. Eur. Respir. J. 2003, 22, 3s–14s. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Depiazzi, J.; Everard, M.L. Dysfunctional breathing and reaching one’s physiological limit as causes of exercise-induced dyspnoea. Breathe 2016, 12, 120–129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Ratnagiri, M.V.; Zhu, Y.; Rahman, T.; Theroux, M.; Tomatsu, S.; Shaffer, T.H. Automated Assessment of Thoracic-Abdominal Asynchrony in Patients with Morquio Syndrome. Diagnostics 2021, 11, 880. [Google Scholar] [CrossRef] [PubMed]
  62. Ryan, L.; Rahman, T.; Strang, A.; Heinle, R.; Shaffer, T.H. Diagnostic differences in respiratory breathing patterns and work of breathing indices in children with Duchenne muscular dystrophy. PLoS ONE 2020, 15, e0226980. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Mccool, F.D.; Hilbert, J.; Wolfe, L.F.; Benditt, J.O. The respiratory system and neuromuscular diseases. In Murray and Nadel’s Textbook of Respiratory Medicine; Elsevier: Philadelphia, PA, USA, 2022; pp. 1812–1828.e5. [Google Scholar]
  64. Praud, J.P.; Redding, G.J. Chest wall and respiratory muscle disorders. In Kendig’s Disorders of the Respiratory Tract in Children; Elsevier: Amsterdam, The Netherlands, 2019; pp. 1044–1061.e2. [Google Scholar] [CrossRef]
  65. Roussos, C.; Zakynthinos, S. Fatigue of the respiratory muscles. Intensive Care Med. 1996, 22, 134–155. [Google Scholar] [CrossRef]
  66. Fregonezi, G.; Sarmento, A.; Pinto, J.; LoMauro, A.; Resqueti, V.; Aliverti, A. Thoracoabdominal asynchrony contributes to exercise limitation in mild asthmatic subjects. Front. Physiol. 2018, 9, 719. [Google Scholar] [CrossRef]
  67. Carroll, T.J.; Taylor, J.L.; Gandevia, S.C. Recovery of central and peripheral neuromuscular fatigue after exercise. J. Appl. Physiol. 2017, 122, 1068–1076. [Google Scholar] [CrossRef] [Green Version]
  68. Davis, J.M. Central and peripheral factors in fatigue. J. Sports Sci. 1995, 13, S49–S53. [Google Scholar] [CrossRef]
  69. Meeusen, R.; Watson, P.; Hasegawa, H.; Roelands, B.; Piacentini, M.F. Central fatigue. Sports Med. 2006, 36, 881–909. [Google Scholar] [CrossRef]
  70. Finsterer, J. Biomarkers of peripheral muscle fatigue during exercise. BMC Musculoskelet. Disord. 2021, 13, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Lowery, M.M.; Vaughan, C.L.; Nolan, P.J.; O’Malley, M.J. Spectral compression of the electromyographic signal due to decreasing muscle fiber conduction velocity. IEEE Trans. Neural Syst. Rehabil. Eng. 2000, 8, 353–361. [Google Scholar] [CrossRef] [PubMed]
  72. De Luca, C.J. The use of surface electromyography in biomechanics. J. Appl. Biomech. 1997, 13, 135–163. [Google Scholar] [CrossRef] [Green Version]
  73. Christensen, H. Muscle activity and fatigue in the shoulder muscles during repetitive work. Eur. J. Appl. Physiol. 1986, 54, 596–601. [Google Scholar] [CrossRef]
  74. McDonald, A.C.; Calvin, T.F.; Keir, P.J. Adaptations to isolated shoulder fatigue during simulated repetitive work. Part II: Recovery. J. Electromyogr. Kinesiol. 2016, 29, 42–49. [Google Scholar] [CrossRef]
  75. Bonato, P.; Ebenbichler, G.R.; Roy, S.H.; Lehr, S.; Posch, M.; Kollmitzer, J.; Della Croce, U. Muscle fatigue and fatigue-related biomechanical changes during a cyclic lifting task. Spine 2003, 28, 1810–1820. [Google Scholar] [CrossRef]
  76. Calvin, T.F.; McDonald, A.C.; Keir, P.J. Adaptations to isolated shoulder fatigue during simulated repetitive work. Part I: Fatigue. J. Electromyogr. Kinesiol. 2016, 29, 34–41. [Google Scholar] [CrossRef]
  77. Goethel, M.F.; Gonçalves, M.; Brietzke, C.; Cardozo, A.C.; Vilas-Boas, J.P.; Ervilha, U.F. A global view on how local muscular fatigue affects human performance. Proc. Natl. Acad. Sci. USA 2020, 117, 19866–19872. [Google Scholar] [CrossRef]
  78. Celli, B.; Criner, G.; Rassulo, J. Ventilatory muscle recruitment during unsupported arm exercise in normal subjects. J. Appl. Physiol. 1988, 64, 1936–1941. [Google Scholar] [CrossRef]
  79. McDonald, J.; Potyk, D.; Fischer, D.; Parmenter, B.; Lillis, T.; Tompkins, L.; Belenky, G. Napping on the night shift: A study of sleep, performance, and learning in physicians-in-training. J. Grad. Med. Educ. 2013, 5, 634–638. [Google Scholar] [CrossRef] [Green Version]
  80. Minelli, A.; Di Palma, M.; Rocchi, M.B.L.; Ponzio, E.; Barbadoro, P.; Bracci, M.; Prospero, E. Cortisol, chronotype, and coping styles as determinants of tolerance of nursing staff to rotating shift work. Chronobiol. Int. 2021, 38, 666–680. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Work cycle sequence: (a) Neutral position at the beginning of the cycle; (b) grabbing the box from the lower shelf; (c) placing the box on the upper shelf; (d) neutral position after returning from the upper shelf placement; (e) grabbing the box from the upper shelf; (f) placing the box on the lower shelf. Fatigue sequence: three trials of biceps curl until muscular failure in each one.
Figure 1. Work cycle sequence: (a) Neutral position at the beginning of the cycle; (b) grabbing the box from the lower shelf; (c) placing the box on the upper shelf; (d) neutral position after returning from the upper shelf placement; (e) grabbing the box from the upper shelf; (f) placing the box on the lower shelf. Fatigue sequence: three trials of biceps curl until muscular failure in each one.
Sensors 22 04247 g001
Figure 2. Example of a participant during rest (left column), baseline (middle column), and after fatigue protocol (right column): raw respiration signals (first row), respiration signals after EMD and normalization (second row), and signal phase (fourth row). In each of these figure rows, the chest movement (light blue) and abdominal movement (dark blue) walls are represented. Correlation and phase synchrony are shown in the third and fifth rows, respectively (both in black).
Figure 2. Example of a participant during rest (left column), baseline (middle column), and after fatigue protocol (right column): raw respiration signals (first row), respiration signals after EMD and normalization (second row), and signal phase (fourth row). In each of these figure rows, the chest movement (light blue) and abdominal movement (dark blue) walls are represented. Correlation and phase synchrony are shown in the third and fifth rows, respectively (both in black).
Sensors 22 04247 g002
Figure 3. Trial mean (red dashed line), means of each decile (blue line), and standard deviation (light blue shadow) of rib cage and abdominal walls correlation for each trial (Baseline, Fatigue 1, and Fatigue 2) across their deciles (trial %). The values of correlation range from “0” to “1”, where 0 corresponds to no correlation and “1” to total correlation.
Figure 3. Trial mean (red dashed line), means of each decile (blue line), and standard deviation (light blue shadow) of rib cage and abdominal walls correlation for each trial (Baseline, Fatigue 1, and Fatigue 2) across their deciles (trial %). The values of correlation range from “0” to “1”, where 0 corresponds to no correlation and “1” to total correlation.
Sensors 22 04247 g003
Figure 4. Trial mean (red dashed line), means of each decile (blue line), and standard deviation (light blue shadow) of rib cage and abdominal walls phase synchronization for each trial (Baseline, Fatigue 1, and Fatigue 2) across their deciles (trial %). The values of correlation range from “0” to “1”, where 0 corresponds to no synchrony at all (antiphase) and “1” to total synchronization (in phase).
Figure 4. Trial mean (red dashed line), means of each decile (blue line), and standard deviation (light blue shadow) of rib cage and abdominal walls phase synchronization for each trial (Baseline, Fatigue 1, and Fatigue 2) across their deciles (trial %). The values of correlation range from “0” to “1”, where 0 corresponds to no synchrony at all (antiphase) and “1” to total synchronization (in phase).
Sensors 22 04247 g004
Figure 5. Trial mean (red dashed line), means of each decile (blue line), and standard deviation (light blue shadow) of both amplitude (upper) and median frequency (bottom) among the three trials (baseline, Fatigue 1, and Fatigue 2) for biceps brachii, considering their deciles (trial %).
Figure 5. Trial mean (red dashed line), means of each decile (blue line), and standard deviation (light blue shadow) of both amplitude (upper) and median frequency (bottom) among the three trials (baseline, Fatigue 1, and Fatigue 2) for biceps brachii, considering their deciles (trial %).
Sensors 22 04247 g005
Figure 6. Trial mean (red dashed line), Means of each decile (blue line), and standard deviation (light blue shadow) of both amplitude (upper) and median frequency (bottom) among the three trials (baseline, fatigue 1, and fatigue 2) for trapezius, considering their deciles (trial %).
Figure 6. Trial mean (red dashed line), Means of each decile (blue line), and standard deviation (light blue shadow) of both amplitude (upper) and median frequency (bottom) among the three trials (baseline, fatigue 1, and fatigue 2) for trapezius, considering their deciles (trial %).
Sensors 22 04247 g006aSensors 22 04247 g006b
Figure 7. Relationship between RIP correlation/phase synchrony and muscle activity that have shown significant effects in the mixed-effects model. The arbitrary units are due to the Yeo-Johnson power transformation.
Figure 7. Relationship between RIP correlation/phase synchrony and muscle activity that have shown significant effects in the mixed-effects model. The arbitrary units are due to the Yeo-Johnson power transformation.
Sensors 22 04247 g007
Figure 8. Number of cycles (a), pitch (b), and roll (c) results for the three trials (baseline, Fatigue 1, and Fatigue 2). * means significant differences.
Figure 8. Number of cycles (a), pitch (b), and roll (c) results for the three trials (baseline, Fatigue 1, and Fatigue 2). * means significant differences.
Sensors 22 04247 g008
Table 1. Mixed effect regression with the rib cage and abdominal walls correlation and synchronization as dependent variable and trials and divisions as independent variables.
Table 1. Mixed effect regression with the rib cage and abdominal walls correlation and synchronization as dependent variable and trials and divisions as independent variables.
RIP CorrelationRIP Synchrony
β Zp β Zp
Intercept0.3731.9280.0540.2751.4140.157
TrialFatigue 1−0.364−6.153<0.001−0.324−5.490<0.001
Fatigue 2−0.348−5.881<0.001−0.290−4.904<0.001
Division2−0.095−0.8780.380−0.082−0.7580.449
3−0.169−1.5630.118−0.041−0.3800.704
4−0.130−1.2070.228−0.068−0.6350.525
5−0.138−1.2780.2010.0180.1630.870
6−0.095−0.8820.378−0.128−1.1830.237
7−0.207−1.9100.056−0.041−0.3760.707
8−0.120−1.1120.266−0.049−0.4530.650
9−0.153−1.4180.156−0.133−1.2330.217
10−0.241−2.230.026−0.179−1.6580.097
SubjectVariation0.637 0.647
Legend: significant values are in bold.
Table 2. Mixed effect regression with biceps brachii and trapezius EMG amplitude and median frequency as dependent variables.
Table 2. Mixed effect regression with biceps brachii and trapezius EMG amplitude and median frequency as dependent variables.
AmplitudeFrequency
Biceps BrachiiTrapeziusBiceps BrachiiTrapezius
β Zp β Zp β Zp β Zp
Intercept−0.096−0.4670.641−0.066−0.3040.7610.0760.3520.725−0.061−0.2780.781
TrialFatigue 10.43310.207<0.0010.1415.727<0.0010.0501.4230.1550.1696.966<0.001
Fatigue 20.2946.93<0.0010.2459.985<0.001−0.094−2.6970.0070.0712.9130.004
Division2−0.058−0.7530.451−0.039−0.8790.379−0.046−0.720.4720.0170.3780.705
3−0.060−0.7760.438−0.030−0.6710.502−0.091−1.4270.154−0.035−0.7850.432
4−0.089−1.150.250−0.036−0.7980.425−0.061−0.9650.335−0.052−1.1820.237
5−0.189−2.4340.015−0.078−1.7450.081−0.068−1.0750.282−0.03−0.6690.503
6−0.195−2.5150.012−0.086−1.9120.056−0.062−0.9680.333−0.042−0.9470.343
7−0.200−2.5760.010−0.078−1.7350.083−0.058−0.9120.362−0.032−0.710.478
8−0.195−2.5170.012−0.094−2.1020.036−0.059−0.9220.356−0.033−0.7470.455
9−0.226−2.9130.004−0.073−1.6240.104−0.092−1.440.1500.0200.4480.654
10−0.251−3.2470.001−0.109−2.4280.015−0.071−1.1180.263−0.005−0.1080.914
SubjectVariation0.814 0.973 0.914 0.980
Legend: significant values are in bold.
Table 3. Mixed effect regression with the rib cage and abdominal walls correlation and synchronization as dependent variables adding muscle activity to the models.
Table 3. Mixed effect regression with the rib cage and abdominal walls correlation and synchronization as dependent variables adding muscle activity to the models.
AmplitudeFrequency
RIP CorrelationRIP SynchronizationRIP CorrelationRIP Synchronization
β Zp β Zp β Zp β Zp
Intercept0.3421.8180.0690.2451.3040.1920.3551.8250.0680.2491.2880.198
TrialFatigue 1−0.283−4.503<0.001−0.238−3.79<0.001−0.345−5.649<0.001−0.304−5.004<0.001
Fatigue 2−0.241−3.863<0.001−0.185−2.9740.003−0.327−5.438<0.001−0.259−4.328<0.001
Division2−0.113−1.0610.289−0.100−0.9370.349−0.087−0.8080.419−0.07−0.6510.515
3−0.184−1.7280.084−0.056−0.5280.598−0.164−1.5160.130−0.03−0.2750.783
4−0.15−1.4040.16−0.088−0.8310.406−0.131−1.2150.225−0.066−0.6140.539
5−0.18−1.6790.093−0.025−0.2350.814−0.135−1.2480.2120.0260.2380.812
6−0.14−1.3060.192−0.173−1.6220.105−0.095−0.8750.381−0.123−1.1460.252
7−0.249−2.3230.020−0.084−0.7890.43−0.205−1.8940.058−0.035−0.3260.745
8−0.168−1.5670.117−0.097−0.9090.363−0.119−1.0980.272−0.043−0.4040.686
9−0.196−1.8260.068−0.178−1.6630.096−0.14−1.2940.196−0.112−1.0390.299
10−0.298−2.7740.006−0.237−2.2130.027−0.234−2.1630.031−0.166−1.5440.122
MuscleBiceps−0.076−1.2830.199−0.101−1.7070.0880.1121.5760.1150.1942.7390.006
Trapezius−0.346−3.709<0.001−0.304−3.2780.001−0.147−1.5300.126−0.178−1.8560.063
SubjectVariation0.601 0.599 0.647 0.640
Legend: significant values are in bold.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Silva, L.; Dias, M.; Folgado, D.; Nunes, M.; Namburi, P.; Anthony, B.; Carvalho, D.; Carvalho, M.; Edelman, E.; Gamboa, H. Respiratory Inductance Plethysmography to Assess Fatigability during Repetitive Work. Sensors 2022, 22, 4247. https://doi.org/10.3390/s22114247

AMA Style

Silva L, Dias M, Folgado D, Nunes M, Namburi P, Anthony B, Carvalho D, Carvalho M, Edelman E, Gamboa H. Respiratory Inductance Plethysmography to Assess Fatigability during Repetitive Work. Sensors. 2022; 22(11):4247. https://doi.org/10.3390/s22114247

Chicago/Turabian Style

Silva, Luís, Mariana Dias, Duarte Folgado, Maria Nunes, Praneeth Namburi, Brian Anthony, Diogo Carvalho, Miguel Carvalho, Elazer Edelman, and Hugo Gamboa. 2022. "Respiratory Inductance Plethysmography to Assess Fatigability during Repetitive Work" Sensors 22, no. 11: 4247. https://doi.org/10.3390/s22114247

APA Style

Silva, L., Dias, M., Folgado, D., Nunes, M., Namburi, P., Anthony, B., Carvalho, D., Carvalho, M., Edelman, E., & Gamboa, H. (2022). Respiratory Inductance Plethysmography to Assess Fatigability during Repetitive Work. Sensors, 22(11), 4247. https://doi.org/10.3390/s22114247

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop