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.     

Differential quadrature method for analyzing nonlinear dynamic characteristics of viscoelastic plates with shear effects

The integro-partial differential equations governing the dynamic behavior of viscoelastic plates taking account of higher-order shear effects and finite deformations are presented. From the matrix formulas of differential quadrature, the special matrix product and the domain decoupled technique presented in this work, the nonlinear governing equations are converted into an explicit matrix form in the spatial domain. The dynamic behaviors of viscoelastic plates are numerically analyzed by introducing new variables in the time domain. The methods in nonlinear dynamics are synthetically applied to reveal plenty and complex dynamical phenomena of viscoelastic plates. The numerical convergence and comparison studies are carried out to validate the present solutions. At the same time, the influences of load and material parameters on dynamic behaviors are investigated. One can see that the system will enter into the chaotic state with a paroxysm form or quasi-periodic bifurcation with changing of parameters.     

Fractional Derivative Viscoelasticity at Large Deformations

A time domain viscoelastic model for large three-dimensional responses underisothermal conditions is presented. Internal variables with fractional orderevolution equations are used to model the time dependent part of the response. By using fractional order rate laws, the characteristics of the timedependency of many polymeric materials can be described using relatively fewparameters. Moreover, here we take into account that polymeric materials are often used in applications where the small deformations approximation does nothold (e.g., suspensions, vibration isolators and rubber bushings). A numerical algorithm for the constitutive response is developed and implemented into a finite element code forstructural dynamics. The algorithm calculates the fractional derivatives by means of the Grünwald–Lubich approach.Analytical and numerical calculations of the constitutive response in the nonlinearregime are presented and compared. The dynamicstructural response of a viscoelastic bar as well as the quasi-static response of athick walled tube are computed, including both geometrically and materiallynonlinear effects. Moreover, it isshown that by applying relatively small load magnitudes, the responses ofthe linear viscoelastic model are recovered.     

Direct numerical simulation of a bubble suspension in small amplitude oscillatory shear flow

Bubble suspensions can be found in many different fields and studying their rheology is crucial in order to improve manufacturing processes. When bubbles are added to a liquid, the magnitude of the viscosity changes and the behavior of the material is modified, giving it viscoelastic properties. For the purpose of this work, the suspended bubbles are considered to be monodisperse. It is assumed that Brownian motion and inertia can be neglected and that the fluid of the matrix is Newtonian and incompressible. The suspension is subject to an oscillatory strain while remaining in the linear regime. The resulting equations are solved in 3D with direct numerical simulation using a finite element discretization. Results of an ordered and random distribution of bubbles of volume fractions up to 40% are presented. The presence of bubbles has an opposite effect on the rheology of the suspension depending on the applied frequency. When the frequency is low, bubbles act as rigid fillers giving a rise to viscosity. On the contrary, when the frequency is high, the strain rate is being accommodated by the gaseous phase. Hence, bubbles deform, leading to a decrease of the viscosity.     

Suspension-based rheological modeling of crystallizing polymer melts

The applicability of suspension models to polymer crystallization is discussed. Although direct numerical simulations of flowing particle-filled melts are useful for gaining understanding about the rheological phenomena involved, they are computationally expensive. A more coarse-grained suspension model, which can relate the parameters in a constitutive equation for the two-phase material to morphological features, such as the volume fractions of differently shaped crystallites and the rheological properties of both phases, will be more practical in numerical polymer processing simulations. General issues, concerning the modeling of linear and nonlinear viscoelastic phenomena induced by rigid and deformable particles, are discussed. A phenomenological extension of linear viscoelastic suspension models into the nonlinear regime is proposed. A number of linear viscoelastic models for deformable particles are discussed, focusing on their possibilities in the context of polymer crystallization. The predictions of the most suitable model are compared to direct numerical simulation results and experimental data.     

非连续变形计算力学模型的粘弹性分析方法Ⅰ:基本理论

基于粘弹性广义有限单元和接触力元,发展了适用于多体相互作用系统非连续变形分析的粘弹性数值分析方法,通过虚功原理,给出了其分区参变量最小势能原理,从而阐明了其理论基础。粘弹性广义有限单元的本构关系可由粘弹性退化为弹性或刚性,因此本文所提出的方法可对由刚体、弹性体和粘弹性体所构成的复杂多体系统在外荷载作用下的力学行为进行数值模拟,同时能够比本文精确地直接得到多体之间的接触应力。     

Numerical simulation of particle migration in asymmetric bifurcation channel

In this work we provide numerical validation of the particle migration during flow of concentrated suspension in asymmetric T-junction bifurcation channel observed in a recent experiment [1]. The mathematical models developed to explain particle migration phenomenon basically fall into two categories, namely, suspension balance model and diffusive flux model. These models have been successfully applied to explain migration behavior in several two-dimensional flows. However, many processes often involve flow in complex 3D geometries. In this work we have carried out numerical simulation of concentrated suspension flow in 3D bifurcation geometry using the diffusive flux model. The simulation method was validated with available experimental and theoretical results for channel flow. After validation of the method we have applied the simulation technique to study the flow of concentrated suspensions through an asymmetric T-junction bifurcation composed of rectangular channels. It is observed that in the span-wise direction inhomogeneous concentration distribution that develops upstream persists throughout the inlet and downstream channels. Due to the migration of particles near the bifurcation section there is almost equal partitioning of flow in the two downstream branches. The detailed comparison of numerical simulation results is made with the experimental data.     

Review of the entry flow problem: Experimental and numerical

This paper is concerned with a review of both experimental and numerical studies of axisymmetric and planar entry flows which have been considered as test problems for the numerical simulation of viscoelastic fluids. The test of the method is usually based upon whether the numerical model predicts vortices in the entry corners. However, it is not clear as to whether one should observe vortices for all viscoelastic fluids. Polyacrylamide solutions and Boger fluids exhibit vortices in axisymmetric flow and the size of the vortex does increase with fluid elasticity. However, the vortex is nearly suppressed in planar entry flow. On the other hand, not all polymer melts are found to exhibit vortices in either axisymmetric or planar entry flow. It is our belief that the origin of vortices is not related to the elasticity based on shear flow propertes but to the behavior of the transient extensional viscosity. Certain polymer melts such as low density polyethylene exhibit vortices in both planar and axisymmetric flow along with unbounded stress growth at the start up of extensional flow. It is believed that the constitutive equations used in the numerical simulation must reflect this extensional behavior if vortices are to be predicted. A review of the numerical simulations concerned with entry flow shows that there is considerable doubt about the accuracy of the predictions for most of the studies. Even for those where the numerical solution is thought to be accurate, the magnitude of the stream function associated with the vortices is usually very low. None of the differential models used to date predicts strain hardening extensional viscosity, but those which are thought to predict vortices do rise more rapidly to the steady-state extensional viscosity values with time. It is recommended that the search of test fluids be widened beyond polymer solutions as there may already exist a number of polymer melts which behave similarly to the predictions of existing constitutive equations.     

The ultratough peeling of elastic tapes from viscoelastic substrates

The peeling of an elastic thin tape from a flat smooth viscoelastic substrate is investigated. Based on a Green function approach and on the translational invariance, a closed form analytical solution is proposed, which takes into account the viscoelastic dissipation in the substrate material.We find that peeling is prevented from taking place, only when the external force is smaller than the one predicted by Kendall's formula for elastic tapes on rigid substrates. However, we also find that, regardless of the value of the applied force, steady state detachment may occur when the elastic tape is sufficiently stiff. In this case, the constant peeling velocity can be modulated by properly defining the geometrical parameters and the material properties of tape and viscoelastic foundation. On the other hand, for relatively high peeling angles or compliant tapes a threshold value of the peeling force is found, above which the steady-state equilibrium is no longer possible and unstable detachment occurs.The present study contributes to shed light on the behavior of pressure sensitive adhesives in contact with viscoelastic substrates like the human skin. At the same time, it can be considered a first step towards a better understanding of the effect of viscoelastic dissipation on the fracture behavior of solids.     

Fracture-Flow-Enhanced Matrix Diffusion in Solute Transport Through Fractured Porous Media

Over the past few decades, significant progress of assessing chemical transport in fractured rocks has been made in laboratory and field investigations as well as in mathematic modeling. In most of these studies, however, matrix diffusion on fracture–matrix surfaces is considered as a process of molecular diffusion only. Mathematical modeling based on this traditional concept often had problems in explaining or predicting tracer transport in fractured rock. In this article, we propose a new conceptual model of fracture-flow-enhanced matrix diffusion, which correlates with fracture-flow velocity. The proposed model incorporates an additional matrix-diffusion process, induced by rapid fluid flow along fractures. According to the boundary-layer theory, fracture-flow-enhanced matrix diffusion may dominate mass-transfer processes at fracture–matrix interfaces, where rapid flow occurs through fractures. The new conceptual model can be easily integrated with analytical solutions, as demonstrated in this article, and numerical models, as we foresee. The new conceptual model is preliminarily validated using laboratory experimental results from a series of tracer breakthrough tests with different velocities in a simple fracture system. Validating of the new model with field experiments in complicated fracture systems and numerical modeling will be explored in future research.     

Extensional behavior influence on viscoelastic turbulent channel flow

In this work we use in the simulation of a viscoelastic turbulent channel flow a modification of the finitely extensible of non-linear elastic dumbbells with the Peterlin approximation (FENE-P) constitutive model for dilute polymer solutions, applicable to high extensional deformations. The new feature introduced by this modification is that the free energy of the polymer (since it is assumed to be entirely entropically driven) remains always bounded (FENE-PB). The characteristics of the model under steady shear flow, pure elongational flow and transient extensional behavior are presented. It is found that the FENE-PB model is more shear thinning than FENE-P. Most importantly, it also shows a higher extensional viscosity than the FENE-P model. Although the steady-state Trouton ratio asymptotically reaches at high extensional rates the same limit as the FENE-P model, the transition from the Newtonian value is sharper and faster. We use the FENE-PB model in direct numerical simulations (DNS) of viscoelastic turbulent channel flow using spectral approximations. The results for various statistics of the flow and the polymer conformation, when compared against those obtained with the original FENE-P model and the same rheological parameters, show an enhanced polymer-induced drag reduction effect and enhanced deformation of the polymer molecules. This indicates that it is not only the asymptotic but also details from the extensional rheological behavior that matter in quantitatively specifying turbulent viscoelastic flow behavior.     

Modeling the anisotropic finite-deformation viscoelastic behavior of soft fiber-reinforced composites

This paper presents constitutive models for the anisotropic, finite-deformation viscoelastic behavior of soft fiber-reinforced composites. An essential assumption of the models is that both the fiber reinforcements and matrix can exhibit distinct time-dependent behavior. As such, the constitutive formulation attributes a different viscous stretch measure and free energy density to the matrix and fiber phases. Separate flow rules are specified for the matrix and the individual fiber families. The flow rules for the fiber families then are combined to give an anisotropic flow rule for the fiber phase. This is in contrast to many current inelastic models for soft fiber-reinforced composites which specify evolution equations directly at the composite level. The approach presented here allows key model parameters of the composite to be related to the properties of the matrix and fiber constituents and to the fiber arrangement. An efficient algorithm is developed for the implementation of the constitutive models in a finite-element framework, and examples are presented examining the effects of the viscoelastic behavior of the matrix and fiber phases on the time-dependent response of the composite.     

Implications of Grain-Scale Mineralogical Heterogeneity for Radionuclide Transport in Fractured Media

The geological disposal of nuclear waste is based on the multi-barrier concept, comprising various engineered and natural barriers, to confine the radioactive waste and isolate it from the biosphere. Some of the planned repositories for high-level nuclear waste will be hosted in fractured crystalline rock formations. The potential of these formations to act as natural transport barriers is related to two coupled processes: diffusion into the rock matrix and sorption onto the mineral surfaces available in the rock matrix. Different in situ and laboratory experiments have pointed out the ubiquitous heterogeneous nature of the rock matrix: mineral surfaces and pore space are distributed in complex microstructures and their distribution is far from being homogeneous (as typically assumed by Darcy-scale coarse reactive transport models). In this work, we use a synthetically generated fracture–matrix system to assess the implications of grain-scale physical and mineralogical heterogeneity on cesium transport and retention. The resulting grain-scale reactive transport model is solved using high-performance computing technologies, and the results are compared with those derived from two alternative models, denoted as upscaled models, where mineral abundance is averaged over the matrix volume. In the grain-scale model, the penetration of cesium into the matrix is faster and the penetration front is uneven and finger-shaped. The analysis of the cesium breakthrough curves computed at two different points in the fracture shows that the upscaled models provide later first-arrival time estimates compared to the grain-scale model. The breakthrough curves computed with the three models converge at late times. These results suggest that spatially averaged upscaled parameters of sorption site distribution can be used to predict the late-time behavior of breakthrough curves but could be inadequate to simulate the early behavior.     

Nonlinear damping in a micromechanical oscillator

Nonlinear elastic effects play an important role in the dynamics of microelectromechanical systems (MEMS). A Duffing oscillator is widely used as an archetypical model of mechanical resonators with nonlinear elastic behavior. In contrast, nonlinear dissipation effects in micromechanical oscillators are often overlooked. In this work, we consider a doubly clamped micromechanical beam oscillator, which exhibits nonlinearity in both elastic and dissipative properties. The dynamics of the oscillator is measured in both frequency and time domains and compared to theoretical predictions based on a Duffing-like model with nonlinear dissipation. We especially focus on the behavior of the system near bifurcation points. The results show that nonlinear dissipation can have a significant impact on the dynamics of micromechanical systems. To account for the results, we have developed a continuous model of a geometrically nonlinear beam-string with a linear Voigt–Kelvin viscoelastic constitutive law, which shows a relation between linear and nonlinear damping. However, the experimental results suggest that this model alone cannot fully account for all the experimentally observed nonlinear dissipation, and that additional nonlinear dissipative processes exist in our devices.     

A time dependent micro-mechanical fiber composite model for inelastic zone growth in viscoelastic matrices

In this work, a fiber composite model is developed to predict the time dependent stress transfer behavior due to fiber fractures, as driven by the viscoelastic behavior of the polymer matrix, and the initiation and propagation of inelastic zones. We validate this model using in situ, room temperature, micro-Raman spectroscopy fiber strain measurements. Multifiber composites were placed under constant load creep tests and the fiber strains were evaluated with time after one fiber break occurred. These composite specimens ranged in fiber volume fraction and strain level. Comparison between prediction and MRS measurements allows us to characterize key in situ material parameters, the critical matrix shear strain for inelastic zones and interfacial frictional slip shear stress. We find that the inelastic zone is predominately either shear yielding or interfacial slipping, and the type depends on the local fiber spacing.     

A model for the viscoelastic behavior of polymers at finite strains

Summary A constitutive model is derived for the isothermal nonlinear viscoelastic response in polymers, which do not possess the separability property. The model is based on the concept of transient networks, and treats a polymer as a system of nonlinear elastic springs (adaptive links), which break and emerge due to micro-Brownian motion of chains. The breakage and reformation rates for adaptive links are assumed to depend on some strain energy density. The viscoelastic behavior is described by an integral constitutive equation, where the relaxation functions satisfy partial differential equations with coefficients depending on the strain history. Adjustable parameters of the model are found by fitting experimental data for a number of polymers in tension at strains up to 400 per cent. To validate the constitutive relations, we consider loading with different strain rates, determine adjustable parameters at one rate of strains, and compare prediction of the model with observations at another rate of strains. Fair agreement between experimental data and results of numerical simulation is demonstrated when the rates of strains differ by more than a decade. Received 1 July 1997; accepted for publication 7 October 1997     

Stress analysis of an elastic cracked layer bonded to a viscoelastic substrate

The effect of a viscoelastic substrate on an elastic cracked layer under an in-plane concentrated load is solved and discussed in this study. Based on a correspondence principle, the viscoelastic solution is directly obtained from the corresponding elastic one. The elastic solution in an anisotropic trimaterial is solved as a rapidly convergent series in terms of complex potentials via the successive iterations of the alternating technique in order to satisfy the continuity condition along the interfaces between dissimilar media. This trimaterial solution is then applied to a problem of a thin layer bonded to a half-plane substrate. Using the standard solid model to formulate the viscoelastic constitutive equation, the real-time stress intensity factors can be directly obtained by performing the numerical calculations. The results obtained in this paper are useful in studying the problem with bone defects where a crack is assumed to exist in an elastic body made of the cortical bone that is bonded to a viscoelastic substrate made of the cancellous bone.     

A general hereditary multimechanism-based deformation model with application to the viscoelastoplastic response of titanium alloys

The formulation of a general model for the hereditary behavior of materials, in the viscoelastic and viscoplastic regimes, is presented. In this, we utilize the complete-potential structure as a general framework, together with the notion of strain- and stress- partitioning in terms of separate contributions of several submechanisms (viscoelastic and viscoplastic) to the thermodynamic functions (stored energy and dissipation). Detailed numerical treatments are given for both (i) the implicit integration algorithm for the governing flow and evolutionary rate equations of the model, and (ii) the automated parameter-estimation methodology (using the software code COMPARE) for characterization. For illustration, a specific form of the model presented is characterized for the TIMETAL 21S material using a very comprehensive test matrix, including creep, relaxation, constant strain-rate tension tests, etc. Discussion of these correlations tests, together with comparisons to several other experimental results, are given to assess the performance and predictive capabilities of the present model as well as the effectiveness and practical utility of the algorithms proposed.     

Efficient Thermo-Chemistry Tabulation for Non-Premixed Combustion at High-Pressure Conditions

Propellant injection and turbulent combustion in high-pressure engines is often dominated by real-gas effects. However, previous studies suggested that the departure of the fluid properties from an ideal gas behavior has only a limited effect on the laminar flame structure. This is due to the fact that chemical reactions take place in the flame zone where the temperature is sufficiently high and molecular interactions are negligible, i.e., the ideal gas assumption is valid. On the other hand, various experimental and numerical studies of injection processes at high-pressure conditions demonstrated that real-gas effects can have a strong impact on the turbulent flow. Mixing is influenced by the rapid change of fluid properties. In this work, we exploit the gap in the fidelity of the thermodynamics model needed to describe the laminar flame structure and that needed to describe the turbulent flow field. We then propose a new real-gas flamelet model with increased numerical performance. The computational cost of the new formulation is not significantly higher than that of an ideal gas simulation. The performance of the method is analyzed and the error that is introduced by our assumptions is assessed by comparison to more complete modeling. Finally, the method is used to simulate a turbulent jet flame emanating from a coaxial injector at supercritical pressure and cryogenic oxidizer temperature. The results are compared with experimental OH^? images giving evidence of the suitability of the present method.