解二维抛物型方程初边值问题的交替方向隐式方法

详细地研究了各种情形的求解二维抛物型方程初边值问题的交替方向隐式差分法。这种方法可以将二维隐式方法归结为求解三对角线性方程组。与一维情况类似,可以继续利用追赶 法求解。它具有运算速度快,存储量小,无条件稳定等优点。该方法是求解二维物型方程的有效方法,必将得到更广泛的应用。     

冻粘土蠕变损伤耦合本构关系研究

以冻土试验为基础,提出了冻粘土服从DP屈服准则的粘塑损伤变量,通过推导得到了相关联流动法则下的冻土粘弹塑性损伤耦合本构方程.为了研究粘弹塑性和损伤的耦合作用,将建立的冻土材料粘弹塑性损伤耦合本构模型,通过用户子程序嵌入到ADINA有限元中.结果表明:粘弹塑损伤耦合本构模型能较好的模拟冻粘土加速蠕变过程,与实测蠕变吻合较好.     

新型边界积分方程的无单元解及其应用

目前隧道及大型地下工程往往在裂隙岩体中开挖,而裂隙与地下空间的距离及裂隙的扩展条件,制约着隧道及地下工程的稳定性。应用能考虑孔洞和裂纹问题的新型边界积分方程与无网格加辽金法结合,建立一种新型的边界无单元法。在该方法中基本的未知量是由边界上的面力和边界上位移密度函数构成的复变量边界函数 。文中应用的边界积分公式和Muskhelishvili的积分公式直接相关。将无网格构造方法引入新型的边界积分方程,建立了新型的边界无单元法。应用该方法详细分析了含隧道和裂纹间相互关系等问题,其数值结果与解析结果吻合很好,说明该方法的正确性和可行性。     

正交各向异性岩体孔径变形法的二维地应力分析方法改进

根据正交各项异性材料的弹性理论,推导出偏轴坐标系下的三维、二维本构方程;利用应力变分法对钻孔周围的应力分布作出修正,进一步得出孔壁径向位移的计算公式。结合孔径变形法测量地应力的特点,建立正交各向异性岩体地应力的分析理论,通过一典型工程与传统的各向同性理论进行了对比。     

TTI介质弹性波频率-空间域有限差分数值模拟

由周期性薄互层引起的VTI介质是研究比较广泛的一类各向异性介质。当VTI介质对称轴偏离垂向,本构坐标系与观测坐标系不重合时,会形成观测坐标系下的TTI介质。引入25点优化差分算子,推导出二维TTI介质频率域弹性波动方程;为压制边界反射,采用完全匹配层法吸收边界条件,并计算出优化差分系数;最后采用集中力源,模拟了弹性波在TTI介质中的传播过程。从波场快照和地面共炮记录可以看出,笔者采用的数值模拟算法能有效压制数值频散。TTI介质中的波场传播比较复杂,纵波传播相对稳定,横波波前的三分叉现象比较明显,并存在振幅奇异性。当VTI介质的对称轴偏转后,还会增加地面地震记录的复杂性。     

结合抛物型缓坡方程计算波浪辐射应力

将波浪辐射应力与抛物型缓坡方程中的待求变量联系起来,提出了一种计算辐射应力的新方法,并用有限差分法对控制方程进行了数值求解。数值结果表明这种方法精度高、编程简单、求解快速,可用于实际大区域波浪辐射应力的计算。     

海浪场非平衡信息论熵的研究初探

鉴于传统熵与时间无关 ,不能反映海 -气耦合机制下海浪场的基本属性———不可逆性 ,文章将非平衡信息论熵的概念引入海浪场的研究。依据海浪场相空间的Fokker Planck方程 ,首次推导出了非平衡信息论熵的熵平衡方程、熵密度、熵流密度、熵产生密度、广义力和广义力密度的具体表达式。     

二维各向异性介质中三分量波动方程有限元法模拟

利用六面体单元和三线性插值函数 ,推导出在任意弹性各向异性介质中三维三分量波动方程所满足的有限元方程 ,同时给出也适用于矩形单元和双线性插值的二维三分量各向异性介质中的波动方程模拟的一般性公式 ,最后讨论两个数值模拟的结果。     

求解代数Riccati方程的一些数值方法评述

在过去的三十多年中 ,代数 Riccati方程的数值求解一直是数值代数界和控制论界关注的重要课题。作者简要评述几种用于求解代数 Riccati方程的有代表性的数值方法。     

波浪导致的海床边坡稳定性分析

利用Biot二维固结理论,计算出波浪引起的海床应力场;在假设海床土体的粘聚力不受波浪作用影响的条件下,推导出了波浪导致弱粘性土或粉土海床失稳的破坏准则,丰富了波浪导致海床失稳的判别方法。并把这一破坏准则用于黄河三角洲埕岛油田海域的海床底坡稳定性分析,讨论了50年一遇和5年一遇两种海况条件下,土体粘聚力对该地区波浪导致海床破坏深度的影响,计算结果表明,土体粘聚力对海床稳定性影响明显。