Stress Relaxation Tests Research Articles

Enhancing Hole Transport and Autonomous Healing Properties of Supramolecular Columns in Unsymmetrical Discotics.

Designing smart autonomous healing soft materials is crucial to attaining cost-efficiency and optimal performance in organic semiconductors. In this context, we design an unsymmetrical thiophene-fused phenazine (TFP)-based discotic liquid crystal (DLC) with the goal of creating an active organic semiconductor that encompasses favorable attributes, such as polarizability, mobility, and processability. Aligned with our objective, we successfully synthesized two unsymmetrical TFP core-based DLCs by linking alkyl chains of variable lengths at the periphery through a coupling reaction. These DLCs show room temperature columnar oblique (Colob) mesophase as evident from temperature-dependent studies employing polarized optical microscopy (POM) and small-angle X-ray scattering (SAXS). We performed space charge-limited current (SCLC) techniques to evaluate the hole mobility of TFP-based DLCs and found that they have high hole mobility (∼10-3 cm2/V s). Additionally, the film morphology and its self-healing nature have been examined using atomic force microscopy (AFM) and stress relaxation test. One of the unsymmetrical DLCs exhibited an ability to relax stress over time while maintaining a constant strain (up to 1-3%). The rational design of these unsymmetrical discotics, exhibiting reduced threshold voltage (as confirmed by conductivity studies) highlights their possible potential in various organic semiconductor devices.

Advancing mechanical testing of biological tissues and hydrogels: A buoy-based approach

Quantifying the material properties of tissues and hydrogels aids in the development of biomedical applications through better understanding of the mechanics and mechanobiological principles at play. This study introduces a mechanical testing platform designed to address challenges in measuring mm-scale tissue and hydrogel material properties. Using a floating buoy design, the platform enables horizontal submerged tensile testing with non-submersible load cells. Buoy drag testing without a sample attached resulted in signal noise (mean ± standard deviation) of −1.6E-4 ± 2.8E-2 mN for stationary recording and 1.3E-2 ± 6.7E-2 mN for maximum buoy displacement speed (1000 µm/s), suggesting the magnitude of drag forces from buoy movement are negligible in comparison to the minimum resolution of force measurement. Validation testing with latex orthodontic bands showed a ∼ 28 × reduction in signal noise with the buoy approach compared to a previously used approach, and similar force displacement recordings using the buoy approach with 2 separate hardware systems. Simultaneous imaging enabled geometrical and microstructural analysis of the sample. Murine bladder tissue was mechanically tested using two different hardware systems and testing protocols. The platform was able to accurately capture the nonlinear stress-stretch response, alongside expected strain-softening and preconditioning behavior. Stress-relaxation testing yielded results consistent with expected microstructural and viscoelastic responses of mouse bladder tissue. The versatility of the platform underscores the potential for it to integrate with various force measurement and actuator-based systems. In conclusion, this platform offers a new avenue for accurate measurement of tissue and hydrogel mechanics, facilitating mm-scale soft material research.

Modeling and Simulation of Curing Kinetics and Stress of a Bonded Mirror with a Room-Temperature Curable Adhesive

Adhesive bonding widely appears in optomechanical products. Nevertheless, the occurrence of curing stress often leads to structure deformation, which significantly compromises the performance of optomechanical systems. This work focuses on the formation and evolution of curing stress and its influence on a multi-point face-bonded mirror. The adhesive’s degree of cure (DOC) was measured using the differential scanning calorimetry (DSC) technique. The functional relationship between DOC and glass-transition temperature ([Formula: see text]) was determined utilizing DiBenedetto’s equation. The time–DOC shift factor of the adhesive was obtained by stress relaxation tests on specimens with different DOCs. Subsequently, a viscoelastic constitutive model concerning the curing process and an accurate finite-element (FE) model of the bonded mirror were developed to simulate the surface-figure error (SFE) evolution. The predicted surface figure was compared with experiments, and the model was verified and validated. Further, simulation analyses were conducted on the performance of the multi-point face-bonded mirror, and the influences of adhesive thickness, bonding area, and spot number on the SFE of the mirror were examined.

Flash Sintering of Bimetallic Assemblies for Turbine Disks in Next‐Generation Jet Engines

René 77 and another advanced powder metallurgy nickel‐based superalloy are flash‐sintered and joined to produce a dual‐alloy assembly. This innovative process allows for the sintering of powders to achieve high relative densities or the formation of solid joints in a shorter time than is feasible with conventional diffusion bonding techniques. Examination of the microstructure after heat treatments reveals that diffusion and recrystallization phenomena can modulate the transition zone. The creep properties of both alloys, evaluated using stress relaxation tests, do not seem to differ by more than an order of magnitude, while the higher thermal stability of René 77 compared to Alloy B at 800 °C tips the scales. In tension, base materials show very high properties, while bimetallic specimens exhibit at least higher properties than René 77. However, failure occurs at the bond at high temperatures due to the formation of an oxide film during the bonding process. These encouraging results illustrate the potential and constraints of the flash‐sintering process for producing dissimilar assemblies, which may enhance the properties of next‐generation aero‐engine turbine disks.

Polymerization-Induced Self-Assembly Providing PEG-Gels with Dynamic Micelle-Crosslinked Hierarchical Structures and Overall Improvement of Their Comprehensive Performances.

Polymer gels are fascinating soft materials and have become excellent candidates for wearable electronics, biomedicine, sensors, etc. Synthetic gels usually suffer from poor mechanical properties, and integrating good mechanical properties, adhesiveness, stability, and self-healing performances in one gel is more difficult. Herein, polymerization-induced self-assembly (PISA) providing PEG-gels with an overall improvement in their comprehensive performances is reported. PISA synthesis is carried out in PEG (solvent) to efficiently produce various nanoparticles, which are used as the nanofillers in the subsequent synthesis of PEG-gels with dynamic micelle-crosslinked hierarchical structures. Compared to hydrogels, PEG-gels show excellent long-term stability due to the nonvolatile feature of PEG solvent. The hierarchical PEG-gels (with nanofillers) exhibit better mechanical and adhesive properties than the homogeneous-gels (without nanofillers). The energy dissipation mechanism of the PEG-gels is analyzed via stress relaxation and cyclic mechanical tests. High-density hydrogen bonds between the micelles and PAA matrix can be broken and reformed, endowing better self-healing properties of the dynamic micelle-crosslinked PEG gels. This work provides a simple strategy for producing hierarchical structural gels with enhanced properties, which offers fundamentals and inspirations for the designing of various advanced functional materials.

Vitrimer-like elastomers with rapid stress-relaxation by high-speed carboxy exchange through conjugate substitution reaction

We report vitrimer-like elastomers that exhibit significantly fast stress relaxation using carboxy exchange via the conjugate substitution reaction of α-(acyloxymethyl) acrylate skeletons. This network design is inspired by a small-molecule model that shows the carboxy exchange reaction even at ambient temperature in the presence of 1,4-diazabicyclo[2.2.2]octane (DABCO). The acrylate and acrylic acid copolymers are cross-linked using bis[α-(bromomethyl)acrylates] and doped with 10 wt% DABCO, exhibiting processability to obtain a transparent film by hot pressing. The high-speed bond exchange in the network, validated by stress-relaxation tests, allows quick molding with household iron. In addition, the material is applied as an adhesion sheet for plastic and metal substrates. Because dynamic cross-linking with the proposed bond exchange mechanism can be implemented for any polymer bearing carboxyl pendants, our approach can be applied to versatile backbones, which must thus be meaningful in the practical sense.

Thermomechanical and Viscoelastic Characterization of Continuous GF/PETG Tape for Extreme Environment Applications

The thermomechanical and viscoelastic properties of a glass fiber polyethylene terephthalate glycol (GF/PETG) continuous unidirectional (UD) tape were investigated using differential scanning calorimetry (DSC), thermomechanical analysis (TMA), and dynamic mechanical analysis (DMA). This study identified five operational conditions based on the Army Regulation 70-38 Standard. The DSC results revealed a glass transition temperature of 78.0 ± 0.3 °C, guiding the selection of temperatures for TMA and DMA tests. TMA provided the coefficient of thermal expansion in three principal directions, consistent with known values for PETG and GF materials. DMA tests, including strain sweep, temperature ramp, frequency sweep, creep, and stress relaxation, defined the material’s linear viscoelastic region and temperature-dependent properties. The frequency sweep indicated an increased modulus with rising frequency, identifying several natural frequency modes. Creep and stress relaxation tests showed time-dependent behavior, with strain increasing under higher loads and stress decreasing over time for all tested input values. Viscoelastic models fitted to the data yielded R2 values of 0.99, demonstrating good agreement. The study successfully measured thermomechanical and viscoelastic properties across various conditions, providing insights into how temperature influences the material’s mechanical response under extreme conditions.

Alterations in mechanical properties of rabbit collateral ligaments eight weeks after anterior cruciate ligament transection

Anterior cruciate ligament (ACL) injury is a common knee ligament injury among young, active adults; however, little is known about its impact on the viscoelastic properties of the knee joint’s collateral ligaments. This study aimed to characterize and compare the viscoelastic properties of rabbit collateral ligaments in healthy control knees, injured knees, and knees contralateral to the injured knees.Unilateral anterior cruciate ligament transection was performed on six New Zealand white rabbits to create an ACL injury model. Medial and lateral collateral ligaments (MCL and LCL) were collected from the injured and contralateral knees eight weeks after ACL transection. Ligaments were also harvested from both knees of four unoperated rabbits. The ligaments underwent tensile stress-relaxation testing at strain levels of 2, 4, 6, and 8 %, and a sinusoidal loading test at 8 % strain with 0.5 % strain amplitude using frequencies of 0.01, 0.05, 0.1, 0.5, 1, and 2 Hz.The results showed that collateral ligaments of ACL-transected knees relaxed slower compared to control knees. Sinusoidal testing revealed that contralateral knee LCLs had significantly higher storage and loss modulus across all test frequencies.The results indicate that contralateral knee LCLs become stiffer compared to LCLs from control and ACL-transected knees, while LCLs from ACL-transected knees become less viscous compared to LCLs from control and contralateral knees. This study suggests that knee ligaments undergo adaptations following an ACL injury that may affect the mechanics of the ACL-transected knee, which should be considered in biomechanical and rehabilitation studies of patients with an ACL injury.

Simple Aseismic Reinforcement of Steel Structures Using Knee Braces with High-Hardness Vises

A novel technique for upgrading the seismic resistance of steel buildings by adding knee braces to existing structures using vises was proposed by researchers in 2022. A feature of this retrofitting method is the easy setup owing to its use of vises made from high-hardness metal. Tests were conducted to investigate two main failure modes: slipping failure at the connection and yielding and buckling failure of the knee brace. The retrofitting design is discussed based on a comparison between the slipping strengths obtained through tests and calculations. Furthermore, an analytical study, using the finite element method (FEM), was conducted to evaluate the test results of retrofitted frames that failed in terms of the yielding and buckling of the knee braces. The findings of the analyses are consistent with the test results. This study included a stress relaxation test to assess the long-term performance of the vises.

Digital Representation of Materials Testing Data for Semantic Web Analytics: Tensile Stress Relaxation Testing Use Case

This study aims to represent an approach for transferring the materials testing datasets to the digital schema that meets the prerequisites of the semantic web. As a use case, the tensile stress relaxation testing method was evaluated and the testing datasets for several copper alloys were prepared. The tensile stress relaxation testing ontology (TSRTO) was modeled following the test standard requirements and by utilizing the appropriate upper-level ontologies. Eventually, mapping the testing datasets into the knowledge graph and converting the data-mapped graphs to the machine-readable Resource Description Framework (RDF) schema led to the preparation of the digital version of testing data which can be efficiently queried on the web.

Modelling the time-dependent mechanical properties of thermoplastic and thermosetting polymers with Gumbel distribution functions

In this paper, we propose a simple modelling method that accurately describes the time-dependent viscoelastic behaviour of thermoplastic and thermosetting solid polymers. We performed short-term creep and stress relaxation tests on representative samples of three major classes of polymers: an amorphous, a semi-crystalline thermoplastic and a thermosetting polymer. Then, to investigate the relationship between the model parameters and the material composition, we also tested thermoplastic matrix composites with different plasticiser contents. From the short-term tests, master curves were constructed with the use of the time–temperature superposition principle. Then, the temperature dependence of the obtained shift factors was described using the Arrhenius model. The resulting creep compliance and relaxation modulus master curves were normalised and modelled with a two-parameter cumulative Gumbel distribution function (CDF) and its complementary distribution function (CCDF). For the creep and stress relaxation master curves, the median error of modelling was less than 16 % and 16.5 % in all cases, respectively. Furthermore, it was less than 2 % for the thermosetting polymer, making this model an excellent tool for describing the creep and stress relaxation behaviour of thermosetting systems. In the case of the plasticised composites, we found a strong relationship between the material composition and the location parameter of the Gumbel model. Based on this, we developed a method that can be used to estimate the evolution of the creep compliance curve for any arbitrary OLA content (between 0 wt% and 16 wt%) at any desired temperature (between 22 °C and 70 °C).

Assessment of various polymeric sulphur as soft bitumen improver for pavement applications

ABSTRACT In this study, five polymeric sulphurs are utilised, including polymeric sulphur based on both aliphatic copolymers and styrene copolymers, namely PS1 to PS5 correlated to amorphous poly alpha olefins (APAO) PS1, olefin residues from MAC carpet company PS2, Ethylene Vinyl Acetate (EVA) copolymer PS3, Styrene Butadiene Styrene (SBS) copolymer PS4 and Dicyclopentadiene monomer (DCPD) PS5. The five polymeric sulphur compounds were mixed by 25wt% with soft bitumen penetration grade (80/100) to produce five sulphur-extended asphalt binders, namely SEA1, SEA2, SEA3, SEA4, and SEA5. The physical properties, including softening point, viscosity, and penetration, are examined to characterise the consistency of SEA binder. The dynamic mechanical behaviour such as complex modulus, temperature sweep, and stress relaxation test was investigated for SEA binder and compared with a reference soft bitumen (Blank) using a dynamic mechanical analyzer (DMA). It is concluded that polymeric sulphur based on styrene compounds SEA2 and SEA4 binders have the highest resistance against cracking of the pavement at low temperatures according to penetration index readings, and have the least permanent (plastic) deformation at high temperatures according to temperature sweep from DMA data. It is noticed that SEA5 binder has a negative effect on the performance of soft bitumen.

Integrating dynamic relaxation with inelastic deformation in metallic glasses: Theoretical insights and experimental validation

Metallic glasses (MGs), a metastable material far from equilibrium, exhibit intricate dynamic relaxation behaviors. The challenge lies in developing a model that accurately describes the dynamics and deformation mechanisms of MG. This paper introduces a model integrating dynamic relaxation with deformation behavior. Validation through dynamic mechanical analysis, stress relaxation, creep, and strain recovery tests confirm the existence of four deformation modes: elasticity, anelasticity from β relaxation, anelasticity from α relaxation at low temperatures, and viscoplasticity from α relaxation at high temperatures. The model captures all of these deformation modes. The dynamical mechanical spectrum and stress relaxation spectrum unveil dynamic features during glass to liquid transition, and a simple and effective experimental method was developed for identifying ultra-low-frequency dynamic relaxation. This work provides new perspectives on the study of relaxation dynamics in glassy states and establishes important connections between dynamic relaxation behavior and deformation mechanisms. These findings lay a theoretical and experimental foundation for the broad application of MGs.

Reduced warpage in semiconductor packages: Optimizing post-cure temperature profile considering cure shrinkage and viscoelasticity of epoxy molding compound

Warpage in semiconductor packages is a critical issue that affects their reliability and performance. This study aims to minimize the warpage of a bi-material dummy package by optimizing the post-mold curing (PMC) temperature profile. A warpage simulation model was developed considering the viscoelastic properties and cure shrinkage behavior of the epoxy molding compound (EMC). A compression molding monitoring system was also developed to analyze the cure behavior of the EMC. The gelation point of the EMC during compression molding was determined using a dielectric sensor, while cure shrinkage behaviors were analyzed using a fiber Bragg grating sensor. Viscoelastic properties were characterized through stress relaxation tests, and thermal properties were measured using three-dimensional digital image correlation. The developed simulation model was used to optimize the PMC temperature profile using a genetic algorithm. Experimental validation showed that the optimized temperature profile reduced the warpage from ∼1519 μm to 486 μm compared to a linear profile. The findings of this study have broader implications for the design and manufacturing of semiconductor packages, as the proposed approach can be applied to minimize warpage and improve package reliability.

Robust Recovery of Optimally Smoothed Polymer Relaxation Spectrum from Stress Relaxation Test Measurements.

The relaxation spectrum is a fundamental viscoelastic characteristic from which other material functions used to describe the rheological properties of polymers can be determined. The spectrum is recovered from relaxation stress or oscillatory shear data. Since the problem of the relaxation spectrum identification is ill-posed, in the known methods, different mechanisms are built-in to obtain a smooth enough and noise-robust relaxation spectrum model. The regularization of the original problem and/or limit of the set of admissible solutions are the most commonly used remedies. Here, the problem of determining an optimally smoothed continuous relaxation time spectrum is directly stated and solved for the first time, assuming that discrete-time noise-corrupted measurements of a relaxation modulus obtained in the stress relaxation experiment are available for identification. The relaxation time spectrum model that reproduces the relaxation modulus measurements and is the best smoothed in the class of continuous square-integrable functions is sought. Based on the Hilbert projection theorem, the best-smoothed relaxation spectrum model is found to be described by a finite sum of specific exponential-hyperbolic basis functions. For noise-corrupted measurements, a quadratic with respect to the Lagrange multipliers term is introduced into the Lagrangian functional of the model's best smoothing problem. As a result, a small model error of the relaxation modulus model is obtained, which increases the model's robustness. The necessary and sufficient optimality conditions are derived whose unique solution yields a direct analytical formula of the best-smoothed relaxation spectrum model. The related model of the relaxation modulus is given. A computational identification algorithm using the singular value decomposition is presented, which can be easily implemented in any computing environment. The approximation error, model smoothness, noise robustness, and identifiability of the polymer real spectrum are studied analytically. It is demonstrated by numerical studies that the algorithm proposed can be successfully applied for the identification of one- and two-mode Gaussian-like relaxation spectra. The applicability of this approach to determining the Baumgaertel, Schausberger, and Winter spectrum is also examined, and it is shown that it is well approximated for higher frequencies and, in particular, in the neighborhood of the local maximum. However, the comparison of the asymptotic properties of the best-smoothed spectrum model and the BSW model a priori excludes a good approximation of the spectrum in the close neighborhood of zero-relaxation time.

Enhanced stress relaxation behavior via basal 〈a〉 dislocation activity in Zircaloy-4 cladding

This work evaluates the stress relaxation behavior of textured Zircaloy-4 cladding to understand how mechanical anisotropy influences pellet-cladding interactions. Uniaxial and biaxial stress relaxation tests are performed using full-tube axial tension and internal pressurization, respectively, aiming to achieve 0.25 %, 1 %, and 2 % equivalent strains in the cladding samples at a temperature of 300 °C. Internal pressure relaxation test results display enhanced stress relaxation compared to axial testing results, particularly for samples loaded beyond yield. Results of electron backscatter diffraction indicate increased deformation microstructure during loading and increased strain homogenization and recovery during relaxation for samples loaded via internal pressurization. Analysis indicates that the increased production and activity of basal 〈a〉 dislocations play a significant role in the enhanced relaxation measured in samples subjected to internal pressurization.

Abstract Tu048: Viscoelastic remodeling of the left ventricular myocardium in myocardial infarction

Introduction. Myocardial infarction (MI) leads to cardiac myocyte death and the formation of a fibrotic scar in the left ventricular myocardium (LV). LV free wall (LVFW) remodeling in MI is known to significantly alter the biomechanical response of the LV and, in turn, LV function. While changes in the elastic behavior of the LVFW in MI have been reported, changes in the viscoelastic behavior of the LVFW remain understudied. Our objective in this study was to determine the changes in the stress relaxation behavior of the LVFW in MI. Methods. MI was induced in rats via the ligation of the LAD artery, and hearts were harvested for three groups: sham, 2-week (2wk), and 4wk post-MI (n=3 for each group). Prior to sacrifice, echocardiography was performed to estimate ejection fraction (EF). LVFW specimens were harvested and subjected to viscoelastic stress relaxation tests. The specimens were stretched to 15% along the circumferential and longitudinal directions and allowed to relax. The stress was normalized against the peak (initial) stress for each specimen for proper comparison across the groups. Results. The EF at 2wk and 4wk specimens was reduced compared to the sham specimens (Fig. 1A). The 4wk post-MI specimens showed slower stress relaxation compared to that in the sham and 2wk specimens in both directions (Figs. 1B, C). Considering a 2-minute stress relaxation span, the normalized stress in the 4wk specimens was significantly higher than the sham specimens in both directions (Fig. 1D). Conclusions. Our results suggest that the viscoelastic relaxation rate of the LVFW diminishes at late MI timepoints. A larger and more mature scar in the LVFW is expected to contribute to slower relaxation at late MI. Future studies can determine the effect of LVFW viscoelastic remodeling on LV diastolic dysfunction and impaired relaxation in MI.

Significant reduction of anisotropy in stress relaxation aging and mechanical properties improvement for 2195 Al-Cu-Li alloy subjected to plastic loading

Although the plastic loading can enhance creep deformation and yield strength, the anisotropic Stress Relaxation Aging (SRA) behavior and mechanism under plastic loading remain unclear, which presents a significant challenge in accurately shaping aluminum alloy panels. In this study, the SRA behavior of 2195-T4 Al-Cu-Li alloys were thoroughly studied under initial loading stresses within the elastic (210/250 MPa) and plastic (380/420 MPa) ranges at 180 ℃ by stress relaxation and tensile tests as well as microstructure characterization. The findings reveal that compared with those under elastic loadings, in-plane anisotropy (IPA) values of the stress relaxation amount, yield strength and fracture elongation under plastic loadings are reduced by 60%-80%, 70%-90% and 72%-89%, respectively. Similarly, IPA values of precipitate size in grains and Precipitation-Free Zones (PFZ) width at grain boundaries under plastic loading decrease by 31.4% and 94.4% respectively. These results indicate plastic loading significantly weakens the anisotropic SRA behavior, owing to numerous uniformly distributed fine T1 phases and small IPA values of both T1 precipitates size and PFZ width in various loading directions. Compared with those of elastic loading-aged alloys, yield strength of plastic loading-aged alloys shows high strength-ductility because of the combined effect of closely dispersed fine T1 precipitates, narrowed PFZ and numerous sheared and rotated T1 phases at different locations during tensile process. The uniformly distributed larger Kernel Average Misorientation (KAM) and Schmidt factor values of the plastic loading-aged alloy, as well as the cross-slip generated, also help to enhance the strength and ductility of the alloy.

Evaluation of the effects of temperature and centrifugation time on elimination of uncured resin from 3D-printed dental aligners

The study investigated the effects of temperature and centrifugation time on the efficacy of removing uncured resin from 3D-printed clear aligners. Using a photo-polymerizable polyurethane resin (Tera Harz TC-85, Graphy Inc., Seoul, Korea), aligners were printed and subjected to cleaning processes using isopropyl alcohol (IPA) or centrifugation (g-force 27.95g) at room temperature (RT, 23 °C) and high temperature (HT, 55 °C) for 2, 4, and 6 min. The control group received no treatment (NT). Cleaning efficiency was assessed through rheological analysis, weight measurement, transparency evaluation, SEM imaging, 3D geometry evaluation, stress relaxation, and cell viability tests. Results showed increased temperature and longer centrifugation times significantly reduced aligner viscosity, weight (P < 0.05), and transmittance. IPA-cleaned aligners exhibited significantly lower transparency and rougher surfaces in SEM images. All groups met ISO biocompatibility standards in cytotoxicity tests. The NT group had higher root mean square (RMS) values, indicating greater deviation from the original design. Stress relaxation tests revealed over 95% recovery in all groups after 60 min. The findings suggest that a 2-min HT centrifugation process effectively removes uncured resin without significantly impacting the aligners’ physical and optical properties, making it a clinically viable option.

On the Schapery nonlinear viscoelastic model: A review

Among the various models presented so far to describe the nonlinear viscoelastic behavior of materials, the model proposed by Schapery in the 1960s stands out for its wide acceptance and significant contributions to the field. Through a detailed inspection of the main articles presented by Schapery and other related studies, the present paper reviews the Schapery nonlinear viscoelastic model. One-dimensional and three-dimensional formulations of the model for both isotropic and anisotropic materials, along with model accuracy, are discussed. The present article highlights the challenges of material characterization methods, such as creep-recovery testing, stress relaxation testing, and dynamic mechanical analysis. The article also addresses the issues of characterizing the material's behavior under non-recoverable strain and varying temperatures. Finally, the future research needs in this field are outlined, emphasizing the potential for further advancements in our understanding of nonlinear viscoelasticity.