Modelling framework for mass-growth II: The general case

The mass-growth modelling framework proposed in Soldatos (2013) is extended to consider an intermediate and, essentially, more general stage during which pseudo-elastic and pseudo-plastic mass-growth take place simultaneously. This “mixed” mass-growth stage is based on the postulation that the strain energy density for mass-growth depends not only on the deformation gradient and the state of pre-stress, but also on the rate of plastic deformation. The rate of total deformation may be thought of as composed by elastic and plastic parts, but the rule according which those parts are assembled is not specified in an a-priori manner. It is instead left to be determined through the solution of the particular mass-growth problem of interest. Still though, either of the purely pseudo-elastic and purely pseudo-plastic mass-growth stages considered in Soldatos (2013) may be obtained as particular cases of the present formulation. Uniform mass-growth of a pre-stressed solid cube (Soldatos, 2013) is again engaged in a pilot, example application where, however, the cube is now assumed to grow in the implied mixed, pseudo-elastic–plastic manner. All unknown quantities involved in this application are identified in a complete manner. It is thus verified that the proposed, mixed stage of growth may be further thought of as an overlap between consecutive pseudo-elastic and pseudo-plastic stages of mass-growth.     

Modelling framework for mass-growth

This study adopts a simple and rather known concept which suggests that stimulus of mass-growth may be built up in a elastic-like stage of mass-growth deformation, followed immediately afterwards by a stage of growth that may involve features of plastic flow. Such an analogy between known solid mechanics and mass-growth concepts suggests that the end of an elastic-like stage of mass growth can be marked through the satisfaction of a certain yield condition for mass-growth. This paper illustrates a manner in which hyper-elasticity and plasticity theory of solids may merge in order to produce such a suitable framework for modelling mass-growth features of matter, in general, or soft material tissue, in particular. In this context, it derives a new type of hyper-elasticity-type of constitutive equation which includes explicit contribution of mass-growth, and illustrates the efficiency of this new modelling framework through a relatively simple example application.     

On a variant of the Maxwell and Oldroyd-B models within the context of a thermodynamic basis

In this paper we develop models within a thermodynamic standpoint that are very similar in form to the classical Maxwell and Oldroyd-B models but differ from them in one important aspect, the manner in which they unload instantaneously from the deformed configuration. As long as the response is not instantaneous, the models that are derived cannot be differentiated from the Maxwell and Oldroyd-B models, respectively. The models can be viewed within the context of materials whose natural configuration evolves, the evolution being determined by the maximization of the rate of entropy production of the material. However, the underpinnings to develop the model are quite different from an earlier development by Rajagopal and Srinivasa [8] in that while the total response of the viscoelastic fluid satisfies the constraint of an incompressible material, the energy storage mechanism associated with the elastic response is allowed to be that for a compressible elastic solid and the dissipative mechanism associated with the viscous response allowed to be that for a compressible fluid, the total deformation however being isochoric. The analysis calls for a careful evaluation of firmly held customs in viscoelasticity wherein it is assumed that it is possible to subject a material to a purely instantaneous elastic response without any dissipation whatsoever. Finally, while the model developed by Rajagopal and Srinivasa [8] arises from the linearization of the non-linear elastic response that they chose and leads to a model wherein the instantaneous elastic response is isochoric, here we develop the model within the context of a different non-linear elastic response that need not be linearized but the instantaneous elastic response not necessarily being isochoric.     

A polygonal finite element formulation for modeling nearly incompressible materials

The objective of the present paper is to develop a finite element formulation for modeling nearly incompressible materials at large strains using polygonal elements. The present finite element formulation is a simplified version of the three-field mixed formulation and, in particular, it reduces the functional of the internal potential energy by expressing the field of the average volume-change in terms of the displacement field, where the latter is discretized using the Wachspress shape functions. The reduced mixed formulation eliminates the volumetric locking in nearly incompressible materials and enhances the computational efficiency as the static condensation is circumvented. A detailed implementation of the finite element formulation is presented in this study. Also, different example problems, including eigenvalue analysis, nonlinear patch test and other benchmark problems are presented for demonstrating the accuracy and the reliability of the developed formulation for polygonal elements.

     

On constitutive modelling of porous neo-Hookean composites

In this paper a hyperelastic constitutive model is developed for neo-Hookean composites with aligned continuous cylindrical pores in the finite elasticity regime. Although the matrix is incompressible, the composite itself is compressible because of the existence of voids. For this compressible transversely isotropic material, the deformation gradient can be decomposed multiplicatively into three parts: an isochoric uniaxial deformation along the preferred direction of the material (which is identical to the direction of the cylindrical pores here); an equi-biaxial deformation on the transverse plane (the plane perpendicular to the preferred direction); and subsequent shear deformation (which includes “along-fibre” shear and transverse shear). Compared to the multiplicative decomposition used in our previous model for incompressible fibre reinforced composites [Guo, Z., Peng, X.Q., Moran, B., 2006, A composites-based hyperelastic constitutive model for soft tissue with application to the human annulus fibrosus. J. Mech. Phys. Solids 54(9), 1952–1971], the equi-biaxial deformation is introduced to achieve the desired volume change. To estimate the strain energy function for this composite, a cylindrical composite element model is developed. Analytically exact strain distributions in the composite element model are derived for the isochoric uniaxial deformation along the preferred direction, the equi-biaxial deformation on the transverse plane, as well as the “along-fibre” shear deformation. The effective shear modulus from conventional composites theory based on the infinitesimal strain linear elasticity is extended to the present finite deformation regime to estimate the strain energy related to the transverse shear deformation, which leads to an explicit formula for the strain energy function of the composite under a general finite deformation state.     

Modelling of stress softening in elastomeric materials: foundations of simple theories

This work is concerned with formulation of constitutive relations for materials exhibiting the stress softening phenomenon (known as the Mullins effect) typical observed in elastomeric and other amorphous materials during loading–reloading cycles. It is assumed that microstructural changes in such materials during the deformation process can be represented by a single scalar-valued softening variable whose evolution is accompanied by microforces satisfying their own law of balance, besides the classical laws of mechanics underlying macroscopic deformation of a material. The constitutive equations are then derived in consistency with thermodynamics of irreversible processes with the restriction to purely mechanical theory. The general form of the derived constitutive equations is subsequently simplified through introduction of additional assumptions leading to various models of the stress softening phenomenon. As an illustration of the general theory, it is shown that the so-called pseudo-elastic model proposed in the literature may be derived without an ad hoc postulate of the variational principle.     

Development of an Artificial Compressibility Methodology with Implicit LU-SGS Method

A numerical method has been developed to solve the steady and unsteady incompressible Navier-Stokes equations in a two-dimensional, curvilinear coordinate system. The solution procedure is based on the method of artificial compressibility and uses a third-order flux-difference splitting upwind differencing scheme for convective terms and second-order center difference for viscous terms. A time-accurate scheme for unsteady incompressible flows is achieved by using an implicit real time discretization and a dual-time approach, which introduces pseudo-unsteady terms into both the mass conservation equation and momentum equations. An efficient fully implicit algorithm LU-SGS, which was originally derived for the compressible Eulur and Navier-Stokes equations by Jameson and Toon [1], is developed for the pseudo-compressibility formulation of the two dimensional incompressible Navier-Stokes equations for both steady and unsteady flows. A variety of computed results are presented to validate the present scheme. Numerical solutions for steady flow in a square lid-driven cavity and over a backward facing step and for unsteady flow in a square driven cavity with an oscillating lid and in a circular tube with a smooth expansion are respectively presented and compared with experimental data or other numerical results.     

Motion of a circular cylinder near a vertical wall

We consider the arbitrary motion of a circular cylinder in an ideal fluid near a vertical wall. This problem is usually solved in the approximate formulation with a degree of error which is difficult to assess, increasing with approach of the cylinder to the wall [1, 2], The exact solution has previously been carried out only for the case of purely circulatory flow about the cylinder [3].     

Integral transform solution of developing laminar duct flow in Navier-Stokes formulation

The generalized integral transform technique is employed in the hybrid numerical-analytical solution of the Navier-Stokes equations in streamfunction-only formulation, which govern the incompressible laminar flow of a Newtonian fluid within a parallel plate channel. Owing to the analytic nature of this approach, the outflow boundary condition for an infinite duct is handled exactly, and the error involved in considering finite duct lengths is investigated. The present error-controlled solutions are used to inspect the relative accuracy of previously reported purely numerical schemes and to compare Navier-Stokes and boundary layer formulations for various combinations of inlet conditions and Reynolds number.     

Computational aspects of a pseudo-elastic constitutive model for muscle properties in a soft-bodied arthropod

In this paper we discuss the computational implementation of a new constitutive model that describes the muscle properties in a soft-bodied arthropod. Qualitatively, the muscle tissues behave similar to particle-reinforced rubber and are capable of large non-linear elastic deformations, show a hysteretic behavior, and display stress softening during the first few cycles of repeated loading. Such behavior can be described by the framework of pseudo-elastic transversely isotropic hyperelasticity. The computational model assumes compressible overall response, and is based upon a multiplicative split of the deformation gradient tensor into volumetric and isochoric parts. Details regarding the implementation of the computational model in the context of an implicit finite element solution procedure are presented. In particular, an explicit expression is provided for the material tangent stiffness tensor. Results obtained utilizing the new implementation are also presented.     

近似不可压软材料动力分析的完全拉格朗日物质点法

物质点法(MPM)在模拟非线性动力问题时具有很好的效果,其已被广泛应用于许多大变形动力问题的分析中.然而传统的MPM在模拟不可压或近似不可压材料的动力学行为时会产生体积自锁,极大地影响模拟精度和收敛性.本文针对近似不可压软材料的大变形动力学行为,提出一种混合格式的显式完全拉格朗日物质点法(TLMPM).首先基于近似不可压软材料的体积部分应变能密度,引入关于静水压力的方程;之后将该方程与动量方程基于显式物质点法框架进行离散,并采用完全拉格朗日格式消除物质点跨网格产生的误差,提升大变形问题的模拟精度;对位移和压强场采用不同阶次的B样条插值函数并通过引入针对体积变形的重映射技术改进了算法,提升算法的准确性.此外,算法通过实施一种交错求解格式在每个时间步对位移场和压强场依次进行求解.最后,给出几个典型数值算例来验证本文所提出的混合格式TLMPM的有效性和准确性,计算结果表明该方法可以有效处理体积自锁,准确地模拟近似不可压软材料的大变形动力学行为.     

无网格Galerkin法与非协调元方法的耦合

无网格Galerkin法(EFGM)处理不可压缩问题时不存在自锁现象，有限元方法(FEM)也常被用来与其耦合以方便地施加边界条件和提高计算效率。在有限元方法中使用等参元，EFGM与FEM的耦合方法在处理不可压缩问题时仍然存在自锁现象。本文在有限元方法中，采用非协调元，将无网格kGalerkin法与非协调元耦合，保留了耦合方法的优点，且避免了求解不可压缩问题时的自锁现象。算例显示本文方法在分析平面应变不可压缩问题时能得到合理的结果。     

Distortion of Anisotropic Hyperelastic Solids Under Pure Pressure Loading: Compressibility, Incompressibility and Near-Incompressibility

The purpose of this note is to examine distortion during pure pressure loading for anisotropic hyperelastic solids. We contrast the corresponding issues in compressible and incompressible hyperelasticity, and then use these results to examine nearly incompressible materials. An anisotropic compressible hyperelastic solid will generally exhibit both volume change and distortion under hydrostatic pressure loading. In contrast, an incompressible hyperelastic solid—both isotropic and anisotropic—exhibits no change to its current state of deformation as the hydrostatic pressure is varied. Nearly incompressible hyperelastic materials are compressible, but approach an incompressible response in an appropriate limit. We examine this limiting process in the context of transverse isotropy. The issue arises as to how to implement a nearly incompressible version of a given truly incompressible material model. Here we examine how certain implementations eliminate distortion under pure pressure loading and why alternative implementations do not eliminate the distortion.     

A new finite element formulation for both compressible and nearly incompressible fluid dynamics

A new finite element formulation designed for both compressible and nearly incompressible viscous flows is presented. The formulation combines conservative and non‐conservative dependent variables, namely, the mass–velocity (density * velocity), internal energy and pressure. The central feature of the method is the derivation of a discretized equation for pressure, where pressure contributions arising from the mass, momentum and energy balances are taken implicitly in the time discretization. The method is applied to the analysis of laminar flows governed by the Navier–Stokes equations in both compressible and nearly incompressible regimes. Numerical examples, covering a wide range of Mach number, demonstrate the robustness and versatility of the new method. Copyright © 2000 John Wiley & Sons, Ltd.     

One-parameter families of isochoric deformations

The purpose of this paper is to show a construction for obtaining smooth one-parameter families of isochoric deformations. Our construction is useful in the variational approach to the equilibrium problem of an incompressible elastic solid.     

Theory of turbulent mixing at the interface of fluids in a gravity field

The theory of turbulent mixing at the interface of two media in accelerated motion was constructed in [1], and an approximate solution was given for incompressible fluids. The time variation of kinetic energy was neglected in the equation of balance for the kinetic energy of the turbulent motion. In [2] the characteristic turbulent velocity is averaged over the mixing region. This allows the initial equations to be solved allowing for the time variation of kinetic energy. It turns out that the resulting density profile roughly coincides with the profile of [1] within a wide range of variation of the initial density differential. In the present paper the equations for the mixing of incompressible fluids are studied in their complete form. It is established that the solutions of [1, 2] are applicable within a limited region, valid for small density ratios. The resulting solution is analyzed qualitatively, and it is shown that the density gradient at the mixing front is discontinuous. The dependence of the solution on two empirical constants is investigated. An approximate choice of the values of these constants is made on the basis of the theoretical considerations of [2, 3], and by comparison with the solution of [1]. The mixing asymmetry is found numerically as a function of the initial density differential. Quantitative characteristics of the solution are illustrated in graphs.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 74–81, July–August, 1976.     

Impact of a figure of revolution in a stream of incompressible fluid with separation of a jet

The impact interaction of bodies with a fluid in a flow with jet separation has been considered in [1–3], for example. This investigation was in the two-dimensional formulation. The present paper considers the three-dimensional problem of impact of a figure of revolution in a stream of an ideal incompressible fluid with separation of a jet in accordance with Kirchhoff's scheme. A boundary-value problem is formulated for the impact flow potential and solved by the Green's function method. A method for constructing the Green's function is described. Expressions are given for the coefficients of the apparent masses. The results are given of computer calculations of these coefficients in the case of a cone using the flow geometry of the corresponding two-dimensional problem.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 176–180, November–December, 1980.     

Parallel finite element computation of unsteady incompressible flows

A parallel semi-explicit iterative finite element computational procedure for modelling unsteady incompressible fluid flows is presented. During the procedure, element flux vectors are calculated in parallel and then assembled into global flux vectors. Equilibrium iterations which introduce some ‘local implicitness’ are performed at each time step. The number of equilibrium iterations is governed by an implicitness parameter. The present technique retains the advantages of purely explicit schemes, namely (i) the parallel speed-up is equal to the number of parallel processors if the small communication overhead associated with purely explicit schemes is ignored and (ii) the computation time as well as the core memory required is linearly proportional to the number of elements. The incompressibility condition is imposed by using the artificial compressibility technique. A pressure-averaging technique which allows the use of equal-order interpolations for both velocity and pressure, this simplifying the formulation, is employed. Using a standard Galerkin approximation, three benchmark steady and unsteady problems are solved to demonstrate the accuracy of the procedure. In all calculations the Reynolds number is less than 500. At these Reynolds numbers it was found that the physical dissipation is sufficient to stabilize the convective term with no need for additional upwind-type dissipation. © 1998 John Wiley & Sons, Ltd.