集中载荷作用下梯度复合材料梁的小波-微分求积法分析

将梯度复合材料梁作为平面应力问题处理,采用小波和微分求积混合法,对集中荷载作用下结构的响应进行了分析.考虑材料特性参数沿高度方向呈梯度分布,在该方向上采用广义微分求积法进行离散;鉴于广义微分求积法求解集中荷载问题精度不高的缺点,在梁的长度方向上引入对突变信号敏感的小波插值函数.数值计算表明,小波-微分求积混合法不仅保留了广义微分求积法高效的优点,而且能够很好地模拟结构局部化特征.     

微分求积单元法在结构工程中的应用

微分求积法（Differential Quadrature Method）是求鳃偏微分方程和积分-微分方程的一种数值方法，该法具有计算简便、精度较高和易于实现等优点。微分求积单元法（Differential Quadrature Element Method）是在微分求积法的基础上结合区域分割和集成规则而形成的一种新的数值计算方法，能通过自适应地选取微分求积网点数目正确模拟构件的刚度和荷载性质，其精度可通过细分单元或增加离散点数目加以提高。微分求积单元法是一种可供选择的、性能优越的数值计算方法。本文将详细论述这一数值方法的基本原理，并通过数值算例说明该方法的应用过程及其优越性，为这一方法在结构工程中的推广应用提供参考。     

微分求积法的特性及其改进

微分求积法已在科学和工程计算中得到了广泛应用。然而,有关时域微分求积法的数值稳定性、计算精度即阶数等基本特性,仍缺乏系统性的分析结论。依据微分求积法的基本原理,推导证明了微分求积法的权系数矩阵满足V-变换这一重要特性;利用微分求积法和隐式Runge-Kutta法的等值性,证明了时域微分求积法是A-稳定、s级s阶的数值方法。在此基础上,为进一步提高传统微分求积法的计算精度,利用待定系数法和Padé逼近,推导出了一类新的s级2s阶的微分求积法。数值计算对比结果验证了所提出的新微分求积法比传统的微分求积法具有更高的计算精度。     

混合微分求积法在功能梯度材料平板柱形弯曲问题中的应用

利用混合微分求积法,对任意荷载作用下不同材料梯度分布的功能梯度材料平板柱形弯曲问题进行了分析。针对广义微分求积法求解集中荷载问题精度不高的缺点,本文利用小波微分求积法进行了改进。由于小波对突变信号具有良好的自适应描述能力,因此在平板宽度方向上,利用小波微分求积法可以有效地处理集中荷载;而在材料梯度变化的板厚方向上,则利用广义微分求积法计算量小且精度高的特点进行离散计算。计算表明,混合微分求积法不仅保留了广义微分求积法高效的特点,而且能有效地求解任意荷载作用的问题。通过算例,分析了在机械荷载作用下,材料不同梯度形式、平板上下表面材料性质差异对功能梯度平板结构响应的影响。     

框架结构P-△效应分析的微分求积单元法

采用一种新的数值方法——微分求积单元法分析框架结构的P-△效应。微分求积单元法采用微分求积法直接求解微分方程的技术，并结合有限分割技术而形成。首先建立考虑剪切变形和轴力二阶效应的框架结构单元平衡微分方程，通过微分求积离散而得到梁单元的一般弹性刚度方程；同时考虑变形后节点的平衡条件和变形协调条件，导出框架结构整体二阶分析的微分求积单元法力学模型。由于该分析模型中包括了单元及结构的所有离散形式的控制方程，因此采用该模型进行结构分析可得出较为精确的解。数值算例的分析比较，表明了该法用于框架结构P-△效应分析的正确性和有效性。本文导出的框架结构二阶分析的微分求积单元法力学模型可用于框架结构剪切变形与几何非线性的耦合效应分析。     

拟谱-微分求积混合方法求解一类双曲电报方程

拟谱方法和微分求积法是两类重要的无网格法，二者都已在科学和工程计算中获得了广泛应用。采用拉格朗日插值多项式作为二者的试函数，且采用同一种网格点分布，指出了在空间域上，微分求积法是拟谱方法的一种特殊形式。在此基础上，结合二者各自的特点，提出了拟谱-微分求积混合方法用于求解一类双曲电报方程。理论分析和数值测试表明，新方法在空间域上具有谱精度收敛性，在时间域上是A-稳定的，比较适合于求解多维电报方程。     

基于小波微分求积法的薄板弯曲分析

利用小波微分求积法(WDQM)对任意荷载作用下的薄板弯曲问题进行了求解分析。数值算例表明,小波微分求积法与一般的DQ法相比具有很好的适用性,特别是薄板受集中荷载或不连续分布荷载作用时,由于小波基函数的紧支撑特性与其对突变信号良好的描述能力,WDQ法的精度明显优于一般的DQ法,具有良好的应用前景。     

基于边界值方法的微分动力系统数值计算方法

对高维非线性初值问题,微分求积法在每一步的积分过程中需要求解一个更高维的非线性方程组,因而计算量巨大。基于微分求积法与边界值方法两者之间的关系,可以将广义向后差分方法和扩展的隐式梯形积分方法看作是经典微分求积法的稀疏表达形式。将广义向后差分方法以及扩展的隐式梯形积分方法这两类边界值方法应用于微分动力系统的数值计算,提出了一类新的数值计算方法。理论分析及算例结果表明,对高维非线性微分初值问题的数值计算,本文方法相对于经典的微分求积法具有更高的计算效率。     

弹性地基上变厚度矩形板自由振动的GDQ法求解

基于广义微分求积法(GDQ法),对弹性地基上变厚度矩形板横向自由振动的控制微分方程及其不同边界条件进行离散,研究了其自由振动的频率特性。数值计算得到了不同长宽比?、不同厚度变化参数?、不同地基参数K条件下以及简支或固定边界条件下弹性地基上变厚度矩形板的量纲为一的振动频率,并与已有文献进行了比较。结果表明:运用广义微分求积法对弹性地基上变厚度矩形板的频率求解结果在退化到K=0时与幂级数解的结果非常吻合;在条件相同的情况下,采用广义微分求积法仅需较少的节点(N=M=13)就能达到满意的求解精度。本文的研究为求解此类问题的低阶、高阶振动频率提供了一种简便有效的数值方法。     

???????????????????????е????

叙述了微分求积法的一般方法，并把该方法应用于求解单摆的非线性振动的周期问题，取得了比常用的L-P法更满意的结果．     

Homotopy analysis and differential quadrature solution of the problem of free-convective magnetohydrodynamic flow over a stretching sheet with the Hall effect and mass transfer taken into account

This paper presents an analytical solution of the problem of free-convective magnetohydrodynamic flow over a stretched sheet with the Hall effect and mass transfer taken into account. A similarity transform reduces the Navier-Stokes, energy, Ohm law, and mass-transfer equations to a system of nonlinear ordinary differential equations. The governing equations are solved analytically using an analytical method for solving nonlinear problems, namely, the homotopy analysis method. The results are compared with the results of a promising numerical method of differential quadrature developed by the authors. It is shown that there is very good agreement between analytical results and those obtained by the differential quadrature method. The differential quadrature method was validated, and the effects of non-dimensional parameters on the velocity, temperature and concentration profiles were studied.     

HOPF BIFURCATION OF A NONLINEAR RESTRAINED CURVED PIPE CONVEYING FLUID BY DIFFERENTIAL QUADRATURE METHOD

This paper proposes a new method for investigating the Hopf bifurcation of a curved pipe conveying fluid with nonlinear spring support. The nonlinear equation of motion is derived by forces equilibrium on microelement of the system under consideration. The spatial coordinate of the system is discretized by the differential quadrature method and then the dynamic equation is solved by the Newton-Raphson method. The numerical solutions show that the inner fluid velocity of the Hopf bifurcation point of the curved pipe varies with different values of the parameter,nonlinear spring stiffness. Based on this, the cycle and divergent motions are both found to exist at specific fluid flow velocities with a given value of the nonlinear spring stiffness. The results are useful for further study of the nonlinear dynamic mechanism of the curved fluid conveying pipe.     

高速小展弦比机翼颤振的微分求积法分析

引入微分求积法，分析高速小展弦比机翼的气动弹性问题。将小展弦比机翼等效为悬臂板，基于一阶活塞气动力理论建立机翼颤振偏微分方程，采用微分求积法将偏微分方程转化为常微分方程，根据频率重合理论对颤振问题进行求解。分析了机翼的固有频率及颤振速度，并与有限元软件计算结果进行比较，误差在2％以内，很好的验证了微分求积法求解小展弦比机翼颤振问题的有效性。分析了机翼面积、展弦比及厚度对颤振速度的影响，结果表明，小展弦比机翼的颤振速度受结构尺寸的影响较大，颤振速度随面积和展弦比的增大而减小，随机翼厚度的增大而增大。     

In-plane vibration analysis of plates in curvilinear domains by a differential quadrature hierarchical finite element method

Free in-plane vibration analysis of plates is carried out by a differential quadrature hierarchical finite element method (DQHFEM). The NURBS (Non-Uniform Rational B-Splines) patches of geometries were first transformed into differential quadrature hierarchical (DQH) patches, and then the elastic field was discretized by the same DQH basis. The DQHFEM solved the compatibility problem caused by different parametrization of neighbouring patches of isogeometric analysis using NURBS. And mesh refinement in DQHFEM does not propagate from patch to patch. The DQHFEM matrices also have the embedding property as the hierarchical finite element method (HFEM). In-plane vibration analyses of plates of several planforms showed that the DQHFEM is similar as the fixed interface mode synthesis method that can analyse a structure using a few nodes on the boundary of substructure elements and only several clamped modes inside each substructure element, but the DQHFEM does not need modal analysis and is of high accuracy. The accuracy and convergence of the DQHFEM were validated through comparison with exact and approximate results in literatures and computed by the authors.     

NONLINEAR RESPONSES OF A FLUID-CONVEYING PIPE EMBEDDED IN NONLINEAR ELASTIC FOUNDATIONS

The nonlinear responses of planar motions of a fluid-conveying pipe embedded in nonlinear elastic foundations are investigated via the differential quadrature method discretization （DQMD） of the governing partial differential equation. For the analytical model, the effect of the nonlinear elastic foundation is modeled by a nonlinear restraining force. By using an iterative algorithm, a set of ordinary differential dynamical equations derived from the equation of motion of the system are solved numerically and then the bifurcations are analyzed. The numerical results, in which the existence of chaos is demonstrated, are presented in the form of phase portraits of the oscillations. The intermittency transition to chaos has been found to arise.     

压电-弹性层合梁自由振动分析的新方法

将状态空间法和微分求积法相结合，分析了压电-弹性层合梁的自由振动. 通过微分求积把状态方程在每一个节点处离散，进而求得解答. 选用不同的节点数目， 分析了方法的收敛性. 计算结果与相关文献的结果能较好地符合. 该方法 对于分析压电-弹性层合梁的工程振动问题非常方便.     

Edge crack opening in an elastic plane with a wedge-like cut in contact with a punch

The equilibrium of an elastic plane with a wedge-like cut and an internal or edge crack on the symmetry axis was studied in [1] in the case of punch indentation in the lateral faces of the cut at a distance from the cut tip. In [1], the systemof singular integral equations of the problemwas solved numerically by the mechanical quadrature method. In this paper, the generalized Wiener-Hopf method [2] is used to obtain the analytic solution of a similar problem in the case of an edge crack under punch pressure on parts of the cut lateral faces adjacent to the cut tip. Some special cases of this problem were considered earlier without a crack [3, 4] or a punch [5, 6].     

ANALYSIS OF ELASTO-PLASTIC POSTBUCKLING AND ENERGY RELEASE RATE FOR DELAMINATED FIBER METAL LAMINATED BEAMS IN THERMAL ENVIRONMENT

The elasto-plastic postbuckling of fiber metal laminated beams with delamination and the energy release rate along the delamination front are discussed in this paper. Considering geometrical nonlinearity, thermal environment and geometrical initial imperfection, the incremental nonlinear equilibrium equations of delaminated fiber metal laminated beams are established,which are solved using the differential quadrature method and iterative method. Based on these,according to the J-integral theory, the elasto-plastic energy release rate is studied. The effects of some important parameters on the elasto-plastic postbuckling behavior and energy release rate of the aramid reinforced aluminum laminated beams are discussed in details.