区间五次Hermite样条多小波Euler-Bernoulli梁单元

首先采用区间五次Hermite样条函数,分别构造了三节点梁的边界和中间节点的多小波尺度函数,然后,基于小波多辨分析思想,构建了梁单元位移多尺度近似空间的基函数系;最后,采用最小势能原理,得到弯曲梁的平衡方程,从而构造了区间五次Hermite样条多小波Euler-Bernoulli梁单元。算例结果表明,该小波单元可通过改变尺度来重新划分网格,从而可自由调节单个小波单元的计算精度,其计算精度与在相同网格划分下采用传统三节点Hermite梁单元计算的完全一致;与其它小波单元相比较,该小波单元具有计算简单明了,物理意义明确,易于理解的特点。     

样条厚薄板通用矩形单元

本文利用样条分段插值表示式,借助于最小势能原理,构造了二次样条及三次样条厚薄梁通用单元的位移模式,并推广到二维问题上去,建立了双二次样条及双三次样条厚薄板通用单元。文中给出的算例表明,这些单元均具有自由度少,连续性强,精度高的优点。     

Delaunay多边形单元的有理函数插值格式

本文提出了基于Delaunay多边形化的多边形单元有理函数插值格式。给出了Delaunay多边形化的概念和Delaunay多边形单元有理函数插值形函数的计算表达式。与Delaunay三角化网格不同，Delaunay多边形化网格形成对区域的唯一剖分。Delaunay多边形单元有理函数插值是以Delaunay多边形的顶点作为插值点，构造的有理函数形式插值。Delaunay多边形单元有理函数插值克服了有限元方法中难以构造边数大于4单元多项式形式位移插值的困难。有理函数插值形函数在多边形单元的内部是无穷次光滑的，在多边形的边界上是线性的。在三角形单元和矩形单元上，有理函数插值分别等价于有限元的三角形面积坐标插值和四边形双线性插值。给出了Delaunay多边形有理函数插值在圆域温度分布插值近似中的两个算例。     

一维区间B样条小波单元的构造研究

基于区间B样条小波及小波有限元理论,提出了一种区间B样条小波有限元方法。传统有限元多项式插值被一维区间B样条小波尺度函数取代,进而构造形状函数和单元。与小波Galer-kin方法不同,本文构造的区间B样条小波单元通过转换矩阵将无明确物理意义的小波插值系数转换到物理空间。转换矩阵在小波单元构造过程中起到关键作用,为了保证求解的稳定性,转换矩阵必须非奇异。构造了以区间B样条尺度函数为插值函数的一系列一维区间B样条小波单元。数值算例表明,本文构造的区间B样条小波单元与传统有限元方法相比,在求解变截面,变载荷等问题时具有收敛快和精度高等优势;有效地丰富了小波有限元法单元库。     

区间B样条小波平面弹性及Mindlin板单元构造研究

基于二维张量积区间B样条小波及小波有限元理论,构造了一类用于分析弹性力学平面问题和中厚板问题的C0型区间B样条小波板单元。在二维小波单元的构造过程中,传统多项式插值被二维区间B样条小波尺度函数取代,进而构造形状函数和单元。与小波Galerkin方法不同,本文构造的区间B样条小波单元通过转换矩阵将无明确物理意义的小波插值系数转换到物理空间。区间B样条小波单元同时具有传统有限元和B样条函数数值逼近精度高及多种用于结构分析的基函数的优点。数值算例表明:与传统有限元和解析解相比,本文构造的二维小波单元具有求解精度高,单元数量和自由度少等优点。     

基于三角网格的四阶Hermite型对流守恒格式

提出一种基于三角网格的求解双曲对流方程的高阶守恒型格式.该格式首先在每个三角单元上重构二元三次Hermite插值多项式,以当前时刻单元节点处解的函数值、一阶空间导数值和该单元的积分平均值为插值条件.然后,利用Semi-Lagrange方法得到单元节点处的下一时刻解的函数值及导数值,而下一时刻的解的单元积分平均值由有限体积方法得到.本文所提出的格式将原始CIP方法从结构网格推广到非结构网格上,使得CIP方法能灵活地用于处理复杂边界问题.该格式为显式紧致格式,计算简单且易于实现.数值实验表明,该格式对于光滑解问题能达到四阶空间精度,而对于非光滑解问题能准确地捕捉激波的位置,改进了原始CIP格式的不守恒性.     

自然单元法研究进展

自然单元法是一种基于Voronoi图和Delaunay三角化几何结构,以自然邻点插值为试函数的一种新型数值方法.其既具有无网格方法和经典有限元方法的优点,又克服了两者的一些缺陷,是一种发展前景广阔的求解微分方程的数值方法.自然单元法的形函数满足插值性质,可以像有限元法一样直接施加本质边界条件,不存在基于移动最小二乘拟合的无网格方法不能直接施加本质边界条件的难题.由于自然单元法是无网格方法,可以方便处理有限元方法较难处理的一些问题,例如移动边界和大变形等问题.自然单元法与其他数值方法的最根本区别于其插值格式的不同.将自然邻点插值用于Galerkin过程,就得到基于Voronoi结构的自然单元Galerkin法.自然邻点插值有自然邻点Sibson插值和Laplace插值(非Sibson插值)两种.Laplace插值比Sibson插值在计算上要简单的多,并且不论对凸的或非凸的区域都能精确施加本质边界条件.以Laplace插值为试函数的自然单元法在数值实施上比以Sibson插值为试函数的自然单元法简单.本文对基于Voronoi结构的自然邻点插值和自然单元法的基本思想作了介绍,综述了国内外关于自然单元法的研究成果,总结了自然单元法的优点和尚需解决的问题.     

一种区间B样条小波平板壳单元及应用

基于二维张量积区间B样条小波,构造了一种件能良好的小波平板壳单元．在小波单元的构造过程中,用二维区间B样条小波尺度函数取代传统多项式插值,在所构造的区间B样条平面弹性单元和平面Mindlin板单元的基础上组合而成．区间B样条小波单元同时具有B样条函数数值逼近精度高和多种用于结构分析的基函数的特点．数值算例表明：与传统有限元和解析解相比,构造的小波平板壳单元具有求解精度高,单元数量和自由度少等优点．     

一种考虑剪切变形的平行四边形厚/薄板通用单元

根据Timoshenko二广义位移梁理论，构造了深梁位移场的插值函数。利用斜坐标系与直角坐标系的变换关系、有限条带思想和深梁位移插值函数，构造了一种考虑剪切变形的平行四边形厚／薄板弯曲通用单元的位移(曲率、剪应变、转角、横向位移)插值函数，导出了刚度矩阵和非结点荷载等效力。并对简支阍支方板、Razzaque斜板、四边简支斜交板弯曲进行了数值计算。算例表明此单元有较好的精度，对于薄板不出现剪切闭锁，可适应于目前桥梁建设中大量采用的斜交板桥结构分析。     

轴对称弹性体扭转问题的无网格自然邻接点Petrov-Galerkin法

本文将无网格自然邻接点Petrov-Galerkin 法应用于轴对称弹性体扭转问题的求解．无网格自然邻接点Petrov-Galerkin 法采用自然邻接点插值构造试函数，并且采用三角形线性单元的形函数作为加权残值法的加权函数．自然邻接点插值构造的试函数满足Kronecker delta 函数性质，因此本质边界条件的施加十分方便．由于几何形状和边界条件的轴对称特点，原来的空间问题简化为二维问题求解，因此计算时只需要横截面上离散节点的信息．数值算例结果表明，所提出的方法对求解轴对称弹性体扭转问题是行之有效的．     

非线性动力方程的三次样条-增维精细算法

提出一种针对非线性动力方程的改进精细积分方法。该方法是在时间步长内采用分段的三次样条函数拟合非齐次项,保持高精度拟合的同时避免了求导运算和高次多项式插值带来的Runge现象。通过引入4×2个变量将动力方程增加四维转化为齐次方程,并建立相应的通解格式,避免了状态空间下系统矩阵求逆。将指数矩阵分为四个子模块,利用各模块的特点分别进行理论推导及基于精细积分法进行分步、分块计算得到相应的理论解和高精度数值解,无需反复计算整个指数矩阵,提高了解算效率。针对含未知状态量的非齐次项,引入预测-校正的方法进行迭代求解。数值计算结果表明了本文方法的有效性。     

Characteristics method with cubic‐spline interpolation for open channel flow computation

In the framework of the specified‐time‐interval scheme, the accuracy of the characteristic method is greatly related to the form of the interpolation. The linear interpolation was commonly used to couple the characteristics method (LI method) in open channel flow computation. The LI method is easy to implement, but it leads to an inevitable smoothing of the solution. The characteristics method with the Hermite cubic interpolation (HP method, originally developed by Holly and Preissmann, 1977) was then proposed to largely reduce the error induced by the LI method. In this paper, the cubic‐spline interpolation on the space line or on the time line is employed to integrate with characteristics method (CS method) for unsteady flow computation in open channel. Two hypothetical examples, including gradually and rapidly varied flows, are used to examine the applicability of the CS method as compared with the LI method, the HP method, and the analytical solutions. The simulated results show that the CS method is comparable to the HP method and more accurate than the LI method. Without tackling the additional equations for spatial or temporal derivatives, the CS method is easier to implement and more efficient than the HP method. Copyright © 2004 John Wiley & Sons, Ltd.     

IMPROVEMENTS OF GEOMETRICALLY NONLINEAR ANALYSIS ALGORITHMS FOR SPATIAL FRAME STRUCTURES

以几何精确梁理论为基础,分别采用高阶拉格朗日插值和埃米特插值构造高精度空间梁单元。提出基于单元层次平衡迭代的自由度凝聚方法,以保证单元的通用性。实现了基于载荷控制或柱面弧长控制的结构几何非线性分析算法。算例研究结果表明,提出的改进方法不但提高了计算效率,而且还具有较高的数值稳定性;特别是基于三次埃米特插值构造的单元表现出较好的性态,适用于结构屈曲后分析。     

高层筒体结构整体稳定及二阶位移分析的改进条元法

根据柱壳理论,构造了一种柱壳曲条,本文结合柱壳曲条和平壳条元求解高层筒体结构的整体稳定及二阶位移。采用三次H erm ite插值函数模拟条元横截面的翘曲位移变化,能较好地反映筒体受力“剪切滞后”效应;采用一族能较好地逼近弯剪型变形曲线的正交多项式作基函数来描述位移沿竖向变化。用最小势能原理建立稳定及二阶位移分析方程。该方法适用于任意平面形状的高层建筑筒体结构及剪力墙结构的稳定及二阶位移分析。与其它方法相比,该方法具有精度高、通用性强、计算量小等优点。     

Numerical solution of Poisson equation with wavelet bases of Hermite cubic splines on the interval

A new wavelet-based finite element method is proposed for solving the Poisson equation. The wavelet bases of Hermite cubic splines on the interval are employed as the multi-scale interpolation basis in the finite element analysis. The lifting scheme of the wavelet-based finite element method is discussed in detail. For the orthogonal characteristics of the wavelet bases with respect to the given inner product, the corresponding multi-scale finite element equation can be decoupled across scales, totally or partially, and suited for nesting approximation. Numerical examples indicate that the proposed method has the higher efficiency and precision in solving the Poisson equation.     

A finite wavelet domain method for wave propagation analysis in thick laminated composite and sandwich plates

An efficient numerical method is developed for the simulation of three dimensional transient dynamic response in thick laminated composite and sandwich plate structures involving very high frequencies and wave numbers. The proposed method incorporates Daubechies wavelet scaling functions for the interpolation of the in-plane displacements with a Galerkin formulation. It further explores the orthonormality and compact support of wavelet scaling functions to produce near diagonal consistent mass matrices and banded stiffness matrices. Hence, an uncoupled equivalent discrete spatial dynamic system is formulated, synthesized and rapidly solved in the wavelet domain using an explicit time integration scheme. The in-plane wavelet interpolation is further combined with an efficient high order layerwise laminate plate theory, that implements Hermite cubic splines for the through-the-thickness approximation of displacement fields. Numerical results are presented on the prediction of guided waves in laminated and thick sandwich composite plates and compared with respective solutions obtained by analytical, semi-analytical and time domain spectral element models. The method yielded higher convergence rates and substantial reductions in computational effort compared to respective time domain spectral finite elements.     

用样条法研究高速滑行体的运动

研究高速滑行体在规则波浪中滑行时，得到了二阶常微分方程组，本文应用样条插值求积方法和三次样条配点计算，求解滑行体在不同波浪条件下的纵倾和升沉的响应．实例计算结果与实验值比较一致．     

Refined Series Methodology for the Fully Coupled Thin-Walled Laminated Beams Considering Foundation Effects

The refined power series solutions are presented for the coupled static analysis of thin-walled laminated beams resting on elastic foundation. For this purpose, the elastic strain energy considering the material and structural coupling effects and the energy including the foundation effects are constructed. The equilibrium equations and the force-displacement relationships are derived from the extended Hamilton's principle, and the explicit expressions for displacement parameters are presented based on power series expansions of displacement components. Finally, the member stiffness matrix is determined by using the force-displacement relationships. For comparison, the finite element model based on the Hermite cubic interpolation polynomial is presented. In order to verify the accuracy and the superiority of the laminated beam element developed by this study, the numerical solutions are presented and compared with results obtained from the regular finite beam elements and the ABAQUS's shell elements. The influences of the fiber angle change and the boundary conditions on the coupled behavior of laminated beams with mono-symmetric I-sections are investigated.     

13加速度计捷联惯导姿态解算方法(英文)

为了减小全加速度计捷联惯导(亦称无陀螺惯导)角速度重构误差以提高姿态解算精度,提出了一种13加速度计无陀螺惯导配置方案。针对合理设计的13加速度计捷联惯导能同时提供角速度和角加速度信息的特点,提出了一种基于角速度/角加速度信息的捷联姿态算法,采用艾尔米特插值方法进行离散角速度插值重构,利用重构后的角速度信息进行姿态矩阵求解。分析了新姿态算法的基本步骤,并进行了典型圆锥运动仿真实验。结果表明,采用新方法求解姿态矩阵的精度提高40%以上。不同于传统姿态算法仅利用角速度信息,新的姿态算法利用更多的信息,能获得更好的精度。     

THE CONSTRUCTION OF WAVELET-BASED TRUNCATED CONICAL SHELL ELEMENT USING B-SPLINE WAVELET ON THE INTERVAL

Based on B-spline wavelet on the interval (BSWI), two classes of truncated conicalshell elements were constructed to solve axisymmetric problems, i.e. BSWI thin truncated conicalshell element and BSWI moderately thick truncated conical shell element with independent slope-deformation interpolation. In the construction of wavelet-based element, instead of traditionalpolynomial interpolation, the scaling functions of BSWI were employed to form the shape functionsthrough the constructed elemental transformation matrix,and then construct BSWI element viathe variational principle. Unlike the process of direct wavelets adding in the wavelet Galerkinmethod, the elemental displacement field represented by the coefficients of wavelets expansionwas transformed into edges and internal modes via the constructed transformation matrix. BSWIelement combines the accuracy of B-spline function approximation and various wavelet-basedelements for structural analysis. Some static and dynamic numerical examples of conical shellswere studied to demonstrate the present element with higher efficiency and precision than thetraditional element.