Flow of elasto-viscoplastic liquids through an axisymmetric expansion–contraction

In this work, we used a recently proposed constitutive model that addresses the elastic behavior of viscoplastic liquids. We used the equation in an expansion–contraction axisymmetric flow, and compared the results with experimental data found in the literature. Steady, inertialess numerical solutions were obtained by solving the balance equations of mass and momentum via the finite element method. We obtained the velocity and stress fields for different combinations of the governing parameters, and observed that elasticity leads to significant changes in the shape and position of the yield surface, affecting both the extra-stress and the rate of deformation fields in the cavity. The trends observed are in qualitative agreement with visualization results available in the literature.     

Characterization of the post-necking strain hardening behavior using the virtual fields method

The present study aims at characterizing the post-necking strain hardening behavior of three sheet metals having different hardening behavior. Standard tensile tests were performed on sheet metal specimens up to fracture and heterogeneous logarithmic strain fields were obtained from a digital image correlation technique. Then, an appropriate elasto-plastic constitutive model was chosen. Von Mises yield criterion under plane stress and isotropic hardening law were considered to retrieve the relationship between stress and strain. The virtual fields method (VFM) was adopted as an inverse method to determine the constitutive parameters by calculating the stress fields from the heterogeneous strain fields. The results show that the choice of a hardening law which can describe the hardening behavior accurately is important to derive the true stress–strain curve. Finally, post-necking hardening behavior was successfully characterized up to the initial stage of localized necking using the VFM with Swift and modified Voce laws.     

Flow of elasto-viscoplastic liquids through a planar expansion–contraction

The steady flow of incompressible elasto-viscoplastic liquids through a planar expansion–contraction is investigated. A novel constitutive model is employed to describe the mechanical behavior of the flowing liquids. Numerical solutions of the constitutive and conservation equations were obtained via a finite element method to investigate the role of elasticity, yield stress, and inertia. The fields of velocity, stress, elastic strain, and rate of strain were obtained for different combinations of the governing parameters. It was observed that these fields, as well as the shape and position of the yield surface, are all strong functions of elasticity, yield stress, and inertia. The trends observed agree well with previous numerical and visualization results available in the literature. The present work offers a detailed study on the effects of elasticity, presenting, in particular, the fields of elastic strain.     

Multiaxial yield surface of transversely isotropic foams: Part I—Modeling

A new yield criterion is proposed for transversely isotropic solid foams. Its derivation is based on the hypothesis that the yielding in foams is driven by the total strain energy density, rather than a completely phenomenological approach. This allows defining the yield surface with minimal number of parameters and does not require complex experiments. The general framework used leads to the introduction of new scalar measures of stress and strain (characteristic stress and strain) for transversely isotropic foams. Furthermore, the central hypothesis that the total strain energy density drives yielding in foams ascribes to the characteristic stress an analogous role of von Mises stress in metal plasticity. Unlike the overwhelming majority of yield models in literature the proposed model recognizes the tension–compression difference in yield behavior of foams through a linear mean stress term. Predictions of the proposed yield model are in excellent agreement with the results of uniaxial, biaxial and triaxial FE analyses implemented on both isotropic and transversely isotropic Kelvin foam models.     

STRUCTURAL DAMAGE MODEL OF UNSATURATED EXPANSIVE SOIL AND ITS APPLICATION IN MULTI-FIELD COUPLE ANALYSIS ON EXPANSIVE SOIL SLOPE

Focused on the sensitivity to climate change and the special mechanical characteristics of undisturbed expansive soil, an elastc-plastic damage constitutive model was proposed based on the mechanics of unsaturated soil and the mechanics of damage. Undisturbed expansive soil was considered as a compound of non-damaged part and damaged part. The behavior of the non-damaged part was described using non-linear constitutive model of unsaturated soil. The property of the damaged part was described using a damage evolution equation and two yield surfaces, i.e., loading yield (LY) and shear yield (SY). Furthermore, a consolidation model for unsaturated undisturbed expansive soil was established and a FEM program named UESEPDC was designed. Numerical analysis on solid-liquid-gas tri-phases and multi-field couple problem was conducted for four stages and fields of stress, displacement, pore water pressure, pore air pressure, water content, suction, and the damage region as well as plastic region in an expansive soil slope were obtained.     

Measuring static yield stress of electrorheological fluids using the slotted plate device

An electrorheological (ER) response is defined as the dramatic change in rheological properties of a suspension of small particles due to the application of a large electric field transverse to the direction of flow. ER fluids are typically composed of nonconducting or semiconducting particles dispersed in a nonconducting continuous phase. A sufficiently large electric field will cause ER fluids to solidify, giving rising to a yield stress. Many applications in torque and stress transfer devices were proposed employing the reversible yielding behavior of ER fluids. Successful applications depend on a large yield stress of ER fluids and therefore accurate measurements of the yield stress of ER fluids are required. Reported experimental yield stresses of ER fluids have been dynamic yield stresses obtained by extrapolating the shear stress–shear rate data to zero-shear rate. It would be very helpful to the understanding of ER behaviors and the applications of ER fluids to be able to measure the static yield stress of ER fluids accurately. The slotted plate technique has been shown to be a successful method to determine the static yield stress of suspensions. The values obtained via the slotted plate method are static yield stress as the platform is designed for extremely low-speed motion. In this study, we modified the slotted plate device for the application of large electric fields and measured the static yield stress of TiO₂ ER fluids under various electric fields. The measured static yield stress values are also compared with the static yield stress values from a commercial rheometer.     

A plasticity model for cellular materials with open-celled structure

Based on a rigid-plastic material model that obeys the von Mises yield criterion, the plastic behavior of foams with an open-celled structure is studied in this paper using a single unit cell. An approximate continuum plasticity model is developed within the framework of the upper bound theorem of plasticity to describe the yield behavior of foams. The microscopic velocity fields are derived for the unit cell, which satisfy the incompressibility and the kinematic boundary conditions, and expressed in macroscopic rate of deformation. From the microscopic velocity fields, a macroscopic yield function is developed for foams under multi-axial stresses and includes the effects of the hydrostatic stress due to the void presence and growth. The dependency of the derived yield surfaces of foams on their relative densities is studied. The plastic behavior of foams is also studied numerically using the finite element method. The newly developed plasticity model is compared with the finite element analysis results and other available foam models and then correlated with the finite element results.     

Modeling the effects of material non-linearity using moving window micromechanics

In the analysis of materials with random heterogeneous microstructure the assumption is often made that material behavior can be represented by homogenized or effective properties. While this assumption yields accurate results for the bulk behavior of composite materials, it ignores the effects of the random microstructure. The spatial variations in these microstructures can focus, initiate and propagate localized non-linear behavior, subsequent damage and failure. In previous work a computational method, moving window micromechanics (MW), was used to capture microstructural detail and characterize the variability of the local and global elastic response. Digital images of material microstructure described the microstructure and a local micromechanical analysis was used to generate spatially varying material property fields. The strengths of this approach are that the material property fields can be consistently developed from digital images of real microstructures, they are easy to import into finite element models (FE) using regular grids, and their statistical characterizations can provide the basis for simulations further characterizing stochastic response. In this work, the moving window micromechanics technique was used to generate material property fields characterizing the non-linear behavior of random materials under plastic yielding; specifically yield stress and hardening slope, post yield. The complete set of material property fields were input into FE models of uniaxial loading. Global stress strain curves from the FE–MW model were compared to a more traditional micromechanics model, the generalized method of cells. Local plastic strain and local stress fields were produced which correlate well to the microstructure. The FE–MW method qualitatively captures the inelastic behavior, based on a non-linear flow rule, of the sample continuous fiber composites in transverse uniaxial loading.     

A non-associated flow rule for sheet metal forming

An improved model of material behavior is proposed that shows good agreement with experimental data for both yield and plastic strain ratios in uniaxial, equi-biaxial, and plane-strain tension under proportional loading for steel, aluminum and possibly other alloys. This model is based on a non-associated flow rule in which the plastic potential and yield surface functions are defined by quadratic functions of the stress tensor. The plastic potential aspect of the model is identical to that proposed by Hill for a quadratic anisotropic plastic potential defined in terms of measured r values. The new model differs in that the yield surface, although also defined by a quadratic function of the stress tensor, is defined independently of the plastic potential in terms of measured yield stresses. The model is developed and implemented in an FEM code that is based on a convected coordinate system. Since the associated flow rule, which assumes equivalency between the plastic potential and yield functions, is commonly accepted as a valid law in the theory of plastic deformation of most metals, the arguments for the associated flow rule are also discussed.     

Rheological studies and optimization of Herschel-Bulkley flow parameters of viscous karaya polymer suspensions using GA and PSO algorithms

Rheological characteristics of gum karaya suspensions which is proposed as a fracturing fluid were investigated with the main objective to determine the yield stress and other rheological parameters using various models. The flow hysteresis confirms the thixotropic behavior of fluid with increased structural breakdown at higher concentration and temperature. An empirical model developed for the studied samples accurately predicts the temperature and polymer concentration sensitivity of the apparent viscosity. The Herschel-Bulkley model showed the best fit to the experimental data; however, the yield stress obtained from some of the samples using nonlinear regression (NL) method reported physically insignificant, negative values. The proposed optimization technique, i.e., “Particle Swarm Optimization” offered the most realistic results with faster convergence over genetic algorithm making it a better choice. The oscillatory study provided more reliable yield stress values and revealed the thermogelation behavior of polymer suspensions making it suitable for high-temperature fracturing application.     

岩土的一种强度准则及其变换应力法

基于岩土摩擦性,假设岩土破坏是由其物理空间内特征面上的应力比决定,提出了等效应力比的概念,即物理空间特征面上的剪应力合力与正应力合力的比值.在二维条件下,等效应力比可表示为σ-τ坐标系下与摩尔圆相切的直线扣除截距正切值;在三维条件下,假设在XYZ空间内存在一三维物理空间平面,此三维空间特征平面的等效应力比为影响材料强度特性的决定性因素,基于上述三维空间特征面建立了强度准则并称之为a准则.SMP准则以及广义Mises准则都是a准则的特例,当二维坐标中的截距为零时,则强度准则退化为SMP(spatially mobilized plane)强度准则,而当正切角为零时,则强度准则退化为广义Mises准则.而当截距与外切角均不为零时,则强度准则为介于上述两者之间的一种强度准则,在偏平面上为介于SMP曲边三角形与广义Mises圆形之间的曲边三角形.在子午面上,采用考虑岩土压剪耦合的屈服准则,破坏准则采用幂函数表达式.在偏平面上提出了基于a准则的形状函数,并采用真三维应力状态表示的破坏强度准则表示在三轴压缩路径下以p,q二维应力变量表达的准则公式,推导得到了基于a准则的变换应力公式,可简单地将一般以p,q为基本变量的二维模型转变为三维应力模型.通过强度以及多种应力路径的测试对比,验证了a准则及基于该准则的变换应力公式的合理性.     

A cyclic anisotropic-plasticity model for metal matrix composites

Based on a six parameter general anisotropic yield surface proposed earlier by Voyiadjis and Thiagarajan (An Anisotropic Yield Surface Model for Directionally Reinforced Metal Matrix Composites, Int. J. Plasticity [1995]), a cyclic plasticity model to model the behavior of directionally reinforced metal matrix composite, has been proposed here. Apart from being able to model different initial yielding behavior along different stress directions, a number of features have been incorporated into the plasticity model. They include the usage of a proposed non-associative flow rule, kinematic hardening rule of Phillips type, a modified form of the bounding surface model for modelling the cyclic behavior, and the usage of a proposed form for evaluating the plastic modulus for anisotropic materials. Previous experimental data have been used for the evaluation of the yield surface parameters as well as those for the determination of the plastic modulus. The stress-strain results generated from the model have then been compared with those from the experiments. The behavior of the model under certain simulated cyclic loading situations has also been presented.     

金属材料在复杂应力状态下的塑性流动特性及本构模型

工程应用中,金属材料和结构往往处于复杂应力状态。材料的塑性行为会受到应力状态的影响,要精确描述材料在复杂应力状态下的塑性流动行为,必须在本构模型中考虑应力状态效应的影响。然而,由于在动态加载下材料的应变率效应和应力状态效应相互耦合、难以分离,给应力状态效应的研究和模型的建立造成很大困难。通过对Ti-6Al-4V钛合金材料开展不同加载条件下的力学性能测试,提出了一个包含应力三轴度和罗德角参数影响的新型本构模型,并通过VUMAT用户子程序嵌入ABAQUS/Explicit软件。分别采用新提出的塑性模型和Johnson-Cook模型对压剪复合试样的动态实验进行了数值模拟。结果表明,新模型不仅在对材料本构曲线的拟合方面具有较强的优势,而且由该模型所得到的透射脉冲和载荷-位移曲线均更加准确。因此,该模型能够更精确地描述和预测金属材料在复杂应力状态下的塑性流变行为。     

Quasi-static bubble in a yield stress fluid: elasto-plastic model

The influence of the yield stress of Carbopol® gel dispersions on the behaviour of quasi-static bubbles was investigated. Many fluids, from many different industrial fields, have yield stress behaviour. Most of them contain gas bubbles. To study bubble behaviour in such suspensions, a transparent model fluid (dispersion of Carbopol® in water) was used. The experimental device allowed to quasi-statically increase bubble internal pressure with small pressure step to reach a maximum target internal pressure and the pressure setpoint was inverted to return to the initial pressure. Hysterical behaviour of the bubbles was highlighted as they did not regain their initial shape because of yield stress. We show that the rheological behaviour is related to the internal pressure, bubble geometry and yield stress in quasi-static conditions. A modification of the Laplace law depending on the yield stress of the fluid and bubble sphericity was proposed and validated.     

K0超固结结构性黏土的本构模型

K₀固结黏土在自然界广泛分布, 其通常同时具有超固结性与天然结构性, 而K₀超固结性又与K₀正常固结性质存在很大差异. 为了有效的描述K₀超固结性质, 在结构性模型基础上, 做了如下三点改进, 使得原模型拓展为同时考虑K₀超固结特性与天然结构性影响的本构模型. (1)引入相对应力比来描述屈服面, 并引入初始各向异性转轴参量ξ来表达初始各向异性对屈服面在p-q空间的位置影响. (2)基于给定的屈服面方程, 推导得到变相应力比参量, 并将变相应力比引入到统一硬化参数中, 利用统一硬化参数可以有效描述初始各向异性固结黏土在剪切加载下的剪缩与剪胀, 应变硬化及软化现象. (3)引入反映结构性胶结强度性质的胶结参量p_e, 并给出p_e随塑性偏应变的衰减演化方程, 利用胶结参量可描述结构性黏土的剪胀特性. 预测与试验结果对比表明, 所提的K₀超固结结构性模型可有效描述K₀超固结黏土的刚度提高效应, 黏土的包辛格效应, 结构性黏土胶结强度的损失现象以及结构性黏土的应变软化现象. 证明了所提模型的适用性以及合理性.      

Using instrumental inertia in controlled stress rheometry

We describe a new method for characterizing the non-linear behavior of complex fluids at both small and large deformations. For creep measurements, we use the coupling between the instrumental inertia and the material‘s elasticity to follow the rheological behavior of a solution of iota carrageenan both above and below the yield stress. It is shown that this coupling selectively excites one particular frequency of the relaxation spectrum. An analytical calculation is used to quantify the non-linear behavior near the yield stress. The “free“ oscillations observed during the first few seconds allow us to choose the most appropriate mechanical model. Comparison with experiment shows that even above the yield stress, a linear model can still give independently reliable information about the changes in each element of the mechanical model. A comparison of free and forced oscillations in controlled stress rheometry shows both experimentally and theoretically the conditions under which the use of free oscillations is advantageous. Received: 4 September1997 Accepted: 13 January 1998     

平面幂律塑性Ⅰ型裂纹尖端应力场全解

在航空航天、船舶、石油管道和核电等领域,服役结构或部件在长期极端条件下运行,不可避免地会产生裂纹,因此,为研究含裂纹结构的准静态断裂行为,必须了解裂纹尖端附近区域的应力应变场特点.对于幂律材料裂纹构元,研究平面应变和平面应力条件下Ⅰ型裂纹尖端应力场的解析分布.基于能量密度等效和量纲分析,推导了能量密度中值点代表性体积单元(representative volume element, RVE)的等效应力解析方程,并定义其为应力因子,进而针对有限平面应变和平面应力紧凑拉伸(compact tension, CT)试样和单边裂纹弯曲(single edge bend, SEB)试样,以应力因子作为应力特征量,并构造用于表征裂尖等效应力等值线的蝶翅轮廓式和扇贝轮廓式三角特殊函数,提出描述幂律塑性条件下平面I型裂纹尖端应力场的半解析模型.该半解析模型形式简单,对CT和SEB试样的裂尖应力场的预测结果与有限元分析的结果比较表明,两者之间均密切吻合,模型公式可直接用于预测Ⅰ型裂纹尖端应力分布,方便于断裂安全评价和理论发展.     

复杂加载下混凝土的弹塑性本构模型

混凝土材料在不同应力路径下或复杂加载条件下会表现出差异性显著的应力应变关系,在小幅循环加载条件下,其应力应变关系会表现出类似于弹性变形的滞回曲线.在不同应力水平下,混凝土的应力应变关系以及破坏特性都具有静水压力相关特点,即随着静水压力增大,各向异性强度特性弱化.此外,混凝土受压及受拉破坏机理不同,因而对应于混凝土硬化损伤亦有不同,即可分为受压硬化损伤,受拉硬化损伤及两者的混合硬化损伤类型.基于Hsieh模型,对该模型进行了三点改进.（1）针对小幅循环加载下混凝土无塑性变形的试验规律,而模型中在应力水平较低的循环加载条件下始终存在塑性变形的预测问题,采用在边界面模型框架下,设置了应力空间的弹性域,初始屈服面与后续临界状态屈服面几何相似的假定.（2）基于广义非线性强度准则将原模型采用变换应力方法将其推广为三维弹塑性本构模型,采用变换后模型可合理的考虑不同应力路径对于子午面以及偏平面上静水压力效应形成的影响,并避免了边界面应力点奇异问题.（3）分别对拉压两种加载损伤模式建议了相应的硬化参数表达式,可分别用于描述上述加载中产生的应变软化及强度退化行为.基于多种加载路径模拟表明:所建立的三维弹塑性本构模型可合理地用于描述混凝土的一般应力应变关系特性.      

Plastic spin associated with a non-normality theory of plasticity

A plasticity model using a vertex-type plastic flow rule on a smooth yield surface for an anisotropic solid has been proposed recently. This model is here completed by incorporating the effect of plastic spin. Simple shear with a large shear strain is one of the hardest tests on finite strain anisotropic plasticity models, and it is here shown which plastic spin expression is needed to avoid unrealistic oscillatory behavior of the shear stress under large shear strains. The idea of using non-normality with a smooth yield surface originates from a recent proposal of using an abrupt strain path change to determine the subsequent yield surface shape. For this method both polycrystal plasticity calculations and experiments have shown a vertex-type response on the apparently smooth yield surface.