构造极坐标中Airy应力函数的观察法

通过引入Airy应力函数，平面问题可以归结为在给定的边界条件下求解一个双调和方程．因此对双调和函数性质的研究将有利于平面问题的求解．首先给出一个有关双调和函数的引理，并分别从复变和微分两种角度提供该引理的证明．借助这个引理，提出了一种构造极坐标中Airy应力函数的观察法．最后，举例说明了该观察法在几个经典平面问题中的应用．这些例子说明，利用本的观察法可以将某些平面问题应力函数构造的过程简单化。     

用Level Set方法求解具有自由面的流动问题

为采用Ｌｅｖｅｌ Ｓｅｔ方法来计算有自由的流动问题提出了一种方案,把自由水面视为水和空气的交界面,两种介质用统一的Ｎ－Ｓ方法求解,在自由面两侧分别采用各自的密度和粘性,并在自由面上给以适当的光滑;采用边界元法求解双调和方程来确定距离函数;Ｎ－Ｓ方程用投影法求解,文中给出了二维水池水面振荡和瞬时溃坝问题的算例,可以看出用ＬｅｖｅｌＳｅｔ方法求解有自由面流动问题是有效的。     

带源参数的二维热传导反问题的无网格方法

利用无网格有限点法求解带源参数的二维热传导反问题，推导了相应的离散方程. 与 其它基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不 需要划分网格，用配点法离散控制方程，可以直接施加边界条件，不需要在区域内部求积分. 用有限点法求解二维热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点. 最后通过算例验证了该方法的有效性.     

弹性力学平面问题的应力函数的选择

弹性力学平面问题的应力解法,归之为求解满足边界条件的双调和方程.要从纯数学上来求出双调和方程的通解是很困难的,也是不必要的.所以弹性力学中不得不采用逆解法和半逆解法,来试凑出一个满足边界条件和双调和方程的解.但是,要从众多的函数中,选择一个既满足边界条件,又满足双调和方程的应力函数,谈何容易,这常使一些初学者感到束手无策.如果我们从边界上的已知应力分布规律出发,就很容易找到所需的应力函数了.例如,在直角坐标解法中,双调和方程为     

弹性力学平面问题的应力函数的选择

弹性力学平面问题的应力解法,归之为求解满足边界条件的双调和方程.要从纯数学上来求出双调和方程的通解是很困难的,也是不必要的.所以弹性力学中不得不采用逆解法和半逆解法,来试凑出一个满足边界条件和双调和方程的解.但是,要从众多的函数中,选择一个既满足边界条件,又满足双调和方程的应力函数,谈何容易,这常使一些初学者感到束手无策.如果我们从边界上的已知应力分布规律出发,就很容易找到所需的应力函数了.例如,在直角坐标解法中,双调和方程为...     

二维双曲型守恒律的一种五阶松弛格式

松弛格式是Jin和Xin提出的无振荡有限差分方法,其主要思想是将守恒律转化为松弛方程组进行求解.本文用逐维五阶WENO重构和显隐式Runge-Kutta方法对松弛方程组的空间和时间进行离散,得到了一种求解二维双曲型守恒律五阶松弛格式.所得格式保持了松弛格式简单的优点,不用求解Riemann问题和计算通量函数的雅可比矩阵.通过二维Burgers方程和二维浅水方程的数值算例验证了格式的有效性.     

弹性力学空间问题的复变函数解

弹性力学空间问题的复变函数解贾普荣（西北工业大学，西安７１００７２）１基本方程对于弹性力学一般的空间问题，位移和应力分量可用伽辽金的３个双调和函数表出￣［１］。设ｆ＿１（ｘ，ｙ，ｚ），ｆ＿２（ｘ，ｙ，ｚ），ｆ＿３（ｘ，ｙ，ｚ）为双调和函数，即位移分量...     

平面问题等价边界积分方程的三次边界轮廓法

基于弹性力学平面问题等的边界积分方程,给出了三次单元的边界轮廓法。根据平面问题解的复变函数表示,构造了三次形函数。给出了对于混合边值问题求解系统方程确定的边界轮廓方程配置和三次单元界轮廓法的实施。     

R-函数理论及准Green函数方法在正交各向异性板振动问题中的应用

首次将R-函数理论及准Green函数方法应用于求解固支正交各向异性薄板的自由振动问题。首先引入参数变换,将正交各向异性薄板的自由振动微分方程转化为双调和算子的边值问题,并应用R-函数理论,以解析函数形式描述复杂边界形状;利用问题的基本解和边界方程构造了一个准Green函数,该函数满足了问题的齐次边界条件;通过R-函数理论构造适当的边界方程,消除了积分方程核的奇异性;再采用Green公式将其化为第二类Fredholm积分方程。数值算例表明:该方法减少了理论计算量,精度较高。本文还证明了其优越性和正确性,是一种新型的数学方法。     

由泛复函构造弹性力学平面问题的特解

本文以泛复变函数（简称泛复函）为工具，通过引入双调和数，构造出直角坐标和极坐标下平面应力函数的一系列特解，其中有些是以往文献中尚未出现的。     

开孔薄板的弯曲问题

本文利用边界型最小二乘配点法,取极座标下的双调和方程的特解序列作为试函数,对开大孔的薄板弯曲问题进行了计算研究。计算中应用了复变函数保角变换的技术,得到了具有正多边形孔的圆板弯曲问题的近似解。这种方法可推广至其它形状的孔与其它形状的外边界。此外,在计算过程中提出了试函数逐项加权的方法,在一定程度上改善了矩阵病态的情况。     

广义平面应变问题

一般情况下平面应变问题归结为一个双调和方程、一个调和方程及轴向应变应满足平面坐标x、y 的线性函数条件下的求解.进一步分析证明,轴向载荷分布必须是x、y 的线性函数类问题都能按广义平面应变公式进行分析.     

矩形薄板弯曲问题的二维广义有限积分变换法

运用二维广义有限积分变换解法，本文推导出不同边界条件下矩形薄板弯曲问题的解析解．在推导过程中，选取满足边界条件的梁振型函数为广义积分变换的积分核，由此构造出广义有限积分变换对，通过对薄板弯曲问题的控制方程进行二维广义积分变换，可以将控制方程转换为易于求解的线性代数方程组．该方法无需预先选取位移函数，无需进行繁琐的叠加过程，求解过程思路清晰，说明该方法更加正确合理．最后通过计算实例对比，验证了该方法的合理性及所推导公式的正确性．     

光弹性分析中主应力分离的新方法

本文利用极坐标形式双调和方程通解得到应力第一不变量函数,结合光弹性实验等差线,用边界配置法求得满足求解对象边界条件的主应力和函数。于是由本文方法得到的主应力和函数和光弹性中的等差线,可求得物体内部各点的主应力值。     

最小二乘边界配点法解二维问题探讨

1.引言薄板弯曲问题和弹性力学平面问题的边界配点法,通常都是选用双调和方程的各类特解序列来构造试函数。在弹性力学教本中出现过的三类特解序列都可用来构造试函数。文献[2]以多项式级数为试函数,使用混合法分析过薄板强度。文[3]在边界配点法中使用了两种试函数——多项式级数和双曲三角函数。文[4]提出了一种产生双调和多项式的方法。文[5,6]使用的是双曲三角函数。文[7]综合了极坐标系中各种     

含有非圆形双孔的无限平板中应力的解析解研究

无限平板中含有任意形状单个孔的问题可以使用复变函数方法获得其应力解析解.对于无限平板中含有两个圆孔或两个椭圆孔的双连通域问题,也可以利用多种方法进行求解,比如双极坐标法、应力函数法、复变函数法以及施瓦茨交替法等.其中复变函数中的保角变换方法是获得应力解析解的一个重要方法.但目前尚未见到用此方法求解无限板中含有一个正方形孔和一个椭圆孔的问题.当板在无穷远处受有均布载荷和孔边作用垂直均布压力时,利用保角变换方法可以求解板中含有两个特定形状孔的问题.该方法将所讨论的区域映射成象平面里的一个圆环,其中最关键的一步是找出相应的映射函数.基于黎曼映射定理,提出了该映射函数一般形式,并利用最优化方法,找到了该问题的具体映射函数,然后通过孔边应力边界条件建立了求解两个解析函数的基本方程,获得了该问题的应力解析解.运用ANSYS有限单元法与结果进行了对比.研究了孔距、椭圆形孔大小和两孔布置方位对边界切向应力的影响,以及不同载荷下两孔中心线上应力分布规律.      

不可压N-S方程高效算法及二维槽道湍流分析

构造了基于非等距网格的迎风紧致格式,并将其与三阶精度的Adams半隐方法相结合,构造了求解不可压N－S方程高效算法。该算法利用基于交错网格的离散形式的压力Poisson方程求解压力项,解决了边界处的残余散度问题;同时还利用快速Fourier变换将方程的隐式部分解耦,离散后的代数方程组利用追赶法求解,大大减少了计算量。通过对二维槽道流动的数值模拟,证实了所构造的数值方法具有精度高,稳定性好,能抑制混淆误差等优点,同时具有很高的计算效率,是进行壁湍流直接数值模拟的有效方法。在数值模拟的基础上对二维槽道流动进行了分析,得到了Reynolds数从6000到15000的二维流动饱和态解（所谓“二维槽道湍流”）;定性及定量结果均与他人的数值计算结果吻合十分理想。对流场进行了分析,指出了“二维湍流”与三维湍流统计特性的区别。     

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

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

含有非圆形双孔的无限平板中应力的解析解研究

无限平板中含有任意形状单个孔的问题可以使用复变函数方法获得其应力解析解.对于无限平板中含有两个圆孔或两个椭圆孔的双连通域问题,也可以利用多种方法进行求解,比如双极坐标法、应力函数法、复变函数法以及施瓦茨交替法等.其中复变函数中的保角变换方法是获得应力解析解的一个重要方法.但目前尚未见到用此方法求解无限板中含有一个正方形孔和一个椭圆孔的问题.当板在无穷远处受有均布载荷和孔边作用垂直均布压力时,利用保角变换方法可以求解板中含有两个特定形状孔的问题.该方法将所讨论的区域映射成象平面里的一个圆环,其中最关键的一步是找出相应的映射函数.基于黎曼映射定理,提出了该映射函数一般形式,并利用最优化方法,找到了该问题的具体映射函数,然后通过孔边应力边界条件建立了求解两个解析函数的基本方程,获得了该问题的应力解析解.运用ANSYS有限单元法与结果进行了对比.研究了孔距、椭圆形孔大小和两孔布置方位对边界切向应力的影响,以及不同载荷下两孔中心线上应力分布规律.     

Blasius边界层的二次稳定性

采用人为中性方法使二维基本流达到有限幅值的平衡态,用小参数展开法求解二次稳定性方程.同时按 Floquet 理论研究了空间模式三维亚谐扰动的二次稳定性及幅值沿流向的演化.数值结果与实验结果一致.