基于单轴数据决定橡胶类材料多轴刚性化有界超弹性势

指出表征橡胶类材料应变刚化效应的现有超弹性模型涉及应变能无穷发散困难,提出新方法解决该困难。基于对数应变不变量的多轴扩张和多轴匹配步骤,建议直接构造橡胶类材料大变形弹性势的显式直接方法。该方法从单轴应力-应变关系直接得到多轴弹性势,所得结果避免了现有各方法决定待定参数组的复杂数值计算,能够准确描述应变刚性化效应,且给出有界弹性应变能,从而避免了前述发散困难。数值结果表明,从单轴数据所得到的弹性势可同时很好的拟合平面应变拉伸(剪切)数据以及等双轴拉伸数据。     

OBTAINING MULTI-AXIAL ELASTIC POTENTIALS FOR RUBBER-LIKE MATERIALS VIA AN EXPLICIT,EXACT APPROACH BASED ON SPLINE INTERPOLATION

An explicit, exact approach is proposed to obtain multi-axial elastic potentials for isotropic rubber-like materials undergoing large incompressible deformations. By means of two direct, explicit procedures, this approach reduces the problem of determining multi-axial poten- tials to that of determining one-dimensional elastic potentials. To this end, two one-dimensional potentials for uniaxial case and simple shear case are respectively determined via spline inter- polation and, then, the two potentials are extended to generate a multi-axial elastic potential using a novel method based on certain logarithmic invariants. Eventually, each of the multi-axial potentials will exactly match the finite strain data from four benchmark tests.     

Dynamic multi-axial behavior of shape memory alloy nanowires with coupled thermo-mechanical phase-field models

The objective of this paper is to provide new insight into the dynamic thermo-mechanical properties of shape memory alloy (SMA) nanowires subjected to multi-axial loadings. The phase-field model with Ginzburg–Landau energy, having appropriate strain based order parameter and strain gradient energy contributions, is used to study the martensitic transformations in the representative 2D square-to-rectangular phase transformations for FePd SMA nanowires. The microstructure and mechanical behavior of martensitic transformations in SMA nanostructures have been studied extensively in the literature for uniaxial loading, usually under isothermal assumptions. The developed model describes the martensitic transformations in SMAs based on the equations for momentum and energy with bi-directional coupling via strain, strain rate and temperature. These governing equations of the thermo-mechanical model are numerically solved simultaneously for different external loadings starting with the evolved twinned and austenitic phases. We observed a strong influence of multi-axial loading on dynamic thermo-mechanical properties of SMA nanowires. Notably, the multi-axial loadings are quite distinct as compared to the uniaxial loading case, and the particular axial stress level is reached at a lower strain. The SMA behaviors predicted by the model are in qualitative agreements with experimental and numerical results published in the literature. The new results reported here on the nanowire response to multi-axial loadings provide new physical insight into underlying phenomena and are important, for example, in developing better SMA-based MEMS and NEMS devices     

Direct determination of multi-axial elastic potentials for incompressible elastomeric solids: an accurate,explicit approach based on rational interpolation

With a novel approach based on certain logarithmic invariants, we demonstrate that a multi-axial elastic potential for incompressible, isotropic rubber-like materials may be obtained directly from two one-dimensional elastic potentials for uniaxial case and simple shear case, in a sense of exactly matching finite strain data for four benchmark tests, including uniaxial extension, simple shear, bi-axial extension, and plane-strain extension. As such, determination of multi-axial elastic potentials may be reduced to that of two one-dimensional elastic potentials. We further demonstrate that the latter two may be obtained by means of rational interpolating procedures for uniaxial data and shear data displaying strain-stiffening effects. Numerical examples are presented in fitting Treloar’s data and other data.     

类橡胶材料变形直到破坏的显式本构模型

本文通过直接、显式的方法提出一个多轴可压缩应变能弹性势来模拟类橡胶材料受载荷直到软化破坏的变形行为．首先，我们提出一个多轴可压缩应变能函数；其次，通过特定的不变量，该多轴应变能函数在单轴拉伸，平面应变和等双轴拉伸三个基准实验的情况下，可以退化为各自的单轴形函数形式；再次，我们显式给出带有软化破坏特性的形函数；最后，模型结果和试验数据可以精确匹配，同时可以预测材料临近破坏以后，接下来的变形行为．     

显式方法精确模拟形状记忆聚合物热力学行为

通过构建一个热耦合的多轴可压缩应变能函数,得到应力-应变、应力-温度和应变-温度之间的函数关系,建立形状记忆聚合物的本构方程.本文引入三个基于对数应变的不变量使得模型（i）可以模拟可压缩情况;（ii）适用于单轴拉伸和等双轴拉伸至少两个基准实验;（iii）多轴有效.通过显式方法（i）给出自由能和熵的具体表达,证明模型热力学定律;（ii）给出应变-应力,温度-应力以及,温度-应变的形函数具体表达.多轴模型在特定的情况下可以自动退化到各自的单轴情况. 通过调节形函数的参数,最终得到的模型结果和实验结果能够精确匹配.新方法建立的本构模型得到的结果能更加准确地指导形状记忆聚合物的工程设计。     

显式方法精确模拟形状记忆聚合物热力学行为

通过构建一个热耦合的多轴可压缩应变能函数,得到应力-应变、应力-温度和应变-温度之间的函数关系,建立形状记忆聚合物的本构方程.本文引入三个基于对数应变的不变量使得模型（i）可以模拟可压缩情况;（ii）适用于单轴拉伸和等双轴拉伸至少两个基准实验;（iii）多轴有效.通过显式方法（i）给出自由能和熵的具体表达,证明模型热力学定律;（ii）给出应变-应力,温度-应力以及,温度-应变的形函数具体表达.多轴模型在特定的情况下可以自动退化到各自的单轴情况. 通过调节形函数的参数,最终得到的模型结果和实验结果能够精确匹配.新方法建立的本构模型得到的结果能更加准确地指导形状记忆聚合物的工程设计。     

显式模拟类橡胶材料Mullins效应滞回圈

通过显式、直接的方法提出一个多轴可压缩应变能函数,用来模拟类橡胶材料在加载——卸载作用下,由于Mullins效应而产生的应力——应变滞回圈. 本文的创新点在于将表征能量耗散的变量引入到应变能函数.新的弹性势具有以下两个特点:第一,在加载情况下,新引入的变量不会对弹性势产生任何影响,因此,只要给出合适的形函数显式表达,3个基准实验,包括单轴拉伸和压缩,等双轴拉伸和压缩,以及平面应变,都可精确模拟;第二,新引入的变量在卸载情况下将被激活.在不同的卸载应力下,变量将发生改变,从而影响弹性势,使其最终产生不同的应力——应变关系卸载曲线,与对应的加载曲线共同构成应力——应变滞回圈.通过对Mullins效应实验数据进行分析和研究,得出了卸载形函数在不同卸载应力下变化的规律,并预测不同卸载应力下的应力——应变关系.最后,我们将得到精确匹配实验数据的数值模拟结果,从而证明本文方法不仅可以精确匹配至少3个基准实验,还可以模拟和预测类橡胶材料在加载——卸载作用下由于Mullins效应而产生的滞回圈.      

The multi-axial deformation behavior of bulk metallic glasses at high homologous temperatures

In this work, we study the behavior of a recently-developed Lanthanum-based bulk metallic glass under uniaxial and multi-axial stress-states using the constitutive model developed by Thamburaja and Ekambaram (2007). The material parameters in the constitutive model are fitted to match the stress–strain responses obtained from a set of simple compression experiments conducted at temperatures within the supercooled liquid region under a variety of strain rates spanning approximately three decades. With the material parameters calibrated, we show that the aforementioned constitutive model is able to accurately predict the force vs. displacement responses of representative experiments conducted under multi-axial stress-states at temperatures within the supercooled liquid region, namely three-point bending and the superplastic forming of a miniature gear component. In particular, the evolution of the specimen geometry during the deformation under multi-axial loading conditions are also well-predicted by the constitutive model.     

多轴加载下混凝土细观破坏模拟的强度准则探讨

实际工程结构中混凝土材料大多处于双轴或三轴的复杂应力状态,已有的细观力学数值研究工作大多针对单轴加载问题,对于双轴或者三轴加载条件下混凝土破坏模拟的研究相对较少。复杂受力条件下的混凝土材料破坏模拟中,细观组分强度准则选取的合理与否将成为混凝土破坏模式及宏观力学性能数值研究准确和成功与否的关键。本文旨在探讨单轴强度准则,如最大拉应变准则在多轴加载条件下混凝土破坏过程研究中运用的合理性。鉴于此,首先在细观尺度上建立了混凝土试件的二维随机骨料模型,分别采用弹性损伤本构关系模型及塑性损伤本构关系模型来描述细观组分(即砂浆基质)的力学性能,对双轴加载条件下混凝土的细观破坏过程进行数值模拟,对比了单轴强度准则和多轴强度准则下混凝土试件破坏路径及宏观应力-应变关系的差异。数值结果表明,简单的单轴强度准则难以反映双轴加载下混凝土内部应力状态的复杂性,不宜采用单轴强度准则来描述多轴加载下混凝土的破坏行为。     

Optimization of the strain energy density in linear anisotropic elasticity

The problem considered here is that of extremizing the strain energy density of a linear anisotropic material by varying the relative orientation between a fixed stress state and a fixed material symmetry. It is shown that the principal axes of stress must coincide with the principal axes of strain in order to minimize or maximize the strain energy density in this situation. Specific conditions for maxima and minima are obtained. These conditions involve the stress state and the elastic constants. It is shown that the symmetry coordinate system of cubic symmetry is the only situation in linear anisotropic elasticity for which a strain energy density extremum can exist for all stress states. The conditions for the extrema of the strain energy density for transversely isotropic and orthotropic materials with respect to uniaxial normal stress states are obtained and illustrated with data on the elastic constants of some composite materials. Not surprisingly, the results show that a uniaxial normal stress in the grain direction in wood minimizes the strain energy in the set of all uniaxial stress states. These extrema are of interest in structural and material optimization.     

A multi-axial, multimechanism based constitutive model for the comprehensive representation of the evolutionary response of SMAs under general thermomechanical loading conditions

We present a fully general, three dimensional, constitutive model for Shape Memory Alloys (SMAs), aimed at describing all of the salient features of SMA evolutionary response under complex thermomechanical loading conditions. In this, we utilize the mathematical formulation we have constructed, along with a single set of the model’s material parameters, to demonstrate the capturing of numerous responses that are experimentally observed in the available SMA literature. This includes uniaxial, multi-axial, proportional, non-proportional, monotonic, cyclic, as well as other complex thermomechanical loading conditions, in conjunction with a wide range of temperature variations. The success of the presented model is mainly attributed to the following two main factors. First, we use multiple inelastic mechanisms to organize the exchange between the energy stored and energy dissipated during the deformation history. Second, we adhere strictly to the well established mathematical and thermodynamical requirements of convexity, associativity, normality, etc. in formulating the evolution equations governing the model behavior, written in terms of the generalized internal stress/strain tensorial variables associated with the individual inelastic mechanisms. This has led to two important advantages: (a) it directly enabled us to obtain the limiting/critical transformation surfaces in the spaces of both stress and strain, as importantly required in capturing SMA behavior; (b) as a byproduct, this also led, naturally, to the exhibition of the apparent deviation from normality, when the transformation strain rate vectors are plotted together with the surfaces in the space of external/global stresses, that has been demonstrated in some recent multi-axial, non-proportional experiments.     

Constitutive relations describing creep deformation for multi-axial time-dependent stress states

A theory of primary and secondary creep deformation in metals is presented, which is based upon the concept of tensor internal state variables and the principles of continuum mechanics and thermodynamics. The theory is able to account for both multi-axial and time-dependent stress and strain states. The wellknown concepts of elastic, anelastic and plastic strains follow naturally from the theory. Homogeneous stress states are considered in detail and a simplified theory is derived by linearizing with respect to the internal state variables. It is demonstrated that the model can be developed in such a way that multi-axial constant-stress creep data can be presented as a single relationship between an equivalent stress and an equivalent strain. It is shown how the theory may be used to describe the multi-axial deformation of metals which are subjected to constant stress states. The multi-axial strain response to a general cyclic stress state is calculated. For uni-axial stress states, square-wave loading and a thermal fatigue stress cycle are analysed.     

Thermo-mechanical response of nylon 101 under uniaxial and multi-axial loadings: Part I,Experimental results over wide ranges of temperatures and strain rates

A comprehensive study of the thermo-mechanical response of a thermoplastic polymer, nylon 101 is presented. Quasi-static and dynamic compression uniaxial and multi-axial experiments (stress states) were performed at a wide range of strain rates (10⁻⁵ to 5000 s⁻¹) and temperatures (−60 to 177 °C or −76 to 350 °F). The material is found to be non-linearly dependent on strain rate and temperature. The change in volume after plastic deformation is investigated and is found to be negligibly small. The relaxation and creep responses at room temperature are found to be dependent on strain rate and the stress–strain level at which these phenomena are initiated. Total deformation is decomposed into visco-elastic and visco-plastic components; these components have been determined at different levels of deformation. Results from non-proportional uniaxial to biaxial compression, and torsion experiments, are also reported for three different strain rates at room temperature. It is shown that nylon 101 has a response dependent on the hydrostatic pressure.     

Compressive failure of fiber composites under multi-axial loading

This paper examines the compressive strength of a fiber reinforced lamina under multi-axial stress states. An equilibrium analysis is carried out in which a kinked band of rotated fibers, described by two angles, is sandwiched between two regions in which the fibers are nominally straight. Proportional multi-axial stress states are examined. The analysis includes the possibility of bifurcation from the current equilibrium state. The compressive strength of the lamina is contingent upon either attaining a load maximum in the equilibrium response or satisfaction of a bifurcation condition, whichever occurs first. The results show that for uniaxial loading a non-zero kink band angle β produces the minimum limit load. For multi-axial loading, different proportional loading paths show regimes of bifurcation dominated and limit load dominated behavior. The present results are able to capture the beneficial effect of transverse compression in raising the composite compressive strength as observed in experiments.     

Large-deformation properties of wheat dough in uni- and biaxial extension. Part I. Flour dough

Rheological and fracture properties of optimally mixed flour doughs from three wheat cultivars which perform differently in cereal products were studied in uniaxial and biaxial extension. Doughs were also tested in small angle sinusoidal oscillation. In accordance with previously published results the linear region was found to be very small. The rheological properties at small deformations hardly depended on the cultivar. A higher water content of the dough resulted in a lower value for the storage modulus and a slightly higher value for tan . For both uniaxial and biaxial extension a more than proportional increase in stress was found with increasing strain, a phenomenon called strain hardening. In uniaxial extension (i) stresses at a certain strain were higher and (ii) the stress was less dependent on the strain rate than in biaxial extension. This indicates that in elongational flow orientational effects are of large importance for the mechanical properties of flour dough. This conclusion is consistent with published data on birefringence of stretched gluten. Fracture stress and strain increased with increasing deformation rate. The observed time-dependency of fracture properties can best be explained by inefficient transport of energy to the crack tip. Presumably, this is caused by energy dissipation due to inhomogeneous deformation because of friction between structural elements, e.g. between dispersed particles and the network. Differences in the rheological properties at large deformations between the cultivars were observed with respect to (i) stress, (ii) strain hardening, (iii) strain rate dependency of the stress, (iv) fracture properties and (v) the stress difference between uniaxial and biaxial extension.     

Size dependent response of large shell elements under in-plane tensile loading

Large-scale thin-walled structures with a low weight-to-stiffness ratio provide the means for cost and energy efficiency in structural design. However, the design of such structures for crash and impact resistance requires reliable FE simulations. Large shell elements are used in those simulations. Simulations require the knowledge of the true stress–strain response of the material until fracture initiation. Because of the size effects, local material relation determined with experiments is not applicable to large shell elements. Therefore, a numerical method is outlined to determine the effect of element size on the macroscopic response of large structural shell elements until fracture initiation. Macroscopic response is determined by introducing averaging unit into the numerical model over which volume averaged equivalent stress and plastic strain are evaluated. Three different stress states are considered in this investigation: uniaxial, plane strain and equi-biaxial tension. The results demonstrate that fracture strain is highly sensitive to size effects in uniaxial tension whereas in plane strain or equi-biaxial tension size effects are much weaker. In uniaxial and plane strain tension the fracture strain for large shell elements approaches the Swift diffuse necking condition.     

基于颗粒离散元法的连接键应变软化模型及其应用

基于颗粒间的有限接触假设,提出了可表述颗粒间力、力矩传递的连接键模型. 为了表征连接键的塑性、损伤及断裂过程,在连接键中引入了考虑应变软化效应的Mohr-Coulomb 准则及最大拉应力准则. 单一连接键的单向拉伸测试及直剪测试表明了上述连接键应变软化模型的计算精度. 研究了颗粒体系的宏观应变能与颗粒平均配位数的对应关系. 通过计算发现,对于二维颗粒体系,当平均配位数为5 时,颗粒体系的宏观应变能与相同参数下连续介质方法(如有限元等) 计算获得的应变能基本一致. 利用上述连接键应变软化模型对岩石的单轴压缩过程进行了模拟,计算结果表明:岩石单轴压缩的应力应变曲线经历了线性上升段、非线性上升段、非线性下降段及缓变段等4 个阶段,并给出了上述4 个阶段与岩石内部损伤破裂状态的内在联系. 计算结果还表明,随着断裂应变的增大,岩石的破裂模式逐渐由拉剪复合型破裂向单一压剪型破裂转化;随着断裂应变的增大,峰值应力及达到峰值应力时的应变均逐渐增大,但峰值时的破裂度及终态时的破裂度将逐渐减小.