Chebyshev谱元法求解含吸收边界的二维均匀稳定流场的声传播

从群速度的角度推导了包含均匀稳定来流的二维波动方程的1阶吸收边界条件,基于Che-byshev谱元法提出了二维均匀稳定来流波动方程的求解方法.在空间上采用谱元方法,在时间上采用隐式Newmark积分法,从而获得了波动方程的离散形式.经具体算例验证表明:与1阶Clay-ton-Engquist-Majda吸收边界条件相比,所推导的吸收边界条件能更有效地削弱边界上的数值反射,避免解的失真,求解方法在空间上具有谱精度,在时间上达到了2阶精度.     

二维各向异性位势薄体问题的虚边界元法

传统边界元法分析各向异性薄体问题时涉及奇异边界积分和拟奇异边界积分的处理,估计这些积分具有相当的难度而且耗时.提出了求解二维各向异性位势薄体问题的虚边界元方法,给出了求解此类问题的新途径,同时拓展了虚边界元法的应用范围.数值算例表明,虚边界元法可有效求解二维各向异性位势薄体问题,且方法简单、精度高、易于程序设计.     

FMBEM求解直圆柱线性波绕射问题

运用快速多极子边界元法(Fast Multipole Boundary Element Method,FMBEM)求得单圆柱在线性波浪中的绕射问题的数值解.所谓的快速多极子边界元法就是采用快速多极子法(Fast Multipole Method,FMM)加速传统边界元法的求解速度.在文中通过求解二维的Helmholtz方程证明FMBEM法具有高精度和高效率,适用于求解大规模的数值问题.另外,给出了单圆柱线性平面波绕射问题中相关水动力学系数的数值计算结果.     

轴向流中平行简支弹性薄板的大挠度流固耦合的数值模拟

分析研究了轴向流中简支弹性薄板大挠度流固耦合系统的振动响应和流场特性.板结构动力学方程采用基于位移的有限元法离散;流场采用二维不可压缩粘性流体N-S方程,并用有限体积法离散;在此基础上结合动网格控制技术,建立模拟双向流固耦合作用下轴向流中简支弹性薄板的二维数值模型.利用该数值模型得到了单块简支板随流速变化流致振动特性,研究了结构大挠度的振动稳定性,分别得到了Pitchfork分岔曲线和非线性系统结构的Hopf分叉曲线.通过轴向流恒定流速下不同间距的平行两块简支弹性薄板流固耦合的数值模拟得到了的流致振动特性.     

轴向流中固支弹性薄板的大挠度流固耦合系统的数值模拟

就轴向流中两端固支大挠度弹性薄板的流固耦合振动特性,固支薄板的结构动力学方程用有限元法离散,流场采用不可压缩的二维粘性流体(N-S方程)用有限体积法离散,结合ADINA中的流体单元划分技术,建立了双向流固耦合作用下轴向流中两端固支薄板的二维仿真模型.通过模拟仿真分析研究了给定不同流速下固支板的流固耦合振动特征和大挠度系统的振动稳定性.分别得出了不同流速下固支板中点的挠度—流速曲线、挠度时程曲线及挠曲线图.结果表明:当流速小于固支板的临界流速时,板将处于稳定的直线平衡状态;当流速大于固支板的临界流速时,板将在新的位置达到弯曲平衡状态,以及在弯曲平衡位置附近发生极限环振动.     

用贴体正交曲线网格求解二维复杂几何域对流扩散问题

文章结合有限差分法和边界元法的优点,利用贴体正交曲线网格和有限控制体积积分来分析二维复杂几何域对流扩散问题.叙述了构作贴体正交网格的边界元技术,推导了通用微分方程形式下粘性可压缩层流流动动量方程和紊流流动时均动量方程的源项表达式,给出了曲壁通道内紊流流场的计算实例和实验测量结果.     

边界积分方程中近边界点几乎奇异积分的计算

该文针对边界元法存在近边界点力学量计算的困难,给出了一个通用性方法,将近边界点到边界单元的距离参数通过分部积分变换到积分式之外,从而计算出二维问题近边界点参量的几乎强奇异和超奇异积分.该法同样适用于板壳问题的边界元法,尤其是对于将超奇异边界积分方程正则化为强奇异边界积分方程的边界元法,求解近边界点参量更加有效.     

快速多极子边界元法在三维势流问题中应用

以Rankin源为基本解,采用快速多极子方法加速后边界元法求解由格林第二公式导出的三维势流边界积分方程,进而计算其势场的分布.无限区域中水流绕射算例的数值计算证明,多极子边界元法能给出满意的结果,与传统边界元方法在运算速度和内存消耗上相比具有明显的优势,表明其适合于在现有的计算条件下求解大尺度多未知量势流问题.     

固体力学中快速多极边界元法研究进展

和快速多极方法相结合，使边界元法处理大规模工程与科学问题变得十分有效,首先概括介绍了快速多极边界元法，接着介绍其精度和效率的验证以及与常规边界元法的比较，并给出在微机机群上的快速多极边界元并行算法,给出了快速多极边界元法的一些应用，其中包括：复合材料的二维、三维模拟，含大量裂纹的二维弹性固体及其疲劳裂纹扩展的模拟．此外还介绍了用于弹塑性问题的快速多极边界元新方法．     

涡方法数值模拟高雷诺数圆柱绕流

圆柱绕流问题是二维物体粘性绕流的基本问题.二维N-S方程的各种差分求解格式往往难以用来实际计算高雷诺数流动问题。80年代中期以来,Stansby等发展了70年代初Chorin提出的随机点涡法,用网格涡方法计算点涡的对流速度,大大减小了计算量,笔者采用Stansby等发展了的随机点涡法,对雷诺数2000的圆柱绕流进行了数值模拟。     

深水中复合材料椭圆柱壳声辐射的数值分析

用有限元／边界元法分析深水中复合材料椭圆柱壳的声辐射．首先利用有限元分析壳体的受迫振动,然后提取壳体表面（耦合面）数据,作为声场分析的边界条件,并计及流-固耦合作用,从而可用边界元计算壳体的辐射功率以及壳体表面、近场和远场声压,结果表明系统的耦合振动模态对于其辐射声场具有重要影响．     

圆柱绕流的三维数值模拟

利用计算流体力学软件CFX－4,对粘性不可压缩流体的圆柱绕流进行了三维数值模拟,采用有限体积法和SIMPLE计算程式,利用不可压缩Navier-Stokes方程,模拟雷诺数在亚临界区内的绕流流动,并计算了流体的水动力特性。为克服数值模拟高雷诺数时的数值不稳定性,计算中采用了QUICK迎风格式,其对流项为三阶精度,其余项如扩散项等为二阶精度,圆柱两端边界采用周期性边界条件。计算结果表明,高雷诺数时圆柱周围的流动具有明显的三维特性,且沿柱长方向不同断面的升力和阻力系数并不相同。同时,对圆柱绕流进行了二维数值模拟,并与三维数值结果进行比较,发现三维模拟的升力和阻力系数均小于二维模拟。     

薄板自由振动虚边界元解法

依据虚边界元法思想,提出一种求解薄板自由振动问题的新算法,通过采用薄板弯曲问题的静力基本解建立了薄板自由振动问题的虚边界积分方程,及满足边界条件和域内点动力位移方程,将薄板自由振动问题转变为代数特征值问题,可直接求解,与边界元直接法相比,本方法无需处理奇异积分,避免了“边界层效应”,而且思想简单,计算省时、方便,算例证实了本方法的可行性和计算精度。     

直接边界元法奇异积分的研究

在用直接边界元法解决弹性问题时,当场点与加载点重合时,边界积分方程将出现奇异.为了减少计算误差,有必要求出奇异积分的解析式.应用高等数学基本理论推导出二维弹性问题直接边界元法奇异积分的解析式,采用线性单元离散边界,求出了奇异的对角线子矩阵元素的解析表达式.根据理论推导的结果,编制了相应的计算程序.可用于分析弹性二维问题的位移场和应力场.     

弹性结构与晃动液体耦合系统的动力特性分析

基于弹性波传播理论,构造液体晃动的等效弹性体模型,将液体晃动问题归结为弹性体动力问题, 建立位移-压力格式的流固耦合方程,采用静态缩聚技术求解流固耦合特征值问题.在此基础上,推导地震条件下的流固耦合系统有限元模型.对储液容器的特征值问题和动态响应问题进行了实例分析和计算,算例结果验证了方法的正确性与有效性.     

求解二维非齐次Helmholtz方程的边界元法

基于Laplace方程的基本解讨论了二维非齐次Helmholtz方程的直接边界元解法.通过将Helmholtz方程变形之后加权Laplace方程的基本解和应用Green公式得到相应的直接积分方程,针对积分方程中同时存在域积分项和边界积分项,在应用边界元法分析求解时采用了耦合关于内点和边界点的积分方程求解,最后,通过数值算例验证方法的有效性.     

基于虚拟弹性体的快速动网格方法

为减少流固耦合的计算时间,在现有弹性体方法的基础上,发展了一种快速动网格方法。首先根据弹性体方法的基本假设,将流场网格所包围的空间区域视为虚拟弹性体,然后将结构与该虚拟弹性体视为一个整体系统,并计算其固有振动的振型及频率,最后以结构受到的流体作用力为激励,通过振型叠加法计算结构网格及流场网格节点的位移。考虑到实际结构的流固耦合振动多为低阶模态的振动,在流固耦合计算中可以通过低阶模态的叠加计算流场网格节点的位移,从而达到快速更新流场网格的目的。采用该快速动网格算法,对某弹性梁颤振问题进行了流固耦合分析,计算结果与已有文献的结果吻合很好,说明了该算法的正确性。与现有的弹性体方法相比,该算法使流固耦合计算时间减少了65.5%。对Wing 445.6模型的颤振问题进行了分析,得到颤振边界与实验值吻合良好,且与现有弹性体法相比,可以减少计算时间54.8%。     

弹性底板矩形贮箱流固耦合系统的自由振动分析

基于Bernoulli方程,通过对弹性体与流场耦合作用机理的分析,建立了弹性底板矩形贮箱流固耦合系统的自由振动方程.将系统的自由振动问题降为一维问题来处理,提出了一种分析此系统自由振动问题的半解析方法,并采用迦辽金方法,求解了系统的自由振动频率.最后讨论了弹性底板的抗弯刚度、结构几何参数以及液体的深度与密度对系统自由振动频率的影响,以便更确切地反映弹性底板矩形贮箱与流体耦合作用的实际机理.     

用边界元法求解分片介质中的Helmholtz方程

该文研究了求解分片介质中的Helmholtz方程的边界元法。边界元求解的思路是将分片介质子区域的公共边界当作子区域的外部边界处理 ,在每个子区域采用边界元法 ,再在公共边界上加衔接条件。该文通过大量数值实验 ,并对比边界元法、有限元法、广义差分法求解效果 ,得出边界元法能很好地克服Helmholtz方程解的震荡性 ,采用边界元法求解Helmholtz方程具有稳定性好 ,精度高的优点。     

透平S_1流面不可压缩流场的边界元法

本文介绍用边界元法计算不可压缩回转面流场,并给出了S_1流面上的运动方程,采用线性单元对边界进行离散化,将库塔条件与边界元方程直接联立求解,并且较好地解决了用边界元法解回转面时所遇到的周期性边界问题。