双重贝寒尔函数积分的数值计算   总被引：2，自引：1，他引：2  

华军  蒋延生  汪文秉 《煤田地质与勘探》2001,29(3):58-62

双重贝塞尔函数积分则由于贝塞尔函数互乘项的强振荡与慢衰减特性而以以应用通常的数值积分算法。本文将被积区间[0,∞]划分为[0，λ0]、[λ0，∞]两部分，应用贝塞尔函数的克尔函数表述式及后者的大宗量渐近特性，区间[λ0，∞]的双重贝寒尔函数积分可被转化为Fourier正（余）弦变换，并可利用各种快速算法对其进行数值计算；区间[0，λ0]上的双重贝塞尔函数积分的计算可直接应用一般的数值积分算法并能获得较高的计算精度；当需大量计算有共同参量的双重贝塞 尔函数积分时，其计算效率仍显不足。此时，可应用贝塞尔函数的导数关系式对[0,λ0]内的双重贝塞尔函数积分进行恒等变换，再用差商近似导数，将其转化为对贝塞尔函数本身的积分，而该积分又仅需计算一次。故本算法对双重贝塞尔函数积分的计算效率有明显提高。  相似文献   

海上平台的船舶静力性能的精确计算   总被引：1，自引：0，他引：1  

周超骏  曹沅 《海洋工程》1996,(2)

本文提出了用双三次贝齐尔（Ｂｅｚｉｅｒ）曲面拟合水上水下工程结构物，包括海上平台、深潜器和各种类型的船舶。用曲面方程，精确计算船舶的静力性能，然后同传统近似的计算静力性能的数值作一比较，指出精确与近似计算结果的误差范围，供设计人员参考应用。  相似文献   

Precise integration method without inverse matrix calculation for structural dynamic equations     

汪梦甫 《地震工程与工程振动(英文版)》2007,6(1):57-64

The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.  相似文献   

Theis解析模型的数值解及其在MatLab中的实现     

李林子  朱国荣  江思珉  闵望 《江苏地质》2007,31(3):242-246

讨论了如何应用数值积分方法求解含水层中地下水非稳定流抽水试验的解析模型以获得含水层的水文地质参数,并以Theis公式为例讨论求解原理。引入了MatLab的数值积分工具和优化工具,轻松实现了对水文地质参数的求解,此举既保证了解的可靠性和唯一性,也大大提高了方法的可操作性。最后以一个具体的承压水非稳定流抽水试验为例验证了这个方法,所得的结果与应用配线法得到的结果基本一致,从而肯定了方法的应用价值。研究结果尽管是针对Theis模型进行的,但求解原理同样适用于其他各种模型。  相似文献   

基于节点计算的自适应数值积分及其程序实现     

孙海涛  王元汉 《岩土力学》2007,28(5):995-1000

数值积分是基于节点计算的数值分析方法的关键技术之一。针对完全无网格法的计算要求，提出了一种能自动适应任意计算域上各种节点分布方式的数值积分算法。该算法能随计算点的位置不同，自动确定积分域及积分域内的求积点，计算出数值积分结果，给出了临空边坡滑移体和圆孔方板2个实例。结果表明，按文中算法编制的程序具有良好的边界再生能力和较高的计算精度，可作为基于节点数值分析方法的一个程序模块，用于实际计算。底层算法思想，使得程序可方便地转换为其他计算语言。  相似文献   

长偏移距瞬变电磁测深甚晚期响应及视电阻率的数值计算   总被引：2，自引：0，他引：2       下载免费PDF全文

翁爱华  王雪秋 《地震地质》2003,25(4):664-670

文中讨论了提高长偏移距瞬变电磁测深甚晚期响应计算精度及相应的全区视电阻率的计算方法。在频率域 ,利用直接数值积分结合连分式展开提高低频段响应的计算精度 ,从而提高时间域甚晚期响应的计算精度。同时 ,利用连分式渐进展开 ,将均匀半空间电阻率表示为该模型长偏移距瞬变电磁测深响应的反函数 ,利用迭代方法求出适合全时间段的全区视电阻率。数值结果表明 ,文中的方法可有效地提高长偏移距瞬变电磁测深 (甚 )晚期数据的模拟与解释精度  相似文献   

饱和土表面在水平集中荷载作用下的瞬态反应   总被引：3，自引：0，他引：3  

周香莲  王建华  陆建飞 《岩土力学》2002,23(3):387-391

分析了水平集中荷载作用在半空间饱和土表面时的瞬态问题。利用Laplace-Hankel变换对非轴对称Biot固结方程进行解耦，利用Laplace-Hankel数值逆变换得到半空间饱和土在时域内的数值解。退化到线弹性中的解与文献中的结果进行比较，验证了文中结果的正确性和数值逆变换的可靠性.可以用于研究地震工程中地震波的传播以及饱和土与结构之间相互作用等问题。  相似文献   

无穷元刚度矩阵及其数值积分     

张伯艳  陈厚群 《地震工程与工程振动》1992,12(3):92-99

本文将二维映射动力无穷元进行了合理分类,并讨论了各类无穷元刚度矩阵计算中的敛散性,推导出更有效、更适宜的计算无穷元刚度矩阵的数值积分公式。  相似文献   

广义索末菲积分的离散复镜像算法   总被引：1，自引：0，他引：1       下载免费PDF全文

王建 《地球物理学报》1992,35(5):644-649

本文对多分层媒质并矢格林函数中的索末菲积分,用Prony方法对其被积函数分式部分用指数函数逼近.利用索末菲恒等式,无穷积分就可化为闭合的级数形式.该闭合形式可看作一串离散镜像在观察点的解.所得结果与数值积分结果吻合很好.  相似文献   

应用数值积分进行重力地形改正的F77程序设计     

向运川 《青海地质》1994,3(1):48-52

  相似文献