瞬态热传导问题的加权最小二乘无网格法

加权最小二乘无网格法是一种基于节点信息的纯无网格法,该方法使用最小二乘法建立系统的变分原理,通过移动最小二乘法构造近似函数,控制方程在节点处的残量使用最小二乘法予以消除,边界条件通过罚函数法引入。本文推导了瞬态热传导问题的加权最小二乘无网格计算格式,编制了相应的计算程序,算例结果表明,该方法具有精度高、前后处理简单的优点,是一种高效的的新型无网格法。     

偶应力理论的无网格法

根据偶应力理论的基本方程,在能量泛函中引入罚函数来满足偶应力理论的几何约束条件,从而得到偶应力理论的无网格法的控制方程.通过算例,分析了偶应力对带中心圆孔的无限平板在单轴拉伸及纯剪状态下的应力集中的影响,验证了该无网格方法的有效性和可行性.     

基于局部Petrov-Galerkin离散方案的无网格法

基于局部Petrov-Galerkin离散方案,选用自然邻近插值构造试函数,用Shepard函数作为权函数,提出了一种无网格方法(MNNPG),这种方法充分发挥了局部Petrov-Galerkin法的优势,并且结合了自然邻近插值的特点,方便引入边界条件,由于以Shepard函数的圆形支集作为积分子域,用分片中点插值来完成区域积分,无需额外背景网格,是一种真正的无网格法。本文将该无网格方法用于求解二维弹性力学边值问题,算例结果很好地吻合了精确解,表明该方法具有良好的数值精度和稳定性。     

The finite element method based on interpolating with wavelet basis function

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｐｒｏｂｌｅｍｓｗｉｔｈｌａｒｇｅｇｒａｄｉｅｎｔａｒｅｃｏｍｍｏｎｉｎｐｒａｃｔｉｃａｌｅｎｇｉｎｅｅｒｉｎｇｆｉｅｌｄｓ,ｅ．ｇ．ｉｎｍａｔｅｒｉａｌｌｏｃａｌｉｚａｔｉｏｎ,ｗｉｔｈｉｎｔｈｅｌｏｃａｌｉｚａｔｉｏ...     

Meshless method based on collocation with consistent compactly supported radial basis functions

Based on our previously study, the accuracy of derivatives of interpolating functions are usually very poor near the boundary of domain when Compactly Supported Radial Basis Functions (CSRBFs) are used, so that it could result in significant error in solving partial differential equations with Neumann boundary conditions. To overcome this drawback, the Consistent Compactly Supported Radial Basis Functions (CCSRBFs) are developed, which satisfy the predetermined consistency conditions. Meshless method based on point collocation with CCSRBFs is developed for solving partial differential equations. Numerical studies show that the proposed method improves the accuracy of approximation significantly. The project supported by the National Natural Science Foundation of China (10172052)     

A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.     

基于单位分解法的无网格数值流形方法

在数值流形方法和单位分解法的基础上，提出了无网格数值流形方法. 无网格数值流形 方法在分析时采用了双重覆盖系统，即数学覆盖和物理覆盖. 数学覆盖提供的节点形成求解 域的有限覆盖和单位分解函数；而物理覆盖描述问题的几何区域及其域内不连续性. 与原有 的数值流形方法相比，无网格数值流形方法的数学覆盖形状更加灵活，可以用一系列节点的 影响域来建立数学覆盖和单位分解函数，具有无网格方法的特性，从而摆脱了传统的数值流 形方法中网格所带来的困难. 与无网格方法相比，由于采用了有限覆盖技术，试函数的构造 不受域内不连续的影响，克服了原有的无网格方法在处理不连续问题时所遇到的困难. 详细推导了无网格数值流形方法的试函数和求解方程，最后给出了算例，验证了该方法的正 确性.     

轴对称热冲击SIF过渡过程分析的热权函数法

热权函数法利用温度与热权函数的乘积的积分来直接计算热冲击过程中裂纹尖端的应力强度因子过渡过程。热权函数与时间τ无关,由于免除了对每一时刻τ所需作的有限元或边界元分析,计算过程大大简化,计算效率得到极大提高,本文将热权函数与有限元法直接耦合,给出了基于刚度阵导数法的轴对称问题的热权函数计算格式,实例计算表明,本文给出的热权函数计算格式具有满意的计算精度。     

Application of multi-dimensional wavelet transform to fluid mechanics

This paper first reviews the application research works of wavelet transform on the fluid mechanics. Then the theories of continuous wavelet transform and multi-dimensional orthogonal(discrete) wavelet transform, including wavelet multiresolution analysis, are introduced. At last the applications of wavelet transform on 2 D and 3 D turbulent wakes and turbulent boundary layer flows are described based on the hot-wire, 2 D particle image velocimetry(PIV) and 3 D tomographic PIV.     

三维热权函数法和多虚拟裂纹扩展技术

热权函数法直接利用温度场与热权函数的乘积的积分来求应力强度因子(SIF)的过渡过程，它可以免除对每一时刻进行热弹性力学有限元或边界元应力分析，计算效率大大提高。本文给出了三维热权函数法的基本方程，并提出了求解三维热权函数法基本方程的多虚拟裂纹扩展法(MVCE法)。在MVCE法中，可以引入无穷多个虚拟裂纹扩展模式；虚拟裂纹扩展模式与应力强度因子的插值直接相联系，所得到的方程组的系数矩阵是一个三对角矩阵，具有良好的计算性能。它对于裂纹前缘SIF分布急剧变化的情况，有良好的数值模拟能力。实例计算表明，MVCE法具有权高的计算效率，并具有很高的计算精度。     

The wavelet analysis method of stationary random processes

I.IntroductionRalldomvibrationoftenoccursinpracticalellgilleeringproblclllsl'l.Forexalllple,theunevellroadexcitesthewheelwithuncertainnZagnitudeandirregulartimelapseandresultsinrandonlvibrationofvehicle,whcllsufferingfromearthquakeor11urricallc,tilestructurewillvibraterandomly.Thewaytostudyrandomvibrationisusillgmathematicaltool,randomprocessestheory,todcalwithvibrationproblem,soitisveryilnporlant1'L,r1.alldomvibrationstudyingtostudytilel.andolllprocesses.Genel'ally.thestudyofrandomprocesses…     

基于无网格Galerkin法求解矩形绝缘管道内的磁流体流动

无量纲磁流体流动控制方程中的哈特曼数较大将导致数值计算发散或误差过大。将无网格Galerkin法引入绝缘管道内的稳定磁流体流动计算中,针对磁流体控制方程中大哈特曼数导致计算误差大的情况,对无网格Galerkin法添加了稳定项。计算结果表明,同等条件下,添加了稳定项的无网格Galerkin法总体相对误差远小于标准无网格Galerkin法的结果,且可以计算哈特曼数最大达50000绝缘管道内的磁流体流动。     

无网格法在三维复杂超音速流场模拟中的应用研究

对最小二乘无网格方法在含复杂外形三维超音速流场中的应用进行了研究.选用分解法求解采用最小二乘法得到的对称方程组,针对最小二乘无网格方法的计算特点生成近似正交均匀分布的离散点,对B1AC2R常规导弹超音速流场采用最小二乘无网格方法进行了无粘数值模拟,计算了B1AC2R常规导弹在不同攻角下的轴向力、法向力及俯仰力矩系数,并将数值结果与实验结果进行了比较.结果表明,最小二乘无网格方法在求解含复杂外形超音速流场时具有较高的准确度,将其应用于三维含复杂外形超音速流场的模拟是完全可行的.     

Multiresolution reproducing kernel particle method for computational fluid dynamics

Multiresolution analysis based on the reproducing kernel particle method (RKPM) is developed for computational fluid dynamics. An algorithm incorporating multiple-scale adaptive refinement is introduced. The concept of using a wavelet solution as an error indicator is also presented. A few representative numerical examples are solved to illustrate the performance of this new meshless method. Results show that the RKPM is a good candidate for tackling the widespread large-scale problems in fluid dynamics. © 1997 John Wiley & Sons, Ltd.     

基于小波变换的爆破地震信号能量分析法的应用研究

根据爆破地震信号具有持时短、突变快等特点,结合工程爆破地震监测资料,利用小波包良好的时频局部化性质对爆破地震信号的能量分布特征进行了分析,研究了爆破地震信号的能量在传播过程中的变化规律,得到了爆破地震信号不同频带上的能量分布。根据爆破地震信号不同频带的特征频率与受控建筑物自振频率之间的关系,确定爆破地震对建筑物的影响。并用工程实例说明了该方法比用单一强度因子作为爆破地震的安全判据更加有效。     

影响无网格方法求解精度的因素分析

基于移动最小二乘法的无网格方法的计算精度除受到节点的分布密度和基底函数的阶次影响外，还受到其它因素的影响，其中权函数的选取、权函数影响域的大小及位移边界条件的引入对计算精度影响较大。本文分析了几种常用权函数在数值计算时的特点，包括计算精度、收敛情况、计算效率等，同时分析了影响域大小及边界条件的引入对计算精度的影响。通过分析给出了确定权函数及其影响域大小的方法。当受约束的自由度较多时，通过配点法引入位移边界条件会引起计算结果的振荡，通过施加稳定项可以消除振荡现象，通过对带孔方板的受力分析证明了其可行性。应用以上结论对J23—10曲柄压力机机身进行了受力分析，应力集中部位的计算结果得到了较高的精度。     

MLPG混合配点法在形状优化中的应用研究

建立了无网格MLPG(Meshless Local Petrov-Galerkin)混合配点法求解二维弹性体位移、应力的数学模型,使用罚函数法添加本质边界条件,并将其应用到结构形状优化,结合遗传算法提出了一种新的连续体结构优化设计方法.对于节点支持域半径的选取进行了重点探讨,提出一种动态支持域选择方法,建立了基于MLP...     

基于小波分析法的光纤陀螺系统标定方法

光纤陀螺的随机噪声很大程度地影响系统的标定精度,通常的解决方法是增加平均次数、延长采样时间等,不能从根本上解决问题。利用小波变换阈值滤波对光纤陀螺系统标定的测试数据进行预先处理,采用Allan方差法分析滤波前后的测试数据,发现滤波后陀螺各项随机误差系数均得到明显抑制。用消噪后的陀螺信号进行标定处理,结果表明仅用原数据量的十分之一即可获得相同的标定精度,减少了测试时间,提高了标定效率。     

Study on spline wavelet finite-element method in multi-scale analysis for foundation

A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D FEM of 7 nodes and 10 nodes are constructed based on the basic formula. Using these proposed elements, the multiscale numerical model for foundation subjected to harmonic periodic load, the foundation model excited by external and internal dynamic load are studied. The results show the pro- posed finite elements have higher precision than the tradi- tional elements with 4 nodes. The proposed finite elements can describe the propagation of stress waves well whenever the foundation model excited by extemal or intemal dynamic load. The proposed finite elements can be also used to con- nect the multi-scale elements. And the proposed finite elements also have high precision to make multi-scale analysis for structure.     

THE WAVELET TRANSFORM OF PERIODIC FUNCTION AND NONSTATIONARY PERIODIC FUNCTION

Some properties of the wavelet transform of trigonometric function, periodic function and nonstationary periodic function have been investigated. The results show that the peak height and width in wavelet energy spectrum of a periodic function are in proportion to its period. At the same time, a new equation, which can truly reconstruct a trigonometric function with only one scale wavelet coefficient, is presented. The reconstructed wave shape of a periodic function with the equation is better than any term of its Fourier series. And the reconstructed wave shape of a class of nonstationary periodic function with this equation agrees well with the function.