Wave propagation and the frequency domain Green's functions in viscoelastic Biot/squirt (BISQ) media

In this paper, we examine the characteristics of elastic wave propagation in viscoelastic porous media, which contain simultaneously both the Biot-flow and the squirt-flow mechanisms (BISQ). The frequency-domain Green's functions for viscoelastic BISQ media are then derived based on the classic potential function methods. Our numerical results show that S-waves are only affected by viscoelasticity, but not by squirt-flows. However, the phase velocity and attenuation of fast P-waves are seriously influenced by both viscoelasticity and squirt-flows; and there exist two peaks in the attenuation-frequency variations of fast P-waves. In the low-frequency range, the squirt-flow characteristic length, not viscoelasticity, affects the phase velocity of slow P-waves, whereas it is opposite in the high-frequency range. As to the contribution of potential functions of two types of compressional waves to the Green's function, the squirt-flow length has a small effect, and the effects of viscoelastic parameter are mainly in the higher frequency range. Crown Copyright (C) 2006 Published by Elsevier Ltd. All rights reserved.

Creep characterization of a fiber reinforced plastic material

Creep behavior of [±45°]_s composite material is characterized by using uniaxial creep and recovery tests. The well-known Schapery nonlinear viscoelastic consti tutive relation was modified to make it suitable for characterizing the creep behavior of this material. Then, using this modified Schapery constitutive equation, by which the vis coplastic and creep damage can be taken into consideration, the creep behavior of [±45°]. glass fiber reinforced epoxy laminate was studied. The constitutive parameters of the material were determined experimentally, and the procedure and method of determination of the material parameters are proved to be valid.

Diffusion in glassy polymers

Fluid diffusion in glassy polymers proceeds in ways that are not explained by the standard diffusion model. Although the reasons for the anomalous effects are not known, much of the observed behavior is attributed to the long times that polymers below their glass transition temperature take to adjust to changes in their condition. The slow internal relaxations of the polymer chains ensure that the material properties are history-dependent, and also allow both local inhomogeneities and differential swelling to occur. Two models are developed in this thesis with the intent of accounting for these effects in the diffusion process.

In Part I, a model is developed to account for both the history dependence of the glassy polymer, and the dual sorption which occurs when gas molecules are immobilized by the local heterogeneities. A preliminary study of a special case of this model is conducted, showing the existence of travelling wave solutions and using perturbation techniques to investigate the effect of generalized diffusion mechanisms on their form. An integral averaging method is used to estimate the penetrant front position.

In Part II, a model is developed for particle diffusion along with displacements in isotropic viscoelastic materials. The nonlinear dependence of the materials on the fluid concentration is taken into account, while pure displacements are assumed to remain in the range of linear viscoelasticity. A fairly general model is obtained for three-dimensional irrotational movements, with the development of the model being based on the assumptions of irreversible thermodynamics. With the help of some dimensional analysis, this model is simplified to a version which is proposed to be studied for Case II behavior.

The short-timescale behavior of glacial ice

Glaciers are often assumed to deform only at slow (i.e., glacial) rates. However, with the advent of high rate geodetic observations of ice motion, many of the intricacies of glacial deformation on hourly and daily timescales have been observed and quantified. This thesis explores two such short timescale processes: the tidal perturbation of ice stream motion and the catastrophic drainage of supraglacial meltwater lakes. Our investigation into the transmission length-scale of a tidal load represents the first study to explore the daily tidal influence on ice stream motion using three-dimensional models. Our results demonstrate both that the implicit assumptions made in the standard two-dimensional flow-line models are inherently incorrect for many ice streams, and that the anomalously large spatial extent of the tidal influence seen on the motion of some glaciers cannot be explained, as previously thought, through the elastic or viscoelastic transmission of tidal loads through the bulk of the ice stream. We then discuss how the phase delay between a tidal forcing and the ice stream’s displacement response can be used to constrain in situ viscoelastic properties of glacial ice. Lastly, for the problem of supraglacial lake drainage, we present a methodology for implementing linear viscoelasticity into an existing model for lake drainage. Our work finds that viscoelasticity is a second-order effect when trying to model the deformation of ice in response to a meltwater lake draining to a glacier’s bed. The research in this thesis demonstrates that the first-order understanding of the short-timescale behavior of naturally occurring ice is incomplete, and works towards improving our fundamental understanding of ice behavior over the range of hours to days.

Extracting material response from simple mechanical tests on hardening-softening-hardening viscoplastic solids

Compliant foams are usually characterized by a wide range of desirable mechanical properties. These properties include viscoelasticity at different temperatures, energy absorption, recoverability under cyclic loading, impact resistance, and thermal, electrical, acoustic and radiation-resistance. Some foams contain nano-sized features and are used in small-scale devices. This implies that the characteristic dimensions of foams span multiple length scales, rendering modeling their mechanical properties difficult. Continuum mechanics-based models capture some salient experimental features like the linear elastic regime, followed by non-linear plateau stress regime. However, they lack mesostructural physical details. This makes them incapable of accurately predicting local peaks in stress and strain distributions, which significantly affect the deformation paths. Atomistic methods are capable of capturing the physical origins of deformation at smaller scales, but suffer from impractical computational intensity. Capturing deformation at the so-called meso-scale, which is capable of describing the phenomenon at a continuum level, but with some physical insights, requires developing new theoretical approaches.

A fundamental question that motivates the modeling of foams is ‘how to extract the intrinsic material response from simple mechanical test data, such as stress vs. strain response?’ A 3D model was developed to simulate the mechanical response of foam-type materials. The novelty of this model includes unique features such as the hardening-softening-hardening material response, strain rate-dependence, and plastically compressible solids with plastic non-normality. Suggestive links from atomistic simulations of foams were borrowed to formulate a physically informed hardening material input function. Motivated by a model that qualitatively captured the response of foam-type vertically aligned carbon nanotube (VACNT) pillars under uniaxial compression [2011,“Analysis of Uniaxial Compression of Vertically Aligned Carbon Nanotubes,” J. Mech.Phys. Solids, 59, pp. 2227–2237, Erratum 60, 1753–1756 (2012)], the property space exploration was advanced to three types of simple mechanical tests: 1) uniaxial compression, 2) uniaxial tension, and 3) nanoindentation with a conical and a flat-punch tip. The simulations attempt to explain some of the salient features in experimental data, like
1) The initial linear elastic response.
2) One or more nonlinear instabilities, yielding, and hardening.

The model-inherent relationships between the material properties and the overall stress-strain behavior were validated against the available experimental data. The material properties include the gradient in stiffness along the height, plastic and elastic compressibility, and hardening. Each of these tests was evaluated in terms of their efficiency in extracting material properties. The uniaxial simulation results proved to be a combination of structural and material influences. Out of all deformation paths, flat-punch indentation proved to be superior since it is the most sensitive in capturing the material properties.

Simulation of shock-induced bubble collapse with application to vascular injury in shockwave lithotripsy

Shockwave lithotripsy is a noninvasive medical procedure wherein shockwaves are repeatedly focused at the location of kidney stones in order to pulverize them. Stone comminution is thought to be the product of two mechanisms: the propagation of stress waves within the stone and cavitation erosion. However, the latter mechanism has also been implicated in vascular injury. In the present work, shock-induced bubble collapse is studied in order to understand the role that it might play in inducing vascular injury. A high-order accurate, shock- and interface-capturing numerical scheme is developed to simulate the three-dimensional collapse of the bubble in both the free-field and inside a vessel phantom. The primary contributions of the numerical study are the characterization of the shock-bubble and shock-bubble-vessel interactions across a large parameter space that includes clinical shockwave lithotripsy pressure amplitudes, problem geometry and tissue viscoelasticity, and the subsequent correlation of these interactions to vascular injury. Specifically, measurements of the vessel wall pressures and displacements, as well as the finite strains in the fluid surrounding the bubble, are utilized with available experiments in tissue to evaluate damage potential. Estimates are made of the smallest injurious bubbles in the microvasculature during both the collapse and jetting phases of the bubble's life cycle. The present results suggest that bubbles larger than 1 μm in diameter could rupture blood vessels under clinical SWL conditions.

Collective behavior of asperities as a model for friction and adhesion

Understanding friction and adhesion in static and sliding contact of surfaces is important in numerous physical phenomena and technological applications. Most surfaces are rough at the microscale, and thus the real area of contact is only a fraction of the nominal area. The macroscopic frictional and adhesive response is determined by the collective behavior of the population of evolving and interacting microscopic contacts. This collective behavior can be very different from the behavior of individual contacts. It is thus important to understand how the macroscopic response emerges from the microscopic one. In this thesis, we develop a theoretical and computational framework to study the collective behavior. Our philosophy is to assume a simple behavior of a single asperity and study the collective response of an ensemble. Our work bridges the existing well-developed studies of single asperities with phenomenological laws that describe macroscopic rate-and-state behavior of frictional interfaces. We find that many aspects of the macroscopic behavior are robust with respect to the microscopic response. This explains why qualitatively similar frictional features are seen for a diverse range of materials. We first show that the collective response of an ensemble of one-dimensional independent viscoelastic elements interacting through a mean field reproduces many qualitative features of static and sliding friction evolution. The resulting macroscopic behavior is different from the microscopic one: for example, even if each contact is velocity-strengthening, the macroscopic behavior can be velocity-weakening. The framework is then extended to incorporate three-dimensional rough surfaces, long- range elastic interactions between contacts, and time-dependent material behaviors such as viscoelasticity and viscoplasticity. Interestingly, the mean field behavior dominates and the elastic interactions, though important from a quantitative perspective, do not change the qualitative macroscopic response. Finally, we examine the effect of adhesion on the frictional response as well as develop a force threshold model for adhesion and mode I interfacial cracks.

Exploring the Kinetics of Domain Switching in Ferroelectrics for Structural Applications

The complex domain structure in ferroelectrics gives rise to electromechanical coupling, and its evolution (via domain switching) results in a time-dependent (i.e. viscoelastic) response. Although ferroelectrics are used in many technological applications, most do not attempt to exploit the viscoelastic response of ferroelectrics, mainly due to a lack of understanding and accurate models for their description and prediction. Thus, the aim of this thesis research is to gain better understanding of the influence of domain evolution in ferroelectrics on their dynamic mechanical response. There have been few studies on the viscoelastic properties of ferroelectrics, mainly due to a lack of experimental methods. Therefore, an apparatus and method called Broadband Electromechanical Spectroscopy (BES) was designed and built. BES allows for the simultaneous application of dynamic mechanical and electrical loading in a vacuum environment. Using BES, the dynamic stiffness and loss tangent in bending and torsion of a particular ferroelectric, viz. lead zirconate titanate (PZT), was characterized for different combinations of electrical and mechanical loading frequencies throughout the entire electric displacement hysteresis. Experimental results showed significant increases in loss tangent (by nearly an order of magnitude) and compliance during domain switching, which shows promise as a new approach to structural damping. A continuum model of the viscoelasticity of ferroelectrics was developed, which incorporates microstructural evolution via internal variables and associated kinetic relations. For the first time, through a new linearization process, the incremental dynamic stiffness and loss tangent of materials were computed throughout the entire electric displacement hysteresis for different combinations of mechanical and electrical loading frequencies. The model accurately captured experimental results. Using the understanding gained from the characterization and modeling of PZT, two applications of domain switching kinetics were explored by using Micro Fiber Composites (MFCs). Proofs of concept of set-and-hold actuation and structural damping using MFCs were demonstrated.

Kautxu birziklatuzko materialek inpaktura duten erantzunaren simulazioa.

[EU]Ondorengo lanean ingeniaritzaren eremuan aurkitzen den eginkizunik garrantzitsuenetako bat egingo da: material berri bat karakterizatu. Hau da, errealitatean ikusten duguna era matematiko batean idatzi, edo beste era batean esanda, modelo matematiko bat lortu. Jorratuko den materiala kautxu birziklatua izango da, pneumatiko zaharretatik abiatuz egiten dena. Material honek, besteak beste, hurrengo abantailak ditu: produzitzeko erraza (beraz, prezio ekonomikoak, lehengaiak hondakin solidoak direla ikusita) eta inpakturako erantzun ona. Kautxu hau, beste erabileren artean, segurtasun-errailetan erabiltzen da, moto-gidariak ebaketetatik babesteko. Modelo matematikoa edo konstitutiboa lortzeko, lege hiperelastiko eta biskoelastikoetatik abiatuz, portaera bisko-hiperelastikorako lege bat proposatuko da. Lege hau deformazio handietarako eta deformazio abiadura handietarako aplikatu egingo da. Honen oinarria aurreko urteetan eskola honetako beste ikasleek eginiko gradu eta karrera amaierako lanak izango dira. Behin hau izanda, elementu finituetako software baten bidez simulatuko da, eta emaitza numerikoak esperimentalekin alderatuko dira. Honen bidez, modeloa hobetu ahal izango dugu, zehaztasun gehiago lortzeko asmoz.

A model for nonlinear viscoelastic mechanical responses of collagenous soft tissues

Experimental observations of the time-dependent mechanical responses of collagenous tissues have demonstrated behavior that deviates from standard treatments of linear or quasi-linear viscoelasticity. In particular, time-dependent deformation can be strongly coupled to strain level, and strain-rate independence can be observed under monotonic loading, even for a tissue with dramatic stress relaxation. It was postulated that this nonlinearity is fundamentally associated with gradual recruitment of individual collagen fibrils during applied mechanical loading. Based on previously observed experimental results for the time-dependent response of collagenous soft tissues, a model is developed to describe the mechanical behavior of these tissues under uniaxial loading. Tissue stresses, under applied strain-controlled loading, are assumed to be a sum of elastic and viscoelastic stress contributions. The relative contributions of elastic and viscoelastic stresses is assumed to vary with strain level, leading to strain- and time-dependent mechanical behavior. The model formulation is examined under conditions of monotonic loading at varying constant strain rates and stress-relaxation at different applied strain levels. The model is compared with experimental data for a membranous biological soft tissue, the amniotic sac, and is found to agree well with experimental results. The limiting behavior of the novel model, at large strains relative to the collagen recruitment, is consistent with the quasi-linear viscoelastic approach. © 2006 Materials Research Society.

Poroelastic indentation analysis for hydrated biological tissues

Indentation techniques are employed for the measurement of mechanical properties of a wide range of materials. In particular, techniques focused at small length-scales, such as nanoindentation and AFM indentation, allow for local characterization of material properties in heterogeneous materials including natural tissues and biomimetic materials. Typical elastic analysis for spherical indentation is applicable in the absence of time-dependent deformation, but is inappropriate for materials with time-dependent responses. Recent analyses for the viscoelastic indentation problem, based on elastic-viscoelastic correspondence, have begun to address the issue of time-dependent deformation during an indentation test. The viscoelastic analysis has been shown to fit experimental indentation data well, and has been demonstrated as useful for characterization of viscoelasticity in polymeric materials and in hydrated mineralized tissues. However, a viscoelastic analysis is not necessarily sufficient for multi-phase materials with fluid flow. In the current work, a poroelastic analysis-based on fluid motion through a porous elastic network-is used to examine spherical indentation creep responses of hydrated biological materials. Both analytical and finite element approaches are considered for the poroelastic Hertzian indentation problem. Modeling results are compared with experimental data from nanoindentation of hydrated bone immersed in water and polar solvents (ethanol, methanol, acetone). Baseline (water-immersed) bone responses are characterized using the poroelastic model and numerical results are compared with altered hydration states due to polar solvents. © 2007 Materials Research Society.

Short- and long-term changes in joint co-contraction associated with motor learning as revealed from surface EMG.

In the field of motor control, two hypotheses have been controversial: whether the brain acquires internal models that generate accurate motor commands, or whether the brain avoids this by using the viscoelasticity of musculoskeletal system. Recent observations on relatively low stiffness during trained movements support the existence of internal models. However, no study has revealed the decrease in viscoelasticity associated with learning that would imply improvement of internal models as well as synergy between the two hypothetical mechanisms. Previously observed decreases in electromyogram (EMG) might have other explanations, such as trajectory modifications that reduce joint torques. To circumvent such complications, we required strict trajectory control and examined only successful trials having identical trajectory and torque profiles. Subjects were asked to perform a hand movement in unison with a target moving along a specified and unusual trajectory, with shoulder and elbow in the horizontal plane at the shoulder level. To evaluate joint viscoelasticity during the learning of this movement, we proposed an index of muscle co-contraction around the joint (IMCJ). The IMCJ was defined as the summation of the absolute values of antagonistic muscle torques around the joint and computed from the linear relation between surface EMG and joint torque. The IMCJ during isometric contraction, as well as during movements, was confirmed to correlate well with joint stiffness estimated using the conventional method, i.e., applying mechanical perturbations. Accordingly, the IMCJ during the learning of the movement was computed for each joint of each trial using estimated EMG-torque relationship. At the same time, the performance error for each trial was specified as the root mean square of the distance between the target and hand at each time step over the entire trajectory. The time-series data of IMCJ and performance error were decomposed into long-term components that showed decreases in IMCJ in accordance with learning with little change in the trajectory and short-term interactions between the IMCJ and performance error. A cross-correlation analysis and impulse responses both suggested that higher IMCJs follow poor performances, and lower IMCJs follow good performances within a few successive trials. Our results support the hypothesis that viscoelasticity contributes more when internal models are inaccurate, while internal models contribute more after the completion of learning. It is demonstrated that the CNS regulates viscoelasticity on a short- and long-term basis depending on performance error and finally acquires smooth and accurate movements while maintaining stability during the entire learning process.

Prediction of tread block forces for a free-rolling tyre in contact with a rough road

This paper describes a new approach to model the forces on a tread block for a free-rolling tyre in contact with a rough road. A theoretical analysis based on realistic tread mechanical properties and road roughness is presented, indicating partial contact between a tread block and a rough road. Hence an asperity-scale indentation model is developed using a semi-empirical formulation, taking into account both the rubber viscoelasticity and the tread block geometry. The model aims to capture the essential details of the contact at the simplest level, to make it suitable as part of a time-domain dynamic analysis of the coupled tyre-road system. The indentation model is found to have a good correlation with the finite element (FE) predictions and is validated against experimental results using a rolling contact rig. When coupled to a deformed tyre belt profile, the indentation model predicts normal and tangential force histories inside the tyre contact patch that show good agreement with FE predictions. © 2012 Elsevier B.V..

How PEDOT:PSS solutions produce satellite-free inkjets

Satellite droplets are unwanted in inkjet printing and various approaches have been suggested for their reduction. Low jetting speeds limit applications of the process. Added surfactants for wetting and conductivity enhancement may help but dynamic surface tension effects may counteract improvements. A higher fluid viscosity delays ligament break-up, but also leads to slower jets, while viscoelasticity reduces satellite formation only in certain cases. We show here that aqueous solutions of PEDOT:PSS (1:2.5 by weight) are strongly shear-thinning. They exhibit low viscosity within the printing nozzle over a wide range of jet speeds, yet rapidly (<100 μs) recover a higher viscosity at the low shear rates applicable once the jet has formed, which give the benefit of delayed satellite formation. The delay over a 0.8 mm stand-off distance can be sufficient to completely suppress satellites, which is significant for many printing applications. © 2012 Elsevier B.V. All rights reserved.

Separating poroviscoelastic deformation mechanisms in hydrogels

Hydrogels have applications in drug delivery, mechanical actuation, and regenerative medicine. When hydrogels are deformed, load-relaxation arising from fluid flow - poroelasticity - and from rearrangement of the polymer network - viscoelasticity - is observed. The physical mechanisms are different in that poroelastic relaxation varies with experimental length-scale while viscoelastic does not. Here, we show that poroviscoelastic load-relaxation is the product of the two individual responses. The difference in length-scale dependence of the two mechanisms can be exploited to uniquely determine poroviscoelastic properties from simultaneous analysis of multi-scale indentation experiments, providing insight into hydrogel physical behavior. © 2013 American Institute of Physics.