Next Article in Journal
Metal Oxide Nanoparticles and Nanotubes: Ultrasmall Nanostructures to Engineer Antibacterial and Improved Dental Adhesives and Composites
Previous Article in Journal
Translational Application of 3D Bioprinting for Cartilage Tissue Engineering
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Viscoelastic Characterization of Parasagittal Bridging Veins and Implications for Traumatic Brain Injury: A Pilot Study

by
Silvia García-Vilana
1,†,
David Sánchez-Molina
1,*,†,
Jordi Llumà
1,
Ignasi Galtés
2,3,
Juan Velázquez-Ameijide
1,
M. Carmen Rebollo-Soria
4 and
Carlos Arregui-Dalmases
5
1
Escola d’Enginyeria de Barcelona Est, Universitat Politècnica de Catalunya, Av. Eduard Maristany, 16, 08019 Barcelona, Spain
2
Institut de Medicina Legal i Ciències Forenses de Catalunya, G.V. Corts Catalanes, 111, 08014 Barcelona, Spain
3
Departament de Psiquiatria i de Medicina Legal, Universitat Autònoma de Barcelona, Cerdanyola del Vallès, 08193 Barcelona, Spain
4
Instituto de Medicina Legal de Aragón, Irene Izárbez, 22005 Huesca, Spain
5
Centro Zaragoza, Crtra. 232, km.273, 50690 Pedrola, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Bioengineering 2021, 8(10), 145; https://doi.org/10.3390/bioengineering8100145
Submission received: 19 August 2021 / Revised: 5 October 2021 / Accepted: 12 October 2021 / Published: 18 October 2021

Abstract

:
Many previous studies on the mechanical properties of Parasagittal Bridging Veins (PSBVs) found that strain rate had a significant effect on some mechanical properties, but did not extensively study the viscoelastic effects, which are difficult to detect with uniaxial simple tensile tests. In this study, relaxation tests and tests under cyclic loading were performed, and it was found that PSBVs do indeed exhibit clear viscoelastic effects. In addition, a complete viscoelastic model for the PSBVs is proposed and data from relaxation, cyclic load and load-unload tests for triangular loads are used to find reference values that characterize the viscoelastic behavior of the PSBVs. Although such models have been proposed for other types of blood vessels, this is the first study that clearly demonstrates the existence of viscoelastic effects from an experimental point of view and also proposes a specific model to explain the data obtained. Finally, this study provides reference values for the usual viscoelastic properties, which would allow more accurate numerical simulation of PSBVs by means of computational models.

Graphical Abstract

1. Introduction

The importance of the TBI at the global level cannot be overlooked. Every year, about 70 million people worldwide are afflicted with some form of TBI, with the global incidence reaching about 940 cases per 100,000 people [1]. Among the factors leading to different types of TBI, a particularly pressing one is SDH produced by the rupture of some types of PSBVs [2,3]. In particular, the superior sagittal sinus-bridging veins, which link the central sagittal sinus to the upper part of the brain mass, have been identified as a highly critical region prone to producing SDHs of some severity [4,5,6,7].
For this reason, detailed and accurate studies of the mechanical properties of the collagenous tissue constituting the PSBVs are important, and should be updated using better and more accurate measurement techniques. For example, one influential paper detected a negative correlation of ultimate stress and strain rate [4]; however, a later paper showed a positive correlation [7], although neither author found the correlation to be statistically significant. The issue was recently closed when the use of more accurate digital measurement systems found a positive and significant correlation [3]. Although strain-rate dependent behavior for some mechanical properties have been studied by several authors [4,6,7,8], there are still few studies in the literature properly investigating the viscoelastic behavior of PSBVs, including the measurement of relaxation times or preconditioning under cyclic loading [3,9]. An important advance is found in [10] which examines preconditioning by cyclic loading and fatigue. However, although some work has been done for other types of blood vessels [11], no specific and complete viscoelastic model has yet been proposed for the PSBVs. The aim of this study is to contribute to filling the gap in viscoelastic characterization, using a more complex series of tests including load-unload cycles, stress relaxation cycles and fast repeated loading cycles, following the line of other works on the characterization of viscoelastic behavior of biological tissues [12,13,14].
This work shows that, the viscoelastic behavior of the PSBVs is clearly observable, measurable and quantifiable when observing the rehology of the tissue in a series of more complex tests, such as cyclic-load and relaxation tests, instead of conventional simple uniaxial tests, where the contribution of the viscoelastic effect to the mechanical behavior is difficult to quantify [3]. Such viscoelastic behavior can be used to improve both the understanding of injury mechanisms [15,16], and the computational models used to assess the probability of severe injuries.

2. Data and Methods

2.1. Material and Specimen Preparation

For this study, different sections of the meningeal-cortex space were obtained from four autopsies of post-mortem human subjects (PMHS) conducted at the Forensic Pathology Service of the Legal Medicine and Forensic Science Institute of Catalonia (FPS/IMLCFC). The study was approved by the Research Committee of the IMLCFC. None of the subjects had been previously diagnosed with any blood vessel pathology.
Once received, the meningeal sections were kept refrigerated for, at most, 96 h in airtight containers that maintained the natural degree of hydration. The sections of the meningeal-cortex space included the meninges, the subarachnoid space, and the upper part of the cerebral cortex. From those four sections, six PSBVs were carefully dissected (more than one PSBV was dissected from two sections). After dissection, the PSBVs were kept in contact with towel tissue soaked in phosphate-buffered saline solution. Twelve to twenty-four hours before testing, the PSBVs were placed in a refrigerator at 2 C and before testing they were allowed to come to room temperature for one hour.

2.2. Cyclic and Relaxation Tests

The tests were performed with an Allround-Table-Top Universal Test Machine (UTM) Zwick-Roell® and load was measured with a 20N load cell HBM®. Special fixtures were used for the PSBVs tests and the displacement was carefully measured with a digital control unit attached to the UTM (due to the small dimensions of the cross-section of the PSBVs, no significant local effects appear in the displacement measures). The specimens were kept in physiological saline solution at 5 C to minimize tissue degradation during the testing.
During the tests, a digitally displacement-controlled maximum strain ε max was applied and the force experienced by each PSBV specimen was measured very accurately. In order not to cause permanent damage to the specimens, ε max was selected as 20% of the average ultimate strain, determined in previous tensile studies of PSBVs [3]. The sequence of stages of the loading process was similar to that used in other soft tissue studies such as [12]. The different stages are as follows: (1) load-hold at ε max , (2) rapid load-unload cycles at ε max , (3) three stages of load-hold at ε max , 2/3 · ε max and 1/3 · ε max , (4) three load-unload stages at ε max at stretch rates of 0.01, 0.1 and 1 s 1 and (5) load to failure. Figure 1b shows a schematic of the above stages. The tests were performed at a speed of 50 mm/s (except for stage 4) and the holding time and the time between load stages were 50 s, allowing the specimen to recover completely.
In this study, each PSBV is modeled as a homogeneous, transversely isotropic, hollow cylinder with constant cross-section [3,4,5,7]. The coordinate axes have their origin upon the static bottom fixture; the X axis is aligned with the (longitudinal) force, while the Y and Z axes are contained in the PSBV cross-section. Thus, the stretch for each time instant is easily determined as λ t = t / 0 = 1 + δ t / 0 where 0 is the initial length, t the instantaneous length and the displacement δ t = t 0 has been digitally measured. Let X t = ( X t , Y t , Z t ) be the coordinates in the initial configuration and x t = ( x t , y t , z t ) in the deformed configuration, the deformation of the cylinder is expressed as
x t = X λ t , y t = Y 1 ν ¯ ( λ t 2 1 ) 1 / 2 , z t = Z 1 ν ¯ ( λ t 2 1 ) 1 / 2
where ν ¯ ( λ t ) represents the Poisson effect. Thus, being F t = x t / X the deformation gradient tensor, the (Green-Lagrangian) strain tensor ε = ( F t T F t 1 ) / 2 is obtained as
ε ( t ) = 1 2 λ t 2 1 0 0 0 ν ¯ ( λ t ) ( λ t 2 1 ) 0 0 0 ν ¯ ( λ t ) ( λ t 2 1 )
On the other hand, the force F t has been measured for each time instant during the whole test. From the force, the (second Piola-Kirchhoff) stress tensor σ is obtained:
σ ( t ) = F t λ t A 0 0 0 0 0 0 0 0 0
being A 0 the initial (undeformed) cross-section of the PSBV.

2.3. Quasi-Linear Viscoelastic Model

Many previous studies had used uniaxial tensile tests on PSBVs, showing that an elastic model reasonably explains the results obtained in this type of test. For this reason, the collagenous tissue of PSBVs has been treated simply as an elastic material. However, the cycle tests and the relaxation tests performed show that there is a clear viscoelastic effect on the mechanical behavior of PSBVs and that, therefore, viscoelastic models would work better for more general situations than simple traction.
In this study, a viscoelastic model of type QLVE [13,17] will be used to model the data. In the QLVE model, the mechanical response is divided into a strain-rate dependent part and a strain-rate independent part: σ x = σ x ( e ) + σ x ( v ) , using a separable relaxation function R ( ε x , t ) = G ( t ) · σ x ( e ) / ε x [13,18], where G ( t ) is the relaxation function, given by the well-known Prony series [19,20,21,22]:
G ( t ) = 1 + k = 1 N g k e t / τ k
with g k being the “weights” of each k-term and τ k being the relaxation times, both to be obtained from experimental trials. The value of N depends on the needs of the setting, in Section 3.1 it is shown that N = 3 is sufficient to represent the data in this study. With these assumptions, the axial strain (using the second Piola-Kirchhoff stress tensor) is expressed as
σ x ( t ) = Ψ ε x = 0 t G ( t τ ) σ x ( e ) ε x ε ˙ x ( τ ) d τ

3. Results

For this study, cyclic and uniaxial relaxation tests were performed on six PSBV specimens. Figure 2 clarifies the five stages of the testing process for each specimen. The same figure shows the variation of the imposed displacements over time, together with the axial force measured by the load cell.
Stages (1) and (3) are stages with constant strain: the first one up to the established strain level ε max , and the second one divided in three stages up to ε max , 2 / 3 · ε max and 1 / 3 · ε max . In both stages, while keeping constant the applied displacement and, therefore, the strain, a clear relaxation of the reaction force of the PSBV specimen is observed, with a drop of exponential type, which is a clear indication of the presence of viscoelastic effects in the mechanical behavior of the specimen.
Likewise, stage (2) in which successive loading and unloading cycles have been applied up to ε max at high speed, the same behavior is perceived, showing a progressive drop in the force as the number of accumulated cycles increases.
Finally, stage (4) is divided into three load-unload stages at speeds of 0.01 , 0.1 and 1 s 1 respectively. It can be seen that, the loads being triangular and applied at a constant velocity, the force response is clearly a concave curve, which also indicates the presence of viscoelasticity. Thus, the results obtained from the cycling and relaxation tests indicate that the bridging vein is a viscoelastic material; therefore, a constitutive model that contemplates such viscoelasticity is required to describe its mechanical behavior.

3.1. Relaxation Tests

The parameters of the viscoelastic model of Equation (5) were fitted from the data of stages (1) and (3) of Figure 2. For this purpose, the set of force-time values from the instant at which the maximum displacement is reached until the end of the stage was considered.
The axial stress ( σ x ) is expressed in terms of the force as σ x ( t ) = F t / λ t A 0 , see Equation (3). Given the low level of deformation, ( ε < 5 %), the material can be considered linear elastic [23] without making large errors and employing the following relationship:
F t E A 0 = λ t ε x ( t ) + k = 1 N g k 0 t e ( t τ ) / τ k ε ˙ x ( τ ) d τ
Note that, during the initial instants of loading, ε ˙ x 0 and that makes the integral not identically zero (see that once the maximum displacement ε ˙ x = 0 is reached, the integral term is null). However, at the end of the stage, the specimen is completely relaxed, so the viscoelastic contribution (VC) ends up being zero. Moreover, from practically the beginning of the trial, both λ t and ε x are constants: λ t = λ and ε x ( t ) = ε max . All this implies that the quotient F t / E A 0 also tends to the constant given by λ ε max . Thus, measuring the force F at the end of the stage provides a way to determine the value of E A 0 from the ratio E A 0 = F / ( λ ε max ) . For that reason, calculations can be made without a direct measurement of the elastic modulus or cross-section of the specimen by simply using the normalization f t = F t / E A 0 , where the latter magnitude f t will be called reduced force.
During the relaxation itself, λ t and, therefore, also ε x are constant. However, before reaching their constant value, there is a very fast sudden stretching, during which the elongation is given by λ t = 1 + δ t / 0 , with δ t being the displacement and the stretching rate being λ ˙ t = δ ˙ t / 0 . Experimentally, it was observed that the fast stretching velocity process can be adequately described by a function of “gamma distribution” type [24], with the form
δ ˙ t = Δ max t a t t a n e t / t a Γ ( n + 1 )
where t a is a parameter related to the time required to reach the maximum displacement and, in practice, is obtained by fitting the measured displacement velocity data. Similarly, several preliminary tests showed that taking n = 3 is sufficient to adequately represent the velocity at the beginning of each relaxation stage. Therefore, the displacement δ t can be obtained by integrating the above expression (see Figure 3).
From the explicit form (7) of the velocity δ ˙ t , we obtain the strain rate as ε ˙ x ( t ) = δ t which, introduced in Equation (6), leads to the formula:
f t = F t E A 0 = λ t ε max + k = 1 N g k e α 1 t / θ k ( 1 θ k ) 4 ( 1 2 θ k ) 7 m = 1 11 θ k m ( μ m + ν m e α 2 / θ k )
In this formula, α 1 and α 2 are real numbers that depend on the initial conditions of the specimen (length) and the test velocity curve (time of load application, etc.). In the same way, the explicit analytical integration yields the real numbers μ m and ν m (the same calculation shows that ν m > 0 , μ m < 0 for even m and ν m < 0 , μ m > 0 if m is odd, and furthermore one has, μ 1 , μ 2 , μ 3 , μ 4 = 0 ). The only adjustable parameters in the above expression are θ k = τ k / t a and g k which are dimensionless. The relaxation times are obtained as τ k = θ k t a . The number of adjustable parameters is 2 N , where N is the number of terms in the Prony series, see Equation (4).
Figure 4 shows the measured force data for stage 1 of one of the specimens, and the corresponding fits with the viscoelastic model for N = 1 , N = 2 and N = 3 in the Prony series (4). The fitted parameters have been set out in Table 1. It can be seen that as N increases, the fit improves markedly as expected, since the number of adjustable g k and g k parameters increases:
  • Notice that, although for N = 1 the coefficient r 2 = 0.889 is quantitatively adequate, qualitatively it is observed that the fit does not adequately represent the vertical asymptote of force drop, nor the curve in general.
  • With N = 2 there is a marked improvement in the fit ( r 2 = 0.990 ).
  • Finally, for N = 3 ( r 2 = 0.995 ), both the asymptote and the curve are adequately represented by the viscoelastic model.
Thus, in this study the parameters g k and τ k have been determined for N = 3 ( g 1 , g 2 , g 3 ; τ 1 , τ 2 , τ 3 ) . The viscoelastic model fits have been repeated for stages (1) and (3), the latter divided into subplots (3.1), (3.2) and (3.3), thus obtaining the parameters. The average for each specimen is presented in Table 2.
It can be observed in the table that the higher the value of g k (from g 1 to g 3 ), the lower the value of τ k , thus increasing the degree of contribution of the model term, and the shorter the relaxation time. The relaxation times decrease by an order of magnitude progressively between τ 1 and τ 3 , obtaining similar values among the different specimens. The average characteristic times for the sample are τ 1 = 17.22 ± 5.17 s, τ 2 = 1.49 ± 0.55 s and τ 3 = 0.13 ± 0.05 s. The average values of the contribution coefficients for the sample are g 1 = 0.17 ± 0.06 , g 2 = 0.17 ± 0.08 and g 3 = 0.87 ± 0.62 .
In addition, to show the importance of the viscoelastic effect in PSBVs, the VC for the reduced force has been calculated, i.e., what percentage of the response corresponds to the elastic effect and what amount to the viscoelastic one. The VC for all specimens, calculated from the parameters g k and τ k in Table 2, is shown in Figure 5. It can be seen that at the initial instant, the VC is between 26 and 35%, decreasing exponentially until it reaches approximately 0% at the end of the relaxation stretch. In fact, after 15 s of relaxation, the specimen has relaxed to a great extent and the VC is less than 10% in all samples.
Likewise, the average VC and viscoelasticity and the confidence interval have been determined from the parameters of the studied specimens (Figure 5). Thus, it can be seen that initially the viscoelastic part corresponds to 29.8 ± 3.35%, decreasing to 11.4 ± 2.40% at 10 s and 4.8 ± 1.06% at 25 s, while the elastic contribution starts from 69.8 ± 4.03%, rising to 88.5 ± 2.64% at 10 s to 94.9 ± 1.10% at 25 s.

3.2. Fast Loading/Unloading Cycle Tests

The viscoelastic effect can also be analyzed in test stage (2), in which fast loading and unloading cycles have been performed. These cycles are fast compared to some of the relaxation times τ k of the material. This leads to the assumption that in the time in which a loading and unloading cycle elapses, the BV may not have enough time to relax the internal stresses and this is reflected in the maximum forces reached in each cycle.
Figure 6 shows the force peaks in successive cycles. It can be seen how in each new cycle the maximum force reached is lower, showing a progressive relaxation of the tissue of the PSBV as the number of cycles increases. For each sample, the USF (unconditioned scale factor), defined in [12], has been calculated as the following ratio:
USF ( n , T ) = F 1 F n > 1
where n is the cumulative cycle number and T is the period of the cycles ( T = 0.55 ). As can be seen in the figure above, the magnitude of the USF grows progressively with the number of cycles. In some samples, it has reached 20 cycles, while in others it has been limited to 10 cycles. From the results, it has been calculated that the USF ( n = 10 ) = 1.14 ± 0.085 and USF ( n = 20 ) = 1.23 ± 0.095 . In Appendix A, it is justified that the USF should be able to be approximated by a function of the type
USF ( n , T ) 1 + g 0 1 + g 0 e α ( T ) n
This equation has the form of Equation (A7) and that prediction is just what is observed in Figure 6, where the fit between the above formula and the experimental data is very good ( r 2 = 0.983 ). The results obtained for the parameters in the sample set are g 0 = 0.206 ± 0.047 and α = 0.127 ± 0.028 .

3.3. Load-Unload Tests

In addition to relaxation tests and fast loading cycle tests, three additional loading and unloading tests were performed in the form of triangular waves at very different strain rates, corresponding to stages 4a, 4b and 4c performed at constant stretching speeds λ ˙ t = 0.01 , 0.1 and 1 s 1 . The main purpose was to verify whether the characteristic times ( τ k ) and contribution coefficients ( g k ) obtained from the relaxation tests could well explain the shape of these curves. The displacement had a very roughly triangular shape (except for brief periods of acceleration and deceleration), each triangular wave included about 350 measurements along the loading and unloading cycle ( N p = 285, 365 and 385 points, respectively), that allowed making a numerical integration of the viscoelastic contribution in (6), using Simpson’s 3/8 rule [25]:
I k : = 0 t = N p h e τ / τ k ε ˙ x ( τ ) f ( τ ) d τ = 1 4 k = 0 N p 3 3 h 8 f ( t k ) + 3 f ( t k + 1 ) + 3 f ( t k + 2 ) + f ( t k + 3 )
where h = t / N p and t k = k h . The result can be seen in Figure 7 where the measured and model-predicted force curves are shown. The fit is reasonable although the medium and low velocity curves show larger divergences from the model.
This could be partly due to the fact that the integration error is larger for the slower loading curves. This happens because the number of points is similar in the three curves, but the duration T w of the triangular wave is different, the numerical error depends on the time step h used for the numerical integration, which is given by the relation h = T w / N p ( N p being similar for the three curves). The error resulting from Simpson’s 3/8 rule of integration is given by:
e 3 h 5 80 N p k max ξ k ( t k 1 , t k ) d 4 f ( ξ k ) d τ 4 < C 3 T w 5 80 N p 4 ε ˙
So, for the faster curve a smaller error was made, while for the slower curve the integration error was larger (having a similar number of points, but the duration of the last wave being longer, which makes the division coarser h = T / N p ).

4. Discussion

The mechanical behavior of PSBVs has been previously studied by some researchers [5,7,9,26], although most of the previous studies do not explicitly consider viscoelastic effects in the stress/force–strain curves [10]. Some strain rate related effects are reported, especially with respect to the influence of strain rate on the mechanical failure of PSBVs [3,4,6,7]. Interestingly, like other biological tissues, PSBVs appear to exhibit failure stresses negatively correlated with strain rate [10,27].
However, the relaxation and fast-cycle tests performed in this study clearly show a viscoelastic effect in the PSBV response, especially in stress relaxation stages under constant strain, where the force on the PSBVs is progressively reduced. This fact had not been investigated in detail in the literature on mechanical properties of PSBV [9], although viscoelastic effects had been measured in rapid load cycles [10].
Regarding the results obtained, this work is possibly the first study of PSBVs estimating the necessary relaxation times of the Prony series, see Equation (4), for modeling the viscoelastic behavior of PSBVs. Interestingly, excellent fits ( r 2 > 0.99 ) are obtained with three relaxation times and these differ by an order of magnitude, as τ 1 = O ( 10 1 ) , τ 2 = O ( 10 0 ) , τ 3 = O ( 10 1 ) , so each of them seems to capture a different time scale. This is in agreement with the findings of Funk et al. (2000) [12], where the viscoelasticity of ankle ligaments is studied and where the observed scales are τ 1 = O ( 10 2 ) , τ 2 = O ( 10 1 ) , τ 3 = O ( 10 0 ) . On the other hand, in Davis & De Vita (2012) [13], it is pointed out that the QLVE model applied to rat tendons only captures the behavior well for t < 10 s, which is why they propose a nonlinear viscoelastic constitutive model that overcomes the shortcomings of QLVE to explain their data.
As for the analysis of the preconditioning, which is based on the USF given by Equation (9), an analysis based on viscoelastic behavior has been presented here, which allows predicting that this magnitude will vary from one cycle to another by means of Equation (10). With only two parameters, this formula predicts a better fit than the exponential type heuristic formula:
USF n = a 0 + k = 1 3 a k e b k n
used in [10] and previously in [28,29,30]. Although the Formula (13) is numerically adequate, it is a heuristic formula and is not directly derived from the viscoelastic behavior equations unlike the Formula (10).
As for the applications of this work, they can be grouped into three different areas:
  • Clinicians are charged with the significant task of distinguishing between accidental and inflicted head trauma. Some times this distinction is straightforward, but in many cases the probabilities of injuries from accidental scenarios are unknown, making the differential diagnosis difficult [31]. A refinement of the knowledge of the tolerance ranges against rupture of PSBVs may provide greater accuracy in the reconstruction of injury mechanisms.
  • Computational biomechanics can simulate many potentially traumatic situations, so that models already allow detailed reconstructions of the sequence of events leading to a severe SDH. Accurate knowledge of the material behavior can improve FEHMs, as their inaccuracy is often not so much a computational problem, but a poor characterization of the biomechanical behavior of brain structures. For example, a good number of FEHMs use a stress-strain response for PSBVs that does not reflect the measured nonlinear behavior [3], e.g., the UDS FEHM (Université de Strasbourg) [32], the KTH FEHM (KTH Royal Institute of Technology) [33], UCDBTM (University College Dublin) [34,35], WSUBIM (Wayne State University) [36] or G/LHM [37] also model PSBVs as elastic beams with a linear stress-strain response [9]. The recognition of the importance of the nonlinear behavior and viscoelasticity of brain structures has been explicitly pointed out by the developers of YEAHM (University of Aveiro) [38] and interestingly, some FEHMs model PSBVs as nonlinear elastic materials [39]. In particular, Equation (4), together with the averages obtained from Table 2, allow us to compute an estimation of the viscoelastic effect, independent of the starting elastic model for the PSBVs, which is being used in the FEHM.
  • The improvement of injury metrics used to assess restraint systems in vehicles or the design of other preventive elements against head trauma. Currently, the estimation is often done by the injury metric called “relative motion damage measure” (RMDM) [40,41], used to predict the probability of a SDH due to the failure of PSBVs [42,43]. However, that metric was developed based on obsolete data [26], and the data from this study can be used to update that injury metric.
Finally, some current limitations of the present study that could be improved in future work:
  • The sample used is consistent but small. So, the effects of age, gender, or other anthropometric characteristics on the viscoelastic mechanical properties could not be determined.
  • In addition, a QLVE model has been used in which the relaxation function is separable, in the sense of [13]. Given the low strain levels used for the tests (since care was taken not to cause irreversible damage to the specimen from one test stage to the next), no effects of non-separability were found. However, further work could build a somewhat more general model on that basis. In any case, the proposed model is a first approximation that even explains the data that were not used for the fits (see Section 3.3).
  • Moreover, a further extension of this work would be to examine whether the relaxation curves could be modeled, using the Prony series of stretched exponential relaxation [44]. This could lead to series with fewer and/or more accurate terms, although it is not clear if this is the case. Further work is needed in order to determine whether the use of stretched exponentials or a more general non-separable viscoelastic model would provide better models.

5. Conclusions

The results reported in the current study are the first complete characterization of the viscoelastic response of PSBVs. Relaxation times and coefficients of the Prony series are computed allowing the use of these values in combination with any hyper-elastic model for PSBVs. The model explains accurately all the load cases considered.
In addition, the fitted constitutive parameters can be used for computational FEHMs; these models and the obtained values could be used to reconstruct the circumstances of head trauma and for assessing the causes of TBI. Finally, the results of this study are highly relevant to assess acute SDH due to the mechanical failure of PSBVs.

Author Contributions

Conceptualization, D.S.-M., S.G.-V. and C.A.-D.; methodology, formal analysis, project administration, supervision, validation, and visualization, D.S.-M. and S.G.-V.; software, data curation, S.G.-V. and J.L.; investigation, D.S.-M., S.G.-V., J.L. and I.G.; funding acquisition, C.A.-D. and D.S.-M.; resources I.G., J.V.-A., J.L. and M.C.R.-S.; writing—original draft preparation, D.S.-M. and S.G.-V.; writing—review and editing, I.G. and J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fundació La Marató de TV3 (Barcelona, No. Grant 289/C/2017, 201704.30). The research was developed within the studies of the research project “Age Effects on Traumatic Brain Injuries”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Due to the restrictions imposed by the legal collaboration agreement between IMLFC and UPC, the data are not publicly available, although they can be accessed upon non-anonymous request to the authors of the study and under the conditions set by the aforementioned legal collaboration agreement.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PSBV(s)Parasagittal Bridging Vein(s)
FEHMFinite Element Head Model
QLVEQuasi-Linear Viscoelastic
SDHSubdural Hematoma
TBITraumatic Brain Injury
VCViscoelastic Contribution

Appendix A. Prediction of the USF

In this appendix. we deduce what behavior is to be expected from USF ( n , T ) , defined in the Equation (9), and we justify mathematically why we expect a behavior given by Formula (10). The displacement δ t in the cyclic test is periodic, therefore, it can be expressed by a trigonometric Fourier series:
δ t = Δ δ a 0 2 + m = 1 a m cos 2 m π T t + m = 1 b m sin 2 m π T t
where T = 0.55 s is the period of the cyclic stretch. The stretch λ t , its corresponding rate λ ˙ t , the strain ε x and the strain rate ε ˙ x can be calculated with the following formulas:
λ t = 1 + δ t 0 , λ ˙ t = δ ˙ t 0 , ε x = λ t 2 1 2 , ε ˙ x = λ t λ ˙ t
To estimate the viscoelastic effect we consider a finite sequence of integrals of the form
J n ( k ) : = 0 t n e t n τ τ k δ ˙ τ 0 1 + δ τ 0 d τ
where t n = ( n 1 + γ ) T are the times for which successive force peaks occur:
J n ( k ) ( T ) = m = 1 a m P 1 , m ( θ k ) + b m P 2 , m ( θ k ) + [ a m P 3 , m ( θ k ) + b m P 4 , m ] e ( n 1 + γ ) T τ k ( 1 + 4 m 2 θ k 2 ) ( 1 + 16 m 2 θ k 2 )
where θ k = π T / τ k has been introduced for brevity, and further we have P j , m ( θ k ) as easily computable fourth-degree polynomials, whose coefficients depend on cos ( m γ π ) and sin ( m γ π ) . The above formula allows us to compute successive force peaks and, thus,
USF ( n , T ) = 1 + k = 1 N g k J 1 ( k ) ( T ) 1 + k = 1 N g k J n ( k ) ( T )
where the g k are the same coefficients that appear in the Equation (6) and j is the index for which | T 1 , 5 τ j | is minimum. It can be seen that the value of USF approaches a maximum when the θ j ( 4 , 5 ) (the exact position of the maximum depends on the g k , although very little). Taking the Equation (A4) into consideration, the expression USF ( n , T ) would have the form
USF ( n , T ) = 1 + k = 1 N g k [ c ( θ k ) + d ( θ k ) ] 1 + k = 1 N g k [ c ( θ k ) + d ( θ k ) e n T τ k ] = 1 + k = 1 N g ˜ k 1 + k = 1 N g ˜ k e n T τ k
where c ( θ k ) and d ( θ k ) are functions which can be computed from Equation (A4), and where the g ˜ k are given by
g ˜ k = g k d ( θ k ) 1 + k = 1 N g k c ( θ k )
The expression (A6) is exact, but inconvenient, in fact, numerically it can be seen that the summand for which it is satisfied that | T 1.5 τ j | is minimal, so a reasonable approximation can be obtained using only one characteristic time:
USF ( n , T ) 1 + g j J 1 ( j ) ( T ) 1 + g j J n ( j ) ( T ) = 1 + g ˜ j 1 + g ˜ j exp n T τ j

References

  1. Dewan, M.C.; Rattani, A.; Gupta, S.; Baticulon, R.E.; Hung, Y.C.; Punchak, M.; Agrawal, A.; Adeleye, A.O.; Shrime, M.G.; Rubiano, A.M.; et al. Estimating the global incidence of traumatic brain injury. J. Neurosurg. 2018, 130, 1080–1097. [Google Scholar] [CrossRef] [Green Version]
  2. Lee, M.C.; Haut, R.C. Insensitivity of tensile failure properties of human bridging veins to strain rate: Implications in biomechanics of subdural hematoma. J. Biomech. 1989, 22, 537–542. [Google Scholar] [CrossRef]
  3. Sánchez-Molina, D.; García-Vilana, S.; Llumà, J.; Galtés, I.; Velázquez-Ameijide, J.; Rebollo-Soria, M.C.; Arregui-Dalmases, C. Mechanical behaviour of blood vessels during a Traumatic Brain Injury: Elastic and viscoelastic contributions. Biology 2021, 10, 831. [Google Scholar] [CrossRef]
  4. Delye, H.; Goffin, J.; Verschueren, P.; Vander Sloten, J.; Van der Perre, G.; Alaerts, H.; Verpoest, I.; Berckmans, D. Biomechanical properties of the superior sagittal sinus-bridging vein complex (No. 2006-22-0024). Sae Tech. Pap. 2006, 50, 1–12. [Google Scholar]
  5. Monson, K.L.; Converse, M.I.; Manley, G.T. Cerebral blood vessel damage in traumatic brain injury. Clin. Biomech. 2019, 64, 98–113. [Google Scholar] [CrossRef]
  6. Monson, K.L.; Goldsmith, W.; Barbaro, N.M.; Manley, G.T. Axial mechanical properties of fresh human cerebral blood vessels. J. Biomech. Eng. 2003, 125, 288–294. [Google Scholar] [CrossRef] [PubMed]
  7. Monea, A.G.; Baeck, K.; Verbeken, E.; Verpoest, I.; Vander Sloten, J.; Goffin, J.; Depreitere, B. The biomechanical behaviour of the bridging vein–superior sagittal sinus complex with implications for the mechanopathology of acute subdural haematoma. J. Mech. Behav. Biomed. Mater. 2014, 32, 155–165. [Google Scholar] [CrossRef] [PubMed]
  8. Lee, M.C.; Haut, R.C. Strain rate effects on tensile failure properties of the common carotid artery and jugular veins of ferrets. J. Biomech. 1992, 25, 925–927. [Google Scholar] [CrossRef]
  9. Famaey, N.; Cui, Z.Y.; Musigazi, G.U.; Ivens, J.; Depreitere, B.; Verbeken, E.; Vander Sloten, J. Structural and mechanical characterisation of bridging veins: A review. J. Mech. Behav. Biomed. Mater. 2015, 41, 222–240. [Google Scholar] [CrossRef] [PubMed]
  10. Pasquesi, S.A.; Margulies, S.S. Failure and fatigue properties of immature human and porcine parasagittal bridging veins. Ann. Biomed. Eng. 2017, 45, 1877–1889. [Google Scholar] [CrossRef]
  11. Wang, Z.; Golob, M.J.; Chesler, N.C. Viscoelastic properties of cardiovascular tissues. Viscoelastic Viscoplastic Mater. 2016, 2, 64. [Google Scholar]
  12. Funk, J.R.; Hall, G.W.; Crandall, J.R.; Pilkey, W.D. Linear and quasi-linear viscoelastic characterization of ankle ligaments. J. Biomech. Eng. 2000, 122, 15–22. [Google Scholar] [CrossRef] [PubMed]
  13. Davis, F.M.; De Vita, R. A nonlinear constitutive model for stress relaxation in ligaments and tendons. Ann. Biomed. Eng. 2012, 40, 2541–2550. [Google Scholar] [CrossRef] [PubMed]
  14. Pascoletti, G.; Di Nardo, M.; Fragomeni, G.; Barbato, V.; Capriglione, T.; Gualtieri, R.; Talevi, R.; Catapano, G.; Zanetti, E.M. Dynamic characterization of the biomechanical behaviour of bovine ovarian cortical tissue and its short-term effect on ovarian tissue and follicles. Materials 2020, 13, 3759. [Google Scholar] [CrossRef]
  15. Miller, J.D.; Nader, R. Acute subdural hematoma from bridging vein rupture: A potential mechanism for growth. J. Neurosurg. 2014, 120, 1378–1384. [Google Scholar] [CrossRef]
  16. Costa, J.M.; Fernandes, F.A.; de Sousa, R.J.A. Prediction of subdural haematoma based on a detailed numerical model of the cerebral bridging veins. J. Mech. Behav. Biomed. Mater. 2020, 111, 103976. [Google Scholar] [CrossRef]
  17. Fung, Y.C. Biomechanics Mechanical Properties of Living Tissues; Springer Science & Business Media: New York, NY, USA, 2013. [Google Scholar]
  18. Jin, Z.H. Some notes on the linear viscoelasticity of functionally graded materials. Math. Mech. Solids 2006, 11, 216–224. [Google Scholar] [CrossRef]
  19. Soussou, J.E.; Moavenzadeh, F.; Gradowczyk, M.H. Application of Prony series to linear viscoelasticity. Trans. Soc. Rheol. 1970, 14, 573–584. [Google Scholar] [CrossRef]
  20. Tzikang, C. Determining a Prony Series for a Viscoelastic Material from Time Varying Strain Data; NASA Center for Aero Space Information (CASI): Hampton, VA, USA, 2000.
  21. Tapia-Romero, M.A.; Dehonor-Gómez, M.; Lugo-Uribe, L.E. Prony series calculation for viscoelastic behavior modeling of structural adhesives from DMA data. Ing. Investig. Tecnol. 2020, 21, 1–10. [Google Scholar] [CrossRef]
  22. Sun, Z.; Gepner, B.D.; Lee, S.H.; Rigby, J.; Cottler, P.S.; Hallman, J.J.; Kerrigan, J.R. Multidirectional mechanical properties and constitutive modeling of human adipose tissue under dynamic loading. Acta Biomater. 2021, 129, 188–198. [Google Scholar] [CrossRef]
  23. Braess, D. Finite Elements: Theory, Fast Solvers, and Applications in Solid Mechanics; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  24. Iriarte, Y.A.; Varela, H.; Gómez, H.J.; Gómez, H.W. A Gamma-Type distribution with applications. Symmetry 2020, 12, 870. [Google Scholar] [CrossRef]
  25. Al-A’asam, J.A. Deriving the composite Simpson rule by using Bernstein polynomials for solving Volterra integral equations. Baghdad Sci. J. 2014, 11, 1274–1282. [Google Scholar]
  26. Löwenhielm, P. Dynamic properties of the parasagittal bridging beins. Z. FüR Rechtsmed. 1974, 74, 55–62. [Google Scholar] [CrossRef]
  27. Sánchez-Molina, D.; García-Vilana, S.; Velázquez-Ameijide, J.; Arregui-Dalmases, C. Probabilistic assessment for clavicle fracture under compression loading: Rate-dependent behavior. Biomed. Eng. Appl. Basis Commun. 2020, 32, 2050040. [Google Scholar] [CrossRef]
  28. Gefen, A.; Margulies, S.S. Are in vivo and in situ brain tissues mechanically similar? J. Biomech. 2004, 37, 1339–1352. [Google Scholar] [CrossRef]
  29. Komatsu, K.; Sanctuary, C.; Shibata, T.; Shimada, A.; Botsis, J. Stress-relaxation and microscopic dynamics of rabbit periodontal ligament. J. Biomech. 2007, 40, 634–644. [Google Scholar] [CrossRef]
  30. Zhang, K.; Siegmund, T.; Chan, R.W. Modeling of the transient responses of the vocal fold lamina propria. J. Mech. Behav. Biomed. Mater. 2009, 2, 93–104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Coats, B.; Eucker, S.A.; Sullivan, S.; Margulies, S.S. Finite element model predictions of intracranial hemorrhage from non-impact, rapid head rotations in the piglet. Int. J. Dev. Neurosci. 2012, 30, 191–200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Raul, J.S.; Roth, S.; Ludes, B.; Willinger, R. Influence of the benign enlargement of the subarachnoid space on the bridging veins strain during a shaking event: A finite element study. Int. J. Leg. Med. 2008, 122, 337–340. [Google Scholar] [CrossRef]
  33. Kleiven, S. Predictors for traumatic brain injuries evaluated through accident reconstructions (No. 2007-22-0003). SAE Tech. Pap. 2007, 51, 1–35. [Google Scholar]
  34. Doorly, M.C. Investigations into Head Injury Criteria Using Numerical Reconstruction of Real Life Accident Cases. Ph.D. Thesis, University College Dublin, Dublin, Ireland, 2007. [Google Scholar]
  35. Yan, W.; Pangestu, O.D. A modified human head model for the study of impact head injury. Comput. Methods Biomech. Biomed. Eng. 2011, 14, 1049–1057. [Google Scholar] [CrossRef] [PubMed]
  36. Viano, D.C.; Casson, I.R.; Pellman, E.J.; Zhang, L.; King, A.I.; Yang, K.H. Concussion in professional football: Brain responses by finite element analysis: Part 9. Neurosurgery 2005, 57, 891–916. [Google Scholar] [CrossRef]
  37. Zoghi-Moghadam, M.; Sadegh, A.M. Global/local head models to analyse cerebral blood vessel rupture leading to ASDH and SAH. Comput. Methods Biomech. Biomed. Eng. 2009, 12, 1–12. [Google Scholar] [CrossRef]
  38. Fernandes, F.A.; Tchepel, D.; de Sousa, R.J.A.; Ptak, M. Development and validation of a new finite element human head model: Yet another head model (YEAHM). Eng. Comput. 2018, 35, 477–496. [Google Scholar] [CrossRef]
  39. Mao, H.; Zhang, L.; Jiang, B.; Genthikatti, V.V.; Jin, X.; Zhu, F.; Yang, K.H. Development of a finite element human head model partially validated with thirty five experimental cases. J. Biomech. Eng. 2013, 135, 111002. [Google Scholar] [CrossRef] [PubMed]
  40. Sanchez-Molina, D.; Velazquez-Ameijide, J.; Arregui-Dalmases, C.; Crandall, J.R.; Untaroiu, C.D. Minimization of analytical injury metrics for head impact injuries. Traffic Inj. Prev. 2012, 13, 278–285. [Google Scholar] [CrossRef]
  41. Sánchez-Molina, D.; Arregui-Dalmases, C.; Velázquez-Ameijide, J.; Angelini, M.; Kerrigan, J.; Crandall, J. Traumatic brain injury in pedestrian–vehicle collisions: Convexity and suitability of some functionals used as injury metrics. Comput. Methods Programs Biomed. 2016, 136, 55–64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Takhounts, E.G.; Ridella, S.A.; Hasija, V.; Tannous, R.E.; Campbell, J.Q.; Malone, D.; Duma, S. Investigation of traumatic brain injuries using the next generation of simulated injury monitor (SIMon) finite element head model (No. 2008-22-0001). SAE Tech. Pap. 2008, 52, 1. [Google Scholar]
  43. Fernandes, F.A.; Sousa, R.J.A.D. Head injury predictors in sports trauma—A state-of-the-art review. Proc. Inst. Mech. Eng. Part J. Eng. Med. 2015, 229, 592–608. [Google Scholar] [CrossRef]
  44. Mauro, J.C.; Mauro, Y.Z. On the Prony series representation of stretched exponential relaxation. Phys. Stat. Mech. Appl. 2018, 506, 75–87. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Experimental setting measurements: the coordinate system is centered on the static fixture with the x axis parallel to the load direction, and a position change from the initial X to the deformed x t configuration is represented. (b) Cyclic and relaxation test stages scheme: (1) load-holding to maximum strain ( ε max ), (2) fast load/unload cycle stage to maximum strain, (3) load-holding to ε max , 2/3 · ε max and 1/3 · ε max , (4) load and unload to ε max at 0.01 , 0.1 , 1 s 1 and (5) load to fracture.
Figure 1. (a) Experimental setting measurements: the coordinate system is centered on the static fixture with the x axis parallel to the load direction, and a position change from the initial X to the deformed x t configuration is represented. (b) Cyclic and relaxation test stages scheme: (1) load-holding to maximum strain ( ε max ), (2) fast load/unload cycle stage to maximum strain, (3) load-holding to ε max , 2/3 · ε max and 1/3 · ε max , (4) load and unload to ε max at 0.01 , 0.1 , 1 s 1 and (5) load to fracture.
Bioengineering 08 00145 g001
Figure 2. Imposed displacements in the different test stages and reaction force of the PSBV for specimen 2632A/20. The stages with constant displacement have been enlarged to show the decrease in force in stages (1) to (3), and the concavity of the load–discharge curve of stage (4).
Figure 2. Imposed displacements in the different test stages and reaction force of the PSBV for specimen 2632A/20. The stages with constant displacement have been enlarged to show the decrease in force in stages (1) to (3), and the concavity of the load–discharge curve of stage (4).
Bioengineering 08 00145 g002
Figure 3. Experimental displacement and experimental velocity applied to a specimen, plotted together with the fit by gamma-type function. The plot also shows the derivatives in Section 1 and Section 3. It can be seen that the fit is in both cases more extreme than 0.997.
Figure 3. Experimental displacement and experimental velocity applied to a specimen, plotted together with the fit by gamma-type function. The plot also shows the derivatives in Section 1 and Section 3. It can be seen that the fit is in both cases more extreme than 0.997.
Bioengineering 08 00145 g003
Figure 4. Fitting of the proposed viscoelastic model, for stage (1) of one of the specimens. Fits with N = 1 , N = 2 and N = 3 terms in the Prony series are considered, with the corresponding parameters from Table 1.
Figure 4. Fitting of the proposed viscoelastic model, for stage (1) of one of the specimens. Fits with N = 1 , N = 2 and N = 3 terms in the Prony series are considered, with the corresponding parameters from Table 1.
Bioengineering 08 00145 g004
Figure 5. Evolution of VC in strength over time for the analyzed specimens: the black lines indicate the average elastic and viscoelastic contributions as a function of time, the shaded area indicates the error interval in the sample.
Figure 5. Evolution of VC in strength over time for the analyzed specimens: the black lines indicate the average elastic and viscoelastic contributions as a function of time, the shaded area indicates the error interval in the sample.
Bioengineering 08 00145 g005
Figure 6. Force peaks in stage (2) as a function of the number of cycles performed and USF computed and predicted by Equation (10).
Figure 6. Force peaks in stage (2) as a function of the number of cycles performed and USF computed and predicted by Equation (10).
Bioengineering 08 00145 g006
Figure 7. Experimental force–time curves compared with model predictions: (a) fast curve ( ε ˙ 1 s 1 ), (b) medium velocity curve ( ε ˙ 0.1 s 1 ) and (c) slow curve ( ε ˙ 0.01 s 1 ). The predicted curve has been obtained by numerical integration. The cumulative integration error is larger for the slow curve and smaller for the fast curve. However, the fast curve seems to have a larger measurement error. Possibly, for that reason, the best fit is obtained for the medium speed curve.
Figure 7. Experimental force–time curves compared with model predictions: (a) fast curve ( ε ˙ 1 s 1 ), (b) medium velocity curve ( ε ˙ 0.1 s 1 ) and (c) slow curve ( ε ˙ 0.01 s 1 ). The predicted curve has been obtained by numerical integration. The cumulative integration error is larger for the slow curve and smaller for the fast curve. However, the fast curve seems to have a larger measurement error. Possibly, for that reason, the best fit is obtained for the medium speed curve.
Bioengineering 08 00145 g007
Table 1. Parameters g k , τ k of the viscoelastic model obtained for the specimen 635/21B for different values of k.
Table 1. Parameters g k , τ k of the viscoelastic model obtained for the specimen 635/21B for different values of k.
N g 1 [–] g 2 [–] g 3 [–] τ 1 [s] τ 2 [s] τ 3 [s]
N = 1 0.36010.756
N = 2 0.2430.39313.3970.799
N = 3 0.2150.2210.39818.1621.8850.211
Table 2. Average ± standard deviation of the parameters g 1 , g 2 , g 3 and τ 1 , τ 2 , τ 3 of the viscoelastic model in the four relaxation stages obtained for each specimen.
Table 2. Average ± standard deviation of the parameters g 1 , g 2 , g 3 and τ 1 , τ 2 , τ 3 of the viscoelastic model in the four relaxation stages obtained for each specimen.
Specimen g 1 [–] g 2 [–] g 3 [–] τ 1 [s] τ 2 [s] τ 3 [s]
2632A0.10 ± 0.050.13 ± 0.041.70 ± 0.2910.96 ± 2.640.88 ± 0.190.07 ± 0.01
617A0.21 ± 0.100.32 ± 0.160.34 ± 0.0615.38 ± 4.901.99 ± 2.620.17 ± 0.19
617B0.24 ± 0.070.13 ± 0.080.30 ± 0.2313.78 ± 9.150.96 ± 0.700.09 ± 0.08
621A0.20 ± 0.070.18 ± 0.100.48 ± 0.4217.08 ± 16.111.20 ± 1.230.09 ± 0.10
635A0.12 ± 0.070.13 ± 0.070.84 ± 1.1320.73 ± 13.342.16 ± 1.460.20 ± 0.10
635B0.14 ± 0.060.16 ± 0.081.60 ± 1.6225.41 ± 8.221.76 ± 1.930.17 ± 0.20
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

García-Vilana, S.; Sánchez-Molina, D.; Llumà, J.; Galtés, I.; Velázquez-Ameijide, J.; Rebollo-Soria, M.C.; Arregui-Dalmases, C. Viscoelastic Characterization of Parasagittal Bridging Veins and Implications for Traumatic Brain Injury: A Pilot Study. Bioengineering 2021, 8, 145. https://doi.org/10.3390/bioengineering8100145

AMA Style

García-Vilana S, Sánchez-Molina D, Llumà J, Galtés I, Velázquez-Ameijide J, Rebollo-Soria MC, Arregui-Dalmases C. Viscoelastic Characterization of Parasagittal Bridging Veins and Implications for Traumatic Brain Injury: A Pilot Study. Bioengineering. 2021; 8(10):145. https://doi.org/10.3390/bioengineering8100145

Chicago/Turabian Style

García-Vilana, Silvia, David Sánchez-Molina, Jordi Llumà, Ignasi Galtés, Juan Velázquez-Ameijide, M. Carmen Rebollo-Soria, and Carlos Arregui-Dalmases. 2021. "Viscoelastic Characterization of Parasagittal Bridging Veins and Implications for Traumatic Brain Injury: A Pilot Study" Bioengineering 8, no. 10: 145. https://doi.org/10.3390/bioengineering8100145

APA Style

García-Vilana, S., Sánchez-Molina, D., Llumà, J., Galtés, I., Velázquez-Ameijide, J., Rebollo-Soria, M. C., & Arregui-Dalmases, C. (2021). Viscoelastic Characterization of Parasagittal Bridging Veins and Implications for Traumatic Brain Injury: A Pilot Study. Bioengineering, 8(10), 145. https://doi.org/10.3390/bioengineering8100145

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