Characterizing the creep of viscoelastic materials by fractal derivative models

In this paper, we make the first attempt to apply the fractal derivative to modeling viscoelastic behavior. The methodology of scaling transformation is utilized to obtain the creep modulus and relaxation compliance for the proposed fractal Maxwell and Kelvin models. Comparing with the fractional derivatives reported in the literature, the fractal derivative as a local operator has lower calculation costs and memory storage requirements. Moreover, numerical results show that the proposed fractal models require fewer parameters, have simpler mathematical expression and result in higher accuracy than the classical integer-order derivative models. Results further confirm that the proposed fractal models can characterize the creep behavior of viscoelastic materials.     

Viscoelasticity of brain corpus callosum in biaxial tension

Computational models of the brain rely on accurate constitutive relationships to model the viscoelastic behavior of brain tissue. Current viscoelastic models have been derived from experiments conducted in a single direction at a time and therefore lack information on the effects of multiaxial loading. It is also unclear if the time-dependent behavior of brain tissue is dependent on either strain magnitude or the direction of loading when subjected to tensile stresses. Therefore, biaxial stress relaxation and cyclic experiments were conducted on corpus callosum tissue isolated from fresh ovine brains. Results demonstrated the relaxation behavior to be independent of strain magnitude, and a quasi-linear viscoelastic (QLV) model was able to accurately fit the experimental data. Also, an isotropic reduced relaxation tensor was sufficient to model the stress-relaxation in both the axonal and transverse directions. The QLV model was fitted to the averaged stress relaxation tests at five strain magnitudes while using the measured strain history from the experiments. The resulting model was able to accurately predict the stresses from cyclic tests at two strain magnitudes. In addition to deriving a constitutive model from the averaged experimental data, each specimen was fitted separately and the resulting distributions of the model parameters were reported and used in a probabilistic analysis to determine the probability distribution of model predictions and the sensitivity of the model to the variance of the parameters. These results can be used to improve the viscoelastic constitutive models used in computational studies of the brain.     

Modeling long-term creep rupture by debonding in unidirectional fibre-reinforced composites

Delayed fracture due to debonding can be observed in many unidirectional fibre-reinforced composites when the fibre/matrix interface experiences creep. The aim of this work is to describe such a phenomenon within the recently proposed modeling framework of transverse isotropy that allows for a neat decomposition of the mechanical behavior into fibre-directional, transverse, and pure shear parts. Specifically, debonding is here chosen to be governed by the tension transverse to the fibres. One can then speak of a mode-I debonding if use is made of the terminology adopted in fracture mechanics. On another hand, the time-dependent response is attributed to the matrix constituent. As the role of this latter is to deform and support stresses primarily in shear, a viscoelastic behavior is introduced that affects solely the pure shear part of the behavior. We show that both characteristics can be easily embedded into the aforementioned formulation. Among others, the occurrence of tertiary creep is made possible to predict. It is otherwise found that the predicted debonding path always propagates along the direction of the fibres in agreement with many experimental observations found in the literature. On the numerical side, the algorithmic treatment of debonding is independent of the one for viscoelasticity. This renders the implementation within the context of the finite element method very easy.     

Viscoelastic analysis of single-component and composite PEG and alginate hydrogels简

Hydrogels, three-dimensional hydrophilic polymer networks, are appealing candidate materials for study- ing the cellular microenvironment as their substantial water content helps to better mimic soft tissue. However, hydrogels can lack mechanical stiffness, strength, and tough- ness. Composite hydrogel systems have been shown to improve upon mechanical properties compared to their single- component counterparts. Poly （ethylene glycol） dimethacrylate （PEGDMA） and alginate are polymers that have been used to form hydrogels for biological applications. Single- component and composite PEGDMA and alginate systems were fabricated with a range of total polymer concentrations. Bulk gels were mechanically characterized using spherical indentation testing and a viscoelastic analysis framework. An increase in shear modulus with increasing polymer con- centration was demonstrated for all systems. Alginate hydro- gels were shown to have a smaller viscoelastic ratio than the PEGDMA gels, indicating more extensive relaxation over time. Composite alginate and PEGDMA hydrogels exhib- ited a combination of the mechanical properties of the con- stituents, as well as a qualitative increase in toughness. Additionally, multiple hydrogel systems were produced that had similar shear moduli, but different viscoelastic behaviors. Accurate measurement of the mechanical properties of hydrogels is necessary in order to determine what parameters are key in modeling the cellular microenvironment.     

Peeling behavior of a viscoelastic thin-film on a rigid substrate

In order to study the adhesion mechanism of a viscoelastic thin-film on a substrate, peeling experiment of a viscoelastic polyvinylchloride (PVC) thin-film on a rigid substrate (glass) is carried out. The effects of peeling rate, peeling angle, film thickness, surface roughness and the interfacial adhesive on the peel-off force are considered. It is found that both the viscoelastic properties of the film and the interfacial adhesive contribute to the rate-dependent peel-off force. For a fixed peeling rate, the peel-off force decreases with the increasing peeling angle. Increasing film thickness or substrate roughness leads to an increase of the peel-off force. Viscoelastic energy release rate in the present experiment can be further predicted by adopting a recently published theoretical model. It is shown that the energy release rate increases with the increase of peeling rates or peeling angles. The results in the present paper should be helpful for understanding the adhesion mechanism of a viscoelastic thin-film.     

Reliability evaluation and control for wideband noise-excited viscoelastic systems

Reliability of first-passage type for wideband noise-excited viscoelastic systems and the quasi-optimal bounded control strategy for maximizing system reliability are investigated. The viscoelastic term is approximately replaced by equivalent damping and stiffness separately. By using the stochastic averaging method based on the generalized harmonic functions, the averaged Itô stochastic differential equation is obtained for the system amplitude. The associated backward Kolmogorov equation is derived and solved to obtain the system reliability. By applying the dynamic programming principle to the averaged system, the quasi-optimal bounded control is devised by maximizing system reliability. The application of the proposed analytical procedures and the effectiveness of the control strategy are illustrated through one example.     

Development of secondary flows in viscoelastic curved ducts and the influence of outlet region

This paper presents numerical simulations of Newtonian and viscoelastic flows through a 180° curved duct of square cross section with a long straight outlet region. A particular attention is paid to the development of the flow in the output rectangular region after the curved part. The viscoelastic fluid is modeled using the constitutive equation proposed by Phan–Thien–Tanner (PTT). The numerical results, obtained with a finite-volume method, are shown for three different Dean numbers (125,137,150)$(125,137,150)$ and for three Deborah numbers (0.1,0.2,0.3)$(0.1,0.2,0.3)$. The necessary outlet length to impose boundary conditions is presented and discussed for these cases. Streamlines and vortex formation are shown to illustrate and analyze the evolution of the secondary flow in this region.     

具有不规则自由表面的粘弹性介质中声波模拟的有限差分算法

声波在具有不规则自由表面的粘弹性介质传播的数值模拟工作在地震勘探,地震预测中非常重要,特别是当模型有起伏的自由界面和较强的衰减特性时更是如此。基于镜像方法和直接方法,本文发展了一种二维有限差分算法,可以模拟不规则自由表面引起的声波散射问题。该方法将自由表面条件与速度应力方程结合求解粘弹性波动方程,在垂直和水平自由段及拐点处设置相应的边界条件。该方法假设自由表面穿过剪切应力和相应参数的网格点。为了提高计算精度,分别计算了水平和垂直方向上应力镜像值,而对质点速度,采用先水平方向后垂直方向分析进行镜像计算和更新。将粘弹性水平自由表面镜像方法和不同倾斜度平滑自由表面的改进方法的数值结果进行对比,并验证了算法的精度。     

Deterministic and stochastic equations of motion arising in Oldroyd fluids of order one: existence,uniqueness, exponential stability and invariant measures

Abstract

In this work, we consider the two-dimensional viscoelastic fluid flow equations, arising from the Oldroyd model for the non-Newtonian fluid flows. We investigate the well-posedness of such models in two-dimensional bounded and unbounded (Poincaré domains) domains, both in deterministic and stochastic settings. The existence and uniqueness of weak solution in the deterministic case is proved via a local monotonicity property of the linear and nonlinear operators and a localized version of the Minty-Browder technique. Some results on the exponential stability of stationary solutions are also established. The global solvability results for the stochastic counterpart are obtained by a stochastic generalization of the Minty-Browder technique. The exponential stability results in the mean square as well as in the pathwise (almost sure) sense are also discussed. Using the exponential stability results, we finally prove the existence of a unique invariant measure, which is ergodic and strongly mixing.     

Dynamic toughness in elastic nonlinear viscous solids

This work addresses the interrelationship among dissipative mechanisms—material separation in the fracture process zone (FPZ), nonelastic deformation in the surrounding background material and kinetic energy—and how they affect the macroscopic dynamic fracture toughness as well as the limiting crack speed in strain rate sensitive materials. To this end, a micromechanics-based model for void growth in a nonlinear viscous solid is incorporated into a microporous strip of cell elements that forms the FPZ. The latter is surrounded by background material described by conventional constitutive relations. In the first part of the paper, the background material is assumed to be purely elastic. Here, the computed dynamic fracture toughness is a convex function of crack velocity. In the second part, the background material as well as the FPZ are described by similar rate-sensitivity parameters. Voids grow in the strain rate strengthened FPZ as the crack velocity increases. Consequently, the work of separation increases. By contrast, the inelastic dissipation in the background material appears to be a concave function of crack velocity. At the lower crack velocity regime, where dissipation in the background material is dominant, toughness decreases as crack velocity increases. At high crack velocities, inelastic deformation enhanced by the inertial force can cause a sharp increase in toughness. Here, the computed toughness increases rapidly with crack velocity. There exist regimes where the toughness is a non-monotonic function of the crack velocity. Two length scales—the width of the FPZ and the ratio of the shear wave speed to the reference strain rate—can be shown to strongly affect the dynamic fracture toughness. Our computational simulations can predict experimental data for fracture toughness vs. crack velocity reported in several studies for amorphous polymeric materials.