薄体结构温度场的高阶边界元分析

分析了二维问题边界元法3节点二次单元的几何特征,区分和定义了源点相对高阶单元的Ⅰ型和Ⅱ型接近度.针对二维位势问题高阶边界元中奇异积分核,构造出具有相同Ⅱ型几乎奇异性的近似核函数,在几乎奇异积分单元上分离出积分核中主导的奇异函数部分.原积分核扣除其近似核函数后消除几乎奇异性,成为正则积分核函数,并采用常规Gauss数值方法计算该正则积分;对奇异核函数的积分推导出解析公式,从而建立了一种新的边界元法高阶单元几乎奇异积分半解析算法.应用该算法计算了二维薄体结构温度场算例,计算结果表明高阶单元半解析算法能充分发挥边界元法优势,显著提高计算精度.     

瞬态热传导的奇异边界法及其MATLAB实现

基于动力学问题时间依赖基本解的奇异边界法是一种无网格边界配点法.该方法引入源点强度因子的概念从而避免了基本解的源点奇异性,具有数学简单、编程容易、精度高等优点.将该方法用于瞬态热传导问题的数值模拟,运用MATLAB实现该问题的数值研究,并创建相应的MATLAB工具箱.针对二维和三维瞬态热传导问题,进行了基于反插值技术和经验公式的奇异边界法MATLAB算例实现.针对支撑圆坯低温瞬态温度场的模拟结果表明,瞬态热传导奇异边界法的MATLAB工具箱具有简单、方便、精确可靠的优点.研究成果有助于发展瞬态热传导的奇异边界法,并为瞬态热传导问题的数值分析和仿真提供了一种简单高效的模拟工具.     

A Novel Localized Meshless Collocation Method for Numerical Simulation of Flume Dynamic Characteristics北大核心CSCD

局部边界节点法是一种基于非奇异半解析基函数和移动最小二乘原理的新型无网格配点技术,该方法把每个节点处的未知变量表示为该点对应的局部子域内节点处物理量的线性组合,该文基于局部边界节点法对数值波浪水槽进行了研究。首先,通过基准算例确定了Laplace算子非奇异半解析基函数的合理形状参数值。进一步,基于合理的参数选取,用较少的离散节点即可成功模拟波浪传播行为,将得到的数值结果与其他文献数值结果比较,可以发现局部边界节点法用更少的局部点即可得到较好的数值结果。最后,以保护近海岸建筑物为目标,模拟了水下防波堤对波浪传播的影响。结果表明,当波浪与梯形防波堤发生作用后,波峰变得比较陡峭,而波谷变得相对比较平坦,为近海岸防波堤的相关研究和设计提供了数值参考。     

利用边界元法求解一类重调和方程

边界元法(BEM)和多重互易法(MRM)相结合求解一类重调和方程.通过重调和基本解序列给出的MRM-方法和BEM, 推导出该类问题的MRM-边界变分方程, 用边界元法求解该变分方程, 从而得到重调和方程的近似解, 并给出了解的存在唯一性证明.通过数值算例说明了MRM-方法具有收敛速度快、计算精度高, 易编程等优点, 为使用边界元法数值求解重调和方程提供了方法和理论依据.适合于工程中的实际运算.     

板弯曲断裂计算的边界配置解法

采用复变函数理论和边界配置方法,分析计算了Kirchhoff板的弯曲断裂问题.假设了位移及内力的复变函数式,它们能满足一系列的基本方程和支配条件,例如域内的平衡方程、裂纹表面的边界条件、裂纹尖端的应力奇异性质.这样,仅板边界的边界条件需要考虑.它们可用边界配置法和最小二乘法近似满足.对不同边界条件和载荷情形进行了分析计算.数值算例表明,本文方法精度较高,计算量小,是一种有效的半解析、半数值计算方法.     

一种改进边界元法在弹性扭转问题中的应用

本文用一种改进边界元法分析与计算了椭圆截面等直杆的扭转问题.并与边界元法的解进行比较,其结果极为符合.然而,改进边界元法较边界元法所需要的数据量少得多,计算时间也将大大减少了.因此,本文方法对求解Poisson方程问题是一种经济而行之有效的数值计算方法.     

多裂纹问题计算分析的本征COD边界积分方程方法

针对多裂纹问题，若采用常规的数值求解技术，计算效率较低.为实现多裂纹问题的大规模数值模拟，建立了本征裂纹张开位移（crack opening displacement, COD）边界积分方程及其迭代算法，并引入Eshelby矩阵的定义，将多裂纹分为近场裂纹和远场裂纹来处理裂纹间的相互影响.以采用常单元作为离散单元的快速多极边界元法为参照，对提出的计算模型和迭代算法进行了数值验证.结果表明，本征COD边界积分方程方法在处理多裂纹问题时取得较大的改进，其计算效率显著高于传统的边界元法和快速多极边界元法.     

平面位势问题中一种非奇异边界积分法

确立平面位势的边界元直接法中边界积分的解析计算框架系统 ,从而避免了传统的高斯近似求积分 .数值算例表明它具有较高的精度和效率 .特别是在边界量和边界附近区域内点物理量的计算可获得较高的精度 .     

平面定常Stokes问题的无奇异第一类边界积分方程

对无奇异边界积分方程归化法的研究,已有的结果都是针对直接变量的,其核心思想是利用刚体位移(包括刚体的转动和平移)或均匀场．然而,对第一类边界积分方程的无奇异边界归化法的研究,至今还未涉足．本文提交一种新方法,归化出平面定常Stokes问题的第一类无奇异边界积分方程,并建立完整的数值求解体系．一个简单的算例表明本文方法可获得理想的数值结果,特别是边界量的数值结果。     

直接边界元法中边界积分的解析处理

确立了平面位势和弹性力学问题的边界元直接法中边界积分的解析计算框架系统,从而避免了传统的主似求积分,数值算例表胆它具有较高的精度和效率,特别是在边界量和边界附近区域内点物理量的计算可获得较高的精度。     

用时变雷诺方程模型模拟孤立波与半圆型防波堤的相互作用

以时变雷诺方程为控制方程,用k-ε模型封闭该方程,采用体积函数(VOF)方法来跟踪波动自由表面,建立了二维垂向波浪数学模型,并用已有的实验资料进行了验证．随后用该模型模拟了半圆型防波堤与孤立波在淹没、平顶水位、完全露顶且不越浪3种典型工况下的相互作用过程．得到了半圆堤附近的流场、压强场和波面的变形过程．结果表明,在淹没状态下,半圆堤背浪面的底部会产生涡旋;平顶水位时,由于越浪的冲击作用,在半圆堤的背浪面将逐渐形成一对较大的涡旋,而半圆堤背浪面的底部,速度始终相对较小;而在露顶不越浪时,半圆堤的迎浪面会出现波浪的二次爬升的现象．为进一步研究结构物附近的污染物的输移扩散和泥沙运动提供基础．     

Self-adaptive FEM numerical modeling of the mild-slope equation

Based on the linear wave theory, the mild-slope equation (MSE) is a preferred mathematical model to simulate nearshore wave propagation. A numerical model to solve the MSE is developed here on the basis of a self-adaptive finite element model (FEM) combined with an iterative method to determine the wave direction angle to the boundary and thus to improve the treatment of the boundary conditions. The numerical resolution of the waves into ideal domains and multidirectional waves through a breakwater gap shows that the numerical model developed here is effective in representing wave absorption at the absorbing boundaries and can be used to simulate multidirectional wave propagation. Finally, the simulated wave distribution in a real harbor shows that the numerical model can be used for engineering practice.     

基于非局部理论的黏弹性纳米杆轴向振动与波传播研究

根据非局部理论和Kelvin黏弹性理论，针对黏弹性纳米杆自由振动和波传播的轴向动力学问题进行研究.首先，推导了黏弹性纳米杆的轴向动力学微分控制方程.然后，通过无量纲化讨论了3种典型边界纳米杆的前三阶振动特性.最后，研究黏弹性纳米杆波的传播问题，导出了圆频率、波速与波数之间的关系.数值结果表明，非局部效应使第一、二阶固有频率持续减小，第三阶频率先增大再减小，出现结构刚度削弱和增强两种趋势.特别地，对于自由端存在集中质量的情形，第二阶频率随着黏性系数增大出现了多值情况，易导致杆件失稳.数值算例还说明了非局部效应的增强可有效降低黏性材料的阻尼效应，产生逃逸频率，使得纵波能够在高波数段传播.另外，黏性系数在低波数段对阻尼比影响可忽略不计，而在高波数段下，黏性系数越大则阻尼比越大.     

Generating Boundary Conditions for the Calculation of Tsunami Propagation on Nested Grids

Some boundary conditions used to numerically simulate tsunami generation and propagation are studied. Special attention is given to generating boundary conditions thatmake it possible to simulate tsunami waves with desired characteristics (amplitude, time period and, in general, waveform). Since the water flow velocity in a propagating tsunami wave is uniquely defined by its height and ocean depth, one can simulate a wave propagating from the boundary into the simulation area. This can be done by specifying the wave height and water flow velocity on the boundary. This method is used to numerically simulate the propagation of a tsunami from the source to the coast on a sequence of refined grids. In this numerical experiment the wave parameters are transferred from the larger area to the subarea via boundary conditions. This method can also generate a wave that has certain characteristics on a specified line.     

Singular boundary method based on time‐dependent fundamental solutions for active noise control

Active noise control is an efficient strategy of noise control. A numerical wave shielding model to inhibit wave propagation, which can be considered as an extension of traditional active noise control, is established using the singular boundary method using time‐dependent fundamental solutions in this study. Two empirical formulas to evaluate the origin intensity factors with Dirichlet and Neumann boundary conditions are derived respectively. In comparison with other similar numerical methods, the method can obtain highly accurate results using very few boundary nodes and small CPU time. These meet the major technical requirements of simulation of active noise control. The subsequent numerical experiments show that the proposed model can shield efficiently from the wave propagation for both inner and exterior problems. By applying the newly derived empirical formulas, the CPU time of the singular boundary method is further reduced significantly, which makes the method a competitive new and efficient meshless method. In addition, the singular boundary method makes active noise control in an online manner via time‐dependent fundamental solutions as its basis functions.     

Dynamic contact between a spherical inclusion and a matrix upon incidence of an elastic wave

Using the boundary element method, we have studied the dynamic displacements and stresses in an infinite elastic matrix with a spherical elastic inclusion, caused by the propagation of an elastic wave. The original problem has been reduced to a system of boundary integral equations for the contact displacements and tractions on the interface between the inclusion and matrix. Based on the numerical solution of these equations, we have analyzed the influence of the direction of wave propagation and frequency on the important physical parameters, depending on the elastic characteristics of composite constituents.     

Numerical study of the stability of breakwater built on a sloped porous seabed under tsunami loading

Strong earthquake induced huge tsunami has occurred for three times in Pacific ocean in recent ten years; for example, the tsunami triggered by the Sumatra earthquake in 2004, Chile earthquake in 2010 and Tohoku earthquake (Japan) in 2011. Tsunami carrying huge energy always would bring high risks to the population living near to coastline. Breakwater is widely used to dissipate the wave energy, and protect coastline and ports. However, they are vulnerable when being attacked by tsunami wave. At present, the interaction mechanism between tsunami, breakwater and its seabed foundation is not fully understood. In this study, the dynamics and stability of a breakwater under the attacking of tsunami wave is investigated by adopting an integrated model PORO-WSSI 2D, in which the VARANS equation for wave motion, and the Biot’s dynamic equation for soil are used. Based on the numerical results, it is found that offshore breakwater interacts intensively with tsunami wave when it overtopping and overflowing over a breakwater. The impact force on the lateral side of breakwater applied by tsunami wave is huge. The shear failure is likely to occur in the seabed foundation of breakwater. The liquefaction is unlikely to occur due to the fact that there is basically no upward seepage force in seabed foundation in the process of tsunami wave passing through the breakwater.     

入水冲击问题变分原理及其它

首先建立入水前后两个衔接阶段的较为严密的场方程.再得到与之对应的各类变分原理,界限定理,第二阶段问题的边界积分方程.证明了解的存在性并提供了求解实施方案.最后以船舶兴波阻力问题的算例,论证了第二阶段问题的一种特殊应用及其正确性.从而为求取较为精确的入水冲击问题基本方程的变分有限元及边界元方法奠定了严密的理论基础.     

Energetic BEM-FEM coupling for the numerical solution of the damped wave equation

Time-dependent problems modeled by hyperbolic partial differential equations can be reformulated in terms of boundary integral equations and solved via the boundary element method. In this context, the analysis of damping phenomena that occur in many physics and engineering problems is a novelty. Starting from a recently developed energetic space-time weak formulation for the coupling of boundary integral equations and hyperbolic partial differential equations related to wave propagation problems, we consider here an extension for the damped wave equation in layered media. A coupling algorithm is presented, which allows a flexible use of finite element method and boundary element method as local discretization techniques. Stability and convergence, proved by energy arguments, are crucial in guaranteeing accurate solutions for simulations on large time intervals. Several numerical benchmarks, whose numerical results confirm theoretical ones, are illustrated and discussed.