势问题的重构核粒子边界无单元法

将重构核粒子法和势问题的边界积分方程方法结合,提出了势问题的重构核粒子边界无单元法. 推导了势问题的重构核粒子边界无单元法的公式,研究其数值积分方案,建立了重构核粒子边界无单元法的离散化边界积分方程,并推导了重构核粒子边界无单元法的内点位势的积分公式. 重构核粒子法形成的形函数具有重构核函数的光滑性,且能再现多项式在插值点的精确值,所以该方法具有更高的精度. 最后给出了数值算例,验证了所提方法的有效性和正确性. }      

复变量重构核粒子法与有限元法耦合解弹性力学问题

提出了弹性力学的复变量重构核粒子法与有限元法的耦合法(CVRKPM/FEM).采用场量耦合试函数法将弹性力学的复变量重构核粒子法与有限元法进行耦合,详细推导了在整个求解域上的耦合公式.最后通过数值算例证实了本文所提弹性力学的复变量重构核粒子法与有限元的耦合法的有效性.本文的耦合法不仅可以很方便地施加本质边界条件,而且可以充分利用无网格方法和有限元法的优势,弥补各自不足以提高计算效率.     

势问题的复变量重构核粒子法

在重构核粒子法的基础上,引入复变量,提出复变量重构核粒子法,在构造形函数时采用一维基函数建立二维问题的修正函数.应用于势问题,具有计算量小、精度高的优点.数值算例证明了方法的有效性.     

弹性力学的重构核粒子边界无单元法与有限元的耦合法

将重构核粒子边界无单元法(RKP-BEFM)与有限元法(FEM)耦合,形成求解具有区域特征的弹性力学问题的重构核粒子边界无单元与有限元的耦合方法RKP-BEF/FE.推导了重构核粒子边界无单元与有限元耦合方法的离散化公式,建立了节点未知量的耦合方程.重构核粒子边界无单元法和有限单元法的较高精度保证了这一直接耦合方法的成功实现与求解精度.最后给出了平面问题的数值算例,验证了提出的耦合方法RKP-BEF/FE的有效性.     

饱和多孔介质动力学大变形分析耦合对流粒子域插值方法

提出了饱和多孔介质大变形动力学响应分析的耦合对流粒子域插值方法(Coupling convectedparticle domain interpolation method,CCPDI)。采用u-p形式的控制方程和超弹性材料本构关系对具有饱和多孔特性的介质进行了大变形动力学模拟,建立了耦合对流粒子域插值方法的弱形式离散求解方程并给出了该方法的计算流程。通过数值算例,验证了所提出的耦合对流粒子域插值方法的正确性。本文工作为生物软组织、肌肉、骨骼和其它一些具有饱和多孔特性的软物质的几何非线性动力学行为分析奠定了基础。     

应力波数值计算中的SPH方法

对一维波动方程的SPH(smoothed particle hydrodynamics)格式和有限差分格式进行比较,并采用SPH法模拟了一维应力/应变波, 获得1个可衡量SPH法模拟应力波准确性的重要指标。结果表明,SPH法模拟应力波传播中采用的光滑长度必须不小于粒子间距;采用B-样条核函数和高斯型核函数能够获得良好的应力波图像,而二次型核函数不能,因此二次型核函数不适用于冲击动力学的数值计算。     

基于核重构思想的最小二乘配点型无网格方法

介绍重构核点法的基本原理和近似函数的构造方法，并基于核重构思想，应用配点法和最小二乘原理，离散微分方程，建立求解的代数方程，提出了一种基于核重构思想的最小二乘配点型无网格方法．与一般配点法相比，该方法的系数矩阵是有对称正定的，计算精度高，稳定性好．该方法的实施不需要背景网格，不需要进行高斯积分，与Galerkin法相比，具有计算量小、边界条件处理简单的特点，是一种真正的无网格法．对该方法构造过程中的近似函数及其导数的计算、修正函数的计算及方法的实现等问题进行了探讨．文中结合若干典型算例，检验了该方法的有效性．     

用局部Petrov-Galerkin法分析薄板自由振动

利用薄板振型方程的等效积分弱形式和对振型函数采用移动最小二乘近似函数进行插值，本文进一步研究了无网格局部Petrov-Galerkin方法在薄板自由振动问题中的应用。它不需要任何形式的网格划分，所有的积分都在规则形状的子域及其边界上进行。在插值近似时，采用虚拟-实际节点值变换方法直接引入本质边界条件。通过数值算例和与其他方法的结果进行比较，表明无网格局部Petrov-Galerkin法求解弹性薄板自由振动问题具有收敛性好、精度高等一系列优点。     

基于表面力模型的液滴冲击弹性基底SPH模拟

采用光滑粒子动力学SPH方法建立液滴冲击弹性基底的流固耦合数值模型,给出描述粘性流体和弹性固体运动的SPH离散方程和数值处理格式,引入人工耗散项来抑制标准SPH方法的数值震荡。为模拟液滴的表面张力效应,通过精确检测边界粒子,采用拉格朗日插值方法计算表面法向量和曲率,结合界面理论中的连续表面力CSF方法,建立了适用于自由表面液滴的表面力模型,方形液滴变形的模拟结果与拉普拉斯理论解吻合较好。随后,采用SPH流固耦合模型模拟1.0 mm直径水滴以不同速度(0.2 m/s～3.0 m/s)冲击两种薄板型基底,分析了基底弹性变形对液滴铺展、收缩以及回弹行为的影响。     

弹性薄板的广义多维微分求积分析方法

多维微分求积DC(differential cubature)法是近年迅速发展起来的一种处理多维微分方程的高精度数值方法。本文围绕多元函数插值适定性这一基本问题，探讨了实现DC法时插值结点组和插值空间应满足的条件，研究了离散结点组选择及其构造适定试函数的方法。提出了一种沿各坐标方向可以布置不同数量结点的广义交错网格DC方法，并给出了确定相应的插值空间及其选择试函数的具体方案。通过弹性薄板变形分析研究了广义交错网格DC法解决实际结构力学问题的可靠性及其潜力。研究结果表明，广义交错网格DC法较传统的交错网格DC格式具有更强的适用性和灵活性，采用数量少得多的离散结点即可达到传统交错网格DC法相同的数值精度，表现出精度高和计算工作量小的优点。     

广义概率密度演化方程的再生核质点加密算法

采用一般质点近似和再生核质点近似表示系统响应量,给出了动力系统响应量的一般表达式。在此基础上,发展了一类求解广义概率密度演化方程的再生核质点加密算法,给出了详细求解步骤。以单自由度系统为例,从响应概率密度的角度考察了再生核质点加密算法的精度。以多自由度框架结构为例,验证了再生核质点加密算法求取非线性随机动力系统响应概率密度的正确性。     

Elastodynamic analysis of a crack at an arbitrary angle to the graded interfacial zone in bonded half-planes under antiplane shear impact

A solution is provided for the elastodynamic problem of a crack at an arbitrary angle to the graded interfacial zone in bonded media under the action of antiplane shear impact. The interfacial zone is modeled by a nonhomogeneous interlayer with the spatially varying shear modulus and mass density in terms of power functions between the two dissimilar, homogeneous half-planes. Based on the use of Laplace and Fourier integral transforms and the coordinate transformations of basic field variables, formulation of the transient crack problem is reduced to solving a Cauchy-type singular integral equation in the Laplace transform domain. The crack-tip response in the physical domain is recovered via the inverse Laplace transform and the values of dynamic mode III stress intensity factors are obtained as a function of time. A comprehensive parametric study is then presented of the effects of crack obliquity on the overshoot behavior of the transient crack-tip response, by plotting the peak values of the dynamic stress intensity factors versus the crack orientation angle for various material and geometric combinations of the bonded system.     

考虑在地基固结和流变的上部结构和地基基础的共同作用问题

本文提出一种能同时考虑地基的瞬时变形,团结变形和流变的地基基础与上部结构共同作用的一种新的计算方法,在这种方法中,土的计算考虑了Biot固结理论和流变,上部结构利用子结构法进行凝聚可得边界刚度矩阵和边界力,利用位移协讯条件则可以得到共同作用的控制方程,对上述方程进行Laplace变换,可以得到Laplace变换域内的控制方程,在Laplace变换域内对上述方程进行数值求解并进行相应的Laplace逆变换则可以得到时间境内任意时刻的解。     

Laplace transform-boundary element coupling method for viscous fluid-structure impact analysis

Based on linearized 2-D Navier-Stokes equation, a Laplace transform-boundary element coupling method for viscous fluid-structure impact analysis is proposed. Under assumption of incompressibility for the fluid, the corresponding equivalent boundary integral equation in terms of the potential function and stream function is first established by Lamb's transform in the Laplace transform domain. It enables us to solve impact water problems in frequency domain by the boundary element method, in which the effect of viscous flow on the dynamic response can be taken into account. Then a complete solution of the problem under consideration in time domain is obtained by means of Durbin's formulas for the numerical inversion of the Laplace transform. Finally, a practical example is given to confirm the validity of the present method. Project supported by the National Defence Foundation of Science & Technology of China (No. J14. 8. 1. JW0515).     

轴向运动弦线横向振动的频域分析

轴向运动弦线是多种工程系统的模型。为明确轴向运动横向振动的频域特性，及探索频域方法的应用特点．本文用频域方法分析轴向运动弦线的横向振动。基于轴向运动弦线横向振动方程和边界条件．通过Laplace变换导出频率域中的控制方程，并将该控制方程和边界条件用状态变量表示。由状态空间中的控制方程导出特征方程，从而求出固有频率。由轴向运动弦线的矩阵函数计算得到系统的传递函数，然后用留数定理计算传递函数的Laplace逆变换．这样就可以得到时域响应。最后分析了轴向运动弦线的横向共振，若简谐外激励的频率与系统固有频率相同，系统响应将随时间无限增加。     

线黏弹性球面发散应力波的频率响应特性

基于线黏弹性球面波Laplace域的理论解, 得到了不同传播距离处粒子速度、粒子位移、应力、应变等力学量的传递函数。以标准线性固体模型为例, 重点讨论了粒子速度频率响应函数的传播特征, 指出随着传播距离的增加, 粒子速度幅频响应函数的高频响应会低于低频响应, 而在理想弹性条件下, 粒子速度幅频响应函数的高频响应一直高于低频响应。以弹性半径为0.025 m的空腔爆炸为例, 采用Laplace数值逆变换方法对粒子速度波形的演化进行了分析, 给出了粒子速度强间断幅值及粒子速度峰值随传播距离变化的衰减规律曲线, 指出黏弹性介质中粒子速度幅值的衰减曲线介于理想弹性介质中粒子速度幅值衰减曲线和黏弹性介质中粒子速度强间断幅值衰减曲线之间。     

求解饱和半空间上弹性圆板固结沉降的积分方程

本文采用解析方法分析了弹性圆板在饮和半空间上的固结沉降。考虑弹性圆板与饮和半空间的接触面上无摩擦力,且饱和半空间表面为全部透水的。运用Ｂｉｏｔ固结理论和积分方程技术,在Ｌａｐｌａｃｅ变换域上建立了弹性圆板固结沉降的对偶积分方程,并化此对偶积分方程为第二类Ｆｒｅｄｈｏｌｍ积分方程。通过对其核函数的有效数值发得到第二类Ｆｒｅｄｈｏｌｍ积分方程的解,再利用Ｌａｐａｃｅ反演技术获得弹性板在时间域中的固结沉     

The elastodynamic solution for a solid sphere and dynamic stress-focusing phenomenon

This paper presents an analytical method of solving the elastodynamic problem of a solid sphere.The basic solution of the elastodynamic problem is decomposed into a quasi-static solution satisfying the inhomogeneous compound boundary conditions and a dynamic solution satisfying the homogeneous compound boundary conditions.By utilizing the variable transform,the dynamic equation may be transformed into Bassel equation.By defining a finite Hankel transform,we can easily obtain the dynamic solution for the inhomogeneous dynamic equation.Thereby,the exact elastodynamic solution for a solid sphere can be obtained.From results carried out,we have observed that there exists the dynamic stress-focusing phenomenon at the center of a solid sphere under shock load and it results in very high dynamic stress-peak.