Random vibration of viscoelastic system under broad-band excitations

An analytical method is proposed to study the response of a viscoelastic system with strongly non-linear stiffness force and under broad-band random excitations. The random excitations can be additive, or multiplicative, or both, and they can be stationary or non-stationary with evolutionary spectra. With the proposed method, contributions of the viscoelastic force to both damping and stiffness are taken into account separately, and then the extended version of the stochastic averaging, called the quasi-conservative averaging, is applied to the system to derive the averaged equation of energy envelope. Probability density functions of system responses, such as the total energy, the amplitude, and the state variables, can then be obtained analytically. The accuracy of the method is substantiated by comparing the analytical results with those from Monte Carlo simulations. Effects of parameters in the viscoelastic force and in the non-linear stiffness force on the system responses are also investigated.     

Spectrum of one-dimensional oscillations in the combined stratified medium consisting of a viscoelasticmaterial and a viscous compressible fluid

Two averaged (effective) models corresponding to transverse and longitudinal oscillations in a periodic combined stratified medium are constructed. The medium considered consists of alternating layers of an isotropic viscoelastic material and a viscous compressible fluid. The dependence of the qualitative characteristics of the spectra of the averaged models on the linear dimensions of the stratified medium specimen is investigated.     

Stress analysis of an orthotropic material under diametral compression

The diametral compression test is commonly used to determine the tensile strength of brittle materials. For isotropic materials a simple relation based on specimen geometry and the applied load at failure is used to calculate the tensile strength. Previous to this work the effect of material orthotropy and material orientation on the specimen stress state had not been completely determined. In this study, both isotropic and orthotropic specimens were analyzed using a finite element analysis and experimentally verified by strain gage and photoelastic measurements. Further, this work investigated the effect of the applied load area on the specimen stress state. Results of this work show that there is a significant difference between the theoretical calculations based on the assumption of material isotropy when compared to an orthotropic material. This difference can be as much as 45 percent depending on the degree of orthotropy and the orientation. It was also determined that the applied load area and material orientation significantly influence the specimen stress state. An applied load area of 8 percent of the circumference was found to reduce the stresses in the applied load region.     

Slow flows of yield stress fluids: yielding liquids or flowing solids?

Yield stress fluids (YSF) exhibit strongly non-linear rheological characteristics. As a consequence, they develop original flow features (as compared to simple fluids) under various boundary conditions. This paper reviews and analyzes the characteristics of a series of slow flows (just beyond yielding) under more or less complex conditions (simple shear flow, flow through a cavity, dip-coating, blade-coating, Rayleigh-Taylor instability, Saffman-Taylor instability) and highlights some of their common original characteristics: (i) a transition from a solid regime to a flowing regime which does not correspond to a true “liquid state,” the flow in this regime may rather be seen as a succession of solid states during very large deformation; (ii) a strong tendency to localization of the yielded regions in some small region of the material while the rest of the material undergoes some deformation in its solid state; (iii) the deformation of YSF interface with another fluid, in the form of fingers tending to penetrate the material via a local liquefaction process. Finally, these observations suggest that slow flows of YSF are a kind of extension of plastic flows for very large deformations and without irreversible changes of the structure. This suggests that the field of plasticity and the field of slow flows of YSF could benefit from each other.     

孔洞缺陷对3D打印Ti-6Al-4V合金疲劳试样应力分布的影响

3D打印金属技术因其个性化及可用于加工复杂零件等显著优点,在医用骨植入体领域得到了快速发展,但3D打印金属材料的孔洞缺陷所引起的应力集中现象严重降低了其疲劳强度,限制了3D打印生物金属材料的运用。本文针对3D打印Ti-6Al-4V合金超声疲劳试样,分析了Micro CT扫描试样得到的三维图像,获得了试样内孔洞缺陷的数量与体积;选择体积分数占比最大的孔洞,采用有限元方法分析了三种不同孔洞分布形式下的局部应力集中现象。研究发现,因空间位置的不同,独立的孔洞、接近自由表面的孔洞、相邻的孔洞三种不同孔洞的分布情况的应力集中系数差异显著。研究结果在一定程度上解释了目前EBM技术打印Ti-6Al-4V合金的孔洞缺陷如何对材料受力后的局部应力情况产生影响。     

软材料粘接结构界面破坏研究综述

软材料已经在软机器人、生物医学及柔性电子等各个领域得到广泛的应用. 实际应用中, 软材料多需要粘附于不同类型的基底上, 与之共同组成工程构件进而实现特定的功能, 粘接界面性能对构件的结构完整性与功能可靠性起着关键性作用. 本文对目前软材料粘接结构界面破坏行为方面的研究进行了系统总结. 首先通过与传统粘接结构的对比, 指出了“软界面”与“软基体”两种软材料粘接结构界面破坏行为的独特性及其物理本质. 接着分别总结了“软界面”与“软基体”两种粘接结构界面破坏行为的实验表征方面的研究成果, 对界面及基体黏弹性耗散对界面破坏机理的影响分别进行了分析. 然后从理论角度, 介绍了针对两种软材料粘接结构界面破坏行为的理论分析方法, 并对已建立的相关理论模型进行了总结. 之后以内聚力模型方法为基础, 介绍了软材料粘接结构界面破坏行为数值模拟方面的相关研究进展. 最后基于已有的研究成果, 提出了目前研究所面临的挑战, 并对可能的软材料粘接结构界面破坏的未来研究方向进行了讨论和展望.      

SHPB试验中粘弹性材料的应力均匀性分析

采用特征线解法，对满足ZWT方程的粘弹性材料在高应变率SHPB试验中的应力均匀性进行了数值研究。着重分析了不同的材料本构粘性(松弛时间)、瞬态波阻抗比和入射波升时对于试样中应力均匀性、应变均匀性和平均应变率等的影响。为今后动态试验的试样设计提供了一定的理论依据。     

Elastic-plastic behavior of elastically homogeneous materials with a random field of inclusions

A multiphase material is considered, which consists of a homogeneous elastic-plastic matrix containing a homogeneous statistically uniform random set of ellipsoidal elastic-plastic inclusions. The elastic properties of the matrix and the inclusions are the same, but the so-called “stress-free strains”, i.e. the strain contributions due to temperature loading, phase transformations, and the plastic strains, fluctuate. A general theory of the yielding for arbitrary loading (by the stress and by the temperature) is employed. The realization of an incremental plasticity scheme is based on averaging over each component of the nonlinear yield criterion. Usually, averaged stresses are used inside each component for this purpose. In distinction from this usual practice physically consistent assumptions about the dependence of these functions on the component's values of the second stress moments are applied. The application of the proposed theory to the prediction of the thermomechanical deformation behavior of a model material is shown.     

Viscoelastic Deformation of a Reinforced Plate with a Crack

The deformation of a viscoelastic reinforced composite is studied. The composite has an axis of elastic symmetry and consists of transversally isotropic fibers and a viscoelastic matrix, which differ by the volume concentration and mechanical characteristics. The material is modeled by a transversely isotropic homogeneous linearly viscoelastic medium with some averaged characteristics. A plate fabricated from the composite in question is weakened by a through mode I crack and is subjected to constant tensile forces. The viscoelastic properties of the matrix material are described by a convolution operator. The Volterra principle is used to derive expressions for the viscoelastic characteristics and crack opening. The irrational function of the integral operator that describes the crack opening is expanded into an operator continued fraction and is represented as the sum of base operators     

Stable growth of penny-shaped crack in viscoelastic composite material under time-dependent loading

Considered is the long-term cracking of the three-dimensional fiber-reinforced viscoelastic composite with a plane penny-shaped crack under time-dependent loading. The composite has a hexagonal structure and consists of elastic isotropic fibers and viscoelastic isotropic matrix. The material is modeled by transversally isotropic homogeneous linearly viscoelastic medium with some averaged characteristics. The crack propagation planecoincides with the plane of isotropy. A ring-shaped yield zone in front of the moving crack is modeled as a Dugdale's zone with time-dependent stresses. Crack growth under deformation of the composite occurs by application of a slowly increasing tensile load; it is normal to the plane of crack propagation. A convolution-type time operator describes the viscoelastic properties of the matrix material. Use is made of the Volterra principle and the theory of long-term cracking of viscoelastic bodies. The irrational function of integral operator associated with the viscoelastic crack opening expression is expanded into a continued fraction of operators. The solution is reduced to the nonlinear integral equations of crack growth. Numerical results are obtained for a specific material. Crack growth kinetics is discussed in connection with the onset of stable crack growth and crack border stress intensity factor.     

Hydraulic impact of “exponential” and nonlinearly viscoplastic media in pipes made of a viscoplastic material

Differential equations are derived and the hydraulic impact process for “exponential” and nonlinearly viscoplastic media in pipes made of a viscoelastic material is analyzed. Hydraulic impact problems for actual media in pipes has been repeatedly treated in the literature [1–4]. The hydraulic impact of a viscous and linearly viscoplastic media in pipes made of an elastic and viscoelastic material was studied in this work. It is well known [5] that many media in the region of low and moderate shear rates reveal a nonlinearity of the flow curve (oil, drilling fluids, polymer solutions and melts, loaded fuels, fuel mixtures, blood, etc.). It should be noted that flexible pipes made of natural materials (pipe boreholes made of polymer materials, membranes of blood vessels, etc.) are described by complicated rheological equations of state for viscoelastic media. Thus a calculation of the influence of nonlinearity of these media and of the viscoelastic properties of the pipe material on the hydraulic impact process is of theoretical and practical interest in many engineering problems.     

Approximate estimates of the high‐temperature creep of structural elements

The creep and longterm strength of structural elements under temperatures much higher than operational temperatures are studied. The actual nonuniform stress–strain state is reduced to a fictitious uniform state by averaging the specific dissipated power over the volume of the body. This allows one to estimate the intensity of creep and longterm strength of structural elements in terms of volumeaveraged energy parameters by using methods of ideal plasticity and considering statically possible stress fields and kinematically possible velocity and strain rate fields.     

考虑损伤的内变量黏弹-黏塑性本构方程

基于Rice 不可逆内变量热力学框架,在约束构型空间中讨论材料的蠕变损伤问题. 通过给定具体的余能密度函数和内变量演化方程推导出考虑损伤的内变量黏弹-黏塑性本构方程. 通过模型相似材料单轴蠕变加卸载试验对一维情况下的本构方程进行参数辨识和模型验证,本构方程能很好地描述黏弹性变形和各蠕变阶段.不同的蠕变阶段具有不同的能量耗散特点. 受应力扰动后,不考虑损伤的材料系统能自发趋于热力学平衡态或稳定态. 在考虑损伤的整个蠕变过程中,材料系统先趋于平衡态再背离平衡态发展. 能量耗散率可作为材料系统热力学状态偏离平衡态的测度;能量耗散率的时间导数可用于表征系统的演化趋势;两者的域内积分值可作为结构长期稳定性的评价指标.      

肌肉类粘弹性超软材料SHPB实验的应力应变均匀性分析

试样的轴向应力和应变均匀问题是霍普金森杆实验行之有效的必要条件,而这对肌肉类粘弹性超软材料尤为重要.本文用量纲分析法,分析出了影响粘弹性材料SHPB实验应力和应变均匀性的无量纲量,结合LS-DYNA数值模拟方法,考察了各无量纲量对粘弹性材料均匀性的影响规律.结果表明,为使肌肉类材料SHPB实验更早达到均匀条件,理想的入...     

On the energy conservation and critical velocities for the propagation of a ＂steady-shock＂ wave in a bar made of cellular material

The propagation of shock waves in a cellular bar is systematically studied in the framework of continuum solids by adopting two idealized material models, viz. the dynamic rigid, perfectly plastic, locking （D-R-PP-L） model and the dynamic rigid, linear hardening plastic, locking （D-R-LHP-L） model, both considering the effects of strain-rate on the material properties. The shock wave speed relevant to these two models is derived. Consider the case of a bar made of one of such material with initial length L 0 and initial velocity v i impinging onto a rigid target. The variations of the stress, strain, particle velocity, specific internal energy across the shock wave and the cease distance of shock wave are all determined analytically. In particular the ＂energy conservation condition＂ and the ＂kinematic existence condition＂ as proposed by Tan et al. （2005） is re-examined, showing that the ＂energy conservation condition＂ and the consequent ＂critical velocity＂, i.e. the shock can only be generated and sustained in R-PP-L bars when the impact velocity is above this critical velocity, is incorrect. Instead, with elastic deformation, strain-hardening and strain-rate sensitivity of the cellular materials being considered, it is appropriate to redefine a first and a second critical impact velocity for the existence and propagation of shock waves in cellular solids. Starting from the basic relations for shock wave propagating in D-R-LHP-L cellular materials, a new method for inversely determining the dynamic stress-strain curve for cellular materials is proposed. By using e.g. a combination of Taylor bar and Hopkinson pressure bar impact experimental technique, the dynamic stress-strain curve of aluminum foam could bedetermined. Finally, it is demonstrated that this new formulation of shock theory in this one-dimensional stress state can be generalized to shocks in a one-dimensional strain state, i.e. for the case of plate impact on cellular materials, by simply making proper replacements of the elastic and plastic constants.     

Scale-dependent constitutive relations and the role of scale on nominal properties

Size effects in strength and fracture energy of heterogeneous materials is considered within a context of scale-dependent constitutive relations. Using tools of wavelet analysis, and considering the failure state of a one-dimensional solid, constitutive relations which include scale as a parameter are derived from a ‘background’ gradient formulation. In the resulting theory, scale is not a fixed quantity independent of deformation, but rather directly dependent on the global deformation field. It is shown that strength or peak nominal stress (maximum point at the engineering stress–strain diagram) decreases with specimen size while toughness or total work to fracture per nominal area (area under the curve in the engineering stress–strain diagram integrated along the length of the considered one-dimensional specimen) increases. This behavior is in agreement with relevant experimental findings on heterogeneous materials where the overall mechanical response is determined by variations in local material properties. The scale-dependent constitutive relations are calibrated from experimental data on concrete specimens.     

软质聚氨酯泡沫的动态压缩力学性能和本构模型

采用Instron 9350落锤试验机研究了中低应变率下软质聚氨酯泡沫的动态压缩力学性能,分析了其应力-应变响应特征和应变率敏感性,讨论了应变率对材料应变率敏感性指数和能量吸收特性的影响,并基于实验结果建立了可准确描述其压缩力学响应的率相关本构模型。结果表明,软质聚氨酯泡沫的静动态压缩应力-应变响应具有典型的三阶段特征,且呈现出明显的应变率强化效应。准静态加载下,材料具有较高的吸能效率但能量吸收值较小,应变率对最大吸能效率和比吸能的影响较小;动态加载下,随着应变率的增加,最大吸能效率显著减小而比吸能明显增大。考虑应变率影响的修正Sherwood-Frost模型和修正Avalle模型都能够很好地表征软质聚氨酯泡沫的静动态压缩应力-应变响应,但修正Avalle模型的参数较少,更便于工程应用。研究结果可为软质聚氨酯泡沫抗冲击结构的设计和优化提供指导。     

Macroscopic modelling of transport phenomena in porous media. 2: Applications to mass,momentum and energy transport

In this second paper, the averaging rules presented in Part 1 are employed in order to develop a general macroscopic balance equation and particular equations for mass, mass of a component, momentum and energy, all of a phase in a porous medium domain. These balance equations involve averaged fluxes. Then macroscopic equations are developed for advective, dispersive and diffusive fluxes, all in terms of averaged state variables of the system. These are combined with the macroscopic balance equations to yield field equations that serve as the core of the mathematical models that describe the transport of extensive quantities in a porous medium domain. It is shown that the methodology of averaging leads to a better understanding of the effective stress concept employed in dealing with transport phenomena in deformable porous media.     

Optimal bounded control of stochastically excited MDOF nonlinear viscoelastic systems

A new procedure for designing optimal bounded control of stochastically excited multi-degree-of-freedom (MDOF) nonlinear viscoelastic systems is proposed based on the stochastic averaging method and the stochastic maximum principle. First, the system is formulated as a quasi-integrable Hamiltonian system with viscoelastic terms and each viscoelastic term is replaced approximately by an elastically restoring force and a visco-damping force based on the randomly periodic behavior of the motion of quasi-integrable Hamiltonian system. Thus, a stochastically excited MDOF nonlinear viscoelastic system is converted to an equivalent quasi-integrable Hamiltonian system without viscoelastic terms. Then, by applying stochastic averaging, the system is further reduced to a partially averaged system of less dimension. The adjoint equation and maximum condition for the optimal control problem of the partially averaged system are derived by using the stochastic maximum principle, and the optimal bounded control force is determined from the maximum condition. Finally, the probability and statistics of the stationary response of optimally controlled system are obtained by solving the Fokker–Plank–Kolmogorov equation (FPK) associated with the fully averaged Itô equation of the controlled system. An example is worked out to illustrate the proposed procedure and its effectiveness.     

Numerical-graphical method for determining characteristics of damaged viscoelastic materials

A method for determining the material functions of nonlinear endochronic theory of aging viscoelastic materials (NETAVEM) with preliminary mechanical damage was developed. The proposed method is based on an analysis of the differences between two graphs of the stress dependence on time obtained in tension with the same constant speed of two specimens made of the same filled polymer material. One of the specimens was not preloaded, and the other was preloaded. The reduced time [1] contained in the NETAVEM constitutive relations and its dependence on the actual time are determined by the distances from the stress axis to two points corresponding to the same stress value and lying on the graphs for the damaged and undamaged specimens. The relaxation kernel is determined in the experiment with the undamaged specimen. These two material functions and the curve obtained for the damaged specimen are used to obtain the NETAVEM aging function, and then the function of viscosity can be calculated. As a result, all characteristics of the damaged material become known, and the strength of structures made of this material can be calculated.