改进的EP法在稳定渗流分析中的应用

为了提高稳定渗流分析中自由面处的模拟精度，通过分析EP(extended pressure)法的控制方程及其离散格式，指出分段连续函数hε(p)中参数ε的选取对计算精度影响显著，提出ε取值为网格点在y方向上间距的一半更合理，得到了改进的EP法，通过数值实验，将此方法的结果与已有结果对比，证实改进后自由面光滑，逸出点位置合理，模拟精度较高。     

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

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

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

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

基于六节点三角形单元优化虚单元法

选取六节点三角形单元优化了虚单元法,弥补了三节点单元求解的自由面是折线的不足,同时解决了等参四边形单元难以寻找节点移动路径的问题.此外,在计算程序中引入节点相关次数矩阵,实现了对虚单元法的进一步优化,提高了计算效率和计算精度.与等效渗透系数法计算结果分别同电模拟试验解比较表明,本文方法在大网格单元划分下仍能得到准确且连续光滑的自由面曲线,计算精度没有因单元数量的减少而降低;在相同的精度条件下,计算效率较三节点单元提高.     

基于稳态随机动力响应的蒙特卡洛模拟精度研究

以高斯白噪声作用下的单质点线性系统为例,探讨了系统稳态随机响应的蒙特卡洛模拟精度.结果表明,模拟的精度与样本数、样本值域分段数和系统参数的取值有关:样本数的选择应综合考虑计算精度与计算量,而合理选择分段数可以有效改善计算结果的精度.     

蒙德法定量安全评价的误差分析与参数选择

为分析与讨论蒙德法(Mond)优缺点、评价精度及适应性,对参数取值的主观偏差影响、参数选取的合理性、对蒙德法中燃烧热的温度、单元中重要物质选取、混合物质影响以及Mond法算法误差传递进行了研究,采用误差传递法对偏差的传递进行了理论分析与论证.结果表明蒙德法中燃烧热、重要物质选取、混合物质等参数选取是合理的,同时阐明了偏差传递过程及对评价结果的影响;最后针对具体问题提出了建议和解决办法,如改进计算模型和量化参数取值等.     

运动潜体尾迹及其与随机海面相互作用的数值模拟

研究了水下运动物体产生的自由面尾迹与随机海面相互作用的问题.以三维非定常不可压缩雷诺平均(RANS)方程为控制方程,采用重整化群(RNG)k-ε湍流模型,应用体积分数(VOF)方法处理自由面,对水下运动球体产生的自由面尾迹特征进行了数值模拟与分析.结合海浪谱模型,数值分析了水下运动球体产生的自由面尾迹与随机海面的相互作用机理及其表现特征.     

相场法模拟过冷熔体枝晶生长的界面厚度参数的取值

研究了用相场法模拟纯物质过冷熔体枝晶生长过程中,界面厚度的取值对模拟结果可靠性的影响及影响界面厚度取值的因素;然后,在选取合理界面厚度取值的基础上,计算了不同过冷度下纯镍过冷熔体枝晶生长的速度并与实验进行了比较.研究表明:为了获得可靠的计算结果,界面厚度的取值应足够小;合理的界面厚度的取值由过冷度、各向异性、界面动力学、热扩散系数等参数综合决定;在选取合理界面厚度的基础上,模拟结果与实验吻合,证明界面厚度足够小时可获得可靠的模拟结果.为相场法模拟过冷熔体枝晶生长时的界面厚度参数的取值提供了依据.     

均质承压含水层稳定渗流的基本解法

基本解法是一种新型的无网格法,仅需边界配点信息即可对待求问题进行数值模拟,具有理论简单、编程容易和精度高等优点.利用地下水稳定渗流问题的数学模型,以基本解法对均质承压含水层的稳定渗流进行了探讨.与解析解和其它数值方法计算结果进行了对比研究,数值结果表明基本解法数值模拟稳定渗流问题具有令人满意的结果.     

FEPG软件在渗流分析中的应用

应用FEPG软件自动生成渗流计算程序,并将该程序应用于三维渗流场分析中.在稳定渗流计算时,采用死活节点法迭代求解自由面,并对该方法进行了改进;在饱和-非饱和渗流计算时,对计算参数加以适当的处理,采用合适的时间差分格式求解方程.算例表明,将FEPG应用于渗流分析是可行的.     

紊动浮射流中的普朗特数取值研究

采用k-ε双方程模型研究存在浮力的传热流动问题时，数学模型中的参数-紊动普朗特数的取值大小对模拟计算的流速和温度分布存在不同程度的影响．本文在已有实验推荐的取值范围的基础上，采用圆柱轴对称坐标系下的k-ε模型，研究了不同紊动普朗特数取值对数值模拟结果逼近真实值的影响程度．研究得到紊动普朗特数取值的变化对轴线流速的大小和分布影响很小，即数值模拟得到的轴线上的流速仍能很好的预报出实验或真实流速值；但对温度分布影响较大，且紊动普朗特数的取值越小，数值模拟出来的轴线上的温度值与真实值相比偏小．得出了紊动普朗特数在大小不同雷诺数下的合理取值变化范围，其结论对采用数学模型研究各种复杂流动具有指导意义．     

文献老化B-K方程的改进式

针对文献老化B-K方程在实际应用中存在系数a,b的取值会超出[0,1],范围的问题,将B-K方程中的e~x或e~(2x)修改为enx,用改进后的方程进行数据模拟.结果表明,对B-K方程的改进,不仅有合理的意义,也提高了模拟的精度,特别是可以解决B-K方程中有权重意义的系数a,b的取值范围问题.     

山体地形下低矮房屋数值风洞模拟的计算域设定

针对山体地形下低矮房屋数值风洞模拟,研究其计算域大小的设定问题,并提出设置的参考方法.通过建立正弦山丘轮廓函数,对4种不同形状因子的山丘进行计算流体动力学数值模拟.采用CFX软件作为计算平台,对影响计算结果精度的各个参数进行分析,在此基础上,提出合理的取值范围.     

基于VOF法的CSG坝非稳定渗流场数值分析

当前国内外对胶凝砂砾石大坝(CSG Dam)渗流场的研究均局限于稳定渗流场,但由于库水位的不断波动,大坝在运行期大都处于非稳定渗流状态,保证这种新型坝在库水位变化时的安全稳定非常重要.采用有限体积法及计算流体力学中的VOF(volume of fluid,VOF)法,对CSG坝两个极端工况下的非稳定渗流场进行模拟,获得...     

海堤非稳定渗流场渗透参数反演分析

针对受到快速涨落潮位作用的海堤进行非稳定渗流有限元分析,在此基础上结合实测渗压序列构建目标函数对海堤的渗透参数进行反演研究.反演中采用分层搜索技术,先由逐步扫描法确定出渗透参数一个合理的小区间,再用一阶梯度优化方法在该小区间上进行迭代优化以获得合适的渗透参数.以反演得到的渗透参数进行海堤非稳定渗流数值模拟,算例表明,计算及反演步骤合理可行,模拟计算得到的渗压结果与实测序列在数值和动态变化规律上都取得了理想的拟合效果.     

混凝土塑性损伤模型参数的研究

基于相关的试验研究成果,提出了一种确定混凝土塑性损伤模型(KC模型)强度参数取值的方法.针对KC模型损伤参数取值依赖于单元尺寸的问题,以我国混凝土规范为基准,阐述了损伤参数值的调整方法,使得数值计算结果更加合理.利用ANSYS/LSDYNA模拟混凝土常规三轴压缩试验,KC模型分别采用本文确定和程序自动生成的强度参数值,并将计算结果与国内外试验结果比较,研究表明:采用本文提出的方法确定KC模型强度参数值比较合理.     

不同湍流模型在列车外流场计算中的比较

通过对简化缩比的列车模型在无侧偏角工况下分别进行模型风洞试验和数值模拟,评估了5种常用湍流模型在列车外流场模拟中的计算精度.数值模拟与风洞试验结果的分析对比表明:Realizable(可实现)k-ε在阻力的计算上精度最高,误差最小为2.5%;SST(剪切应力输运)k-w在升力的计算上更为精确,误差最小为0.5%.Standard(标准)k-w和SST k-w对准二维分离的捕捉与油流显示试验最接近,Realizable k-ε则对三维分离位置的预测更准确.     

流固耦合法模拟岩石爆破时耦合范围的确定

为研究岩石爆破问题模拟中的流固耦合算法以及提高计算效率,文章采用LS-DYNA软件开展大量计算,分别讨论了不同装药半径(r)、不耦合系数(K)和网格精度条件下流固耦合区域半径(R)对单元峰值应力衰减曲线的影响,并结合应力云图和应力时程曲线分析,确定出合理的流固耦合范围。结果表明:随着耦合装药半径r的增大,最小合理流固耦合半径R_(min)与r比值呈减小趋势;相同装药半径下随着K增大,R_(min)/r逐渐增大;网格精度对流固耦合半径取值影响较小;相同条件下,当R较小时,应力云图不稳定且应力时程曲线后期会迅速衰减,适当增大流固耦合半径会逐渐消除该现象。对于一般的岩石爆破问题,建议R/r取5～10。     

嵌岩桩水平承载力地基反力Matlock法分析

嵌岩桩在海上风电、跨海大桥等结构中的应用十分广泛,但在嵌岩桩水平承载特性方面的研究仍有较多不足。弹性地基反力法中的m法在我国是用于计算水平承载桩承载力时应用最多的方法,但地基反力m法针对嵌岩桩缺少合理的取值方法。本文基于现场试验数据建立并验证有限差分数值模型,通过数值模拟与参数分析,研究了岩体力学参数和桩基参数等因素对桩岩相互作用的影响,提出了m法中地基反力系数比例系数的拟合计算公式。通过现场试验数据对文中提出的方法进行了验证,计算结果与实际情况较为相符,通过与其它方法比较,本文提出的修正m法相对非线性的p-y曲线法具有取值计算较为简便的优势。     

近自由面三维振动水翼的水动力分析

采用非定常面元法对近自由面三维振动水翼进行势流数值分析.在数值计算中,将Rankine源和偶极子置于边界面上,用时间步进法模拟水翼的势流场和自由表面波形.在自由面采用线性自由面边界条件,在尾涡面上布置偶极子及采用尾涡存储效应以满足库塔条件.确定了数值计算模型和相关参数,对不同展长的水翼在不同浸深、不同航速和不同振动频率...