利用边界单元法求平面问题应力强度因子K_I的应力-位移杂交法

本文提出用裂尖附近2点或3点的应力和位移计算应力强度因子K_I的杂交方法.这种方法充分利用了边界单元法的计算结果,考虑了裂尖应力场和位移场渐近展开式的高阶项,使用远离裂尖的点算出的K_I也有较好的精度,拟合线十分平坦.用算例的结果将杂交法与一般的位移法和应力法进行了比较,同时,对常量单元和线性单元也进行了比较.     

一种改进格林元方法及在渗流问题中的应用

采用格林公式和基本解推导出直接边界积分方程来求解渗流问题.边界积分方程数值离散基于格林元方法(Green element methond),改进了原方法中压力和压力导数的求解方法,命名为混合边界元方法(Mixed boundary element method).相较于格林元类方法,该方法显式考虑了求解节点的外法向流量值和压力值,并使求得的数值解在求解区域上能够连续,符合实际的物理过程,在不增加额外未知数的情况下提高了计算精度.分析了不同网格类型对模拟计算结果的影响,并对稳定渗流问题、非稳定(瞬态)渗流问题和非稳态问题进行了实例计算,结果显示改进方法提高了计算精度,并对各类渗流问题有较好的适应性.     

快速求解叶栅流动的边界元方法

本文利用边界元方法，通过ｋｒｉｃｈｈｏｆｆ变换将描述叶栅流动的控制方程转换成线性方程。并将广义ｋ－Ｊ条件与边界积分方程联立求解，避免了非线性项和叶片出气角的迭代计算。完成了一种快速求解任意迥转面叶栅流场的计算程序，实用表明与其它数值方法及实验结果符合较好，具有快带、简明、实用的特点。     

Boundary element analysis of contact problems using artificial boundary node approach

An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.     

黏弹性人工边界等效荷载计算的改进方法

黏弹性人工边界在场地地震反应和结构-地基动力相互作用等问题的计算中已得到了广泛的应用.地震波在黏弹性人工边界中的输入是通过将地震波转化为作用于人工边界处的等效载荷来实现的.计算等效节点载荷的常规方法默认边界节点对应区域内的应力为均布力,但实际上该节点对应区域内的应力分布通常是不均匀的.本文在有限元方法结合黏弹性局部人工边界的显式时域波动方法的基础上,建立了无限域散射问题地震波等效载荷计算的一种改进方法.该方法采用细化网格与应力积分相结合的方法计算人工边界等效节点力,有效地降低了人工边界上等效节点力的计算误差.以不同角度入射地震波的二维算例为例,算例给出的波场位移云图和节点位移时程曲线验证了本文方法的有效性,其计算精度与网格尺寸和地震波入射角度密切相关,且网格越小、入射角度越小,计算精度越高.对于相同的网格尺寸,本文采用方法的计算精度明显高于常规方法,尤其是对于斜入射问题优势更为明显.     

敷设多孔吸声材料声腔特征值分析的径向积分边界元法

由于Helmholtz方程的基本解是频率的函数,因此传统边界元法在处理声场特征值问题时具有天生的缺陷。本文采用Laplace方程基本解生成积分方程,通过径向积分法将在此过程中产生的域积分项转化为边界积分。此方法克服了传统边界元法系数矩阵对频率的依赖,同时克服了特解积分法对特解的依赖,并通过对表面声导纳的多项式逼近,将敷设多孔吸声材料声腔特征值问题转化为矩阵多项式,从而避免了复杂的非线性求解。通过数值算例验证了算法的有效性。     

SOLVING THE FREE BOUNDARY PROBLEM IN CONTINUOUS CASTING BY USING BOUNDARY ELEMENT METHOD

ＳＯＬＶＩＮＧＴＨＥＦＲＥＥＢＯＵＮＤＡＲＹＰＲＯＢＬＥＭＩＮＣＯＮＴＩＮＵＯＵＳＣＡＳＴＩＮＧＢＹＵＳＩＮＧＢＯＵＮＤＡＲＹＥＬＥＭＥＮＴＭＥＴＨＯＤＬｉＹａｏｙｏｎｇ（李耀勇）；ＺｈａｎｇＺｈｉｌｉ（张自立）（ＲｅｃｅｉｖｅｄＪｕｎｅ，１８，１９...     

动态断裂力学的无限相似边界元法

对弹性动力学的相似边界元法进行了进一步研究,推导了相应的计算公式,并在此基础上提出了动态断裂力学的无限相似边界元法.与传统的边界元法相比,相似边界元法由于只需在少数单元上进行数值积分,大大减少了计算量.对动态断裂力学问题,无限相似边界元法由于在裂纹尖端的边界上设置了逼近于裂纹尖端的无限个相似边界单元,可直接得到裂纹尖端具有奇异性的应力,而不需要设置奇异单元,从而突破了奇异单元对应力奇异性阶次的局限.另外,还讨论了无限相似边界元法得到的无限阶的线性代数方程组的求解方法.     

三维势流场的比例边界有限元求解方法

比例边界有限元法(SBFEM)是线性偏微分方程的一种新的数值求解方法。该方法只对计算域边界利用Galerkin方法进行数值离散,相对于有限元方法(FEM)减少了一个空间坐标的维数,而在减少的空间坐标方向利用解析方法进行求解;相对于边界元法(BEM),比例边界有限元方法不需要基本解,避免了奇异积分的计算,所以它结合了有限元和边界元方法的优点。本文建立了利用比例边界有限元法求解三维Laplace方程的数值模型并用于计算三维物体周围的水流场,将计算结果与解析解和边界元方法进行了对比,结果表明此方法可以很好地模拟水流场,且具有较高的计算精度。     

Research on the companion solution for a thin plate in the meshless local boundary integral equation method

The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method ( GFEM ), boundary element method (BEM) and element free Galerkin method (EFGM), and is a truly meshless method possessing wide prospects in engineeringapplications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.     

冲击波对核反应堆安全壳的动力响应研究

本文用时域边界积分方程提出研究半球形核反应堆安全壳在爆炸冲击波作用下动力响应分析的计算模型。从空中爆炸和地面爆炸两方面探求壳体位移和变形过程,分析冲击波在壳体内的传播规律。计算了在水平冲击波作用下壳体的内力分布,其计算结果可以作为我国核电站设计的参考分析依据。     

正交各向异性平面断裂问题应力强度因子的边界元解法

本文研究单向复合材料或正交各向异性体平面断裂问题,构造了一个1/4节点混合参数应力奇异边界元,综合运用该元素与1/4节点等参边界元,提出了求解应力强度因子的混合边界元解法,用所述方法计算了含中心裂纹无限大与有限大正交各向异性板的应力强度因子,算例表明,本文所述方法不仅计算精度高,而且适应性强,便于工程应用。     

模拟相变问题的精细积分边界元法

将精细积分边界元法和界面追踪法相结合求解相变问题。因为边界元法只需要将待求解空间域的边界离散,方便连续追踪移动界面位置和重构网格,所以边界元法适合应用于移动边界问题的模拟。首先,利用精细积分边界元法在固相区域和液相区域分别求解相应的瞬态热传导控制方程,从而求得温度场和边界热流密度。然后,根据固-液相变界面上的能量平衡方程,利用热流密度求得相变界面的移动速度,再采用界面追踪法预测移动相变界面的位置变化。最后,给出了几个数值算例,并通过与参考解的对比验证本文方法的准确性。     

NUMERICAL SIMULATION OF 2D FIBER-REINFORCED COMPOSITES USING BOUNDARY ELEMENT METHOD

Introduction Withthedevelopmentofmodernindustry,compositesareincreasinglybeingappliedto agreatnumberofimportantstructures.Todeterminethemacroscopicaleffective characteristicsofcompositesisanessentialprobleminmanyengineeringapplications.The macroscopicalef…     

压电材料平面问题的虚边界元-等额配点解法

利用压电材料平面问题的基本解和弹性力学虚边界元方法的基本思想,提出了压电材料平面问题的虚边界元-等额配点解法。该解法继承了传统边界元方法的优点,而避免了传统边界元方法遇到的边界积分奇异性问题。最后给出了压电材料平面问题的一些具体算例,并与解析解作了比较。结果表明本文的方法有很高的精度,是该问题一个十分有效的数值求解方法。     

裂缝分析的比例边界有限元与有限元耦合的虚拟结构面模型

比例边界有限元法SBFEM(Scaled Boundary Finite Element Method)是一种半解析数值方法,在裂缝分析特别是强度因子计算上具有相当高的精度。本文提出了一种用于裂缝分析的基于虚拟结构面的SBFEM与常规FEM的耦合分析方法。首先选取裂缝周边一定范围的计算域,并将结构分成不含裂缝区域和含裂缝区域两部分。然后,对不含裂缝区域,采用FEM进行网格离散;对含裂缝区域,采用SBFEM进行网格离散;两者相互独立,在这两个域内,分别采用各自相应的位移模式。最后通过在SBFEM网格的外边界设置虚拟耦合结构面的模式,实现有限元网格和比例边界有限元网格的耦合。通过两个经典的含裂缝平板的算例研究,探讨了本文方法在I型开裂和混合型开裂分析中,影响应力强度因子精度的因素。算例表明,SBFEM具有的降维和半解析性质,使本文方法在裂缝分析中的前处理简单易行,且计算结果具有相当高的计算精度。     

圆外Stokes问题的边界积分公式

针对圆外区域的Stokes问题,利用Fourier展开法,通过自然边界归化得到了一个只与边界速度有关的Stokes问题的边界积分公式.根据此公式及边界速度值,求得区域内速度及压力分布的解析表达式,并通过数值积分的方法进行求解,计算量小,所得速度及压力分布图曲线光滑.最后借助流体软件进行数值计算,结果验证了边界积分公式的正确性、可行性.     

双材料基本解弹塑性边界元法

提出并研究采用双材料基本解的弹塑性边界元法,得到了内点应力公式中有关奇点塑性应变自由项的完整表达式,并利用非连续边界单元和非连续区域单元解决了当奇点位于界面上时该自由项难于确定,以及计算区域Cauchy主值积分的常塑性应变场法在与界面相连的奇异区域单元上无法实施的困难．采用双材料基本解的弹塑性边界无法针对双材料的结构特点,特别适于分析有关弹塑性双材料界面及界面裂纹问题．     

边界元法计算近边界点参量的一个通用算法

针对边界元法存在近力界点参量计算的困难,给出了一个通用性方法,将近边界点到边界单元的距离参数通过分部积分变换到积分式之外,从而计算出二维问题近边界点参量的几乎强奇异和超奇异积分,由此,对任何近边界点参量,提出了整套计算方案,算例证明了本法的有效性。     

几何精确NURBS有限元中边界条件施加方式对精度影响的三维计算分析

非均匀有理B样条(NURBS)有限元法把计算机辅助几何设计(CAGD)中的NURBS几何构形方法与有限元方法有机结合起来,有效消除了有限元离散模型的几何误差,提高了计算精度。但是由于NURBS基函数不是插值函数,直接在控制节点上施加位移边界条件会引起较大误差。本文详细讨论了NURBS基函数的插值特性,在NURBS有限元分析中采用罚函数法施加位移边界条件,提高了收敛率和计算精度。结合典型三维弹性力学问题,对两种施加位移边界条件的方法进行了对比和分析。计算结果表明,直接施加位移边界条件会导致收敛率和精度的明显降低,而基于罚函数法的NURBS有限元分析则能达到最优收敛率,并具有更高的精度。