有自由面渗流的无单元法

为解决传统有限元方法很难处理复杂渗透介质的计算域这一问题,提出了有自由面渗流的无单元法。无单元法利用滑动最小二乘法来建立在全域高阶连续可导的插值函数,具有积分网格和节点相互独立的优点,可以避免有限元法中网格在迭代过程中变化的问题,实现了真正意义上的网格固定。计算实例表明,无单元法计算结果的精度较其他方法有所提高,并具有简单、灵活的优点。     

弦单元及其在悬挂结构动力分析中的应用

为了合理地对工程中包含有悬挂系统的整体结构进行动力有限元分析， 根据有限单元法的基本原理推导了弦单元的刚度方程，求得了其等效横向 刚度矩阵；最后研究了自由悬挂系统的动力特性．     

圆形隧洞各向异性渗流的研究

基于饱和均质介质中,圆形隧洞各向同性稳定渗流的一个解析解,利用等效渗透系数法、坐标变换及流量不变法则,导出饱和成层介质中,圆形隧洞各向异性稳定渗流的解析解.结合有限单元法计算精度的要求,以及Polubarinova-Kochina隧洞渗流流量计算假设,对圆形隧洞各向异性稳定渗流孔隙水压分布的解析式进行替换和简化,推导出进行渗流场分析时确定有限元分析合理计算区域的理论公式,并对圆形隧洞各向异性稳定渗流孔隙水压和洞周流量的理论解和有限元解进行数值比对和误差分析.     

无压渗流渗出面确定的等效结点流量法

结合溢出点的特性及初流量法求解具有自由面渗流的优越性,以可能溢出点等效结点流量控制为外循环,非饱和区高斯点初流量为内循环,以渗出面只有流量流出,没有流量流入原则进行渗出面判定,编制了相应三维渗流分析程序并应用于某除险加固土石坝渗流分析中,结果表明：渗出面的确定方法可行．土石坝防渗加固效果明显,墙体不会产生渗透破坏．     

基于可靠度的黄河下游堤防工程渗流稳定分析

将可靠性理论应用于黄河下游堤防工程的渗流稳定分析,根据与Monte-Carlo方法相结合的有限单元法求解单元堤段渗流稳定的基本原理,基于黄河下游堤防渗透破坏的特点,建立了黄河堤防渗透破坏计算模型,介绍了基于Monte-Carlo方法的有限单元法求解典型堤防断面渗透破坏概率的计算步骤.最后利用基于Monte-Carlo方...     

结构任意截面内力求解的有限元内力法

在水工结构配筋设计中常需要事先确定截面上的内力值,有限元内力法是求解指定截面上内力值的有效方法,但实际工程中指定的内力求解截面往往不满足有限元内力法"截面由单元面组成"的要求,而且商业有限元软件无法导出内力求解公式中需要的结构整体劲度矩阵.针对这些问题,从有限元基本原理出发,提出 "采用等效变换将指定截面转换为由单元面组成" 和"采用单元劲度矩阵求解结点等效荷载"的方法,从而不需形成整体劲度矩阵,使得利用商业软件的计算结果来求解任意指定截面内力成为可能.最后通过一个简单算例对文中提出方法的正确性与可行性进行了验证.     

岩体结构面网格渗透张量分析方法

工程岩体的渗透特性研究是其渗流分析和稳定性设计的技术关键·作者以现场岩体渗透结构面概率模型统计资料为依据,采用离散介质方法建立典型裂隙网络模型,提出计算岩体结构面网络的等效渗透系数张量方法·应用等效连续介质模型分析方法研究区域岩体的渗流问题·针对现场的矿山边坡水文地质实体,运用该技术建立了反映岩体渗透各向异性的等效连续介质模型·实例分析表明,该技术途径不但适合工程岩体宏观渗流分析的特点,而且能够较为深刻地认识裂隙渗流的本质·     

渗流问题数值模拟计算的新方法

采用边界拟合坐标系，取渗流区域的边界为坐标线（二维问题）或坐标面（三维问题），用 Ｔｈｏｍｐｓｏｎ法将不规则渗流域转变为规则域，从而在规则域上求解转变后的渗流方程．避免了在不规则域上划分单元及迭代运算中生成新网格的繁琐过程，尤其对三维有自由面的无压渗流问题，本方法具有精度高而工作量小等优点。     

非饱和土弹塑性模型在土坝蓄水变形中的应用

对非饱和土弹塑性模型（BBM）进行了总结,推导出弹塑性矩阵,采用Galerkin加权余量法以及向后欧拉法对控制方程进行空间域和时间域上的离散,编制了用于分析非饱和土水一应力耦合问题的非线性有限元程序．采用该本构模型,对一均质堤坝在蓄水之后的流固耦合问题进行了有限元分析,结果表明堤坝内的渗流和变形的分布符合规律．     

等效连续裂隙岩体渗流与应力全耦合分析

推导了等效连续裂隙岩体渗透性与应力变形的关系,根据应变等效原则推导了等效连续裂隙岩体的弹性矩阵,首次采用四自由度全耦合法,建立了等效连续裂隙岩体渗流与应力的全耦合有限元方程组,进行了等效连续裂隙岩体渗流与应力的全耦合分析。     

求解渗流自由面的有限元混合不动点法

基于有自由面渗流分析的高斯点,建立了求解渗流问题的非光滑非线性方程组模型和求解此类问题的有限元混合不动点算法,此方法属于固定网格法,只需划分一次网格,不需要对数据做任何近似处理,完全利用程序迭代计算渗流自由面.讨论了非光滑方程组解的存在性和该不动点算法的收敛性,通过节点压强插值绘制出渗流自由面.算例结果表明,该方法简单且收敛速度快.对不动点法的收敛性分析为迭代法的收敛提供了理论依据.     

非饱和土壤渗流问题的混合元-变网格有限元方法

非饱和土壤渗流运动是多孔介质流体运动的一种重要形式，单纯运用有限元等方法不易得到最优误差估计，为更好地研究该问题，在混合元方法基础上提出了混合元-变网格有限元格式，并应用微分方程先验估计的理论和技巧，得到最优阶的L^2-模误差估计结果，从理论上证明了该方法在研究多孔介质流体运动中的实用性，从而为更好地研究非饱和土壤渗流问题提供了直接的理论依据．     

改进复合单元渗透系数调整法渗流自由面分析

在复合单元渗透系数调整法基础上,应用连续的罚函数对其改进,优化了此类算法的稳定性,提高了收敛速度和精度,并将此改进算法编制成相应的有限元分析程序,通过算例验证,可应用于稳定和非稳定渗流自由面边界复杂渗流场的高效精细模拟.     

曲面均匀网格轮式行进生成法

提出一种新的曲面有限元网格生成方法，将曲面边界离散相邻节点连接构造线段序列，并从序列中的线段出发，以线段的空间位置为依据，构造圆，再利用该圆在曲面上布点筛选求取网格节点，生成单元，使节点一单元一线段交替生成，该方法实现了曲面均匀网格的自动生成，网格质量良好，适用于不同形式的曲面，具有很大程度的通用性，有利于解决有限元分析建模的瓶颈之一——曲面网格自动生成问题。     

基于非稳定渗流有限元的坝坡稳定分析

对库水位连续下降的非稳定渗流问题 ,建立与定解问题相对应的泛函 ,再根据有限元原理 ,以变动的渗流自由面为流量 ,补给边界条件来模拟非稳定渗流过程 ,确定各个时段的自由面及等势线 .在渗流分析基础上进行坝坡稳定分析 ,采用BISHOP法 ,并搜寻最危险滑弧 ,计算相应的安全系数 .通过工程实例分析 ,得出了一些有益的结论     

无网格伽辽金法在施工导流二维数值模拟中的应用

针对有限元法等传统数值计算方法存在受单元网格限制、前后处理工作复杂的问题,提出应用一种数值计算方法--无网格伽辽金法,对具有简单边界条件的水利水电工程施工导流的恒定二维浅水流动问题进行了分析、计算.同时利用有限元法进行了对比计算,从流速、水位的计算结果来看,两种计算方法结果相近、误差较小,表明采用无网格伽辽金法解决此类问题是可行的.     

有限单元法的非结构化网格及其自动生成

本文介绍了有限单元法非结构化网格的基本原理及其自动形成方法。由于有限元法是一种离散的数值求解方法，其近似求解方法的精确度，很大程度上取决于所形成网格的质量。另外，对于工程中一些形状复杂的问题，一般的网格生成方法很难对其进行离散。非结构化网格及其自动生成，使复杂形状的工程问题容易地进行离散，改善所形成网格的质量，提高近似计算的精确度，并且在有限元网格修正自适应分析中具有重要作用。     

岩体三维主干裂隙网络渗流模型

渗流分析是水电工程中一项重要的研究内容,大坝周转岩体的渗流通道大多为裂隙网络,用连续介质渗流模型难以解决这一问题,针对岩体主干裂隙网络渗流特征,建立了岩体三维主干裂隙风络渗流模型,运用有限元数值方法,结合算例分析了坝基岩体主干裂隙网络渗流问题,理论解与有限元解元解的比较结果表明,岩体三维主干裂隙网络渗流模型及其有限元算法是有效而实用的。     

基于有限元自动生成系统(FEPG)的组合网格算法

针对工程计算中经常出现的局部特性(特别是奇性)问题,和以往解决此类问题的算法的局限性,提出一种基于有限元自动生成系统(FEPG)的组合网格算法.该算法采用两套网格求解,在整个求解区域采用较粗网格,不考虑奇异的影响;而在奇异附近区域采用较细的网格,考虑奇异的影响;整体粗网格求解和局部细网格求解反复迭代,求得最终结果.该算法用于实际工程计算的迭代次数少,与常用的有限元方法所求得的解相符合,为求解大型实际复杂问题提供了一个好的算法和思路.     

A-?方法在高频电磁场计算中的应用

基于节点的有限元方法具有网格剖分、构造高阶基函数容易的优点。由于在节点上定义场量,节点有限元更适用于多物理问题。但节点有限元方法直接求解电磁场会出现伪解,场量在不均匀介质中不连续等问题。基于节点的A-φ方法可以有效避免传统节点有限元方法存在的问题。本文研究A-φ方法的工程应用,研究开域和闭域问题中如何设置关于矢势和标势函数的边界条件,特别是波导问题和理想导体球的散射问题,讨论了端口边界条件,辐射边界条件的使用方法。对于理想导体边界条件采用了阻抗边界条件,与端口条件配合,克服方程的奇异性。数值卖验比较分析了A-φ法节点有限元和棱边法有限元的计算结果,验证了A-φ法节点有限元的正确性和有效性。