非均质中厚板的无网格LRPIM动力学分析

用局部加权残值法建立了非均质中厚板的局部径向点插值离散系统方程,采用无网格局部径向点插值法分析了非均质中厚板的自由振动和强迫振动问题。用径向基函数耦合多项式基函数来近似试函数,用四次样条函数做为加权残值法中的权函数。所构造的形函数具有Kronecker delta性质,可以很方便地施加本质边界条件。该方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行。在计算过程中,取积分中的高斯点的材料参数来模拟问题域材料特性的变化。计算结果表明,利用该方法计算非均质中厚板的自由振动和强迫振动问题可以得到具有较高精度的解。     

用无网格局部径向点插值法分析非均质中厚板的弯曲问题

用无网格局部径向点插值法分析了非均质中厚板的弯曲问题.利用虚位移原理推导了中厚板的离散系统方程.采用径向基函数耦合多项式基函数来近似试函数,用四次样条函数作为加权残值公式中的权函数.所构造成的形函数具有Kronecker delta性质,可以很方便地施加本质边界条件.此方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行,是一种真正的无网格方法.在计算过程中,取积分中的高斯点的材料参数来模拟问题域材料特性的变化.算例结果表明这种无网格方法具有效率高、精度高和易于实现等优点.     

用无网格LRPIM分析中厚板弯曲问题

用径向基函数构造无网格点插值法的形函数,插值函数具有Kronecker delta函数性质,因此可以很方便地施加本质边界条件.利用无网格局部径向点插值方法分别对一个对边固支另对边简支中厚板和一个悬臂中厚板的弯曲进行了分析计算.该方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行,是一种真正的无网格方法.算例表明:将无网格局部径向点插值法应用于计算中厚板的弯曲问题,所求得的位移场和应力场都是光滑的;在径向基函数的基础上,附加多项式大大提高了插值精度;所得结果与弹性力学理论解以及有限元解都十分吻合.     

弹塑性扭转问题具多项式基的径向点插值无网格法

对于弹塑性扭转问题描述的椭圆变分不等式,采用具多项式基的径向点插值法无网格方法与Uzawa方法耦舍,得到了带松弛因子的离散迭代算法,并给出了数值算例,分析了参数对结果的影响.通过与有限元法比较,表明该方法是求解弹塑性扭转问题的有效的方法之一.     

中厚板弯曲分析的无网格弱-强式法(MWS)

基于局部弱式和强式配点相结合的无网格弱-强式法(meshfree weak-strong method,MWS)求解中厚板问题.MWS法对问题域使用整体离散节点表征和强形式配点法进行计算,在自然边界条件上或靠近自然边界条件的区域采用局部弱形式Petrov-Cralerkin法计算,用移动最小二乘法或径向点插值法来构造形函数,是一种理想的真正无网格法.采取MWS法,文中计算了中厚板的弯曲问题和能量误差.算例结果和对比分析表明,无网格弱-强式法(MWS)可以自然协调处理两类边界条件,计算效率高、数值结果稳定;对计算域采用规则节点布置,其解与弹性力学理论解以及有限元解都吻合很好.     

大跨桥梁参数识别响应面方法中的近似函数及样本选取

以某大跨斜拉桥实验室物理模型设计参数与动力特性之间复杂的隐式关系为对象,分析响应面模型中近似函数和试验样本对回归模型精度的影响。首先推导多项式和径向基函数响应面模型的建模方法,然后基于斜拉桥三维空间模型的有限元动力分析,比较分析完全二阶多项式与几种径向基函数模型在设计域和扩展设计域内以及噪声污染情况下的回归精度。最后对几种少量样本的试验设计方法和响应面试验设计方法进行了分析。结果表明,径向基函数模型各方面性能优于多项式模型,响应面实验设计方法优于其他方法,少量样本可采用D优化设计,较多样本情况下以中心复合设计较佳。     

基于分区径向基函数配点法的大变形分析

无网格法因为不需要划分网格, 可以避免网格畸变问题,使得其广泛应用于大变形和一些复杂问题. 径向基函数配点法是一种典型的强形式无网格法,这种方法具有完全不需要任何网格、求解过程简单、精度高、收敛性好以及易于扩展到高维空间等优点,但是由于其采用全域的形函数, 在求解高梯度问题时 存在精度较低和无法很好地反应局部特性的缺点. 针对这个问题,本文引入分区径向基函数配点法来求解局部存在高梯度的大变形问题. 基于完全拉格朗日格式,采用牛顿迭代法建立了分区径向基函数配点法在大变形分析中的增量求解模式.这种方法将求解域根据其几何特点划分成若干个子域, 在子域内构建径向基函数插值, 在界面上施加所有的界面连续条件,构建分块稀疏矩阵统一求解. 该方法仍然保持超收敛性, 且将原来的满阵转化成了稀疏矩阵, 降低了存储空间,提高了计算效率. 相比较于传统的径向基函数配点法和有限元法, 这种方法能够更好地反应局部特性和求解高梯度问题.数值分析表明该方法能够有效求解局部存在高梯度的大变形问题.      

基于分区径向基函数配点法的大变形分析

无网格法因为不需要划分网格,可以避免网格畸变问题,使得其广泛应用于大变形和一些复杂问题.径向基函数配点法是一种典型的强形式无网格法,这种方法具有完全不需要任何网格、求解过程简单、精度高、收敛性好以及易于扩展到高维空间等优点,但是由于其采用全域的形函数,在求解高梯度问题时存在精度较低和无法很好地反应局部特性的缺点.针对这个问题,本文引入分区径向基函数配点法来求解局部存在高梯度的大变形问题.基于完全拉格朗日格式,采用牛顿迭代法建立了分区径向基函数配点法在大变形分析中的增量求解模式.这种方法将求解域根据其几何特点划分成若干个子域,在子域内构建径向基函数插值,在界面上施加所有的界面连续条件,构建分块稀疏矩阵统一求解.该方法仍然保持超收敛性,且将原来的满阵转化成了稀疏矩阵,降低了存储空间,提高了计算效率.相比较于传统的径向基函数配点法和有限元法,这种方法能够更好地反应局部特性和求解高梯度问题.数值分析表明该方法能够有效求解局部存在高梯度的大变形问题.     

基于径向基函数配点法的梁板弯曲问题分析

采用径向基函数配点法分析考虑剪切效应的梁板弯曲问题,该方法利用径向基函数作为近似函数,基于配点法离散方程,通过最小二乘法求解。径向基函数配点法在离散和计算过程中不需要任何形式的网格划分,是一种真正的无网格法;径向基函数可以用一元函数来描述多元函数,存在明显的储存和运算简单的特点;而基于配点法求解不需要积分,提高了计算效率。分析考虑剪切效应的薄梁板问题时,传统的有限元法或无网格法求解均会存在剪切锁闭问题,而径向基函数在全域内存在无限连续性,能够准确地满足Kirchhoff约束条件,因此径向基函数配点法能够消除剪切锁闭现象,而且不会出现应力波动。该方法的优势在于,其不仅易于离散、精度高,而且具有指数收敛率,计算效率高。数值算例验证了上述结论和该方法的稳定性。     

无网格法求解一类分段连续型延迟偏微分方程

无网格法是基于散点信息求解偏微分方程问题的数值方法,无网格法可减少或完全消除对网格的依赖,数值实施更加灵活.因此,考虑采用基于径向基函数的无网格插值法求解一类分段连续型延迟偏微分方程.首先,利用θ-加权有限差分法得到方程时间上的离散格式,利用基于径向基函数的无网格插值法近似空间导数,得到了全离散数值格式.采用的基函数是Multiquadric (MQ)径向基函数,MQ径向基函数在精度及稳定性等方面都优于其他径向基函数.其次,采用傅里叶分析方法对该方法进行稳定性分析,得到了该方法稳定的条件,且该条件只与时间步长有关.最后,通过数值算例验证了方法的收敛性和稳定性,从而说明了方法的有效性和适用性.     

ELASTIC DYNAMIC ANALYSIS OF MODERATELY THICK PLATE USING MESHLESS LRPIM

A meshless local radial point interpolation method （LRPIM） for solving elastic dy-namic problems of moderately thick plates is presented in this paper. The discretized system equation of the plate is obtained using a locally weighted residual method. It uses a radial basis function （RBF） coupled with a polynomial basis function as a trial function,and uses the quartic spline function as a test function of the weighted residual method. The shape function has the properties of the Kronecker delta function,and no additional treatment is done to impose essen-tial boundary conditions. The Newmark method for solving the dynamic problem is adopted in computation. Effects of sizes of the quadrature sub-domain and influence domain on the dynamic properties are investigated. The numerical results show that the presented method can give quite accurate results for the elastic dynamic problem of the moderately thick plate.     

A smoothed Hermite radial point interpolation method for thin plate analysis

A smoothed Hermite radial point interpolation method using gradient smoothing operation is formulated for thin plate analysis. The radial basis functions augmented with polynomial basis are used to construct the shape functions that have the important Delta function property. The smoothed Galerkin weakform is adopted to discretize the governing partial differential equations, and a curvature smoothed operation is developed to relax the continuity requirement and achieve accurate bending solutions. The approximation based on both deflection and rotation variables make the proposed method very effective in enforcing the essential boundary conditions. The effects of different numbers of sub-smoothing-domains created based on the triangular background cell are investigated in detail. A number of numerical examples have been studied and the results show that the present method is very stable and accurate even for extremely irregular background cells.     

基于径向基点插值法的旋转Mindlin板高次刚柔耦合动力学模型

将无网格径向基点插值法(radial point interpolation method,RPIM)用于中心刚体?旋转柔性板的动力学分析.基于浮动坐标系方法和一阶剪切变形理论即Mindlin板理论,考虑剪切变形的影响,并计入板面内变形的非线性耦合变形项,采用径向基点插值法描述板的变形场,保留动能中有关非线性耦合变形项...     

四边简支矩形中厚板的弯曲

本文采用Reissner中厚板理论求解了四边简支矩形中厚板的弯曲问题。文中首先对Reissner中厚板理论的控制方程进行了适当的变更,使之成为非耦联的二阶偏微分方程组,然后利用有限积分变换法求解所得新的控制方程,得到了四边简支矩形中厚板受均布载荷作用下的解析解。文中所述方法可用以求解具有其它边界条件和载荷的矩形中厚板的弯曲问题,同时还可移植应用于其它中厚板理论。     

配点型广义节点无网格法的基本原理及其应用

借鉴流形方法思想,引入广义节点的概念,对传统的无网格法进行了改进,建立了可具有任意高阶多项式插值函数的广义节点无网格方法。同时采用径向插值函数构造具有插值特性的逼近函数;采用配点法建立系统的离散方程。在阐述了这种方法基本原理的同时,针对线弹性力学问题给出了这种方法的数值计算列式。与传统无网格方法相比,这种方法更具有一般性;同时由于采用了配点法而不需要背景积分网格,所以可以认为这种方法是某种真正意义上的无网格法。当选取0阶广义节点位移插值函数时便可得到传统的无网格法;在不增加支持域内节点数目的条件下,通过选取高阶广义节点位移插值函数可以提高计算精度。最后通过算例分析,对0阶、1阶及2阶广义节点无网格法与现有的有关解答进行了对比,论证了其合理性。     

高速公路沉降的非线性等阶径向点插值法解

无单元法一个突出的优点在于其只需要结点信息而不需要单元信息。先介绍等阶径向点插值法这种新型无单元的形函数构造思路，接着给出了它非线性求解平面比奥固结问题的主要方程，然后对一软基高速公路的断面沉降进行了计算，并与非线性有限元法结果进行了对比。可以看出该法不但计算精度高，而且在解路堤分级施工的这类移动边界问题的沉降时，比有限元法更方便，具有较好的应用前景。     

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

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