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

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

样条厚薄板通用矩形单元

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

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

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

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

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

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

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

样条矩形单元

本文讨论了样条分段(分片)插值,给出了插值基函数的显式及B样条表达式,使得样条插值在局部单元上完成,并用于有限元分析,建立了样条矩形单元。用样条矩形单元求解问题对,可套用有限元的现成计算程序,处理各类边界条件及区域内部的约束条件,此外对于非均匀划分的单元网格,阶梯形边界形状的使用也较方便灵活。     

基于多边形网格的变系数问题边界单元法

提出一种用多边形网格计算二维变系数问题域积分的新型边界单元法。首先,构造了由任意多边组成的多边形网格形函数,用于几何与物理量的插值;其次,用径向积分法将多边形域积分转换成沿多边形周边的线积分,有效解决了各类非规则多边形网格的单元积分难题;最后,三个有关功能梯度材料与结构的数值算例结果显示本文提出的算法和常规有限元相比误差小于1%,说明本文方法具有很高的精度,且由于其单元积分时无需对积分函数或者积分域进行三角化等额外处理,该方法具有很高的效率。     

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

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

Kirchhoff-Love壳的等几何分析

论文利用等几何分析研究了基于Kirchhoff-Love理论的薄壳的静态问题.等几何分析采用等参思想,将精确描述几何形状的NURBS基函数同时作为场变量的插值函数,保证了在分析和网格优化过程中模型的几何精确性,并可以轻易地构造任意高阶连续的单元.该方法具有很高的数值精度.计算结果表明,在等几何分析中,NURBS单元的阶次越高,网格数越多,计算结果越精确.     

建立基于力的变截面梁单元质量矩阵的新方法

基于位移的有限梁单元中三次Hermite插值函数不能有效地描述变截面梁单元内部位移变化,只能通过加密网格增加单元数解决,会造成计算量增大。基于力的有限梁单元由于使用的力插值函数不受截面形状变化的影响,在处理变截面梁时有很大优势,可以得到精确的位移插值函数,利用较少的单元可以达到很高的精度,解决了基于位移的有限梁单元在处理变截面梁时的不足。本文得到了考虑剪切变形的位移插值函数和考虑转动惯量的一致质量矩阵。利用算例验证了本文理论的正确性和高效性。     

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

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

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.     

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

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

IMPROVEMENTS OF GEOMETRICALLY NONLINEAR ANALYSIS ALGORITHMS FOR SPATIAL FRAME STRUCTURES

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

爆炸荷载作用下地下拱形结构动力响应样条小波有限元研究

采用小波有限元方法研究爆炸荷载作用下地下结构的动力响应,克服在模拟过程中由于材料的奇异性、地质条件的复杂性和加载的快速性出现的应力集中和计算效率低下,根据样条有限点法,构造了单向和双向区间B样条圆环扇形小波单元,用尺度函数作为插值函数;结合工程实例,通过Matlab软件编程对爆炸荷载作用地下拱形结构的动力响应进行了数值模拟,并与应用传统有限元程序模拟结果进行对比。结果表明,小波有限元用很少单元取得了较高的精度,计算效率比传统有限元提高一倍。     

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.     

A numerical method for the evaluation of Fourier integrals

This paper suggests the use of spline function interpolation in the evaluation of Fourier integrals. At the same time, the numerical results of some common functions by various interpolation methods and a simplified method of construction of spline function for various boundary conditions are also presented.     

具脱层轴对称层合圆柱壳的振动模态分析

对具环向贯穿脱层的轴对称层合圆柱壳进行振动模态分析.首先,采用Heaviside阶梯函数,构造了一种适合于脱层壳的位移模式.通过对脱层壳的能量分析,应用瑞利--里兹法后,得到用时间函数表示的系统振动控制方程,然后对其求解,得到脱层壳模态分析的特征方程式.算例中,讨论了不同的脱层位置、脱层大小和脱层深度对脱层壳振动模态的影响.     

h-ADAPTIVE ANALYSIS BASED ON MESHLESS LOCAL PETROV-G ALERKIN METHOD WITH B SPLINE WAVELET FOR PLATES AND SHELLS

Using the two-scale decomposition technique, the h-adaptive meshless local Petrov- Galerkin method for solving Mindlin plate and shell problems is presented. The scaling functions of B spline wavelet are employed as the basis of the moving least square method to construct the meshless interpolation function. Multi-resolution analysis is used to decompose the field variables into high and low scales and the high scale component can commonly represent the gradient of the solution according to inherent characteristics of wavelets. The high scale component in the present method can directly detect high gradient regions of the field variables. The developed adaptive refinement scheme has been applied to simulate actual examples, and the effectiveness of the present adaptive refinement scheme has been verified.     

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.