1. 引言
微分求积法(DQM) [1] 具有计算逻辑简明、程序过程紧凑、计算量少、计算精度高等优点，近年来已经被广泛应用于众多研究领域 [2] [3] 。但是，DQM在求解具有局域特殊性问题时，其计算收敛性和准确性受到很大影响，为了改善这一问题，一些研究者将DQM与势能原理相结合 [4] ，使得求解过程中同时存在数值微分和积分。对于势能泛函中的微分和积分计算，现有做法是利用微分求积法离散微分形式，利用Gauss-Lobatto求积、积分求积法 [5] 及其他求积方法来离散积分形式 [6] [7] 。由于微分和积分需要分开离散，且易受微积分计算节点不统一的约束，导致计算流程较为繁琐，计算效率降低。有鉴于此，本文提出了基于DQM和积分求积法(Integral Quadrature Method，简称IQM)思想的统一化的微分–积分求积法(Differential-Integral Quadrature Method，简称D-IQM)，建立计算微分求积法和积分求积法的权系数的统一方式，进一步提高了势能原理框架下的新型数值方法D-IQM的计算效率，从而有利于解决局域特殊性问题。
2. 原理简介
2.1. 微分求积法
微分求积法(DQM)的基本原理：函数
在区间[a,b]上连续可微，则其在给定节点处的导数值可以表示为域内全部离散节点的函数值的加权和，即
(1)
其中，
是插值基函数，函数的高阶导数可表示为：
(2)
式中，
为k阶导数加权系数。
若基函数采用拉格朗日插值多项式，则一阶导数权系数表达式如下：
(3)
其中，
。已知一阶权系数，各高阶权系数
可由下式推出：
(4)
2.2. 统一化的微分–积分求积法(D-IQM)
下面，基于DQM和积分求积法(Integral Quadrature Method，简称IQM)思想，本文提出统一化的微分–积分求积法(Differential-Integral Quadrature Method，简称D-IQM)，建立计算微分求积法和积分求积法的权系数的统一方式。
设定在区间[a,b]上连续可微的函数
、
，对
进行n阶求导，有
(5)
式(5)结合DQM，且设
，
，则
(6)
一般的，设
(7)
则式(6)可以表达为如下矩阵形式：
(8)
(9)
的一次积分可以表达 [2] 为：
(10)
其中，c为积分区域外任意一点，由上述关系可得：
(11)
(12)
式(10)可进一步表示为：
(13)
一般的，设
(14)
式(13)的矩阵表达式如下
(15)
因此，可以利用微积分统一形式的权系数计算方法解决微分计算和积分计算。
对于上述方法中n = 1的特殊情况，参看式(5)可知
的积分为原函数
，设插值基函数
，
为插值函数的系数矩阵，则有
(16)
因此，由式(13)可知
(17)
结合式(9)，
(18)
式中
为统一化的集成权系数。
至此，我们得到了一套计算微分求积法和积分求积法的权系数的统一方式，即根据微分求积法方式所确定的权系数，在求解积分问题时可以采用式(15)进行计算；特别是当n = 1时可以采用式(18)进行计算。上述方式有助于提高新型数值方法D-IQM的计算效率，从而进一步拓展微分求积法的应用范围。
2.3. 基于势能原理的微分–积分求积法(D-IQM)
势能原理以变分形式表示物理定律 [8] ，即在满足一定约束条件的所有可能物体运动状态中，真实的运动状态使某物理量(如势能泛函)取极值或驻值。对于复杂结构，特别是具有局部特殊性的问题，引进势能泛函，应用D-IQM进行离散计算，有助于求解这类问题，改善计算结果的收敛性与准确性。
一般来讲，一维梁结构的总势能泛函表达式 [9] 如下：
(19)
根据D-IQM，上式离散成如下表达形式：
(20)
定义下列矩阵及向量：
(21)
那么势能泛函式（20）可以写成矩阵表达形式：
(22)
根据最小势能原理，结合式（21）中定义的权系数矩阵有
(23)
其中，在一维问题中，梁的刚度矩阵K为
(24)
梁的载荷矩阵Q为
(25)
在微分–积分求积法中，仍然可以采用方程替代法将边界条件加入到刚度矩阵中去，在此不再详述。
无论一维或是二维的局域性问题，结合势能原理的微分求积法具有良好的计算稳定性、收敛性及准确性，这主要是因为在势能原理的框架下，不仅方程中需要求导的阶次降低，而且泛函的形式更有利于局域载荷在求解中的处理，这些都有利于DQM更好地解决局域特殊性问题。
3. 实例结果及分析
3.1. 简单函数
为了检验上述微分–积分求积法(D-IQM)的正确性，现以若干简单函数(参见表1)为例进行对比。
所求结果与精确解对比如图1所示，可以发现以D-IQM所求结果与精确解近似，并随着节点的增多，呈现出较好的准确性、收敛性，计算结果有着良好的稳定性。
3.2. 力学问题
本文提出的D-IQM沿承了传统DQM的计算逻辑，在势能原理的框架下同样可以更好的处理应用性问题。为了进一步验证D-IQM的应用性，将其与势能原理结合(参见2.3节有关推导)求解具有局域特殊性的实际问题，下面以梁为例。
如图2所示，梁的材料弹性模量为
，泊松比
，梁长
，梁横截面高
，宽
，坐标轴x以梁的左端为原点，梁的载荷类型及边界条件如表2所示，其中q为分布于全梁的均布载荷。
表2中“精确解”和“本文解”均指梁的最大挠度，其中，“本文解”即采用本文提出的D-IQM得到的解。从表2可以看出，上述两种解对于表2所列各种梁问题都取得了较为一致的结果。图3显示了
Table 1. Integrable function and integral interval
表1. 被积函数及积分区间
Table 2. Load type and boundary condition of beam
表2. 梁的载荷类型及边界条件
Figure 2. Simply supported beam of single concentrated force
图2. 单集中力作用下的简支梁
Figure 3. Curve: maximum deflection of beam
图3. 梁最大挠度收敛情况
本文D-IQM求解的收敛趋势，由该图可知，D-IQM解随着节点数增多逐渐收敛于有限元值，并保持了良好的稳定性。这些结果表明，统一化的微分–积分求积法在保持原有DQM计算逻辑的简明性前提下，具有良好的数值计算性能，以及良好的应用性；同时，微分权系数计算方法和积分求积法计算方法的统一有助于简化计算逻辑，提高计算效率。
4. 结论
本文基于DQM和积分求积法，本文提出统一化的微分–积分求积法(简称D-IQM)，建立了计算微分求积法和积分求积法的权系数的统一方式，根据微分求积法的权系数计算积分权系数，从而简化了计算逻辑。文中对基础函数积分进行计算，通过对比讨论证明了该方法在收敛速度、求解精度、稳定性、应用范围方面的良好性能。进一步通过力学问题实例的求解，证明新型数值方法D-IQM可以有效提高势能原理框架下相关问题的计算效率，从而有利于解决局域特殊性问题。上述应用实例从不同侧面说明了本文提出的D-IQM的应用性和有效性。