基于拉氏变换的弹性动力学插值型重构核粒子法

插值型重构核粒子法的形函数具有离散点插值特性和不低于核函数的高阶光滑性,因而不仅可以直接施加本质边界条件,同时也保证了较高的计算精度．本文将弹性动力学方程作拉氏变换后,在变换域内用插值型重构核粒子法求解,最后再借助Durbin数值反演方法求得时间域的解．针对典型的弹性动力学问题,给出了插值型重构核粒子法的数值算例,并验证了本文方法的有效性．     

光滑粒子法与有限元的耦合算法及其在冲击动力学中的应用

扼要讨论了光滑粒子法的离散思想,充分利用光滑粒子法和有限元方法各自的优点,提出了一种初始时刻用有限元建模,计算过程中大变形单元自动转换为光滑粒子的耦合算法。高速碰撞的系列算例说明,耦合算法不但适宜于计算大变形冲击动力学问题,而且由于集两种方法的优点于一身,可以更高效地模拟一些高速碰撞问题,提高计算效率。     

拉氏变换求解梁的挠曲线方程

运用拉普拉斯变换求解梁的挠曲线近似微分方程, 并利用坐标系平移变换导出了分段梁挠曲线方程的一般形式, 通过算例验证简述了用此方法可方便地根据弯矩方程和边界条件求出梁各段挠曲线方程的表达式.     

光滑粒子动力学方法的发展与应用

光滑粒子动力学(smoothed particle hydrodynamics,SPH)是一种拉格朗日型无网格粒子方法,已经成功地应用到了工程和科学的众多领域.SPH使用粒子离散及代表所模拟的介质,并且基于粒子体系估算和近似介质运动的控制方程.本文分析和综述了SPH模拟方法的发展历程、数值方法与应用进展.介绍了SPH方法的基本思想;从连续性、边界处理、稳定性和计算效率4个方面阐述了SPH方法的研究现状;介绍了SPH方法近年来在可压缩流动、不可压缩流动以及弹塑性材料高速变形与失效方面的一些典型应用;并对SPH方法的发展与应用进行了预测与展望.      

光滑粒子模拟方法在超高速碰撞现象中的应用

简要介绍了基于黎曼解的光滑粒子法,并将改进的SPH方法应用于超高速碰撞,对二维轴对称条件下的弹丸超高速碰撞薄板问题进行了数值模拟,研究了靶板厚度、弹丸速度、弹丸形状等因素对形成碎片云的影响。通过与实验数据比较,该算法模拟的碎片云的形状及特征与实验相吻合,验证了光滑粒子法对冲击动力学问题数值模拟的有效性。     

插值摄动法及其在力学中的应用

本文研究了最高阶导数乘以小参数，或出现奇点的微分方程的定解问题，用插值摄动法求得了一级近似解，它和通常的奇异摄动法（匹配法、多尺度法）的一级近似解的精度相同。     

特征线法在三维弹／粘塑性波数值解析中的应用

根据广义特征理论建立了求解三维非线性应力波传播问题的特征关系式,并采用特征线法对纵向冲击载荷下弹／粘塑性矩形截面棒的三维应力波传播过程进行了数值模拟。     

改进特征线法在水下爆炸正反演问题中的应用

介绍了改进特征线法及其在水下爆炸正反演问题中的应用。将改进特征线法用于TNT炸药球水下爆炸的模拟,并与实验数据和AUTODYN程序的结果比较,结果表明,改进特征线法求解水下爆炸问题具有较好的准确性,且能够捕捉到流场中弱冲击的传播。而后应用改进的特征线法将正问题中的冲击波参数作为初值条件,利用逆序差分格式,对水气界面做反演求解;将反演的水气界面与正问题中的界面作比较,结果显示两者在近场吻合良好。     

抗差自适应插值粒子滤波及其在组合导航中的应用

针对粒子滤波存在的粒子退化和重要性密度函数难以选取的问题,在吸收抗差自适应滤波、二阶插值滤波和粒子滤波算法优点的基础上,提出了一种新的抗差自适应插值粒子滤波算法。该算法利用二阶插值滤波算法得到重要性密度函数,通过抗差自适应因子实时控制动力学模型误差及观测异常对导航解的影响。将该算法应用于SINS/CNS/SAR组合导航系统进行计算仿真,并与经典的粒子滤波算法进行比较分析。结果表明,提出的滤波算法得到的姿态误差控制在[-0.3′,+0.3′],速度误差控制在[-0.4 m/s,+0.4 m/s],位置误差控制在[-5 m,+5 m],性能明显优于经典的粒子滤波算法。新的滤波算法不但能够有效地抑制粒子退化,而且能够有效地控制动力学模型误差及观测异常的影响,提高了组合导航的滤波精度。     

复合材料层合板自由振动的插值型重构核粒子法

基于一阶剪切理论,利用插值型重构核粒子法对复合材料层合板的自由振动进行了数值分析。由于插值型重构核粒子法的形函数具有离散点插值性和全域高阶光滑性,可直接施加本质边界条件,保证计算精度。通过本文给出的方法,对不同边界条件、不同材料参数、不同边厚比的层合板的振动频率进行了计算,并与已有文献的结果进行了对比,其相对误差均在4%以内,典型算例的计算验证了本文方法的可行性和有效性。     

AUGMENTED LAGRANGIAN METHOD AND OPEN BOUNDARY CONDITIONS IN 2D SIMULATION OF POISEUILLE–BÉNARD CHANNEL FLOW

The main objective of this study is to compare the influence of different boundary conditions upon the incompressible Poiseuille –Bénard channel flow (PBCF) in a 2D rectangular duct heated from below. In a first technical part the algorithm used to carry out this work, based on the augmented Lagrangian method, is presented. The implementation details of the five different open boundary conditions (OBCs) and the periodic boundary conditions (PBCs) tested in the present paper are also given. The study is then carried out for 1800<Ra≤ 10,000, 0<Re≤10 and 0·67≤Pr≤6·4. The five selected OBCs, applied at the outlet of the computational domain, respectively express the following conditions: a square profile for the velocity (OBC1), mass conservation (OBC2), zero second derivative of the horizontal velocity component (OBC3), a mixed boundary condition combining Dirichlet and Neumann conditions (OBC4) and an Orlanski-type boundary condition (OBC5). A good estimation of the perturbation amplitude and of the length of the perturbed zone at the outlet boundary is proposed. It is shown that OBC5 causes very little perturbation in the recirculating flow compared with the other OBCs. © 1997 John Wiley & Sons, Ltd.     

NEW TREATMENT OF ESSENTIAL BOUNDARY CONDITIONS IN EFG METHOD BY COUPLING WITH RPIM

One of major difficulties in the implementation of meshfree methods using the moving least square (MLS) approximation, such as element-free Galerkin method (EFG), is the imposition of essential boundary conditions as the approximations do not pass through the nodal parameter values. Another class of meshfree methods based on the radial basis point interpolation can satisfy the essential boundary conditions exactly since its approximation function passes through each node in an influence domain and thus its shape functions possess the properties of delta function. In this paper, a coupled element-free Galerkin(EFG)-radial point interpolation method (RPIM) is proposed to enhance their advantages and avoid their disadvantages. Discretized equations of equilibrium are obtained in the RPIM region and the EFG region, respectively. Then a collocation approach is introduced to couple the RPIM and the EFG method. This method satisfies the linear consistency exactly and can maintain the stiffness matrix symmetric. Numerical tests show that this method gives reasonably accurate results consistent with the theory.     

光滑节点插值法:计算固有频率下界值的新方法

将光滑节点插值法用于悬臂梁的静力学,并首次用于旋转柔性梁的频率分析. 采用梯度光滑技术,用线性插值形函数描述梁的位移场,求解4 阶微分方程. 在静力学分析中,将该方法所得梁中各点位移与假设模态法、有限元法及解析解的结果对比,可知该方法虽用简单的线性插值形函数描述梁的位移场,但精度却很高. 进一步研究表明,采用模态高于9 阶的假设模态法会使刚度阵条件数变差,导致结果发散. 在频率分析中,与有限元法、假设模态法和解析解对比,表明该方法一个重要特性:能提供固有频率的下界值,而有限元法和假设模态法只能提供固有频率的上界值,说明该方法结合有限元法在处理无解析解的问题时可以从上下界最大程度的逼近真实解,提高精度. 光滑节点插值法具有形函数结构简单、独立变量少且能提供固有频率下界值的特性,因此,具有较高的推广及应用价值.      

非齐次边界条件下轴向运动梁的非线性振动

轴向运动系统的横向非线性振动一直是国内外研究的热点课题之一.目前相关研究大都是针对齐次边界条件的.但是在工程实际中,非齐次边界条件更为常见,而针对非齐次边界条件的研究相对较少.为深入研究非齐次边界条件对轴向运动系统横向非线性振动的影响,本文以轴向变速运动黏弹性Euler梁为例,引入由黏弹性引起的非齐次边界条件,同时还引入由轴向加速度引起的径向变化张力,建立梁横向振动的积分-偏微分型运动方程,并导出了相应的非齐次边界条件.采用直接多尺度法分析了梁的次谐波参数共振.由可解性条件得到了梁的稳态响应,并根据Routh-Hurvitz判据确定了系统稳态响应的稳定性.通过数值例子讨论了黏弹性系数,轴向运动速度,轴向速度脉动幅值和非线性系数对幅频响应的影响,并详细对比分析了非齐次边界条件和齐次边界条件对幅频响应的影响.结果表明:随着黏弹性系数的增大,非齐次边界条件下的零解失稳区域和稳态响应幅值比齐次边界条件下的失稳区域和幅值大,非齐次边界条件对高阶次谐波参数共振的影响更加显著.最后,引入微分求积法来验证直接多尺度法的近似解结果.      

格子Boltzmann方法模拟多孔介质惯性流的边界条件改进

格子Boltzmann方法可以有效地模拟水动力学问题,边界处理方法的选择对于可靠的模拟计算至关重要.本文基于多松弛时间格子Boltzmann模型开展了不同边界条件下,周期对称性结构和不规则结构中流体流动模拟,阐述了不同边界条件的精度和适用范围. 此外,引入一种混合式边界处理方法来模拟多孔介质惯性流, 结果表明:对于周期性对称结构流动模拟,体力格式边界条件和压力边界处理方法是等效的,两者都能精确地捕捉流体流动特点; 而对于非周期性不规则结构,两种边界处理方法并不等价,体力格式边界条件只适用于周期性结构;由于广义化周期性边界条件忽略了垂直主流方向上流体与固体格点的碰撞作用,同样不适合处理不规则模型;体力-压力混合式边界格式能够用来模拟周期性或非周期性结构流体流动,在模拟多孔介质流体惯性流时,比压力边界条件有更大的应用优势,可以获得更大的雷诺数且能保证计算的准确性.      

CRACK FACE CONTACT PROBLEM ANALYSIS USING THE SCALED BOUNDARY FINITE ELEMENT METHOD

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.     

????????????????????????е????

提出了将杂交边界点法和双重互易法结合求解势问题的一种新的算法. 将势问题的解分为通解和特解两部分，通解使用 杂交边界点方法求解，特解则利用局部径向基函数近似. 该方法输入数据只是求解域上离散 的点，不需要额外的方程来计算域内物理量，后处理十分简便. 数值算例表明了该方法的稳 定性和有效性.     

STABILITY AND ACCURACY OF NUMERICAL BOUNDARY CONDITIONS IN AEROELASTIC ANALYSIS

This paper analyses the accuracy and numerical stability of coupling procedures in aeroelastic modelling. A two-dimensional model problem assuming unsteady inviscid flow past an oscillating wall leads to an even simpler one-dimensional model problem. Analysis of different numerical algorithms shows that in general the coupling procedures are numerically stable, but care is required to achieve accuracy when using very few time steps per period of natural oscillation of the structure. The relevance of the analysis to fully three-dimensional applications is discussed. © 1997 by John Wiley & Sons, Ltd.     

压电材料中裂纹面电边界条件的适用性

先分析理想裂纹面电边界条件(包括导通型裂纹和不导通型裂纹),然后分析缺陷的厚度及其对尖端场的影响,用封闭形式的解表示出了缺陷厚度的影响.结果指出,裂纹长厚比对计算结果是有影响的,除非裂纹极扁,不可通裂纹边界条件可以给出合理的结果,而可通边界条件只有在一些限制条件下才适用.     

比例边界坐标插值方法在谱元法中的应用——无穷域Euler方程的数值模拟

将比例边界坐标插值方法引入谱元法, 构成比例边界谱单元, 对无穷域Euler方程进行数值模拟.阐述了比例边界谱单元的基本使用方法以及基于比例边界谱元的Runge-Kutta间断Galerkin方法求解Euler方程的过程;计算了无穷域圆柱和NACA0012翼型绕流问题, 并与已有结果进行了比较, 显示了计算结果的正确性.用基于比例边界谱元的间断Galerkin方法求解无穷域Euler方程时, 最多只需将求解域划分为2个子域, 避免了一般谱方法将求解域划分为9个或者27个子域的麻烦. 比例边界谱单元为无穷域Euler方程的直接求解提供了一个可供参考的方法.