非线性节块展开法

非线性节块展开法是一种快速计算反应堆物理参数的节块方法。与有限差分法相比，它具有速度快、精度高的特点；与另一种节块法—Ｇｒｅｅｎ函数展开法相比，它又保留了有限差分法的优点：不但可以得到有效倍增因子、功率分布等参数，还可以得到差分方程的系数矩阵，因而十分适合与广义微扰理论（ＧＰＴ）相结合以获得更快的计算速度，应用于堆芯燃料管理优化中。通过两个二维基准问题的检验计算，比较了非线性节块展开法和Ｇｒｅｅｎ函数展开法程序的速度与精度，结果表明两个程序基本相当，非线性节块展开法可以满足实际应用的需要。     

一种节约内存的局部一维减缩时域有限差分方法

为了提高电大尺寸目标电磁计算效率,提出了一种用于克服时域有限差分算法稳定条件限制、降低计算所需内存的新方法。证明了即使在无源区域,局部一维时域有限差分法所给出的电磁场量不满足零散度关系,并推导了该散度关系的具体表达式。基于该非零散度关系和麦克斯韦旋度方程,将三维局部一维时域有限差分法与减缩时域有限差分法相结合,从而得到三维的局部一维减缩时域有限差分法。通过仿真计算证明该方法与LOD-FDTD方法的计算结果一致,具有良好的计算稳定性。     

用有限差分法分析列车对隧道中电磁波传播特性的影响

文章推导出复杂边界条件下二维Helmholtz方程的差分公式，实例分析了列车对矩形隧道和半圆形隧道中电磁波截止频率的影响、列车-矩形隧道中主模的衰减常数。实例计算结果表明：有限差分法的计算结果与有限元法计算结果完全一致，有限差分法是分析列车-隧道中电磁波传播问题的一种有效方法。     

水泥混凝土路面的有限差分法

在修改《刚性路面》设计规范中,为了验算有限单元法解,本文利用有限差分法算出水泥混凝土路面在车轮荷载作用下的应力和变形.与有限单元法比较,有限差分法解具有同样的精度;所用数学力学模型及计算方法较为简单,易于掌握.有限差分法也是计算水泥混凝土路面的可行方法.     

有限元-有限差分法二维波动逆时偏移初探

提出了采用有限元-有限差分实现二维波动方程的逆时偏移算法。该方法在空间上，联合采用有限元法和有限差分法；对于地表(水平)方向，使用有限元法进行离散，将原方程转化为一个一维(深度和时间)问题的方程组；在深度和时间方向上，采用有限差分法来求解。介绍了算法的基本原理，给出了计算实例并与使用F-K(频率-波数)域相移法、频率-空间域有限差分法的结果进行了比较。与采用有限元的偏移方法相比，本方法可以节省大量内存；与采用有限差分的偏移方法相比，可以在一定程度上提高计算精度。本算法有可能在地震勘探数据处理中发挥一定的作用。     

运用局部有限差分法确定变截面梁最大挠度的位置

在工程实际中,通常使用有限差分法来确定变截面梁最大挠度及其所在位置,然而这种方法十分繁杂.提出了局部有限差分法,采用逐次减小差分区域和步长的方法,使计算精度层层递进,计算过程简单易行.     

波导截止特性和衰减常数的有限差分分析

文章提出用有限差分法计算波导截止频率、截止波数、衰减常数,推导出用有限差分法计算相关参数的计算公式,编程计算了矩行波导、L型波导和单脊型波导的TE、TM波的截止频率、截止波数及衰减常数,计算结果与解析解非常吻合,证明文章提出的计算方法是正确和可行的。     

径向薛定谔方程的有限差分解法

本文应用有限差分法对谐振子势与非谐振子势进行了计算。结果表明，该方法具有广泛的应用前景。     

格林函数节块法求解中子伴随通量方程

将中子伴随通量方程变换处理以后，利用格林函数节块法的原理求解中子伴随通量，以此代替过去采用最多的细网有限差分法。对ＴＷＩＧＬ二维模型进行了校算，并与有限差分法ＣＩＴＡＴＩＯＮ程序结果进行了比较。两种方法的计算结果十分吻合，这表明采用粗网格的格林函数节块法求解中子伴随通量是可行的，而且这种粗网格的方法具有空间网格大，计算精度高，计算速度快的特点，因而具有较高的计算效率。     

线性微分方程的有限差分法求解周期正解

为了提高线性微分方程在求解计算中的性能，提出了线性微分方程的有限差分法求解周期正解。证明了线性微分方程的有限差分格式与能量守恒的离散规则是等价关系，分析了线性微分方程差分格式的正性和周期性，实现了基于有限差分法的线性微分方程求解周期正解。最后通过数值实验的方式，得出与传统差分法求解线性微分方程的周期正解相比，不仅满足了正性要求，还满足了周期性要求。     

异形网格差分在三维场计算中的应用

论述了一种用于三维场计算的异形同格差分方法.此法不但具有差分法要求计算机内存小、计算速度快等优点,而且具有有限元法的单元剖分灵活、适宜任意形状边界等优点,并可自动生成网格.在已有的二维数学模型的基础上,作者开发了三维数值计算差分公式,以四极透镜为例编制程序并进行了计算,计算结果表明:理论值与计算值相符,最大相对误差不超过10~(-3).     

基于光滑有限元的含裂纹复合材料的虚拟裂纹闭合法

为提高求解断裂参数的精度和效率,将光滑有限元法与虚拟裂纹闭合法相结合,提出光滑有限元-虚拟裂纹闭合法.对含不同长度和角度的倾斜裂纹复合材料圆板的断裂参数进行了求解,并与有限元-虚拟裂纹闭合法计算结果进行了对比.数值算例结果验证了该方法具有高精度,裂纹尖端处单元不需特殊处理,对网格尺寸要求低等优点,是分析断裂问题简洁高效的数值计算方法.     

基于自适应移动网格的Cahn-Hilliard方程的有限元法

针对二维Cahn-Hilliard方程,使用自适应移动网格,建立有限元数值模型。由于Cahn-Hilliard方程在初期变换迅速,且在后期变化缓慢,使用基于移动网格偏微分方程(moving mesh partial differential equation,MMPDE)的移动网格准则能够更好地捕捉相变的过程。在移动网格上,对空间方向使用线性有限元离散,对时间方向使用五阶Radau IIA格式离散。数值结果表明在移动网格下的数值解能够很好地保持原方程固有的质量守恒与能量稳定定律,提高计算效率,验证了该方法的有效性和可行性。     

无网格伽辽金法在施工导流二维数值模拟中的应用

针对有限元法等传统数值计算方法存在受单元网格限制、前后处理工作复杂的问题,提出应用一种数值计算方法--无网格伽辽金法,对具有简单边界条件的水利水电工程施工导流的恒定二维浅水流动问题进行了分析、计算.同时利用有限元法进行了对比计算,从流速、水位的计算结果来看,两种计算方法结果相近、误差较小,表明采用无网格伽辽金法解决此类问题是可行的.     

典型大型商用飞机撞击问题数值建模及计算

基于对空客A320的结构分析,开展了大型商用飞机结构简化和建模方法研究,建立了结构较为详细的飞机有限元计算模型,确定了适当的网格尺寸;通过将飞机撞击刚性靶体的数值模拟结果与Riera理论计算结果进行对比,验证了所建立的大型商用飞机有限元计算模型的合理性以及数值计算方法的可靠性.该大型商用飞机有限元计算模型可用于飞机冲击载荷特性的进一步分析,也可直接用于飞机撞击作用下结构响应安全性评价的数值模拟.     

基于BIM+数字孪生的高边坡施工检算动态反馈方法

为有效进行高边坡施工过程中不利工况的施工检算及动态反馈,以二次开发为依托,提出了基于参数化BIM模型及有限元方法快速进行高边坡不同工况稳定性检算的整套方案,并通过实际工程案例进行验证。首先,采用Revit二次开发技术实现不利工况高精度正向反馈至BIM模型;其次,将该BIM模型导入Abaqus中进行快速检算,基于Abaqus二次开发将有限元计算结果的自定义结构化输出;最后,采用Mesh重构方法将快速检算的结果向BIM环境反馈。结果表明：所提方案实现了数字孪生体与施工现场之间的同步关联,并且实现了数值计算结果与BIM环境互通。可见,BIM技术通过开发后可真实高效反映施工情况,并作为中间枢纽实现施工检算与成果管理,可为相关施工检算及BIM应用提供新思路。     

Finite volume method for investigating anisotropic conductivity in EEG

A novel finite volume method is presented to investigate the effect of anisotropic conductivity on the potential distribution of the scalp. The second-order interpolation of the tetrahedral mesh is used to avoid employing the secondary element. The calculation precision is enhanced by using the Gaussian integration. To avoid the geometric singularity, the triangular prism is employed in place of the conventional hexahedron mesh. With the method, the spherical models as well as the realistic head models are simulated. The calculation results indicate that the anisotropic ratio and the position of dipole sources have great influence on the potential distribution in the electroencephalogram.     

斜齿轮时变啮合刚度算法修正及影响因素研究

针对势能法计算斜齿轮时变啮合刚度精度不足问题,提出一种刚度修正算法.考虑端面重合度大于或小于轴向重合度两种情况下单齿接触线长度的不同表达形式,建立齿根圆与基圆不重合时的变截面悬臂梁模型,采用切片法和积分思想推导并计算了斜齿轮啮合刚度,通过与ISO算法和有限元法对比分析,验证了该修正算法的可行性.在此基础上,探讨了螺旋角、模数、齿数、齿宽和压力角等参数对啮合刚度的影响.计算与分析表明,啮入段的相对时间与端面重合度和轴向重合度大小及比重有关;齿轮基本参数的变化引起重合度和单齿啮合刚度的改变,进而影响综合啮合刚度波动值和均值;当端面重合度或轴向重合度在整数附近时,啮合刚度波动值较小,而总重合度在整数附近时,啮合刚度波动值较大.与传统势能法相比,修正算法提高了斜齿轮时变啮合刚度的计算精度,在斜齿轮刚度激励的准确计算方面具有较强的实用性.     

Finite volume method for investigating anisotrooic conductivitv in EEG

A novel finite volume method is presented to investigate the effect of anisotropic conductivity on the potential distribution of the scalp. The second-order interpolation of the tetrahedral mesh is used to avoid employing the secondary element. The calculation precision is enhanced by using the Gaussian integration. To avoid the geometric singularity, the triangular prism is employed in place of the conventional hexahedron mesh. With the method, the spherical models as well as the realistic head models are simulated. The calculation results indicate that the anisotropic ratio and the position of dipole sources have great influence on the potential distribution in the electroencephalogram.     

用高阶有限元素法模拟双侧向测井仪器在层状介质中的响应

利用有限元法模拟双侧向测井响应,基函数的特性对有限元计算效率和精度影响较大.当基函数为线性函数时,要达到理想的仿真精度需要采用十分稠密的网格剖分,这将导致计算量增加从而影响仿真效率.采用高阶基函数插值,模拟双侧向测井响应,结果表明,稀疏的网格剖分在采用高阶基函数时仍然可以保证较好的仿真精度,同比线性插值基函数,高阶基函数具有高效率和高精度特点.针对深侧向截断误差较大的问题,通过优化网格划分和边界处理方法可提高深侧向结果精度.