考虑缠结效应的超弹性本构模型

经典熵弹性模型, 如 Neo-Hookean模型和Arruda-Boyce八链模型, 被广泛应用于预测橡胶等软材料的超弹性力学行为. 然而, 大量实验结果也显示仅采用一套模型参数, 这类模型不能同时准确地描述橡胶在多种加载模式下的应力响应. 为了克服上述模型的不足, 本文在熵弹性的模型基础上引入缠结约束效应. 微观上, 采用Langevin统计模型来表征熵弹性变形自由能, 通过管模型(tube model)引入缠结约束自由能, 并基于仿射假设, 建立微观变形与宏观变形之间的映射关系. 在宏观上, 所建立的超弹性模型的Helmholtz自由能同时包含熵弹性和缠结约束两部分, 其中熵弹性自由能与经典的Arruda-Boyce八链模型一致, 依赖于柯西-格林应变张量的第一不变量, 而缠结约束自由能依赖于柯西-格林应变张量的第二不变量. 与文献中的实验结果对比发现, 该三参数模型能准确地预测实验中所测得的橡胶材料在单轴拉伸、纯剪切和等双轴拉伸变形条件下的应力响应, 也能较好地描述不同预拉伸比条件下双轴拉伸实验结果. 最后, 本文比较了所建立的基于应变不变量的缠结约束模型与文献中相关的缠结约束模型在多种加载模式下自由能的异同. 总的来说, 本文所建立的本构理论能准确模拟橡胶等软材料的大变形力学行为, 对其工程应用有促进作用.      

Physical mechanism of the compressive response of F-actin networks: significance of crosslinker unbinding events

A model of crosslinker unbinding is implemented in a highly coarse-grained granular model of F-actin cytoskeleton. We employ this specific granular model to study the mechanisms of the compressive responses of F-actin networks. It is found that the compressive response of F-actin cytoskeleton has dependency on the strain rate. The evolution of deformation energy in the network indicates that crosslinker unbinding events can induce the remodelling of F-actin cytoskeleton in response to external loadings. The internal stress in F-actin cytoskeleton can efficiently dissipate with the help of crosslinker unbinding, which could lead to the spontaneous relaxation of living cells.     

Dilatational plastic constitutive equations and their application to ductile fracture analysis of axisymmetric specimens

In many ductile materials voids nucleate and grow under large strain and triaxial stress, which yield volumetric plastic expansion. A constitutive equation is presented, which accounts for this plastic dilatancy. The plastic moduli involved in this equation can be calibrated by using necking tests of axisymmetric bars, void model analysis and computer simulation. To verify the rationality of such a constitutive equation and adjust the values of plastic moduli, the constitutive equation with its moduli to be determined is applied to analyse the ductile fracture behaviour of axisymmetric bars.     

Shear band analysis and shear moduli calibration

Strain localization is a well known phenomenon, generally associated with plastic deformation and rupture in solids, especially in geomaterials. In this process, deformation is observed to concentrate in narrow zones called shear bands. This phenomenon has been studied extensively in the last 20 years by different researchers, experimentally, theoretically and numerically. A criterion for the onset of localization can be predicted solely on the basis of the constitutive law of the material, using the so-called shear band analysis. This criterion gives the critical orientation, and the critical stress state and strain for a given loading history. An important point, already stressed by Vardoulakis in 1980, is that in particular, out-of-axes shear moduli play a central role in the criterion. These are the moduli involved in the response to a deviatoric stress increment with principal axes oriented at 45° from total stress principal axes. Out-of-axes shear moduli are difficult parameters to calibrate; common tests, with fixed principal stress and strain directions, do not provide any information on these moduli, as long as they remain homogeneous. Still, real civil engineering and environmental problems are definitely not simple axisymmetric triaxial tests; practical modeling involves complex stress paths, and need complex parameters to be calibrated. Only special tests, like compression–torsion on hollow cylinder tests, or even more complex tests can be used for shear moduli calibration. However, shear band initiation in homogeneous, fixed-axes tests does activate out-of-axes shear. Hence, it is natural that shear band analysis makes shear moduli enter into the analysis.Then, a typical inverse analysis approach can be used here: experimental observation of strain localization in triaxial tests can be used together with a proper shear band analysis for the model considered, in order to determine out-of-axes shear moduli.This approach has been used for a stiff marl in the framework of a calibration study on a set of triaxial tests. The steps of the method are presented, and the bifurcation surface in the stress space is exhibited.     

Constitutive equations from Gaussian slip-link network theories in polymer melt rheology

Allowing for flow-dependent slip in the junctions of a temporary junction network, we derive the constitutive equations of temporary slip-link networks. The stress tensor is determined by three material functions, namely, the time-dependent linear-viscoelastic memory function, and two strain-dependent functions describing slip and disentanglement of network strands. Slip and disentanglement are related via a mass balance for network strands.By specifying slip and disentanglement, the constitutive equations of Lodge, Wagner, Doi-Edwards, and Marrucci are shown to be special temporary slip-link constitutive equations. To demonstrate the predictive power of temporary slip-link network theories, we compare predictions and extensional flow data with step change in flow direction.Dedicated to Professor Arthur S. Lodge on the occasion of his 70th birthday and his retirement from the University of Wisconsin.     

柔性硅泡沫薄片的短时松弛行为研究

根据多孔硅泡沫材料的单轴压缩和应力松弛实验,利用最小二乘法的LM法拟合得到硅泡沫材料的本构关系.基于上述模型开展了组合结构中多孔泡沫薄片应力松弛行为的数值模拟,得到了短时松弛过程中硅泡沫结构件的应力变化规律.     

Multi-cycle deformation of silicone elastomer: observations and constitutive modeling with finite strains

Observations are reported on a medical grade of silicone elastomer in uniaxial tensile tests up to breakage of specimens, short-term relaxation tests, and cyclic tests with monotonically increasing maximum elongation ratios. Experimental data in cyclic tests demonstrate the fading memory phenomenon: stress–strain diagrams for two specimens with different deformation histories along the first n?1 cycles and coinciding loading programs for the other cycles become identical starting from the nth cycle. A constitutive model is developed in cyclic viscoplasticity of elastomers with finite strains, and its adjustable parameters are found by fitting the experimental data. Ability of the stress–strain relations to predict the mechanical response in cyclic tests with various deformation programs is confirmed by numerical simulation.     

Self-consistent modeling of large plastic deformation, texture and morphology evolution in semi-crystalline polymers

A self-consistent model for semi-crystalline polymers is proposed to study their constitutive behavior, texture and morphology evolution during large plastic deformation. The material is considered as an aggregate of composite inclusions, each representing a stack of crystalline lamellae with their adjacent amorphous layers. The deformation within the inclusions is volume-averaged over the phases. The interlamellar shear is modeled as an additional slip system with a slip direction depending on the inclusion's stress. Hardening of the amorphous phase due to molecular orientation and, eventually, coarse slip, is introduced via Arruda-Boyce hardening law for the corresponding plastic resistance. The morphology evolution is accounted for through the change of shape of the inclusions under the applied deformation gradient. The overall behavior is obtained via a viscoplastic tangent self-consistent scheme. The model is applied to high density polyethylene (HDPE). The stress-strain response, texture and morphology changes are simulated under different modes of straining and compared to experimental data as well as to the predictions of other models.     

Finite viscoplasticity of amorphous glassy polymers in the logarithmic strain space

The paper outlines a new constitutive model and experimental results of rate-dependent finite elastic–plastic behavior of amorphous glassy polymers. In contrast to existing kinematical approaches to finite viscoplasticity of glassy polymers, the formulation proposed is constructed in the logarithmic strain space and related to a six-dimensional plastic metric. Therefore, it a priori avoids difficulties concerning with the uniqueness of a plastic rotation. The constitutive framework consists of three major steps: (i) A geometric pre-processing defines a total and a plastic logarithmic strain measures determined from the current and plastic metrics, respectively. (ii) The constitutive model describes the stresses and the consistent moduli work-conjugate to the logarithmic strain measures in an analogous structure to the geometrically linear theory. (iii) A geometric post-processing maps the stresses and the algorithmic tangent moduli computed in the logarithmic strain space to their nominal, material or spatial counterparts in the finite deformation space. The analogy between the formulation of finite plasticity in the logarithmic strain space and the geometrically linear theory of plasticity makes this framework very attractive, in particular regarding the algorithmic implementation. The flow rule for viscoplastic strains in the logarithmic strain space is adopted from the celebrated double-kink theory. The post-yield kinematic hardening is modeled by different network models. Here, we compare the response of the eight chain model with the newly proposed non-affine micro-sphere model. Apart from the constitutive model, experimental results obtained from both the homogeneous compression and inhomogeneous tension tests on polycarbonate are presented. Besides the load–displacement data acquired from inhomogeneous experiments, quantitative three-dimensional optical measurements of the surface strain fields are carried out. With regard to these experimental data, the excellent predictive quality of the theory proposed is demonstrated by means of representative numerical simulations.     

PBX炸药本构关系的实验研究

在SHPB装置上开展了某典型PBX炸药的单轴压缩、间接拉伸实验。利用入射波整形技术,结合 半导体应变片和石英晶体联合测试的方法,实现了PBX炸药的应力平衡和恒应变率加载,得到了不同应变率 下某PBX炸药的压缩应力应变曲线,初步建立了该炸药的修正Sargin唯象模型;开展了3种PBX炸药的动 态巴西实验,通过高速摄影获得了炸药的破坏过程,结合数字相关技术,获得了试样裂纹附近的应变场分布, 初步建立了描述3种炸药动态拉伸行为的修正Johnson-Cook模型。模型曲线与实验结果符合较好。     

Cyclic stress-strain behaviour of alloy AA6060 T4, part II: Biaxial experiments and modelling

Laboratory tests have been conducted to investigate the inelastic behaviour of aluminium alloy AA6060 T4 subjected to non-proportional cyclic loading. The results of four tests with variable strain path shapes and strain amplitudes are reported in this paper. The tests were carried out by applying combined axial force and torque to thin-walled tubular specimens, using effective strain amplitudes in the range 0.4–0.8%. Major emphasis has been put on the two important material properties: plastic anisotropy and influence of strain range and strain path shapes on cyclic hardening. A constitutive model for cyclic plasticity is used to predict the stress response of the alloy for the non-proportional strain paths applied in the experiments. The model adopts a quadratic yield function and multi-component non-linear isotropic and kinematic hardening rules to describe plastic anisotropy, the shape of the hysteresis loops and the evolution of cyclic hardening. Good agreement is obtained between the physical and correlated stress response of the alloy.     

Anisotropic Mullins stress softening of a deformed silicone holey plate

Rubber like materials parts are designed using finite element code in which more and more precise and robust constitutive equations are implemented. In general, constitutive equations developed in literature to represent the anisotropy induced by the Mullins effect present analytical forms that are not adapted to finite element implementation. The present paper deals with the development of a constitutive equation that represents the anisotropy of the Mullins effect using only strain invariants. The efficiency of the modeling is first compared to classical homogeneous experimental tests on a filled silicone rubber. Second, the model is tested on a complex structure. In this aim, a silicone holey plate is molded and tested in tension, its local strain fields are evaluated by means of digital image correlation. The experimental results are compared to the simulations from the constitutive equation implemented in a finite element code. Global measurements (i.e. force and displacement) and local strain fields are successfully compared to experimental measurements to validate the model.     

On the Correlation of FEM and Experiments for Hyperelastic Elastomers

Correlation of modern finite element methods (FEM) with advanced experimental techniques for elastomers, biomedical materials, and living organs requires study and modification of the behavior of these materials. In this study, the mechanical behavior of a commonly-used elastomer, silicone rubber, which provides excellent biocompatibility, was examined under different applied loading configurations, and large deformations were investigated through both experiment and simulation. The stress-strain behaviors of silicone rubber were tested, using multiple homogeneous experiments, including uniaxial extension and equibiaxial tension, the load-apex displacement response, and digitized deformed shapes of two of the most-used structures for nonlinear hyperelasticity—the inflation of a clamped circular membrane, and indentation of the membrane by a spherical indenter. Uniaxial and equibiaxial data were evaluated simultaneously, characterized by various constitutive models for implementation in the FE simulation. These constitutive models examined the prediction of the FE simulations for the inflation and indentation tests in comparison to the results of experiments at various load-apex displacement levels. The results showed that the constitutive models calibrated with the uniaxial and equibiaxial tests, predicted nearly the same results as the actual experimental results, particularly for the applied loads that generated moderate strain. However, when the FE simulations based on the constitutive models were adjusted, employing only uniaxial or equibiaxial tests, they predicted different results, where the degree of their correlations with experimental results was incomplete or in some states simply poor. The simulations suggested that the inverse FE procedure should not be restricted to the choice of material models, while more attention should be given to the choice of ranges of deformation.     

高应变率下航空透明聚氨酯的动态本构模型

采用低阻抗分离式霍普金森压杆对航空透明聚氨酯进行了高应变率下的动态力学性能测试,得到的应力应变曲线表现出了显著的非线性黏弹性特征。基于本构理论和实验数据,构建了航空透明聚氨酯的松弛时间应变相关的超黏弹性本构形式。该本构模型由2部分组成:一部分表征准静态下的超弹性行为,另一部分描述非线性应变率的相关特性。利用超黏弹性本构模型对不同应变率下航空透明聚氨酯的动态应力应变曲线进行拟合,拟合曲线与实验曲线一致性良好。     

Modelling and experimental characterisation of the Lüders strain in complex loaded ferritic steel compact tension specimens

The use of 3D digital image correlation (DIC) has been used to capture the Lüders strains in a low carbon ferritic steel. Results were used to calibrate and compare with finite element (FE) results based on a constitutive plasticity model, capable of yield drop behaviour and therefore Lüders strains, by Zhang et al. (2001). Tensile tests were carried out at several strain rates to characterise the material behaviour. The results of these tests were used to fit parameters in the constitutive plasticity model. The FE model was then tested on a complex loading situation of in-plane compression of a compact tension (CT) specimen. The FE model predicts the shape and formation of the Lüders bands well. This FE model, using Zhang’s constitutive plasticity model, was used to predict the residual stress profile to compare with standard elastic–plastic isotropic hardening models with no yield point. The yield point reduced both the predicted peak tensile stress, at the notch root, and the amount of plastic strain. In regions where the plastic strain was of a similar size to the Lüders strain the stress profiles were perturbed from flat profiles predicted by the standard elastic–plastic hardening models.     

基于深度学习和细观力学的颗粒材料本构关系研究

颗粒材料的本构关系对岩土工程等众多领域至关重要. 不同于传统的唯象本构理论, 本文基于机器学习模型探索了一种细观力学理论指导下的数据驱动型颗粒材料本构关系预测方法. 根据Vogit均质化假设, 建立了小应变条件下颗粒材料应力?应变解析关系, 此关系唯一地确定了一组与颗粒材料本构行为相关的细观组构变量. 这些变量与反应颗粒材料宏观性质的主应力和主应变信息通过一系列离散元三轴压缩数值试验获得. 考虑到细观组构变量为内变量, 不能直接作为本构模型的输入. 本文基于有向图方法将颗粒材料微观结构信息隐式地包含在应力?应变的预测当中, 并采用门控循环单元(GRU)循环神经网络作为基础深度学习模型描述有向图中结点之间的映射关系. 通过将有向图从目标节点沿源节点展开, 整个应力?应变预测模型可由两个神经网络分别训练并组装而成. 将训练后的深度学习模型在全新的数据集上进行测试, 结果表明该训练策略能有效捕捉到颗粒材料在常规三轴任意加卸载, 等中主应力系数b的真三轴加载, 和等平均有效应力p的真三轴加卸载等复杂多轴加载工况下的应力?应变响应关系, 模型具有良好的内插和外推预测能力. 考虑到深度学习模型捕捉颗粒材料力学响应的能力及其开放式学习的特点, 充分结合数据驱动方法和理论本构模型可能是颗粒材料本构研究的一个重要方向.      

An investigation on microscopic and macroscopic stability phenomena of composite solids with periodic microstructure

An analysis of the effects of microscopic instabilities on the homogenized response of heterogeneous solids with periodic microstructure and incrementally linear constitutive law is here carried out. In order to investigate the possibility to obtain a conservative prediction of microscopic primary instability in terms of homogenized properties, novel macroscopic constitutive stability measures are introduced, corresponding to the positive definiteness of the homogenized moduli tensors relative to a class of conjugate stress–strain pairs.Numerical simulations, addressed to hyperelastic microstructural models representing cellular solids and reinforced composites, are worked out through the implementation of an innovative one-way coupled finite element formulation able to determine sequentially the principal equilibrium solution, the incremental equilibrium solutions providing homogenized moduli and the stability eigenvalue problem solution, for a given monotonic macrostrain path. Both uniaxial and equibiaxial loading conditions are considered.The exact microscopic stability region in the macrostrain space, obtained by taking into account microstructural details, is compared with the macroscopic stability regions determined by means of the introduced macroscopic constitutive measures. These results highlight how the conservativeness of the adopted macroscopic constitutive stability measure with respect to microscopic primary instability, strictly depends on the type of loading condition (tensile or compressive) and the kind of microstructure.     

多孔硅橡胶材料有限元计算的切线本构矩阵

为进行多孔硅橡胶材料非线性有限元计算，针对用伸长比表示的应变能函数，推导了多孔硅橡胶材料非线性有限元计算中切线本构关系的具体形式。并采用便于编程的矩阵记法，给出了非线性有限元计算中克希荷夫应力和格林应变之间的切线关系式。     

Membrane Curvatures and Stress-strain Full Fields of Axisymmetric Bulge Tests from 3D-DIC Measurements. Theory and Validation on Virtual and Experimental results

The bulge test is mostly used to analyze equibiaxial tensile stress state at the pole of inflated isotropic membranes. Three-dimensional digital image correlation (3D-DIC) technique allows the determination of three-dimensional surface displacements and strain fields. In this paper, a method is proposed to determine also the membrane stress tensor fields for in-plane isotropic materials, independently of any constitutive equation. Stress-strain state is then known at any surface point which enriches greatly experimental data deduced from the axisymmetric bulge tests. Our method consists, first in calculating from the 3D-DIC experimental data the membrane curvature tensor at each surface point of the bulge specimen. Then, curvature tensor fields are used to investigate axisymmetry of the test. Finally in the axisymmetric case, membrane stress tensor fields are determined from meridional and circumferential curvatures combined with the measurement of the inflating pressure. Our method is first validated for virtual 3D-DIC data, obtained by numerical simulation of a bulge test using a hyperelastic material model. Afterward, the method is applied to an experimental bulge test performed using as material a silicone elastomer. The stress-strain fields which are obtained using the proposed method are compared with results of the finite element simulation of this overall bulge test using a neo-Hookean model fitted on uniaxial and equibiaxial tensile tests.