弹性力学问题的一种新解法

弹性力学问题的一种新解法邹广德，沈玉凤（山东工程学院，淄博２５５０１２）对于弹性力学问题，常常采用边界积分方程的间接解法，将假想作用在边界上的虚拟力作为未知量，故也称虚拟力法。这种方法与边界积分方程的直接解法一样，存在奇异积分，且边界附近精度较低。为...     

弹性力学问题的一阶有限元解法

求解弹性力学问题的应力时,如果采用常规的位移有限元法,需要先求得单元的节点位移,再经过求导运算得到。为了解决这种求解方式引起的应力精度下降的问题,提出了弹性力学问题的一阶多变量形式,使得应力与位移精度同阶,并推导了弱形式。采用有限元方法,对弹性力学问题给出了一阶解法的二维、三维数值算例,并且将一阶解法的结果与常规位移有限元法的解进行了比较。数值计算结果表明,一阶解法有效提高了应力的精度,并且应力的误差和节点位移的误差具有相同的收敛阶,验证了本文方法的有效性,为提高有限元法的应力精度提供了新的思路。     

弹性力学平面问题的新解法—应力分量法

本文提出了一种求解弹性力学平面问题的新方法－应力分量法，并对经典的平面问题进行了求解，解答正确，方法便捷，与传统的应力函数法相比，难度小，费时少，速度快。     

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

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

弹性力学平面问题的新解法─—应力分量法

本文提出了一种求解弹性力学平面问题的新方法─—应力分量法,并对经典的平面问题进行了求解,解答正确、方法便捷,与传统的应力函数法相比,难度小、费时少、速度快。     

Kelvin问题的一种新解法

1.引言在静弹性力学的边界元法中,Kelvin问题的解是一种非常重要的基本解.自1848年Lord Kelvin首先用分量形式的Helmhoty分解求出作用在一无限大域内的集中力产生的位移和应力场的特解以来,文献[2]、[3]又用矢量分解的形式进行了求解.其实矢量分解的方法只不过是Lord Kelvin提出的解法的一种冠以矢量算子的简洁书写形式而已.本文将要采用一种全新的方法,直接推导出二维和三维Kelvin问题与边界元法相应的基本解函数.与其它方法相比,该法新颖,且尤为简洁、方便.     

弹性力学的一种正交关系

在弹性力学求解新体系中，将对偶向量进行重新排序后，提出了一种新的对偶微分矩阵，对于有一个方向正交的各向异性材料的三维弹性力学问题发现了一种新的正交关系．将材料的正交方向取为z轴，证明了这种正交关系的成立．对于z方向材料正交的各向异性弹性力学问题，新的正交关系包含弹性力学求解新体系提出的正交关系。     

弹性地基梁的一种解法

由于弹性地基梁与地基连续接触,因而成为一无穷多次超静定结构,且又必须考虑地基的变形,这就使问题的求解相当冗繁.尤其在受载复杂时,计算更感不便.本文通过积分变换,提出一种解法,使其在公式推导与解题程序上均有所简化,并能对长梁、短梁进行统一描述.(一)基本公式推导考察图1所示弹性地基梁,假定全梁抗弯刚度 EJ为常数,并认为地基符合温克尔假设.这样,该梁的地 ....     

利用计算机辅助求解弹性力学问题的一种新方法

采用Matlab编程与多项式应力函数相结合的方式直接求解弹性力学平面问题. 该方法一方面使选择应力函数具有非常强的规律性,另一方面简化了确定待定常数的过程. 该方法改进了用传统的逆解法、半逆解法选择应力函数,需要经过多次尝试、不断修改,才能得到满意结果的缺点. 同时指出了在一些情况下弹性力学平面问题是不具有多项式的解答.     

功能梯度材料平面问题的辛弹性力学解法

将辛弹性力学解法推广用于功能梯度材料平面问题的分析,考虑沿长度方向弹性模量为指数函数变化而泊松比为常数的矩形域平面弹性问题,给出了具体的求解步骤. 提出了移位Hamilton矩阵的新概念,建立起相应的辛共轭正交关系;导出了对应特殊本征值的本征解,发现材料的非均匀特性使特殊本征解的形式发生明显的变化.      

弹性力学问题复变函数解法的应用与发展

对利用复变函数法求解弹性力学问题在过去100年中的发展进行了简要综述,介绍了自1909年Kolosov提出复应力函数法以来的弹性力学复变函数解法相关的几本具有影响的经典专著或教科书、复变函数方法在弹性问题壳体和空间问题的拓展,以及近些年来在解决更多力学问题上的新发展和新应用等.     

球面各向同性弹性力学的位移解法

本文引入三个位移函数（ｗ，Ｇ，ψ），将球面各向同性弹性力学运动方程，简化为关于ψ的二阶偏微分方程，和关于Ｗ和Ｇ的联立方程。在静力学问题中，联立方程可进一步简化，ｗ和Ｇ可用另一位移函数Ｆ表示，而Ｆ满足一个四阶偏微分方程。在球壳固有振动问题中，则简化为一个独立的二阶常微分方程，和另两个二阶的联立的常微分方程，证明了在多层球壳中，它们分别对应独立的两类振动。改进了常微分方程的解法，并计算了一个二层球壳的频率。     

弹性力学混合状态方程的小波解法

应用小波理论求解弹性力学混合状态方程，讨论了解的收敛性。从文中的数值算例不难看出，该方法不失为混合状态方程一种新的求解途径。     

弹性力学空间轴对称问题的一种复变函数法

本文引用了复变量广义解析函数概念,证明了空间轴对称问题的拉甫应力函数可以用两个适当选择的广义解析函数表示,据此即可导出应力分量、位移分量及边界条件的复变函数表示式,进而就可利用复变函数法求解空间轴对称问题。本文用这种方法求解了含有一个球形空腔的圆轴在两端受拉时的解答,表明了以此法分析求解空间轴对称问题的可能性。     

一种弹性波反演的快速迭代解法

本文研究弹性波反演的数值解法,通过求低频响应值的残差最小和利用有限元及振型分解法计算弹性波正问题的低频响应值,给出了一种弹性波反演的快速迭代解法。该算法允许引入反演参数的约束条件,计算效率很高;尽管一般只能得到弹性波反问题的近似解,但收敛性很好,因此可应用于为其它弹性波反演解法提供较好的初始值。     

Kelvin问题的一种解法

1.Kelvin问题所对应的定解问题在无界弹性体内,把集中力的作用点选为坐标系Ox_1x_2x_3的原点.设集中力的大小为P(常数),其方向沿单位矢量n,则此集中力可表示为Pδ(x)n,其中δ(x)为Dirac的δ函数,x为空间中一点的位置矢.弹性体内各点的位移矢量u是点的坐标的函数,表示为u=u(x).位移u在无界域中所应满足静     

弹性力学的一种边界无单元法

首先对移动最小二乘副近法进行了研究，针对其容易形成病态方程的缺点，提出了以带权的正交函数作为基函数的方法－改进的移动最小二乘副近法，改进的移动最小二乘逼近法比原方法计算量小，精度高，且不会形成病态方程组，然后，将弹性力学的边界积分方程方法与改进的移动最小二乘逼近法结合，提出了弹性力学的一种边界无单元法，这种边界无单元法法是边界积分方程的无网格方法，与原有的边界积分方程的无网格方法相比，该方法直接采用节点变量的真实解为基本未知量，是边界积分方程无网格方法的直接解法，更容易引入界条件，且具有更高的精度，最后给出了弹性力学的边界无单元法的数值算例，并与原有的边界积分方程的无网格方法进行了较为详细的比较和讨论。     

薄壳问题的一种数值解法

本文直接从三维弹性力学微分方程出发,依据三维的Kelvin解,应用最小二乘法建立了三维虚边界元法解薄壳问题的一般方法。本方法的显著优点是：不论求解何种壳体问题,思想是不变的,均以三维的Kelvin解来建立方程,而勿需对不同几何形状的壳本采用不同的基本解。文中给出了数值算例,以作为本方法的应用。本文方法与边界元直接法相比,优点在于无需处理奇异积分,且系数阵是对称的;再者,本文方法思想简单,程序实现容易。     

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

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