无网格方法在平面粘弹性力学问题中的应用

本文综合应用无网格方法(EFGM)、线性粘弹性与弹性力学之间的对应原理,Laplace变换和逆变换等方法求解了拟静态平面弹性和粘弹性力学问题。首先,利用Laplace变换和逆变换推导了平面问题的粘弹性本构关系,建立了拟静态粘弹性平面问题的边值问题;其次,利用粘弹性与弹性力学之间的对应原理得到了Laplace变换域中平面问题的基本方程,在Laplace变换域中建立了相应的泛函,并得到了用无网格方法离散的控制方程;同时,求解了几个拟静态弹性和粘弹性平面问题,给出了它们的表达式和数值结果;最后,采用Laplace逆变换和数值逆变换,得到了粘弹性力学平面问题在物理空间中的解,并比较了由解析解和无网格数值方法所得到的数值结果,可以看到它们是非常吻合的。说明本文方法的正确性和有效性。     

非线性粘弹性平面问题的数学模型及其失记效应

本文以位移为基本未知量,利用非线性粘弹性力学中的Leaderman本构关系和线性几何假设,建立了非线性粘弹性平面问题的数学模型;在粘弹性泊松比为常数的情况下,利用Titchmarsh定理和Laplace变换法证明了非线性粘弹性平面问题与非线性弹性平面问题之间存在着某些对应关系,对应关系为粘弹性问题的求解提供了一种新的思路,利用这些关系可直接从相应弹性问题获得粘弹性平面问题的部分响应,与传统的时域有限差分法相比,计算时间明显缩短,另外,对应关系也揭示了粘弹性结构的失记效应,即结构的部分响应仅与外部输入的现时值有关,而与其历史无关。     

粘弹性力学的对应原理及其数值反演方法

积分变换是处理粘弹性混合边值问题的重要数学工具,积分变换的应用使粘弹性混合边值问题在象空间与相应弹性混合边值问题对应起来,从而使粘弹性混合边值问题的求解可以继承和借鉴弹性问题的求解方法,再利用积分反演方法就可求得时间域粘弹性边值问题的解．本文结合国内外的研究成果,就粘弹性力学中存在的各种对应原理及数值反演方法进行了归类和总结．结合在求解粘弹性边值问题中的应用,对各类方法的特点进行了评述,并指出存在的问题及发展新的数值方法的研究重点．     

考虑损伤的粘弹性梁的纯弯曲

根据粘弹性损伤理论,分析了带损伤粘弹性矩形梁在受纯弯曲时损伤对应力的影响,得到了在Laplace变换域内损伤场和应力场的分布．利用Laplace数值逆变换,分别得到了损伤弹性梁和损伤粘弹性梁的最大应力和最大损伤值,分析了材料的粘性对梁内应力和损伤的影响。     

辛体系下平面热黏弹性圣维南问题的解析解

借助积分变换，将辛体系引入平面热黏弹性问题，建立了基本问题的对偶方程，并将全 部圣维南问题归结为满足共轭辛正交关系的零本征值本征解问题. 同时，利用变量代换和本 征解展开等技术给出了一套求解边界条件问题的具体方法. 算例讨论了几种典型边界条件问 题，描述了热黏弹性材料的蠕变和松弛特征，体现了这种辛方法的有效性.     

一维流体饱和粘弹性多孔介质层的动力响应

本文研究了不可压流体饱和粘弹性多孔介质层的一维动力响应问题。基于粘弹性理论和多孔介质理论，在流相和固相微观不可压、固相骨架服从粘弹性积分型本构关系和小变形的假定下，建立了不可压流体饱和粘弹性多孔介质层一维动力响应的数学模型，利用Laplace变换，求得了原初边值问题在变换空间中的解析解，并利用Laplace逆变换的Crump数值反演方法，得到原动力响应问题的数值解。数值研究了饱和标准线性粘弹性多孔介质层的动力响应，分析了固相位移、渗流速度、孔隙压力及固相有效应力等的响应特征。结果表明，与不可压流体饱和弹性多孔介质相同，不可压流体饱和粘弹性多孔介质中亦只存在一个纵波，并且固相骨架的粘性对动力行为有显著的影响。     

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

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

考虑损伤时带微孔粘弹性体的应力分布

本文根据粘弹性理论和带空沿材料的线性理论，利用Laplace变换及其逆变换，给出了粘弹性损伤材料的一种本构模型。应用这种本构关系。讨论了带损伤的圆环形板或圆筒体在内外压力作用下的平面问题，得到了圆孔边缘附近的应力场和损伤增量场的分布情况，同时根据Laplace变换的终值定理，得到圆环或圆筒的终态应力和损伤增量的分布，并分析了工程中的一个典型的例子。得到了最终损伤的一种近似分布。可为工程设计提供一定的参考。     

矩形梁受分布荷载的辛解答

在弹性力学平面直角坐标辛体系中,采用分离变量法,放弃齐次边界条件,得到了矩形梁侧边受幂函数形式分布荷载问题的辛解答,给出了这类问题在辛体系中的一般解法,分别对矩形梁受法向和切向分布荷载的问题进行了求解,显示了此方法的有效性.辛解法采用对偶的二类变量进行求解,可同时给出位移和应力;由于辛解法能较好地处理各种边界条件,因此不仅能求解静定问题,也能直接求解静不定问题.     

哈密顿体系与弹性楔体问题

将哈密体系引入到级坐标下的弹性力学楔体问题,利用该体系辛空间的性质,将问题化为本征值和本征向量求解上,得到了完备的解空间,从而改变了弹性力学传统的拉格朗日体系以应力函数为特征的半逆法的讨论去解决该类问题的思路,给出了一条求解该类问题的直接法。     

Stresses in rotating heterogeneous viscoelastic composite cylinders with variable thickness

An analytical solution is presented for the rotation problem of a two-layer composite elastic cylinder under a plane strain assumption. The external cylinder has variable-thickness formulation, and is made of a heterogeneous orthotropic material. It contains a fiber-reinforced viscoelastic homogeneous isotropic solid core of uniform thickness. The thickness and elastic properties of the external cylinder are taken as power functions of the radial direction. By the boundary and continuity conditions, the radial displacement and stresses for the rotating composite cylinder are determined. The effective moduli and Illyushin’s approximation methods are used to obtain the viscoelastic solution to the problem. The effects of heterogeneity, thickness variation, constitutive, time parameters on the radial displacement, and stresses are investigated.     

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

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

A combined approach of the Laplace transform and Padé approximation solving viscoelasticity problems

In this paper the Laplace transform method is combined with Padé approximations to solve linear viscoelastic problems. This approach allows to avoid the usual difficulties of original function determination. An algorithm is given to find solution with arbitrary precision. As an example the solution for problem of viscoelastic orthotropic half-plane stress state under concentrated normal force is given.     

Inverse problem of deformation of a physically nonlinear inhomogeneous medium

An isotropic linearelastic (viscoelastic) plane containing various physically nonlinear elliptic inclusions is considered. It is assumed that the distances between the centers of the inclusions are much greater than their dimensions. The problem is to determine the orientation of the inclusions and the loads applied at infinity which ensure a specified value of the principal shear stress in each inclusion. Necessary and sufficient conditions of existence of the solution of the problem are formulated for a plane strain of an incompressible inhomogeneous medium.     

无限粘弹性平面中孔洞扩展的时变力学解析解

从粘弹性时变力学基本方程出发,针对Maxwell本构模型,导出了双向等压下粘弹性平面中圆孔半径任意规律扩展时的时变力学解析解.在内径线性时变时,其结果与用对应原理法得出的解相同,其解的可靠性得到验证.与对应原理法相比,文中的结果适用于半径任意时变情况,更具一般性.     

Exact Solution of the Problem on a Six‐Constant Jeffreys Model of Fluid in a Plane Channel

An exact analytic solution of the problem of a generalized viscoelastic Jeffreys fluid flow in a plane channel under the action of a pressure gradient is found. The velocity profiles are obtained in a parametric form with a velocity gradient taken as a parameter. The critical values of the pressure gradient are determined, which, when exceeded, lead to weak tangential discontinuities in the longitudinal velocity profile. When the pressure gradient changes smoothly over some range of parameters, a hysteresis loop emerges on the graph of the flow rate versus the pressure gradient.     

An Edge Crack Problem in a Semi-Infinite Plane Subjected to Concentrated Forces

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｅｄｇｅｃｒａｃｋｐｒｏｂｌｅｍｉｎａｓｅｍｉ＿ｉｎｆｉｎｉｔｅｐｌａｎｅｗａｓｃｏｎｓｉｄｅｒｅｄｂｙｍａｎｙｉｎｖｅｓｔｉｇａｔｏｒｓ.ＵｓｉｎｇｔｈｅＦｒｅｄｈｏｌｍｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎａｎｄｔｈｅａｌｔｅｒｎａｔｉｎｇｍｅｔｈｏｄ ,ａｎｅｄｇｅｃｒａｃｋｐｒｏｂｌｅｍｉｎａｓｅｍｉ＿ｉｎｆｉｎｉｔｅｐｌａｎｅｗａｓｓｏｌｖｅｄ[1].Ｔｈｅｏｂｔａｉｎｅｄｒｅｓｕｌｔｓａｒｅｌｉｍｉｔｅｄｔｏｔｈｅｃａｓｅｔｈａｔｔｈｅｅｄｇｅｃｒａｃｋｉｓｎｏｒｍ…     

Identification of viscoelastic properties by means of nanoindentation taking the real tip geometry into account

Motivated by recent progress in viscoelastic indentation analysis, the identification of viscoelastic properties from nanoindentation test data taking the real tip geometry into account is presented in this paper. Based on the elastic solution of the indentation problem, the corresponding viscoelastic solution is obtained by the application of the method of functional equations. This general solution, which accounts for the real geometric properties of the indenter tip, is specialized for the case of a trapezoidal load history, commonly employed in nanoindentation testing. Three deviatoric creep models, the single dash-pot, the Maxwell, and the three-parameter model are considered. The so-obtained expressions allow us to determine viscoelastic model parameters via back calculation from the measured load–penetration history. The presented approach is illustrated by the identification of short-term viscoelastic properties of bitumen. Hereby, the influence of loading rate, maximum load, and temperature on the model parameters is investigated.