粘弹性随机有限元

以近似不可压缩粘弹增量有限元和摄动法为基础，利用增量法处理遗传积分，应用参数摄动考虑随机性，采用局部平均方法对随机场进行离散，通过相关结构分解减少计算量，发展了一种粘弹性随机有限元方法。研究表明，尽管粘弹性材料本构关系具有时间相依性，其随机摄动格式并不存在“长期项”的影响。应用该方法进行粘弹性结构的随机模拟，程序实施简单，计算效率较高，精度较高。     

线性随机结构的非平稳随机响应变异性分析

对于具有随机参数的结构受到非平稳随机激励的问题，给出了结构随机响应变异系数的虚拟激励摄动算法。它应用虚拟激励法先将随机荷载转化为确定性荷载，以使随机问题精确地转化为仅结构参数具有随机性的问题，从而将问题归结为应用随机摄动法求解单随机问题。求解过程简单高效，且有较高的精度。     

粘滞和粘弹性阻尼器减震结构的随机响应特性

利用积分型本构关系,针对带支撑任意线性粘滞和粘弹性阻尼器单自由度耗能结构,建立了微分和积分混合地震响应方程;基于随机平均分析法,推导出耗能结构振幅与相位瞬态联合概率密度函数、位移与速度瞬态联合概率密度函数、位移与速度瞬态响应方差、振幅动力可靠性、振幅首超时间任意阶统计矩的一般解析解;给出了带支撑广义Maxwell阻尼器和广义微分模型阻尼器耗能结构上述各种随机响应特性,从而建立了带支撑任意线性粘滞和粘弹性阻尼器单自由度减震结构的各种随机响应特性分析的统一解析解法。     

树脂基复合材料板的粘弹性损伤本构关系

1 引言一般树脂基复合材料板具有相当强烈的各向异性和非均匀性,受载后很容易发生基体裂纹群等损伤.同时,即使在常温下,这类材料也显示粘弹性.而且,损伤与粘弹性都是各向异性的.因此,复合材料力学响应的分析比均匀无损的粘弹性材料要复杂得多,困难得多.实际上,在复合材料的结构强度和尺寸稳定性设计中,它的时间相关性和存在损伤是两个不能回避的重要问题.为此,特别需要复合材料粘弹性损伤本构关系的知识.最近一个时期,复合材料的粘弹性本构关系已得到一定的研究.作者曾提出适用于复合材料分析的弹脆性损伤模型以及考虑损伤的粘弹性本构关系.在此基础     

分数导数型本构关系描述粘弹性梁的振动分析

本文研究粘弹性梁在周期激励作用下的受迫振动问题。梁的材料满足Kelvin-Volgt分数导数型本构关系。基于动力学方程、本构关系和应变－位移关系建立了小变形粘弹性梁的振动方程。采用分离变量法分析粘弹性梁的自由振动,导出模态坐标满足的常微分－积分方程和模态函数满足的常微分方程,对于两端简支的截面梁给出了固有频率和模态函数。对于简谐激励作用下粘弹性梁的受迫振动,利用模态叠加得到了稳态响应。最后给出数值算例说明本文方法的应用。     

固体推进剂非线性拉伸行为的实验研究

固体推进剂非线性拉伸行为的实验研究王晓明，沈亚鹏（西安交通大学，７１００４９）关键词本构关系，粘弹性，非线性，推进剂１引言Ｈａｄｄａｄ（１）和杨挺青（２）对粘弹性非线性本构方程作了比较全面的综述，其中最简单，而又有利于实验确定材料常数的本构方程是其中...     

设置连接阻尼器的相邻结构随机地震动响应的简明封闭解

针对两相邻结构间设置连接阻尼器对结构的减震影响问题,研究了基于Kanai-Tajimi谱地震动激励下的Kelvin型粘弹性阻尼器与相邻结构形成的组合体系的随机地震动系列响应(绝对位移及层间位移)的简明封闭解。首先,利用Kelvin型粘弹性阻尼器本构关系及Kanai-Tajimi谱的滤波方程,将组合体系基于复杂地震动激励精确转化为基于简明白噪声激励的运动方程;其次,利用复模态法获得了组合结构相对于地面的绝对位移、层间位移等系列响应方差及0阶~2阶谱矩的统一简明封闭解。最后,通过算例及与虚拟激励法进行对比,证明本文方法的正确性和简明性;通过与未设置阻尼装置结构体系的动力响应对比,说明了阻尼装置对相邻结构具有良好的减震效果,但局部楼层的层间位移及层间剪力会有所增加。     

刚性微粒填充的高聚物非线性粘弹性本构关系的微力学分析

在应变具有可加性条件下,本文分析表明Eshelby线性微力学理论,包括其等效包含体方法和相应的平均化方法可推广应用于研究高刚度微粒填充高聚物的非线性粘弹性本构关系。     

分数阶导数粘弹性模型的有限元法

在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数值解。并以二维沥青路面结构为例进行了路面动态粘弹性响应分析。算例分析表明,该方法能够正确有效地进行路面动态粘弹性分析。     

粘弹性有限元法研究

以粘弹性材料的积分蠕变型本构关系为基础，建立了粘弹结构有限元分析模型，大量数值计算结果表明，本模型是行之有效，精度很高的。     

A non‐iterative implicit algorithm for the solution of advection–diffusion equation on a sphere

A numerical algorithm for the solution of advection–diffusion equation on the surface of a sphere is suggested. The velocity field on a sphere is assumed to be known and non‐divergent. The discretization of advection–diffusion equation in space is carried out with the help of the finite volume method, and the Gauss theorem is applied to each grid cell. For the discretization in time, the symmetrized double‐cycle componentwise splitting method and the Crank–Nicolson scheme are used. The numerical scheme is of second order approximation in space and time, correctly describes the balance of mass of substance in the forced and dissipative discrete system and is unconditionally stable. In the absence of external forcing and dissipation, the total mass and L₂‐norm of solution of discrete system is conserved in time. The one‐dimensional periodic problems arising at splitting in the longitudinal direction are solved with Sherman–Morrison's formula and Thomas's algorithm. The one‐dimensional problems arising at splitting in the latitudinal direction are solved by the bordering method that requires a prior determination of the solution at the poles. The resulting linear systems have tridiagonal matrices and are solved by Thomas's algorithm. The suggested method is direct (without iterations) and rapid in realization. It can also be applied to linear and nonlinear diffusion problems, some elliptic problems and adjoint advection–diffusion problems on a sphere. Copyright © 2015 John Wiley & Sons, Ltd.     

Solving the Stress Problem for a Finite Cylinder by the Schwarz Method with Hybrid Approximations

An axisymmetric elastic problem is solved by a heterogeneous numerical scheme of the alternating Schwarz method based on the finite- and boundary-element methods. The Schwarz method and the modified scheme are compared. A relaxation parameter is used to accelerate the iterative process     

自由单元法及其在结构分析中的应用

通过吸收有限元与无网格法的优点,提出了一种新的数值方法------自由单元法.此方法在离散方面,采用有限元法中的等参单元,表征几何形状和进行物理量的插值;在算法方面,采用单元配点技术,逐点产生系统方程.主要特点是,在每个配置点只需要一个和周围自由选择的节点而形成的一个独立的等参单元,因而不需要考虑物理量在单元之间的相互连接关系与导数连续性问题. 本文介绍强形式与弱形式两种自由单元法,前者直接由控制方程和边界条件直接产生系统方程,后者通过在自由单元上建立控制方程的加权余量式产生弱形式积分式,并通过像传统有限元法中的积分过程建立系统方程组.本文提出的方法是一种单元配点法,对于域内点为了获得较高的导数精度,需要采用至少具有一个内部点的等参单元,为此除了可使用各阶次的拉格朗日四边形单元外, 还 给出了七节点三角形等参单元,用于模拟较为复杂的几何形状问题.      

A new modification of false position method based on homotopy analysis method

A new modification of false position method for solving nonlinear equations is presented by applying homotopy analysis method （HAM）. Some numerical illustrations are given to show the efficiency of algorithm.     

Preconditioned Iterative Methods for Linear Equations Arising from the Boundary Element Method

A precondition for the Gauss–Seidel iterative method to solve a linear system of equations arising from the boundary element method for the Laplace and convective diffusion with first-order reaction problems is presented in this paper. The present precondition is based on the elementary matrix operation. We discuss the effect of the precondition in comparison with the Gauss elimination (GE) method in some numerical experiments.     

无网格MPM法三维前处理系统设计

无网格法在解决爆炸冲击、裂纹扩展、高速碰撞以及大变形等问题时,由于与网格无关,能够避免有限元法因网格畸变计算精度降低或计算终止的矛盾,是目前求解固体力学问题的有效方法之一,也是近年来研究的热点。前处理对于无网格法数值计算是非常重要的环节,基于面向对象程序设计思想,应用C 和OpenGL开发了三维图形交互方式的MPM法前处理系统,提出了应用递归四叉树法、扫描表示法和CSG法划分任意形状非等间距的MPM法质点单元,为建立通用的无网格前处理平台提出了有效的解决方案。     

参数多项式方法的扩展——参数分布形状特征式法

在深入分析参数多项式方法数学原理的基础上，提出了参数分布形状特征式法。根据流动问题的物理本质及大量试验数据，总结出参数分布的形状特征式，将其代入流动基本方程，就可以快速而准确地求解流场。它是参数多项式方法的一种扩展，具有重要的应用前景。     

基于核重构思想的配点型无网格方法的研究--一维问题

无网格方法按其离散原理可分为Galerkin型、配点型等。其中Galerkin型无网格方法的实施需要背景网格，不属于真正的无网格法；配点型无网格方法的实施不需要背景网格，是真正的无网格法。本文首先介绍了重构核点法的基本原理，然后基于核重构思想，与配点法相结合，以一维问题为例，研究了配点型无网格方法，对该方法构造过程中的近似函数及其导数的计算、修正函数的计算及方法的实现等问题进行了探讨。并结合若干典型算例，检验了其计算精度与收敛姓。     

无单元Galerkin法和边界元耦合法

无网格方法与有限元法或边界元法耦合是无网格方法处理边界条件的方法之一,在无网格方法中研究无网格方法与有限元法或边界元法耦合的研究显得非常重要.本文在无单元Galerkin法和边界元法的基础上,基于无单元Galerkin法子域和边界元法子域的界面上位移连续和面力平衡条件,提出了一种新的无单元Galerkin法和边界元法的直接耦合方法,对弹性力学问题详细推导了在整个求解域上的耦合公式.与以往的耦合法相比,这种方法简单直观,不需要增加新的耦合区域,也不需要建立新的逼近函数来保证界面位移的连续性.算例结果表明,该方法具有较好的计算精度.     

离散元数值模拟中查找邻居元关系的改进算法

离散元方法是近年发展起来的有前途的新概念数值模拟方法,然而查找邻居元关系耗费大量机时和内存。本文引入一种改进算法-“区域法”,并对其原理、算法、运行时间和空间进行了分析。和原算法中的“窗口法”的比较结果表明,新算法在运行时间和空间方面均明显优于原算法。