基于边界元法的弹性动力问题并行求解

为了扩大结构弹性动力分析的规模和提高分析速度,在微机机群环境下给出了两种基于边界元法的瞬态问题并行求解算法,即并行拉普拉斯变换求解算法和并行时域求解算法.并行拉氏变换法通过拉氏变换隐去时间变量,由各结点机独立求解各自负责的变换边界元问题.并行时域法采用与时间有关的基本解,使得边界元系统矩阵可以实现时间域上的并行形成.系数矩阵采用卷帘存储,以保持负载平衡.通过矩阵向量运算的并行化实现时间步进算法的并行化.理论分析和数值试验结果表明:两种算法都具有较好的并行性能.可以用于大型问题的高效求解.     

物性值随温度变化热弹性问题的摄动—边界元分析

本文将摄动法和边界元法相结合求解物性值随温度变化的热弹性问题,简述了基本方程和积分方程的建立,导出了有关计算公式。算例表明本文方法简便、有效。     

一种处理弹性动力边界元法中奇异积分的方法

利用边界元法求解瞬态弹性动力学问题时,时域基本解函数的分段连续性和奇异性为该问题的求解带来很大的困难。为了解决时域基本解中的奇异性问题,本文依据柯西主值的定义,对经过时间解析积分之后的时域基本解进行奇异值分解,将其分成奇异和正则积分两部分;其中正则部分可通过采用常规高斯积分方法来计算,而奇异部分具有简单的形式,可以利用解析积分计算。经过上述操作之后,就可以达到直接消除时域基本解中奇异积分的目的。和传统方法相比,本文方法并不依赖静力学基本解来消除奇异性,是一种直接求解方法。最后给定两个数值算例来验证本文提出方法的正确性和可行性,结果表明使用本文算法可以解决弹性动力学边界积分方程中的奇异性问题。     

弹性力学中一种新的边界轮廓法

利用基本解的特性,将面力积分方程化成仅含有Ｃａｕｃｈｙ主值积分的形式,基于这种边界积分方程,提出了一种新的边界轮廓法,对于三维问题,该方法只须计算沿边界单元界线的线积分,对二维问题,则只需计算边界单元两点的热函数之差,无须进行数值积分计算,实例计算说明该方法是有效的。     

一种求解薄板稳定及振动问题的边界元法

本文利用薄板小挠度弯曲问题的基本解,引入等效荷载的概念,用以求解薄板稳定问题及振动问题,取得良好的结果。本文所提出的方法也可用于其他的薄板问题。     

边界元法求解瞬态弹性动力学问题的一种新解法

本文首次将Newmark法引入到边界元法之中,从而使得边界元法与有限元法有机地结合起来,使边界元法的通用性在求解瞬态弹性动力学问题上大大加强,在程序的实现过程中应用一定的技巧,提出了多重子单元划分法,解决了时间步长与单元网格之间的耦合关系。提高了计算精度,缩短了计算时间,从所给出的算例可知,此算法是可行的,程序也是可靠的。     

变边界粘弹性轴对称问题的复变函数法

对边界几何形状、位置随时间变化的变边界结构,给出了用复变函数求解粘弹问题的解析方法.文中用拉普拉斯变换结合平面弹性复变方法,对内外边界变化时粘弹性轴对称问题进行求解.引入两个与时间、空间相关的解析函数,给出了变边界情况下应力、位移以及边界条件与解析函数的关系.当解析函数形式部分确定,则可用边界条件求解其中与时间相关的待定函数.求解待定函数的方程一般情况下为一系列积分方程,特殊情况可求得解析解.对轴对称问题中应力边值问题、位移边值问题以及混合边值问题,分别利用边界条件求得相关系数,从而得到了应力与位移的解析表达.当取Boltzmann粘弹模型时,进行不同边值问题的分析.分析显示,应力、位移的形态与大小均与边界变化过程相关,与固定边界粘弹性问题有较大不同.本文解答可用于粘弹性轴对称问题内外边界任意变化及各种边值问题的力学分析.此外,该法可进一步进行荷载非对称、复杂孔型变边界问题的求解.     

三维弹性快速多极边界元法

将静电场多极展开法和广义极小残值法结合于三维弹性问题的边界元法,使其求解的计算量及所需内存量同节点的自由度总数成正比,变革计算结构,加快求解速度以适应大规模数值计算。两者结合的关键点在于边界元法基本解的合理分解,并用广义极小残值法(GMRES)求解方程。轧机支承辊变形场大规模数值算例的总自由度数首次达N=34008并获得成功。清晰地描述了支承辊和工作辊接触区的辊型。     

断裂力学的相似边界元法及其应用

首先对弹性力学的相似边界元法进行了研究,推导了相应的计算公式。与传统的边界元法相比,相似边界元法由于只需在少数单元上进行数值积分,当边界单元数目较多时大大减少了计算量。在此基础上,将相似边界元法应用于断裂力学,对路面断裂力学问题进行了计算,与有限元法的结果比较,说明了本文方法在减少计算量的情况下仍能较好地保证精度。     

边界元法边界层处理的一种新方法—拟多连域法

提出了处理边界元法边界层问题的新方法。这种方法在求得边界未知量后，在距边界较远的域中构造一工作边界，将实际边界与工作边界看作一多连域问题求得工作边界的未知量。再将工作边界与边界层边界视为一多连域问题求解，可得到满足精度要求的边界附近点的位移与应力。这种方法理论简洁、计算方便、有效、精度高，对于需求边界上多点的未知量问题很具优越性     

结动振动摄动分析的新方法

提出了一种用于结构动力修改的设计灵敏度分析的新方法。它发展了Ｎｅｌｓｏｎ［１］陈［２］等人的方法，较好地解决了结构修改量大时计算精度低之间的矛盾。数值算例表明，本文新方法计算精度高，易于计算机实施。     

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

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

Laplace变换法及其在粘弹性波中传播的应用

针对波传播分析理论的发展历程进行了简要的综述,详细介绍了几种处理粘弹性波传播问题的分析方法,重点讲解Laplace变换法以及Laplace变换在粘弹性波中的应用,对比分析几种方法在各自应用上的优劣,由于Laplace变换法能准确地描述应力波在任意时刻、任意点的波动情况,在处理大尺寸混凝土类构件中应力波传播问题时具有其独特的优势。     

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

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

?????????????????????????е??й????

 系统综述了计算机透析成像技术(CT)的基本原理及其在岩土工程中的应用现状,包 括CT的图像重建理论、CT图像分析的进展、岩土材料损伤的探测方法、岩土材料结构变 化的动态监测、岩土材料本构模型的建立等方面.     

基于等效黏性弹簧的黏弹性Timoshenko裂纹梁弯曲解析解

考虑裂纹的缝隙和黏性效应,将梁中横向裂纹等效为黏弹性扭转弹簧,利用广义Delta函数,给出了Laplace变换域内裂纹梁的等效抗弯刚度,得到了具有任意开闭裂纹数目且满足标准线性固体黏弹性本构的Timoshenko梁在时间域内的弯曲变形显式解析通解．在此基础上,通过两个数值算例,分析了时间、梁跨高比和裂纹深度等参数对黏弹性Timoshenko开裂纹梁弯曲变形的影响．结果表明：裂纹黏性对Timoshenko裂纹梁的弯曲具有显著的影响．相比于裂纹的弹性扭转弹簧模型,考虑裂纹黏性效应的黏弹性Timoshenko裂纹梁在裂纹处挠度尖点和转角跳跃现象十分明显．另外,由于横向剪切引起的附加变形,Timoshenko裂纹梁的稳态挠度与Euler-Bernoulli梁挠度的差值为常数,其大小与裂纹模型、梁跨高比或裂纹深度无关,这些结果对梁裂纹无损检测具有指导意义．     

Some remarks on the propagation and non-propagation of discontinuities in linearly viscoelastic liquids

The equations of linear viscoelasticity with a bounded memory kernel have been shown to propagate singularities in a similar way as hyperbolic equations. In this paper, we investigate a model problem for a certain class of unbounded memory kernels. It is shown thatC -solutions are obtained, although there is a discontinuity in the boundary conditions.     

快速小波边界元的矩阵后压缩方法

介绍了一种基于传统边界元单元划分的小波 Galerkin 边界元法,该方法具有几乎线性(即 ON,N 为自由度)的求解复杂度。在准消失矩小波的框架下介绍了非标准型系数矩阵的压缩问题,提出了一种后压缩算法以降低小波边界元法的内存消耗。求解 Stokes 方程的算例表明,后压缩算法在保证结果收敛特性的情况下可以将系数矩阵的内存占用量降低 5 倍以上。     

任意外型薄板非线性后屈曲问题的边界元解

本文就薄板后屈曲问题建立一组新型的边界元计算公式，用这组公式求解能方便处理各种边界问题，另外文中将面内应力分解成基本部份和附加部份，并利用微分算子分解理论导得了挠度的一个不同形式的基本解，由于计算公式中，实现了面内位移和挠度的解耦，从而使迭代过程得到简化，文末还对圆板后屈曲路径进行了计算，得到了满意的结果。