首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
本给出了用边界元法积分方程的算法。通过处例验算表明,这种算法综合了边界元法精度高和用矩最法解积分方程速度快的优点,是一种求解积分方程十分有效的算法。  相似文献   

2.
本文给出了用边界元法求积分方程的算法。通过算例验算表明,这种算法综合了边界元法精度高和用矩量法解积分方程速度快的优点,是一种求解积分方程十分有效的算法。  相似文献   

3.
吴大志  张振营 《岩土力学》2015,36(12):3393-3399
考虑到地基在竖直方向上的非均匀性,结合扭转振动的特点,建立了简谐扭转动荷载作用时剪切模量随深度增大的广义Gibson饱和地基的动力方程,通过积分变换求解了动力方程。考虑到半空间地基表面处应力自由、埋置圆板所在平面为混合边界和无穷深度处为波的辐射边界等条件,得到了广义Gibson饱和地基中刚性圆板扭转振动时的对偶积分方程,通过合适的变换转化成了一个第2类Fredholm积分方程,求解了相应的动力响应问题。对比静扭距作用时的荷载-位移关系,给出了动力柔度系数和扭转角位移幅值的表达式,并把所研究的问题进行退化且与前人成果进行了对比。数值研究表明:当基础的埋置深度小于5倍基础半径时,广义Gibson饱和地基中埋置基础的扭转振动存在明显的边界层现象,且埋置深度越小,边界层现象越明显。  相似文献   

4.
针对图解法求取孔隙结构特征值使用条件严格,同时误差较大,而矩法饱和度范围又受限于试验。采用积分法求取特征值,即对孔隙半径进行积分求取特征值。积分法必须先获取孔隙半径数学表达式,即毛管压力数学模型。对毛管压力数学模型进行了探究,在双对数坐标系下推导并建立了毛管压力幂函数模型,采用最小二乘法求取模型参数,并通过对比Sami论文数据验证了推导的模型。深入分析矩法与积分法的联系可知,代入孔隙半径数学表达式即可用积分法求取岩石孔隙结构特征值。最后,对3种方法求取结果进行分析讨论,得出积分法算法严谨,更加接近真实值。这一研究加深了对岩石微观孔隙结构认识。  相似文献   

5.
讨论了平面上广义解析函数的非线性Riemann问题,把它转化为奇异积分方程,并利用压缩原理证明解的存在性。  相似文献   

6.
系数矩阵存储和线性方程组求解是限制三维电磁积分方程方法发展的主要因素。Zhdanov提出准线性(QL)近似技术,建立了复杂散射场与背景场的线性关系,有效地避免了积分方程中大型线性方程组的求解,但是该算法用于多源问题航空电磁正演模拟时精度不高。因此,本文提出一种基于多重网格准线性(MGQL)近似的算法,并利用系数矩阵的Toeplitz性质存储矩阵和快速傅里叶变换,实现了矩阵与向量的快速乘积、降低了计算复杂度,采用多重网格结合了积分方程方法和准线性近似解法的优点,在保证精度的条件下提高计算速度、减少存储量。针对不同类型网格的模拟实验表明,相比于传统积分方程方法,本文算法在保证计算精度的同时,可以将计算速度极大地提高(>10倍)。  相似文献   

7.
利用传输线理论、Sommerfeld积分快速计算以及最小平方拟合技术研究建立多分量感应测井数据的一种新的快速参数化迭代反演算法,同时重构水平层状横向同性地层的纵、横向电阻率以及水平层界面深度。首先,通过Fourier变换与传输线理论给出频率波数域中电磁场并矢Green函数在各个地层中的解析解,并利用三次样条插值和贝塞尔函数递推公式建立Sommerfeld积分的半解析算法,快速计算多分量感应的测井响应。然后在此基础上,利用摄动理论建立磁场并矢Green函数与模型向量间变化关系的摄动方程,并将摄动方程中各个积分转化为Sommerfeld积分,实现正演模拟的同时用半解析算法快速确定多分量感应测井响应的Fréchet导数。最后,利用归一化处理和奇异值分解技术,同时反演所有地层的纵、横向电阻率和层界面深度,实现输入数据和反演模型的模拟数据优化拟合。理论模型的数值结果验证了该反演算法的有效性及抗噪性。  相似文献   

8.
二重网格方法用于曲面位场转换   总被引:1,自引:0,他引:1  
本文介绍了二重网格方法及在求解曲面位场转换的积分方程中的应用。给出该积分方程的二重网格方法的具体算法和二重网格迭代算子的结构图。并提出了三点注意事项。  相似文献   

9.
利用格林公式及点源场的边界条件,将地下半空间位函数的体积积分变为地面位函教的积分方程以后,用边界元方法解此积分方程便可得到地面各节点的电位值。其计算精度和速度取决于剖分方式和节点密度。本文采用变网格的剖分方式,用较少的点实现了计算,并兼顾了精度和速度。  相似文献   

10.
本文利用脉冲技术研究了地震勘探反问题的算法。同时还应用积分方程的解核表出示了正演问题的解。从而得到了点距、频率、反射层深度之间的级数关系,对于反演问所导出的第一类积分方程,文中指出了应用作者在位场反演中所采用的方法,这就保证了反问题解的稳定性。  相似文献   

11.
高源  刘海笑  李洲 《岩土力学》2019,(10):3951-3958
基于适用于饱和砂土循环动力分析的边界面模型,利用修正广义Mises方法将其推广至三维应力空间中,引入与状态相关的剪胀函数,采用历史最大加载面硬化准则,以反映饱和砂土排水条件下刚度变化、剪胀剪缩等特性。结合子增量步显式积分算法的思想,建立了适用于饱和砂土循环动力分析边界面模型的显式积分算法流程。借助有限元软件子程序接口,将该算法流程开发到有限元软件中,通过建立土体单元数值模型,对Toyoura砂在单调荷载和循环荷载下的三轴试验工况进行模拟。结果表明,利用子增量步显式积分算法对模型进行积分,能够准确有效地模拟砂土的应力-应变曲线、以及砂土在单调和循环荷载作用下均展现的剪胀剪缩现象,验证了显式积分算法的有效性。通过设置不同的增量步长,验证了子增量步显式积分算法对于较小应变增量步的高度稳定性和收敛性。  相似文献   

12.
韩泽军  林皋  周小文  杨林青 《岩土力学》2018,39(6):2287-2294
动力响应问题的求解对于地基在外荷载作用下引起的弹性波动问题研究有重要的意义。本文提出了一种求解横观各向同性层状地基在施加时间简谐荷载作用下任意点的应力响应的算法。此算法利用傅里叶变换将广义平面应变问题频率-空间域的动力方程转化到频率-波数域内,结合对偶变量的引入,利用高精度的精细积分算法对状态方程进行求解,在得到频率-波数域内的位移响应的基础上,利用傅里叶逆变换得到任意点的动应力响应。简谐荷载不仅可以施加在地基表面,而且可以施加在地基内部。对比算例验证了本文算法的准确性,同时对地基各向异性特性、激励频率和阻尼比对动应力响应的影响进行了参数分析,为工程实际提供可靠的数值依据。  相似文献   

13.
一类断层效应反演问题的条件稳定性和数值算法   总被引:2,自引:2,他引:2  
考虑一类垂直断层效应反演问题,其数学模型是具有解析核的第一类积分方程。因其不适定性,获得它的稳定解非常困难。为获得其数值解,有必要先揭示它的条件稳定性。由积分方程方法、解析延拓方法和积分变换性质,获得了它的条件稳定性。基于条件稳定性设计了一种正则化算法,获得了近似解的误差估计。  相似文献   

14.
肖洪天  王小敏 《岩土力学》2004,25(Z1):91-94
在边界积分方程离散的基础上,根据远场位移的特性,用衰减函数描述远场位移,阐述了自适应高斯积分算法,并结合前人成果—用解析法求出在远场上的强奇异积分,编制了适合于半无限域线弹性问题的计算机程序,并进行了事例验证,验证表明,该法比传统边界元具有更高的精度。  相似文献   

15.
为提高电测深曲线反演过程的稳定性和自动性,基于广义交叉验证自适应正则化方法,给出了电测深曲线反演的全自动迭代方法.首先介绍了引入光滑约束的线性反演方程的构建方法,综合奇异值分解法与广义交叉验证自适应正则化方法求解线性反演方程,接着给出了反演过程中初始模型的构建方法和反演偏导数矩阵的差分计算方法,并采用C++语言实现了电...  相似文献   

16.
本文介绍识别多类异常的句法模式识别方法。它将异常值分段(或分小区间)并提取特征,根据特征值确定基元,然后将异常表示成链。用已知异常的链设计最近邻分类器,通过它对未知异常进行分类和识别。文章中仅对异常识别给出了特征值的提取,链的形成以及求不同长度链之间距离的算法。  相似文献   

17.
从目标函数的线性化出发,采用奇异值分解法求广义逆矩阵并适当加阻尼,用迭代法直接解超定方程实现改进广义逆矩阵反演。进一步从理论与实际结合上阐明改进广义逆矩阵反演法的灵活、稳定,它涵盖了多种反演方法的优点,又具有自身特点,在实际应中适应性强,并能提供一些辅助信息,可更好地评价解释结果。   相似文献   

18.
修正剑桥模型的隐式积分算法在ABAQUS中的数值实施   总被引:5,自引:0,他引:5  
范庆来  栾茂田  杨庆 《岩土力学》2008,29(1):269-273
利用大型有限元软件ABAQUS所提供的用户材料子程序UMAT接口,针对修正的剑桥本构模型开发了隐式积分算法,并且与自动选择时间步长的增量有限元方程迭代解法相结合,对正常固结土与超固结土的三轴排水与不排水试验进行了数值模拟。结果表明,所发展的隐式本构积分算法与时间步长自动选择方法具有较好的稳定性和较高的计算精度,能够得到比较合理的数值分析结果。  相似文献   

19.
第三代相干体算法(C3算法)具有分辨率高、压制噪声能力强的优点,但是该算法需要计算协方差矩阵的特征值,所以耗时较多。为提高大矩阵特征值分解的运算速度,提出了一种估计最大特征值的快速收敛算法,通过运用A2k矩阵的迹来构造一种收敛速度更快的特征值估计方法,并给出了这种方法的误差估计表达式,实现计算精度的有效控制。文中还讨论了地震数据体边缘补地震道方法,解决了将递推算法应用到平面方向时空间窗越过边界的问题。实例表明,本算法的计算效率明显高于常规的C3算法。  相似文献   

20.
常规的逆时偏移算法都是基于双程波动方程的数值解法,实际上利用Kirchhoff积分进行波场延拓,也可以实现逆时偏移。以Kirchhoff积分为基础,采用高斯束叠加积分表示的格林函数来表示正反向延拓波场,将正反向延拓波场进行保幅的互相关得到地下成像。保幅高斯束逆时偏移方法结合了射线类偏移方法的高计算效率和逆时波动方程偏移的高精度,能很好地处理焦散点、陡倾角等问题,并且具有面向目标成像的能力,并消除了偏移过程中的几何扩散效应和入射角变化引起的振幅误差。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号