Mathematical modeling of electrorheological materials

Received May 11, 2000 / Published online January 23, 2001     

On constitutive equations for electrorheological materials

Constitutive equations for electrorheological (ER) fluids have been based on experimental results for steady shearing flows and constant electric fields. The fluids have been modeled as being rigid until a yield stress is reached. Additional stress is then proportional to the shear rate. Recent experimental results indicate that ER materials have a regime of solid-like response when deformed from a rest state. They behave in a viscoelastic-like manner under sinusoidal shearing and exhibit time-dependent response under sudden changes in shear rate or electric field. In this work, a constitutive theory for ER materials is presented which accounts for these recent experimental observations. The stress is given by a functional of the deformation gradient history and the electric field vector. Using the methods of continuum mechanics, a general three-dimensional constitutive equation is obtained. A sample constitutive equation is introduced which is then used to determine the response of an ER material for different shear histories. The calculated shear response is shown to be qualitatively similar to that observed experimentally.     

On the use of universal relations in the modeling of transversely isotropic materials

A material is of coaxial type   if the Cauchy stress tensor T$\mathit{T}$ and the strain tensor B$\mathit{B}$ are coaxial for all deformations. Clearly a hyperelastic material is of coaxial type if and only if it is isotropic. Here we present a weaker definition of materials of coaxial type. Anisotropic materials may be of a coaxial type in a weak sense if for a given specific  B$\mathit{B}$ we have that TB=BT$\mathit{TB}=\mathit{BT}$. We denote these materials B$\mathit{B}$-coaxial. We show that for transverse isotropic materials weak coaxial constitutive equations may be characterized using universal relations. We discuss the impact of B$\mathit{B}$-coaxial materials in the modeling of soft tissues. We conclude that B$\mathit{B}$-coaxial materials are a strong evidence that in real world materials two anisotropic invariants are always necessary to model in a meaningful and correct way single fiber reinforced materials.     

Nonlinear multiscale modeling of polymer materials

In this study, a hyperelastic multiscale modeling technique is used to predict elastic properties of polycarbonate and polyimide polymer systems using a set of widely accepted atomistic force fields. The model incorporates molecular simulations and a nonlinear, continuum mechanics-based, constitutive formulation that incorporates the behavior of the polymer materials as predicted from molecular simulations. The predicted properties of the polymers using multiple force fields are compared to experimentally measured values. Both static and dynamic molecular simulations are performed using molecular mechanics energy minimizations and molecular dynamics simulation techniques, respectively. The results of this study indicate that static molecular simulation is a useful tool to predict the bulk-level nonlinear mechanical behavior of polymers for finite deformations. It is found that the AMBER force field yields the most accurate predicted mechanical and physical properties of the modeled polymer systems compared to the other force fields used in this study.     

Homogenized modeling for vascularized poroelastic materials

A new mathematical model for the macroscopic behavior of a material composed of a poroelastic solid embedding a Newtonian fluid network phase (also referred to as vascularized poroelastic material), with fluid transport between them, is derived via asymptotic homogenization. The typical distance between the vessels/channels (microscale) is much smaller than the average size of a whole domain (macroscale). The homogeneous and isotropic Biot’s equation (in the quasi-static case and in absence of volume forces) for the poroelastic phase and the Stokes’ problem for the fluid network are coupled through a fluid-structure interaction problem which accounts for fluid transport between the two phases; the latter is driven by the pressure difference between the two compartments. The averaging process results in a new system of partial differential equations that formally reads as a double poroelastic, globally mass conserving, model, together with a new constitutive relationship for the whole material which encodes the role of both pore and fluid network pressures. The mathematical model describes the mutual interplay among fluid filling the pores, flow in the network, transport between compartments, and linear elastic deformation of the (potentially compressible) elastic matrix comprising the poroelastic phase. Assuming periodicity at the microscale level, the model is computationally feasible, as it holds on the macroscale only (where the microstructure is smoothed out), and encodes geometrical information on the microvessels in its coefficients, which are to be computed solving classical periodic cell problems. Recently developed double porosity models are recovered when deformations of the elastic matrix are neglected. The new model is relevant to a wide range of applications, such as fluid in porous, fractured rocks, blood transport in vascularized, deformable tumors, and interactions across different hierarchical levels of porosity in the bone.     

Three-dimensional micromechanical modeling of voided polymeric materials

A three-dimensional micromechanical unit cell model for particle-filled materials is presented. The cell model is based on a Voronoi tessellation of particles arranged on a body-centered cubic (BCC) array. The three-dimensionality of the present cell model enables the study of several deformation modes, including uniaxial, plane strain and simple shear deformations, as well as arbitrary principal stress states.The unit cell model is applied to studies on the micromechanical and macromechanical behavior of rubber-toughened polycarbonate. Different load cases are examined, including plane strain deformation, simple shear deformation and principal stress states. For a constant macroscopic strain rate, the different load cases show that the macroscopic flow strength of the blend decreases with an increase in void volume fraction, as expected. The main mechanism for plastic deformation is broad shear banding across inter-particle ligaments. The distributed nature of plastic straining acts to reduce the amount of macroscopic strain softening in the blend as the initial void volume fraction is increased. In the case of plane strain deformation, the plastic flow is observed to initiate across inter-particle ligaments in the direction of constraint. This particular mode of deformation could not have been captured using a two-dimensional, plane strain idealization of cylindrical voids in a matrix.The potential for localized crazing and/or cavitation in the matrix is addressed. It is observed that the introduction of voids acts to relieve hydrostatic stress in the matrix material, compared to the homopolymer. It is also seen that the predicted peak hydrostatic stress in the matrix is higher under plane strain deformation than under triaxial tension (with equal lateral stresses), for the same macroscopic stress triaxiality.The effect of void volume fraction on the macroscopic uniaxial tension behavior of the different blends is examined using a Considère construction for dilatant materials. The natural draw ratio was predicted to decrease with an increase in void volume fraction.     

Modulus-switching viscoelasticity of electrorheological networks

To form an electrorheological network (ERN), semiconducting nanoparticles were embedded in a polymer that was cross-linked to restrict particle motion. The microstructure ranged from random to aligned, depending on the degree of field-induced particle alignment during chemical network formation. We investigated in detail the softness effects of the matrix, having a relatively low storage modulus, on the dynamic rheological behavior of the ERN and analyzed its anisotropy. The anisotropy of the microstructure was probed rheologically with the modes of small-amplitude oscillatory shear (loading perpendicular to the field direction) and compression (loading in the field direction). The storage shear modulus was found to be a function of the applied electric field, particle volume fraction, and the pre-alignment electric field strength during the cross-linking reaction of the matrix, which governs the thickness of particle columns and intercolumn distance. Nonlinear behavior at small strain (below 0.1%) was conspicuous; this nonlinear viscoelasticity was accompanied by only a limited deformation of ordered connectivity. Throughout this study, we fabricated the ERN with the highly controllable modulus-switching effect acting in a shear-mode operation. Managing this anisotropy of an ERN by the electrical and chemical process is important in the design of smart materials that will provide improved stability and mechanical strength compared with fluid-type electrorheological materials and faster response time compared with that of conventional charged polymer gel.     

Computational modeling of surface phenomena in soft-wet materials

We consider the problem of extracting tribological information from experimental observations of contact with soft-wet materials. Particular attention is placed on simulating the response of two rotating cylinders of soft specimens placed in frictional contact, with a variable coefficient of friction dependent on the relative sliding velocity. The bulk behavior is modeled by means of a finite deformation viscoelasticity formulation, with constitutive parameters taken to be representative of hydrogels. We focus on the modeling of the surface behavior and employ a mortar-finite element contact formulation. Through a series of numerical studies, we demonstrate the strong sensitivity of the results to the choice of interfacial constitutive parameters. The difficulties of extracting such parameters using only experimental data and approximate analytical expressions are also examined.     

Physicomathematical modeling of the processes of capillary impregnation of porous materials

The process of capillary impregnation of porous materials is studied numerically. A physicomathematical model of liquid diffusion in a porous sample is proposed. The model involves an analytical presentation of the diffusion coefficient, which describes available experimental data. A method of solving one-dimensional unsteady problems of impregnation is developed and tested on a self-similar solution of the corresponding boundary-value problem of impregnation. If the impregnation process is sufficiently long, the motion of the liquid in the sample is described by a stable self-similar solution. A classification of moisture diffusion on the basis of the initial humidity on the sample boundary is proposed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 1, pp. 42–51, January–February, 2009.     

A phenomenological model of electrorheological fluids

The fundamental assumption of the paper is that the extra stress tensor of an electrorheological fluid is an isotropic tensor valued function of the rate of strain tensor D and the vector n (which characterizes the orientation and length N of the fibers formed by application of an electric field). The resulting constitutive equation for is supplemented by the solution of the previously studied time evolution equation for n. Plastic behavior for the shear and normal stresses is predicted. Anticipating that the action of increasing shear rate is i) to orient the fibers more and more in the direction of flow and ii) simultaneously to break up the fibers leads to the conclusion that for the same behavior is encountered as without an electric field. Using realistically possible approximation formulas for the dependence of and N on leads to the Bingham behavior for and power law behavior for large shear rates.

Multiscale modeling for ferroelectric materials: a transition from the atomic level to phase-field modeling

To link the atomic level and the mesoscale within a knowledge-based multiscale modeling approach for ferroelectric materials, a method is suggested to transfer results from first-principles calculations into a phase-field model. DFT calculations and atomistic simulations are applied and provide a set of intrinsic and extrinsic material properties for PbTiO₃ and tetragonal Pb(Zr_0.5Ti_0.5)O₃. The Helmholtz free energy of the phase-field model that contains all crystallographic and domain wall information is discussed in detail, and a sensitivity analysis is performed to identify the coefficients of the energy function. Then, a method is developed to adjust the coefficients of the Helmholtz free energy solely based on results from first-principles calculations. Full sets of adjusted energy coefficients for PbTiO₃ and Pb(Zr_0.5Ti_0.5)O₃ are presented and discussed, as well the limits of the suggested adjustment method.     

Constitutive modeling and the trousers test for fracture of rubber-like materials

The classical constitutive modeling of incompressible hyperelastic materials such as vulcanized rubber involves strain-energy densities that depend on the first two invariants of the strain tensor. The most well-known of these is the Mooney-Rivlin model and its specialization to the neo-Hookean form. While each of these models accurately predicts the mechanical behavior of rubber at moderate stretches, they fail to reflect the severe strain-stiffening and effects of limiting chain extensibility observed in experiments at large stretch. In recent years, several constitutive models that capture the effects of limiting chain extensibility have been proposed. Here we confine attention to two such phenomenological models. The first, proposed by Gent in 1996, depends only on the first invariant and involves just two material parameters. Its mathematical simplicity has facilitated the analytic solution of a wide variety of basic boundary-value problems. A modification of this model that reflects dependence on the second invariant has been proposed recently by Horgan and Saccomandi. Here we discuss the stress response of the Gent and HS models for some homogeneous deformations and apply the results to the fracture of rubber-like materials. Attention is focused on a particular fracture test, namely the trousers test where two legs of a cut specimen are pulled horizontally apart. It is shown that the cut position plays a key role in the fracture analysis, and that the effect of the cut position depends crucially on the constitutive model employed. For stiff rubber-like or biological materials, it is shown that the influence of the cut position is diminished. In fact, for linearly elastic materials, the critical driving force for fracture is independent of the cut position. It is also shown that the limiting chain extensibility models predict finite fracture toughness as the cut position approaches the edge of the specimen whereas classical hyperelastic models predict unbounded toughness in this limit. The results are relevant to the structural integrity of rubber components such as vibration isolators, vehicle tires, earthquake bearings, seals and flexible joints.