首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 281 毫秒
1.
三维位势问题的边界元分析中,关于坐标变量的边界位势梯度的计算是一个困难的问题.已有一些方法着手解决这个问题,然而,这些方法需要复杂的理论推导和大量的数值计算.本文提出求解一般边界位势梯度边界积分方程的辅助边值问题法.该方法构造了与原边界值问题具有相同解域的辅助边值问题,该辅助边值问题具有已知解,因此通过求解此辅助边值问题,可获得梯度边界积分方程对应的系统矩阵,然后将此系统矩阵应用于求解原边值问题,求解过程非常简单,只需求解一个线性系统即可获得原边值问题的解.值得注意的是,在求解原边值问题时,不再需要重新计算系统矩阵,因此辅助边值问题法的效率并不很差.辅助边值问题法避免了强奇异积分的计算,具有数学理论简单、程序设计容易、计算精度高等优点,为坐标变量梯度边界积分方程的求解提供了一个新的途径. 3个标准的数值算例验证了方法的有效性.  相似文献   

2.
基于转化域方程为边界积分方程的极限定理及一个新颖的基本解分解技术, 建立间接变量规则化边界积分方程, 它有效地避免了奇异积分的直接计算. 与已有方法比,该方法不将问题变换为各向同性的问题去处理, 因而无需反演运算, 也有别于Galerkin方法, 无需计算重积分. 可计算任意边界位势梯度, 而不仅限于法向通量. 针对椭圆边界的边值问题, 提交一种精确单元来描述边界几何. 数值算例表明, 所提算法稳定且效率高, 所得数值结果与精确解吻合较好.   相似文献   

3.
基于移动相似中心的比例边界有限元方法   总被引:1,自引:0,他引:1  
在传统的比例边界有限元中,相似中心是固定的,难以用其求解关于偏心域的场问题。本文引入移动相似中心的概念,建立新的比例边界坐标变换,并利用加权余量法将控制方程半弱化为关于径向坐标的二阶常微分方程,引入对偶变量,将其降为系数矩阵为Hamilton矩阵的一阶常微分方程。对Hamilton矩阵进行Schur分解,得到微分方程的通解,代入边值条件可得关于积分常数的代数方程。此方法将比例边界有限元扩展到偏心域的边值问题,同时在径向是半解析的,解的精度高;仅需要离散求解域的一个边界,数据量小;在计算中仅需要对Hamilton矩阵进行Schur分解以及求解关于积分常数的代数方程,运算量少。将偏心环形域静电场边值问题的算例与解析解或其他数值方法计算结果的比较,表明此方法具有精度高、数据量小及运算量小的优点。  相似文献   

4.
弹性力学平面问题中一类无奇异边界积分方程   总被引:6,自引:2,他引:6  
从理论上提出一种新的方法,归化出间接变量无奇异边界积分方 程. 采用Lagrange二次单元,建立一个数值求解框架系统. 此外,基于问题的计算区域的 特殊性,给出一种边界近似方法. 数值算例表明该方法所取得的数值结果与精确解相当接 近,特别是边界量的数值结果. 此外,该方法容易被推广到三维问题. 和已有的直接变量的情形相比较,具有优点:1)无需处理HFP积分. 大大降低处理问 题的复杂性,并提高了计算效率和解的精度;2)摆脱了问题的具体形式,进入纯代数操作. 这样做的好处是从理论上建立一种普遍适用的方法,不仅适用于弹性力学问题,同样可应用 于其它问题,如位势问题, Stokes问题等. 3)提供了一种计算CPV积分的方法.  相似文献   

5.
本文以调和函数的边值问题为例,探讨了边界积分方程的充要条件.文中首次提出了超定问题的概念,并建立了超定问题有解的一个充要条件,它也就是直接变量边界积分方程的一个充要条件.文中首次阐明了边界积分方程与变分原理的内在的联系,还指出了间接变量与直接变量两类边界积分方程之间存在着一一对应的关系.文中的慨念、思路和论点不难用于其它有变分原理的问题的边界积分方程.  相似文献   

6.
非连续边界元积分的精确表达式及相关问题   总被引:5,自引:0,他引:5  
以二维位势问题边界元分析为例,给出了利用线性非连续边界元离散边界积分方程时系数矩阵积分计算的精确表达式,通过和利用Gauss积分方法计算系数矩阵所得数值结果的比较表明:配位点选择不同对数值计算结果精度影响的主要原因是积分计算的精度,尤其当配位因子选择较大时,存在的准奇异积分(Nearly Singular Integrals)很难利用常规Gauss积分方法准确求得。  相似文献   

7.
本文讨论了一类简化的Signorini问题。首先将原问题和一个边值问题建立联系,其次将原问题的解分解为不带不等边界条件的变分方程的解和一个变分不等式的解。然后利用边值问题的边界积分方程将变分不等式等价地化解为边界变分不等式。这样原求区域上的第一类椭圆变分不等式问题化解为求一个区域上的变分方程和一个边界变分不等式。最后说明了边界变分不等式解的存在唯一性。文末计算了柱面和半无限刚性基础的摩擦接触问题。结论表明文中方法具有较好的精度。  相似文献   

8.
对于热弹性平面问题,过去广泛集中在直接变量边界元法研究,本文研究间接变量规则化边界元法,建立了间接变量规则化边界积分方程。和直接边界元法相比,间接法具有降低密度函数的连续性要求、位移梯度方程中的热载荷体积分具有较弱奇异性等优点。数值实施中,用精确单元描述边界几何,不连续插值函数逼近边界量。算例表明,本文方法效率高,所得数值结果与精确解相当吻合。  相似文献   

9.
功能梯度材料动态断裂力学的径向积分边界元法   总被引:1,自引:0,他引:1  
高效伟  郑保敬  刘健 《力学学报》2015,47(5):868-873
采用径向积分边界元法分析功能梯度材料动态断裂力学问题. 该方法使用与弹性模量无关的弹性静力学开尔文基本解作为问题的基本解,在导出的边界-域积分方程中含有由材料的非均质性和惯性项引起的域积分,通过径向积分法将域积分转化为等效的边界积分,得到只含边界积分的纯边界积分方程;从而建立只需边界离散的无内部网格边界元算法. 采用候博特方法求解关于时间二阶导数的系统离散的常微分方程组. 最后通过数值算例验证本文方法的精度和有效性.   相似文献   

10.
薄体位势问题边界元法中的解析积分算法   总被引:1,自引:0,他引:1  
薄体结构的数值分析是边界元法的难点问题之一。该文导出了一种完全解析积分算法,用这种算法计算了薄体平面位势问题边界元法中出现的几乎弱奇异、强奇异和超奇异积分。当边界离散为一系列线性单元,边界积分方程离散计算的积分可归纳为三种形式。对薄体问题,源点与积分单元距离通常相距很近,这些积分产生显著几乎奇异性,直接采用常规高斯积分不能有效计算。为此该文导出了这些几乎奇异积分的全解析计算公式。按源点与单元的距离是否为零,公式分两种情况。新算法采用全解析积分公式处理几乎奇异积分,首先精确计算出薄体问题边界未知位势和法向位势梯度,然后再进一步计算了域内点的物理参量。算例表明该文算法可处理狭长比为1.E-08的薄体问题,显示了边界元法分析薄体问题具有独特的优势。  相似文献   

11.
伸缩虚拟边界元法解二维Helmholtz外问题   总被引:1,自引:0,他引:1  
向宇  黄玉盈 《力学学报》2003,35(3):272-279
以位势理论为基础,提出了求解Helmholtz外问题的伸缩虚拟边界元法.给出了该方法在全波数域内获得唯一解的严格数学证明,其核心是通过伸缩虚拟边界使对偶内问题的特征频率(本征值)避开与波数重合,从而保证了解的唯一性,同以往前人提出的几种解法途径相比,该法简单得多;通过诸多边界曲线形状和不同边界量的声辐射算例,从计算精度、稳定性以及克服解的非唯一性等方面,对该方法进行了检验.计算结果表明:对远场或近场辐射声压,该方法都具有非常高的效率和精度.  相似文献   

12.
An analytical approach is proposed here to study scattering of deep water waves by a submerged or a surface piercing vertical porous barrier. It involves a connection between two wave potentials of which one is the solution of a boundary value problem associated with wave scattering by the porous barrier and the other is the solution of a complementary type problem where barrier and gap positions are interchanged and solid barrier takes the position of the porous barrier. The connection also involves an auxiliary or a connection wave potential. The potential for the solid barrier problem involves incident wave forcing while the auxiliary potential describes a solid barrier type problem that involves a non-physical forcing. The solution procedure of Ursell (Ursell, 1947) is chosen to solve these boundary value problems explicitly in the case of normal wave incidence as it also determines necessary exact behavior of the potential at the barrier edge. The reflection coefficients are also connected and the reflection amplitudes of the normally incident wave against the vertical porous barriers are obtained analytically. Numerical results for reflection and transmission coefficients are presented.  相似文献   

13.
提出了将杂交边界点法和双重互易法结合求解势问题的一种新的算法. 将势问题的解分为通解和特解两部分,通解使用 杂交边界点方法求解,特解则利用局部径向基函数近似. 该方法输入数据只是求解域上离散 的点,不需要额外的方程来计算域内物理量,后处理十分简便. 数值算例表明了该方法的稳 定性和有效性.  相似文献   

14.
在同一界面的不同区域具有多种边界条件, 称之为混合边界, 这是一个熟知的力学问题. 对这类问题进行精确分析时, 必须要进行混合边值问题的求解. 而对于一般的三维非轴对称情形, 混合边值问题的求解往往存在数学困难. 本文利用Hilbert定理和双重Fourier变换, 给出了一种求解三维非轴对称混合边值问题的解析方法, 利用该方法对具有混合透水边界的饱和多孔地基上矩形板的振动弯曲进行了解析研究(板与地基接触面为不透水边界, 其余为透水边界). 首先, 基于Kirchhoff理论和Biot多孔介质理论建立矩形板与饱和多孔地基的动力控制方程, 进行耦合求解. 针对板土接触面和非接触面的混合边值问题, 采用双重Fourier变换构造出两对二维对偶积分方程, 以接触应力和接触面孔隙压力为基本未知量, 用Jacobi正交多项式将未知量展开, 再利用Schmidt法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.  相似文献   

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

16.
王熙  龚育宁 《力学学报》1992,24(1):93-101
本文给出了弹性动力学轴对称问题基本方程的一种理论解。它由满足非齐次边界条件的准静态解和满足齐次边界条件的动态解的叠加构成。在求得准静态解后,代入基本方程,得到动态解所需满足的非齐次方程。由相应的齐次方程的特征值问题,定义了有限Hankel变换。通过这种变换及Laplace变换,求得动态解,从而得到了一个完整的理论解。文中通过对一个实例求解,表明该方法求解过程简便,实用,求解结果精确。  相似文献   

17.
The paper combines a self-adaptive precise algorithm in the time domain with Meshless Element Free Galerkin Method (EFGM) for solving viscoelastic problems with rotationally periodic symmetry. By expanding variables at a discretized time interval, the variations of variables can be described more precisely, and iteration is not required for non-linear cases. A space-time domain coupled problem with initial and boundary values can be converted into a series of linear recursive boundary value problems, which are solved by a group theory based on EFGM. It has been proved that the coefficient matrix of the global EFG equation for a rotationally periodic system is block-circulant so long as a kind of symmetry-adapted reference coordinate system is adopted, and then a partitioning algorithm for facilitating parallel processing was proposed via a completely orthogonal group transformation. Therefore instead of solving the original system, only a series of independent small sub-problems need to be solved, leading to computational convenience and a higher computing efficiency. Numerical examples are given to illustrate the full advantages of the proposed algorithm. The project supported by the National Natural Science Foundation of China ((10421002, 10472019 and 10172024); NKBRSF (2005CB321704) and the Fund of Disciplines Leaders of Young and Middle Age Faculty in Colleges of Liaoning Province. The English text was polished by Yunming Chen.  相似文献   

18.
薛冰寒  林皋  胡志强  庞林 《力学学报》2016,48(3):615-623
摩擦接触问题是计算力学领域最具挑战性的问题之一,接触系统的泛函具有非线性、非光滑的特点,导致接触算法的收敛性与精确性难以保证.因此将比例边界等几何分析(scaled boundary isogeometric analysis,SBIGA)与B可微方程组(B dierential equation,BDE)相结合,提出了求解二维摩擦接触问题的比例边界等几何B可微方程组方法.在比例边界等几何坐标变换的基础上,通过虚功原理推导了关于边界控制点变量的接触平衡方程,表示成B可微方程组形式的接触条件可被严格满足,求解B可微方程组的算法的收敛性有理论保证.此比例边界等几何B可微方程组方法(SBIGA-BDE)只需在接触体边界进行等几何离散,使问题降低一维,能精确描述接触边界,并可通过节点插入算法进行真实接触区域的识别.此外,由于几何建模和数值分析使用相同的基函数,节约了划分网格的时间.以赫兹接触问题和悬臂梁摩擦接触问题为例,通过与解析解及数值计算软件ANSYS计算结果进行对比,验证了该方法求解二维摩擦接触问题的有效性及高精度等特点.   相似文献   

19.
李珺璞  陈文 《力学学报》2018,50(4):961-969
大规模高频声场的数值模拟是一项非常有计算挑战性的课题. 为了解决传统边界型离散方法由于全局支撑的满阵限制, 不易应用于大规模高频声场模拟的计算瓶颈, 本文提出了一种用于模拟大规模高频声场的双层奇异边界法. 在该方法中, 通过引入双层结构, 细网格上的全局支撑的满阵被转化为局部支撑的大规模稀疏矩阵, 传统奇异边界法模拟大规模问题时所面临的高计算量以及过度存储需求遂得以解决. 其次, 双层奇异边界法仅通过粗网格评估远场作用, 且独立于特定的插值核函数. 相较于快速多级方法, 该方法具有更强的适应性和灵活性, 且多层结构使该方法具有一定的预调节作用, 非常适合求解具有大规模、高秩、高条件数特点的高频波矩阵. 在其后的散射球模型算例中, 双层奇异边界法配置10万个节点, 成功模拟了无量纲波数高达160的声散射问题. 在对于人头模型的声散射特性分析中, 双层奇异边界法比COMSOL软件计算速度快了约78.13%. 当配置8万个节点时, 双层奇异边界法成功模拟了频率高达25 kHz 的工况, 该频率已远远超出了人耳的听力极限.   相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号