粘弹性双材料界面裂缝的缝端位移场及动态断裂分析

本文研究的是经常在实际工程中遇到的粘弹性双材料界面裂缝的动断裂问题.由于粘弹性自身的复杂性,使得粘弹性双材料界面裂缝缝端应力的奇异性较弹性呈现出更为复杂的形式,从而使动断裂问题的分析变得更为困难.根据此情况,本文采用复阻尼理论反映粘弹性体的运动规律,用复势理论和平面问题复变函数解答的科洛索夫公式推导了粘弹性双材料界面裂缝缝端位移场及动态应力强度因子的求解公式,利用特解边界元进行了粘弹性双域耦合动力响应计算,按求得的公式用位移外推法计算了单边裂纹板在动荷载作用下的动态应力强度因子.分析了粘性,弹模比和缝长对动态应力强度因子的影响,得出了一些有益的结论.     

基于Ogden公式的充压伸缩式水封非线性计算研究

考虑橡胶类止水材料性能对水封设计的重要作用,论文采用Mooney-Rivlin公式和Ogden公式对橡胶类止水材料的单向拉伸、单轴压缩和纯剪实验数据进行曲线拟合与研究.结合工程实例将该两种公式的拟合参数应用于充压伸缩式水封自由外伸量的非线性有限元计算中,并将计算结果与该水封自由外伸量模型实验进行对比研究.结果表明:Og...     

SHPB实验中粘弹性试件内部应力波的传播

建立描述SHPB实验中线性粘弹性试件内部应力波传播的控制方程组,根据试件两端与入射杆及透射杆接触的应力波特征关系给出耦合边界条件.对方程组和定解条件进行Laplace变换,求得试件内部应力在变换域像函数的表达式.采用数值反变换技术进行反Laplace变换,获得试件两端的应力时程曲线.对现有的固定Tal-bot反变换算法进行改进:将入射波像函数分解为基本部分和延迟部分,利用固定Talbot算法对基本部分入射波作用下的波动问题求解,其他部分的解通过延迟定理得到,最终解为两部分的叠加.采用这种改进算法得到的不同入射波下粘弹性试件的内部应力解与传统的基于特征线数值模拟方法的结果吻合.在此基础上探讨了粘弹性试件的几何参数和材料本构参数对透射波波形的影响.     

热载荷作用下功能梯度材料梁的热粘弹性弯曲

由记忆型非均匀热粘弹性材料的积分型本构关系出发,在时空可分离松弛函数假设和平截面几何假设下,通过引进"结构热函数",建立了FGM梁热粘弹性弯曲问题的数学模型及其简化Gurtin型变分原理.在热弹性参数沿厚度方向呈幂律形式变化和热粘弹性松弛函数空域部分沿厚度方向呈指数形式变化的情况下,借助Ritz解和解析解,研究了热载荷作用下材料组分对热弹性/热粘弹性挠度响应和应力分布的影响,发现了热应力反向分布现象.     

各向异性弹性体中矩形孔边缘应力分布

推导了各向异性平面弹性体中矩形孔洞边缘的应力分布公式并给出了数值计算结果。工程中会遇到矩形孔洞的问题，而矩形孔导致的应力集中增加了工程问题的复杂性。目前尚没有简单实用的公式计算矩形孔在各向异性体中应力的分布。本文采用复变函数保角变换的方法推导了矩形孔边缘的应力公式。在利用本文给出的方法，可以计算矩形孔在不同性质的各向异性材料中，承受任意的双轴荷载下孔边缘的应力分布情况。     

测定粘弹性材料动态特性的新方法探讨

本文用近年来兴起的数字谱分析技术对粘弹性材料动态特性的测定问题进行了较为深入的探讨和研究.提出用数字谱分析的方法高效率地进行粘弹性材料动态特性的测试.文中对实际粘弹性材料NYJlA-5-37进行了实验测试.并与以前的结果进行了对比,结果比较令人满意.本文还对试验中存在的静重去除问题以及粘接问题进行了比较深入的分析.     

带式输送机压陷阻力理论研究

为了减小带式输送机的压陷阻力,基于粘弹性材料在动载荷下的能耗理论,覆盖层采用三元件固体模型,对输送机压陷阻力的能耗机理进行了理论研究,推导出输送带通过托辊时的压陷阻力及压陷阻力系数理论公式.并结合具体输送机技术参数、利用公式和Matlab软件对各种影响因素进行了分析和探讨,结果表明公式能够正确反映各种因素对压陷阻力系数的影响.输送带的粘弹性特性是造成输送机压陷阻力的主要因素,可以通过选配合适的输送带橡胶材料获得合适的输送带粘弹性参数,并结合调整输送带和输送机的结构参数达到降低输送机压陷阻力的目的.     

不可压缩和近乎不可压缩粘弹性问题的有限单元法

本文首先从Herrmann泛函出发,导出了不可压缩和近乎不可压缩粘弹性问题的新的本构关系,然后根据虚功原理建立了相应的有限元公式。最后编制了平面问题的有限元分析程序VFAPINP,并计算了实例。计算结果表明所介绍的分析方法和计算程序解决了泊松比趋近或等于0.5时粘弹性应力分析问题,特别适用于固体推进剂的应力分析。     

无限粘弹性平面中椭圆孔口变边界过程的复变函数解答

工程中存在一类几何边界随时间变化的变边界结构,例如土木工程中处于施工阶段的结构。本文以粘弹性岩体中隧道开挖为背景,尝试用变边界问题对应关系和平面弹性复变方法求取无限平面中椭圆孔口自相似变边界情况下的解析解答。首先建立了复变函数法求解变边界粘弹性问题的基本步骤和公式。然后通过建立逆映射函数将已知?平面复位势转至z平面,从而解耦参与拉普拉斯变换的时间与孔口映射函数所带来的时间,从而导出了粘弹性类材料的应力与位移的统一表达。作为一个例子,本文选择Boltzmann粘弹性模型,代入模型参数后得到积分形式的位移、应力解析解,通过与数值解的比较验证了该解答的可靠性,并通过一个算例分析了变边界过程对位移、应力的影响。分析结果显示,采用不同变边界过程的位移、应力变化形态和数值均有差别。本文解答可用于进行地下椭圆孔型隧道在开挖过程中的力学分析,为实际工程提供初步设计的手段。此外,本文给出的方法可用于推导任意形状孔型变边界问题的解答。     

多孔硅橡胶有限变形的粘弹性行为

针对孔隙度较大（孔隙度大于50％）的硅橡胶材料在有限变形时的粘弹性行为，从建立描述材料粘弹性特征的松驰函数和变形特征的应变能函数出发，提出了适合多孔隙、可压硅橡胶材料的非线性粘弹性力学行为的本构关系，松驰函数和应变能函数可解耦为等容和体积变形两部分，并引入了拟时间的概念来反映变形对材料特征时间的影响，利用硅橡胶材料的单轴压缩松驰实验与材料模型进行了对比，讨论了多孔硅橡胶的等容变形和体积变形对应力松驰的影响。     

Waves in nonlinear fluid filled tubes

Wave propagation and shock formation in nonlinear elastic and viscoelastic fluid filled tubes is discussed. For a Mooney-Rivlin material a simple exact solution exhibiting distortionless propagation is found.     

Universal motions for a class of viscoelastic materials of differential type

Universal quasi-static motions for a class of incompressible, viscoelastic materials of differential type are examined. These time dependent motions are similar to corresponding static universal deformations well-known for incompressible, isotropic elastic materials. General details are illustrated for the pure torsion problem, and specific results and physical effects are provided for the viscoelastic Mooney-Rivlin model.     

丙烯酸弹性体的率相关分数阶黏弹性模型研究

丙烯酸弹性体VHB 4910作为一种重要的介电弹性体,在软体机器人、致动器、俘能器和智能隔振器等领域有很好的应用前景.但材料的非线性黏弹性对其力学行为有显著影响.近来分数阶模型在复杂材料的建模中取得了成功.本文基于分数阶有限变形Kelvin-Voigt流变学模型建立弹性体的三维张量本构,并进一步推导单向拉伸情况下的本构...     

Rheological model optimization using advanced evolutionary computation for the analysis of the influence of recycled rubber on rubber blend dynamical behaviour

Rheological models of tested rubber blend samples, with different ratios of smut/recycled material, were developed, analyzed and discussed according to experimental data. The experimental measurements show the dynamical behavior of these samples after a short-term axial load. The rheological model of samples was assumed as a viscoelastic Voigt-Maxwell and hyperelastic Mooney-Rivlin model. Real coded genetic algorithms were implemented for optimal model parameters identification in order to minimize the difference between the real sample behavior and dynamical model output. Numerical results obtained for analyzed mixture were compared with the original sample to verify the presented method. In purpose of predicting the dynamical characteristics for desirable rubber blends, an interrelation was formed based on the concentration of recycled material in sample blends and dynamical model parameters dependence.     

Finite ramp time correction method for non-linear viscoelastic material model

In this paper, a simple formula for the prediction of stress (strain) in the relaxation (creep) period is derived for a non-linear viscoelastic material model which takes into account the finite ramp time. Usually, it is assumed that the ramp time is small and, therefore, loading can be described via Heaviside step function. This assumption, when applied to the material parameters identification process, can lead to a large errors in the values of the approximated material parameters. Especially, for the materials which undergo significant stress decay in the beginning of relaxation the assumption of infinite small ramp time can induce severe errors. With the help of the derived formula more reliable material parameter identification can be accomplished. The proposed method is tested with numerical simulations and compared with analytical results, Heaviside step loading case and method described by Nordin and Varna. Simulations show good agreement with analytical results.     

110甲基乙烯基硅橡胶材料参数的确定

应用INSTRON5544对110甲基乙烯基硅橡胶材料进行了单轴大变形压缩的力学特性测试。本文根据已有的Mooney-Rivilin模型和Yeoh模型的本构关系,采用泰勒级数展开的办法获得了基于小变形的本构关系表达式。利用本文所得本构关系,将实验数据进行了非线性拟合处理,得到了相应的Mooney-Rivilin模型和Yeoh模型的材料参数。拟合所得的名义应力-应变曲线与实测的名义应力-应变曲线吻合较好。利用非线性拟合所得的Mooney-Rivilin模型和Yeoh模型的材料参数,应用ABAQUS软件进行有限元分析计算。通过有限元的分析计算,验证了Mooney-Rivilin模型和Yeoh模型材料参数的正确性与可行性。     

Response of frictional contact problems in thermo-rheologically complex structures

This paper presents a comprehensive computational model for predicting the nonlinear response of frictional viscoelastic contact systems under thermo-mechanical loading and experience geometrical nonlinearity. The nonlinear viscoelastic constitutive model is expressed by an integral form of a creep function, whose elastic and time-dependent properties change with stresses and temperatures. The thermo-viscoelastic behavior of the contacting bodies is assumed to follow a class of thermo-rheologically complex materials. An incremental-recursive formula for solving the nonlinear viscoelastic integral equation is derived. Such formula necessitates data storage only from the previous time step. The contact problem as a variational inequality constrained model is handled using the Lagrange multiplier method for exact satisfaction of the inequality contact constraints. A local nonlinear friction law is adopted to model friction at the contact interface. The material and geometrical nonlinearities are modeled in the framework of the total Lagrangian formulation. The developed model is verified using available benchmarks. The effectiveness and accuracy of the developed computational model is validated by solving two thermo-mechanical contact problems with different natures. Moreover, obtained results show that the mechanical properties and the class of thermo-rheological behavior of the contacting bodies as well as the coefficient of friction have significant effects on the contact response of nonlinear thermo-viscoelastic materials.     

On the radial oscillations of incompressible, isotropic, elastic and limited elastic thick-walled tubes

The finite amplitude, free radial oscillations of a thick-walled circular cylindrical tube are studied for an arbitrary incompressible, isotropic and homogeneous rubber-like material having limiting molecular chain extensibility. First, based on classical results for hyperelastic tubes, some results for thick-walled Mooney-Rivlin tubes are described graphically in the phase plane. Then the periodicity of the finite amplitude, free oscillations of a general limited elastic, thick-walled tube is studied; and some analytical results for the Gent model are illustrated in several numerical examples. Results for thick-walled Gent tubes are compared with those for corresponding Mooney-Rivlin tubes; and the motion of thin-walled Gent tubes is illustrated in the phase plane. Physical conclusions are presented. The period of small amplitude oscillations of an arbitrary elastic or limited elastic tube is derived from relations obtained by a linearization of a general class of equations of which the tube problem is a special case. Classical results of the linear theory are thereby recovered and compared with results for Mooney-Rivlin and Gent tubes.     

The ultratough peeling of elastic tapes from viscoelastic substrates

The peeling of an elastic thin tape from a flat smooth viscoelastic substrate is investigated. Based on a Green function approach and on the translational invariance, a closed form analytical solution is proposed, which takes into account the viscoelastic dissipation in the substrate material.We find that peeling is prevented from taking place, only when the external force is smaller than the one predicted by Kendall's formula for elastic tapes on rigid substrates. However, we also find that, regardless of the value of the applied force, steady state detachment may occur when the elastic tape is sufficiently stiff. In this case, the constant peeling velocity can be modulated by properly defining the geometrical parameters and the material properties of tape and viscoelastic foundation. On the other hand, for relatively high peeling angles or compliant tapes a threshold value of the peeling force is found, above which the steady-state equilibrium is no longer possible and unstable detachment occurs.The present study contributes to shed light on the behavior of pressure sensitive adhesives in contact with viscoelastic substrates like the human skin. At the same time, it can be considered a first step towards a better understanding of the effect of viscoelastic dissipation on the fracture behavior of solids.