非线性粘弹性拟静态问题与非线性弹性静力问题对应原理

本文应用多重单边拉氏变换导出了非线性粘弹性拟静态问题与非线性弹性静力问题的对应关系,从而把过去认为只有线性条件下存在的粘弹性——弹性对应原理拓展到了非线性范畴     

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

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

一种适合橡胶类材料的非线性粘弹性本构模型

借助非线性流变模型建立大变形情况非线性粘弹材料的本构关系,考虑到大多数橡胶类材料具有的几乎不可压缩性,以及体积响应和剪切响应的流变性能不同,将变形梯度乘法分解为等容部分和体积变形部分,给出了一种适合橡胶类材料的非线性粘弹性本构模型,并模拟了粘滞效应。对于极快或极慢的过程,该模型退化为橡胶弹性理论;在小变形情况下退化为经典的广义Maxwell粘弹性材料。模型与热力学第二定律相容,适合于大规模数值分析。     

平面粘弹性问题的辛求解方法

将弹性力学辛对偶求解方法与Laplace变换相结合,提出了一个求解粘弹性平面问题的新方法。首先利用Laplace变换,将粘弹性平面问题转化为一个准弹性问题,在辛弹性力学的框架下,利用分离变量和辛本征展开法对其进行求解,然后由逆变换得到原问题的解。为证明方法的有效性,求解分析了矩形域平面粘弹性圣维南问题,得到了令人满意的结果。     

幂率型非线性粘弹性裂纹尖端场

研究幂率型非线性粘弹性裂纹尖端场．为了推导的需要,首先列出了幂率型硬化材料的HRR奇异场和高阶渐近场．论证了它们实质上是各向同性、不可压缩、幂率型、非线性弹性裂纹尖端场．回颐了求解非线性粘弹性问题的弹性回复对应原理之后,给出了在第一类边界条件下求解幂率型非线性粘弹性材料裂纹问题的对应原理．得到了幂率型非线性粘弹性材料,特别是改性聚丙稀的裂纹尖端应力、应变和位移场的解答．     

粘弹性结合材料界面端奇异场

在电子封装等结构中存在大量的粘弹性界面问题,其破坏一般均始于界面端,但目前尚无关于粘弹性界面端奇异场的解.粘弹性问题在拉普拉斯域内与弹性问题有对应关系,理论上可以利用对应性原理由弹性解经拉氏逆变换得到粘弹性问题的解.但是,对于粘弹性界面端,由于奇异场的奇异指数也是与时间有关的,因此进行严密的拉氏逆变换是非常困难的.本文借鉴弹性界面端奇异场,近似地给出了线性粘弹性体界面端奇异场的具体形式,并通过数值计算验证了近似理论解的有效性.     

粘弹性界面裂纹奇异场

对于许多粘弹性问题,通常可以利用对应性原理,即由弹性问题的结果得到对应的粘弹性问题在拉普拉斯变换域内的解,再通过反演变换求得最终时域中的解.但是,由于界面裂纹场存在着振荡奇异性,弹性问题解的形式就已经非常复杂,对应的粘弹性问题要通过反演变换直接求得准确的解析解几乎是不可能的.本文在利用对应性原理时做了更简单的准静态处理,即将弹性结果中的材料参数用粘弹性材料参数做对应替代,得到了粘弹性界面裂纹场近似的经典解,并与有限元分析结果作了比较.同时,利用Comninou接触模型,对粘弹性界面裂纹在远场拉剪混合加载情况下的裂尖应力场和接触区做了考察,并与经典解作了比较.     

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

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

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

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

几何缺陷浅拱的动力稳定性分析

研究了几何缺陷对粘弹性铰支浅拱动力稳定性能的影响。从达朗贝尔原理和欧拉-贝努利假定出发推导了粘弹性铰支浅拱在正弦分布突加荷载作用下的动力学控制方程,并采用Galerkin截断法得到了可用龙格-库塔法求解的无量纲化非线性微分方程组。同时引入能有效追踪结构动力后屈曲路径的广义位移控制法,对含几何缺陷浅拱的响应曲线进行几何、材料双重非线性有限元分析。用这两种方法分析了前三阶谐波缺陷对浅拱动力稳定性能的影响,其中动力临界荷载由B-R准则判定。主要结论有:材料粘弹性使浅拱动力临界荷载增大且结构响应曲线与弹性情况差别很大;二阶谐波缺陷影响显著,它使动力临界荷载明显下降且使得浅拱粘弹性动力临界荷载可能低于弹性动力临界荷载。     

含随机参数的非线性黏弹性有限元计算

进行了粗粒土与结构接触面单调和循环加载试验，基于宏细观测量结果, 扩展了 损伤概念以 描述该类接触面在受载过程中的物态演化, 及由于物态演化导致的力学特性从初始状态到最终 稳定状态的连续变化过程. 揭示了接触面损伤的细观物理基础主要是接触面内土的颗粒破碎 和剪切压密这两种物态演化；指出接触面的剪胀体应变可以划分为可逆性和不可逆性剪胀体 应变两部分，其中不可逆性剪胀体应变可作为接触面损伤发展的宏观量度，因此其归一化 形式可作为一种损伤因子的定义；提出了建立粗粒土与结构接触面一种损伤本构关系的基本思路.     

Large deformations of nonlinear viscoelastic and multi-responsive beams

This study presents analyses of deformations in nonlinear viscoelastic beams that experience large displacements and rotations due to mechanical, thermal, and electrical stimuli. The studied beams are relatively thin so that the effect of the transverse shear deformation is neglected, and the stretch along the transverse axis of the beams is also ignored. It is assumed that the plane that is perpendicular to the longitudinal axis of the undeformed beam remains plane during the deformations. The nonlinear kinematics of the finite strain beam theory presented by Reissner [27] is adopted, and a nonlinear viscoelastic constitutive relation based on a quasi-linear viscoelastic (QLV) model is considered for the beams. Deformation in beams due to mechanical, thermal, and electric field inputs are incorporated through the use of time integral functions, by separating the time-dependent function and nonlinear measures of field variables. The nonlinear measures are formulated by including higher order terms of the field variables, i.e. strain, temperature, and electric field. Responses of beams under mechanical, thermal, and electrical stimuli are illustrated and the effects of nonlinear constitutive relations on the overall deformations of the beams are highlighted.     

Numerical analysis of the stress concentration near holes originating in previously loaded viscoelastic bodies at finite strains

The state of stress and strain of previously loaded viscoelastic bodies with holes originating in them, successively or simultaneously, is analyzed under finite plane deformations. The problem statement and solution are based on the theory of repeatedly superimposed large deformations. The material mechanical properties are described using integral relations of the convolution type over time with a weakly singular kernel. The problem solving is based on the finite-element method. To calculate the integral of the convolution type, a recurrence formula is used that can be obtained by approximating the initial kernel with a linear combination of exponential functions (the truncated Prony’s series). The nonlinear effects and the effect of the interaction between holes on the stress concentration are analyzed. For the dynamic problems, the results for incompressible and weakly compressible materials are compared.     

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.     

粘弹性大变形动力响应的有限元分析

本文基于ＴｏｔａｌＬａｇｒａｎｇｉａｎ增量叠加方法，采用Ｋｉｒｃｈｈｏｆｆ应力增量和Ｇｒｅｅｎ应变增量表示的动力虚功方程和Ｋｉｒｃｈｈｏｆｆ应力－Ｇｒｅｅｎ应变的单积分型本构关系，导出粘弹性大变形的动力变分方程。依此采用Ｎｅｗｍａｒｋ法和八节点轴对称等参数元与二十节点三维等参数元编制了轴对称及三维问题的动力响应计算程序，典型例题的计算结果表明分析符合结构的物理性质。     

SOME IDENTITY RELATIONS BETWEEN PLANE PROBLEMS FOR VISCO-AND ELASTICITY

I.IntroductionThereareimportantapplicationsfortheor}:ofplanea'iscoelasticit}'Inthefieldsofgeology,miningandconstructingetc.,butformostproblemsofviscoelastici[}'.theirsolutionsareobtainedfromthecorrespondingelasticsolutionsb}'"leansofthecorrespondenceprinc…     

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

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

Development of the large increment method for elastic perfectly plastic analysis of plane frame structures under monotonic loading

The displacement-based finite element method dominates current practice for material nonlinear analysis of structures. However, there are several characteristics that may limit the effectiveness of this approach. In particular, for elastoplastic analysis, the displacement method relies upon a step-by-step incremental approach stemming from flow theory and also requires significant mesh refinement to resolve behavior in plastic zones. This leads to computational inefficiencies that, in turn, encourage the reconsideration of force-based approaches for elastoplastic problems.One of these force algorithms that has been recently developed is the large increment method. The main advantage of the flexibility-based large increment method (LIM) over the displacement method is that it separates the global equilibrium and compatibility equations from the local constitutive relations. Consequently, LIM can reach the solution in one large increment or in a few large steps, thus, avoiding the development of cumulative errors. This paper discusses the extension of the large increment methodology for the nonlinear analysis of plane frame structures controlled by an elastic, perfectly plastic material model. The discussion focuses on the power of LIM to handle these nonlinear problems, especially when plastic hinges form in the frame and ultimately as the structure approaches the collapse stage. Illustrative planar frame examples are presented and the results are compared with those obtained from a standard displacement method.     

Using Fourier differential quadrature method to analyze transverse nonlinear vibrations of an axially accelerating viscoelastic beam

In this paper, a Fourier expansion-based differential quadrature (FDQ) method is developed to analyze numerically the transverse nonlinear vibrations of an axially accelerating viscoelastic beam. The partial differential nonlinear governing equation is discretized in space region and in time domain using FDQ and Runge–Kutta–Fehlberg methods, respectively. The accuracy of the proposed method is represented by two numerical examples. The nonlinear dynamical behaviors, such as the bifurcations and chaotic motions of the axially accelerating viscoelastic beam, are investigated using the bifurcation diagrams, Lyapunov exponents, Poincare maps, and three-dimensional phase portraits. The bifurcation diagrams for the in-plane responses to the mean axial velocity, the amplitude of velocity fluctuation, and the frequency of velocity fluctuation are, respectively, presented when other parameters are fixed. The Lyapunov exponents are calculated to further identify the existence of the periodic and chaotic motions in the transverse nonlinear vibrations of the axially accelerating viscoelastic beam. The conclusion is drawn from numerical simulation results that the FDQ method is a simple and efficient method for the analysis of the nonlinear dynamics of the axially accelerating viscoelastic beam.     

黏弹性传动带1:3内共振时的周期和混沌运动

研究了参数激励作用下黏弹性传动带在1:3内共振时的周期解分岔和混沌动力学. 同时考虑传动带的线性外阻尼因素和材料内阻尼因素. 首先建立了具有线性外阻尼情况下的黏弹性传动带平面运动时的非线性动力学方程, 黏弹性材料的本构关系用Kelvin模型描述. 然后考虑黏弹性传动带的横向振动问题, 利用多尺度法和Galerkin离散法得到黏弹性传动带系统在1:3内共振时的平均方程. 最后利用数值模拟方法研究了黏弹性传动带系统的周期振动和混沌动力学, 得到了系统在不同参数下的混沌运动. 数值模拟结果说明黏弹性传动带系统存在周期分岔, 概周期运动及混沌运动.