基于重叠边界的有限元与离散元耦合异步长并行计算方法

为了减少有限元和离散元耦合领域问题的计算时间,本文建立了有限元-离散元耦合异步长并行计算方法。整个计算模型分为有限元和离散元两类计算区域,根据每个子域单元特性选择合适的积分时间步长。每个子分区均采用显式时间积分格式进行求解,采用重叠边界法处理边界耦合问题,不涉及子循环过程的插值和截断过程。数值算例表明该方法在保证高精度的同时有效降低了计算时间。

用边界元-有限元法研究高架结构辐射噪声

高架结构是城市轨道交通的重要部分,其引起的结构辐射噪声问题引起了越来越多的关注。基于边界元Helmholtz积分公式和有限元流固耦合系统理论,推导了结构振动加速度与辐射声压之间的声传递向量,该传递向量可通过单元辐射叠加法求得,结合采用直线电机车辆的北京首都机场快轨线某段高架桥梁段现场试验结果进行了比较分析。结果表明:使用该方法计算的声压结果与实测结果的声压频谱分布及变化趋势吻合较好,该方法可用于城市轨道交通高架结构辐射噪声的预测研究。

一种计算轴对称磁场的边界元-有限元混合法的研究

由静磁场的基本方程和有限元的原理出发,研究了一种新的计算带非饱和磁介质轴对称磁场的边界元-有限元混合法,编制了计算电磁聚焦成象系统聚焦磁场的程序。这种方法解决了单纯使用有限元法的边界封闭问题,为成象系统及其它磁透镜的工程设计提供了一种有效方法。

有限元-边界元耦合法在二维电磁场数值计算中的应用及优化措施

在二维电磁场数值分析中，有限元-边界元耦合法将有限元法与边界元法结合，可以同时发挥两种方法的优势．笔者在介绍有限元-边界元耦合法原理与算法的基础上，分别研究子域逼近有限元法、Mortar元法和并行计算与有限元-边界元耦合法结合使用分别为关心区域相对整体区域尺寸较小、求解区域中包含运动导体和大规模高自由度电磁问题的数值计算提供优化途径．

基于比例边界有限元法的结构-地基动力相互作用时域算法的改进

本文基于比例边界有限元法(Scaled Boundary Finite Element Method)的缩减基函数解法,对结构-无限地基动力相互作用的时域算法进行了改进。通过选择适当的基函数数目,缩减结构-地基接触面上的自由度,以减小卷积积分所带来的计算工作量,推导了缩减自由度后运动方程的表达式。通过重力坝和拱坝加速度频响曲线的算例,对比了不同程度缩减的基函数和全部基函数对计算精度和效率的影响。结果表明,基函数的缩减可使计算效率明显提高,但精度损失不大。当采用60%的基函数时,计算效率提高5倍,而峰值频响曲线的精度损失却不超过4%。因此,该算法比较适合于大型结构-地基动力相互作用问题的时域分析。

流-固耦合问题边界元-有限元耦合方法分析

利用边界元 -有限元耦合方法对流 -固耦合振动问题进行了分析 ,假设流体控制方程为 L aplace方程 ,利用非连续边界元对流体域进行离散 ,从而有效地解决了边界元分析中的“角点效应”问题。固体以平面梁为模型 ,采用有限元进行离散 ,对非连续边界元和有限元的耦合问题进行了分析 ,通过对悬臂梁在一侧受液体作用时的瞬态响应分析的数值解同解析解的比较 ,表明了本文所给方法的有效性 ,同时为利用边界元 -有限元耦合方法对流

三维有限元-边界元耦合法的并行计算及优化措施

在应用有限元-边界元耦合法分析大规模三维电磁场数值计算问题时,采用并行计算方式可以有效地节省计算时间和提高求解精度.在介绍并行计算在线性单元和高阶单元问题实施原理的基础上,分别论述了应用自适应交叉逼近技术降低边界元求解区域内存消耗、应用分布式预处理共轭梯度法求解器对复杂且包含易变化部件求解区域进行基于相对编号的建模以提高并行机群的运行效果、应用Borland C＋＋Builder软件实现分析结果导出自动化以降低人工工作量这三种对有限元-边界元耦合法并行计算的优化措施.

有限元-边界元法计算电磁响应(英文)

提出了用有限元边界元法计算三维电磁场感应问题 .首先阐述了棱边单元法的基本性质 ,并且根据基函数的几何平面特性推导出对于任意位置面的表面棱边单元法基本公式 ,然后用棱边单元法离散数学模型 .用此方法模拟地下金属军事目标的电磁响应 ,计算了目标的最低谐振频率及表面电流密度 ,和单一目标的频率特性及位移特性 ,计算结果和实验结果一致 ,因此这种数值方法对计算三维电磁场问题行之有效 .

结构声学耦合随机性分析的等几何有限元-边界元法研究

采用有限元法(FEM)进行水下壳结构振动响应分析,边界元法(BEM)进行结构振动声学分析;组合有限元法与边界元法构成耦合FEM-BEM方法进行水下薄壳结构声振强耦合分析。为了克服传统拉格朗日函数近似几何模型与物理场插值计算时的不连续与低精度问题,采用Loop细分曲面等几何法构建几何模型,并采用相同的样条函数进行物理场高阶插值计算,实现水下声振强耦合系统的CAD/CAE的集成分析。随机性分析致力于研究系统输入的不确定造成的输出不确定。蒙特卡罗模拟(MCs)因简单直接被认为是解决复杂多维不确定性问题的通用工具,然而巨大的计算成本降低了其适用性。采用本征正交分解(POD)和径向基函数(RBF)可降低计算成本,提高计算效率,实现基于MCs的快速随机性分析。考虑结构材料属性参数以及结构形状参数的不确定性对计算结果的影响,采用MCs分析随机变量下的结构声学响应的统计特征。最后通过若干算例验证该算法的正确性与有效性。

航天器结构声振耦合问题的有限元-边界元方法数值仿真研究

针对高超声速、结构轻量化带来的航天器复杂结构声振耦合动力学环境难以精准预示问题,提出一种基于有限元-边界元(FEM-BEM)方法的声振耦合问题数值仿真方法。分析国内外求解声振耦合问题的主要数值仿真方法,从影响航天器性能设计的低中频声振耦合动力学特性入手,采用有限元-边界元方法,以航天器典型结构截锥壳为研究对象,开展混响噪声作用下的声振耦合研究。建立了声振耦合模型,获得了主体结构、主要支撑结构及设备的声振耦合响应。提出的数值仿真方法和仿真分析过程,为航天器声振耦合分析提供了理论基础和仿真手段。

有限元–边界元耦合法在运动导体涡流场中的应用

基于A,φ-A法和库伦规范,推导了导体区域和非导体区域的有限元方程及自由空间的边界元方程,通过引入交界面条件,实现了将边界元矩阵等效为有限元矩阵求解的有限元-边界元耦合法(finite element and boundary element coupling method,FE-BECM)。将FE-BECM应用于TEAM-7问题的计算,验证了该方法处理开域涡流问题的有效性。当FE-BECM应用于运动导体涡流场(moving conductor eddy current,MCEC)问题时,用有限元离散源电流区域和运动部件,用边界元离散自由空间并关联相互独立的有限元区域。该方法克服了常规有限元法使用1套网格处理运动问题所遇到的麻烦。使用有限元-边界元耦合法对单级线圈炮问题进行了计算,验证了算法处理运动导体涡流场问题的有效性。

基于二次场的可控源电磁法三维有限元-无限元数值模拟

本文提出了一种新的混合有限元-无限元三维可控源电磁法(CSEM)问题快速高精度正演模拟算法.首先从电场双旋度方程出发,推导了水平电偶极子源的二次场边值问题,采用无限元代替截断边界条件和有限元离散内部计算区域的新策略,达到减小计算区域的目的,基于并行直接求解技术,实现多源CSEM问题的快速精确求解.其次,通过层状解析模型测试,一方面验证了新算法的正确性,另一方面通过与其他三种已知CSEM问题求解策略进行对比,表明了本文提出的基于二次场有限元-无限元算法具有离散区域小、求解速度快和计算精度高等优点.最后,通过3D模型计算,清晰直观地模拟了场源阴影效应,为野外数据的处理与解释提供指导.

有限元—边界元耦合法在3维开域涡流场中的应用

基于A、φ-A法和库伦规范,推导了导体区域和非导体区域的有限元方程及自由空间的边界元方程,通过引入交界面条件,实现了将边界元矩阵等效为有限元矩阵求解的有限元—边界元耦合法。分别使用单纯有限元和有限元—边界元耦合法对TEAM 7问题进行了计算,通过对比3种有限元模型的计算结果指出了截断边界和远场单元在处理开域问题上的缺陷。为了探讨有限元—边界元耦合法的建模方式,根据包裹空气层的不同建立了3种模型并进行了对比分析。通过对数值仿真结果和实测数据进行对比,证明了有限元—边界元耦合法在处理3维开域涡流场问题上的有效性。

无限元与有限元区域边界位移值的D—N迭代求解法

无限元法是基于有限元方法之上，将无限剖分思想与有限元方法相结合，在应力、形变等规则区域若使用有限剖分计算可达到精度要求时，则使用有限元法计算，而在形变剧烈、有应力集中的区域则采用无限相似单元剖分的无限元法计算，以求得应力集中区域确切应力分布状况，以弥...

地铁振动预测的周期性有限元-边界元耦合模型

针对地铁列车运行引起的隧道结构和自由场中的振动响应问题,提出了一个在频率-波数域内的三维周期性有限元-边界元耦合的数值模型,此模型中隧道结构采用有限元法计算,自由场采用边界元法模拟。此模型采用Floquet变换,利用隧道和自由场在隧道轴线方向上的周期性,把无限长的隧道及自由场的网格划分限制在一个基本元内,这样使动力学数值计算的效率大大提高。利用此模型计算了在隧道底板上施加固定单位谐振荷载情况下隧道-自由场相互作用系统的动力响应,结果表明此模型可应用于地铁列车运行引起的隧道和自由场中的动力响应预测。

基于2.5维有限元-边界元分析轨道随机不平顺影响下的铁路地基振动

建立了完善的2.5维有限元-边界元耦合模型分析移动谐荷载作用下钢轨、扣件、轨枕、道砟、路基和地基等各部分的振动响应;利用既有的车辆-轨道-路基-地基耦合系统垂向振动解析模型得到轨道谐波高低不平顺引起的垂向动态轮轨力;在此基础上,结合轨道随机不平顺功率谱密度,提出了列车运行引起的地基振动功率谱计算方法。对比分析了地基表面测点垂向振动加速度级的理论计算与现场实测结果,证明了该模型的合理性。模型能有效地分析具有复杂横截面形状的轨道-路基-地基系统的振动响应以及多种轨道、地基减振隔振措施的影响,且具有较高的计算效率,适用于铁路线路设计阶段的方案比较研究。

用时域有限元边界元耦合法计算三维瞬态涡流场

文中首次应用时域有限元-边界元耦合法(T-FEBEM) 研究了涡流检测中有限厚平板导体中含有窄裂缝时的瞬态涡流场。采用伽辽金加权余量法,推导出A-φ位函数的时域有限元-边界元耦合算法,讨论了确定时间步长的基本方法, 用不完全分解共轭梯度法(ICCG)求解线性方程组。研究了平板导体中不同深度裂缝的时域响应,并与频域有限元- 边界元耦合法(FEBEM)和付里叶变换方法所得结果进行了对比。计算结果表明,时域有限元-边界元耦合法效率更高,是一种求解瞬态开域涡流场的有效方法。

有限元与边界元耦合法在船舶结构直接设计分析系统中的应用

讨论了船舶结构分析中的结构模型化问题，提出一种适用于船体舭部、肘板等特殊部位的有限元与边界元耦合模型，解决单纯用有限元法难于处理的问题。此方法已用于船舶结构直接设计分析系统（DDASSV1.0），在初步设计中可一次获得较满意的纵横强度，多数情况下不需要进行再分析。因此，本方法可大大减轻前后处理的工作量，并能保证良好的精度。

一类各项异性半线性椭圆方程自然边界元与有限元耦合法

将冯康和余德浩提出的自然边界归化方法用于研究一类半线性椭圆方程外区域问题,提出一种自然边界元与有限元的耦合算法.针对某一类半线性椭圆方程,应用变分原理,研究其弱解性及Galerkin逼近,得到有限元解的误差估计及收敛阶O(hn),最后给出相应数值例子.

土、结构相互作用问题的边界元、有限元耦合法

本文简要介绍了土、结构相互作用边界元、有限元耦合法的发展；着重对近年来应用边界元、有限元耦合法解动力相互作用问题的代表性的研究工作做了介绍，对边界元、有限元耦合法的几个问题的进展进行了综合评述，最后提出了作者对几个问题的见解。