Thermodynamics of viscoelastic rate-type fluids with stress diffusion

NASA Astrophysics Data System (ADS)

Málek, Josef; Průša, Vít; Skřivan, Tomáš; Süli, Endre

2018-02-01

We propose thermodynamically consistent models for viscoelastic fluids with a stress diffusion term. In particular, we derive variants of compressible/incompressible Maxwell/Oldroyd-B models with a stress diffusion term in the evolution equation for the extra stress tensor. It is shown that the stress diffusion term can be interpreted either as a consequence of a nonlocal energy storage mechanism or as a consequence of a nonlocal entropy production mechanism, while different interpretations of the stress diffusion mechanism lead to different evolution equations for the temperature. The benefits of the knowledge of the thermodynamical background of the derived models are documented in the study of nonlinear stability of equilibrium rest states. The derived models open up the possibility to study fully coupled thermomechanical problems involving viscoelastic rate-type fluids with stress diffusion.

Modelling Viscoelastic Behaviour of Polymer by A Mixed Velocity, Displacement Formulation - Numerical and Experimental Results

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pham, VT.; Silva, L.; Digonnet, H.

2011-05-04

The objective of this work is to model the viscoelastic behaviour of polymer from the solid state to the liquid state. With this objective, we perform experimental tensile tests and compare with simulation results. The chosen polymer is a PMMA whose behaviour depends on its temperature. The computation simulation is based on Navier-Stokes equations where we propose a mixed finite element method with an interpolation P1+/P1 using displacement (or velocity) and pressure as principal variables. The implemented technique uses a mesh composed of triangles (2D) or tetrahedra (3D). The goal of this approach is to model the viscoelastic behaviour ofmore » polymers through a fluid-structure coupling technique with a multiphase approach.« less

Fluid-structure interaction with pipe-wall viscoelasticity during water hammer

NASA Astrophysics Data System (ADS)

Keramat, A.; Tijsseling, A. S.; Hou, Q.; Ahmadi, A.

2012-01-01

Fluid-structure interaction (FSI) due to water hammer in a pipeline which has viscoelastic wall behaviour is studied. Appropriate governing equations are derived and numerically solved. In the numerical implementation of the hydraulic and structural equations, viscoelasticity is incorporated using the Kelvin-Voigt mechanical model. The equations are solved by two different approaches, namely the Method of Characteristics-Finite Element Method (MOC-FEM) and full MOC. In both approaches two important effects of FSI in fluid-filled pipes, namely Poisson and junction coupling, are taken into account. The study proposes a more comprehensive model for studying fluid transients in pipelines as compared to previous works, which take into account either FSI or viscoelasticity. To verify the proposed mathematical model and its numerical solutions, the following problems are investigated: axial vibration of a viscoelastic bar subjected to a step uniaxial loading, FSI in an elastic pipe, and hydraulic transients in a pressurised polyethylene pipe without FSI. The results of each case are checked with available exact and experimental results. Then, to study the simultaneous effects of FSI and viscoelasticity, which is the new element of the present research, one problem is solved by the two different numerical approaches. Both numerical methods give the same results, thus confirming the correctness of the solutions.

Finite element method for viscoelastic medium with damage and the application to structural analysis of solid rocket motor grain

NASA Astrophysics Data System (ADS)

Deng, Bin; Shen, ZhiBin; Duan, JingBo; Tang, GuoJin

2014-05-01

This paper studies the damage-viscoelastic behavior of composite solid propellants of solid rocket motors (SRM). Based on viscoelastic theories and strain equivalent hypothesis in damage mechanics, a three-dimensional (3-D) nonlinear viscoelastic constitutive model incorporating with damage is developed. The resulting viscoelastic constitutive equations are numerically discretized by integration algorithm, and a stress-updating method is presented by solving nonlinear equations according to the Newton-Raphson method. A material subroutine of stress-updating is made up and embedded into commercial code of Abaqus. The material subroutine is validated through typical examples. Our results indicate that the finite element results are in good agreement with the analytical ones and have high accuracy, and the suggested method and designed subroutine are efficient and can be further applied to damage-coupling structural analysis of practical SRM grain.

Implications of the earthquake cycle for inferring fault locking on the Cascadia megathrust

USGS Publications Warehouse

Pollitz, Fred; Evans, Eileen

2017-01-01

GPS velocity fields in the Western US have been interpreted with various physical models of the lithosphere-asthenosphere system: (1) time-independent block models; (2) time-dependent viscoelastic-cycle models, where deformation is driven by viscoelastic relaxation of the lower crust and upper mantle from past faulting events; (3) viscoelastic block models, a time-dependent variation of the block model. All three models are generally driven by a combination of loading on locked faults and (aseismic) fault creep. Here we construct viscoelastic block models and viscoelastic-cycle models for the Western US, focusing on the Pacific Northwest and the earthquake cycle on the Cascadia megathrust. In the viscoelastic block model, the western US is divided into blocks selected from an initial set of 137 microplates using the method of Total Variation Regularization, allowing potential trade-offs between faulting and megathrust coupling to be determined algorithmically from GPS observations. Fault geometry, slip rate, and locking rates (i.e. the locking fraction times the long term slip rate) are estimated simultaneously within the TVR block model. For a range of mantle asthenosphere viscosity (4.4 × 1018 to 3.6 × 1020 Pa s) we find that fault locking on the megathrust is concentrated in the uppermost 20 km in depth, and a locking rate contour line of 30 mm yr−1 extends deepest beneath the Olympic Peninsula, characteristics similar to previous time-independent block model results. These results are corroborated by viscoelastic-cycle modelling. The average locking rate required to fit the GPS velocity field depends on mantle viscosity, being higher the lower the viscosity. Moreover, for viscosity ≲ 1020 Pa s, the amount of inferred locking is higher than that obtained using a time-independent block model. This suggests that time-dependent models for a range of admissible viscosity structures could refine our knowledge of the locking distribution and its epistemic uncertainty.

Pressure induced ageing of polymers

NASA Technical Reports Server (NTRS)

Emri, I.; Knauss, W. G.

1988-01-01

The nonlinearly viscoelastic response of an amorphous homopolymer is considered under aspects of time dependent free volume behavior. In contrast to linearly viscoelastic solids, this model couples shear and volume deformation through a shift function which influences the rate of molecular relaxation or creep. Sample computations produce all those qualitative features one observes normally in uniaxial tension including the rate dependent formation of a yield point as a consequence of the history of an imposed pressure.

Coupled dynamics of a viscoelastically supported infinite string and a number of discrete mechanical systems moving with uniform speed

NASA Astrophysics Data System (ADS)

Roy, Soumyajit; Chakraborty, G.; DasGupta, Anirvan

2018-02-01

The mutual interaction between a number of multi degrees of freedom mechanical systems moving with uniform speed along an infinite taut string supported by a viscoelastic layer has been studied using the substructure synthesis method when base excitations of a common frequency are given to the mechanical systems. The mobility or impedance matrices of the string have been calculated analytically by Fourier transform method as well as wave propagation technique. The above matrices are used to calculate the response of the discrete mechanical systems. Special attention is paid to the contact forces between the discrete and the continuous systems which are estimated by numerical simulation. The effects of phase difference, the distance between the systems and different base excitation amplitudes on the collective behaviour of the mechanical systems are also studied. The present study has relevance to the coupled dynamic problem of more than one railway pantographs and an overhead catenary system where the pantographs are modelled as discrete systems and the catenary is modelled as a taut string supported by continuous viscoelastic layer.

Study of dynamic fluid-structure coupling with application to human phonation

NASA Astrophysics Data System (ADS)

Saurabh, Shakti; Faber, Justin; Bodony, Daniel

2013-11-01

Two-dimensional direct numerical simulations of a compressible, viscous fluid interacting with a non-linear, viscoelastic solid are used to study the generation of the human voice. The vocal fold (VF) tissues are modeled using a finite-strain fractional derivative constitutive model implemented in a quadratic finite element code and coupled to a high-order compressible Navier-Stokes solver through a boundary-fitted fluid-solid interface. The viscoelastic solver is validated through in-house experiments using Agarose Gel, a human tissue simulant, undergoing static and harmonic deformation measured with load cell and optical diagnostics. The phonation simulations highlight the role tissue nonlinearity and viscosity play in the glottal jet dynamics and in the radiated sound. Supported by the National Science Foundation (CAREER award number 1150439).

Strong interseismic coupling, fault afterslip, and viscoelastic flow before and after the Oct. 9, 1995 Colima-Jalisco earthquake: continuous GPS measurements from Colima, Mexico

USGS Publications Warehouse

Azua, B.M.; DeMets, C.; Masterlark, Timothy

2002-01-01

Continuous GPS measurements from Colima, Mexico during 4/93-6/01, bracketing the Oct. 9, 1995 M = 8.0 Colima-Jalisco earthquake, provide new constraints on Rivera plate subduction mechanics. Modeling of margin-normal strain accumulation before the earthquake suggests the Rivera-North America subduction interface was fully locked. Transient postseismic motion from 10/ 95-6/97 is well fit by a model that includes logarithmically-decaying fault afterslip, elastic strain from shallow fault relocking, and possibly a minor viscoelastic response, but is fit poorly by models that assume a dominant Maxwell viscoelastic response of the lower crust and upper mantle, independent of the assumed viscosities. Landward, margin-normal motion since mid-1997 is parallel to but ??? 75% slower than the pre-seismic velocity. Afterslip alone fails to account for this slowdown. The viscoelastic response predicted by a FEM correctly resolves the remaining velocity difference within the uncertainties. Both processes thus offset elastic strain accumulating from the relocked subduction interface.

Jeans Instability of the Self-Gravitating Viscoelastic Ferromagnetic Cylinder with Axial Nonuniform Rotation and Magnetic Field

NASA Astrophysics Data System (ADS)

Dhiman, Joginder Singh; Sharma, Rajni

2017-12-01

The effects of nonuniform rotation and magnetic field on the instability of a self gravitating infinitely extending axisymmetric cylinder of viscoelastic ferromagnetic medium have been studied using the Generalised Hydrodynamic (GH) model. The non-uniform magnetic field and rotation are acting along the axial direction of the cylinder and the propagation of the wave is considered along the radial direction, while the ferrofluid magnetization is taken collinear with the magnetic field. A general dispersion relation representing magnetization, magnetic permeability and viscoelastic relaxation time parameters is obtained using the normal mode analysis method in the linearized perturbation equation system. Jeans criteria which represent the onset of instability of self gravitating medium are obtained under the limits; when the medium behaves like a viscous liquid (strongly coupled limit) and a Newtonian liquid (weakly coupled limit). The effects of various parameters on the Jeans instability criteria and on the growth rate of self gravitating viscoelastic ferromagnetic medium have been discussed. It is found that the magnetic polarizability due to ferromagnetization of medium marginalizes the effect of non-uniform magnetic field on the Jeans instability, whereas the viscoelasticity of the medium has the usual stabilizing effect on the instability of the system. Further, it is found that the cylindrical geometry is more stable than the Cartesian one. The variation of growth rate against the wave number and radial distance has been depicted graphically.

Nonlinear Viscoelastic Mechanism for Aftershock Triggering and Decay

NASA Astrophysics Data System (ADS)

Shcherbakov, R.; Zhang, X.

2016-12-01

Aftershocks are ubiquitous in nature. They are the manifestation of relaxation phenomena observed in various physical systems. In one prominent example, they typically occur after large earthquakes. They also occur in other natural or experimental systems, for example, in solar flares, in fracture experiments on porous materials and acoustic emissions, after stock market crashes, in the volatility of stock prices returns, in internet traffic variability and e-mail spamming, to mention a few. The observed aftershock sequences usually obey several well defined non-trivial empirical laws in magnitude, temporal, and spatial domains. In many cases their characteristics follow scale-invariant distributions. The occurrence of aftershocks displays a prominent temporal behavior due to time-dependent mechanisms of stress and/or energy transfer. In this work, we consider a slider-block model to mimic the behavior of a seismogenic fault. In the model, we introduce a nonlinear viscoelastic coupling mechanism to capture the essential characteristics of crustal rheology and stress interaction between the blocks and the medium. For this purpose we employ nonlinear Kelvin-Voigt elements consisting of an elastic spring and a dashpot assembled in parallel to introduce viscoelastic coupling between the blocks and the driving plate. By mapping the model into a cellular automaton we derive the functional form of the stress transfer mechanism in the model. We show that the nonlinear viscoelasticity plays a critical role in triggering of aftershocks. It explains the functional form of the Omori-Utsu law and gives physical interpretation of its parameters. The proposed model also suggests that the power-law rheology of the fault gauge and underlying lower crust and upper mantle control the decay rate of aftershocks. To verify this, we analyze several prominent aftershock sequences to estimate their decay rates and correlate with the rheological properties of the underlying lower crust and mantle.

Process Modelling of Curing Process-Induced Internal Stress and Deformation of Composite Laminate Structure with Elastic and Viscoelastic Models

NASA Astrophysics Data System (ADS)

Li, Dongna; Li, Xudong; Dai, Jianfeng

2018-06-01

In this paper, two kinds of transient models, the viscoelastic model and the linear elastic model, are established to analyze the curing deformation of the thermosetting resin composites, and are calculated by COMSOL Multiphysics software. The two models consider the complicated coupling between physical and chemical changes during curing process of the composites and the time-variant characteristic of material performance parameters. Subsequently, the two proposed models are implemented respectively in a three-dimensional composite laminate structure, and a simple and convenient method of local coordinate system is used to calculate the development of residual stresses, curing shrinkage and curing deformation for the composite laminate. Researches show that the temperature, degree of curing (DOC) and residual stresses during curing process are consistent with the study in literature, so the curing shrinkage and curing deformation obtained on these basis have a certain referential value. Compared the differences between the two numerical results, it indicates that the residual stress and deformation calculated by the viscoelastic model are more close to the reference value than the linear elastic model.

Sea-level and solid-Earth deformation feedbacks in ice sheet modelling

NASA Astrophysics Data System (ADS)

Konrad, Hannes; Sasgen, Ingo; Klemann, Volker; Thoma, Malte; Grosfeld, Klaus; Martinec, Zdeněk

2014-05-01

The interactions of ice sheets with the sea level and the solid Earth are important factors for the stability of the ice shelves and the tributary inland ice (e.g. Thomas and Bentley, 1978; Gomez et al, 2012). First, changes in ice extent and ice thickness induce viscoelastic deformation of the Earth surface and Earth's gravity field. In turn, global and local changes in sea level and bathymetry affect the grounding line and, subsequently, alter the ice dynamic behaviour. Here, we investigate these feedbacks for a synthetic ice sheet configuration as well as for the Antarctic ice sheet using a three-dimensional thermomechanical ice sheet and shelf model, coupled to a viscoelastic solid-Earth and gravitationally self-consistent sea-level model. The respective ice sheet undergoes a forcing from rising sea level, warming ocean, and/or changing surface mass balance. The coupling is realized by exchanging ice thickness, Earth surface deformation and sea level periodically. We apply several sets of viscoelastic Earth parameters to our coupled model, e.g. simulating a low-viscous upper mantle present at the Antarctic Peninsula (Ivins et al., 2011). Special focus of our study lies on the evolution of Earth surface deformation and local sea level changes, as well as on the accompanying grounding line evolution. N. Gomez, D. Pollard, J. X. Mitrovica, P. Huybers, and P. U. Clark 2012. Evolution of a coupled marine ice sheet-sea level model, J. Geophys. Res., 117, F01013, doi:10.1029/2011JF002128. E. R. Ivins, M. M. Watkins, D.-N. Yuan, R. Dietrich, G. Casassa, and A. Rülke 2011. On-land ice loss and glacial isostatic adjustment at the Drake Passage: 2003-2009, J. Geophys. Res. 116, B02403, doi: 10.1029/2010JB007607 R. H. Thomas and C. R. Bentley 1978. A model for Holocene retreat of the West Antarctic Ice Sheet, Quaternary Research, 10 (2), pages 150-170, doi: 10.1016/0033-5894(78)90098-4.

Rate dependent constitutive behavior of dielectric elastomers and applications in legged robotics

NASA Astrophysics Data System (ADS)

Oates, William; Miles, Paul; Gao, Wei; Clark, Jonathan; Mashayekhi, Somayeh; Hussaini, M. Yousuff

2017-04-01

Dielectric elastomers exhibit novel electromechanical coupling that has been exploited in many adaptive structure applications. Whereas the quasi-static, one-dimensional constitutive behavior can often be accurately quantified by hyperelastic functions and linear dielectric relations, accurate predictions of electromechanical, rate-dependent deformation during multiaxial loading is non-trivial. In this paper, an overview of multiaxial electromechanical membrane finite element modeling is formulated. Viscoelastic constitutive relations are extended to include fractional order. It is shown that fractional order viscoelastic constitutive relations are superior to conventional integer order models. This knowledge is critical for transition to control of legged robotic structures that exhibit advanced mobility.

Postseismic deformation following the 2015 Gorkha earthquake and implications for rheology

NASA Astrophysics Data System (ADS)

Rollins, C.; Gualandi, A.; Avouac, J. P.; Liu, J.; Zhang, Z.

2017-12-01

The 2015 Mw 7.9 Gorkha earthquake ruptured the lower, northern edge of the interseismically locked section of the Main Himalayan Thrust (MHT). Independent Component Analysis of location timeseries at GPS stations in Nepal and Tibet reveals significant transient postseismic motion following the mainshock. In order to probe the frictional properties of the MHT and the viscoelastic properties of the crust and upper mantle, we compare the extracted postseismic motions to those predicted by forward models of afterslip and viscoelastic relaxation. Postseismic displacements are minimal south of the coseismic rupture, suggesting that minimal afterslip occurred there and that the upper MHT remains mostly locked. North of the rupture, postseismic displacements feature south-southwest horizontal motion and uplift, each on the order of a few cm in the first postseismic year. A model of stress-driven afterslip extending 100 km north of the coseismic rupture reproduces the horizontal postseismic timeseries and the general pattern of uplift and subsidence; however, this model significantly overpredicts the uplift at stations overlying the rupture, and the down-dip extent of afterslip may be unrealistic. Viscoelastic relaxation in the high-temperature Tibetan crust reproduces the observed SSW motion without overpredicting the uplift; viscoelastic relaxation in the downgoing Indian mantle, however, produces northward motion and subsidence north of the rupture, i.e. opposite to the observed motions. We argue that models of coupled afterslip (confined close to the rupture) and viscoelastic relaxation can reproduce the postseismic timeseries with physically plausible parameters.

Seismic waves and earthquakes in a global monolithic model

NASA Astrophysics Data System (ADS)

Roubíček, Tomáš

2018-03-01

The philosophy that a single "monolithic" model can "asymptotically" replace and couple in a simple elegant way several specialized models relevant on various Earth layers is presented and, in special situations, also rigorously justified. In particular, global seismicity and tectonics is coupled to capture, e.g., (here by a simplified model) ruptures of lithospheric faults generating seismic waves which then propagate through the solid-like mantle and inner core both as shear (S) or pressure (P) waves, while S-waves are suppressed in the fluidic outer core and also in the oceans. The "monolithic-type" models have the capacity to describe all the mentioned features globally in a unified way together with corresponding interfacial conditions implicitly involved, only when scaling its parameters appropriately in different Earth's layers. Coupling of seismic waves with seismic sources due to tectonic events is thus an automatic side effect. The global ansatz is here based, rather for an illustration, only on a relatively simple Jeffreys' viscoelastic damageable material at small strains whose various scaling (limits) can lead to Boger's viscoelastic fluid or even to purely elastic (inviscid) fluid. Self-induced gravity field, Coriolis, centrifugal, and tidal forces are counted in our global model, as well. The rigorous mathematical analysis as far as the existence of solutions, convergence of the mentioned scalings, and energy conservation is briefly presented.

Numerical modeling of bubble dynamics in viscoelastic media with relaxation

NASA Astrophysics Data System (ADS)

Warnez, M. T.; Johnsen, E.

2015-06-01

Cavitation occurs in a variety of non-Newtonian fluids and viscoelastic materials. The large-amplitude volumetric oscillations of cavitation bubbles give rise to high temperatures and pressures at collapse, as well as induce large and rapid deformation of the surroundings. In this work, we develop a comprehensive numerical framework for spherical bubble dynamics in isotropic media obeying a wide range of viscoelastic constitutive relationships. Our numerical approach solves the compressible Keller-Miksis equation with full thermal effects (inside and outside the bubble) when coupled to a highly generalized constitutive relationship (which allows Newtonian, Kelvin-Voigt, Zener, linear Maxwell, upper-convected Maxwell, Jeffreys, Oldroyd-B, Giesekus, and Phan-Thien-Tanner models). For the latter two models, partial differential equations (PDEs) must be solved in the surrounding medium; for the remaining models, we show that the PDEs can be reduced to ordinary differential equations. To solve the general constitutive PDEs, we present a Chebyshev spectral collocation method, which is robust even for violent collapse. Combining this numerical approach with theoretical analysis, we simulate bubble dynamics in various viscoelastic media to determine the impact of relaxation time, a constitutive parameter, on the associated physics. Relaxation time is found to increase bubble growth and permit rebounds driven purely by residual stresses in the surroundings. Different regimes of oscillations occur depending on the relaxation time.

Viscoelastic and elastomeric active matter: linear instability and nonlinear dynamics

NASA Astrophysics Data System (ADS)

Hemingway, Ewan J.; Cates, M. E.; Marchetti, M. C.; Fielding, S. M.

We consider a continuum model of active viscoelastic matter, whereby a model of an active nematic liquid-crystal is coupled to a minimal model of polymer dynamics with a viscoelastic relaxation time τc. To explore the resulting interplay between active and polymeric dynamics, we first generalise a linear stability analysis (from earlier studies without polymer) to derive criteria for the onset of spontaneous flow. Perhaps surprisingly, our results show that the spontaneous flow instability persists even for divergent polymer relaxation times. We explore the novel dynamical states to which these instabilities lead by means of nonlinear numerical simulations. This reveals oscillatory shear-banded states in 1D, and activity-driven turbulence in 2D, even in the limit τc --> ∞ . Adding polymer can also have calming effects, increasing the net throughput of spontaneous flow along a channel in a new type of ''drag-reduction'', an effect that may have implications for cytoplasmic streaming processes within the cell.

Broadscale postseismic gravity change following the 2011 Tohoku-Oki earthquake and implication for deformation by viscoelastic relaxation and afterslip.

PubMed

Han, Shin-Chan; Sauber, Jeanne; Pollitz, Fred

2014-08-28

The analysis of GRACE gravity data revealed postseismic gravity increase by 6 μGal over a 500 km scale within a couple of years after the 2011 Tohoku-Oki earthquake, which is nearly 40-50% of the coseismic gravity change. It originates mostly from changes in the isotropic component corresponding to the M rr moment tensor element. The exponential decay with rapid change in a year and gradual change afterward is a characteristic temporal pattern. Both viscoelastic relaxation and afterslip models produce reasonable agreement with the GRACE free-air gravity observation, while their Bouguer gravity patterns and seafloor vertical deformations are distinctly different. The postseismic gravity variation is best modeled by the biviscous relaxation with a transient and steady state viscosity of 10 18 and 10 19  Pa s, respectively, for the asthenosphere. Our calculated higher-resolution viscoelastic relaxation model, underlying the partially ruptured elastic lithosphere, yields the localized postseismic subsidence above the hypocenter reported from the GPS-acoustic seafloor surveying.

Broadscale postseismic gravity change following the 2011 Tohoku-Oki earthquake and implication for deformation by viscoelastic relaxation and afterslip

PubMed Central

Han, Shin-Chan; Sauber, Jeanne; Pollitz, Fred

2014-01-01

The analysis of GRACE gravity data revealed postseismic gravity increase by 6 μGal over a 500 km scale within a couple of years after the 2011 Tohoku-Oki earthquake, which is nearly 40–50% of the coseismic gravity change. It originates mostly from changes in the isotropic component corresponding to the Mrr moment tensor element. The exponential decay with rapid change in a year and gradual change afterward is a characteristic temporal pattern. Both viscoelastic relaxation and afterslip models produce reasonable agreement with the GRACE free-air gravity observation, while their Bouguer gravity patterns and seafloor vertical deformations are distinctly different. The postseismic gravity variation is best modeled by the biviscous relaxation with a transient and steady state viscosity of 1018 and 1019 Pa s, respectively, for the asthenosphere. Our calculated higher-resolution viscoelastic relaxation model, underlying the partially ruptured elastic lithosphere, yields the localized postseismic subsidence above the hypocenter reported from the GPS-acoustic seafloor surveying. PMID:25821272

Nonlinear Viscoelastic Rheology and the Occurrence of Aftershocks

NASA Astrophysics Data System (ADS)

Shcherbakov, R.; Zhang, X.

2017-12-01

Aftershocks are ubiquitous in nature. They are the manifestation of relaxation phenomena observed in various physical systems. In one prominent example, they typically occur after large earthquakes. The observed aftershock sequences usually obey several well defined non-trivial empirical laws in magnitude, temporal, and spatial domains. In many cases their characteristics follow scale-invariant distributions. The occurrence of aftershocks displays a prominent temporal behavior due to time-dependent mechanisms of stress and/or energy transfer. There are compelling evidences that the lower continental crust and upper mantle are governed by various solid state creep mechanisms. Among those mechanisms a power-law viscous flow was suggested to explain the postseismic surface deformation after large earthquakes. In this work, we consider a slider-block model to mimic the behavior of a seismogenic fault. In the model, we introduce a nonlinear viscoelastic coupling mechanism to capture the essential characteristics of crustal rheology and stress interaction between the blocks and the medium. For this purpose we employ nonlinear Kelvin-Voigt elements consisting of an elastic spring and a dashpot assembled in parallel to introduce viscoelastic coupling between the blocks and the driving plate. By mapping the model into a cellular automaton we derive the functional form of the stress transfer mechanism in the model. We show that the nonlinear viscoelasticity plays a critical role in triggering of aftershocks. It explains the functional form of the Omori-Utsu law and gives physical interpretation of its parameters. The proposed model also suggests that the power-law rheology of the fault gauge and underlying lower crust and upper mantle controls the decay rate of aftershocks. To verify this, we analyze several prominent aftershock sequences to estimate their decay rates and correlate with the rheological properties of the underlying lower crust and mantle, which were estimated from the postseismic surface deformation. Our modelling suggests that the power-law rheology exponent n controls the decay rate of aftershocks and is related to the parameter p of the Omori-Utsu law.

Finite volume multigrid method of the planar contraction flow of a viscoelastic fluid

NASA Astrophysics Data System (ADS)

Moatssime, H. Al; Esselaoui, D.; Hakim, A.; Raghay, S.

2001-08-01

This paper reports on a numerical algorithm for the steady flow of viscoelastic fluid. The conservative and constitutive equations are solved using the finite volume method (FVM) with a hybrid scheme for the velocities and first-order upwind approximation for the viscoelastic stress. A non-uniform staggered grid system is used. The iterative SIMPLE algorithm is employed to relax the coupled momentum and continuity equations. The non-linear algebraic equations over the flow domain are solved iteratively by the symmetrical coupled Gauss-Seidel (SCGS) method. In both, the full approximation storage (FAS) multigrid algorithm is used. An Oldroyd-B fluid model was selected for the calculation. Results are reported for planar 4:1 abrupt contraction at various Weissenberg numbers. The solutions are found to be stable and smooth. The solutions show that at high Weissenberg number the domain must be long enough. The convergence of the method has been verified with grid refinement. All the calculations have been performed on a PC equipped with a Pentium III processor at 550 MHz. Copyright

Viscoelastic coupling model of the San Andreas fault along the big bend, southern California

USGS Publications Warehouse

Savage, J.C.; Lisowski, M.

1997-01-01

The big bend segment of the San Andreas fault is the 300-km-long segment in southern California that strikes about N65??W, roughly 25?? counterclockwise from the local tangent to the small circle about the Pacific-North America pole of rotation. The broad distribution of deformation of trilateration networks along this segment implies a locking depth of at least 25 km as interpreted by the conventional model of strain accumulation (continuous slip on the fault below the locking depth at the rate of relative plate motion), whereas the observed seismicity and laboratory data on fault strength suggest that the locking depth should be no greater than 10 to 15 km. The discrepancy is explained by the viscoelastic coupling model which accounts for the viscoelastic response of the lower crust. Thus the broad distribution of deformation observed across the big bend segment can be largely associated with the San Andreas fault itself, not subsidiary faults distributed throughout the region. The Working Group on California Earthquake Probabilities [1995] in using geodetic data to estimate the seismic risk in southern California has assumed that strain accumulated off the San Andreas fault is released by earthquakes located off the San Andreas fault. Thus they count the San Andreas contribution to total seismic moment accumulation more than once, leading to an overestimate of the seismicity for magnitude 6 and greater earthquakes in their Type C zones.

Viscoelastic and elastomeric active matter: Linear instability and nonlinear dynamics.

PubMed

Hemingway, E J; Cates, M E; Fielding, S M

2016-03-01

We consider a continuum model of active viscoelastic matter, whereby an active nematic liquid crystal is coupled to a minimal model of polymer dynamics with a viscoelastic relaxation time τ(C). To explore the resulting interplay between active and polymeric dynamics, we first generalize a linear stability analysis (from earlier studies without polymer) to derive criteria for the onset of spontaneous heterogeneous flows (strain rate) and/or deformations (strain). We find two modes of instability. The first is a viscous mode, associated with strain rate perturbations. It dominates for relatively small values of τ(C) and is a simple generalization of the instability known previously without polymer. The second is an elastomeric mode, associated with strain perturbations, which dominates at large τ(C) and persists even as τ(C)→∞. We explore the dynamical states to which these instabilities lead by means of direct numerical simulations. These reveal oscillatory shear-banded states in one dimension and activity-driven turbulence in two dimensions even in the elastomeric limit τ(C)→∞. Adding polymer can also have calming effects, increasing the net throughput of spontaneous flow along a channel in a type of drag reduction. The effect of including strong antagonistic coupling between the nematic and polymer is examined numerically, revealing a rich array of spontaneously flowing states.

Viscoelastic and elastomeric active matter: Linear instability and nonlinear dynamics

NASA Astrophysics Data System (ADS)

Hemingway, E. J.; Cates, M. E.; Fielding, S. M.

2016-03-01

We consider a continuum model of active viscoelastic matter, whereby an active nematic liquid crystal is coupled to a minimal model of polymer dynamics with a viscoelastic relaxation time τC. To explore the resulting interplay between active and polymeric dynamics, we first generalize a linear stability analysis (from earlier studies without polymer) to derive criteria for the onset of spontaneous heterogeneous flows (strain rate) and/or deformations (strain). We find two modes of instability. The first is a viscous mode, associated with strain rate perturbations. It dominates for relatively small values of τC and is a simple generalization of the instability known previously without polymer. The second is an elastomeric mode, associated with strain perturbations, which dominates at large τC and persists even as τC→∞ . We explore the dynamical states to which these instabilities lead by means of direct numerical simulations. These reveal oscillatory shear-banded states in one dimension and activity-driven turbulence in two dimensions even in the elastomeric limit τC→∞ . Adding polymer can also have calming effects, increasing the net throughput of spontaneous flow along a channel in a type of drag reduction. The effect of including strong antagonistic coupling between the nematic and polymer is examined numerically, revealing a rich array of spontaneously flowing states.

Bubble dynamics in viscoelastic soft tissue in high-intensity focal ultrasound thermal therapy.

PubMed

Zilonova, E; Solovchuk, M; Sheu, T W H

2018-01-01

The present study is aimed to investigate bubble dynamics in a soft tissue, to which HIFU's continuous harmonic pulse is applied by introducing a viscoelastic cavitation model. After a comparison of some existing cavitation models, we decided to employ Gilmore-Akulichev model. This chosen cavitation model should be coupled with the Zener viscoelastic model in order to be able to simulate soft tissue features such as elasticity and relaxation time. The proposed Gilmore-Akulichev-Zener model was investigated for exploring cavitation dynamics. The parametric study led us to the conclusion that the elasticity and viscosity both damp bubble oscillations, whereas the relaxation effect depends mainly on the period of the ultrasound wave. The similar influence of elasticity, viscosity and relaxation time on the temperature inside the bubble can be observed. Cavitation heat source terms (corresponding to viscous damping and pressure wave radiated by bubble collapse) were obtained based on the proposed model to examine the cavitation significance during the treatment process. Their maximum values both overdominate the acoustic ultrasound term in HIFU applications. Elasticity was revealed to damp a certain amount of deposited heat for both cavitation terms. Copyright © 2017 Elsevier B.V. All rights reserved.

A spatially nonlocal model for polymer-penetrant diffusion

NASA Astrophysics Data System (ADS)

Edwards, D. A.

Diffusion of a penetrant in a polymer entanglement network cannot be described by Fick's Law alone; rather, one must incorporate other nonlocal effects. In contrast to previous viscoelastic models which have modeled these effects through hereditary integrals in time, a new model is presented exploiting the disparate lengths of the polymer in the glassy (dry) and rubbery (saturated) states. This model leads to a partial integrodifferential equation which is nonlocal in space. The system is recast as a moving boundary-value problem between sets of coupled partial differential equations. Using singular perturbation techniques, sorption in a semi-infinite polymer is studied on several time scales with varying exposed interface conditions. Though some of the results match with those from viscoelastic models, new physically relevant behaviors also appear. These include the formation of stopping fronts and overshoot in the pseudostress.

Numerical modeling of bubble dynamics in viscoelastic media with relaxation

PubMed Central

Warnez, M. T.; Johnsen, E.

2015-01-01

Cavitation occurs in a variety of non-Newtonian fluids and viscoelastic materials. The large-amplitude volumetric oscillations of cavitation bubbles give rise to high temperatures and pressures at collapse, as well as induce large and rapid deformation of the surroundings. In this work, we develop a comprehensive numerical framework for spherical bubble dynamics in isotropic media obeying a wide range of viscoelastic constitutive relationships. Our numerical approach solves the compressible Keller–Miksis equation with full thermal effects (inside and outside the bubble) when coupled to a highly generalized constitutive relationship (which allows Newtonian, Kelvin–Voigt, Zener, linear Maxwell, upper-convected Maxwell, Jeffreys, Oldroyd-B, Giesekus, and Phan-Thien-Tanner models). For the latter two models, partial differential equations (PDEs) must be solved in the surrounding medium; for the remaining models, we show that the PDEs can be reduced to ordinary differential equations. To solve the general constitutive PDEs, we present a Chebyshev spectral collocation method, which is robust even for violent collapse. Combining this numerical approach with theoretical analysis, we simulate bubble dynamics in various viscoelastic media to determine the impact of relaxation time, a constitutive parameter, on the associated physics. Relaxation time is found to increase bubble growth and permit rebounds driven purely by residual stresses in the surroundings. Different regimes of oscillations occur depending on the relaxation time. PMID:26130967

Actuation of flexoelectric membranes in viscoelastic fluids with applications to outer hair cells

PubMed Central

Herrera-Valencia, E. E.; Rey, Alejandro D.

2014-01-01

Liquid crystal flexoelectric actuation uses an imposed electric field to create membrane bending, and it is used by the outer hair cells (OHCs) located in the inner ear, whose role is to amplify sound through generation of mechanical power. Oscillations in the OHC membranes create periodic viscoelastic flows in the contacting fluid media. A key objective of this work on flexoelectric actuation relevant to OHCs is to find the relations and impact of the electromechanical properties of the membrane, the rheological properties of the viscoelastic media, and the frequency response of the generated mechanical power output. The model developed and used in this work is based on the integration of: (i) the flexoelectric membrane shape equation applied to a circular membrane attached to the inner surface of a circular capillary and (ii) the coupled capillary flow of contacting viscoelastic phases, such that the membrane flexoelectric oscillations drive periodic viscoelastic capillary flows, as in OHCs. By applying the Fourier transform formalism to the governing equation, analytical expressions for the transfer function associated with the curvature and electrical field and for the power dissipation of elastic storage energy were found. PMID:25332388

Broadscale Postseismic Gravity Change Following the 2011 Tohoku-Oki Earthquake and Implication for Deformation by Viscoelastic Relaxation and Afterslip

NASA Technical Reports Server (NTRS)

Han, Shin-Chan; Sauber, Jeanne; Pollitz, Fred

2014-01-01

The analysis of GRACE gravity data revealed post-seismic gravity increase by 6 micro-Gal over a 500 km scale within a couple of years after the 2011 Tohoku-Oki earthquake, which is nearly 40-50% of the co-seismic gravity change. It originates mostly from changes in the isotropic component corresponding to the M(sub rr) moment tensor element. The exponential decay with rapid change in a year and gradual change afterward is a characteristic temporal pattern. Both viscoelastic relaxation and afterslip models produce reasonable agreement with the GRACE free-air gravity observation, while their Bouguer gravity patterns and seafloor vertical deformations are distinctly different. The post-seismic gravity variation is best modeled by the bi-viscous relaxation with a transient and steady state viscosity of 10(exp 18) and 10(exp 19) Pa s, respectively, for the asthenosphere. Our calculated higher-resolution viscoelastic relaxation model, underlying the partially ruptured elastic lithosphere, yields the localized post-seismic subsidence above the hypocenter reported from the GPS-acoustic seafloor surveying.

Three dimensional modelling of earthquake rupture cycles on frictional faults

NASA Astrophysics Data System (ADS)

Simpson, Guy; May, Dave

2017-04-01

We are developing an efficient MPI-parallel numerical method to simulate earthquake sequences on preexisting faults embedding within a three dimensional viscoelastic half-space. We solve the velocity form of the elasto(visco)dynamic equations using a continuous Galerkin Finite Element Method on an unstructured pentahedral mesh, which thus permits local spatial refinement in the vicinity of the fault. Friction sliding is coupled to the viscoelastic solid via rate- and state-dependent friction laws using the split-node technique. Our coupled formulation employs a picard-type non-linear solver with a fully implicit, first order accurate time integrator that utilises an adaptive time step that efficiently evolves the system through multiple seismic cycles. The implementation leverages advanced parallel solvers, preconditioners and linear algebra from the Portable Extensible Toolkit for Scientific Computing (PETSc) library. The model can treat heterogeneous frictional properties and stress states on the fault and surrounding solid as well as non-planar fault geometries. Preliminary tests show that the model successfully reproduces dynamic rupture on a vertical strike-slip fault in a half-space governed by rate-state friction with the ageing law.

Linear oscillation of gas bubbles in a viscoelastic material under ultrasound irradiation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hamaguchi, Fumiya; Ando, Keita, E-mail: kando@mech.keio.ac.jp

2015-11-15

Acoustically forced oscillation of spherical gas bubbles in a viscoelastic material is studied through comparisons between experiments and linear theory. An experimental setup has been designed to visualize bubble dynamics in gelatin gels using a high-speed camera. A spherical gas bubble is created by focusing an infrared laser pulse into (gas-supersaturated) gelatin gels. The bubble radius (up to 150 μm) under mechanical equilibrium is controlled by gradual mass transfer of gases across the bubble interface. The linearized bubble dynamics are studied from the observation of spherical bubble oscillation driven by low-intensity, planar ultrasound driven at 28 kHz. It follows frommore » the experiment for an isolated bubble that the frequency response in its volumetric oscillation was shifted to the high frequency side and its peak was suppressed as the gelatin concentration increases. The measurement is fitted to the linearized Rayleigh–Plesset equation coupled with the Voigt constitutive equation that models the behavior of linear viscoelastic solids; the fitting yields good agreement by tuning unknown values of the viscosity and rigidity, indicating that more complex phenomena including shear thinning, stress relaxation, and retardation do not play an important role for the small-amplitude oscillations. Moreover, the cases for bubble-bubble and bubble-wall systems are studied. The observed interaction effect on the linearized dynamics can be explained as well by a set of the Rayleigh–Plesset equations coupled through acoustic radiation among these systems. This suggests that this experimental setup can be applied to validate the model of bubble dynamics with more complex configuration such as a cloud of bubbles in viscoelastic materials.« less

A Comparison of Curing Process-Induced Residual Stresses and Cure Shrinkage in Micro-Scale Composite Structures with Different Constitutive Laws

NASA Astrophysics Data System (ADS)

Li, Dongna; Li, Xudong; Dai, Jianfeng; Xi, Shangbin

2018-02-01

In this paper, three kinds of constitutive laws, elastic, "cure hardening instantaneously linear elastic (CHILE)" and viscoelastic law, are used to predict curing process-induced residual stress for the thermoset polymer composites. A multi-physics coupling finite element analysis (FEA) model implementing the proposed three approaches is established in COMSOL Multiphysics-Version 4.3b. The evolution of thermo-physical properties with temperature and degree of cure (DOC), which improved the accuracy of numerical simulations, and cure shrinkage are taken into account for the three models. Subsequently, these three proposed constitutive models are implemented respectively in a 3D micro-scale composite laminate structure. Compared the differences between these three numerical results, it indicates that big error in residual stress and cure shrinkage generates by elastic model, but the results calculated by the modified CHILE model are in excellent agreement with those estimated by the viscoelastic model.

Transport behaviors of locally fractional coupled Brownian motors with fluctuating interactions

NASA Astrophysics Data System (ADS)

Wang, Huiqi; Ni, Feixiang; Lin, Lifeng; Lv, Wangyong; Zhu, Hongqiang

2018-09-01

In some complex viscoelastic mediums, it is ubiquitous that absorbing and desorbing surrounding Brownian particles randomly occur in coupled systems. The conventional method is to model a variable-mass system driven by both multiplicative and additive noises. In this paper, an improved mathematical model is created based on generalized Langevin equations (GLE) to characterize the random interaction with locally fluctuating number of coupled particles in the elastically coupled factional Brownian motors (FBM). By the numerical simulations, the effect of fluctuating interactions on collective transport behaviors is investigated, and some abnormal phenomena, such as cooperative behaviors, stochastic resonance (SR) and anomalous transport, are observed in the regime of sub-diffusion.

Seismic Full Waveform Modeling & Imaging in Attenuating Media

NASA Astrophysics Data System (ADS)

Guo, Peng

Seismic attenuation strongly affects seismic waveforms by amplitude loss and velocity dispersion. Without proper inclusion of Q parameters, errors can be introduced for seismic full waveform modeling and imaging. Three different (Carcione's, Robertsson's, and the generalized Robertsson's) isotropic viscoelastic wave equations based on the generalized standard linear solid (GSLS) are evaluated. The second-order displacement equations are derived, and used to demonstrate that, with the same stress relaxation times, these viscoelastic formulations are equivalent. By introducing separate memory variables for P and S relaxation functions, Robertsson's formulation is generalized to allow different P and S wave stress relaxation times, which improves the physical consistency of the Qp and Qs modelled in the seismograms.The three formulations have comparable computational cost. 3D seismic finite-difference forward modeling is applied to anisotropic viscoelastic media. The viscoelastic T-matrix (a dynamic effective medium theory) relates frequency-dependent anisotropic attenuation and velocity to reservoir properties in fractured HTI media, based on the meso-scale fluid flow attenuation mechanism. The seismic signatures resulting from changing viscoelastic reservoir properties are easily visible. Analysis of 3D viscoelastic seismograms suggests that anisotropic attenuation is a potential tool for reservoir characterization. To compensate the Q effects during reverse-time migration (RTM) in viscoacoustic and viscoelastic media, amplitudes need to be compensated during wave propagation; the propagation velocity of the Q-compensated wavefield needs to be the same as in the attenuating wavefield, to restore the phase information. Both amplitude and phase can be compensated when the velocity dispersion and the amplitude loss are decoupled. For wave equations based on the GSLS, because Q effects are coupled in the memory variables, Q-compensated wavefield propagates faster than the attenuating wavefield, and introduce unwanted phase shift. Numerical examples show that there are phase (depth) shifts in the Q-compensated RTM images from the GSLS equation. An adjoint-based least-squares reverse-time migration is proposed for viscoelastic media (Q-LSRTM), to compensate the attenuation losses in P and S images. The viscoelastic adjoint operator, and the P and S modulus perturbation imaging conditions are derived using the adjoint-state method and an augmented Lagrangian functional. Q-LSRTM solves the viscoelastic linearized modeling operator for synthetic data, and the adjoint operator is used for back propagating the data residual. Q-LSRTM is capable of iteratively updating the P and S modulus perturbations,in the direction of minimizing data residuals, and attenuation loss is iteratively compensated. A novel Q compensation approach is developed for adjoint seismic imaging by pseudodifferential scaling. With a correct Q model included in the migration algorithm, propagation effects, including the Q effects, can be compensated with the application of the inverse Hessian to the RTM image. Pseudodifferential scaling is used to efficiently approximate the action of the inverse Hessian. Numerical examples indicate that the adjoint RTM images with pseudodifferential scaling approximate the true model perturbation, and can be used as well-conditioned gradients for least-squares imaging.

FEM simulation of the die compaction of pharmaceutical products: influence of visco-elastic phenomena and comparison with experiments.

PubMed

Diarra, Harona; Mazel, Vincent; Busignies, Virginie; Tchoreloff, Pierre

2013-09-10

This work studies the influence of visco-elastic behavior in the finite element method (FEM) modeling of die compaction of pharmaceutical products and how such a visco-elastic behavior may improve the agreement between experimental and simulated compression curves. The modeling of the process was conducted on a pharmaceutical excipient, microcrystalline cellulose (MCC), by using Drucker-Prager cap model coupled with creep behavior in Abaqus(®) software. The experimental data were obtained on a compaction simulator (STYLCAM 200R). The elastic deformation of the press was determined by performing experimental tests on a calibration disk and was introduced in the simulation. Numerical optimization was performed to characterize creep parameters. The use of creep behavior in the simulations clearly improved the agreement between the numerical and experimental compression curves (stresses, thickness), mainly during the unloading part of the compaction cycle. For the first time, it was possible to reproduce numerically the fact that the minimum tablet thickness is not obtained at the maximum compression stress. This study proves that creep behavior must be taken into account when modeling the compaction of pharmaceutical products using FEM methods. Copyright © 2013 Elsevier B.V. All rights reserved.

How Enzymes Work: A Look through the Perspective of Molecular Viscoelastic Properties

NASA Astrophysics Data System (ADS)

Qu, Hao; Zocchi, Giovanni

2013-01-01

We present nanorheology measurements on the folded state of an enzyme that show directly that the (ensemble-averaged) stress-strain relations are nonlinear and frequency dependent beyond 1-Å deformation. We argue that this frequency dependence allows for opening a nonequilibrium cycle in the force-deformation plane if the forward and backward conformational changes of the enzyme during catalysis happen at different speeds. Using a heuristic model for the experimentally established viscoelastic properties of the enzyme, we examine a number of general features of enzymatic action. We find that the proposed viscoelastic cycle is consistent with the linear decrease of the speed of motor proteins with load. We find a relation between the stall force and the maximum rate for enzymes (in general) and motors (in particular). We estimate the stall force of the motor protein kinesin from thermodynamic quantities and estimate the maximum rate of enzymes from purely mechanical quantities. We propose that the viscoelastic cycle provides a framework for considering mechanochemical coupling in enzymes on the basis of possibly universal materials properties of the folded state of proteins.

Acoustic radiation force induced resonance elastography of coagulating blood: theoretical viscoelasticity modeling and ex vivo experimentation

NASA Astrophysics Data System (ADS)

Bhatt, Manish; Montagnon, Emmanuel; Destrempes, François; Chayer, Boris; Kazemirad, Siavash; Cloutier, Guy

2018-03-01

Deep vein thrombosis is a common vascular disease that can lead to pulmonary embolism and death. The early diagnosis and clot age staging are important parameters for reliable therapy planning. This article presents an acoustic radiation force induced resonance elastography method for the viscoelastic characterization of clotting blood. The physical concept of this method relies on the mechanical resonance of the blood clot occurring at specific frequencies. Resonances are induced by focusing ultrasound beams inside the sample under investigation. Coupled to an analytical model of wave scattering, the ability of the proposed method to characterize the viscoelasticity of a mimicked venous thrombosis in the acute phase is demonstrated. Experiments with a gelatin-agar inclusion sample of known viscoelasticity are performed for validation and establishment of the proof of concept. In addition, an inversion method is applied in vitro for the kinetic monitoring of the blood coagulation process of six human blood samples obtained from two volunteers. The computed elasticity and viscosity values of blood samples at the end of the 90 min kinetics were estimated at 411  ±  71 Pa and 0.25  ±  0.03 Pa · s for volunteer #1, and 387  ±  35 Pa and 0.23  ±  0.02 Pa · s for volunteer #2, respectively. The proposed method allowed reproducible time-varying thrombus viscoelastic measurements from samples having physiological dimensions.

Spatio-temporal dynamics of an active, polar, viscoelastic ring.

PubMed

Marcq, Philippe

2014-04-01

Constitutive equations for a one-dimensional, active, polar, viscoelastic liquid are derived by treating the strain field as a slow hydrodynamic variable. Taking into account the couplings between strain and polarity allowed by symmetry, the hydrodynamics of an active, polar, viscoelastic body include an evolution equation for the polarity field that generalizes the damped Kuramoto-Sivashinsky equation. Beyond thresholds of the active coupling coefficients between the polarity and the stress or the strain rate, bifurcations of the homogeneous state lead first to stationary waves, then to propagating waves of the strain, stress and polarity fields. I argue that these results are relevant to living matter, and may explain rotating actomyosin rings in cells and mechanical waves in epithelial cell monolayers.

Thermo-Mechanical Analyses of Dynamically Loaded Rubber Cylinders

NASA Technical Reports Server (NTRS)

Johnson, Arthur R.; Chen, Tzi-Kang

2002-01-01

Thick rubber components are employed by the Army to carry large loads. In tanks, rubber covers road wheels and track systems to protect roadways. It is difficult for design engineers to simulate the details of the hysteretic heating for large strain viscoelastic deformations. In this study, an approximation to the viscoelastic energy dissipated per unit time is investigated for use in estimating mechanically induced viscoelastic heating. Coupled thermo-mechanical simulations of large cyclic deformations of rubber cylinders are presented. The cylinders are first compressed axially and then cyclically loaded about the compressed state. Details of the algorithm and some computational issues are discussed. The coupled analyses are conducted for tall and short rubber cylinders both with and without imbedded metal disks.

Influence of internal viscoelastic modes on the Brownian motion of a λ-DNA coated colloid.

PubMed

Yanagishima, Taiki; Laohakunakorn, Nadanai; Keyser, Ulrich F; Eiser, Erika; Tanaka, Hajime

2014-03-21

We study the influence of grafted polymers on the diffusive behaviour of a colloidal particle. Our work demonstrates how such additional degrees of freedom influence the Brownian motion of the particle, focusing on internal viscoelastic coupling between the polymer and colloid. Specifically, we study the mean-squared displacements (MSDs) of λ-DNA grafted colloids using Brownian dynamics simulation. Our simulations reveal the non-trivial effect of internal modes, which gives rise to a crossover from the short-time viscoelastic to long-time diffusional behaviour. We also show that basic features can be captured by a simple theoretical model considering the relative motion of a colloid to a part of the polymer corona. This model describes well a MSD calculated from an extremely long trajectory of a single λ-DNA coated colloid from experiment and allows characterisation of the λ-DNA hairs. Our study suggests that the access to the internal relaxation modes via the colloid trajectory offers a novel method for the characterisation of soft attachments to a colloid.

Computing dispersion curves of elastic/viscoelastic transversely-isotropic bone plates coupled with soft tissue and marrow using semi-analytical finite element (SAFE) method.

PubMed

Nguyen, Vu-Hieu; Tran, Tho N H T; Sacchi, Mauricio D; Naili, Salah; Le, Lawrence H

2017-08-01

We present a semi-analytical finite element (SAFE) scheme for accurately computing the velocity dispersion and attenuation in a trilayered system consisting of a transversely-isotropic (TI) cortical bone plate sandwiched between the soft tissue and marrow layers. The soft tissue and marrow are mimicked by two fluid layers of finite thickness. A Kelvin-Voigt model accounts for the absorption of all three biological domains. The simulated dispersion curves are validated by the results from the commercial software DISPERSE and published literature. Finally, the algorithm is applied to a viscoelastic trilayered TI bone model to interpret the guided modes of an ex-vivo experimental data set from a bone phantom. Copyright © 2017 Elsevier Ltd. All rights reserved.

Rheology of the lithosphere inferred from postseismic uplift following the 1959 Hebgen Lake earthquake

USGS Publications Warehouse

Nishimura, T.; Thatcher, W.

2003-01-01

We have modeled the broad postseismic uplift measured by geodetic leveling in the epicentral area of the 1959 Mw = 7.3 Hebgen Lake, Montana earthquake, a normal faulting event in the northern Basin and Range province. To fit the observed uplift we calculate synthetic postseismic deformation using the relaxation response of a gravitational viscoelastic Earth to the earthquake. For a model with an elastic plate overlying a viscoelastic half-space, we find that the elastic thickness is 38 ?? 8 km, which isclose to the local crustal thickness. The half-space viscosity is estimated at 4 ?? 1018??0.5 Pa s. The leveling data do not require a viscous lower crust but permit a lower bound viscosity of 1020 Pa s. The observed broad uplift cannot be explained by physically plausible afterslip on and below the coseismic fault. However, local deformation across the coseismic surface rupture requires shallow afterslip reaching the surface. The postseismic deformation induced by the estimated viscoelastic structure decays exponentially with a time constant of ???15 years. Because of coupling between the elastic layer and the viscoelastic substrate, this relaxation time is significantly longer than the 2 year Maxwell relaxation time of the viscous half-space itself. Our result suggests the importance of postseismic relaxation in interpreting high-precision global positioning system velocities. For example, our model results suggest that postseismic transient velocities from both the 1959 Hebgen Lake and the 1983 Mw = 6.9 Borah Peak earthquakes are currently as large as 1-2 mm/yr.

Phonation threshold pressure predictions using viscoelastic properties up to 1,400 Hz of injectables intended for Reinke's space.

PubMed

Klemuk, Sarah A; Lu, Xiaoying; Hoffman, Henry T; Titze, Ingo R

2010-05-01

Viscoelastic properties of numerous vocal fold injectables have been reported but not at speaking frequencies. For materials intended for Reinke's space, ramifications of property values are of great concern because of their impact on ease of voice onset. Our objectives were: 1) to measure viscoelastic properties of a new nonresorbing carbomer and well-known vocal fold injectables at vocalization frequencies using established and new instrumentation, and 2) to predict phonation threshold pressures using a computer model with intended placement in Reinke's space. Rheology and phonation threshold pressure calculations. Injectables were evaluated with a traditional rotational rheometer and a new piezo-rotary vibrator. Using these data at vocalization frequencies, phonation threshold pressures (PTP) were calculated for each biomaterial, assuming a low dimensional model with supraglottic coupling and adjusted vocal fold length and thickness at each frequency. Results were normalized to a nominal PTP value. Viscoelastic data were acquired at vocalization frequencies as high as 363 to 1,400 Hz for six new carbomer hydrogels, Hylan B, and Extracel intended for vocal fold Reinke's space injection and for Cymetra (lateral injection). Reliability was confirmed with good data overlap when measuring with either rheometer. PTP predictions ranged from 0.001 to 16 times the nominal PTP value of 0.283 kPa. Accurate viscoelastic measurements of vocal fold injectables are now possible at physiologic frequencies. Hylan B, Extracel, and the new carbomer hydrogels should generate easy vocal onset and sustainable vocalization based on their rheologic properties if injected into Reinke's space. Applications may vary depending on desired longevity of implant. Laryngoscope, 2010.

Fractional viscoelasticity in fractal and non-fractal media: Theory, experimental validation, and uncertainty analysis

NASA Astrophysics Data System (ADS)

Mashayekhi, Somayeh; Miles, Paul; Hussaini, M. Yousuff; Oates, William S.

2018-02-01

In this paper, fractional and non-fractional viscoelastic models for elastomeric materials are derived and analyzed in comparison to experimental results. The viscoelastic models are derived by expanding thermodynamic balance equations for both fractal and non-fractal media. The order of the fractional time derivative is shown to strongly affect the accuracy of the viscoelastic constitutive predictions. Model validation uses experimental data describing viscoelasticity of the dielectric elastomer Very High Bond (VHB) 4910. Since these materials are known for their broad applications in smart structures, it is important to characterize and accurately predict their behavior across a large range of time scales. Whereas integer order viscoelastic models can yield reasonable agreement with data, the model parameters often lack robustness in prediction at different deformation rates. Alternatively, fractional order models of viscoelasticity provide an alternative framework to more accurately quantify complex rate-dependent behavior. Prior research that has considered fractional order viscoelasticity lacks experimental validation and contains limited links between viscoelastic theory and fractional order derivatives. To address these issues, we use fractional order operators to experimentally validate fractional and non-fractional viscoelastic models in elastomeric solids using Bayesian uncertainty quantification. The fractional order model is found to be advantageous as predictions are significantly more accurate than integer order viscoelastic models for deformation rates spanning four orders of magnitude.

Active viscoelastic matter: from bacterial drag reduction to turbulent solids.

PubMed

Hemingway, E J; Maitra, A; Banerjee, S; Marchetti, M C; Ramaswamy, S; Fielding, S M; Cates, M E

2015-03-06

A paradigm for internally driven matter is the active nematic liquid crystal, whereby the equations of a conventional nematic are supplemented by a minimal active stress that violates time-reversal symmetry. In practice, active fluids may have not only liquid-crystalline but also viscoelastic polymer degrees of freedom. Here we explore the resulting interplay by coupling an active nematic to a minimal model of polymer rheology. We find that adding a polymer can greatly increase the complexity of spontaneous flow, but can also have calming effects, thereby increasing the net throughput of spontaneous flow along a pipe (a "drag-reduction" effect). Remarkably, active turbulence can also arise after switching on activity in a sufficiently soft elastomeric solid.

Active Viscoelastic Matter: From Bacterial Drag Reduction to Turbulent Solids

NASA Astrophysics Data System (ADS)

Hemingway, E. J.; Maitra, A.; Banerjee, S.; Marchetti, M. C.; Ramaswamy, S.; Fielding, S. M.; Cates, M. E.

2015-03-01

A paradigm for internally driven matter is the active nematic liquid crystal, whereby the equations of a conventional nematic are supplemented by a minimal active stress that violates time-reversal symmetry. In practice, active fluids may have not only liquid-crystalline but also viscoelastic polymer degrees of freedom. Here we explore the resulting interplay by coupling an active nematic to a minimal model of polymer rheology. We find that adding a polymer can greatly increase the complexity of spontaneous flow, but can also have calming effects, thereby increasing the net throughput of spontaneous flow along a pipe (a "drag-reduction" effect). Remarkably, active turbulence can also arise after switching on activity in a sufficiently soft elastomeric solid.

Experimental and Computational Studies of Sound Transmission in a Branching Airway Network Embedded in a Compliant Viscoelastic Medium

PubMed Central

Dai, Zoujun; Peng, Ying; Mansy, Hansen A.; Sandler, Richard H.; Royston, Thomas J.

2015-01-01

Breath sounds are often used to aid in the diagnosis of pulmonary disease. Mechanical and numerical models could be used to enhance our understanding of relevant sound transmission phenomena. Sound transmission in an airway mimicking phantom was investigated using a mechanical model with a branching airway network embedded in a compliant viscoelastic medium. The Horsfield self-consistent model for the bronchial tree was adopted to topologically couple the individual airway segments into the branching airway network. The acoustics of the bifurcating airway segments were measured by microphones and calculated analytically. Airway phantom surface motion was measured using scanning laser Doppler vibrometry. Finite element simulations of sound transmission in the airway phantom were performed. Good agreement was achieved between experiments and simulations. The validated computational approach can provide insight into sound transmission simulations in real lungs. PMID:26097256

Experimental and computational studies of sound transmission in a branching airway network embedded in a compliant viscoelastic medium

NASA Astrophysics Data System (ADS)

Dai, Zoujun; Peng, Ying; Mansy, Hansen A.; Sandler, Richard H.; Royston, Thomas J.

2015-03-01

Breath sounds are often used to aid in the diagnosis of pulmonary disease. Mechanical and numerical models could be used to enhance our understanding of relevant sound transmission phenomena. Sound transmission in an airway mimicking phantom was investigated using a mechanical model with a branching airway network embedded in a compliant viscoelastic medium. The Horsfield self-consistent model for the bronchial tree was adopted to topologically couple the individual airway segments into the branching airway network. The acoustics of the bifurcating airway segments were measured by microphones and calculated analytically. Airway phantom surface motion was measured using scanning laser Doppler vibrometry. Finite element simulations of sound transmission in the airway phantom were performed. Good agreement was achieved between experiments and simulations. The validated computational approach can provide insight into sound transmission simulations in real lungs.

Coupled Thermo-Mechanical Analyses of Dynamically Loaded Rubber Cylinders

NASA Technical Reports Server (NTRS)

Johnson, Arthur R.; Chen, Tzi-Kang

2000-01-01

A procedure that models coupled thermo-mechanical deformations of viscoelastic rubber cylinders by employing the ABAQUS finite element code is described. Computational simulations of hysteretic heating are presented for several tall and short rubber cylinders both with and without a steel disk at their centers. The cylinders are compressed axially and are then cyclically loaded about the compressed state. The non-uniform hysteretic heating of the rubber cylinders containing a steel disk is presented. The analyses performed suggest that the coupling procedure should be considered for further development as a design tool for rubber degradation studies.

Micromechanical Fatigue Visco-Damage Model for Short Glass Fiber Reinforced Polyamide-66

NASA Astrophysics Data System (ADS)

Despringre, N.; Chemisky, Y.; Robert, G.; Meraghni, F.

This work presents a micromechanical fatigue damage model developed for short glass fiber reinforced PA66. It has been developed to predict the high cycle fatigue behavior of PA66/GF30. The model is based on an extended Mori-Tanaka method which includes coated inclusions, matrix viscoelasticity and the evolution of micro-scale damage. The developed model accounts for the nonlinear matrix viscoelasticity and the reinforcement orientation. The description of the damage processes is based on the experimental investigation of damage mechanisms previously performed through in-situ SEM tests and X-ray micro-computed tomography observations. Damage chronologies have been proposed involving three different processes: interface debonding/coating, matrix micro-cracking and fiber breakages. Their occurrence strongly depends on the microstructure and the relative humidity. Each damage mechanism is introduced through an evolution law coupled to local stress fields. The developed model is implemented using a UMAT subroutine. Its experimental validation is achieved under stress or strain controlled fatigue tests.

A dynamic unilateral contact problem with adhesion and friction in viscoelasticity

NASA Astrophysics Data System (ADS)

Cocou, Marius; Schryve, Mathieu; Raous, Michel

2010-08-01

The aim of this paper is to study an interaction law coupling recoverable adhesion, friction and unilateral contact between two viscoelastic bodies of Kelvin-Voigt type. A dynamic contact problem with adhesion and nonlocal friction is considered and its variational formulation is written as the coupling between an implicit variational inequality and a parabolic variational inequality describing the evolution of the intensity of adhesion. The existence and approximation of variational solutions are analysed, based on a penalty method, some abstract results and compactness properties. Finally, some numerical examples are presented.

Intracellular Fluid Mechanics: Coupling Cytoplasmic Flow with Active Cytoskeletal Gel

NASA Astrophysics Data System (ADS)

Mogilner, Alex; Manhart, Angelika

2018-01-01

The cell is a mechanical machine, and continuum mechanics of the fluid cytoplasm and the viscoelastic deforming cytoskeleton play key roles in cell physiology. We review mathematical models of intracellular fluid mechanics, from cytoplasmic fluid flows, to the flow of a viscous active cytoskeletal gel, to models of two-phase poroviscous flows, to poroelastic models. We discuss application of these models to cell biological phenomena, such as organelle positioning, blebbing, and cell motility. We also discuss challenges of understanding fluid mechanics on the cellular scale.

Nanoscale effects in the characterization of viscoelastic materials with atomic force microscopy: coupling of a quasi-three-dimensional standard linear solid model with in-plane surface interactions.

PubMed

Solares, Santiago D

2016-01-01

Significant progress has been accomplished in the development of experimental contact-mode and dynamic-mode atomic force microscopy (AFM) methods designed to measure surface material properties. However, current methods are based on one-dimensional (1D) descriptions of the tip-sample interaction forces, thus neglecting the intricacies involved in the material behavior of complex samples (such as soft viscoelastic materials) as well as the differences in material response between the surface and the bulk. In order to begin to address this gap, a computational study is presented where the sample is simulated using an enhanced version of a recently introduced model that treats the surface as a collection of standard-linear-solid viscoelastic elements. The enhanced model introduces in-plane surface elastic forces that can be approximately related to a two-dimensional (2D) Young's modulus. Relevant cases are discussed for single- and multifrequency intermittent-contact AFM imaging, with focus on the calculated surface indentation profiles and tip-sample interaction force curves, as well as their implications with regards to experimental interpretation. A variety of phenomena are examined in detail, which highlight the need for further development of more physically accurate sample models that are specifically designed for AFM simulation. A multifrequency AFM simulation tool based on the above sample model is provided as supporting information.

Experimental validation of solid rocket motor damping models

NASA Astrophysics Data System (ADS)

Riso, Cristina; Fransen, Sebastiaan; Mastroddi, Franco; Coppotelli, Giuliano; Trequattrini, Francesco; De Vivo, Alessio

2017-12-01

In design and certification of spacecraft, payload/launcher coupled load analyses are performed to simulate the satellite dynamic environment. To obtain accurate predictions, the system damping properties must be properly taken into account in the finite element model used for coupled load analysis. This is typically done using a structural damping characterization in the frequency domain, which is not applicable in the time domain. Therefore, the structural damping matrix of the system must be converted into an equivalent viscous damping matrix when a transient coupled load analysis is performed. This paper focuses on the validation of equivalent viscous damping methods for dynamically condensed finite element models via correlation with experimental data for a realistic structure representative of a slender launch vehicle with solid rocket motors. A second scope of the paper is to investigate how to conveniently choose a single combination of Young's modulus and structural damping coefficient—complex Young's modulus—to approximate the viscoelastic behavior of a solid propellant material in the frequency band of interest for coupled load analysis. A scaled-down test article inspired to the Z9-ignition Vega launcher configuration is designed, manufactured, and experimentally tested to obtain data for validation of the equivalent viscous damping methods. The Z9-like component of the test article is filled with a viscoelastic material representative of the Z9 solid propellant that is also preliminarily tested to investigate the dependency of the complex Young's modulus on the excitation frequency and provide data for the test article finite element model. Experimental results from seismic and shock tests performed on the test configuration are correlated with numerical results from frequency and time domain analyses carried out on its dynamically condensed finite element model to assess the applicability of different equivalent viscous damping methods to describe damping properties of slender launch vehicles in payload/launcher coupled load analysis.

Experimental validation of solid rocket motor damping models

NASA Astrophysics Data System (ADS)

Riso, Cristina; Fransen, Sebastiaan; Mastroddi, Franco; Coppotelli, Giuliano; Trequattrini, Francesco; De Vivo, Alessio

2018-06-01

In design and certification of spacecraft, payload/launcher coupled load analyses are performed to simulate the satellite dynamic environment. To obtain accurate predictions, the system damping properties must be properly taken into account in the finite element model used for coupled load analysis. This is typically done using a structural damping characterization in the frequency domain, which is not applicable in the time domain. Therefore, the structural damping matrix of the system must be converted into an equivalent viscous damping matrix when a transient coupled load analysis is performed. This paper focuses on the validation of equivalent viscous damping methods for dynamically condensed finite element models via correlation with experimental data for a realistic structure representative of a slender launch vehicle with solid rocket motors. A second scope of the paper is to investigate how to conveniently choose a single combination of Young's modulus and structural damping coefficient—complex Young's modulus—to approximate the viscoelastic behavior of a solid propellant material in the frequency band of interest for coupled load analysis. A scaled-down test article inspired to the Z9-ignition Vega launcher configuration is designed, manufactured, and experimentally tested to obtain data for validation of the equivalent viscous damping methods. The Z9-like component of the test article is filled with a viscoelastic material representative of the Z9 solid propellant that is also preliminarily tested to investigate the dependency of the complex Young's modulus on the excitation frequency and provide data for the test article finite element model. Experimental results from seismic and shock tests performed on the test configuration are correlated with numerical results from frequency and time domain analyses carried out on its dynamically condensed finite element model to assess the applicability of different equivalent viscous damping methods to describe damping properties of slender launch vehicles in payload/launcher coupled load analysis.

Pressure-induced emergence of unusually high-frequency transverse excitations in a liquid alkali metal: Evidence of two types of collective excitations contributing to the transverse dynamics at high pressures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bryk, Taras; Lviv Polytechnic National University, 12 S. Bandera Street, UA-79013 Lviv; Ruocco, G.

Unlike phonons in crystals, the collective excitations in liquids cannot be treated as propagation of harmonic displacements of atoms around stable local energy minima. The viscoelasticity of liquids, reflected in transition from the adiabatic to elastic high-frequency speed of sound and in absence of the long-wavelength transverse excitations, results in dispersions of longitudinal (L) and transverse (T) collective excitations essentially different from the typical phonon ones. Practically, nothing is known about the effect of high pressure on the dispersion of collective excitations in liquids, which causes strong changes in liquid structure. Here dispersions of L and T collective excitations inmore » liquid Li in the range of pressures up to 186 GPa were studied by ab initio simulations. Two methodologies for dispersion calculations were used: direct estimation from the peak positions of the L/T current spectral functions and simulation-based calculations of wavenumber-dependent collective eigenmodes. It is found that at ambient pressure, the longitudinal and transverse dynamics are well separated, while at high pressures, the transverse current spectral functions, density of vibrational states, and dispersions of collective excitations yield evidence of two types of propagating modes that contribute strongly to transverse dynamics. Emergence of the unusually high-frequency transverse modes gives evidence of the breakdown of a regular viscoelastic theory of transverse dynamics, which is based on coupling of a single transverse propagating mode with shear relaxation. The explanation of the observed high-frequency shift above the viscoelastic value is given by the presence of another branch of collective excitations. With the pressure increasing, coupling between the two types of collective excitations is rationalized within a proposed extended viscoelastic model of transverse dynamics.« less

Time-Dependent Testing Evaluation and Modeling for Rubber Stopper Seal Performance.

PubMed

Zeng, Qingyu; Zhao, Xia

2018-01-01

Sufficient rubber stopper sealing performance throughout the entire sealed product life cycle is essential for maintaining container closure integrity in the parenteral packaging industry. However, prior publications have lacked systematic considerations for the time-dependent influence on sealing performance that results from the viscoelastic characteristics of the rubber stoppers. In this paper, we report results of an effort to study these effects by applying both compression stress relaxation testing and residual seal force testing for time-dependent experimental data collection. These experiments were followed by modeling fit calculations based on the Maxwell-Wiechert theory modified with the Kohlrausch-Williams-Watts stretched exponential function, resulting in a nonlinear, time-dependent sealing force model. By employing both testing evaluations and modeling calculations, an in-depth understanding of the time-dependent effects on rubber stopper sealing force was developed. Both testing and modeling data show good consistency, demonstrating that the sealing force decays exponentially over time and eventually levels off because of the viscoelastic nature of the rubber stoppers. The nonlinearity of stress relaxation derives from the viscoelastic characteristics of the rubber stoppers coupled with the large stopper compression deformation into restrained geometry conditions. The modeling fit with capability to handle actual testing data can be employed as a tool to calculate the compression stress relaxation and residual seal force throughout the entire sealed product life cycle. In addition to being time-dependent, stress relaxation is also experimentally shown to be temperature-dependent. The present work provides a new, integrated methodology framework and some fresh insights to the parenteral packaging industry for practically and proactively considering, designing, setting up, controlling, and managing stopper sealing performance throughout the entire sealed product life cycle. LAY ABSTRACT: Historical publications in the parenteral packaging industry have lacked systematic considerations for the time-dependent influence on the sealing performance that results from effects of viscoelastic characteristic of the rubber stoppers. This study applied compression stress relaxation testing and residual seal force testing for time-dependent experimental data collection. These experiments were followed by modeling fit calculations based on the Maxwell-Wiechert theory modified with the Kohlrausch-Williams-Watts stretched exponential function, resulting in a nonlinear, time-dependent sealing force model. Experimental and modeling data show good consistency, demonstrating that sealing force decays exponentially over time and eventually levels off. The nonlinearity of stress relaxation derives from the viscoelastic characteristics of the rubber stoppers coupled with the large stopper compression deformation into restrained geometry conditions. In addition to being time-dependent stress relaxation, it is also experimentally shown to be temperature-dependent. The present work provides a new, integrated methodology framework and some fresh insights to the industry for practically and proactively considering, designing, setting up, controlling, and managing of the stopper sealing performance throughout the entire sealed product life cycle. © PDA, Inc. 2018.

A theoretical and computational framework for mechanics of the cortex

NASA Astrophysics Data System (ADS)

Torres-SáNchez, Alejandro; Arroyo, Marino

The cell cortex is a thin network of actin filaments lying beneath the cell surface of animal cells. Myosin motors exert contractile forces in this network leading to active stresses, which play a key role in processes such as cytokinesis or cell migration. Thus, understanding the mechanics of the cortex is fundamental to understand the mechanics of animal cells. Due to the dynamic remodeling of the actin network, the cortex behaves as a viscoelastic fluid. Furthermore, due to the difference between its thickness (tens of nanometers) and its dimensions (tens of microns), the cortex can be regarded a surface. Thus, we can model the cortex as a viscoelastic fluid, confined to a surface, that generates active stresses. Interestingly, geometric confinement results in the coupling between shape generation and material flows. In this work we present a theoretical framework to model the mechanics of the cortex that couples elasticity, hydrodynamics and force generation. We complement our theoretical description with a computational setting to simulate the resulting non-linear equations. We use this methodology to understand different processes such as asymmetric cell division or experimental probing of the rheology of the cortex We acknowledge the support of the Europen Research Council through Grant ERC CoG-681434.

Limitations in interpretation of Quartz Crystal Microbalance (QCM) beyond the rigid (Sauerbrey) to viscoelastic (lossy) transition

NASA Astrophysics Data System (ADS)

Wiener, Clinton; Weiss, Robert; White, Christopher; Vogt, Bryan

2014-03-01

Since Sauerbrey's 1959 discovery of the mass-frequency relationship in quartz, the QCM has been utilized to probe deposited mass layers. The mass to frequency (imaginary component of the impedance) relationship breaks down when the added mass is not rigidly coupled to the sensor surface and viscous dissipation of the quartz occurs. This dissipation is important in the deposition of soft materials such as polymers or biological molecules. By using a viscoelastic model for frequency and dissipation; the mass, viscosity, and shear modulus can be accurately determined. Here, we demonstrate an additional breakdown in the coupling of the imaginary component of the impedance to the mass by simultaneous QCM-D and spectroscopic ellipsometry (SE) measurements by examination of the swelling behavior of thin physically crosslinked poly-n-isopropylacrylamide films. A film swollen beyond 3 times its dry thickness shows a frequency increase (mass loss) and dissipation increases (increasing lossy film character) on cooling, but SE results show increased swelling of the film. This behavior was found to be thickness invariant for dry thicknesses of 32 nm and greater. Modeling of this QCM-D data shows non-physical results. Scaling concepts associated with this high loss limit will be discussed.

Ultra wideband (0.5-16 kHz) MR elastography for robust shear viscoelasticity model identification.

PubMed

Liu, Yifei; Yasar, Temel K; Royston, Thomas J

2014-12-21

Changes in the viscoelastic parameters of soft biological tissues often correlate with progression of disease, trauma or injury, and response to treatment. Identifying the most appropriate viscoelastic model, then estimating and monitoring the corresponding parameters of that model can improve insight into the underlying tissue structural changes. MR Elastography (MRE) provides a quantitative method of measuring tissue viscoelasticity. In a previous study by the authors (Yasar et al 2013 Magn. Reson. Med. 70 479-89), a silicone-based phantom material was examined over the frequency range of 200 Hz-7.75 kHz using MRE, an unprecedented bandwidth at that time. Six viscoelastic models including four integer order models and two fractional order models, were fit to the wideband viscoelastic data (measured storage and loss moduli as a function of frequency). The 'fractional Voigt' model (spring and springpot in parallel) exhibited the best fit and was even able to fit the entire frequency band well when it was identified based only on a small portion of the band. This paper is an extension of that study with a wider frequency range from 500 Hz to 16 kHz. Furthermore, more fractional order viscoelastic models are added to the comparison pool. It is found that added complexity of the viscoelastic model provides only marginal improvement over the 'fractional Voigt' model. And, again, the fractional order models show significant improvement over integer order viscoelastic models that have as many or more fitting parameters.

Correlation between the viscoelastic heterogeneity and the domain wall motion of Fe-based metallic glass

NASA Astrophysics Data System (ADS)

Ouyang, S.; Song, L. J.; Liu, Y. H.; Huo, J. T.; Wang, J. Q.; Xu, W.; Li, J. L.; Wang, C. T.; Wang, X. M.; Li, R. W.

2018-06-01

The soft magnetic properties of Fe-based metallic glasses are reduced significantly by external and residual stresses, e.g., the susceptibility decreases and coercivity increases, which limits their application severely. Unraveling the micromechanism of how the stress influences the soft magnetic properties is of great help for enhancing the performance of Fe-based metallic glasses. In this work, we investigate the effect of viscoelastic heterogeneity on the motion of magnetic domain wall surrounding nanoindentations. Compared to the matrix, dissipation of the viscoelastic heterogeneity increases toward the nanoindentation. Meanwhile, the motion of domain wall under external magnetic field becomes more difficult toward the nanoindentations. A correlation between the viscoelastic dissipation and the moving ability of magnetic domain walls is observed, which can be well fitted using magnetoelastic coupling theory. This suggests that manipulating the microscale viscoelastic heterogeneity is probably a helpful strategy for enhancing the soft magnetic properties of metallic glasses.

Multiplex Particle Focusing via Hydrodynamic Force in Viscoelastic Fluids

NASA Astrophysics Data System (ADS)

Lee, Doo Jin; Brenner, Howard; Youn, Jae Ryoun; Song, Young Seok

2013-11-01

We introduce a multiplex particle focusing phenomenon that arises from the hydrodynamic interaction between the viscoelastic force and the Dean drag force in a microfluidic device. In a confined microchannel, the first normal stress difference of viscoelastic fluids results in a lateral migration of suspended particles. Such a viscoelastic force was harnessed to focus different sized particles in the middle of a microchannel, and spiral channel geometry was also considered in order to take advantage of the counteracting force, Dean drag force that induces particle migration in the outward direction. For theoretical understanding, we performed a numerical analysis of viscoelastic fluids in the spiral microfluidic channel. From these results, a concept of the `Dean-coupled Elasto-inertial Focusing band (DEF)' was proposed. This study provides in-depth physical insight into the multiplex focusing of particles that can open a new venue for microfluidic particle dynamics for a concrete high throughput platform at microscale.

Post-seismic relaxation following the great 2004 Sumatra-Andaman earthquake on a compressible self-gravitating Earth

USGS Publications Warehouse

Pollitz, F.F.; Burgmann, R.; Banerjee, P.

2006-01-01

he Mw ??? 9.0 2004 December 26 Sumatra-Andaman and Mw =8.7 2005 March 28 Nias earthquakes, which collectively ruptured approximately 1800 km of the Andaman and Sunda subduction zones, are expected to be followed by vigorous viscoelastic relaxation involving both the upper and lower mantle. Because of these large spatial dimensions it is desirable to fully account for gravitational coupling effects in the relaxation process. We present a stable method of computing relaxation of a spherically-stratified, compressible and self-gravitating viscoelastic Earth following an impulsive moment release event. The solution is cast in terms of a spherical harmonic expansion of viscoelastic normal modes. For simple layered viscoelastic models, which include a low-viscosity oceanic asthenosphere, we predict substantial post-seismic effects over a region several 100s of km wide surrounding the eastern Indian Ocean. We compare observed GPS time-series from ten regional sites (mostly in Thailand and Indonesia), beginning in 2004 December, with synthetic time-series that include the coseismic and post-seismic effects of the 2004 December 26 and 2005 March 28 earthquakes. A viscosity structure involving a biviscous (Burgers body) rheology in the asthenosphere explains the pattern and amplitude of post-seismic offsets remarkably well. ?? 2006 The Authors Journal compilation ?? 2006 RAS.

Nonlinear viscoelasticity and generalized failure criterion for biopolymer gels

NASA Astrophysics Data System (ADS)

Divoux, Thibaut; Keshavarz, Bavand; Manneville, Sébastien; McKinley, Gareth

2016-11-01

Biopolymer gels display a multiscale microstructure that is responsible for their solid-like properties. Upon external deformation, these soft viscoelastic solids exhibit a generic nonlinear mechanical response characterized by pronounced stress- or strain-stiffening prior to irreversible damage and failure, most often through macroscopic fractures. Here we show on a model acid-induced protein gel that the nonlinear viscoelastic properties of the gel can be described in terms of a 'damping function' which predicts the gel mechanical response quantitatively up to the onset of macroscopic failure. Using a nonlinear integral constitutive equation built upon the experimentally-measured damping function in conjunction with power-law linear viscoelastic response, we derive the form of the stress growth in the gel following the start up of steady shear. We also couple the shear stress response with Bailey's durability criteria for brittle solids in order to predict the critical values of the stress σc and strain γc for failure of the gel, and how they scale with the applied shear rate. This provides a generalized failure criterion for biopolymer gels in a range of different deformation histories. This work was funded by the MIT-France seed fund and by the CNRS PICS-USA scheme (#36939). BK acknowledges financial support from Axalta Coating Systems.

A coupled sharp-interface immersed boundary-finite-element method for flow-structure interaction with application to human phonation.

PubMed

Zheng, X; Xue, Q; Mittal, R; Beilamowicz, S

2010-11-01

A new flow-structure interaction method is presented, which couples a sharp-interface immersed boundary method flow solver with a finite-element method based solid dynamics solver. The coupled method provides robust and high-fidelity solution for complex flow-structure interaction (FSI) problems such as those involving three-dimensional flow and viscoelastic solids. The FSI solver is used to simulate flow-induced vibrations of the vocal folds during phonation. Both two- and three-dimensional models have been examined and qualitative, as well as quantitative comparisons, have been made with established results in order to validate the solver. The solver is used to study the onset of phonation in a two-dimensional laryngeal model and the dynamics of the glottal jet in a three-dimensional model and results from these studies are also presented.

Resolving Dynamic Properties of Polymers through Coarse-Grained Computational Studies

DOE Office of Scientific and Technical Information (OSTI.GOV)

Salerno, K. Michael; Agrawal, Anupriya; Perahia, Dvora

2016-02-05

Coupled length and time scales determine the dynamic behavior of polymers and underlie their unique viscoelastic properties. To resolve the long-time dynamics it is imperative to determine which time and length scales must be correctly modeled. In this paper, we probe the degree of coarse graining required to simultaneously retain significant atomistic details and access large length and time scales. The degree of coarse graining in turn sets the minimum length scale instrumental in defining polymer properties and dynamics. Using linear polyethylene as a model system, we probe how the coarse-graining scale affects the measured dynamics. Iterative Boltzmann inversion ismore » used to derive coarse-grained potentials with 2–6 methylene groups per coarse-grained bead from a fully atomistic melt simulation. We show that atomistic detail is critical to capturing large-scale dynamics. Finally, using these models we simulate polyethylene melts for times over 500 μs to study the viscoelastic properties of well-entangled polymer melts.« less

Ice age True Polar Wander: raising debates and new analyses

NASA Astrophysics Data System (ADS)

Sabadini, Roberto; Cambiotti, Gabriele; Ricard, Yanick

2010-05-01

Issues related to long time scale instability in the Earth's rotation, named True Polar Wander (TPW), have continuously been debated, after the pioneering works of the sixties. Since Maxwell Earth models with elastic or high viscosity viscoelastic lithospheres predict different ice-age TPW in the lower mantle viscosity range 1021 - 1022 Pa s, it has been recently suggested that the observed fluid Love number should be used to describe the initial equatorial bulge rather than the tidal fluid limit resulting from the viscoelastic modelling itself. We show that different ice-age TPW predictions have to be expected due to the dependence of TPW on the Earth's initial state, characterized by a larger and stress-free equatorial bulge for the viscoelastic lithosphere, compared to the elastic one, and that there is no shortcomings or errors in the traditional approach based on the use of tidal Love number from the model. The use of the observed fluid Love number represents in fact a simplified attempt to couple the effects on TPW from mantle convection and glacial forcing, by including the non-hydrostatic flattening due to mantle convection but not its driving part. This partial coupling freezes in space the non-hydrostatic contribution due to mantle convection, thus damping the present-day ice-age TPW and forcing the axis of instantaneous rotation to come back to its initial position when ice ages started. In this perspective, we discuss the implication of self-consistent convection calculations of the non-hydrostatic contribution and its impact on the long-term Earth's rotation stability during ice-age. We develop a full compressible model, based on the numerical integration in the radial variable of the momentum and Poisson equations and on the contour integration in the Laplace domain, which allows us to deal with the non-modal contribution from continuous radial rheological variations. We quantify the effects of the compressible rheology, compared to the widely used incompressible ones

Nanoscale effects in the characterization of viscoelastic materials with atomic force microscopy: Coupling of a quasi-three-dimensional standard linear solid model with in-plane surface interactions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Solares, Santiago D.

Significant progress has been accomplished in the development of experimental contact-mode and dynamic-mode atomic force microscopy (AFM) methods designed to measure surface material properties. However, current methods are based on one-dimensional (1D) descriptions of the tip-sample interaction forces, thus neglecting the intricacies involved in the material behavior of complex samples (such as soft viscoelastic materials) as well as the differences in material response between the surface and the bulk. In order to begin to address this gap, a computational study is presented where the sample is simulated using an enhanced version of a recently introduced model that treats the surfacemore » as a collection of standard-linear-solid viscoelastic elements. The enhanced model introduces in-plane surface elastic forces that can be approximately related to a two-dimensional (2D) Young's modulus. Relevant cases are discussed for single-and multifrequency intermittent-contact AFM imaging, with focus on the calculated surface indentation profiles and tip-sample interaction force curves, as well as their implications with regards to experimental interpretation. A variety of phenomena are examined in detail, which highlight the need for further development of more physically accurate sample models that are specifically designed for AFM simulation. As a result, a multifrequency AFM simulation tool based on the above sample model is provided as supporting information.« less

Nanoscale effects in the characterization of viscoelastic materials with atomic force microscopy: Coupling of a quasi-three-dimensional standard linear solid model with in-plane surface interactions

DOE PAGES

Solares, Santiago D.

2016-04-15

Significant progress has been accomplished in the development of experimental contact-mode and dynamic-mode atomic force microscopy (AFM) methods designed to measure surface material properties. However, current methods are based on one-dimensional (1D) descriptions of the tip-sample interaction forces, thus neglecting the intricacies involved in the material behavior of complex samples (such as soft viscoelastic materials) as well as the differences in material response between the surface and the bulk. In order to begin to address this gap, a computational study is presented where the sample is simulated using an enhanced version of a recently introduced model that treats the surfacemore » as a collection of standard-linear-solid viscoelastic elements. The enhanced model introduces in-plane surface elastic forces that can be approximately related to a two-dimensional (2D) Young's modulus. Relevant cases are discussed for single-and multifrequency intermittent-contact AFM imaging, with focus on the calculated surface indentation profiles and tip-sample interaction force curves, as well as their implications with regards to experimental interpretation. A variety of phenomena are examined in detail, which highlight the need for further development of more physically accurate sample models that are specifically designed for AFM simulation. As a result, a multifrequency AFM simulation tool based on the above sample model is provided as supporting information.« less

Active mechanical stabilization of the viscoplastic intracellular space of Dictyostelia cells by microtubule-actin crosstalk.

PubMed

Heinrich, Doris; Sackmann, Erich

2006-11-01

The micro-viscoelasticity of the intracellular space of Dictyostelium discoideum cells is studied by evaluating the intracellular transport of magnetic force probes and their viscoelastic responses to force pulses of 20-700 pN. The role of the actin cortex, the microtubule (MT) aster and their crosstalk is explored by comparing the behaviour of wild-type cells, myosin II null mutants, latrunculin A and benomyl treated cells. The MT coupled beads perform irregular local and long range directed motions which are characterized by measuring their velocity distributions (P(v)). The correlated motion of the MT and the centrosome are evaluated by microfluorescence of GFP-labelled MTs. P(v) can be represented by log-normal distributions with long tails and it is determined by random sweeping motions (v approximately 0.5 microm/s) of the MTs (caused by tangential forces on the filament ends coupled to the actin cortex) and by intermittent bead transports parallel to the MTs (v(max) approximately 1.5 microm/s). The tails are due to spontaneous filament deflections (with speeds up to 10 microm/s) attributed to pre-stressing of the MT by local cortical tensions, generated by dynactin motors generating plus-end directed forces in the MTs. The viscoelastic responses are strongly non-linear and are mostly directed opposite or perpendicular to the force, showing that the cytoplasm behaves as an active viscoplastic body with time and force dependent drag coefficients. Nano-Newton loads exerted on the soft MT are balanced by traction forces arising at the MT ends coupled to the actin cortex and the centrosome, respectively. The mechanical coupling between the soft microtubules and the viscoelastic actin cortex provides cells with high mechanical stability despite the softness of the cytoplasm.

New developments in isotropic turbulent models for FENE-P fluids

NASA Astrophysics Data System (ADS)

Resende, P. R.; Cavadas, A. S.

2018-04-01

The evolution of viscoelastic turbulent models, in the last years, has been significant due to the direct numeric simulation (DNS) advances, which allowed us to capture in detail the evolution of the viscoelastic effects and the development of viscoelastic closures. New viscoelastic closures are proposed for viscoelastic fluids described by the finitely extensible nonlinear elastic-Peterlin constitutive model. One of the viscoelastic closure developed in the context of isotropic turbulent models, consists in a modification of the turbulent viscosity to include an elastic effect, capable of predicting, with good accuracy, the behaviour for different drag reductions. Another viscoelastic closure essential to predict drag reduction relates the viscoelastic term involving velocity and the tensor conformation fluctuations. The DNS data show the high impact of this term to predict correctly the drag reduction, and for this reason is proposed a simpler closure capable of predicting the viscoelastic behaviour with good performance. In addition, a new relation is developed to predict the drag reduction, quantity based on the trace of the tensor conformation at the wall, eliminating the need of the typically parameters of Weissenberg and Reynolds numbers, which depend on the friction velocity. This allows future developments for complex geometries.

Electrostatics of the protein-water interface and the dynamical transition in proteins.

PubMed

Matyushov, Dmitry V; Morozov, Alexander Y

2011-07-01

Atomic displacements of hydrated proteins are dominated by phonon vibrations at low temperatures and by dissipative large-amplitude motions at high temperatures. A crossover between the two regimes is known as a dynamical transition. Recent experiments indicate a connection between the dynamical transition and the dielectric response of the hydrated protein. We analyze two mechanisms of the coupling between the protein atomic motions and the protein-water interface. The first mechanism considers viscoelastic changes in the global shape of the protein plasticized by its coupling to the hydration shell. The second mechanism involves modulations of the local motions of partial charges inside the protein by electrostatic fluctuations. The model is used to analyze mean-square displacements of iron of metmyoglobin reported by Mössbauer spectroscopy. We show that high displacement of heme iron at physiological temperatures is dominated by electrostatic fluctuations. Two onsets, one arising from the viscoelastic response and the second from electrostatic fluctuations, are seen in the temperature dependence of the mean-square displacements when the corresponding relaxation times enter the instrumental resolution window.

Electrostatics of the protein-water interface and the dynamical transition in proteins

NASA Astrophysics Data System (ADS)

Matyushov, Dmitry V.; Morozov, Alexander Y.

2011-07-01

Atomic displacements of hydrated proteins are dominated by phonon vibrations at low temperatures and by dissipative large-amplitude motions at high temperatures. A crossover between the two regimes is known as a dynamical transition. Recent experiments indicate a connection between the dynamical transition and the dielectric response of the hydrated protein. We analyze two mechanisms of the coupling between the protein atomic motions and the protein-water interface. The first mechanism considers viscoelastic changes in the global shape of the protein plasticized by its coupling to the hydration shell. The second mechanism involves modulations of the local motions of partial charges inside the protein by electrostatic fluctuations. The model is used to analyze mean-square displacements of iron of metmyoglobin reported by Mössbauer spectroscopy. We show that high displacement of heme iron at physiological temperatures is dominated by electrostatic fluctuations. Two onsets, one arising from the viscoelastic response and the second from electrostatic fluctuations, are seen in the temperature dependence of the mean-square displacements when the corresponding relaxation times enter the instrumental resolution window.

The Singularity Mystery Associated with a Radially Continuous Maxwell Viscoelastic Structure

NASA Technical Reports Server (NTRS)

Fang, Ming; Hager, Bradford H.

1995-01-01

The singularity problem associated with a radially continuous Maxwell viscoclastic structure is investigated. A special tool called the isolation function is developed. Results calculated using the isolation function show that the discrete model assumption is no longer valid when the viscoelastic parameter becomes a continuous function of radius. Continuous variations in the upper mantle viscoelastic parameter are especially powerful in destroying the mode-like structures. The contribution to the load Love numbers of the singularities is sensitive to the convexity of the viscoelastic parameter models. The difference between the vertical response and the horizontal response found in layered viscoelastic parameter models remains with continuous models.

Automated palpation for breast tissue discrimination based on viscoelastic biomechanical properties.

PubMed

Tsukune, Mariko; Kobayashi, Yo; Miyashita, Tomoyuki; Fujie, G Masakatsu

2015-05-01

Accurate, noninvasive methods are sought for breast tumor detection and diagnosis. In particular, a need for noninvasive techniques that measure both the nonlinear elastic and viscoelastic properties of breast tissue has been identified. For diagnostic purposes, it is important to select a nonlinear viscoelastic model with a small number of parameters that highly correlate with histological structure. However, the combination of conventional viscoelastic models with nonlinear elastic models requires a large number of parameters. A nonlinear viscoelastic model of breast tissue based on a simple equation with few parameters was developed and tested. The nonlinear viscoelastic properties of soft tissues in porcine breast were measured experimentally using fresh ex vivo samples. Robotic palpation was used for measurements employed in a finite element model. These measurements were used to calculate nonlinear viscoelastic parameters for fat, fibroglandular breast parenchyma and muscle. The ability of these parameters to distinguish the tissue types was evaluated in a two-step statistical analysis that included Holm's pairwise [Formula: see text] test. The discrimination error rate of a set of parameters was evaluated by the Mahalanobis distance. Ex vivo testing in porcine breast revealed significant differences in the nonlinear viscoelastic parameters among combinations of three tissue types. The discrimination error rate was low among all tested combinations of three tissue types. Although tissue discrimination was not achieved using only a single nonlinear viscoelastic parameter, a set of four nonlinear viscoelastic parameters were able to reliably and accurately discriminate fat, breast fibroglandular tissue and muscle.

Postseismic viscoelastic surface deformation and stress. Part 1: Theoretical considerations, displacement and strain calculations

NASA Technical Reports Server (NTRS)

Cohen, S. C.

1979-01-01

A model of viscoelastic deformations associated with earthquakes is presented. A strike-slip fault is represented by a rectangular dislocation in a viscoelastic layer (lithosphere) lying over a viscoelastic half-space (asthenosphere). Deformations occur on three time scales. The initial response is governed by the instantaneous elastic properties of the earth. A slower response is associated with viscoelastic relaxation of the lithosphere and a yet slower response is due to viscoelastic relaxation of the asthenosphere. The major conceptual contribution is the inclusion of lithospheric viscoelastic properties into a dislocation model of earthquake related deformations and stresses. Numerical calculations using typical fault parameters reveal that the postseismic displacements and strains are small compared to the coseismic ones near the fault, but become significant further away. Moreover, the directional sense of the deformations attributable to the elastic response, the lithospheric viscoelastic softening, and the asthenospheric viscoelastic flow may differ and depend on location and model details. The results and theoretical arguments suggest that the stress changes accompanying lithospheric relaxation may also be in a different sense than and be larger than the strain changes.

A Viscoelastic Constitutive Model Can Accurately Represent Entire Creep Indentation Tests of Human Patella Cartilage

PubMed Central

Pal, Saikat; Lindsey, Derek P.; Besier, Thor F.; Beaupre, Gary S.

2013-01-01

Cartilage material properties provide important insights into joint health, and cartilage material models are used in whole-joint finite element models. Although the biphasic model representing experimental creep indentation tests is commonly used to characterize cartilage, cartilage short-term response to loading is generally not characterized using the biphasic model. The purpose of this study was to determine the short-term and equilibrium material properties of human patella cartilage using a viscoelastic model representation of creep indentation tests. We performed 24 experimental creep indentation tests from 14 human patellar specimens ranging in age from 20 to 90 years (median age 61 years). We used a finite element model to reproduce the experimental tests and determined cartilage material properties from viscoelastic and biphasic representations of cartilage. The viscoelastic model consistently provided excellent representation of the short-term and equilibrium creep displacements. We determined initial elastic modulus, equilibrium elastic modulus, and equilibrium Poisson’s ratio using the viscoelastic model. The viscoelastic model can represent the short-term and equilibrium response of cartilage and may easily be implemented in whole-joint finite element models. PMID:23027200

The Sedimentation of Particles under Orthogonal Shear in Viscoelastic Fluids

NASA Astrophysics Data System (ADS)

Murch, William L.; Krishnan, Sreenath; Shaqfeh, Eric S. G.

2016-11-01

Many engineering applications, including oil and gas recovery, require the suspension of particles in viscoelastic fluids during fluid transport and processing. A topic of specific importance involves such particle suspensions experiencing an applied shear flow in a direction perpendicular to gravity (referred to as orthogonal shear). Previously, it has been shown that particle sedimentation coupled with an orthogonal shear flow can reduce the particle settling rate in elastic fluids. The underlying mechanism of this enhanced coupling drag is not fully understood, particularly at finite Weissenberg numbers. This talk examines the role of fluid elasticity on a single, non-Brownian, rigid sphere settling in orthogonal shear using experiments and numerical simulations. New experiments were performed in a Taylor-Couette flow cell using Boger fluids to study the coupling drag as a function of the shear and sedimentation Weissenberg numbers as well as particle confinement. The elastic effect was also studied with fully 3D simulations of flow past a rigid sphere, using the FENE-P constitutive model to describe the polymeric fluid rheology. These simulations show good agreement with the experiments and allow for further insight into the mechanism of elasticity-enhanced drag. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship.

An equivalent viscoelastic model for rock mass with parallel joints

NASA Astrophysics Data System (ADS)

Li, Jianchun; Ma, Guowei; Zhao, Jian

2010-03-01

An equivalent viscoelastic medium model is proposed for rock mass with parallel joints. A concept of "virtual wave source (VWS)" is proposed to take into account the wave reflections between the joints. The equivalent model can be effectively applied to analyze longitudinal wave propagation through discontinuous media with parallel joints. Parameters in the equivalent viscoelastic model are derived analytically based on longitudinal wave propagation across a single rock joint. The proposed model is then verified by applying identical incident waves to the discontinuous and equivalent viscoelastic media at one end to compare the output waves at the other end. When the wavelength of the incident wave is sufficiently long compared to the joint spacing, the effect of the VWS on wave propagation in rock mass is prominent. The results from the equivalent viscoelastic medium model are very similar to those determined from the displacement discontinuity method. Frequency dependence and joint spacing effect on the equivalent viscoelastic model and the VWS method are discussed.

The poro-viscoelastic properties of trabecular bone: a micro computed tomography-based finite element study.

PubMed

Sandino, Clara; McErlain, David D; Schipilow, John; Boyd, Steven K

2015-04-01

Bone is a porous structure with a solid phase that contains hydroxyapatite and collagen. Due to its composition, bone is often represented either as a poroelastic or as a viscoelastic material; however, the poro-viscoelastic formulation that allows integrating the effect of both the fluid flow and the collagen on the mechanical response of the tissue, has not been applied yet. The objective of this study was to develop a micro computed tomography (µCT)-based finite element (FE) model of trabecular bone that includes both the poroelastic and the viscoelastic nature of the tissue. Cubes of trabecular bone (N=25) from human distal tibia were scanned with µCT and stress relaxation experiments were conducted. The µCT images were the basis for sample specific FE models, and the stress relaxation experiments were simulated applying a poro-viscoelastic formulation. The model considers two scales of the tissue: the intertrabecular pore and the lacunar-canalicular pore scales. Independent viscoelastic and poroelastic models were also developed to determine their contribution to the poro-viscoelastic model. All the experiments exhibited a similar relaxation trend. The average reaction force before relaxation was 9.28 × 10(2)N (SD ± 5.11 × 10(2)N), and after relaxation was 4.69 × 10(2)N (SD ± 2.88 × 10(2)N). The slope of the regression line between the force before and after relaxation was 1.92 (R(2)=0.96). The poro-viscoelastic models captured 49% of the variability of the experimental data before relaxation and 33% after relaxation. The relaxation predicted with viscoelastic models was similar to the poro-viscoelastic ones; however, the poroelastic formulation underestimated the reaction force before relaxation. These data suggest that the contribution of viscoelasticity (fluid flow-independent mechanism) to the mechanical response of the tissue is significantly greater than the contribution of the poroelasticity (fluid flow-dependent mechanism). Copyright © 2015 Elsevier Ltd. All rights reserved.

Micromechanics of transformation fields in ageing linear viscoelastic composites: effects of phase dissolution or precipitation

NASA Astrophysics Data System (ADS)

Honorio, Tulio

2017-11-01

Transformation fields, in an affine formulation characterizing mechanical behavior, describe a variety of physical phenomena regardless their origin. Different composites, notably geomaterials, present a viscoelastic behavior, which is, in some cases of industrial interest, ageing, i.e. it evolves independently with respect to time and loading time. Here, a general formulation of the micromechanics of prestressed or prestrained composites in Ageing Linear Viscoelasticity (ALV) is presented. Emphasis is put on the estimation of effective transformation fields in ALV. The result generalizes Ageing Linear Thermo- and Poro-Viscoelasticity and it can be used in approaches coping with a phase transformation. Additionally, the results are extended to the case of locally transforming materials due to non-coupled dissolution and/or precipitation of a given (elastic or viscoelastic) phase. The estimations of locally transforming composites can be made with respect to different morphologies. As an application, estimations of the coefficient of thermal expansion of a hydrating alite paste are presented.

Shock-induced thermal wave propagation and response analysis of a viscoelastic thin plate under transient heating loads

NASA Astrophysics Data System (ADS)

Li, Chenlin; Guo, Huili; Tian, Xiaogeng

2018-04-01

This paper is devoted to the thermal shock analysis for viscoelastic materials under transient heating loads. The governing coupled equations with time-delay parameter and nonlocal scale parameter are derived based on the generalized thermo-viscoelasticity theory. The problem of a thin plate composed of viscoelastic material, subjected to a sudden temperature rise at the boundary plane, is solved by employing Laplace transformation techniques. The transient responses, i.e. temperature, displacement, stresses, heat flux as well as strain, are obtained and discussed. The effects of time-delay and nonlocal scale parameter on the transient responses are analyzed and discussed. It can be observed that: the propagation of thermal wave is dynamically smoothed and changed with the variation of time-delay; while the displacement, strain, and stress can be rapidly reduced by nonlocal scale parameter, which can be viewed as an important indicator for predicting the stiffness softening behavior for viscoelastic materials.

Enabling large-scale viscoelastic calculations via neural network acceleration

NASA Astrophysics Data System (ADS)

Robinson DeVries, P.; Thompson, T. B.; Meade, B. J.

2017-12-01

One of the most significant challenges involved in efforts to understand the effects of repeated earthquake cycle activity are the computational costs of large-scale viscoelastic earthquake cycle models. Deep artificial neural networks (ANNs) can be used to discover new, compact, and accurate computational representations of viscoelastic physics. Once found, these efficient ANN representations may replace computationally intensive viscoelastic codes and accelerate large-scale viscoelastic calculations by more than 50,000%. This magnitude of acceleration enables the modeling of geometrically complex faults over thousands of earthquake cycles across wider ranges of model parameters and at larger spatial and temporal scales than have been previously possible. Perhaps most interestingly from a scientific perspective, ANN representations of viscoelastic physics may lead to basic advances in the understanding of the underlying model phenomenology. We demonstrate the potential of artificial neural networks to illuminate fundamental physical insights with specific examples.

Viscoelastic Properties of Human Tracheal Tissues.

PubMed

Safshekan, Farzaneh; Tafazzoli-Shadpour, Mohammad; Abdouss, Majid; Shadmehr, Mohammad B

2017-01-01

The physiological performance of trachea is highly dependent on its mechanical behavior, and therefore, the mechanical properties of its components. Mechanical characterization of trachea is key to succeed in new treatments such as tissue engineering, which requires the utilization of scaffolds which are mechanically compatible with the native human trachea. In this study, after isolating human trachea samples from brain-dead cases and proper storage, we assessed the viscoelastic properties of tracheal cartilage, smooth muscle, and connective tissue based on stress relaxation tests (at 5% and 10% strains for cartilage and 20%, 30%, and 40% for smooth muscle and connective tissue). After investigation of viscoelastic linearity, constitutive models including Prony series for linear viscoelasticity and quasi-linear viscoelastic, modified superposition, and Schapery models for nonlinear viscoelasticity were fitted to the experimental data to find the best model for each tissue. We also investigated the effect of age on the viscoelastic behavior of tracheal tissues. Based on the results, all three tissues exhibited a (nonsignificant) decrease in relaxation rate with increasing the strain, indicating viscoelastic nonlinearity which was most evident for cartilage and with the least effect for connective tissue. The three-term Prony model was selected for describing the linear viscoelasticity. Among different models, the modified superposition model was best able to capture the relaxation behavior of the three tracheal components. We observed a general (but not significant) stiffening of tracheal cartilage and connective tissue with aging. No change in the stress relaxation percentage with aging was observed. The results of this study may be useful in the design and fabrication of tracheal tissue engineering scaffolds.

A finite nonlinear hyper-viscoelastic model for soft biological tissues.

PubMed

Panda, Satish Kumar; Buist, Martin Lindsay

2018-03-01

Soft tissues exhibit highly nonlinear rate and time-dependent stress-strain behaviour. Strain and strain rate dependencies are often modelled using a hyperelastic model and a discrete (standard linear solid) or continuous spectrum (quasi-linear) viscoelastic model, respectively. However, these models are unable to properly capture the materials characteristics because hyperelastic models are unsuited for time-dependent events, whereas the common viscoelastic models are insufficient for the nonlinear and finite strain viscoelastic tissue responses. The convolution integral based models can demonstrate a finite viscoelastic response; however, their derivations are not consistent with the laws of thermodynamics. The aim of this work was to develop a three-dimensional finite hyper-viscoelastic model for soft tissues using a thermodynamically consistent approach. In addition, a nonlinear function, dependent on strain and strain rate, was adopted to capture the nonlinear variation of viscosity during a loading process. To demonstrate the efficacy and versatility of this approach, the model was used to recreate the experimental results performed on different types of soft tissues. In all the cases, the simulation results were well matched (R 2 ⩾0.99) with the experimental data. Copyright © 2018 Elsevier Ltd. All rights reserved.

Analyzing the Deposition of Titanium Dioxide Nanoparticles at Model Rough Mineral Surfaces Using a Quartz Crystal Microbalance with Dissipation Monitoring

NASA Astrophysics Data System (ADS)

Li, Y.; Kananizadeh, N.; Rodenhausen, K. B.; Schubert, M.; Bartelt-Hunt, S.

2015-12-01

Titanium dioxide nanoparticles (nTiO2) is the most extensively manufactured engineered materials. nTiO2 from sunscreens was found to enter sediments after released into a lake. nTiO2 may also enter the subsurface via irrigation using effluents from wastewater treatment plants. Interaction of nTiO2 with soils and sediments will largely influence their fate, transport, and ecotoxicity. Measuring the interaction between nTiO2 and natural substrates (e.g. such as sands) is particularly challenging due to highly heterogeneous and rough natural sand surfaces. In this study, an engineered controllable rough surface known as three dimensional nanostructured sculptured columnar thin films (SCTFs) has been used to mimic surface roughness. SCTFs were fabricated by glancing angle deposition (GLAD), a physical vapor deposition technique facilitated by electron beam evaporation. Interaction between nTiO2 and SCTF coated surfaces was investigated using a quartz crystal microbalance with dissipation monitoring (QCM-D). In parallel, a Generalized Ellipsometry (GE) was coupled with the QCM-D to measure the deposition of nTiO2. We found that the typical QCM-D modeling approach, e.g. viscoelastic model, would largely overestimate the mass of deposited nTiO2, because the frequency drops due to particle deposition or water entrapment in rough areas were not differentiated. Here, we demonstrate a new approach to model QCM-D data for nTiO2 deposition on rough surfaces, which couples the viscoelastic model with a model of flow on the non-uniform surface.

Estimating material viscoelastic properties based on surface wave measurements: A comparison of techniques and modeling assumptions

PubMed Central

Royston, Thomas J.; Dai, Zoujun; Chaunsali, Rajesh; Liu, Yifei; Peng, Ying; Magin, Richard L.

2011-01-01

Previous studies of the first author and others have focused on low audible frequency (<1 kHz) shear and surface wave motion in and on a viscoelastic material comprised of or representative of soft biological tissue. A specific case considered has been surface (Rayleigh) wave motion caused by a circular disk located on the surface and oscillating normal to it. Different approaches to identifying the type and coefficients of a viscoelastic model of the material based on these measurements have been proposed. One approach has been to optimize coefficients in an assumed viscoelastic model type to match measurements of the frequency-dependent Rayleigh wave speed. Another approach has been to optimize coefficients in an assumed viscoelastic model type to match the complex-valued frequency response function (FRF) between the excitation location and points at known radial distances from it. In the present article, the relative merits of these approaches are explored theoretically, computationally, and experimentally. It is concluded that matching the complex-valued FRF may provide a better estimate of the viscoelastic model type and parameter values; though, as the studies herein show, there are inherent limitations to identifying viscoelastic properties based on surface wave measurements. PMID:22225067

Non-integer viscoelastic constitutive law to model soft biological tissues to in-vivo indentation.

PubMed

Demirci, Nagehan; Tönük, Ergin

2014-01-01

During the last decades, derivatives and integrals of non-integer orders are being more commonly used for the description of constitutive behavior of various viscoelastic materials including soft biological tissues. Compared to integer order constitutive relations, non-integer order viscoelastic material models of soft biological tissues are capable of capturing a wider range of viscoelastic behavior obtained from experiments. Although integer order models may yield comparably accurate results, non-integer order material models have less number of parameters to be identified in addition to description of an intermediate material that can monotonically and continuously be adjusted in between an ideal elastic solid and an ideal viscous fluid. In this work, starting with some preliminaries on non-integer (fractional) calculus, the "spring-pot", (intermediate mechanical element between a solid and a fluid), non-integer order three element (Zener) solid model, finally a user-defined large strain non-integer order viscoelastic constitutive model was constructed to be used in finite element simulations. Using the constitutive equation developed, by utilizing inverse finite element method and in vivo indentation experiments, soft tissue material identification was performed. The results indicate that material coefficients obtained from relaxation experiments, when optimized with creep experimental data could simulate relaxation, creep and cyclic loading and unloading experiments accurately. Non-integer calculus viscoelastic constitutive models, having physical interpretation and modeling experimental data accurately is a good alternative to classical phenomenological viscoelastic constitutive equations.

Tension-compression viscoelastic behaviors of the periodontal ligament.

PubMed

Wang, Chen-Ying; Su, Ming-Zen; Chang, Hao-Hueng; Chiang, Yu-Chih; Tao, Shao-Huan; Cheng, Jung-Ho; Fuh, Lih-Jyh; Lin, Chun-Pin

2012-09-01

Although exhaustively studied, the mechanism responsible for tooth support and the mechanical properties of the periodontal ligament (PDL) remain a subject of considerable controversy. In the past, various experimental techniques and theoretical analyses have been employed to tackle this intricate problem. The aim of this study was to investigate the viscoelastic behaviors of the PDL using three-dimensional finite element analysis. Three dentoalveolar complex models were established to simulate the tissue behaviors of the PDL: (1) deviatoric viscoelastic model; (2) volumetric viscoelastic model; and (3) tension-compression volumetric viscoelastic model. These modified models took into consideration the presence of tension and compression along the PDL during both loading and unloading. The inverse parameter identification process was developed to determine the mechanical properties of the PDL from the results of previously reported in vitro and in vivo experiments. The results suggest that the tension-compression volumetric viscoelastic model is a good approximation of normal PDL behavior during the loading-unloading process, and the deviatoric viscoelastic model is a good representation of how a damaged PDL behaves under loading conditions. Moreover, fluid appears to be the main creep source in the PDL. We believe that the biomechanical properties of the PDL established via retrograde calculation in this study can lead to the construction of more accurate extra-oral models and a comprehensive understanding of the biomechanical behavior of the PDL. Copyright © 2012. Published by Elsevier B.V.

Damped response of shells by a constrained viscoelastic layer

NASA Technical Reports Server (NTRS)

El-Raheb, M.; Wagner, P.

1986-01-01

Vibration absorbers are introduced into an asymmetric configuration of thin cylinders and tori enclosing an acoustic medium. The absorbers consist of thin axial strips bonded to the cylinder with a thin viscoelastic layer. The constrained layer dissipates the energy of relative motions between strip and cylinder. The absorber is most effective on response modes with two or more circumferential waves. The use of transfer matrices is extended to the coupled cylinder-absorber system.

A nonlinear and fractional derivative viscoelastic model for rail pads in the dynamic analysis of coupled vehicle-slab track systems

NASA Astrophysics Data System (ADS)

Zhu, Shengyang; Cai, Chengbiao; Spanos, Pol D.

2015-01-01

A nonlinear and fractional derivative viscoelastic (FDV) model is used to capture the complex behavior of rail pads. It is implemented into the dynamic analysis of coupled vehicle-slab track (CVST) systems. The vehicle is treated as a multi-body system with 10 degrees of freedom, and the slab track is represented by a three layer Bernoulli-Euler beam model. The model for the rail pads is one dimensional, and the force-displacement relation is based on a superposition of elastic, friction, and FDV forces. This model takes into account the influences of the excitation frequency and of the displacement amplitude through a fractional derivative element, and a nonlinear friction element, respectively. The Grünwald representation of the fractional derivatives is employed to numerically solve the fractional and nonlinear equations of motion of the CVST system by means of an explicit integration algorithm. A dynamic analysis of the CVST system exposed to excitations of rail harmonic irregularities is carried out, pointing out the stiffness and damping dependence on the excitation frequency and the displacement amplitude. The analysis indicates that the dynamic stiffness and damping of the rail pads increase with the excitation frequency while they decrease with the displacement amplitude. Furthermore, comparisons between the proposed model and ordinary Kelvin model adopted for the CVST system, under excitations of welded rail joint irregularities and of random track irregularities, are conducted in the time domain as well as in the frequency domain. The proposed model is shown to possess several modeling advantages over the ordinary Kelvin element which overestimates both the stiffness and damping features at high frequencies.

The relationship between the instantaneous velocity field and the rate of moment release in the lithosphere

USGS Publications Warehouse

Pollitz, F.F.

2003-01-01

Instantaneous velocity gradients within the continental lithosphere are often related to the tectonic driving forces. This relationship is direct if the forces are secular, as for the case of loading of a locked section of a subduction interface by the downgoing plate. If the forces are static, as for the case of lateral variations in gravitational potential energy, then velocity gradients can be produced only if the lithosphere has, on average, zero strength. The static force model may be related to the long-term velocity field but not the instantaneous velocity field (typically measured geodetically over a period of several years) because over short time intervals the upper lithosphere behaves elastically. In order to describe both the short- and long-term behaviour of an (elastic) lithosphere-(viscoelastic) asthenosphere system in a self-consistent manner, I construct a deformation model termed the expected interseismic velocity (EIV) model. Assuming that the lithosphere is populated with faults that rupture continually, each with a definite mean recurrence time, and that the Earth is well approximated as a linear elastic-viscoelastic coupled system, I derive a simple relationship between the instantaneous velocity field and the average rate of moment release in the lithosphere. Examples with synthetic fault networks demonstrate that velocity gradients in actively deforming regions may to a large extent be the product of compounded viscoelastic relaxation from past earthquakes on hundreds of faults distributed over large ( ≥106 km2) areas.

A K-BKZ Formulation for Soft-Tissue Viscoelasticity

NASA Technical Reports Server (NTRS)

Freed, Alan D.; Diethelm, Kai

2005-01-01

A viscoelastic model of the K-BKZ (Kaye 1962; Bernstein et al. 1963) type is developed for isotropic biological tissues, and applied to the fat pad of the human heel. To facilitate this pursuit, a class of elastic solids is introduced through a novel strain-energy function whose elements possess strong ellipticity, and therefore lead to stable material models. The standard fractional-order viscoelastic (FOV) solid is used to arrive at the overall elastic/viscoelastic structure of the model, while the elastic potential via the K-BKZ hypothesis is used to arrive at the tensorial structure of the model. Candidate sets of functions are proposed for the elastic and viscoelastic material functions present in the model, including a regularized fractional derivative that was determined to be the best. The Akaike information criterion (AIC) is advocated for performing multi-model inference, enabling an objective selection of the best material function from within a candidate set.

Water Touch-and-Bounce from a Soft Viscoelastic Substrate: Wetting, Dewetting, and Rebound on Bitumen.

PubMed

Lee, Jae Bong; Dos Santos, Salomé; Antonini, Carlo

2016-08-16

Understanding the interaction between liquids and deformable solid surfaces is a fascinating fundamental problem, in which interaction and coupling of capillary and viscoelastic effects, due to solid substrate deformation, give rise to complex wetting mechanisms. Here we investigated as a model case the behavior of water drops on two smooth bitumen substrates with different rheological properties, defined as hard and soft (with complex shear moduli in the order of 10(7) and 10(5) Pa, respectively, at 1 Hz), focusing both on wetting and on dewetting behavior. By means of classical quasi-static contact angle measurements and drop impact tests, we show that the water drop behavior can significantly change from the quasi-static to the dynamic regime on soft viscoelastic surfaces, with the transition being defined by the substrate rheological properties. As a result, we also show that on the hard substrate, where the elastic response is dominant under all investigated conditions, classical quasi-static contact angle measurements provide consistent results that can be used to predict the drop dynamic wetting behavior, such as drop deposition or rebound after impact, as typically observed for nondeformable substrates. Differently, on soft surfaces, the formation of wetting ridges did not allow to define uniquely the substrate intrinsic advancing and receding contact angles. In addition, despite showing a high adhesion to the soft surface in quasi-static measurements, the drop was surprisingly able to rebound and escape from the surface after impact, as it is typically observed for hydrophobic surfaces. These results highlight that measurements of wetting properties for viscoelastic substrates need to be critically used and that wetting behavior of a liquid on viscoelastic surfaces is a function of the characteristic time scales.

A viscoelastic-stochastic model of the effects of cytoskeleton remodelling on cell adhesion.

PubMed

Li, Long; Zhang, Wenyan; Wang, Jizeng

2016-10-01

Cells can adapt their mechanical properties through cytoskeleton remodelling in response to external stimuli when the cells adhere to the extracellular matrix (ECM). Many studies have investigated the effects of cell and ECM elasticity on cell adhesion. However, experiments determined that cells are viscoelastic and exhibiting stress relaxation, and the mechanism behind the effect of cellular viscoelasticity on the cell adhesion behaviour remains unclear. Therefore, we propose a theoretical model of a cluster of ligand-receptor bonds between two dissimilar viscoelastic media subjected to an applied tensile load. In this model, the distribution of interfacial traction is assumed to follow classical continuum viscoelastic equations, whereas the rupture and rebinding of individual molecular bonds are governed by stochastic equations. On the basis of this model, we determined that viscosity can significantly increase the lifetime, stability and dynamic strength of the adhesion cluster of molecular bonds, because deformation relaxation attributed to the viscoelastic property can increase the rebinding probability of each open bond and reduce the stress concentration in the adhesion area.

Two-dimensional analysis of post-seismic deformation of the 2011 Tohoku-Oki earthquake with rate-and-state friction and non-linear rock rheology

NASA Astrophysics Data System (ADS)

Muto, J.; Moore, J. D. P.; Barbot, S.; Iinuma, T.; Ohta, Y.; Horiuchi, S.; Hikaru, I.

2017-12-01

We conduct a two-dimensional (2D) analysis of the post-seismic deformation of the 2011 Tohoku-Oki earthquake with the nonlinear coupling between frictional afterslip and viscoelastic flow. We consider slip on the plate boundary and distributed viscous flow of the lower crust and mantle. We created 2D transects across the Miyagi-Yamagata area where the largest coseismic slip was observed. We use the stress change by the coseismic slip model of Iinuma et al. (2012) to drive the post-seismic relaxation. The simulation is performed by the integral method (Lambert & Barbot, 2016) expanded to plane strain (Barbot, Moore, & Lambert, 2017). Despite the simple 2D approximation, we look for a realistic model compatible with mineral physics to explain geodetic observations including 5 years of seafloor observations (Tomita et al., 2017). In the ductile regions, the model employs a bi-viscous Burgers rheology with power-law flow (Masuti et al., 2016). The steady-state viscosity is estimated based on a thermal structure obtained by thermal-flow model including the wedge corner flow (Horiuchi & Iwamori, 2016). We model afterslip by the regularized rate-strengthening approximation of the rate-and-state dependent friction law (Barbot et al., 2009). The combination of power-law rheology with stress-driven afterslip explains the observed 2D displacement fields well during the 5-year post-seismic period. We also find that the model requires a low viscosity ( 1018 Pas) body beneath the quaternary volcano (Mt. Naruko) to reproduce the localized subsidence detected in the 9-month post-seismic period (Muto et al., 2016). The introduction of the low-viscosity body also reproduces quick recovery of the subsidence in the 5-year period. Equipped with a reference model that fits available geodetic observations, we discuss the importance of the mechanical coupling between afterslip and viscoelastic flow. We find that ignoring the traction change on the fault by viscoelastic flow introduces variations of the order of 20% on the amplitude of afterslip. This effect is most pronounced late in the post-seismic relaxation. Our model reconciles laboratory constraints on rock rheology and geophysical observations after the earthquake and serves as a first-order reference to better understand the dynamics of subduction at the Japan trench.

Size-dependent geometrically nonlinear free vibration analysis of fractional viscoelastic nanobeams based on the nonlocal elasticity theory

NASA Astrophysics Data System (ADS)

Ansari, R.; Faraji Oskouie, M.; Gholami, R.

2016-01-01

In recent decades, mathematical modeling and engineering applications of fractional-order calculus have been extensively utilized to provide efficient simulation tools in the field of solid mechanics. In this paper, a nonlinear fractional nonlocal Euler-Bernoulli beam model is established using the concept of fractional derivative and nonlocal elasticity theory to investigate the size-dependent geometrically nonlinear free vibration of fractional viscoelastic nanobeams. The non-classical fractional integro-differential Euler-Bernoulli beam model contains the nonlocal parameter, viscoelasticity coefficient and order of the fractional derivative to interpret the size effect, viscoelastic material and fractional behavior in the nanoscale fractional viscoelastic structures, respectively. In the solution procedure, the Galerkin method is employed to reduce the fractional integro-partial differential governing equation to a fractional ordinary differential equation in the time domain. Afterwards, the predictor-corrector method is used to solve the nonlinear fractional time-dependent equation. Finally, the influences of nonlocal parameter, order of fractional derivative and viscoelasticity coefficient on the nonlinear time response of fractional viscoelastic nanobeams are discussed in detail. Moreover, comparisons are made between the time responses of linear and nonlinear models.

Energy dissipation in quasi-linear viscoelastic tissues, cells, and extracellular matrix.

PubMed

Babaei, Behzad; Velasquez-Mao, A J; Pryse, Kenneth M; McConnaughey, William B; Elson, Elliot L; Genin, Guy M

2018-05-21

Characterizing how a tissue's constituents give rise to its viscoelasticity is important for uncovering how hidden timescales underlie multiscale biomechanics. These constituents are viscoelastic in nature, and their mechanics must typically be assessed from the uniaxial behavior of a tissue. Confounding the challenge is that tissue viscoelasticity is typically associated with nonlinear elastic responses. Here, we experimentally assessed how fibroblasts and extracellular matrix (ECM) within engineered tissue constructs give rise to the nonlinear viscoelastic responses of a tissue. We applied a constant strain rate, "triangular-wave" loading and interpreted responses using the Fung quasi-linear viscoelastic (QLV) material model. Although the Fung QLV model has several well-known weaknesses, it was well suited to the behaviors of the tissue constructs, cells, and ECM tested. Cells showed relatively high damping over certain loading frequency ranges. Analysis revealed that, even in cases where the Fung QLV model provided an excellent fit to data, the the time constant derived from the model was not in general a material parameter. Results have implications for design of protocols for the mechanical characterization of biological materials, and for the mechanobiology of cells within viscoelastic tissues. Copyright © 2018. Published by Elsevier Ltd.

Crustal deformation, the earthquake cycle, and models of viscoelastic flow in the asthenosphere

NASA Technical Reports Server (NTRS)

Cohen, S. C.; Kramer, M. J.

1983-01-01

The crustal deformation patterns associated with the earthquake cycle can depend strongly on the rheological properties of subcrustal material. Substantial deviations from the simple patterns for a uniformly elastic earth are expected when viscoelastic flow of subcrustal material is considered. The detailed description of the deformation pattern and in particular the surface displacements, displacement rates, strains, and strain rates depend on the structure and geometry of the material near the seismogenic zone. The origin of some of these differences are resolved by analyzing several different linear viscoelastic models with a common finite element computational technique. The models involve strike-slip faulting and include a thin channel asthenosphere model, a model with a varying thickness lithosphere, and a model with a viscoelastic inclusion below the brittle slip plane. The calculations reveal that the surface deformation pattern is most sensitive to the rheology of the material that lies below the slip plane in a volume whose extent is a few times the fault depth. If this material is viscoelastic, the surface deformation pattern resembles that of an elastic layer lying over a viscoelastic half-space. When the thickness or breath of the viscoelastic material is less than a few times the fault depth, then the surface deformation pattern is altered and geodetic measurements are potentially useful for studying the details of subsurface geometry and structure. Distinguishing among the various models is best accomplished by making geodetic measurements not only near the fault but out to distances equal to several times the fault depth. This is where the model differences are greatest; these differences will be most readily detected shortly after an earthquake when viscoelastic effects are most pronounced.

A new insight into the dependence of relaxation time on frequency in viscoelastic surfactant solutions: From experimental to modeling study.

PubMed

García, Brayan F; Saraji, Soheil

2018-05-01

The relaxation time in viscoelastic surfactant solutions is a function of temperature, salt/surfactant concentrations, resting conditions, as well as shear frequency. The simplistic assumption of a single and constant relaxation time is not representative of all relaxation modes in these solutions especially at high frequencies. Steady-state and oscillatory measurements are carried out to study the effects of high temperature, concentration and resting condition on the rheology of surfactants/salt mixtures including a non-ionic and a zwitterionic/anionic surfactant system. Furthermore, a novel semi-empirical rheological model is deducted based on Cates theory.This model introduces, for the first time, a frequency-dependence for the continuous relaxation time spectrum. At high temperatures, the non-ionic surfactant become more viscoelastic and the zwitterionic/anionic system loses its viscoelasticity. The addition of surfactant/salt improves the viscoelasticity of both systems, and, for the zwitterionic/anionic mixture, increasing the resting temperature improves its viscoelasticity. In addition, the proposed model significantly improves predictions of traditional Maxwell model for different viscoelastic surfactant solutions (using data from this study and the literature) for a considerable range of surfactant and salt combinations at a wide range of temperature. Copyright © 2018 Elsevier Inc. All rights reserved.

Viscoelastic relaxations of high alcohols and alkanes: Effects of heterogeneous structure and translation-orientation coupling

NASA Astrophysics Data System (ADS)

Yamaguchi, Tsuyoshi

2017-03-01

The frequency-dependent shear viscosity of high alcohols and linear alkanes, including 1-butanol, 1-octanol, 1-dodecanol, n-hexane, n-decane, and n-tetradecane, was calculated using molecular dynamics simulation. The relaxation of all the liquids was bimodal. The correlation functions of the collective orientation were also evaluated. The analysis of these functions showed that the slower relaxation mode of alkanes is assigned to the translation-orientation coupling, while that of high alcohols is not. The X-ray structure factors of all the alcohols showed prepeaks, as have been reported in the literature, and the intermediate scattering functions were calculated at the prepeak. Comparing the intermediate scattering function with the frequency-dependent shear viscosity based on the mode-coupling theory, it was demonstrated that the slower viscoelastic relaxation of the alcohols is assigned to the relaxation of the heterogeneous structure described by the prepeak.

Geodetic estimates of fault slip rates in the San Francisco Bay area

USGS Publications Warehouse

Savage, J.C.; Svarc, J.L.; Prescott, W.H.

1999-01-01

Bourne et al. [1998] have suggested that the interseismic velocity profile at the surface across a transform plate boundary is a replica of the secular velocity profile at depth in the plastosphere. On the other hand, in the viscoelastic coupling model the shape of the interseismic surface velocity profile is a consequence of plastosphere relaxation following the previous rupture of the faults that make up the plate boundary and is not directly related to the secular flow in the plastosphere. The two models appear to be incompatible. If the plate boundary is composed of several subparallel faults and the interseismic surface velocity profile across the boundary known, each model predicts the secular slip rates on the faults which make up the boundary. As suggested by Bourne et al., the models can then be tested by comparing the predicted secular slip rates to those estimated from long-term offsets inferred from geology. Here we apply that test to the secular slip rates predicted for the principal faults (San Andreas, San Gregorio, Hayward, Calaveras, Rodgers Creek, Green Valley and Greenville faults) in the San Andreas fault system in the San Francisco Bay area. The estimates from the two models generally agree with one another and to a lesser extent with the geologic estimate. Because the viscoelastic coupling model has been equally successful in estimating secular slip rates on the various fault strands at a diffuse plate boundary, the success of the model of Bourne et al. [1998] in doing the same thing should not be taken as proof that the interseismic velocity profile across the plate boundary at the surface is a replica of the velocity profile at depth in the plastosphere.

A finite deformation viscoelastic-viscoplastic constitutive model for self-healing materials

NASA Astrophysics Data System (ADS)

Shahsavari, H.; Naghdabadi, R.; Baghani, M.; Sohrabpour, S.

2016-12-01

In this paper, employing the Hencky strain, viscoelastic-viscoplastic response of self-healing materials is investigated. Considering the irreversible thermodynamics and using the effective configuration in the Continuum Damage-Healing Mechanics (CDHM), a phenomenological finite strain viscoelastic-viscoplastic constitutive model is presented. Considering finite viscoelastic and viscoplastic deformations, total deformation gradient is multiplicatively decomposed into viscoelastic and viscoplastic parts. Due to mathematical advantages and physical meaning of Hencky strain, this measure of strain is employed in the constitutive model development. In this regard, defining the damage and healing variables and employing the strain equivalence hypothesis, the strain tensor is determined in the effective configuration. Satisfying the Clausius-Duhem inequality, the evolution equations are introduced for the viscoelastic and viscoplastic strains. The damage and healing variables also evolve according to two different prescribed functions. To employ the proposed model in different loading conditions, the model is discretized in the semi-implicit form. Material parameters of the model are identified employing experimental tests on asphalt mixes available in the literature. Finally, capability of the model is demonstrated comparing the model predictions in the creep-recovery and repeated creep-recovery with the experimental results available in the literature and a good agreement between predicted and test results is revealed.

Investigation of transient cavitating flow in viscoelastic pipes

NASA Astrophysics Data System (ADS)

Keramat, A.; Tijsseling, A. S.; Ahmadi, A.

2010-08-01

A study on water hammer in viscoelastic pipes when the fluid pressure drops to liquid vapour pressure is performed. Two important concepts including column separation and the effects of retarded strains in the pipe wall on the fluid response have been investigated separately in recent works, but there is some curiosity as to how the results for pressure and discharge are when column separation occurs in a viscoelastic pipe. For pipes made of plastic such as polyethylene (PE) and polyvinyl chloride (PVC), viscoelasticity is a crucial mechanical property which changes the hydraulic and structural transient responses. Based on previous developments in the analysis of water hammer, a model which is capable of analysing column separation in viscoelastic pipes is presented and used for solving the selected case studies. For the column-separation modelling the Discrete Vapour Cavity Model (DVCM) is utilised and the viscoelasticity property of the pipe wall is modelled by Kelvin-Voigt elements. The effects of viscoelasticity play an important role in the column separation phenomenon because it changes the water hammer fundamental frequency and so affects the time of opening or collapse of the cavities. Verification of the implemented computer code is performed for the effects of viscoelasticity and column separation - separately and simultaneously - using experimental results from the literature. In the provided examples the focus is placed on the simultaneous effect of viscoelasticity and column separation on the hydraulic transient response. The final conclusions drawn are that if rectangular grids are utilised the DVCM gives acceptable predictions of the phenomenon and that the pipe wall material's retarded behaviour strongly dampens the pressure spikes caused by column separation.

Mechanical properties of multifunctional structure with viscoelastic components based on FVE model

NASA Astrophysics Data System (ADS)

Hao, Dong; Zhang, Lin; Yu, Jing; Mao, Daiyong

2018-02-01

Based on the models of Lion and Kardelky (2004) and Hofer and Lion (2009), a finite viscoelastic (FVE) constitutive model, considering the predeformation-, frequency- and amplitude-dependent properties, has been proposed in our earlier paper [1]. FVE model is applied to investigating the dynamic characteristics of the multifunctional structure with the viscoelastic components. Combing FVE model with the finite element theory, the dynamic model of the multifunctional structure could be obtained. Additionally, the parametric identification and the experimental verification are also given via the frequency-sweep tests. The results show that the computational data agree well with the experimental data. FVE model has made a success of expressing the dynamic characteristics of the viscoelastic materials utilized in the multifunctional structure. The multifunctional structure technology has been verified by in-orbit experiments.

Coefficient of restitution in fractional viscoelastic compliant impacts using fractional Chebyshev collocation

NASA Astrophysics Data System (ADS)

Dabiri, Arman; Butcher, Eric A.; Nazari, Morad

2017-02-01

Compliant impacts can be modeled using linear viscoelastic constitutive models. While such impact models for realistic viscoelastic materials using integer order derivatives of force and displacement usually require a large number of parameters, compliant impact models obtained using fractional calculus, however, can be advantageous since such models use fewer parameters and successfully capture the hereditary property. In this paper, we introduce the fractional Chebyshev collocation (FCC) method as an approximation tool for numerical simulation of several linear fractional viscoelastic compliant impact models in which the overall coefficient of restitution for the impact is studied as a function of the fractional model parameters for the first time. Other relevant impact characteristics such as hysteresis curves, impact force gradient, penetration and separation depths are also studied.

Modelling water hammer in viscoelastic pipelines: short brief

NASA Astrophysics Data System (ADS)

Urbanowicz, K.; Firkowski, M.; Zarzycki, Z.

2016-10-01

The model of water hammer in viscoelastic pipelines is analyzed. An appropriate mathematical model of water hammer in polymer pipelines is presented. An additional term has been added to continuity equation to describe the retarded deformation of the pipe wall. The mechanical behavior of viscoelastic material is described by generalized Kelvin-Voigt model. The comparison of numerical simulation and experimental data from well known papers is presented. Short discussion about obtained results are given.

Pearling Instabilities of a Viscoelastic Thread

NASA Astrophysics Data System (ADS)

Deblais, A.; Velikov, K. P.; Bonn, D.

2018-05-01

Pearling instabilities of slender viscoelastic threads have received much attention, but remain incompletely understood. We study the instabilities in polymer solutions subject to uniaxial elongational flow. Two distinctly different instabilites are observed: beads on a string and blistering. The beads-on-a-string structure arises from a capillary instability whereas the blistering instability has a different origin: it is due to a coupling between stress and polymer concentration. By varying the temperature to change the solution properties we elucidate the interplay between flow and phase separation.

Non-linear vibrations of sandwich viscoelastic shells

NASA Astrophysics Data System (ADS)

Benchouaf, Lahcen; Boutyour, El Hassan; Daya, El Mostafa; Potier-Ferry, Michel

2018-04-01

This paper deals with the non-linear vibration of sandwich viscoelastic shell structures. Coupling a harmonic balance method with the Galerkin's procedure, one obtains an amplitude equation depending on two complex coefficients. The latter are determined by solving a classical eigenvalue problem and two linear ones. This permits to get the non-linear frequency and the non-linear loss factor as functions of the displacement amplitude. To validate our approach, these relationships are illustrated in the case of a circular sandwich ring.

Viscoelastic lower crust and mantle relaxation following the 14-16 April 2016 Kumamoto, Japan, earthquake sequence

NASA Astrophysics Data System (ADS)

Pollitz, Fred F.; Kobayashi, Tomokazu; Yarai, Hiroshi; Shibazaki, Bunichiro; Matsumoto, Takumi

2017-09-01

The 2016 Kumamoto, Japan, earthquake sequence, culminating in the Mw=7.0 16 April 2016 main shock, occurred within an active tectonic belt of central Kyushu. GPS data from GEONET reveal transient crustal motions from several millimeters per year up to ˜3 cm/yr during the first 8.5 months following the sequence. The spatial pattern of horizontal postseismic motions is shaped by both shallow afterslip and viscoelastic relaxation of the lower crust and upper mantle. We construct a suite of 2-D regional viscoelastic structures in order to derive an optimal joint afterslip and viscoelastic relaxation model using forward modeling of the viscoelastic relaxation. We find that afterslip dominates the postseismic relaxation in the near field (within 30 km of the main shock epicenter), while viscoelastic relaxation dominates at greater distance. The viscoelastic modeling strongly favors a very weak lower crust below a ˜65 km wide zone coinciding with the Beppu-Shimabara graben and the locus of central Kyushu volcanism. Inferred uppermost mantle viscosity is relatively low beneath southern Kyushu, consistent with independent inferences of a hydrated mantle wedge within the Nankai trough fore -arc.

Viscoelastic lower crust and mantle relaxation following the 14–16 April 2016 Kumamoto, Japan, earthquake sequence

USGS Publications Warehouse

Pollitz, Fred; Kobayashi, Tomokazu; Yarai, Hiroshi; Shibazaki, Bunichiro; Matsumoto, Takumi

2017-01-01

The 2016 Kumamoto, Japan, earthquake sequence, culminating in the Mw=7.0 16 April 2016 main shock, occurred within an active tectonic belt of central Kyushu. GPS data from GEONET reveal transient crustal motions from several millimeters per year up to ∼3 cm/yr during the first 8.5 months following the sequence. The spatial pattern of horizontal postseismic motions is shaped by both shallow afterslip and viscoelastic relaxation of the lower crust and upper mantle. We construct a suite of 2-D regional viscoelastic structures in order to derive an optimal joint afterslip and viscoelastic relaxation model using forward modeling of the viscoelastic relaxation. We find that afterslip dominates the postseismic relaxation in the near field (within 30 km of the main shock epicenter), while viscoelastic relaxation dominates at greater distance. The viscoelastic modeling strongly favors a very weak lower crust below a ∼65 km wide zone coinciding with the Beppu-Shimabara graben and the locus of central Kyushu volcanism. Inferred uppermost mantle viscosity is relatively low beneath southern Kyushu, consistent with independent inferences of a hydrated mantle wedge within the Nankai trough fore -arc.

Viscoelastic/damage modeling of filament-wound spherical pressure vessels

NASA Technical Reports Server (NTRS)

Hackett, Robert M.; Dozier, Jan D.

1987-01-01

A model of the viscoelastic/damage response of a filament-wound spherical vessel used for long-term pressure containment is developed. The matrix material of the composite system is assumed to be linearly viscoelastic. Internal accumulated damage based upon a quadratic relationship between transverse modulus and maximum circumferential strain is postulated. The resulting nonlinear problem is solved by an iterative routine. The elastic-viscoelastic correspondence is employed to produce, in the Laplace domain, the associated elastic solution for the maximum circumferential strain which is inverted by the method of collocation to yield the time-dependent solution. Results obtained with the model are compared to experimental observations.

Approximations of thermoelastic and viscoelastic control systems

NASA Technical Reports Server (NTRS)

Burns, J. A.; Liu, Z. Y.; Miller, R. E.

1990-01-01

Well-posed models and computational algorithms are developed and analyzed for control of a class of partial differential equations that describe the motions of thermo-viscoelastic structures. An abstract (state space) framework and a general well-posedness result are presented that can be applied to a large class of thermo-elastic and thermo-viscoelastic models. This state space framework is used in the development of a computational scheme to be used in the solution of a linear quadratic regulator (LQR) control problem. A detailed convergence proof is provided for the viscoelastic model and several numerical results are presented to illustrate the theory and to analyze problems for which the theory is incomplete.

Finite element analysis of left ventricle during cardiac cycles in viscoelasticity.

PubMed

Shen, Jing Jin; Xu, Feng Yu; Yang, Wen An

2016-08-01

To investigate the effect of myocardial viscoeslasticity on heart function, this paper presents a finite element model based on a hyper-viscoelastic model for the passive myocardium and Hill's three-element model for the active contraction. The hyper-viscoelastic model considers the myocardium microstructure, while the active model is phenomenologically based on the combination of Hill's equation for the steady tetanized contraction and the specific time-length-force property of the myocardial muscle. To validate the finite element model, the end-diastole strains and the end-systole strain predicted by the model are compared with the experimental values in the literature. It is found that the proposed model not only can estimate well the pumping function of the heart, but also predicts the transverse shear strains. The finite element model is also applied to analyze the influence of viscoelasticity on the residual stresses in the myocardium. Copyright © 2016 Elsevier Ltd. All rights reserved.

A 3D moisture-stress FEM analysis for time dependent problems in timber structures

NASA Astrophysics Data System (ADS)

Fortino, Stefania; Mirianon, Florian; Toratti, Tomi

2009-11-01

This paper presents a 3D moisture-stress numerical analysis for timber structures under variable humidity and load conditions. An orthotropic viscoelastic-mechanosorptive material model is specialized on the basis of previous models. Both the constitutive model and the equations needed to describe the moisture flow across the structure are implemented into user subroutines of the Abaqus finite element code and a coupled moisture-stress analysis is performed for several types of mechanical loads and moisture changes. The presented computational approach is validated by analyzing some wood tests described in the literature and comparing the computational results with the reported experimental data.

Jeans self gravitational instability of strongly coupled quantum plasma

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sharma, Prerana, E-mail: preranaiitd@rediffmail.com; Chhajlani, R. K.

2014-07-15

The Jeans self-gravitational instability is studied for quantum plasma composed of weakly coupled degenerate electron fluid and non-degenerate strongly coupled ion fluid. The formulation for such system is done on the basis of two fluid theory. The dynamics of weakly coupled degenerate electron fluid is governed by inertialess momentum equation. The quantum forces associated with the quantum diffraction effects and the quantum statistical effects act on the degenerate electron fluid. The strong correlation effects of ion are embedded in generalized viscoelastic momentum equation including the viscoelasticity and shear viscosities of ion fluid. The general dispersion relation is obtained using themore » normal mode analysis technique for the two regimes of propagation, i.e., hydrodynamic and kinetic regimes. The Jeans condition of self-gravitational instability is also obtained for both regimes, in the hydrodynamic regime it is observed to be affected by the ion plasma oscillations and quantum parameter while in the kinetic regime in addition to ion plasma oscillations and quantum parameter, it is also affected by the ion velocity which is modified by the viscosity generated compressional effects. The Jeans critical wave number and corresponding critical mass are also obtained for strongly coupled quantum plasma for both regimes.« less

Unified viscoelasticity: Applying discrete element models to soft tissues with two characteristic times.

PubMed

Anssari-Benam, Afshin; Bucchi, Andrea; Bader, Dan L

2015-09-18

Discrete element models have often been the primary tool in investigating and characterising the viscoelastic behaviour of soft tissues. However, studies have employed varied configurations of these models, based on the choice of the number of elements and the utilised formation, for different subject tissues. This approach has yielded a diverse array of viscoelastic models in the literature, each seemingly resulting in different descriptions of viscoelastic constitutive behaviour and/or stress-relaxation and creep functions. Moreover, most studies do not apply a single discrete element model to characterise both stress-relaxation and creep behaviours of tissues. The underlying assumption for this disparity is the implicit perception that the viscoelasticity of soft tissues cannot be described by a universal behaviour or law, resulting in the lack of a unified approach in the literature based on discrete element representations. This paper derives the constitutive equation for different viscoelastic models applicable to soft tissues with two characteristic times. It demonstrates that all possible configurations exhibit a unified and universal behaviour, captured by a single constitutive relationship between stress, strain and time as: σ+Aσ̇+Bσ¨=Pε̇+Qε¨. The ensuing stress-relaxation G(t) and creep J(t) functions are also unified and universal, derived as [Formula: see text] and J(t)=c2+(ε0-c2)e(-PQt)+σ0Pt, respectively. Application of these relationships to experimental data is illustrated for various tissues including the aortic valve, ligament and cerebral artery. The unified model presented in this paper may be applied to all tissues with two characteristic times, obviating the need for employing varied configurations of discrete element models in preliminary investigation of the viscoelastic behaviour of soft tissues. Copyright © 2015 Elsevier Ltd. All rights reserved.

Investigation of pitchfork bifurcation phenomena effects on heat transfer of viscoelastic flow inside a symmetric sudden expansion

NASA Astrophysics Data System (ADS)

Shahbani-Zahiri, A.; Hassanzadeh, H.; Shahmardan, M. M.; Norouzi, M.

2017-11-01

In this paper, the inertial and non-isothermal flows of the viscoelastic fluid through a planar channel with symmetric sudden expansion are numerically simulated. Effects of pitchfork bifurcation phenomena on the heat transfer rate are examined for the thermally developing and fully developed flow of the viscoelastic fluid inside the expanded part of the planar channel with an expansion ratio of 1:3. The rheological model of exponential Phan Thien-Tanner is used to include both the effects of shear-thinning and elasticity in fluid viscosity. The properties of fluids are temperature-dependent, and the viscous dissipation and heat stored by fluid elasticity are considered in the heat transfer equation. For coupling the governing equations, the PISO algorithm (Pressure Implicit with Splitting of Operator) is applied and the system of equations is linearized using the finite volume method on the collocated grids. The main purpose of this study is to examine the pitchfork bifurcation phenomena and its influences on the temperature distribution, the local and mean Nusselt numbers, and the first and second normal stress differences at different Reynolds, elasticity, and Brinkman numbers. The results show that by increasing the Brinkman number for the heated flow of the viscoelastic fluid inside the expanded part of the channel, the value of the mean Nusselt number is almost linearly decreased. Also, the maximum values of the local Nusselt number for the thermally developing flow and the local Nusselt number of the thermally fully developed flow are decremented by enhancing the Brinkman number.

Effect of a viscoelastic target on the impact response of a flat-nosed projectile

NASA Astrophysics Data System (ADS)

Liu, Hu; Yang, Jialing; Liu, Hua

2018-02-01

Taylor impact is a widely used strategy in which a flat-nosed projectile is fired onto a rigid anvil directly to determine the dynamic strength of rod specimens. Nowadays, the rigid anvil is often replaced by an output target bar to ensure the accuracy of measurement via recording strain signals in the output bar. For testing the dynamic strength of low-density materials, a low-impedance target bar, which exhibits viscoelastic characteristics is often employed. In this paper, an extended Taylor model is proposed to improve the idealization of treating the target bar as perfectly rigid material in the classic Taylor model, and the viscoelastic effect of the target bar is incorporated. The viscoelastic target bar is depicted by two elastic springs and one dashpot. Based on the plastic shock wave theory in the flat-nosed projectile associated with the viscoelastic wave analysis in the target bar, the viscoelastic effect of the target bar on the impact response of the flat-nosed projectile is investigated. The finite element simulation is also carried out to verify the theoretical model, and good agreement is found. The present theoretical model is also called the Taylor-cylinder Hopkinson impact, which provides a more accurate way to identify the dynamic material parameters. The dynamic responses of the present model are further compared with previous elastic and rigid target bar models. It is found that the viscoelastic effect of the target bar should be taken into consideration in the Taylor-cylinder Hopkinson impact test for low-impedance materials.

The development and validation of a numerical integration method for non-linear viscoelastic modeling

PubMed Central

Ramo, Nicole L.; Puttlitz, Christian M.

2018-01-01

Compelling evidence that many biological soft tissues display both strain- and time-dependent behavior has led to the development of fully non-linear viscoelastic modeling techniques to represent the tissue’s mechanical response under dynamic conditions. Since the current stress state of a viscoelastic material is dependent on all previous loading events, numerical analyses are complicated by the requirement of computing and storing the stress at each step throughout the load history. This requirement quickly becomes computationally expensive, and in some cases intractable, for finite element models. Therefore, we have developed a strain-dependent numerical integration approach for capturing non-linear viscoelasticity that enables calculation of the current stress from a strain-dependent history state variable stored from the preceding time step only, which improves both fitting efficiency and computational tractability. This methodology was validated based on its ability to recover non-linear viscoelastic coefficients from simulated stress-relaxation (six strain levels) and dynamic cyclic (three frequencies) experimental stress-strain data. The model successfully fit each data set with average errors in recovered coefficients of 0.3% for stress-relaxation fits and 0.1% for cyclic. The results support the use of the presented methodology to develop linear or non-linear viscoelastic models from stress-relaxation or cyclic experimental data of biological soft tissues. PMID:29293558

Effect of microstructure on the coupled electromagnetic-thermo-mechanical response of cyclotrimethylenetrinitramine-estane energetic aggregates to infrared laser radiation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Brown, Judith A.; Zikry, M. A., E-mail: zikry@ncsu.edu

2015-09-28

The coupled electromagnetic (EM)-thermo-mechanical response of cyclotrimethylenetrinitramine-estane energetic aggregates under laser irradiation and high strain rate loads has been investigated for various aggregate sizes and binder volume fractions. The cyclotrimethylenetrinitramine (RDX) crystals are modeled with a dislocation density-based crystalline plasticity formulation and the estane binder is modeled with finite viscoelasticity through a nonlinear finite element approach that couples EM wave propagation with laser heat absorption, thermal conduction, and inelastic deformation. Material property and local behavior mismatch at the crystal-binder interfaces resulted in geometric scattering of the EM wave, electric field and laser heating localization, high stress gradients, dislocation density, andmore » crystalline shear slip accumulation. Viscous sliding in the binder was another energy dissipation mechanism that reduced stresses in aggregates with thicker binder ligaments and larger binder volume fractions. This investigation indicates the complex interactions between EM waves and mechanical behavior, for accurate predictions of laser irradiation of heterogeneous materials.« less

Extending Maxwell's equations for dielectric materials using analytical principles from viscoelasticity based on the fractional calculus

NASA Astrophysics Data System (ADS)

Wharmby, Andrew William

Existing fractional calculus models having a non-empirical basis used to describe constitutive relationships between stress and strain in viscoelastic materials are modified to employ all orders of fractional derivatives between zero and one. Parallels between viscoelastic and dielectric theory are drawn so that these modified fractional calculus based models for viscoelastic materials may be used to describe relationships between electric flux density and electric field intensity in dielectric materials. The resulting fractional calculus based dielectric relaxation model is tested using existing complex permittivity data in the radio-frequency bandwidth of a wide variety of homogeneous materials. The consequences that the application of this newly developed fractional calculus based dielectric relaxation model has on Maxwell's equations are also examined through the effects of dielectric dissipation and dispersion.

Fractional calculus in biomechanics: a 3D viscoelastic model using regularized fractional derivative kernels with application to the human calcaneal fat pad.

PubMed

Freed, A D; Diethelm, K

2006-11-01

A viscoelastic model of the K-BKZ (Kaye, Technical Report 134, College of Aeronautics, Cranfield 1962; Bernstein et al., Trans Soc Rheol 7: 391-410, 1963) type is developed for isotropic biological tissues and applied to the fat pad of the human heel. To facilitate this pursuit, a class of elastic solids is introduced through a novel strain-energy function whose elements possess strong ellipticity, and therefore lead to stable material models. This elastic potential - via the K-BKZ hypothesis - also produces the tensorial structure of the viscoelastic model. Candidate sets of functions are proposed for the elastic and viscoelastic material functions present in the model, including two functions whose origins lie in the fractional calculus. The Akaike information criterion is used to perform multi-model inference, enabling an objective selection to be made as to the best material function from within a candidate set.

Capture mechanism in Palaeotropical pitcher plants (Nepenthaceae) is constrained by climate

PubMed Central

Moran, Jonathan A.; Gray, Laura K.; Clarke, Charles; Chin, Lijin

2013-01-01

Background and Aims Nepenthes (Nepenthaceae, approx. 120 species) are carnivorous pitcher plants with a centre of diversity comprising the Philippines, Borneo, Sumatra and Sulawesi. Nepenthes pitchers use three main mechanisms for capturing prey: epicuticular waxes inside the pitcher; a wettable peristome (a collar-shaped structure around the opening); and viscoelastic fluid. Previous studies have provided evidence suggesting that the first mechanism may be more suited to seasonal climates, whereas the latter two might be more suited to perhumid environments. In this study, this idea was tested using climate envelope modelling. Methods A total of 94 species, comprising 1978 populations, were grouped by prey capture mechanism (large peristome, small peristome, waxy, waxless, viscoelastic, non-viscoelastic, ‘wet’ syndrome and ‘dry’ syndrome). Nineteen bioclimatic variables were used to model habitat suitability at approx. 1 km resolution for each group, using Maxent, a presence-only species distribution modelling program. Key Results Prey capture groups putatively associated with perhumid conditions (large peristome, waxless, viscoelastic and ‘wet’ syndrome) had more restricted areas of probable habitat suitability than those associated putatively with less humid conditions (small peristome, waxy, non-viscoelastic and‘dry’ syndrome). Overall, the viscoelastic group showed the most restricted area of modelled suitable habitat. Conclusions The current study is the first to demonstrate that the prey capture mechanism in a carnivorous plant is constrained by climate. Nepenthes species employing peristome-based and viscoelastic fluid-based capture are largely restricted to perhumid regions; in contrast, the wax-based mechanism allows successful capture in both perhumid and more seasonal areas. Possible reasons for the maintenance of peristome-based and viscoelastic fluid-based capture mechanisms in Nepenthes are discussed in relation to the costs and benefits associated with a given prey capture strategy. PMID:23975653

Capture mechanism in Palaeotropical pitcher plants (Nepenthaceae) is constrained by climate.

PubMed

Moran, Jonathan A; Gray, Laura K; Clarke, Charles; Chin, Lijin

2013-11-01

Nepenthes (Nepenthaceae, approx. 120 species) are carnivorous pitcher plants with a centre of diversity comprising the Philippines, Borneo, Sumatra and Sulawesi. Nepenthes pitchers use three main mechanisms for capturing prey: epicuticular waxes inside the pitcher; a wettable peristome (a collar-shaped structure around the opening); and viscoelastic fluid. Previous studies have provided evidence suggesting that the first mechanism may be more suited to seasonal climates, whereas the latter two might be more suited to perhumid environments. In this study, this idea was tested using climate envelope modelling. A total of 94 species, comprising 1978 populations, were grouped by prey capture mechanism (large peristome, small peristome, waxy, waxless, viscoelastic, non-viscoelastic, 'wet' syndrome and 'dry' syndrome). Nineteen bioclimatic variables were used to model habitat suitability at approx. 1 km resolution for each group, using Maxent, a presence-only species distribution modelling program. Prey capture groups putatively associated with perhumid conditions (large peristome, waxless, viscoelastic and 'wet' syndrome) had more restricted areas of probable habitat suitability than those associated putatively with less humid conditions (small peristome, waxy, non-viscoelastic and'dry' syndrome). Overall, the viscoelastic group showed the most restricted area of modelled suitable habitat. The current study is the first to demonstrate that the prey capture mechanism in a carnivorous plant is constrained by climate. Nepenthes species employing peristome-based and viscoelastic fluid-based capture are largely restricted to perhumid regions; in contrast, the wax-based mechanism allows successful capture in both perhumid and more seasonal areas. Possible reasons for the maintenance of peristome-based and viscoelastic fluid-based capture mechanisms in Nepenthes are discussed in relation to the costs and benefits associated with a given prey capture strategy.

Viscoelastic effect on acoustic band gaps in polymer-fluid composites

NASA Astrophysics Data System (ADS)

Merheb, B.; Deymier, P. A.; Muralidharan, K.; Bucay, J.; Jain, M.; Aloshyna-Lesuffleur, M.; Greger, R. W.; Mohanty, S.; Berker, A.

2009-10-01

In this paper, we present a theoretical analysis of the propagation of acoustic waves through elastic and viscoelastic two-dimensional phononic crystal structures. Numerical calculations of transmission spectra are conducted by extending the finite-difference-time-domain method to account for linear viscoelastic materials with time-dependent moduli. We study a phononic crystal constituted of a square array of cylindrical air inclusions in a solid viscoelastic matrix. The elastic properties of the solid are those of a silicone rubber. This system exhibits very wide band gaps in its transmission spectrum that extend to frequencies in the audible range of the spectrum. These gaps are characteristic of fluid matrix/air inclusion systems and result from the very large contrast between the longitudinal and transverse speeds of sound in rubber. By treating the matrix as a viscoelastic medium within the standard linear solid (SLS) model, we demonstrate that viscoelasticity impacts the transmission properties of the rubber/air phononic crystal not only by attenuating the transmitted acoustic waves but also by shifting the passing bands frequencies toward lower values. The ranges of frequencies exhibiting attenuation or frequency shift are determined by the value of the relaxation time in the SLS model. We show that viscoelasticity can be used to decrease the frequency of pass bands (and consequently stop bands) in viscoelastic/air phononic crystals.

Viscoelastic love-type surface waves

USGS Publications Warehouse

Borcherdt, Roger D.

2008-01-01

The general theoretical solution for Love-Type surface waves in viscoelastic media provides theoreticalexpressions for the physical characteristics of the waves in elastic as well as anelastic media with arbitraryamounts of intrinsic damping. The general solution yields dispersion and absorption-coefficient curves for the waves as a function of frequency and theamount of intrinsic damping for any chosen viscoelastic model.Numerical results valid for a variety of viscoelastic models provide quantitative estimates of the physicalcharacteristics of the waves pertinent to models of Earth materials ranging from small amounts of damping in the Earth’s crust to moderate and large amounts of damping in soft soils and water-saturated sediments. Numerical results, presented herein, are valid for a wide range of solids and applications.

Time-Varying Upper-Plate Deformation during the Megathrust Subduction Earthquake Cycle

NASA Astrophysics Data System (ADS)

Furlong, Kevin P.; Govers, Rob; Herman, Matthew

2015-04-01

Over the past several decades of the WEGENER era, our abilities to observe and image the deformational behavior of the upper plate in megathrust subduction zones has dramatically improved. Several intriguing inferences can be made from these observations including apparent lateral variations in locking along subduction zones, which differs from interseismic to coseismic periods; the significant magnitude of post-earthquake deformation (e.g. following the 20U14 Mw Iquique, Chile earthquake, observed on-land GPS post-EQ displacements are comparable to the co-seismic displacements); and incompatibilities between rates of slip deficit accumulation and resulting earthquake co-seismic slip (e.g. pre-Tohoku, inferred rates of slip deficit accumulation on the megathrust significantly exceed slip amounts for the ~ 1000 year recurrence.) Modeling capabilities have grown from fitting simple elastic accumulation/rebound curves to sparse data to having spatially dense continuous time series that allow us to infer details of plate boundary coupling, rheology-driven transient deformation, and partitioning among inter-earthquake and co-seismic displacements. In this research we utilize a 2D numerical modeling to explore the time-varying deformational behavior of subduction zones during the earthquake cycle with an emphasis on upper-plate and plate interface behavior. We have used a simplified model configuration to isolate fundamental processes associated with the earthquake cycle, rather than attempting to fit details of specific megathrust zones. Using a simple subduction geometry, but realistic rheologic layering we are evaluating the time-varying displacement and stress response through a multi-earthquake cycle history. We use a simple model configuration - an elastic subducting slab, an elastic upper plate (shallower than 40 km), and a visco-elastic upper plate (deeper than 40 km). This configuration leads to an upper plate that acts as a deforming elastic beam at inter-earthquake loading times and rates with a viscously relaxed regime at depths greater than 40 km. Analyses of our preliminary model results lead to the following: 1. Co-seismic stress transfer from the unloading elastic layer (shallow) into an elastically loading visco-elastic layer (deeper) - extends ~ 100 km inboard of locked zone. This stress transfer affects both coseismic and post-seismic surface displacements. 2. Post-seismic response of upper plate involves seaward motion for initial 10-20 years (~ 2 Maxwell times) after EQ. This occurs in spite of there being no slip on locked plate boundary - i.e. this is not plate boundary after-slip but rather is a consequence of stress relaxation in co-seismically loaded visco-elastic layer. However standard inversions of the surface displacement field would indicate significant after-slip along the locked plate interface. 3. By approximately 80 years (8 Maxwell times) system has returned to simple linear displacement pattern - the expected behavior for a shortening elastic beam. Prior to that time, the surface (observable) displacement pattern changes substantially over time and would result in an apparent temporal variation in coupling - from near-zero coupling to fully locked over ~ 80 years post-earthquake. These preliminary results indicate that care is needed in interpreting observed surface displacement fields from GPS, InSAR, etc. during the interseismic period. temporal variations in crustal deformation observed in regions such as the recent Tohoku, Maule, and Iquique megathrust events which are ascribed to fault plane after-slip may in fact reflect processes associated with re-equilibration of the visco-elastic subduction system.

Mapping Viscoelastic and Plastic Properties of Polymers and Polymer-Nanotube Composites using Instrumented Indentation

PubMed Central

Gayle, Andrew J.; Cook, Robert F.

2016-01-01

An instrumented indentation method is developed for generating maps of time-dependent viscoelastic and time-independent plastic properties of polymeric materials. The method is based on a pyramidal indentation model consisting of two quadratic viscoelastic Kelvin-like elements and a quadratic plastic element in series. Closed-form solutions for indentation displacement under constant load and constant loading-rate are developed and used to determine and validate material properties. Model parameters are determined by point measurements on common monolithic polymers. Mapping is demonstrated on an epoxy-ceramic interface and on two composite materials consisting of epoxy matrices containing multi-wall carbon nanotubes. A fast viscoelastic deformation process in the epoxy was unaffected by the inclusion of the nanotubes, whereas a slow viscoelastic process was significantly impeded, as was the plastic deformation. Mapping revealed considerable spatial heterogeneity in the slow viscoelastic and plastic responses in the composites, particularly in the material with a greater fraction of nanotubes. PMID:27563168

Temperature-dependent poroelastic and viscoelastic effects on microscale-modelling of seismic reflections in heavy oil reservoirs

NASA Astrophysics Data System (ADS)

Ciz, Radim; Saenger, Erik H.; Gurevich, Boris; Shapiro, Serge A.

2009-03-01

We develop a new model for elastic properties of rocks saturated with heavy oil. The heavy oil is represented by a viscoelastic material, which at low frequencies and/or high temperatures behaves as a Newtonian fluid, and at high frequencies and/or low temperatures as a nearly elastic solid. The bulk and shear moduli of a porous rock saturated with such viscoelastic material are then computed using approximate extended Gassmann equations of Ciz and Shapiro by replacing the elastic moduli of the pore filling material with complex and frequency-dependent moduli of the viscoelastic pore fill. We test the proposed model by comparing its predictions with numerical simulations based on a direct finite-difference solution of equations of dynamic viscoelasticity. The simulations are performed for the reflection coefficient from an interface between a homogeneous fluid and a porous medium. The numerical tests are performed both for an idealized porous medium consisting of alternating solid and viscoelastic layers, and for a more realistic 3-D geometry of the pore space. Both sets of numerical tests show a good agreement between the predictions of the proposed viscoelastic workflow and numerical simulations for relatively high viscosities where viscoelastic effects are important. The results confirm that application of extended Gassmann equations in conjunction with the complex and frequency-dependent moduli of viscoelastic pore filling material, such as heavy oil, provides a good approximation for the elastic moduli of rocks saturated with such material. By construction, this approximation is exactly consistent with the classical Gassmann's equation for sufficiently low frequencies or high temperature when heavy oil behaves like a fluid. For higher frequencies and/or lower temperatures, the predictions are in good agreement with the direct numerical solution of equations of dynamic viscoelasticity on the microscale. This demonstrates that the proposed methodology provides realistic estimates of elastic properties of heavy oil rocks.

Pure axial flow of viscoelastic fluids in rectangular microchannels under combined effects of electro-osmosis and hydrodynamics

NASA Astrophysics Data System (ADS)

Reshadi, Milad; Saidi, Mohammad Hassan; Ebrahimi, Abbas

2018-02-01

This paper presents an analysis of the combined electro-osmotic and pressure-driven axial flows of viscoelastic fluids in a rectangular microchannel with arbitrary aspect ratios. The rheological behavior of the fluid is described by the complete form of Phan-Thien-Tanner (PTT) model with the Gordon-Schowalter convected derivative which covers the upper convected Maxwell, Johnson-Segalman and FENE-P models. Our numerical simulation is based on the computation of 2D Poisson-Boltzmann, Cauchy momentum and PTT constitutive equations. The solution of these governing nonlinear coupled set of equations is obtained by using the second-order central finite difference method in a non-uniform grid system and is verified against 1D analytical solution of the velocity profile with less than 0.06% relative error. Also, a parametric study is carried out to investigate the effect of channel aspect ratio (width to height), wall zeta potential and the Debye-Hückel parameter on 2D velocity profile, volumetric flow rate and the Poiseuille number in the mixed EO/PD flows of viscoelastic fluids with different Weissenberg numbers. Our results show that, for low channel aspect ratios, the previous 1D analytical models underestimate the velocity profile at the channel half-width centerline in the case of favorable pressure gradients and overestimate it in the case of adverse pressure gradients. The results reveal that the inapplicability of the Debye-Hückel approximation at high zeta potentials is more significant for higher Weissenberg number fluids. Also, it is found that, under the specified values of electrokinetic parameters, there is a threshold for velocity scale ratio in which the Poiseuille number is approximately independent of channel aspect ratio.

Micro-Macro Simulation of Viscoelastic Fluids in Three Dimensions

NASA Astrophysics Data System (ADS)

Rüttgers, Alexander; Griebel, Michael

2012-11-01

The development of the chemical industry resulted in various complex fluids that cannot be correctly described by classical fluid mechanics. For instance, this includes paint, engine oils with polymeric additives and toothpaste. We currently perform multiscale viscoelastic flow simulations for which we have coupled our three-dimensional Navier-Stokes solver NaSt3dGPF with the stochastic Brownian configuration field method on the micro-scale. In this method, we represent a viscoelastic fluid as a dumbbell system immersed in a three-dimensional Newtonian liquid which leads to a six-dimensional problem in space. The approach requires large computational resources and therefore depends on an efficient parallelisation strategy. Our flow solver is parallelised with a domain decomposition approach using MPI. It shows excellent scale-up results for up to 128 processors. In this talk, we present simulation results for viscoelastic fluids in square-square contractions due to their relevance for many engineering applications such as extrusion. Another aspect of the talk is the parallel implementation in NaSt3dGPF and the parallel scale-up and speed-up behaviour.

Measurements of vocal fold tissue viscoelasticity: Approaching the male phonatory frequency range

NASA Astrophysics Data System (ADS)

Chan, Roger W.

2004-06-01

Viscoelastic shear properties of human vocal fold tissues have been reported previously. However, data have only been obtained at very low frequencies (<=15 Hz). This necessitates data extrapolation to the frequency range of phonation based on constitutive modeling and time-temperature superposition. This study attempted to obtain empirical measurements at higher frequencies with the use of a controlled strain torsional rheometer, with a design of directly controlling input strain that introduced significantly smaller system inertial errors compared to controlled stress rheometry. Linear viscoelastic shear properties of the vocal fold mucosa (cover) from 17 canine larynges were quantified at frequencies of up to 50 Hz. Consistent with previous data, results showed that the elastic shear modulus (G'), viscous shear modulus (G''), and damping ratio (ζ) of the vocal fold mucosa were relatively constant across 0.016-50 Hz, whereas the dynamic viscosity (ɛ') decreased monotonically with frequency. Constitutive characterization of the empirical data by a quasilinear viscoelastic model and a statistical network model demonstrated trends of viscoelastic behavior at higher frequencies generally following those observed at lower frequencies. These findings supported the use of controlled strain rheometry for future investigations of the viscoelasticity of vocal fold tissues and phonosurgical biomaterials at phonatory frequencies.

Quantitative sonoelastography for the in vivo assessment of skeletal muscle viscoelasticity

NASA Astrophysics Data System (ADS)

Hoyt, Kenneth; Kneezel, Timothy; Castaneda, Benjamin; Parker, Kevin J.

2008-08-01

A novel quantitative sonoelastography technique for assessing the viscoelastic properties of skeletal muscle tissue was developed. Slowly propagating shear wave interference patterns (termed crawling waves) were generated using a two-source configuration vibrating normal to the surface. Theoretical models predict crawling wave displacement fields, which were validated through phantom studies. In experiments, a viscoelastic model was fit to dispersive shear wave speed sonoelastographic data using nonlinear least-squares techniques to determine frequency-independent shear modulus and viscosity estimates. Shear modulus estimates derived using the viscoelastic model were in agreement with that obtained by mechanical testing on phantom samples. Preliminary sonoelastographic data acquired in healthy human skeletal muscles confirm that high-quality quantitative elasticity data can be acquired in vivo. Studies on relaxed muscle indicate discernible differences in both shear modulus and viscosity estimates between different skeletal muscle groups. Investigations into the dynamic viscoelastic properties of (healthy) human skeletal muscles revealed that voluntarily contracted muscles exhibit considerable increases in both shear modulus and viscosity estimates as compared to the relaxed state. Overall, preliminary results are encouraging and quantitative sonoelastography may prove clinically feasible for in vivo characterization of the dynamic viscoelastic properties of human skeletal muscle.

Identifying viscoelastic parameters of tissue specimens using Hertz contact mechanics

NASA Astrophysics Data System (ADS)

Namiri, Nikan K.; Maccabi, Ashkan; Bajwa, Neha; Badran, Karam W.; St. John, Maie A.; Taylor, Zachary D.; Grundfest, Warren S.; Saddik, George N.

2018-02-01

The unique viscoelastic properties of tissues throughout the human body can be utilized in a variety of clinical applications. Palpation techniques, for instance, enable surgeons to distinguish malignancies in tissue composition during surgical procedures. Additionally, imaging devices have begun utilizing the viscoelastic properties of tissue to delineate tumor margins. Vibroacoustography (VA), a non-invasive, high resolution imaging modality, has the ability to detect sub-millimeter differences in tissue composition. VA images tissue using a low frequency acoustic radiation force, which perturbs the target and causes an acoustic response that is dependent on the target's viscoelastic properties. Given the unique properties specific to human and animal tissues, there are far-reaching clinical applications of VA. To date, however, a comprehensive model that relates viscoelasticity to VA tissue response has yet to be developed. Utilizing tissue-mimicking phantoms (TMPs) and fresh ex vivo tissues, a mechanical stress relaxation model was developed to compare the viscoelastic properties of known and unknown specimens. This approach was conducted using the Hertz theory of contact mechanics. Fresh hepatic tissue was obtained from porcine subjects (n=10), while gelatin and agar TMPs (n=12) were fabricated from organic extracts. Each specimen's elastic modulus (E), long term shear modulus (η), and time constant (τ) were found to be unique. Additionally, each specimen's stress relaxation profiles were analyzed using Weichert-Maxwell viscoelastic modeling, and retained high precision (R2>0.9) among all samples.

Postseismic viscoelastic deformation and stress. Part 2: Stress theory and computation; dependence of displacement, strain, and stress on fault parameters

NASA Technical Reports Server (NTRS)

Cohen, S. C.

1979-01-01

A viscoelastic model for deformation and stress associated with earthquakes is reported. The model consists of a rectangular dislocation (strike slip fault) in a viscoelastic layer (lithosphere) lying over a viscoelastic half space (asthenosphere). The time dependent surface stresses are analyzed. The model predicts that near the fault a significant fraction of the stress that was reduced during the earthquake is recovered by viscoelastic softening of the lithosphere. By contrast, the strain shows very little change near the fault. The model also predicts that the stress changes associated with asthenospheric flow extend over a broader region than those associated with lithospheric relaxation even though the peak value is less. The dependence of the displacements, stresses on fault parameters studied. Peak values of strain and stress drop increase with increasing fault height and decrease with fault depth. Under many circumstances postseismic strains and stresses show an increase with decreasing depth to the lithosphere-asthenosphere boundary. Values of the strain and stress at distant points from the fault increase with fault area but are relatively insensitive to fault depth.

Magnetic susceptibility, nanorheology, and magnetoviscosity of magnetic nanoparticles in viscoelastic environments

NASA Astrophysics Data System (ADS)

Ilg, Patrick; Evangelopoulos, Apostolos E. A. S.

2018-03-01

While magnetic nanoparticles suspended in Newtonian solvents (ferrofluids) have been intensively studied in recent years, the effects of viscoelasticity of the surrounding medium on the nanoparticle dynamics are much less understood. Here we investigate a mesoscopic model for the orientational dynamics of isolated magnetic nanoparticles subject to external fields, viscous and viscoelastic friction, as well as the corresponding random torques. We solve the model analytically in the overdamped limit for weak viscoelasticity. By comparison to Brownian dynamics simulations we establish the limits of validity of the analytical solution. We find that viscoelasticity not only slows down the magnetization relaxation, shifts the peak of the imaginary magnetic susceptibility χ″ to lower frequencies, and increases the magnetoviscosity but also leads to nonexponential relaxation and a broadening of χ″. The model we study also allows us to test a recent proposal for using magnetic susceptibility measurements as a nanorheological tool using a variant of the Germant-DiMarzio-Bishop relation. We find for the present model and certain parameter ranges that the relation of the magnetic susceptibility to the shear modulus is satisfied to a good approximation.

A viscoelastic–stochastic model of the effects of cytoskeleton remodelling on cell adhesion

PubMed Central

Li, Long; Zhang, Wenyan

2016-01-01

Cells can adapt their mechanical properties through cytoskeleton remodelling in response to external stimuli when the cells adhere to the extracellular matrix (ECM). Many studies have investigated the effects of cell and ECM elasticity on cell adhesion. However, experiments determined that cells are viscoelastic and exhibiting stress relaxation, and the mechanism behind the effect of cellular viscoelasticity on the cell adhesion behaviour remains unclear. Therefore, we propose a theoretical model of a cluster of ligand–receptor bonds between two dissimilar viscoelastic media subjected to an applied tensile load. In this model, the distribution of interfacial traction is assumed to follow classical continuum viscoelastic equations, whereas the rupture and rebinding of individual molecular bonds are governed by stochastic equations. On the basis of this model, we determined that viscosity can significantly increase the lifetime, stability and dynamic strength of the adhesion cluster of molecular bonds, because deformation relaxation attributed to the viscoelastic property can increase the rebinding probability of each open bond and reduce the stress concentration in the adhesion area. PMID:27853571

A new analytical method for estimating lumped parameter constants of linear viscoelastic models from strain rate tests

NASA Astrophysics Data System (ADS)

Mattei, G.; Ahluwalia, A.

2018-04-01

We introduce a new function, the apparent elastic modulus strain-rate spectrum, E_{app} ( \\dot{ɛ} ), for the derivation of lumped parameter constants for Generalized Maxwell (GM) linear viscoelastic models from stress-strain data obtained at various compressive strain rates ( \\dot{ɛ}). The E_{app} ( \\dot{ɛ} ) function was derived using the tangent modulus function obtained from the GM model stress-strain response to a constant \\dot{ɛ} input. Material viscoelastic parameters can be rapidly derived by fitting experimental E_{app} data obtained at different strain rates to the E_{app} ( \\dot{ɛ} ) function. This single-curve fitting returns similar viscoelastic constants as the original epsilon dot method based on a multi-curve global fitting procedure with shared parameters. Its low computational cost permits quick and robust identification of viscoelastic constants even when a large number of strain rates or replicates per strain rate are considered. This method is particularly suited for the analysis of bulk compression and nano-indentation data of soft (bio)materials.

Mathematical model of Rayleigh-Taylor and Richtmyer-Meshkov instabilities for viscoelastic fluids

NASA Astrophysics Data System (ADS)

Rollin, Bertrand; Andrews, Malcolm J.

2011-04-01

We extended the Goncharov model [V. N. Goncharov, Phys. Rev. Lett.PRLTAO0031-900710.1103/PhysRevLett.88.134502 88, 134502 (2002)] for nonlinear Rayleigh-Taylor instability of perfect fluids to the case of Rivlin-Ericksen viscoelastic fluids [R. S. Rivlin and J. L. Ericksen, Rat. Mech. Anal. 4, 323 (1955)], with surface tension. For Rayleigh-Taylor instability, viscosity, surface tension, and viscoelasticity decrease the exponential growth rate predicted by linear stability analysis. In particular, we find that viscosity and surface tension decrease the terminal bubble velocity, whereas viscoelasticity is found to have no effect. All three properties increase the saturation height of the bubble. In Richmyer-Meshkov instability, the decay of the asymptotic velocity depends on the balance between viscosity and surface tension, and viscoelasticity tends to slow the asymptotic velocity decay.

A multiscale model for predicting the viscoelastic properties of asphalt concrete

NASA Astrophysics Data System (ADS)

Garcia Cucalon, Lorena; Rahmani, Eisa; Little, Dallas N.; Allen, David H.

2016-08-01

It is well known that the accurate prediction of long term performance of asphalt concrete pavement requires modeling to account for viscoelasticity within the mastic. However, accounting for viscoelasticity can be costly when the material properties are measured at the scale of asphalt concrete. This is due to the fact that the material testing protocols must be performed recursively for each mixture considered for use in the final design.

Viscoelastic propellant effects on Space Shuttle Dynamics

NASA Technical Reports Server (NTRS)

Bugg, F.

1981-01-01

The program of solid propellant research performed in support of the space shuttle dynamics modeling effort is described. Stiffness, damping, and compressibility of the propellant and the effects of many variables on these properties are discussed. The relationship between the propellant and solid rocket booster dynamics during liftoff and boost flight conditions and the effects of booster vibration and propellant stiffness on free free solid rocket booster modes are described. Coupled modes of the shuttle system and the effect of propellant stiffness on the interfaces of the booster and the external tank are described. A finite shell model of the solid rocket booster was developed.

Viscoelastic material inversion using Sierra-SD and ROL

DOE Office of Scientific and Technical Information (OSTI.GOV)

Walsh, Timothy; Aquino, Wilkins; Ridzal, Denis

2014-11-01

In this report we derive frequency-domain methods for inverse characterization of the constitutive parameters of viscoelastic materials. The inverse problem is cast in a PDE-constrained optimization framework with efficient computation of gradients and Hessian vector products through matrix free operations. The abstract optimization operators for first and second derivatives are derived from first principles. Various methods from the Rapid Optimization Library (ROL) are tested on the viscoelastic inversion problem. The methods described herein are applied to compute the viscoelastic bulk and shear moduli of a foam block model, which was recently used in experimental testing for viscoelastic property characterization.

Viscoelastic modeling of deformation and gravity changes induced by pressurized magmatic sources

NASA Astrophysics Data System (ADS)

Currenti, Gilda

2018-05-01

Gravity and height changes, which reflect magma accumulation in subsurface chambers, are evaluated using analytical and numerical models in order to investigate their relationships and temporal evolutions. The analysis focuses mainly on the exploration of the time-dependent response of gravity and height changes to the pressurization of ellipsoidal magmatic chambers in viscoelastic media. Firstly, the validation of the numerical Finite Element results is performed by comparison with analytical solutions, which are devised for a simple spherical source embedded in a homogeneous viscoelastic half-space medium. Then, the effect of several model parameters on time-dependent height and gravity changes is investigated thanks to the flexibility of the numerical method in handling complex configurations. Both homogeneous and viscoelastic shell models reveal significantly different amplitudes in the ratio between gravity and height changes depending on geometry factors and medium rheology. The results show that these factors also influence the relaxation characteristic times of the investigated geophysical changes. Overall, these temporal patterns are compatible with time-dependent height and gravity changes observed on Etna volcano during the 1994-1997 inflation period. By modeling the viscoelastic response of a pressurized prolate magmatic source, a general agreement between computed and observed geophysical variations is achieved.

Structural Model for Viscoelastic Properties of Pericardial Bioprosthetic Valves.

PubMed

Rassoli, Aisa; Fatouraee, Nasser; Guidoin, Robert

2018-03-30

The benefit of bioprosthetic aortic valve over mechanical valve replacements is the release of thromboembolism and digression of long-term anticoagulation treatment. The function of bioprostheses and their efficiency is known to depend on the mechanical properties of the leaflet tissue. So it is necessary to select a suitable tissue for the bioprosthesis. The purpose of the present study is to clarify the viscoelastic behavior of bovine, equine, and porcine pericardium. In this study, pericardiums were compared mechanically from the viscoelastic aspect. After fixation of the tissues in glutaraldehyde, first uniaxial tests with different extension rates in the fiber direction were performed. Then, the stress relaxation tests in the fiber direction were done on these pericardial tissues by exerting 20, 30,40, and 50% strains. After evaluation of viscoelastic linearity, the Prony series, quasilinear viscoelastic (QLV) and modified superposition theory were applied to the stress relaxation data. Finally, the parameters of these constitutive models were extracted for each pericardium tissue. All three tissues exhibited a decrease in relaxation rate with elevating strain, indicating the nonlinear viscoelastic behavior of these tissues. The three-term Prony model was selected for describing the linear viscoelasticity. Among different models, the QLV model was best able to capture the relaxation behavior of the pericardium tissues. More stiffness of porcine pericardium was observed in comparison to the two other pericardium tissues. The relaxation percentage of porcine pericardium was less than the two others. It can be concluded that porcine pericardium behaves more as an elastic and less like a viscous tissue in comparison to the bovine and equine pericardium. © 2018 International Center for Artificial Organs and Transplantation and Wiley Periodicals, Inc.

Airway mechanics and lung tissue viscoelasticity: effects of altered blood hematocrit in the pulmonary circulation.

PubMed

Peták, Ferenc; Fodor, Gergely H; Babik, Barna; Habre, Walid

2016-07-01

The contribution of the hematocrit (Hct) of the blood in the pulmonary vasculature to the overall lung mechanics has not been characterized. We therefore set out to establish how changes of the Hct level in the pulmonary circulation affect the airway and lung tissue viscoelastic properties. The Hct level of the blood in an isolated perfused rat lung model was randomly altered. Intermediate (26.5%), followed by low (6.6%) or normal (43.7%), Hct was set in two consecutive sequences. The pulmonary capillary pressure was maintained constant throughout the experiment, and the pulmonary hemodynamic parameters were monitored continuously. The airway resistance (Raw), the viscous (G) and elastic (H) parameters, and the hysteresivity (η = G/H) of the lung tissues were obtained from measurements of forced oscillatory input impedance data. Raw was not affected by the alterations of the Hct levels. As concerns the lung tissues, the decrease of Hct to intermediate or low levels resulted in close to proportional decreases in the viscoelastic parameters G [16.5 ± 7.7% (SD), 12.1 ± 9.5%, P < 0.005] and H (13.2 ± 8.6%, 10.8 ± 4.7%, P < 0.001). No significant changes in η were detected in a wide range of Hct, which indicates that coupled processes cause alterations in the resistive and elastic properties of the lungs following Hct changes in the pulmonary circulation. The diminishment of the viscous and elastic parameters of the pulmonary parenchyma following a reduction of blood Hct demonstrates the significant contribution of the red blood cells to the overall lung viscoelasticity. Copyright © 2016 the American Physiological Society.

[Research on adaptive quasi-linear viscoelastic model for nonlinear viscoelastic properties of in vivo soft tissues].

PubMed

Wang, Heng; Sang, Yuanjun

2017-10-01

The mechanical behavior modeling of human soft biological tissues is a key issue for a large number of medical applications, such as surgery simulation, surgery planning, diagnosis, etc. To develop a biomechanical model of human soft tissues under large deformation for surgery simulation, the adaptive quasi-linear viscoelastic (AQLV) model was proposed and applied in human forearm soft tissues by indentation tests. An incremental ramp-and-hold test was carried out to calibrate the model parameters. To verify the predictive ability of the AQLV model, the incremental ramp-and-hold test, a single large amplitude ramp-and-hold test and a sinusoidal cyclic test at large strain amplitude were adopted in this study. Results showed that the AQLV model could predict the test results under the three kinds of load conditions. It is concluded that the AQLV model is feasible to describe the nonlinear viscoelastic properties of in vivo soft tissues under large deformation. It is promising that this model can be selected as one of the soft tissues models in the software design for surgery simulation or diagnosis.

Modelling nonlinear viscoelastic behaviours of loudspeaker suspensions-like structures

NASA Astrophysics Data System (ADS)

Maillou, Balbine; Lotton, Pierrick; Novak, Antonin; Simon, Laurent

2018-03-01

Mechanical properties of an electrodynamic loudspeaker are mainly determined by its suspensions (surround and spider) that behave nonlinearly and typically exhibit frequency dependent viscoelastic properties such as creep effect. The paper aims at characterizing the mechanical behaviour of electrodynamic loudspeaker suspensions at low frequencies using nonlinear identification techniques developed in recent years. A Generalized Hammerstein based model can take into account both frequency dependency and nonlinear properties. As shown in the paper, the model generalizes existing nonlinear or viscoelastic models commonly used for loudspeaker modelling. It is further experimentally shown that a possible input-dependent law may play a key role in suspension characterization.

A model for longitudinal and shear wave propagation in viscoelastic media

PubMed

Szabo; Wu

2000-05-01

Relaxation models fail to predict and explain loss characteristics of many viscoelastic materials which follow a frequency power law. A model based on a time-domain statement of causality is presented that describes observed power-law behavior of many viscoelastic materials. A Hooke's law is derived from power-law loss characteristics; it reduces to the Hooke's law for the Voigt model for the specific case of quadratic frequency loss. Broadband loss and velocity data for both longitudinal and shear elastic types of waves agree well with predictions. These acoustic loss models are compared to theories for loss mechanisms in dielectrics based on isolated polar molecules and cooperative interactions.

Fractal Folding and Medium Viscoelasticity Contribute Jointly to Chromosome Dynamics

NASA Astrophysics Data System (ADS)

Polovnikov, K. E.; Gherardi, M.; Cosentino-Lagomarsino, M.; Tamm, M. V.

2018-02-01

Chromosomes are key players of cell physiology, their dynamics provides valuable information about its physical organization. In both prokaryotes and eukaryotes, the short-time motion of chromosomal loci has been described with a Rouse model in a simple or viscoelastic medium. However, little emphasis has been put on the influence of the folded organization of chromosomes on the local dynamics. Clearly, stress propagation, and thus dynamics, must be affected by such organization, but a theory allowing us to extract such information from data, e.g., on two-point correlations, is lacking. Here, we describe a theoretical framework able to answer this general polymer dynamics question. We provide a scaling analysis of the stress-propagation time between two loci at a given arclength distance along the chromosomal coordinate. The results suggest a precise way to assess folding information from the dynamical coupling of chromosome segments. Additionally, we realize this framework in a specific model of a polymer whose long-range interactions are designed to make it fold in a fractal way and immersed in a medium characterized by subdiffusive fractional Langevin motion with a tunable scaling exponent. This allows us to derive explicit analytical expressions for the correlation functions.

Viscoelastic coupling of nanoelectromechanical resonators.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Simonson, Robert Joseph; Staton, Alan W.

2009-09-01

This report summarizes work to date on a new collaboration between Sandia National Laboratories and the California Institute of Technology (Caltech) to utilize nanoelectromechanical resonators designed at Caltech as platforms to measure the mechanical properties of polymeric materials at length scales on the order of 10-50 nm. Caltech has succeeded in reproducibly building cantilever resonators having major dimensions on the order of 2-5 microns. These devices are fabricated in pairs, with free ends separated by reproducible gaps having dimensions on the order of 10-50 nm. By controlled placement of materials that bridge the very small gap between resonators, the mechanicalmore » devices become coupled through the test material, and the transmission of energy between the devices can be monitored. This should allow for measurements of viscoelastic properties of polymeric materials at high frequency over short distances. Our work to date has been directed toward establishing this measurement capability at Sandia.« less

Ultrasonic wave propagation in viscoelastic cortical bone plate coupled with fluids: a spectral finite element study.

PubMed

Nguyen, Vu-Hieu; Naili, Salah

2013-01-01

This work deals with the ultrasonic wave propagation in the cortical layer of long bones which is known as being a functionally graded anisotropic material coupled with fluids. The viscous effects are taken into account. The geometrical configuration mimics the one of axial transmission technique used for evaluating the bone quality. We present a numerical procedure adapted for this purpose which is based on the spectral finite element method (FEM). By using a combined Laplace-Fourier transform, the vibroacoustic problem may be transformed into the frequency-wavenumber domain in which, as radiation conditions may be exactly introduced in the infinite fluid halfspaces, only the heterogeneous solid layer needs to be analysed using FEM. Several numerical tests are presented showing very good performance of the proposed approach. We present some results to study the influence of the frequency on the first arriving signal velocity in (visco)elastic bone plate.

Dynamics and Stability of Rolling Viscoelastic Tires

DOE Office of Scientific and Technical Information (OSTI.GOV)

Potter, Trevor

2013-04-30

Current steady state rolling tire calculations often do not include treads because treads destroy the rotational symmetry of the tire. We describe two methodologies to compute time periodic solutions of a two-dimensional viscoelastic tire with treads: solving a minimization problem and solving a system of equations. We also expand on work by Oden and Lin on free spinning rolling elastic tires in which they disovered a hierachy of N-peak steady state standing wave solutions. In addition to discovering a two-dimensional hierarchy of standing wave solutions that includes their N-peak hiearchy, we consider the eects of viscoelasticity on the standing wavemore » solutions. Finally, a commonplace model of viscoelasticity used in our numerical experiments led to non-physical elastic energy growth for large tire speeds. We show that a viscoelastic model of Govindjee and Reese remedies the problem.« less

A Galerkin least squares approach to viscoelastic flow.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rao, Rekha R.; Schunk, Peter Randall

2015-10-01

A Galerkin/least-squares stabilization technique is applied to a discrete Elastic Viscous Stress Splitting formulation of for viscoelastic flow. From this, a possible viscoelastic stabilization method is proposed. This method is tested with the flow of an Oldroyd-B fluid past a rigid cylinder, where it is found to produce inaccurate drag coefficients. Furthermore, it fails for relatively low Weissenberg number indicating it is not suited for use as a general algorithm. In addition, a decoupled approach is used as a way separating the constitutive equation from the rest of the system. A Pressure Poisson equation is used when the velocity andmore » pressure are sought to be decoupled, but this fails to produce a solution when inflow/outflow boundaries are considered. However, a coupled pressure-velocity equation with a decoupled constitutive equation is successful for the flow past a rigid cylinder and seems to be suitable as a general-use algorithm.« less

Viscoelastic analysis of a dental metal-ceramic system

NASA Astrophysics Data System (ADS)

Özüpek, Şebnem; Ünlü, Utku Cemal

2012-11-01

Porcelain-fused-to-metal (PFM) restorations used in prosthetic dentistry contain thermal stresses which develop during the cooling phase after firing. These thermal stresses coupled with the stresses produced by mechanical loads may be the dominant reasons for failures in clinical situations. For an accurate calculation of these stresses, viscoelastic behavior of ceramics at high temperatures should not be ignored. In this study, the finite element technique is used to evaluate the effect of viscoelasticity on stress distributions of a three-point flexure test specimen, which is the current international standard, ISO 9693, to characterize the interfacial bond strength of metal-ceramic restorative systems. Results indicate that the probability of interfacial debonding due to normal tensile stress is higher than that due to shear stress. This conclusion suggests modification of ISO 9693 bond strength definition from one in terms of the shear stress only to that accounting for both normal and shear stresses.

A micro-macro constitutive model for finite-deformation viscoelasticity of elastomers with nonlinear viscosity

NASA Astrophysics Data System (ADS)

Zhou, Jianyou; Jiang, Liying; Khayat, Roger E.

2018-01-01

Elastomers are known to exhibit viscoelastic behavior under deformation, which is linked to the diffusion processes of the highly mobile and flexible polymer chains. Inspired by the theories of polymer dynamics, a micro-macro constitutive model is developed to study the viscoelastic behaviors and the relaxation process of elastomeric materials under large deformation, in which the material parameters all have a microscopic foundation or a microstructural justification. The proposed model incorporates the nonlinear material viscosity into the continuum finite-deformation viscoelasticity theories which represent the polymer networks of elastomers with an elastic ground network and a few viscous subnetworks. The developed modeling framework is capable of adopting most of strain energy density functions for hyperelastic materials and thermodynamics evolution laws of viscoelastic solids. The modeling capacity of the framework is outlined by comparing the simulation results with the experimental data of three commonly used elastomeric materials, namely, VHB4910, HNBR50 and carbon black (CB) filled elastomers. The comparison shows that the stress responses and some typical behaviors of filled and unfilled elastomers can be quantitatively predicted by the model with suitable strain energy density functions. Particularly, the strain-softening effect of elastomers could be explained by the deformation-dependent (nonlinear) viscosity of the polymer chains. The presented modeling framework is expected to be useful as a modeling platform for further study on the performance of different type of elastomeric materials.

Microstructure-sensitive Crystal Viscoelasticity for Ni-base Superalloys Targeting Long-term Creep-Fatigue Interaction Modeling

DOE Office of Scientific and Technical Information (OSTI.GOV)

Neu, Richard W

The aim of this project is to develop a microstructure-sensitive crystal viscoplasticity (CVP) model for single-crystal Ni-base superalloys to model the behavior of the material and components in the hot gas path sections of industrial gas turbines (IGT). Microstructure degradation associated with aging critical to predicting long-term creep-fatigue interactions will be embedded into the model through the γ' precipitate morphology evolution by coupling the coarsening drivers and kinetics into the constitutive equations of the CVP model. Model parameters will be determined using new experimental protocols that involve systematically artificially aging the alloy under different stress conditions to determine the relationshipmore » between the size and morphology g' precipitates on the creep and thermomechanical fatigue response.« less

Instability behaviour of cosmic gravito-coupled correlative complex bi-fluidic admixture

NASA Astrophysics Data System (ADS)

Das, Papari; Karmakar, Pralay Kumar

2017-10-01

The gravitational instability of an unbounded infinitely extended composite gravitating cloud system composed of gravito-coupled neutral gaseous fluid (NGF) and dark matter fluid (DMF) is theoretically investigated in a classical framework. It is based on a spatially-flat geometry approximation (1D, sheet-like, boundless) at the backdrop that the radius of curvature of the gravito-confined bi-fluidic-boundary is much larger than all the hydro-characteristic scale lengths of interest. The relevant collective correlative dynamics, via the lowest-order mnemonic viscoelasticity, is mooted. We apply a standard formalism of normal mode analysis to yield a unique brand of generalized quadratic dispersion relation having variable multi-parametric coefficients dependent on the diversified equilibrium properties. It is parametrically seen that the DMF flow speed and the DMF viscoelasticity introduce stabilizing effects against the composite cloud collapse. The instability physiognomies, as specialized extreme corollaries, are in good accord with the previously reported predictions. The analysis may be widely useful to see the gravito-thermally coupled wave dynamics leading to the formation of large-scale hierarchical non-homologous structures in dark-matter-dominated dwarf galaxies.

Nonlinear Viscoelastic Mechanics of Cross-linked Rubbers

NASA Technical Reports Server (NTRS)

Freed, Alan D.; Leonov, Arkady I.; Gray, Hugh R. (Technical Monitor)

2002-01-01

The paper develops a general theory for finite rubber viscoelasticity, and specifies it in the form, convenient for solving problems important for rubber, tire and space industries. Based on the quasi-linear approach of non-equilibrium thermodynamics, a general nonlinear theory has been developed for arbitrary nonisothermal deformations of viscoelastic solids. In this theory, the constitutive equations are presented as the sum of known equilibrium (rubber elastic) and non-equilibrium (liquid polymer viscoelastic) terms. These equations are then simplified using several modeling arguments. Stability constraints for the proposed constitutive equations are also discussed. It is shown that only strong ellipticity criteria are applicable for assessing stability of the equations governing viscoelastic solids.

A Thermodynamic Theory of Solid Viscoelasticity. Part II:; Nonlinear Thermo-viscoelasticity

NASA Technical Reports Server (NTRS)

Freed, Alan D.; Leonov, Arkady I.; Gray, Hugh R. (Technical Monitor)

2002-01-01

This paper, second in the series of three papers, develops a general, nonlinear, non-isothermal, compressible theory for finite rubber viscoelasticity and specifies it in a form convenient for solving problems important to the rubber, tire, automobile, and air-space industries, among others. Based on the quasi-linear approach of non-equilibrium thermodynamics, a general nonlinear theory of differential type has been developed for arbitrary non-isothermal deformations of viscoelastic solids. In this theory, the constitutive equations were presented as the sum of a rubber elastic (equilibrium) and a liquid type viscoelastic (non-equilibrium) terms. These equations have then been simplified using several modeling and simplicity arguments.

Fractional viscoelasticity of soft elastomers and auxetic foams

NASA Astrophysics Data System (ADS)

Solheim, Hannah; Stanisauskis, Eugenia; Miles, Paul; Oates, William

2018-03-01

Dielectric elastomers are commonly implemented in adaptive structures due to their unique capabilities for real time control of a structure's shape, stiffness, and damping. These active polymers are often used in applications where actuator control or dynamic tunability are important, making an accurate understanding of the viscoelastic behavior critical. This challenge is complicated as these elastomers often operate over a broad range of deformation rates. Whereas research has demonstrated success in applying a nonlinear viscoelastic constitutive model to characterize the behavior of Very High Bond (VHB) 4910, robust predictions of the viscoelastic response over the entire range of time scales is still a significant challenge. An alternative formulation for viscoelastic modeling using fractional order calculus has shown significant improvement in predictive capabilities. While fractional calculus has been explored theoretically in the field of linear viscoelasticity, limited experimental validation and statistical evaluation of the underlying phenomena have been considered. In the present study, predictions across several orders of magnitude in deformation rates are validated against data using a single set of model parameters. Moreover, we illustrate the fractional order is material dependent by running complementary experiments and parameter estimation on the elastomer VHB 4949 as well as an auxetic foam. All results are statistically validated using Bayesian uncertainty methods to obtain posterior densities for the fractional order as well as the hyperelastic parameters.

Pulsatile flow of non-Newtonian blood fluid inside stenosed arteries: Investigating the effects of viscoelastic and elastic walls, arteriosclerosis, and polycythemia diseases.

PubMed

Nejad, A Abbas; Talebi, Z; Cheraghali, D; Shahbani-Zahiri, A; Norouzi, M

2018-02-01

In this study, the interaction of pulsatile blood flow with the viscoelastic walls of the axisymmetric artery is numerically investigated for different severities of stenosis. The geometry of artery is modeled by an axisymmetric cylindrical tube with a symmetric stenosis in a two-dimensional case. The effects of stenosis severity on the axial velocity profile, pressure distribution, streamlines, wall shear stress, and wall radial displacement for the viscoelastic artery are also compared to the elastics artery. Furthermore, the effects of atherosclerosis and polycythemia diseases on the hemodynamics and the mechanical behavior of arterial walls are investigated. The pulsatile flow of non-Newtonian blood is simulated inside the viscoelastic artery using the COMSOL Multiphysics software (version 5) and by employing the fluid-structure interaction (FSI) method and the arbitrary Lagrangian-Eulerian (ALE) method. Moreover, finite element method (FEM) is used to solve the governing equations on the unstructured grids. For modeling the non-Newtonian blood fluid and the viscoelastic arterial wall, the modified Casson model, and generalized Maxwell model are used, respectively. According to the results, with stenosis severity increasing from 25% to 75% at the time of maximum volumetric flow rate, the maximum value of axial velocity and its gradient increase 7.9 and 19.6 times, and the maximum wall shear stress of viscoelastic wall increases 24.2 times in the constriction zone. With the progression of the atherosclerosis disease (fivefold growth of arterial elastic modulus), the wall radial displacement of viscoelastic arterial walls decreases nearly 40%. In this study, axial velocity profile, pressure distribution, streamlines, wall radial displacement, and wall shear stress were examined for different percentages of stenosis (25%, 50%, and 75%). The atherosclerosis disease was investigated by the fivefold growth of viscoelastic arterial elastic modulus and polycythemia disease was examined by the 21-fold increase in the yield stress of the blood fluid. Furthermore, the comparison of results between the elastic and viscoelastic arterial walls shows that the wall radial displacement for viscoelastic artery is lower than that for the elastic artery as much as 21.7% for the severe stenosis of 75%. Copyright © 2017 Elsevier B.V. All rights reserved.

Probing viscoelastic surfaces with bimodal tapping-mode atomic force microscopy: Underlying physics and observables for a standard linear solid model.

PubMed

Solares, Santiago D

2014-01-01

This paper presents computational simulations of single-mode and bimodal atomic force microscopy (AFM) with particular focus on the viscoelastic interactions occurring during tip-sample impact. The surface is modeled by using a standard linear solid model, which is the simplest system that can reproduce creep compliance and stress relaxation, which are fundamental behaviors exhibited by viscoelastic surfaces. The relaxation of the surface in combination with the complexities of bimodal tip-sample impacts gives rise to unique dynamic behaviors that have important consequences with regards to the acquisition of quantitative relationships between the sample properties and the AFM observables. The physics of the tip-sample interactions and its effect on the observables are illustrated and discussed, and a brief research outlook on viscoelasticity measurement with intermittent-contact AFM is provided.

Improvement of sound insulation performance of double-glazed windows by using viscoelastic connectors

NASA Astrophysics Data System (ADS)

Takahashi, D.; Sawaki, S.; Mu, R.-L.

2016-06-01

A new method for improving the sound insulation performance of double-glazed windows is proposed. This technique uses viscoelastic materials as connectors between the two glass panels to ensure that the appropriate spacing is maintained. An analytical model that makes it possible to discuss the effects of spacing, contact area, and viscoelastic properties of the connectors on the performance in terms of sound insulation is developed. The validity of the model is verified by comparing its results with measured data. The numerical experiments using this analytical model showed the importance of the ability of the connectors to achieve the appropriate spacing and their viscoelastic properties, both of which are necessary for improving the sound insulation performance. In addition, it was shown that the most effective factor is damping: the stronger the damping, the more the insulation performance increases.

Lattice Boltzmann model for three-phase viscoelastic fluid flow

NASA Astrophysics Data System (ADS)

Xie, Chiyu; Lei, Wenhai; Wang, Moran

2018-02-01

A lattice Boltzmann (LB) framework is developed for simulation of three-phase viscoelastic fluid flows in complex geometries. This model is based on a Rothman-Keller type model for immiscible multiphase flows which ensures mass conservation of each component in porous media even for a high density ratio. To account for the viscoelastic effects, the Maxwell constitutive relation is correctly introduced into the momentum equation, which leads to a modified lattice Boltzmann evolution equation for Maxwell fluids by removing the normal but excess viscous term. Our simulation tests indicate that this excess viscous term may induce significant errors. After three benchmark cases, the displacement processes of oil by dispersed polymer are studied as a typical example of three-phase viscoelastic fluid flow. The results show that increasing either the polymer intrinsic viscosity or the elastic modulus will enhance the oil recovery.

Viscoelastic shear properties of human vocal fold mucosa: theoretical characterization based on constitutive modeling.

PubMed

Chan, R W; Titze, I R

2000-01-01

The viscoelastic shear properties of human vocal fold mucosa (cover) were previously measured as a function of frequency [Chan and Titze, J. Acoust. Soc. Am. 106, 2008-2021 (1999)], but data were obtained only in a frequency range of 0.01-15 Hz, an order of magnitude below typical frequencies of vocal fold oscillation (on the order of 100 Hz). This study represents an attempt to extrapolate the data to higher frequencies based on two viscoelastic theories, (1) a quasilinear viscoelastic theory widely used for the constitutive modeling of the viscoelastic properties of biological tissues [Fung, Biomechanics (Springer-Verlag, New York, 1993), pp. 277-292], and (2) a molecular (statistical network) theory commonly used for the rheological modeling of polymeric materials [Zhu et al., J. Biomech. 24, 1007-1018 (1991)]. Analytical expressions of elastic and viscous shear moduli, dynamic viscosity, and damping ratio based on the two theories with specific model parameters were applied to curve-fit the empirical data. Results showed that the theoretical predictions matched the empirical data reasonably well, allowing for parametric descriptions of the data and their extrapolations to frequencies of phonation.

A new mixed subgrid-scale model for large eddy simulation of turbulent drag-reducing flows of viscoelastic fluids

NASA Astrophysics Data System (ADS)

Li, Feng-Chen; Wang, Lu; Cai, Wei-Hua

2015-07-01

A mixed subgrid-scale (SGS) model based on coherent structures and temporal approximate deconvolution (MCT) is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation (LES) of turbulent drag-reducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence (FHIT) with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation (DNS) results. Compared with the LES results using the temporal approximate deconvolution model (TADM) for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number. For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives. Project supported by the China Postdoctoral Science Foundation (Grant No. 2011M500652), the National Natural Science Foundation of China (Grant Nos. 51276046 and 51206033), and the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20112302110020).

Viscoelastic flow past mono- and bidisperse random arrays of cylinders: flow resistance, topology and normal stress distribution.

PubMed

De, S; Kuipers, J A M; Peters, E A J F; Padding, J T

2017-12-13

We investigate creeping viscoelastic fluid flow through two-dimensional porous media consisting of random arrangements of monodisperse and bidisperse cylinders, using our finite volume-immersed boundary method introduced in S. De, et al., J. Non-Newtonian Fluid Mech., 2016, 232, 67-76. The viscoelastic fluid is modeled with a FENE-P model. The simulations show an increased flow resistance with increase in flow rate, even though the bulk response of the fluid to shear flow is shear thinning. We show that if the square root of the permeability is chosen as the characteristic length scale in the determination of the dimensionless Deborah number (De), then all flow resistance curves collapse to a single master curve, irrespective of the pore geometry. Our study reveals how viscoelastic stresses and flow topologies (rotation, shear and extension) are distributed through the porous media, and how they evolve with increasing De. We correlate the local viscoelastic first normal stress differences with the local flow topology and show that the largest normal stress differences are located in shear flow dominated regions and not in extensional flow dominated regions at higher viscoelasticity. The study shows that normal stress differences in shear flow regions may play a crucial role in the increase of flow resistance for viscoelastic flow through such porous media.

Application of a Numerical Inverse Laplace Integration Method to Surface Loading on a Viscoelastic Compressible Earth Model

NASA Astrophysics Data System (ADS)

Tanaka, Yoshiyuki; Klemann, Volker; Okuno, Jun'ichi

2009-09-01

Normal mode approaches for calculating viscoelastic responses of self-gravitating and compressible spherical earth models have an intrinsic problem of determining the roots of the secular equation and the associated residues in the Laplace domain. To bypass this problem, a method based on numerical inverse Laplace integration was developed by T anaka et al. (2006, 2007) for computations of viscoelastic deformation caused by an internal dislocation. The advantage of this approach is that the root-finding problem is avoided without imposing additional constraints on the governing equations and earth models. In this study, we apply the same algorithm to computations of viscoelastic responses to a surface load and show that the results obtained by this approach agree well with those obtained by a time-domain approach that does not need determinations of the normal modes in the Laplace domain. Using the elastic earth model PREM and a convex viscosity profile, we calculate viscoelastic load Love numbers ( h, l, k) for compressible and incompressible models. Comparisons between the results show that effects due to compressibility are consistent with results obtained by previous studies and that the rate differences between the two models total 10-40%. This will serve as an independent method to confirm results obtained by time-domain approaches and will usefully increase the reliability when modeling postglacial rebound.

Mechanisms governing the visco-elastic responses of living cells assessed by foam and tensegrity models.

PubMed

Cañadas, P; Laurent, V M; Chabrand, P; Isabey, D; Wendling-Mansuy, S

2003-11-01

The visco-elastic properties of living cells, measured to date by various authors, vary considerably, depending on the experimental methods and/or on the theoretical models used. In the present study, two mechanisms thought to be involved in cellular visco-elastic responses were analysed, based on the idea that the cytoskeleton plays a fundamental role in cellular mechanical responses. For this purpose, the predictions of an open unit-cell model and a 30-element visco-elastic tensegrity model were tested, taking into consideration similar properties of the constitutive F-actin. The quantitative predictions of the time constant and viscosity modulus obtained by both models were compared with previously published experimental data obtained from living cells. The small viscosity modulus values (10(0)-10(3) Pa x s) predicted by the tensegrity model may reflect the combined contributions of the spatially rearranged constitutive filaments and the internal tension to the overall cytoskeleton response to external loading. In contrast, the high viscosity modulus values (10(3)-10(5) Pa x s) predicted by the unit-cell model may rather reflect the mechanical response of the cytoskeleton to the bending of the constitutive filaments and/or to the deformation of internal components. The present results suggest the existence of a close link between the overall visco-elastic response of micromanipulated cells and the underlying architecture.

Study of non-linear deformation of vocal folds in simulations of human phonation

NASA Astrophysics Data System (ADS)

Saurabh, Shakti; Bodony, Daniel

2014-11-01

Direct numerical simulation is performed on a two-dimensional compressible, viscous fluid interacting with a non-linear, viscoelastic solid as a model for the generation of the human voice. The vocal fold (VF) tissues are modeled as multi-layered with varying stiffness in each layer and using a finite-strain Standard Linear Solid (SLS) constitutive model implemented in a quadratic finite element code and coupled to a high-order compressible Navier-Stokes solver through a boundary-fitted fluid-solid interface. The large non-linear mesh deformation is handled using an elliptic/poisson smoothening technique. Supra-glottal flow shows asymmetry in the flow, which in turn has a coupling effect on the motion of the VF. The fully compressible simulations gives direct insight into the sound produced as pressure distributions and the vocal fold deformation helps study the unsteady vortical flow resulting from the fluid-structure interaction along the full phonation cycle. Supported by the National Science Foundation (CAREER Award Number 1150439).

Nanoscale Viscoelasticity of Extracellular Matrix Proteins in Soft Tissues: a Multiscale Approach

PubMed Central

Miri, Amir K.; Heris, Hossein K.; Mongeau, Luc; Javid, Farhad

2013-01-01

We propose that the bulk viscoelasticity of soft tissues results from two length-scale-dependent mechanisms: the time-dependent response of extracellular matrix proteins (ECM) at the nanometer scale and the biophysical interactions between the ECM solid structure and interstitial fluid at the micrometer scale. The latter was modeled using the poroelasticity theory with an assumption of free motion of the interstitial fluid within the porous ECM structure. Following a recent study (Heris, H.K., Miri, A.K., Tripathy, U., Barthelat, F., Mongeau, L., 2013. Journal of the Mechanical Behavior of Biomedical Materials), atomic force microscopy was used to perform creep loading and 50-nm sinusoidal oscillations on porcine vocal folds. The proposed model was calibrated by a finite element model to accurately predict the nanoscale viscoelastic moduli of ECM. A linear correlation was observed between the in-depth distribution of the viscoelastic moduli and that of hyaluronic acids in the vocal fold tissue. We conclude that hyaluronic acids may regulate the vocal fold viscoelasticity at nanoscale. The proposed methodology offers a characterization tool for biomaterials used in vocal fold augmentations. PMID:24317493

Effects of viscoelasticity on drop impact and spreading on a solid surface

NASA Astrophysics Data System (ADS)

Izbassarov, Daulet; Muradoglu, Metin

2016-06-01

The effects of viscoelasticity on drop impact and spreading on a flat solid surface are studied computationally using a finite-difference-front-tracking method. The finitely extensible nonlinear elastic-Chilcott-Rallison model is used to account for the fluid viscoelasticity. It is found that viscoelasticity favors advancement of contact line during the spreading phase, leading to a slight increase in the maximum spreading, in agreement with experimental observations [Huh, Jung, Seo, and Lee, Microfluid. Nanofluid. 18, 1221 (2015), 10.1007/s10404-014-1518-4]. However, in contrast with the well-known antirebound effects of polymeric additives, the viscoelasticity is found to enhance the tendency of the drop rebound in the receding phase. These results suggest that the antirebound effects are mainly due to the polymer-induced modification of wetting properties of the substrate rather than the change in the material properties of the drop fluid. A model is proposed to test this hypothesis. It is found that the model results in good qualitative agreement with the experimental observations and the antirebound behavior can be captured by the modification of surface wetting properties in the receding phase.

Finite element formulation of viscoelastic sandwich beams using fractional derivative operators

NASA Astrophysics Data System (ADS)

Galucio, A. C.; Deü, J.-F.; Ohayon, R.

This paper presents a finite element formulation for transient dynamic analysis of sandwich beams with embedded viscoelastic material using fractional derivative constitutive equations. The sandwich configuration is composed of a viscoelastic core (based on Timoshenko theory) sandwiched between elastic faces (based on Euler-Bernoulli assumptions). The viscoelastic model used to describe the behavior of the core is a four-parameter fractional derivative model. Concerning the parameter identification, a strategy to estimate the fractional order of the time derivative and the relaxation time is outlined. Curve-fitting aspects are focused, showing a good agreement with experimental data. In order to implement the viscoelastic model into the finite element formulation, the Grünwald definition of the fractional operator is employed. To solve the equation of motion, a direct time integration method based on the implicit Newmark scheme is used. One of the particularities of the proposed algorithm lies in the storage of displacement history only, reducing considerably the numerical efforts related to the non-locality of fractional operators. After validations, numerical applications are presented in order to analyze truncation effects (fading memory phenomena) and solution convergence aspects.

Single Tracking Location Acoustic Radiation Force Impulse Viscoelasticity Estimation (STL-VE): A Method for Measuring Tissue Viscoelastic Parameters

PubMed Central

Langdon, Jonathan H; Elegbe, Etana; McAleavey, Stephen A

2015-01-01

Single Tracking Location (STL) Shear wave Elasticity Imaging (SWEI) is a method for detecting elastic differences between tissues. It has the advantage of intrinsic speckle bias suppression compared to Multiple Tracking Location (MTL) variants of SWEI. However, the assumption of a linear model leads to an overestimation of the shear modulus in viscoelastic media. A new reconstruction technique denoted Single Tracking Location Viscosity Estimation (STL-VE) is introduced to correct for this overestimation. This technique utilizes the same raw data generated in STL-SWEI imaging. Here, the STL-VE technique is developed by way of a Maximum Likelihood Estimation (MLE) for general viscoelastic materials. The method is then implemented for the particular case of the Kelvin-Voigt Model. Using simulation data, the STL-VE technique is demonstrated and the performance of the estimator is characterized. Finally, the STL-VE method is used to estimate the viscoelastic parameters of ex-vivo bovine liver. We find good agreement between the STL-VE results and the simulation parameters as well as between the liver shear wave data and the modeled data fit. PMID:26168170

Linear and Nonlinear Viscoelastic Modeling of Aorta and Carotid Pressure-Area Dynamics under in vivo and ex vivo Conditions

PubMed Central

Valdez-Jasso, Daniela; Bia, Daniel; Zócalo, Yanina; Armentano, Ricardo L.; Haider, Mansoor A.; Olufsen, Mette S.

2013-01-01

A better understanding of the biomechanical properties of the arterial wall provides important insight into arterial vascular biology under normal (healthy) and pathological conditions. This insight has potential to improve tracking of disease progression and to aid in vascular graft design and implementation. In this study, we use linear and nonlinear viscoelastic models to predict biomechanical properties of the thoracic descending aorta and the carotid artery under ex vivo and in vivo conditions in ovine and human arteries. Models analyzed include a four-parameter (linear) Kelvin viscoelastic model and two five-parameter nonlinear viscoelastic models (an arctangent and a sigmoid model) that relate changes in arterial blood pressure to the vessel cross-sectional area (via estimation of vessel strain). These models were developed using the framework of Quasilinear Viscoelasticity (QLV) theory and were validated using measurements from the thoracic descending aorta and the carotid artery obtained from human and ovine arteries. In vivo measurements were obtained from ten ovine aortas and ten human carotid arteries. Ex vivo measurements (from both locations) were made in eleven male Merino sheep. Biomechanical properties were obtained through constrained estimation of model parameters. To further investigate the parameter estimates we computed standard errors and confidence intervals and we used analysis of variance to compare results within and between groups. Overall, our results indicate that optimal model selection depends on the arterial type. Results showed that for the thoracic descending aorta (under both experimental conditions) the best predictions were obtained with the nonlinear sigmoid model, while under healthy physiological pressure loading the carotid arteries nonlinear stiffening with increasing pressure is negligible, and consequently, the linear (Kelvin) viscoelastic model better describes the pressure-area dynamics in this vessel. Results comparing biomechanical properties show that the Kelvin and sigmoid models were able to predict the zero-pressure vessel radius; that under ex vivo conditions vessels are more rigid, and comparatively, that the carotid artery is stiffer than the thoracic descending aorta; and that the viscoelastic gain and relaxation parameters do not differ significantly between vessels or experimental conditions. In conclusion, our study demonstrates that the proposed models can predict pressure-area dynamics and that model parameters can be extracted for further interpretation of biomechanical properties. PMID:21203846

Experimental study and constitutive modeling of the viscoelastic mechanical properties of the human prolapsed vaginal tissue.

PubMed

Peña, Estefania; Calvo, B; Martínez, M A; Martins, P; Mascarenhas, T; Jorge, R M N; Ferreira, A; Doblaré, M

2010-02-01

In this paper, the viscoelastic mechanical properties of vaginal tissue are investigated. Using previous results of the authors on the mechanical properties of biological soft tissues and newly experimental data from uniaxial tension tests, a new model for the viscoelastic mechanical properties of the human vaginal tissue is proposed. The structural model seems to be sufficiently accurate to guarantee its application to prediction of reliable stress distributions, and is suitable for finite element computations. The obtained results may be helpful in the design of surgical procedures with autologous tissue or prostheses.

DEM code-based modeling of energy accumulation and release in structurally heterogeneous rock masses

NASA Astrophysics Data System (ADS)

Lavrikov, S. V.; Revuzhenko, A. F.

2015-10-01

Based on discrete element method, the authors model loading of a physical specimen to describe its capacity to accumulate and release elastic energy. The specimen is modeled as a packing of particles with viscoelastic coupling and friction. The external elastic boundary of the packing is represented by particles connected by elastic springs. The latter means introduction of an additional special potential of interaction between the boundary particles, that exercises effect even when there is no direct contact between the particles. On the whole, the model specimen represents an element of a medium capable of accumulation of deformation energy in the form of internal stresses. The data of the numerical modeling of the physical specimen compression and the laboratory testing results show good qualitative consistency.

Hemodynamics in a Pediatric Ascending Aorta Using a Viscoelastic Pediatric Blood Model

PubMed Central

Good, Bryan C.; Deutsch, Steven; Manning, Keefe B.

2015-01-01

Congenital heart disease is the leading cause of infant death in the United States with over 36,000 newborns affected each year. Despite this growing problem there are few mechanical circulatory support devices designed specifically for pediatric and neonate patients. Previous research has been done investigating pediatric ventricular assist devices (PVADs) assuming blood to be a Newtonian fluid in computational fluid dynamics (CFD) simulations, ignoring its viscoelastic and shear-thinning properties. In contrast to adult VADs, PVADs may be more susceptible to hemolysis and thrombosis due to altered flow into the aorta, and therefore, a more accurate blood model should be used. A CFD solver that incorporates a modified Oldroyd-B model designed specifically for pediatric blood is used to investigate important hemodynamic parameters in a pediatric aortic model under pulsatile flow conditions. These results are compared to Newtonian blood simulations at three physiological pediatric hematocrits. Minor differences are seen in both velocity and WSS during early stages of the cardiac systole between the Newtonian and viscoelastic models. During diastole, significant differences are seen in the velocities in the descending aorta (up to 12%) and in the aortic branches (up to 30%) between the two models. Additionally, peak wall shear stress (WSS) differences are seen between the models throughout the cardiac cycle. At the onset of diastole, peak WSS differences of 43% are seen between the Newtonian and viscoelastic model and between the 20 and 60% hematocrit viscoelastic models at peak systole of 41%. PMID:26159560

Generalization of exponential based hyperelastic to hyper-viscoelastic model for investigation of mechanical behavior of rate dependent materials.

PubMed

Narooei, K; Arman, M

2018-03-01

In this research, the exponential stretched based hyperelastic strain energy was generalized to the hyper-viscoelastic model using the heredity integral of deformation history to take into account the strain rate effects on the mechanical behavior of materials. The heredity integral was approximated by the approach of Goh et al. to determine the model parameters and the same estimation was used for constitutive modeling. To present the ability of the proposed hyper-viscoelastic model, the stress-strain response of the thermoplastic elastomer gel tissue at different strain rates from 0.001 to 100/s was studied. In addition to better agreement between the current model and experimental data in comparison to the extended Mooney-Rivlin hyper-viscoelastic model, a stable material behavior was predicted for pure shear and balance biaxial deformation modes. To present the engineering application of current model, the Kolsky bars impact test of gel tissue was simulated and the effects of specimen size and inertia on the uniform deformation were investigated. As the mechanical response of polyurea was provided over wide strain rates of 0.0016-6500/s, the current model was applied to fit the experimental data. The results were shown more accuracy could be expected from the current research than the extended Ogden hyper-viscoelastic model. In the final verification example, the pig skin experimental data was used to determine parameters of the hyper-viscoelastic model. Subsequently, a specimen of pig skin at different strain rates was loaded to a fixed strain and the change of stress with time (stress relaxation) was obtained. The stress relaxation results were revealed the peak stress increases by applied strain rate until the saturated loading rate and the equilibrium stress with magnitude of 0.281MPa could be reached. Copyright © 2017 Elsevier Ltd. All rights reserved.

Dynamics of multiple viscoelastic carbon nanotube based nanocomposites with axial magnetic field

DOE Office of Scientific and Technical Information (OSTI.GOV)

Karličić, Danilo; Cajić, Milan; Murmu, Tony

2014-06-21

Nanocomposites and magnetic field effects on nanostructures have received great attention in recent years. A large amount of research work was focused on developing the proper theoretical framework for describing many physical effects appearing in structures on nanoscale level. Great step in this direction was successful application of nonlocal continuum field theory of Eringen. In the present paper, the free transverse vibration analysis is carried out for the system composed of multiple single walled carbon nanotubes (MSWCNT) embedded in a polymer matrix and under the influence of an axial magnetic field. Equivalent nonlocal model of MSWCNT is adopted as viscoelasticallymore » coupled multi-nanobeam system (MNBS) under the influence of longitudinal magnetic field. Governing equations of motion are derived using the Newton second low and nonlocal Rayleigh beam theory, which take into account small-scale effects, the effect of nanobeam angular acceleration, internal damping and Maxwell relation. Explicit expressions for complex natural frequency are derived based on the method of separation of variables and trigonometric method for the “Clamped-Chain” system. In addition, an analytical method is proposed in order to obtain asymptotic damped natural frequency and the critical damping ratio, which are independent of boundary conditions and a number of nanobeams in MNBS. The validity of obtained results is confirmed by comparing the results obtained for complex frequencies via trigonometric method with the results obtained by using numerical methods. The influence of the longitudinal magnetic field on the free vibration response of viscoelastically coupled MNBS is discussed in detail. In addition, numerical results are presented to point out the effects of the nonlocal parameter, internal damping, and parameters of viscoelastic medium on complex natural frequencies of the system. The results demonstrate the efficiency of the suggested methodology to find the closed form solutions for the free vibration response of multiple nanostructure systems under the influence of magnetic field.« less

Determining the Viscosity Coefficient for Viscoelastic Wave Propagation in Rock Bars

NASA Astrophysics Data System (ADS)

Niu, Leilei; Zhu, Wancheng; Li, Shaohua; Guan, Kai

2018-05-01

Rocks with microdefects exhibit viscoelastic behavior during stress wave propagation. The viscosity coefficient of the wave can be used to characterize the attenuation as the wave propagates in rock. In this study, a long artificial bar with a readily adjustable viscosity coefficient was fabricated to investigate stress wave attenuation. The viscoelastic behavior of the artificial bar under dynamic loading was investigated, and the initial viscoelastic coefficient was obtained based on the amplitude attenuation of the incident harmonic wave. A one-dimensional wave propagation program was compiled to reproduce the time history of the stress wave measured during the experiments, and the program was well fitted to the Kelvin-Voigt model. The attenuation and dispersion of the stress wave in long artificial viscoelastic bars were quantified to accurately determine the viscoelastic coefficient. Finally, the method used to determine the viscoelastic coefficient of a long artificial bar based on the experiments and numerical simulations was extended to determine the viscoelastic coefficient of a short rock bar. This study provides a new method of determining the viscosity coefficient of rock.

Wormhole Formation in RSRM Nozzle Joint Backfill

NASA Technical Reports Server (NTRS)

Stevens, J.

2000-01-01

The RSRM nozzle uses a barrier of RTV rubber upstream of the nozzle O-ring seals. Post flight inspection of the RSRM nozzle continues to reveal occurrence of "wormholes" into the RTV backfill. The term "wormholes", sometimes called "gas paths", indicates a gas flow path not caused by pre-existing voids, but by a little-understood internal failure mode of the material during motor operation. Fundamental understanding of the mechanics of the RSRM nozzle joints during motor operation, nonlinear viscoelastic characterization of the RTV backfill material, identification of the conditions that predispose the RTV to form wormholes, and screening of candidate replacement materials is being pursued by a joint effort between Thiokol Propulsion, NASA, and the Army Propulsion & Structures Directorate at Redstone Arsenal. The performance of the RTV backfill in the joint is controlled by the joint environment. Joint movement, which applies a tension and shear load on the material, coupled with the introduction of high pressure gas in combination create an environment that exceeds the capability of the material to withstand the wormhole effect. Little data exists to evaluate why the material fails under the modeled joint conditions, so an effort to characterize and evaluate the material under these conditions was undertaken. Viscoelastic property data from characterization testing will anchor structural analysis models. Data over a range of temperatures, environmental pressures, and strain rates was used to develop a nonlinear viscoelastic model to predict material performance, develop criteria for replacement materials, and quantify material properties influencing wormhole growth. Three joint simulation analogs were developed to analyze and validate joint thermal barrier (backfill) material performance. Two exploratory tests focus on detection of wormhole failure under specific motor operating conditions. A "validation" test system provides data to "validate" computer models and predictions. Finally, two candidate replacement materials are being screened and "validated" using the developed test systems.

Dynamic assessment of nonlinear typical section aeroviscoelastic systems using fractional derivative-based viscoelastic model

NASA Astrophysics Data System (ADS)

Sales, T. P.; Marques, Flávio D.; Pereira, Daniel A.; Rade, Domingos A.

2018-06-01

Nonlinear aeroelastic systems are prone to the appearance of limit cycle oscillations, bifurcations, and chaos. Such problems are of increasing concern in aircraft design since there is the need to control nonlinear instabilities and improve safety margins, at the same time as aircraft are subjected to increasingly critical operational conditions. On the other hand, in spite of the fact that viscoelastic materials have already been successfully used for the attenuation of undesired vibrations in several types of mechanical systems, a small number of research works have addressed the feasibility of exploring the viscoelastic effect to improve the behavior of nonlinear aeroelastic systems. In this context, the objective of this work is to assess the influence of viscoelastic materials on the aeroelastic features of a three-degrees-of-freedom typical section with hardening structural nonlinearities. The equations of motion are derived accounting for the presence of viscoelastic materials introduced in the resilient elements associated to each degree-of-freedom. A constitutive law based on fractional derivatives is adopted, which allows the modeling of temperature-dependent viscoelastic behavior in time and frequency domains. The unsteady aerodynamic loading is calculated based on the classical linear potential theory for arbitrary airfoil motion. The aeroelastic behavior is investigated through time domain simulations, and subsequent frequency transformations, from which bifurcations are identified from diagrams of limit cycle oscillations amplitudes versus airspeed. The influence of the viscoelastic effect on the aeroelastic behavior, for different values of temperature, is also investigated. The numerical simulations show that viscoelastic damping can increase the flutter speed and reduce the amplitudes of limit cycle oscillations. These results prove the potential that viscoelastic materials have to increase aircraft components safety margins regarding aeroelastic stability.

Lamb Wave Dispersion Ultrasound Vibrometry (LDUV) Method for Quantifying Mechanical Properties of Viscoelastic Solids

PubMed Central

Nenadic, Ivan Z.; Urban, Matthew W.; Mitchell, Scott A.; Greenleaf, James F.

2011-01-01

Diastolic dysfunction is the inability of the left ventricle to supply sufficient stroke volumes under normal physiological conditions and is often accompanied by stiffening of the left-ventricular myocardium. A noninvasive technique capable of quantifying viscoelasticity of the myocardium would be beneficial in clinical settings. Our group has been investigating the use of Shearwave Dispersion Ultrasound Vibrometry (SDUV), a noninvasive ultrasound based method for quantifying viscoelasticity of soft tissues. The primary motive of this study is the design and testing of viscoelastic materials suitable for validation of the Lamb wave Dispersion Ultrasound Vibrometry (LDUV), an SDUV-based technique for measuring viscoelasticity of tissues with plate-like geometry. We report the results of quantifying viscoelasticity of urethane rubber and gelatin samples using LDUV and an embedded sphere method. The LDUV method was used to excite antisymmetric Lamb waves and measure the dispersion in urethane rubber and gelatin plates. An antisymmetric Lamb wave model was fitted to the wave speed dispersion data to estimate elasticity and viscosity of the materials. A finite element model of a viscoelastic plate submerged in water was used to study the appropriateness of the Lamb wave dispersion equations. An embedded sphere method was used as an independent measurement of the viscoelasticity of the urethane rubber and gelatin. The FEM dispersion data were in excellent agreement with the theoretical predictions. Viscoelasticity of the urethane rubber and gelatin obtained using the LDUV and embedded sphere methods agreed within one standard deviation. LDUV studies on excised porcine myocardium sample were performed to investigate the feasibility of the approach in preparation for open-chest in vivo studies. The results suggest that the LDUV technique can be used to quantify mechanical properties of soft tissues with a plate-like geometry. PMID:21403186

Lamb wave dispersion ultrasound vibrometry (LDUV) method for quantifying mechanical properties of viscoelastic solids.

PubMed

Nenadic, Ivan Z; Urban, Matthew W; Mitchell, Scott A; Greenleaf, James F

2011-04-07

Diastolic dysfunction is the inability of the left ventricle to supply sufficient stroke volumes under normal physiological conditions and is often accompanied by stiffening of the left-ventricular myocardium. A noninvasive technique capable of quantifying viscoelasticity of the myocardium would be beneficial in clinical settings. Our group has been investigating the use of shear wave dispersion ultrasound vibrometry (SDUV), a noninvasive ultrasound-based method for quantifying viscoelasticity of soft tissues. The primary motive of this study is the design and testing of viscoelastic materials suitable for validation of the Lamb wave dispersion ultrasound vibrometry (LDUV), an SDUV-based technique for measuring viscoelasticity of tissues with plate-like geometry. We report the results of quantifying viscoelasticity of urethane rubber and gelatin samples using LDUV and an embedded sphere method. The LDUV method was used to excite antisymmetric Lamb waves and measure the dispersion in urethane rubber and gelatin plates. An antisymmetric Lamb wave model was fitted to the wave speed dispersion data to estimate elasticity and viscosity of the materials. A finite element model of a viscoelastic plate submerged in water was used to study the appropriateness of the Lamb wave dispersion equations. An embedded sphere method was used as an independent measurement of the viscoelasticity of the urethane rubber and gelatin. The FEM dispersion data were in excellent agreement with the theoretical predictions. Viscoelasticity of the urethane rubber and gelatin obtained using the LDUV and embedded sphere methods agreed within one standard deviation. LDUV studies on excised porcine myocardium sample were performed to investigate the feasibility of the approach in preparation for open-chest in vivo studies. The results suggest that the LDUV technique can be used to quantify the mechanical properties of soft tissues with a plate-like geometry.

Spectral modelling of ice-induced wave decay: implementation of a new viscoelastic theory in WAVEWATCH III

NASA Astrophysics Data System (ADS)

Liu, Q.; Rogers, W. E.; Babanin, A. V.; Squire, V. A.; Mosig, J. E. M.; Li, J.; Guan, C.

2017-12-01

A new viscoelastic ice layer model is implemented in the third generation spectral wave model WAVEWATCH III to estimate the ice-induced, frequency-dependent wave attenuation rate. Two case studies are then conducted with this viscoelastic model: one is the hindcast of waves in the autumn Beaufort Sea, 2015, and the other is the modelling of wave fields in the Antarctic marginal ice zone (MIZ), 2012. It is demonstrated that the viscoelastic model is capable of reproducing the measured significant wave heights (Ηs) in these two different geophysical regions. The sensitivity of the simulated wave height on different source terms -- ice-induced decay Sice and other physical processes Sother such as wind input Sin, nonlinear four-wave interaction Snl -- is also investigated in this study. For the Antarctic MIZ experiment, Sother is found to be much less than Sice and thus contributes little to the simulated Hs. The trend of the wave height decay (dHs/dx) discovered recently -- saturating at large wave heights -- is well reproduced by the standalone linear viscoelastic model. The flattening of dHs/dx is most likely due to the only presence of longer waves, with the shorter waves having been already low-pass filtered. Nonetheless, Sother should not be disregarded within a more general modelling perspective as Sin and Snl is shown to be comparable or even much higher than Sice in the Beaufort Sea case.

Effect of Longitudinal Magnetic Field on Vibration Characteristics of Single-Walled Carbon Nanotubes in a Viscoelastic Medium

NASA Astrophysics Data System (ADS)

Zhang, D. P.; Lei, Y.; Shen, Z. B.

2017-12-01

The effect of longitudinal magnetic field on vibration response of a sing-walled carbon nanotube (SWCNT) embedded in viscoelastic medium is investigated. Based on nonlocal Euler-Bernoulli beam theory, Maxwell's relations, and Kelvin viscoelastic foundation model, the governing equations of motion for vibration analysis are established. The complex natural frequencies and corresponding mode shapes in closed form for the embedded SWCNT with arbitrary boundary conditions are obtained using transfer function method (TFM). The new analytical expressions for the complex natural frequencies are also derived for certain typical boundary conditions and Kelvin-Voigt model. Numerical results from the model are presented to show the effects of nonlocal parameter, viscoelastic parameter, boundary conditions, aspect ratio, and strength of the magnetic field on vibration characteristics for the embedded SWCNT in longitudinal magnetic field. The results demonstrate the efficiency of the proposed methods for vibration analysis of embedded SWCNTs under magnetic field.

A multilayer model of time dependent deformation following an earthquake on a strike-slip fault

NASA Technical Reports Server (NTRS)

Cohen, S. C.

1981-01-01

A multilayer model of the Earth to calculate finite element of time dependent deformation and stress following an earthquake on a strike slip fault is discussed. The model involves shear properties of an elastic upper lithosphere, a standard viscoelastic linear solid lower lithosphere, a Maxwell viscoelastic asthenosphere and an elastic mesosphere. Systematic variations of fault and layer depths and comparisons with simpler elastic lithosphere over viscoelastic asthenosphere calculations are analyzed. Both the creep of the lower lithosphere and astenosphere contribute to the postseismic deformation. The magnitude of the deformation is enhanced by a short distance between the bottom of the fault (slip zone) and the top of the creep region but is less sensitive to the thickness of the creeping layer. Postseismic restressing is increased as the lower lithosphere becomes more viscoelastic, but the tendency for the width of the restressed zone to growth with time is retarded.

Probing viscoelastic surfaces with bimodal tapping-mode atomic force microscopy: Underlying physics and observables for a standard linear solid model

PubMed Central

2014-01-01

Summary This paper presents computational simulations of single-mode and bimodal atomic force microscopy (AFM) with particular focus on the viscoelastic interactions occurring during tip–sample impact. The surface is modeled by using a standard linear solid model, which is the simplest system that can reproduce creep compliance and stress relaxation, which are fundamental behaviors exhibited by viscoelastic surfaces. The relaxation of the surface in combination with the complexities of bimodal tip–sample impacts gives rise to unique dynamic behaviors that have important consequences with regards to the acquisition of quantitative relationships between the sample properties and the AFM observables. The physics of the tip–sample interactions and its effect on the observables are illustrated and discussed, and a brief research outlook on viscoelasticity measurement with intermittent-contact AFM is provided. PMID:25383277

Persistence length of collagen molecules based on nonlocal viscoelastic model.

PubMed

Ghavanloo, Esmaeal

2017-12-01

Persistence length is one of the most interesting properties of a molecular chain, which is used to describe the stiffness of a molecule. The experimentally measured values of the persistence length of the collagen molecule are widely scattered from 14 to 180 nm. Therefore, an alternative approach is highly desirable to predict the persistence length of a molecule and also to explain the experimental results. In this paper, a nonlocal viscoelastic model is developed to obtain the persistence length of the collagen molecules in solvent. A new explicit formula is proposed for the persistence length of the molecule with the consideration of the small-scale effect, viscoelastic properties of the molecule, loading frequency, and viscosity of the solvent. The presented model indicates that there exists a range of molecule lengths in which the persistence length strongly depends on the frequency and spatial mode of applied loads, small-scale effect, and viscoelastic properties of the collagen.

Efficiency of perfectly matched layers for seismic wave modeling in second-order viscoelastic equations

NASA Astrophysics Data System (ADS)

Ping, Ping; Zhang, Yu; Xu, Yixian; Chu, Risheng

2016-12-01

In order to improve the perfectly matched layer (PML) efficiency in viscoelastic media, we first propose a split multi-axial PML (M-PML) and an unsplit convolutional PML (C-PML) in the second-order viscoelastic wave equations with the displacement as the only unknown. The advantage of these formulations is that it is easy and efficient to revise the existing codes of the second-order spectral element method (SEM) or finite-element method (FEM) with absorbing boundaries in a uniform equation, as well as more economical than the auxiliary differential equations PML. Three models which are easily suffered from late time instabilities are considered to validate our approaches. Through comparison the M-PML with C-PML efficiency of absorption and stability for long time simulation, it can be concluded that: (1) for an isotropic viscoelastic medium with high Poisson's ratio, the C-PML will be a sufficient choice for long time simulation because of its weak reflections and superior stability; (2) unlike the M-PML with high-order damping profile, the M-PML with second-order damping profile loses its stability in long time simulation for an isotropic viscoelastic medium; (3) in an anisotropic viscoelastic medium, the C-PML suffers from instabilities, while the M-PML with second-order damping profile can be a better choice for its superior stability and more acceptable weak reflections than the M-PML with high-order damping profile. The comparative analysis of the developed methods offers meaningful significance for long time seismic wave modeling in second-order viscoelastic wave equations.

Viscoelastic Properties of Confluent MDCK II Cells Obtained from Force Cycle Experiments.

PubMed

Brückner, Bastian Rouven; Nöding, Helen; Janshoff, Andreas

2017-02-28

The local mechanical properties of cells are frequently probed by force indentation experiments carried out with an atomic force microscope. Application of common contact models provides a single parameter, the Young's modulus, to describe the elastic properties of cells. The viscoelastic response of cells, however, is generally measured in separate microrheological experiments that provide complex shear moduli as a function of time or frequency. Here, we present a straightforward way to obtain rheological properties of cells from regular force distance curves collected in typical force indentation measurements. The method allows us to record the stress-strain relationship as well as changes in the weak power law of the viscoelastic moduli. We derive an analytical function based on the elastic-viscoelastic correspondence principle applied to Hertzian contact mechanics to model both indentation and retraction curves. Rheological properties are described by standard viscoelastic models and the paradigmatic weak power law found to interpret the viscoelastic properties of living cells best. We compare our method with atomic force microscopy-based active oscillatory microrheology and show that the method to determine the power law coefficient is robust against drift and largely independent of the indentation depth and indenter geometry. Cells were subject to Cytochalasin D treatment to provoke a drastic change in the power law coefficient and to demonstrate the feasibility of the approach to capture rheological changes extremely fast and precisely. The method is easily adaptable to different indenter geometries and acquires viscoelastic data with high spatiotemporal resolution. Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.

Finite element implementation of a thermo-damage-viscoelastic constitutive model for hydroxyl-terminated polybutadiene composite propellant

NASA Astrophysics Data System (ADS)

Xu, Jinsheng; Han, Long; Zheng, Jian; Chen, Xiong; Zhou, Changsheng

2017-11-01

A thermo-damage-viscoelastic model for hydroxyl-terminated polybutadiene (HTPB) composite propellant with consideration for the effect of temperature was implemented in ABAQUS. The damage evolution law of the model has the same form as the crack growth equation for viscoelastic materials, and only a single damage variable S is considered. The HTPB propellant was considered as an isotropic material, and the deviatoric and volumetric strain-stress relations are decoupled and described by the bulk and shear relaxation moduli, respectively. The stress update equations were expressed by the principal stresses σ_{ii}R and the rotation tensor M, the Jacobian matrix in the global coordinate system J_{ijkl} was obtained according to the fourth-order tensor transformation rules. Two models having complex stress states were used to verify the accuracy of the constitutive model. The test results showed good agreement with the strain responses of characteristic points measured by a contactless optical deformation test system, which illustrates that the thermo-damage-viscoelastic model perform well at describing the mechanical properties of an HTPB propellant.

An anomalous subdiffusion model for calcium spark in cardiac myocytes

NASA Astrophysics Data System (ADS)

Tan, Wenchang; Fu, Chaoqi; Fu, Ceji; Xie, Wenjun; Cheng, Heping

2007-10-01

The elementary events of excitation-contraction coupling in heart muscle are Ca2+ sparks, which arise from ryanodine receptors in the sarcoplasmic reticulum (SR). Here, an anomalous subdiffusion model is developed to explore Ca2+ spark formation in cardiac myocytes. Numerical simulations reproduce the brightness, the time course, and spatial size of a typical cardiac Ca2+ spark. It is suggested that the diffusion of Ca2+ spark in the cytoplasm may no longer obey Fickian second law, but the anomalous space subdiffusion. The physical reason is perhaps due to the effects of the electric field of the calcium ions and the viscoelasticity of the cytoplasm and its complex structures.

Acoustics of marine sediment under compaction: binary grain-size model and viscoelastic extension of Biot's theory.

PubMed

Leurer, Klaus C; Brown, Colin

2008-04-01

This paper presents a model of acoustic wave propagation in unconsolidated marine sediment, including compaction, using a concept of a simplified sediment structure, modeled as a binary grain-size sphere pack. Compressional- and shear-wave velocities and attenuation follow from a combination of Biot's model, used as the general framework, and two viscoelastic extensions resulting in complex grain and frame moduli, respectively. An effective-grain model accounts for the viscoelasticity arising from local fluid flow in expandable clay minerals in clay-bearing sediments. A viscoelastic-contact model describes local fluid flow at the grain contacts. Porosity, density, and the structural Biot parameters (permeability, pore size, structure factor) as a function of pressure follow from the binary model, so that the remaining input parameters to the acoustic model consist solely of the mass fractions and the known mechanical properties of each constituent (e.g., carbonates, sand, clay, and expandable clay) of the sediment, effective pressure, or depth, and the environmental parameters (water depth, salinity, temperature). Velocity and attenuation as a function of pressure from the model are in good agreement with data on coarse- and fine-grained unconsolidated marine sediments.

Laser Interferometer Measurements of the Viscoelastic Properties of Tectorial Membrane Mutants

NASA Astrophysics Data System (ADS)

Jones, Gareth; Russell, Ian; Lukashkin, Andrei

2011-11-01

The visco-elastic properties of the tectorial membrane (TM) can be determined by measuring the propagation velocity of travelling waves over a range of frequencies. This study presents a new method using laser interferometry and compares the TM's material properties (sheer storage modulus, G' and viscosity, η) at basal and apical locations in wild-type mice and basal locations of three mutant groups (TectaY1870C/+, Tectb-/- and Otoa-/-). The G' and η values calculated for the wild-type mice are similar to estimates derived using other methods whereas the mutant groups all exhibit slower wave propagation velocities and reduced longitudinal coupling.

Structural signature of a brittle-to-ductile transition in self-assembled networks.

PubMed

Ramos, Laurence; Laperrousaz, Arnaud; Dieudonné, Philippe; Ligoure, Christian

2011-09-30

We study the nonlinear rheology of a novel class of transient networks, made of surfactant micelles of tunable morphology reversibly linked by block copolymers. We couple rheology and time-resolved structural measurements, using synchrotron radiation, to characterize the highly nonlinear viscoelastic regime. We propose the fluctuations of the degree of alignment of the micelles under shear as a probe to identify a fracture process. We show a clear signature of a brittle-to-ductile transition in transient gels, as the morphology of the micelles varies, and provide a parallel between the fracture of solids and the fracture under shear of viscoelastic fluids.

Random walks of colloidal probes in viscoelastic materials

NASA Astrophysics Data System (ADS)

Khan, Manas; Mason, Thomas G.

2014-04-01

To overcome limitations of using a single fixed time step in random walk simulations, such as those that rely on the classic Wiener approach, we have developed an algorithm for exploring random walks based on random temporal steps that are uniformly distributed in logarithmic time. This improvement enables us to generate random-walk trajectories of probe particles that span a highly extended dynamic range in time, thereby facilitating the exploration of probe motion in soft viscoelastic materials. By combining this faster approach with a Maxwell-Voigt model (MVM) of linear viscoelasticity, based on a slowly diffusing harmonically bound Brownian particle, we rapidly create trajectories of spherical probes in soft viscoelastic materials over more than 12 orders of magnitude in time. Appropriate windowing of these trajectories over different time intervals demonstrates that random walk for the MVM is neither self-similar nor self-affine, even if the viscoelastic material is isotropic. We extend this approach to spatially anisotropic viscoelastic materials, using binning to calculate the anisotropic mean square displacements and creep compliances along different orthogonal directions. The elimination of a fixed time step in simulations of random processes, including random walks, opens up interesting possibilities for modeling dynamics and response over a highly extended temporal dynamic range.

A nonlocal strain gradient model for dynamic deformation of orthotropic viscoelastic graphene sheets under time harmonic thermal load

NASA Astrophysics Data System (ADS)

Radwan, Ahmed F.; Sobhy, Mohammed

2018-06-01

This work presents a nonlocal strain gradient theory for the dynamic deformation response of a single-layered graphene sheet (SLGS) on a viscoelastic foundation and subjected to a time harmonic thermal load for various boundary conditions. Material of graphene sheets is presumed to be orthotropic and viscoelastic. The viscoelastic foundation is modeled as Kelvin-Voigt's pattern. Based on the two-unknown plate theory, the motion equations are obtained from the dynamic version of the virtual work principle. The nonlocal strain gradient theory is established from Eringen nonlocal and strain gradient theories, therefore, it contains two material scale parameters, which are nonlocal parameter and gradient coefficient. These scale parameters have two different effects on the graphene sheets. The obtained deflection is compared with that predicted in the literature. Additional numerical examples are introduced to illustrate the influences of the two length scale coefficients and other parameters on the dynamic deformation of the viscoelastic graphene sheets.

Discrete quasi-linear viscoelastic damping analysis of connective tissues, and the biomechanics of stretching.

PubMed

Babaei, Behzad; Velasquez-Mao, Aaron J; Thomopoulos, Stavros; Elson, Elliot L; Abramowitch, Steven D; Genin, Guy M

2017-05-01

The time- and frequency-dependent properties of connective tissue define their physiological function, but are notoriously difficult to characterize. Well-established tools such as linear viscoelasticity and the Fung quasi-linear viscoelastic (QLV) model impose forms on responses that can mask true tissue behavior. Here, we applied a more general discrete quasi-linear viscoelastic (DQLV) model to identify the static and dynamic time- and frequency-dependent behavior of rabbit medial collateral ligaments. Unlike the Fung QLV approach, the DQLV approach revealed that energy dissipation is elevated at a loading period of ∼10s. The fitting algorithm was applied to the entire loading history on each specimen, enabling accurate estimation of the material's viscoelastic relaxation spectrum from data gathered from transient rather than only steady states. The application of the DQLV method to cyclically loading regimens has broad applicability for the characterization of biological tissues, and the results suggest a mechanistic basis for the stretching regimens most favored by athletic trainers. Copyright © 2017 Elsevier Ltd. All rights reserved.

An inverse method for determining the spatially resolved properties of viscoelastic-viscoplastic three-dimensional printed materials.

PubMed

Chen, X; Ashcroft, I A; Wildman, R D; Tuck, C J

2015-11-08

A method using experimental nanoindentation and inverse finite-element analysis (FEA) has been developed that enables the spatial variation of material constitutive properties to be accurately determined. The method was used to measure property variation in a three-dimensional printed (3DP) polymeric material. The accuracy of the method is dependent on the applicability of the constitutive model used in the inverse FEA, hence four potential material models: viscoelastic, viscoelastic-viscoplastic, nonlinear viscoelastic and nonlinear viscoelastic-viscoplastic were evaluated, with the latter enabling the best fit to experimental data. Significant changes in material properties were seen in the depth direction of the 3DP sample, which could be linked to the degree of cross-linking within the material, a feature inherent in a UV-cured layer-by-layer construction method. It is proposed that the method is a powerful tool in the analysis of manufacturing processes with potential spatial property variation that will also enable the accurate prediction of final manufactured part performance.

Discrete quasi-linear viscoelastic damping analysis of connective tissues, and the biomechanics of stretching

PubMed Central

Babaei, Behzad; Velasquez-Mao, Aaron J.; Thomopoulos, Stavros; Elson, Elliot L.; Abramowitch, Steven D.; Genin, Guy M.

2017-01-01

The time- and frequency-dependent properties of connective tissue define their physiological function, but are notoriously difficult to characterize. Well-established tools such as linear viscoelasticity and the Fung quasi-linear viscoelastic (QLV) model impose forms on responses that can mask true tissue behavior. Here, we applied a more general discrete quasi-linear viscoelastic (DQLV) model to identify the static and dynamic time- and frequency-dependent behavior of rabbit medial collateral ligaments. Unlike the Fung QLV approach, the DQLV approach revealed that energy dissipation is elevated at a loading period of ~10 seconds. The fitting algorithm was applied to the entire loading history on each specimen, enabling accurate estimation of the material's viscoelastic relaxation spectrum from data gathered from transient rather than only steady states. The application of the DQLV method to cyclically loading regimens has broad applicability for the characterization of biological tissues, and the results suggest a mechanistic basis for the stretching regimens most favored by athletic trainers. PMID:28088071

Ultrasonic measurements of breast viscoelasticity.

PubMed

Sridhar, Mallika; Insana, Michael F

2007-12-01

In vivo measurements of the viscoelastic properties of breast tissue are described. Ultrasonic echo frames were recorded from volunteers at 5 fps while applying a uniaxial compressive force (1-20 N) within a 1 s ramp time and holding the force constant for up to 200 s. A time series of strain images was formed from the echo data, spatially averaged viscous creep curves were computed, and viscoelastic strain parameters were estimated by fitting creep curves to a second-order Voigt model. The useful strain bandwidth from this quasi-static ramp stimulus was 10(-2) < or = omega < or = 10(0) rad/s (0.0016-0.16 Hz). The stress-strain curves for normal glandular tissues are linear when the surface force applied is between 2 and 5 N. In this range, the creep response was characteristic of biphasic viscoelastic polymers, settling to a constant strain (arrheodictic) after 100 s. The average model-based retardance time constants for the viscoelastic response were 3.2 +/- 0.8 and 42.0 +/- 28 s. Also, the viscoelastic strain amplitude was approximately equal to that of the elastic strain. Above 5 N of applied force, however, the response of glandular tissue became increasingly nonlinear and rheodictic, i.e., tissue creep never reached a plateau. Contrasting in vivo breast measurements with those in gelatin hydrogels, preliminary ideas regarding the mechanisms for viscoelastic contrast are emerging.

Ultrasonic measurements of breast viscoelasticity

PubMed Central

Sridhar, Mallika; Insana, Michael F.

2009-01-01

In vivo measurements of the viscoelastic properties of breast tissue are described. Ultrasonic echo frames were recorded from volunteers at 5 fps while applying a uniaxial compressive force (1–20 N) within a 1 s ramp time and holding the force constant for up to 200 s. A time series of strain images was formed from the echo data, spatially averaged viscous creep curves were computed, and viscoelastic strain parameters were estimated by fitting creep curves to a second-order Voigt model. The useful strain bandwidth from this quasi-static ramp stimulus was 10−2 ≤ ω ≤ 100 rad/s (0.0016–0.16 Hz). The stress-strain curves for normal glandular tissues are linear when the surface force applied is between 2 and 5 N. In this range, the creep response was characteristic of biphasic viscoelastic polymers, settling to a constant strain (arrheodictic) after 100 s. The average model-based retardance time constants for the viscoelastic response were 3.2±0.8 and 42.0±28 s. Also, the viscoelastic strain amplitude was approximately equal to that of the elastic strain. Above 5 N of applied force, however, the response of glandular tissue became increasingly nonlinear and rheodictic, i.e., tissue creep never reached a plateau. Contrasting in vivo breast measurements with those in gelatin hydrogels, preliminary ideas regarding the mechanisms for viscoelastic contrast are emerging. PMID:18196803

Thickness-shear mode quartz crystal resonators in viscoelastic fluid media

NASA Astrophysics Data System (ADS)

Arnau, A.; Jiménez, Y.; Sogorb, T.

2000-10-01

An extended Butterworth-Van Dyke (EBVD) model to characterize a thickness-shear mode quartz crystal resonator in a semi-infinite viscoelastic medium is derived by means of analysis of the lumped elements model described by Cernosek et al. [R. W. Cernosek, S. J. Martin, A. R. Hillman, and H. L. Bandey, IEEE Trans. Ultrason. Ferroelectr. Freq. Control 45, 1399 (1998)]. The EBVD model parameters are related to the viscoelastic properties of the medium. A capacitance added to the motional branch of the EBVD model has to be included when the elastic properties of the fluid are considered. From this model, an explicit expression for the frequency shift of a quartz crystal sensor in viscoelastic media is obtained. By combining the expressions for shifts in the motional series resonant frequency and in the motional resistance, a simple equation that relates only one unknown (the loss factor of the fluid) to those measurable quantities, and two simple explicit expressions for determining the viscoelastic properties of semi-infinite fluid media have been derived. The proposed expression for the parameter Δf/ΔR is compared with the corresponding ratio obtained with data computed from the complete admittance model. Relative errors below 4.5%, 3%, and 1.2% (for the ratios of the load surface mechanical impedance to the quartz shear characteristic impedance of 0.3, 0.25, and 0.1, respectively), are obtained in the range of the cases analyzed. Experimental data from the literature are used to validate the model.

Viscoelastic shear lag model to predict the micromechanical behavior of tendon under dynamic tensile loading.

PubMed

Wu, Jiayu; Yuan, Hong; Li, Longyuan; Fan, Kunjie; Qian, Shanguang; Li, Bing

2018-01-21

Owing to its viscoelastic nature, tendon exhibits stress rate-dependent breaking and stiffness function. A Kelvin-Voigt viscoelastic shear lag model is proposed to illustrate the micromechanical behavior of the tendon under dynamic tensile conditions. Theoretical closed-form expressions are derived to predict the deformation and stress transfer between fibrils and interfibrillar matrix while tendon is dynamically stretched. The results from the analytical solutions demonstrate that how the fibril overlap length and fibril volume fraction affect the stress transfer and mechanical properties of tendon. We find that the viscoelastic property of interfibrillar matrix mainly results in collagen fibril failure under fast loading rate or creep rupture of tendon. However, discontinuous fibril model and hierarchical structure of tendon ensure relative sliding under slow loading rate, helping dissipate energy and protecting fibril from damage, which may be a key reason why regularly staggering alignment microstructure is widely selected in nature. According to the growth, injury, healing and healed process of tendon observed by many researchers, the conclusions presented in this paper agrees well with the experimental findings. Additionally, the emphasis of this paper is on micromechanical behavior of tendon, whereas this analytical viscoelastic shear lag model can be equally applicable to other soft or hard tissues, owning the similar microstructure. Copyright © 2017 Elsevier Ltd. All rights reserved.

A theoretical and experimental technique to measure fracture properties in viscoelastic solids

NASA Astrophysics Data System (ADS)

Freitas, Felipe Araujo Colares De

Prediction of crack growth in engineering structures is necessary for better analysis and design. However, this prediction becomes quite complex for certain materials in which the fracture behavior is both rate and path dependent. Asphaltic materials used in pavements have that intrinsic complexity in their behavior. A lot of research effort has been devoted to better understanding viscoelastic behavior and fracture in such materials. This dissertation presents a further refinement of an experimental test setup, which is significantly different from standard testing protocols, to measure viscoelastic and fracture properties of nonlinear viscoelastic solids, such as asphaltic materials. The results presented herein are primarily for experiments with asphalt, but the test procedure can be used for other viscoelastic materials as well. Even though the test is designed as a fracture test, experiments on the investigated materials have uncovered very complex phenomena prior to fracture. Viscoelasticity and micromechanics are used to explain some of the physical phenomena observed in the tests. The material behavior prior to fracture includes both viscoelastic behavior and a necking effect, which is further discussed in the appendix of the present study. The dissertation outlines a theoretical model for the prediction of tractions ahead of the crack tip. The major contribution herein lies in the development of the experimental procedure for evaluating the material parameters necessary for deploying the model in the prediction of ductile crack growth. Finally, predictions of crack growth in a double cantilever beam specimens and asphalt concrete samples are presented in order to demonstrate the power of this approach for predicting crack growth in viscoelastic media.

Modeling and simulation of viscoelastic biological particles' 3D manipulation using atomic force microscopy

NASA Astrophysics Data System (ADS)

Korayem, M. H.; Habibi Sooha, Y.; Rastegar, Z.

2018-05-01

Manipulation of the biological particles by atomic force microscopy is used to transfer these particles inside body's cells, diagnosis and destruction of the cancer cells and drug delivery to damaged cells. According to the impossibility of simultaneous observation of this process, the importance of modeling and simulation can be realized. The contact of the tip with biological particle is important during manipulation, therefore, the first step of the modeling is choosing appropriate contact model. Most of the studies about contact between atomic force microscopy and biological particles, consider the biological particle as an elastic material. This is not an appropriate assumption because biological cells are basically soft and this assumption ignores loading history. In this paper, elastic and viscoelastic JKR theories were used in modeling and simulation of the 3D manipulation for three modes of tip-particle sliding, particle-substrate sliding and particle-substrate rolling. Results showed that critical force and time in motion modes (sliding and rolling) for two elastic and viscoelastic states are very close but these magnitudes were lower in the viscoelastic state. Then, three friction models, Coulomb, LuGre and HK, were used for tip-particle sliding mode in the first phase of manipulation to make results closer to reality. In both Coulomb and LuGre models, critical force and time are very close for elastic and viscoelastic states but in general critical force and time prediction of HK model was higher than LuGre and the LuGre model itself had higher prediction than Coulomb.

Performance evaluation of traveling wave ultrasonic motor based on a model with visco-elastic friction layer on stator.

PubMed

Qu, Jianjun; Sun, Fengyan; Zhao, Chunsheng

2006-12-01

A new visco-elastic contact model of traveling wave ultrasonic motor (TWUSM) is proposed. In this model, the rotor is assumed to be rigid body and the friction material on stator teeth surface to be visco-elastic body. Both load characteristics of TWUSM, such as rotation speed, torque and efficiency, and effects of interface parameters between stator and rotor on output characteristic of TWUSM can be calculated and simulated numerically by using MATLAB method based on this model. This model is compared with that one of compliant slider and rigid stator. The results show that this model can obtain bigger stall torque. The simulated results are compared with test results, and found that their load characteristics have good agreement.

The viscoelastic characterization of polymer materials exposed to the low-Earth orbit environment

NASA Technical Reports Server (NTRS)

Strganac, Thomas; Letton, Alan

1992-01-01

Recent accomplishments in our research efforts have included the successful measurement of the thermal mechanical properties of polymer materials exposed to the low-earth orbit environment. In particular, viscoelastic properties were recorded using the Rheometrics Solids Analyzer (RSA 2). Dynamic moduli (E', the storage component of the elastic modulus, and E'', the loss component of the elastic modulus) were recorded over three decades of frequency (0.1 to 100 rad/sec) for temperatures ranging from -150 to 150 C. Although this temperature range extends beyond the typical use range of the materials, measurements in this region are necessary in the development of complete viscoelastic constitutive models. The experimental results were used to provide the stress relaxation and creep compliance performance characteristics through viscoelastic correspondence principles. Our results quantify the differences between exposed and control polymer specimens. The characterization is specifically designed to elucidate a constitutive model that accurately predicts the change in behavior of these materials due to exposure. The constitutive model for viscoelastic behavior reflects the level of strain, the rate of strain, and the history of strain as well as the thermal history of the material.

The psychomechanics of simulated sound sources: Material properties of impacted bars

NASA Astrophysics Data System (ADS)

McAdams, Stephen; Chaigne, Antoine; Roussarie, Vincent

2004-03-01

Sound can convey information about the materials composing an object that are often not directly available to the visual system. Material and geometric properties of synthesized impacted bars with a tube resonator were varied, their perceptual structure was inferred from multidimensional scaling of dissimilarity judgments, and the psychophysical relations between the two were quantified. Constant cross-section bars varying in mass density and viscoelastic damping coefficient were synthesized with a physical model in experiment 1. A two-dimensional perceptual space resulted, and the dimensions were correlated with the mechanical parameters after applying a power-law transformation. Variable cross-section bars varying in length and viscoelastic damping coefficient were synthesized in experiment 2 with two sets of lengths creating high- and low-pitched bars. In the low-pitched bars, there was a coupling between the bar and the resonator that modified the decay characteristics. Perceptual dimensions again corresponded to the mechanical parameters. A set of potential temporal, spectral, and spectrotemporal correlates of the auditory representation were derived from the signal. The dimensions related to mass density and bar length were correlated with the frequency of the lowest partial and are related to pitch perception. The correlate most likely to represent the viscoelastic damping coefficient across all three stimulus sets is a linear combination of a decay constant derived from the temporal envelope and the spectral center of gravity derived from a cochlear representation of the signal. These results attest to the perceptual salience of energy-loss phenomena in sound source behavior.

The fractional derivative Kelvin–Voigt model of viscoelasticity with and without volumetric relaxation

NASA Astrophysics Data System (ADS)

Rossikhin, Yu A.; Shitikova, M. V.

2018-04-01

The fractional derivative Kelvin–Voigt model of viscoelasticity involving the time-dependent Poisson’s operator has been studied not only for the case of a time-independent bulk modulus, but also when the volumetric relaxation is taken into account. It has been shown that such a model could describe the features of auxetic materials.

Crustal deformation in great California earthquake cycles

NASA Technical Reports Server (NTRS)

Li, Victor C.; Rice, James R.

1986-01-01

Periodic crustal deformation associated with repeated strike slip earthquakes is computed for the following model: A depth L (less than or similiar to H) extending downward from the Earth's surface at a transform boundary between uniform elastic lithospheric plates of thickness H is locked between earthquakes. It slips an amount consistent with remote plate velocity V sub pl after each lapse of earthquake cycle time T sub cy. Lower portions of the fault zone at the boundary slip continuously so as to maintain constant resistive shear stress. The plates are coupled at their base to a Maxwellian viscoelastic asthenosphere through which steady deep seated mantle motions, compatible with plate velocity, are transmitted to the surface plates. The coupling is described approximately through a generalized Elsasser model. It is argued that the model gives a more realistic physical description of tectonic loading, including the time dependence of deep slip and crustal stress build up throughout the earthquake cycle, than do simpler kinematic models in which loading is represented as imposed uniform dislocation slip on the fault below the locked zone.

An immersed-boundary method for flow–structure interaction in biological systems with application to phonation

PubMed Central

Luo, Haoxiang; Mittal, Rajat; Zheng, Xudong; Bielamowicz, Steven A.; Walsh, Raymond J.; Hahn, James K.

2008-01-01

A new numerical approach for modeling a class of flow–structure interaction problems typically encountered in biological systems is presented. In this approach, a previously developed, sharp-interface, immersed-boundary method for incompressible flows is used to model the fluid flow and a new, sharp-interface Cartesian grid, immersed boundary method is devised to solve the equations of linear viscoelasticity that governs the solid. The two solvers are coupled to model flow–structure interaction. This coupled solver has the advantage of simple grid generation and efficient computation on simple, single-block structured grids. The accuracy of the solid-mechanics solver is examined by applying it to a canonical problem. The solution methodology is then applied to the problem of laryngeal aerodynamics and vocal fold vibration during human phonation. This includes a three-dimensional eigen analysis for a multi-layered vocal fold prototype as well as two-dimensional, flow-induced vocal fold vibration in a modeled larynx. Several salient features of the aerodynamics as well as vocal-fold dynamics are presented. PMID:19936017

Modeling and Simulation of Viscous Electro-Active Polymers

PubMed Central

Vogel, Franziska; Göktepe, Serdar; Steinmann, Paul; Kuhl, Ellen

2014-01-01

Electro-active materials are capable of undergoing large deformation when stimulated by an electric field. They can be divided into electronic and ionic electro-active polymers (EAPs) depending on their actuation mechanism based on their composition. We consider electronic EAPs, for which attractive Coulomb forces or local re-orientation of polar groups cause a bulk deformation. Many of these materials exhibit pronounced visco-elastic behavior. Here we show the development and implementation of a constitutive model, which captures the influence of the electric field on the visco-elastic response within a geometrically non-linear finite element framework. The electric field affects not only the equilibrium part of the strain energy function, but also the viscous part. To adopt the familiar additive split of the strain from the small strain setting, we formulate the governing equations in the logarithmic strain space and additively decompose the logarithmic strain into elastic and viscous parts. We show that the incorporation of the electric field in the viscous response significantly alters the relaxation and hysteresis behavior of the model. Our parametric study demonstrates that the model is sensitive to the choice of the electro-viscous coupling parameters. We simulate several actuator structures to illustrate the performance of the method in typical relaxation and creep scenarios. Our model could serve as a design tool for micro-electro-mechanical systems, microfluidic devices, and stimuli-responsive gels such as artificial skin, tactile displays, or artificial muscle. PMID:25267881

Comparison of different constitutive models to characterize the viscoelastic properties of human abdominal adipose tissue. A pilot study.

PubMed

Calvo-Gallego, Jose L; Domínguez, Jaime; Gómez Cía, Tomás; Gómez Ciriza, Gorka; Martínez-Reina, Javier

2018-04-01

Knowing the mechanical properties of human adipose tissue is key to simulate surgeries such as liposuction, mammoplasty and many plastic surgeries in which the subcutaneous fat is present. One of the most important surgeries, for its incidence, is the breast reconstruction surgery that follows a mastectomy. In this case, achieving a deformed shape similar to the healthy breast is crucial. The reconstruction is most commonly made using autologous tissue, taken from the patient's abdomen. The amount of autologous tissue and its mechanical properties have a strong influence on the shape of the reconstructed breast. In this work, the viscoelastic mechanical properties of the human adipose tissue have been studied. Uniaxial compression stress relaxation tests were performed in adipose tissue specimens extracted from the human abdomen. Two different viscoelastic models were used to fit to the experimental tests: a quasi-linear viscoelastic (QLV) model and an internal variables viscoelastic (IVV) model; each one with four different hyperelastic strain energy density functions to characterise the elastic response: a 5-terms polynomial function, a first order Ogden function, an isotropic Gasser-Ogden-Holzapfel function and a combination of a neoHookean and an exponential function. The IVV model with the Ogden function was the best combination to fit the experimental tests. The viscoelastic properties are not important in the simulation of the static deformed shape of the breast, but they are needed in a relaxation test performed under finite strain rate, particularly, to derive the long-term behaviour (as time tends to infinity), needed to estimate the static deformed shape of the breast. The so obtained stiffness was compared with previous results given in the literature for adipose tissue of different regions, which exhibited a wide dispersion. Copyright © 2018 Elsevier Ltd. All rights reserved.

An experimentally based nonlinear viscoelastic model of joint passive moment.

PubMed

Esteki, A; Mansour, J M

1996-04-01

Previous investigations have not converged on a generally accepted model of the dissipative part of joint passive moment. To provide a basis for developing a model, a series of measurements were performed to characterize the passive moment at the metacarpophalangeal joint of the index finger. Two measurement procedures were used, one in moment relaxation over a range of fixed joint angles and the other at a series of constant joint velocities. Fung's quasi-linear viscoelastic theory motivated the development of the passive moment model. Using this approach, it was not necessary to make restrictive assumptions regarding the viscoelastic behavior of the passive moment. The generality of the formulation allowed specific functions to be chosen based on experimental data rather than finding coefficients which attempted to fit a preselected model of the data. It was shown that a nonlinear viscoelastic model described the passive stiffness. No significant frictional effects were found. Of particular importance was the nonlinear behavior of the dissipative part of the passive moment which was modeled by joint speed raised to a power less than one. This result could explain the differing findings among previous investigations, and may have important implications for control of limb movement.

A Brief Review of Elasticity and Viscoelasticity

DTIC Science & Technology

2010-05-27

through electromagnetic or acoustic means. Creating a model that accurately describes these Rayleigh waves is key to modeling and understanding the...technology to be feasible, a mathematical model that describes the propagation of the acoustic wave from the stenosis to the chest wall will be necessary...viscoelastic model is simpler to use than poroelastic models but yields similar results for a wide range of soils and dynamic 30 loadings. In addition

Numerical design and test on an assembled structure of a bolted joint with viscoelastic damping

NASA Astrophysics Data System (ADS)

Hammami, Chaima; Balmes, Etienne; Guskov, Mikhail

2016-03-01

Mechanical assemblies are subjected to many dynamic loads and modifications are often needed to achieve acceptable vibration levels. While modifications on mass and stiffness are well mastered, damping modifications are still considered difficult to design. The paper presents a case study on the design of a bolted connection containing a viscoelastic damping layer. The notion of junction coupling level is introduced to ensure that sufficient energy is present in the joints to allow damping. Static performance is then addressed and it is shown that localization of metallic contact can be used to meet objectives, while allowing the presence of viscoelastic materials. Numerical prediction of damping then illustrates difficulties in optimizing for robustness. Modal test results of three configurations of an assembled structure, inspired by aeronautic fuselages, are then compared to analyze the performance of the design. While validity of the approach is confirmed, the effect of geometric imperfections is shown and stresses the need for robust design.

Modeling the Coupled Chemo-Thermo-Mechanical Behavior of Amorphous Polymer Networks.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zimmerman, Jonathan A.; Nguyen, Thao D.; Xiao, Rui

2015-02-01

Amorphous polymers exhibit a rich landscape of time-dependent behavior including viscoelasticity, structural relaxation, and viscoplasticity. These time-dependent mechanisms can be exploited to achieve shape-memory behavior, which allows the material to store a programmed deformed shape indefinitely and to recover entirely the undeformed shape in response to specific environmental stimulus. The shape-memory performance of amorphous polymers depends on the coordination of multiple physical mechanisms, and considerable opportunities exist to tailor the polymer structure and shape-memory programming procedure to achieve the desired performance. The goal of this project was to use a combination of theoretical, numerical and experimental methods to investigate themore » effect of shape memory programming, thermo-mechanical properties, and physical and environmental aging on the shape memory performance. Physical and environmental aging occurs during storage and through exposure to solvents, such as water, and can significantly alter the viscoelastic behavior and shape memory behavior of amorphous polymers. This project – executed primarily by Professor Thao Nguyen and Graduate Student Rui Xiao at Johns Hopkins University in support of a DOE/NNSA Presidential Early Career Award in Science and Engineering (PECASE) – developed a theoretical framework for chemothermo- mechanical behavior of amorphous polymers to model the effects of physical aging and solvent-induced environmental factors on their thermoviscoelastic behavior.« less

Viscoelastic flow modeling in the extrusion of a dough-like fluid

NASA Technical Reports Server (NTRS)

Dhanasekharan, M.; Kokini, J. L.; Janes, H. W. (Principal Investigator)

2000-01-01

This work attempts to investigate the effect of viscoelasticity and three-dimensional geometry in screw channels. The Phan-Thien Tanner (PTT) constitutive equation with simplified model parameters was solved in conjunction with the flow equations. Polyflow, a commercially available finite element code was used to solve the resulting nonlinear partial differential equations. The PTT model predicted one log scale lower pressure buildup compared to the equivalent Newtonian results. However, the velocity profile did not show significant changes for the chosen PTT model parameters. Past Researchers neglected viscoelastic effects and also the three dimensional nature of the flow in extruder channels. The results of this paper provide a starting point for further simulations using more realistic model parameters, which may enable the food engineer to more accurately scale-up and design extrusion processes.

Near-field postseismic deformation associated with the 1992 Landers and 1999 Hector Mine, California, earthquakes

USGS Publications Warehouse

Savage, J.C.; Svarc, J.L.; Prescott, W.H.

2003-01-01

After the Landers earthquake (Mw = 7.3, 1992.489) a linear array of 10 monuments extending about 30 km N50??E on either side of the earthquake rupture plus a nearby offtrend reference monument were surveyed frequently by GPS until 2003.2. The array also spans the rupture of the subsequent Hector Mine earthquake (Mw = 7.1, 1999.792 . The pre-Landers velocities of monuments in the array relative to interior North America were estimated from earlier trilateration and very long baseline interferometry measurements. Except at the reference monument, the post-Landers velocities of the individual monuments in the array relaxed to their preseismic values within 4 years. Following the Hector Mine earthquake the velocities of the monuments relaxed to steady rates within 1 year. Those steady rates for the east components are about equal to the pre-Landers rates as is the steady rate for the north component of the one monument east of the Hector Mine rupture. However, the steady rates for the north components of the 10 monuments west of the rupture are systematically ???10 mm yr1 larger than the pre-Landers rates. The relaxation to a steady rate is approximately exponential with decay times of 0.50 ?? 0.10 year following the Landers earthquake and 0.32 ?? 0.18 year following the Hector Mine earthquake. The postearthquake motions of the Landers array following the Landers earthquake are not well approximated by the viscoelastic-coupling model of Pollitz et al. [2000]. A similar viscoelastic-coupling model [Pollitz et al., 2001] is more successful in representing the deformation after the Hector Mine earthquake.

Pseudospectral modeling and dispersion analysis of Rayleigh waves in viscoelastic media

USGS Publications Warehouse

Zhang, K.; Luo, Y.; Xia, J.; Chen, C.

2011-01-01

Multichannel Analysis of Surface Waves (MASW) is one of the most widely used techniques in environmental and engineering geophysics to determine shear-wave velocities and dynamic properties, which is based on the elastic layered system theory. Wave propagation in the Earth, however, has been recognized as viscoelastic and the propagation of Rayleigh waves presents substantial differences in viscoelastic media as compared with elastic media. Therefore, it is necessary to carry out numerical simulation and dispersion analysis of Rayleigh waves in viscoelastic media to better understand Rayleigh-wave behaviors in the real world. We apply a pseudospectral method to the calculation of the spatial derivatives using a Chebyshev difference operator in the vertical direction and a Fourier difference operator in the horizontal direction based on the velocity-stress elastodynamic equations and relations of linear viscoelastic solids. This approach stretches the spatial discrete grid to have a minimum grid size near the free surface so that high accuracy and resolution are achieved at the free surface, which allows an effective incorporation of the free surface boundary conditions since the Chebyshev method is nonperiodic. We first use an elastic homogeneous half-space model to demonstrate the accuracy of the pseudospectral method comparing with the analytical solution, and verify the correctness of the numerical modeling results for a viscoelastic half-space comparing the phase velocities of Rayleigh wave between the theoretical values and the dispersive image generated by high-resolution linear Radon transform. We then simulate three types of two-layer models to analyze dispersive-energy characteristics for near-surface applications. Results demonstrate that the phase velocity of Rayleigh waves in viscoelastic media is relatively higher than in elastic media and the fundamental mode increases by 10-16% when the frequency is above 10. Hz due to the velocity dispersion of P and S waves. ?? 2011 Elsevier Ltd.

Constitutive Modeling of Porcine Liver in Indentation Using 3D Ultrasound Imaging

PubMed Central

Jordan, P.; Socrate, S.; Zickler, T.E.; Howe, R.D.

2009-01-01

In this work we present an inverse finite-element modeling framework for constitutive modeling and parameter estimation of soft tissues using full-field volumetric deformation data obtained from 3D ultrasound. The finite-element model is coupled to full-field visual measurements by regularization springs attached at nodal locations. The free ends of the springs are displaced according to the locally estimated tissue motion and the normalized potential energy stored in all springs serves as a measure of model-experiment agreement for material parameter optimization. We demonstrate good accuracy of estimated parameters and consistent convergence properties on synthetically generated data. We present constitutive model selection and parameter estimation for perfused porcine liver in indentation and demonstrate that a quasilinear viscoelastic model with shear modulus relaxation offers good model-experiment agreement in terms of indenter displacement (0.19 mm RMS error) and tissue displacement field (0.97 mm RMS error). PMID:19627823

Viscoelastic crustal deformation by magmatic intrusion: A case study in the Kutcharo caldera, eastern Hokkaido, Japan

NASA Astrophysics Data System (ADS)

Yamasaki, Tadashi; Kobayashi, Tomokazu; Wright, Tim J.; Fukahata, Yukitoshi

2018-01-01

Geodetic signals observed at volcanoes, particularly their temporal patterns, have required us to make the correlation between the surface displacement and magmatic process at depth in terms of viscoelastic crustal rheology. Here we use a parallelized 3-D finite element model to examine the response of the linear Maxwell viscoelastic crust and mantle to the inflation of a sill in order to show the characteristics of a long-term volcano deformation. In the model, an oblate-spheroidal sill is instantaneously or gradually inflated in a two-layered medium that consists of an elastic layer underlain by a viscoelastic layer. Our numerical experiments show that syn-inflation surface uplift is followed by post-inflation surface subsidence as the viscoelastic substrate relaxes. For gradual inflation events, the magnitude of inflation-induced uplift is reduced by the relaxation, through which the volume of a magma inferred by matching the prediction of an elastic model with observed surface uplift could be underestimated. For a given crustal viscosity, sill depth is the principal factor controlling subsidence caused by viscoelastic relaxation. The subsidence rate is highest when the inflation occurs at the boundary between the elastic and the viscoelastic layers. The mantle viscosity has an insignificant impact unless the depth of the inflation is greater than a half the crustal thickness. We apply the viscoelastic model to the interferometric synthetic aperture radar (InSAR) data in the Kutcharo caldera, eastern Hokkaido, Japan, where the surface has slowly subsided over a period of approximately three years following about a two-year period of inflation. The emplacement of a magmatic sill is constrained to occur at a depth of 4.5 km, which is significantly shallower than the geophysically imaged large-scale magma chamber. The geodetically detected deformation in the caldera reflects the small-scale emplacement of a magma that ascended from the deeper chamber, but not the inflation of the chamber itself. The observed ground displacement is controlled by a lower-crustal viscosity of 4 × 1017 Pa s, which is lower than that inferred from some studies of post-seismic deformation, perhaps due to higher temperatures beneath the active caldera. Our results suggest that geodetic signals observed during and following magmatic intrusions need to be revisited. Uzs‧ is the uplift at t‧ = Δt‧ for models with Δt‧ > 0.

Surface loading of a viscoelastic earth-I. General theory

NASA Astrophysics Data System (ADS)

Tromp, Jeroen; Mitrovica, Jerry X.

1999-06-01

We present a new normal-mode formalism for computing the response of an aspherical, self-gravitating, linear viscoelastic earth model to an arbitrary surface load. The formalism makes use of recent advances in the theory of the Earth's free oscillations, and is based upon an eigenfunction expansion methodology, rather than the tradi-tional Love-number approach to surface-loading problems. We introduce a surface-load representation theorem analogous to Betti's reciprocity relation in seismology. Taking advantage of this theorem and the biorthogonality of the viscoelastic modes, we determine the complete response to a surface load in the form of a Green's function. We also demonstrate that each viscoelastic mode has its own unique energy partitioning, which can be used to characterize it. In subsequent papers, we apply the theory to spherically symmetric and aspherical earth models.

Fractional Derivative Models for Ultrasonic Characterization of Polymer and Breast Tissue Viscoelasticity

PubMed Central

Coussot, Cecile; Kalyanam, Sureshkumar; Yapp, Rebecca; Insana, Michael F.

2009-01-01

The viscoelastic response of hydropolymers, which include glandular breast tissues, may be accurately characterized for some applications with as few as 3 rheological parameters by applying the Kelvin-Voigt fractional derivative (KVFD) modeling approach. We describe a technique for ultrasonic imaging of KVFD parameters in media undergoing unconfined, quasi-static, uniaxial compression. We analyze the KVFD parameter values in simulated and experimental echo data acquired from phantoms and show that the KVFD parameters may concisely characterize the viscoelastic properties of hydropolymers. We then interpret the KVFD parameter values for normal and cancerous breast tissues and hypothesize that this modeling approach may ultimately be applied to tumor differentiation. PMID:19406700

Surface waves on a soft viscoelastic layer produced by an oscillating microbubble.

PubMed

Tinguely, Marc; Hennessy, Matthew G; Pommella, Angelo; Matar, Omar K; Garbin, Valeria

2016-05-14

Ultrasound-driven bubbles can cause significant deformation of soft viscoelastic layers, for instance in surface cleaning and biomedical applications. The effect of the viscoelastic properties of a boundary on the bubble-boundary interaction has been explored only qualitatively, and remains poorly understood. We investigate the dynamic deformation of a viscoelastic layer induced by the volumetric oscillations of an ultrasound-driven microbubble. High-speed video microscopy is used to observe the deformation produced by a bubble oscillating at 17-20 kHz in contact with the surface of a hydrogel. The localised oscillating pressure applied by the bubble generates surface elastic (Rayleigh) waves on the gel, characterised by elliptical particle trajectories. The tilt angle of the elliptical trajectories varies with increasing distance from the bubble. Unexpectedly, the direction of rotation of the surface elements on the elliptical trajectories shifts from prograde to retrograde at a distance from the bubble that depends on the viscoelastic properties of the gel. To explain these behaviours, we develop a simple three-dimensional model for the deformation of a viscoelastic solid by a localised oscillating force. By using as input for the model the values of the shear modulus obtained from the propagation velocity of the Rayleigh waves, we find good qualitative agreement with the experimental observations.

Nonlinear viscoelastic characterization of bovine trabecular bone.

PubMed

Manda, Krishnagoud; Wallace, Robert J; Xie, Shuqiao; Levrero-Florencio, Francesc; Pankaj, Pankaj

2017-02-01

The time-independent elastic properties of trabecular bone have been extensively investigated, and several stiffness-density relations have been proposed. Although it is recognized that trabecular bone exhibits time-dependent mechanical behaviour, a property of viscoelastic materials, the characterization of this behaviour has received limited attention. The objective of the present study was to investigate the time-dependent behaviour of bovine trabecular bone through a series of compressive creep-recovery experiments and to identify its nonlinear constitutive viscoelastic material parameters. Uniaxial compressive creep and recovery experiments at multiple loads were performed on cylindrical bovine trabecular bone samples ([Formula: see text]). Creep response was found to be significant and always comprised of recoverable and irrecoverable strains, even at low stress/strain levels. This response was also found to vary nonlinearly with applied stress. A systematic methodology was developed to separate recoverable (nonlinear viscoelastic) and irrecoverable (permanent) strains from the total experimental strain response. We found that Schapery's nonlinear viscoelastic constitutive model describes the viscoelastic response of the trabecular bone, and parameters associated with this model were estimated from the multiple load creep-recovery (MLCR) experiments. Nonlinear viscoelastic recovery compliance was found to have a decreasing and then increasing trend with increasing stress level, indicating possible stiffening and softening behaviour of trabecular bone due to creep. The obtained parameters from MLCR tests, expressed as second-order polynomial functions of stress, showed a similar trend for all the samples, and also demonstrate stiffening-softening behaviour with increasing stress.

Development of a stress-mode sensitive viscoelastic constitutive relationship for asphalt concrete: experimental and numerical modeling

NASA Astrophysics Data System (ADS)

Karimi, Mohammad M.; Tabatabaee, Nader; Jahanbakhsh, H.; Jahangiri, Behnam

2017-08-01

Asphalt binder is responsible for the thermo-viscoelastic mechanical behavior of asphalt concrete. Upon application of pure compressive stress to an asphalt concrete specimen, the stress is transferred by mechanisms such as aggregate interlock and the adhesion/cohesion properties of asphalt mastic. In the pure tensile stress mode, aggregate interlock plays a limited role in stress transfer, and the mastic phase plays the dominant role through its adhesive/cohesive and viscoelastic properties. Under actual combined loading patterns, any coordinate direction may experience different stress modes; therefore, the mechanical behavior is not the same in the different directions and the asphalt specimen behaves as an anisotropic material. The present study developed an anisotropic nonlinear viscoelastic constitutive relationship that is sensitive to the tension/compression stress mode by extending Schapery's nonlinear viscoelastic model. The proposed constitutive relationship was implemented in Abaqus using a user material (UMAT) subroutine in an implicit scheme. Uniaxial compression and indirect tension (IDT) testing were used to characterize the viscoelastic properties of the bituminous materials and to calibrate and validate the proposed constitutive relationship. Compressive and tensile creep compliances were calculated using uniaxial compression, as well as IDT test results, for different creep-recovery loading patterns at intermediate temperature. The results showed that both tensile creep compliance and its rate were greater than those of compression. The calculated deflections based on these IDT test simulations were compared with experimental measurements and were deemed acceptable. This suggests that the proposed viscoelastic constitutive relationship correctly demonstrates the viscoelastic response and is more accurate for analysis of asphalt concrete in the laboratory or in situ.

Viscoelastic optical nonlocality of doped cadmium oxide epsilon-near-zero thin films

DOE Office of Scientific and Technical Information (OSTI.GOV)

Luk, Ting S.; De Ceglia, Domenico; Scalora, Michael

Optical nonlocalities are elusive and hardly observable in traditional plasmonic materials like noble and alkali metals. Here we experimentally observe and theoretically model viscoelastic nonlocalities in the infrared optical response of a doped, cadmium oxide epsilon-near-zero thin film. The nonlocality is clearly detectable thanks to the low damping rate of conduction electrons and the virtual absence of interband transitions at infrared wavelengths. We describe the motion of conduction electrons using a hydrodynamic model for a viscoelastic fluid, and find excellent agreement with experimental results. The electrons’ elasticity blue-shifts the infrared plasmonic resonance associated with the main epsilon-near-zero mode, and triggersmore » the onset of higher-order resonances due to the excitation of electron-pressure modes above the bulk plasma frequency. We also provide evidence of the existence of nonlocal damping, i.e., viscosity, in the motion of optically-excited conduction electrons using a combination of spectroscopic ellipsometry data and predictions based on the viscoelastic hydrodynamic model.« less

Viscoelastic performance of dielectric elastomer subject to different voltage stimulation

NASA Astrophysics Data System (ADS)

Sheng, Junjie; Zhang, Yuqing; Liu, Lei; Li, Bo; Chen, Hualing

2017-04-01

Dielectric elastomer (DE) is capable of giant deformation subject to an electric field, and demonstrates significant advantages in the potentially application of soft machines with muscle-like characteristics. Due to an inherent property of all macromolecular materials, DE exhibits strong viscoelastic properties. Viscoelasticity could cause a time-dependent deformation and lower the response speed and energy conversion efficiency of DE based actuators, thus strongly affect its electromechanical performance and applications. Combining with the rheological model of viscoelastic relaxation, the viscoelastic performance of a VHB membrane in a circular actuator configuration undergoing separately constant, ramp and sinusoidal voltages are analyzed both theoretically and experimentally. The theoretical results indicated that DE could attain a big deformation under a small constant voltage with a longer time or under a big voltage with a shorter time. The model also showed that a higher critical stretch could be achieved by applying ramping voltage with a lower rate and the stretch magnitude under sinusoidal voltage is much larger at a relatively low frequency. Finally, experiments were designed to validate the simulation and show well consistent with the simulation results.

Modeling electrical response of polymer-coated SAW resonators by equivalent circuit representation.

PubMed

Kshetrimayum, Roshan; Yadava, R D S; Tandon, R P

2011-07-01

The paper presents an equivalent circuit model of the polymer coated surface acoustic wave (SAW) resonators by combining coupling-of-mode (COM) description of SAW resonators and perturbation calculation of SAW propagation under polymer loading. An expression for the motional load produced by polymer coating is deduced in terms of COM parameters and polymer characteristics. In addition, expressions for the shifts in resonance frequency and attenuation due to polymer loading are obtained. Simulation results are presented for one-port and two-port resonator devices coated with viscoelastic thin polymer film. The influence of polymer film on resonator response is studied with regard to variations in film thickness and shear modulus. The model simplifies understanding of polymer-coated SAW sensors. Copyright © 2010 Elsevier B.V. All rights reserved.

2.5-D frequency-domain viscoelastic wave modelling using finite-element method

NASA Astrophysics Data System (ADS)

Zhao, Jian-guo; Huang, Xing-xing; Liu, Wei-fang; Zhao, Wei-jun; Song, Jian-yong; Xiong, Bin; Wang, Shang-xu

2017-10-01

2-D seismic modelling has notable dynamic information discrepancies with field data because of the implicit line-source assumption, whereas 3-D modelling suffers from a huge computational burden. The 2.5-D approach is able to overcome both of the aforementioned limitations. In general, the earth model is treated as an elastic material, but the real media is viscous. In this study, we develop an accurate and efficient frequency-domain finite-element method (FEM) for modelling 2.5-D viscoelastic wave propagation. To perform the 2.5-D approach, we assume that the 2-D viscoelastic media are based on the Kelvin-Voigt rheological model and a 3-D point source. The viscoelastic wave equation is temporally and spatially Fourier transformed into the frequency-wavenumber domain. Then, we systematically derive the weak form and its spatial discretization of 2.5-D viscoelastic wave equations in the frequency-wavenumber domain through the Galerkin weighted residual method for FEM. Fixing a frequency, the 2-D problem for each wavenumber is solved by FEM. Subsequently, a composite Simpson formula is adopted to estimate the inverse Fourier integration to obtain the 3-D wavefield. We implement the stiffness reduction method (SRM) to suppress artificial boundary reflections. The results show that this absorbing boundary condition is valid and efficient in the frequency-wavenumber domain. Finally, three numerical models, an unbounded homogeneous medium, a half-space layered medium and an undulating topography medium, are established. Numerical results validate the accuracy and stability of 2.5-D solutions and present the adaptability of finite-element method to complicated geographic conditions. The proposed 2.5-D modelling strategy has the potential to address modelling studies on wave propagation in real earth media in an accurate and efficient way.

Use of segmented constrained layer damping treatment for improved helicopter aeromechanical stability

NASA Astrophysics Data System (ADS)

Liu, Qiang; Chattopadhyay, Aditi; Gu, Haozhong; Liu, Qiang; Chattopadhyay, Aditi; Zhou, Xu

2000-08-01

The use of a special type of smart material, known as segmented constrained layer (SCL) damping, is investigated for improved rotor aeromechanical stability. The rotor blade load-carrying member is modeled using a composite box beam with arbitrary wall thickness. The SCLs are bonded to the upper and lower surfaces of the box beam to provide passive damping. A finite-element model based on a hybrid displacement theory is used to accurately capture the transverse shear effects in the composite primary structure and the viscoelastic and the piezoelectric layers within the SCL. Detailed numerical studies are presented to assess the influence of the number of actuators and their locations for improved aeromechanical stability. Ground and air resonance analysis models are implemented in the rotor blade built around the composite box beam with segmented SCLs. A classic ground resonance model and an air resonance model are used in the rotor-body coupled stability analysis. The Pitt dynamic inflow model is used in the air resonance analysis under hover condition. Results indicate that the surface bonded SCLs significantly increase rotor lead-lag regressive modal damping in the coupled rotor-body system.

Interrogating the viscoelastic properties of tissue using viscoelastic response (VISR) ultrasound

NASA Astrophysics Data System (ADS)

Selzo, Mallory Renee

Affecting approximately 1 in 3,500 newborn males, Duchenne muscular dystrophy (DMD) is one of the most common lethal genetic disorders in humans. Boys with DMD suffer progressive loss of muscle strength and function, leading to wheelchair dependence, cardiac and respiratory compromise, and death during young adulthood. There are currently no treatments that can halt or reverse the disease progression, and translating prospective treatments into clinical trials has been delayed by inadequate outcome measures. Current outcome measures, such as functional and muscle strength assessments, lack sensitivity to individual muscles, require subjective effort of the child, and are impacted by normal childhood growth and development. The goal of this research is to develop Viscoelastic Response (VisR) ultrasound which can be used to delineate compositional changes in muscle associated with DMD. In VisR, acoustic radiation force (ARF) is used to produce small, localized displacements within the muscle. Using conventional ultrasound to track the motion, the displacement response of the tissue can be evaluated against a mechanical model. In order to develop signal processing techniques and assess mechanical models, finite element method simulations are used to model the response of a viscoelastic material to ARF excitations. Results are then presented demonstrating VisR differentiation of viscoelastic changes with progressive dystrophic degeneration in a dog model of DMD. Finally, clinical feasibility of VisR imaging is demonstrated in two boys with DMD.

The Audible Human Project: Modeling Sound Transmission in the Lungs and Torso

NASA Astrophysics Data System (ADS)

Dai, Zoujun

Auscultation has been used qualitatively by physicians for hundreds of years to aid in the monitoring and diagnosis of pulmonary diseases. Alterations in the structure and function of the pulmonary system that occur in disease or injury often give rise to measurable changes in lung sound production and transmission. Numerous acoustic measurements have revealed the differences of breath sounds and transmitted sounds in the lung under normal and pathological conditions. Compared to the extensive cataloging of lung sound measurements, the mechanism of sound transmission in the pulmonary system and how it changes with alterations of lung structural and material properties has received less attention. A better understanding of sound transmission and how it is altered by injury and disease might improve interpretation of lung sound measurements, including new lung imaging modalities that are based on an array measurement of the acoustic field on the torso surface via contact sensors or are based on a 3-dimensional measurement of the acoustic field throughout the lungs and torso using magnetic resonance elastography. A long-term goal of the Audible Human Project (AHP ) is to develop a computational acoustic model that would accurately simulate generation, transmission and noninvasive measurement of sound and vibration within the pulmonary system and torso caused by both internal (e.g. respiratory function) and external (e.g. palpation) sources. The goals of this dissertation research, fitting within the scope of the AHP, are to develop specific improved theoretical understandings, computational algorithms and experimental methods aimed at transmission and measurement. The research objectives undertaken in this dissertation are as follows. (1) Improve theoretical modeling and experimental identification of viscoelasticity in soft biological tissues. (2) Develop a poroviscoelastic model for lung tissue vibroacoustics. (3) Improve lung airway acoustics modeling and its coupling to the lung parenchyma; and (4) Develop improved techniques in array acoustic measurement on the torso surface of sound transmitted through the pulmonary system and torso. Tissue Viscoelasticity. Two experimental identification approaches of shear viscoelasticity were used. The first approach is to directly estimate the frequency-dependent surface wave speed and then to optimize the coefficients in an assumed viscoelastic model type. The second approach is to measure the complex-valued frequency response function (FRF) between the excitation location and points at known radial distances. The FRF has embedded in it frequency-dependent information about both surface wave phase speed and attenuation that can be used to directly estimate the complex shear modulus. The coefficients in an assumed viscoelastic tissue model type can then be optimized. Poroviscoelasticity Model for Lung Vibro-acoustics. A poroviscoelastic model based on Biot theory of wave propagation in porous media was used for compression waves in the lungs. This model predicts a fast compression wave speed close to the one predicted by the effective medium theory at low frequencies and an additional slow compression wave due to the out of phase motion of the air and the lung parenchyma. Both compression wave speeds vary with frequency. The fast compression wave speed and attenuation were measured on an excised pig lung under two different transpulmonary pressures. Good agreement was achieved between the experimental observation and theoretical predictions. Sound Transmission in Airways and Coupling to Lung Parenchyma. A computer generated airway tree was simplified to 255 segments and integrated into the lung geometry from the Visible Human Male for numerical simulations. Acoustic impedance boundary conditions were applied at the ends of the terminal segments to represent the unmodeled downstream airway segments. Experiments were also carried out on a preserved pig lung and similar trends of lung surface velocity distribution were observed between the experiments and simulations. This approach provides a feasible way of simplifying the airway tree and greatly reduces the computation time. Acoustic Measurements of Sound Transmission in Human Subjects. Scanning laser Doppler vibrometry (SLDV) was used as a gold standard for transmitted sound measurements on a human subject. A low cost piezodisk sensor array was also constructed as an alternative to SLDV. The advantages and disadvantages of each technique are discussed.

Nonlinear viscoelastic characterization of polycarbonate

NASA Technical Reports Server (NTRS)

Caplan, E. S.; Brinson, H. F.

1982-01-01

Uniaxial tensile creep and recovery data from polycarbonate at six temperatures and six stress levels are analyzed for nonlinear viscoelastic constitutive modeling. A theory to account for combined effects of two or more accelerating factors is presented.

Prediction of water loss and viscoelastic deformation of apple tissue using a multiscale model.

PubMed

Aregawi, Wondwosen A; Abera, Metadel K; Fanta, Solomon W; Verboven, Pieter; Nicolai, Bart

2014-11-19

A two-dimensional multiscale water transport and mechanical model was developed to predict the water loss and deformation of apple tissue (Malus × domestica Borkh. cv. 'Jonagold') during dehydration. At the macroscopic level, a continuum approach was used to construct a coupled water transport and mechanical model. Water transport in the tissue was simulated using a phenomenological approach using Fick's second law of diffusion. Mechanical deformation due to shrinkage was based on a structural mechanics model consisting of two parts: Yeoh strain energy functions to account for non-linearity and Maxwell's rheological model of visco-elasticity. Apparent parameters of the macroscale model were computed from a microscale model. The latter accounted for water exchange between different microscopic structures of the tissue (intercellular space, the cell wall network and cytoplasm) using transport laws with the water potential as the driving force for water exchange between different compartments of tissue. The microscale deformation mechanics were computed using a model where the cells were represented as a closed thin walled structure. The predicted apparent water transport properties of apple cortex tissue from the microscale model showed good agreement with the experimentally measured values. Deviations between calculated and measured mechanical properties of apple tissue were observed at strains larger than 3%, and were attributed to differences in water transport behavior between the experimental compression tests and the simulated dehydration-deformation behavior. Tissue dehydration and deformation in the high relative humidity range ( > 97% RH) could, however, be accurately predicted by the multiscale model. The multiscale model helped to understand the dynamics of the dehydration process and the importance of the different microstructural compartments (intercellular space, cell wall, membrane and cytoplasm) for water transport and mechanical deformation.

Viscoelastic shock wave in ballistic gelatin behind soft body armor.

PubMed

Liu, Li; Fan, Yurun; Li, Wei

2014-06-01

Ballistic gelatins are widely used as a surrogate of biological tissue in blunt trauma tests. Non-penetration impact tests of handgun bullets on the 10wt% ballistic gelatin block behind soft armor were carried out in which a high-speed camera recorded the crater׳s movement and pressure sensors imbedded in the gelatin block recorded the pressure waves at different locations. The observed shock wave attenuation indicates the necessity of considering the gelatin׳s viscoelasticity. A three-element viscoelastic constitutive model was adopted, in which the relevant parameters were obtained via fitting the damping free oscillations at the beginning of the creep-mode of rheological measurement, and by examining the data of published split Hopkinson pressure bar (SHPB) experiments. The viscoelastic model is determined by a retardation time of 5.5×10(-5)s for high oscillation frequencies and a stress relaxation time of 2.0-4.5×10(-7)s for shock wave attenuation. Using the characteristic-line method and the spherical wave assumption, the propagation of impact pressure wave front and the subsequent unloading profile can be simulated using the experimental velocity boundary condition. The established viscoelastic model considerably improves the prediction of shock wave attenuation in the ballistic gelatin. Copyright © 2014 Elsevier Ltd. All rights reserved.

Measurement of viscoelastic properties of in vivo swine myocardium using Lamb Wave Dispersion Ultrasound Vibrometry (LDUV)

PubMed Central

Urban, Matthew W.; Pislaru, Cristina; Nenadic, Ivan Z.; Kinnick, Randall R.; Greenleaf, James F.

2012-01-01

Viscoelastic properties of the myocardium are important for normal cardiac function and may be altered by disease. Thus, quantification of these properties may aid with evaluation of the health of the heart. Lamb Wave Dispersion Ultrasound Vibrometry (LDUV) is a shear wave-based method that uses wave velocity dispersion to measure the underlying viscoelastic material properties of soft tissue with plate-like geometries. We tested this method in eight pigs in an open-chest preparation. A mechanical actuator was used to create harmonic, propagating mechanical waves in the myocardial wall. The motion was tracked using a high frame rate acquisition sequence, typically 2500 Hz. The velocities of wave propagation were measured over the 50–400 Hz frequency range in 50 Hz increments. Data were acquired over several cardiac cycles. Dispersion curves were fit with a viscoelastic, anti-symmetric Lamb wave model to obtain estimates of the shear elasticity, μ1, and viscosity, μ2 as defined by the Kelvin-Voigt rheological model. The sensitivity of the Lamb wave model was also studied using simulated data. We demonstrated that wave velocity measurements and Lamb wave theory allow one to estimate the variation of viscoelastic moduli of the myocardial walls in vivo throughout the course of the cardiac cycle. PMID:23060325

Analysis of Transient Shear Wave in Lossy Media.

PubMed

Parker, Kevin J; Ormachea, Juvenal; Will, Scott; Hah, Zaegyoo

2018-07-01

The propagation of shear waves from impulsive forces is an important topic in elastography. Observations of shear wave propagation can be obtained with numerous clinical imaging systems. Parameter estimations of the shear wave speed in tissues, and more generally the viscoelastic parameters of tissues, are based on some underlying models of shear wave propagation. The models typically include specific choices of the spatial and temporal shape of the impulsive force and the elastic or viscoelastic properties of the medium. In this work, we extend the analytical treatment of 2-D shear wave propagation in a biomaterial. The approach applies integral theorems relevant to the solution of the generalized Helmholtz equation, and does not depend on a specific rheological model of the tissue's viscoelastic properties. Estimators of attenuation and shear wave speed are derived from the analytical solutions, and these are applied to an elastic phantom, a viscoelastic phantom and in vivo liver using a clinical ultrasound scanner. In these samples, estimated shear wave group velocities ranged from 1.7 m/s in the liver to 2.5 m/s in the viscoelastic phantom, and these are lower-bounded by independent measurements of phase velocity. Copyright © 2018 World Federation for Ultrasound in Medicine and Biology. Published by Elsevier Inc. All rights reserved.

Measuring nanoscale viscoelastic parameters of cells directly from AFM force-displacement curves.

PubMed

Efremov, Yuri M; Wang, Wen-Horng; Hardy, Shana D; Geahlen, Robert L; Raman, Arvind

2017-05-08

Force-displacement (F-Z) curves are the most commonly used Atomic Force Microscopy (AFM) mode to measure the local, nanoscale elastic properties of soft materials like living cells. Yet a theoretical framework has been lacking that allows the post-processing of F-Z data to extract their viscoelastic constitutive parameters. Here, we propose a new method to extract nanoscale viscoelastic properties of soft samples like living cells and hydrogels directly from conventional AFM F-Z experiments, thereby creating a common platform for the analysis of cell elastic and viscoelastic properties with arbitrary linear constitutive relations. The method based on the elastic-viscoelastic correspondence principle was validated using finite element (FE) simulations and by comparison with the existed AFM techniques on living cells and hydrogels. The method also allows a discrimination of which viscoelastic relaxation model, for example, standard linear solid (SLS) or power-law rheology (PLR), best suits the experimental data. The method was used to extract the viscoelastic properties of benign and cancerous cell lines (NIH 3T3 fibroblasts, NMuMG epithelial, MDA-MB-231 and MCF-7 breast cancer cells). Finally, we studied the changes in viscoelastic properties related to tumorigenesis including TGF-β induced epithelial-to-mesenchymal transition on NMuMG cells and Syk expression induced phenotype changes in MDA-MB-231 cells.

TEMPERATURE-DEPENDENT VISCOELASTIC PROPERTIES OF THE HUMAN SUPRASPINATUS TENDON

PubMed Central

Huang, Chun-Yuh; Wang, Vincent M.; Flatow, Evan L.; Mow, Van C.

2009-01-01

Temperature effects on the viscoelastic properties of the human supraspinatus tendon were investigated using static stress-relaxation experiments and Quasi-Linear Viscoelastic (QLV) theory. Twelve supraspinatus tendons were randomly assigned to one of two test groups for tensile testing using the following sequence of temperatures: (1) 37°C, 27°C, and 17°C (Group I, n=6), or (2) 42°C, 32°C, and 22°C (Group II, n=6). QLV parameter C was found to increase at elevated temperatures, suggesting greater viscous mechanical behavior at higher temperatures. Elastic parameters A and B showed no significant difference among the six temperatures studied, implying that the viscoelastic stress response of the supraspinatus tendon is not sensitive to temperature over shorter testing durations. Using regression analysis, an exponential relationship between parameter C and test temperature was implemented into QLV theory to model temperature-dependent viscoelastic behavior. This modified approach facilitates the theoretical determination of the viscoelastic behavior of tendons at arbitrary temperatures. PMID:19159888

Accelerated viscoelastic characterization of T300-5208 graphite-epoxy laminates

NASA Technical Reports Server (NTRS)

Tuttle, M. E.; Brinson, H. F.

1985-01-01

A viscoelastic response scheme for the accelerated characterization of polymer-based composite laminates in applied to T300/5208 graphite/epoxy. The response of uni-directional specimens is modeled. The transient component of the viscoelastic creep compliance is assumed to follow a power law approximation. A recursive relationship is developed, based upon the Schapery single-integral equation, which allows approximation of a continuous time-varying uniaxial load using discrete steps in stress. The viscoelastic response of T300/5208 to transverse normal and shear stresses is determined unsing 90 deg and 10 deg off-axis tensile specimens. In each case the seven viscoelastic material parameters required in the analysis are determined experimentally using short-term creep and creep recovery tests. It is shown that an accurate measure of the power law exponent is crucial for accurate long-term prediction. A short term test cycle selection procedure is proposed, which should provide useful guidelines for the evaluation of other viscoelastic materials.

Helmholtz-Smoluchowski velocity for viscoelastic electroosmotic flows.

PubMed

Park, H M; Lee, W M

2008-01-15

Many biofluids such as blood and DNA solutions are viscoelastic and exhibit extraordinary flow behaviors, not existing in Newtonian fluids. Adopting appropriate constitutive equations these exotic flow behaviors can be modeled and predicted reasonably using various numerical methods. However, the governing equations for viscoelastic flows are not easily solvable, especially for electroosmotic flows where the streamwise velocity varies rapidly from zero at the wall to a nearly uniform velocity at the outside of the very thin electric double layer. In the present investigation, we have devised a simple method to find the volumetric flow rate of viscoelastic electroosmotic flows through microchannels. It is based on the concept of the Helmholtz-Smoluchowski velocity which is widely adopted in the electroosmotic flows of Newtonian fluids. It is shown that the Helmholtz-Smoluchowski velocity for viscoelastic fluids can be found by solving a simple cubic algebraic equation. The volumetric flow rate obtained using this Helmholtz-Smoluchowski velocity is found to be almost the same as that obtained by solving the governing partial differential equations for various viscoelastic fluids.

Computational modeling of electromechanical instabilities in dielectric elastomers (Conference Presentation)

NASA Astrophysics Data System (ADS)

Park, Harold

2016-04-01

Dielectric elastomers are a class of soft, active materials that have recently gained significant interest due to the fact that they can be electrostatically actuated into undergoing extremely large deformations. An ongoing challenge has been the development of robust and accurate computational models for elastomers, particularly those that can capture electromechanical instabilities that limit the performance of elastomers such as creasing, wrinkling, and snap-through. I discuss in this work a recently developed finite element model for elastomers that is dynamic, nonlinear, and fully electromechanically coupled. The model also significantly alleviates volumetric locking due that arises due to the incompressible nature of the elastomers, and incorporates viscoelasticity within a finite deformation framework. Numerical examples are shown that demonstrate the performance of the proposed method in capturing electromechanical instabilities (snap-through, creasing, cratering, wrinkling) that have been observed experimentally.

Maxwell: A semi-analytic 4D code for earthquake cycle modeling of transform fault systems

NASA Astrophysics Data System (ADS)

Sandwell, David; Smith-Konter, Bridget

2018-05-01

We have developed a semi-analytic approach (and computational code) for rapidly calculating 3D time-dependent deformation and stress caused by screw dislocations imbedded within an elastic layer overlying a Maxwell viscoelastic half-space. The maxwell model is developed in the Fourier domain to exploit the computational advantages of the convolution theorem, hence substantially reducing the computational burden associated with an arbitrarily complex distribution of force couples necessary for fault modeling. The new aspect of this development is the ability to model lateral variations in shear modulus. Ten benchmark examples are provided for testing and verification of the algorithms and code. One final example simulates interseismic deformation along the San Andreas Fault System where lateral variations in shear modulus are included to simulate lateral variations in lithospheric structure.

Modeling of viscoelastic properties of nonpermeable porous rocks saturated with highly viscous fluid at seismic frequencies at the core scale

NASA Astrophysics Data System (ADS)

Wang, Zizhen; Schmitt, Douglas R.; Wang, Ruihe

2017-08-01

A core scale modeling method for viscoelastic properties of rocks saturated with viscous fluid at low frequencies is developed based on the stress-strain method. The elastic moduli dispersion of viscous fluid is described by the Maxwell's spring-dash pot model. Based on this modeling method, we numerically test the effects of frequency, fluid viscosity, porosity, pore size, and pore aspect ratio on the storage moduli and the stress-strain phase lag of saturated rocks. And we also compared the modeling results to the Hashin-Shtrikman bounds and the coherent potential approximation (CPA). The dynamic moduli calculated from the modeling are lower than the predictions of CPA, and both of these fall between the Hashin-Shtrikman bounds. The modeling results indicate that the frequency and the fluid viscosity have similar effects on the dynamic moduli dispersion of fully saturated rocks. We observed the Debye peak in the phase lag variation with the change of frequency and viscosity. The pore structure parameters, such as porosity, pore size, and aspect ratio affect the rock frame stiffness and result in different viscoelastic behaviors of the saturated rocks. The stress-strain phase lags are larger with smaller stiffness contrasts between the rock frame and the pore fluid. The viscoelastic properties of saturated rocks are more sensitive to aspect ratio compared to other pore structure parameters. The results suggest that significant seismic dispersion (at about 50-200 Hz) might be expected for both compressional and shear waves passing through rocks saturated with highly viscous fluids.