MLPG法与配点法耦合求解地下水非均质承压稳定流问题

给出了无网格局部Petrov-Galerkin(MLPG)法与配点法耦合求解导水系数为分片常数的非均质承压稳定流问题的方法.在各子区分界线布置的节点上,应用相容条件建立配点方程组,在除分界线外的其他节点上建立MLPG方程组,联立得到求解水头函数数值解的耦合方程组.编写了相应MATLAB程序,进行了具体模型计算,并与MLPG法和边界元法的计算结果进行了比较,结果表明该方法求解问题有效,精度较MLPG法计算精度显著提高,且明显优于边界元法.     

采用断续附设边界的间接边界单元法解位势问题

本文提出了在间接边界单元法界采用断续的附设边界来求解位势问题.该方法避免了奇异积分,简化了计算,同通常的奇异间接法相比,它大大地提高了计算精度.     

导体角域静电问题的直接法与保角变换法

从几个典型导体角域静电问题示例出发阐明直接法与保角变换法的关系,指出保变换法仅仅是一种特殊的变 量代换法;保角变换具有局限性,不宜使用保角变换求解导体角域点电荷静电问题.     

有限元与边界元耦合法在波浪力计算中的应用

利用有限元与边界元耦合法对三维无界区域中直立圆柱所受的波浪力进行计算,把整个求解区域分成内域和外域两部分,在内域采用有限元法,对外域采用边界元法计算,由加权余量法的理论知这种耦合在理论上是可行的,根据此耦合方法编制了完整的计算机程序,对海洋中直立圆柱的波浪力进行了分析.数值计算的结果与理论解吻合良好,表明该方法有效.     

求解类拉普拉斯方程的径向基函数配点法

建立了径向基函数配点法求解类拉普拉斯方程定解问题的方法.将求解域依据系数张量为分片常量来分解为若干子域.在每个子域上分别利用所布置的中心点建立用径向基函数表达的近似待解函数.在每个子域内及子域边界与外边界重合部分的配置节点上分别利用类拉普拉斯方程和定解条件建立近似解函数的待定系数满足的配点方程组,在相邻子域的分界线的配置节点上利用相容条件建立待定系数满足的配点方程组.方程组联立求解就得到了整个域上的近似解.算例计算表明该方法求解问题简捷有效,算法具有很高的精度.     

无侧移结构的广义力矩分配法与位移法联合法

在广义力矩分配法的基础上,引入广义转动刚度系数的概念,建立广义力矩分配法和位移法的联合法,求解得到精确的计算结果.首先,在所有未知转角位移刚节点处施加刚臂,计算得到杆端的固端弯矩,释放尽量多的结点,计算得到未分配结点处的不平衡力矩;然后,以未分配结点处的转角位移为基本未知量,建立位移法典型方程,再利用广义力矩分配法计算得到广义转动刚度系数;最后,求解计算典型方程,利用叠加法绘制最终弯矩图.以多跨连续梁和多层框架作为算例,演示联合法的具体计算步骤,并对比分析矩阵位移法、联合法和分层法的计算结果,验证了联合法的准确性和有效性.     

非线性热传导问题的基于滑动Kriging插值的MLPG法

利用基于滑动Kriging插值的无网格局部Petrov-Galerkin(MLPG)法来求解二维非线性稳态和瞬态热传导问题,Heaviside分段函数作为局部弱形式的权函数,并通过加权余量法推导相应的离散方程.该问题考虑了材料热传导系数随温度的线性变化,并通过拟线性法来求解非线性问题的解,时间域的离散通过向后差分法来实现.基于滑动Kriging插值构造MLPG中的形函数由于满足克罗内克δ性质,因此可以直接准确地施加本质边界条件.在构造刚度矩阵过程中,只涉及边界积分,不涉及区域积分和奇异积分.将数值计算结果与有限元法得到的结果加以对比可以看出,基于滑动Kriging插值的MLPG法能够很好地解决此类热传导问题.     

用分域边界元计算微带内导体同轴传输线的特性阻抗

文章根据明晰的物理概念，提出用分域边界元法计算微带内导体同轴传输线的特性阻抗，将无法用直接边界法进行直接求解的微带内异导体传输线问题，转化为可以用边界元法求解的简单边值问题，并通过三个工程实例检验其计算精度。实例计算结果表明，分域边界元法具有较高的计算精度，可以很方便地应用于微带内导体传输线问题的工程设计与计算。     

轴向运动Timoshenko固支梁固有频率的数值仿真

运用Galerkin方法和微分求积法求解固支边界轴向运动Timoshenko梁的固有频率.讨论系统的前两阶固有频率随轴向速度、刚度系数变化的情况,并将这2种方法得到的数值计算结果与复模态分析方法得到的精确解进行比较,发现用微分求积法和复模态分析法得到的结果几乎重合,而用Galerkin方法得到的结果在随刚度系数的增加和速度的增大时有所差异.     

追赶法并行求解循环三对角方程组

给出了求解循环三对角线性方程组的一种并行算法.在系数矩阵满足对角占优的条件下,利用该方法能够快速、稳定地求解循环三对角线性方程组,在单个进程上的计算量仅为○(17n).与传统算法求解循环三对角线性方程组的计算量相同.而且,本算法可以方便地实施分布式并行计算,各进程仅需向主进程传递8个实数,而主进程向各子进程传递2个实数,通讯量较小.数值实验结果表明:对于大规模的循环三对角线性方程组.利用16个进程计算的并行效率均在0_75以上.求解三对角线性方程组的传统追赶法实则是本文算法的一种特例,因此.该算法也可用于求解三对角线性方程组.     

关于边界元法及其解的精度问题探讨

本文阐述了边界元直接法与间接法的基本原理及步骤,提出了改进边界元解在边界附近精度的方法。     

间接边界元法的物理表述

根据两类位势层的物理意义,得到了边界元法的直接法与间接法的统一公式.不同学科中采用的间接法在表达式的形式上是一致的.因此,有共同的原理。对于同一物理现象,采用不同性质的点源,可以得到不同的基本解,从而,形成间接法中的两种不同方法.经过统一处理,还可编制实用的计算机程序.     

弹性力学问题的第二类边界积分方程——边界元法及其二维基本解

本文导出了以单位集中不连续位移（第二类边界）产生的场为基本解的两种边界积 分方程——直接法和间接法边界积分方程，按 Ｈａｄａｍａｒｄ主值存在，进行高阶奇异积 分并给出了二维问题的两个基本解，与通常的第一类直接法和间接法相比，第二类方法 对处理裂纹等不连续结构面较为有利。     

侧向支柱对冲压增程弹用进气道性能影响研究

结合分区对接网格技术和二阶精度区域分解算法,对某冲压增程弹丸进气道在不同来流攻角和不同侧向支柱形状工况下的内外复杂流场进行了数值模拟.得到了临界工况下超声速进气道内外粘性流场复杂的波系结构,分析了流场结构特性和激波波系结构.当进气道扩压段有侧向支柱存在时,进气道总压恢复系数和动能系数均有所降低,而流场畸变指数则显著增大.在三种形状侧向支柱中,采用两端削尖支柱形状的进气道的性能最优.攻角的存在也在一定程度上降低了进气道性能.     

Comparison study on spherical wave superposition method and spherical wave source boundary point method for realizing nearfield acoustic holography

In the light of the concept of spherical wave source, the theoretical model of nearfield acoustic holography (NAH) based on the spherical wave superposition method (SWSM), including reconstruction of expansion coefficients, prediction of acoustic field, error sensitivity analysis, regularization method and a searching method with dual measurement surfaces for determining the optimal number of expansion terms, is established. Subsequently, the spherical wave source boundary point method (SWSBPM) and its application in the NAH are introduced briefly. Considering the similarity of the SWSM and the SWSBPM for realizing the NAH, they are compared. The similarities and differences of the two methods are illuminated by a rigorous mathematical justification and two experiments on a single source and two coherent sources in the semi-free acoustic field. And, the superiority of the NAH based on the SWSBPM is demonstrated.     

试析直接法与间接法的关系

直接法与间接法是企业在现金流量表中列报现金流量的方法.本文阐述了最新会计准则中直接法与间接法的内在关系.直接法和间接法都反映经营活动的现金流量,两种方法都需要把权责发生制转换为收付实现制,两种方法编制现金流量所需的主要依据相同,但两种方法涵义不同,使用条件和特点也不同.     

各向异性固体层的声反射和透射谱

考虑到由流体斜入射到各向异性固体时超声波模式转换的一般情况 ,研究了浸没于流体中一般各向异性固体层状材料的声反射和透射问题 ,给出了问题的求解方法 ,得到各向异性固体层的声反射和透射系数是介质参数 ,固体层取向和厚度 ,入射角及入射波频率的函数 ,给出了从流体斜入射到正交复合材料板上的声反射和透射谱的数值解例     

快速多极边界元方法在二维声散射问题中的应用

快速多极算法(FMM)是求解边界元方法(BEM)在大尺度情况下的一种非常有效的算法.研究了快速多极算法在二维声散射问题的边界积分方程求解中的应用.给出了积分核函数以及其共轭积分算子核函数的多极展开式,局部展开式以及相应展开系数之间的转化关系.分别应用两种不同的层级树结构的FMM来进行求解,并对两种树结构下的求解效率进行了对比.数值算例表明用快速多极算法求解该问题时在存储量和计算量上比直接求解方法效率更高.     

跳水"压水花"动作的计算机模拟--跳水"压水花"技术研究系列论文之一

用计算机模拟方法研究了跳水运动中,运动员入水时,其手型与水花高度之间的关系。研究中把人和水分别简化为不同角度的楔形体和理想流体,建立了满足控制函数和初、边值条件的流、固耦合运动方程,用有限元方法模拟计算了不同斜升角度(从47—80?)楔形体撞入水面的过程。模拟结果显示水花高度与楔形体斜升角度成反比,据此提出水花控制理论和技术,并运用于专业运动员的训练,取得明显的效果。     

立体角的计算与棱角的识别

在边界元方法中,尖棱角处必须做特殊处理以克服此处的法向量的不连续性。通过立体角的值来判断结构中的尖棱角,对实现边界元算法具有重要的作用。现详细地介绍立体角的数值求解方法,并以两个实例说明应用立体角识别尖棱角的可行性。