基于B样条物质点法的溃坝流模拟研究.pdf
《基于B样条物质点法的溃坝流模拟研究.pdf》由会员分享,可在线阅读,更多相关《基于B样条物质点法的溃坝流模拟研究.pdf(10页珍藏版)》请在咨信网上搜索。
1、 应用数学和力学编委会,:基于 样条物质点法的溃坝流模拟研究徐云卿,周晓敏,赵世一,徐聖飞,孙 政,(江西理工大学 土木与测绘工程学院,江西 赣州;江西理工大学 江西省环境岩土与工程灾害控制重点实验室,江西 赣州)摘要:溃坝流是水利工程中常见的一种自由表面流动问题,准确模拟溃坝流问题具有重要的工程意义 样条物质点法()作为一种物质点法(,)的改进算法,其提高了物质点的计算精度和收敛性,且在自由表面流动问题中有独特的算法优势 基于 样条物质点法,通过引入人工状态方程,发展了一种弱可压缩 样条物质点法();开展溃坝流问题的模拟研究,分析 样条插值基函数阶数对模拟结果的影响 结果表明:模拟所得流体波
2、前位置、波前流速及给定位置处的高程变化与已有实验结果吻合较好;同时,随着 样条基函数阶数的增加,模拟结果与试验结果吻合度逐渐提高 随着基函数阶数的增加,计算耗时呈约 倍增长;不同阶次 样条物质点法的计算耗时随背景网格尺寸的增长率基本一致,约呈线性增长 验证了弱可压缩 样条物质点法模拟溃坝流问题的有效性,为模拟溃坝流问题提供了一种新的思路和方法 关 键 词:溃坝流;样条物质点法;自由表面流;人工状态方程中图分类号:文献标志码:,(,;,):,(),(),(),应用数学和力学 卷 期 年 月 ,收稿日期:;修订日期:基金项目:国家自然科学基金项目(;);江西省主要学科学术和技术带头人青年项目();
3、国家级大学生创新创业训练项目();江西省大学生创新创业训练项目()作者简介:徐云卿(),男,硕士(:);孙政(),男,副教授,博士,硕士生导师(通讯作者:)引用格式:徐云卿,周晓敏,赵世一,徐聖飞,孙政 基于 样条物质点法的溃坝流模拟研究 应用数学和力学,():,:;引 言水库大坝是土木工程中常见的工程结构,在防洪、能源生产、蓄水等方面都发挥着重要作用 但大坝在具有社会经济效益的同时,也存在着溃坝的潜在风险,一旦发生溃坝,将会对人民生命财产安全造成极大的威胁 因此,开展水坝的溃决过程及其演进规律研究,对于溃坝洪水灾害的预测和防治具有重要意义 溃坝流是一种典型的自由表面流动问题,研究内容主要包括
4、波前到达时间、波前流速、波前位置和沿程各处的水位等 目前,研究方法主要有理论分析、物理试验、数值模拟 随着数学理论和计算机技术的高速发展,数值模拟广泛应用于溃坝水流的研究 数值模拟研究相比较于物理试验具有更高的灵活性,不局限于场地设置,计算周期相对较短,能够大幅度降低试验成本等优点 等采用 法的一维模型耦合采用有限差分法的溃坝二维模型求解溃坝扩散波动情况;胡四一和谭维炎以 格式预测了高坝瞬溃的溃坝波演进过程;王大国等采用 有限元法建立了水波模型,精确模拟和分析了下游河床无水时溃坝模型的自由水面运动特征;张建伟等采用光滑粒子流体动力学方法()建立了包含下游水位的溃坝模型,对于溃坝流冲击过程中波前
5、到达时间、波前流速、波前位置能够做到精确模拟 但溃坝流涉及到水波演进大变形问题,有网格类的有限差分法()、有限体积法()、有限元法()等方法在模拟时会出现网格畸变并引起误差 尽管无网格化的 避免了网格畸变,但由于追踪的是网格边界上质量、动量和能量的通量流动,需求解非线性对流项,增加了求解难度,且不易于追踪各质点的运动时间历程,使得计算效率低 因此,需要开发一种可以减少误差和提高计算效率的数值模拟方法 物质点法()作为一种相对新型的粒子型算法,其具有 粒子和 背景网格双重描述,既包含了 粒子型无网格算法的优势,又可避免网格畸变和求解非线性对流项,可方便地追踪自由表面和粒子信息的时间历程;同时利用
6、 型背景网格求解动量方程和梯度,具有计算效率和处理基本边界条件的优势 物质点法已广泛应用各个工程领域问题中 然而,传统的物质点法采用线性插值形函数来实现物质点和背景网格节点之间的信息映射 由于网格边界处的形函数梯度不连续,当物质点穿越背景网格边界时,将产生网格穿越误差,引起数值振荡 等采用高阶 样条基函数替换线性插值基函数,节点自由度空间代替背景网格节点,减少了网格穿越误差,提高了物质点法的计算精度和收敛性 目前 样条物质点法()已经用于模拟自由表面流体问题,如具有自由表面的复杂流动问题,流体结构相互作用问题等 本文主要研究两个方面:一方面,运用弱可压缩 样条物质点法(),对溃坝流问题进行模拟
7、研究,并将模拟的流体波前位置、波前速度及给定位置的高程变化与 等的实验结果进行对比分析,验证运用 样条物质点法模拟溃坝流问题的有效性;另一方面,基于不同阶次 样条插值形函数对溃坝流问题进行模拟,并与已有实验结果进行对比,分析阶数对模拟精度和计算效率的影响 计 算 方 法 弱可压缩 样条物质点法 样条基函数 样条基函数的定义有很多种,尽管所得的表达式不同,但其本质是一致的 本文选用计算稳定、易于理解的 递归公式来计算 样条基函数 在 递归公式中,零阶基函数,()如下应 用 数 学 和 力 学 年 第 卷所示:,()对于 ,基函数由 递归公式定义为,(),(),(),()式中,是第 个节点,是基函
8、数的总数目,是基函数的阶数 由式()可以看出,如果基函数的阶数为零(),则基函数呈阶梯分布,即在第 节点区间,)上基函数,()等于,而其他节点区间等于 采用,()和,(),通过式()可以得到,();当所有,()计算完成,采取同样的方法可以得到,(),持续这个计算过程直到所有需要的,()计算完毕,其递归过程如图 所示 图 递归示意图 ()参数网格空间()节点自由度空间()()图 网格空间示意图 假设,()的导数也可以用递归公式定义,即,(),(),()()本文研究的是二维溃坝流模拟实验,基函数是双变量基函数,假定节点向量 ,第 期 徐云卿,等:基于 样条物质点法的溃坝流模拟研究 和 ,和 是参数
9、的方向,和 分别为沿图()中,方向的基函数数目 利用张量积结构双变量基函数定义为(,),(,)(,),(),(),()式中,;,三变量 样条基函数也可以类似于双变量基函数推导得到 图()参数网格中的空心圆和实心圆分别表示节点向量和物质点 图()给出的是节点自由度空间,图中(,)表示节点自由度空间中对应的(,)自由度;物质点(,)在节点自由度(,)上对应的 样条基函数记作(,),(,)(,),其中(,)分别表示在节点矢量 和方向上 样条基函数的阶数 样条物质点法 样条物质点法是基于物质点法改进而来的,样条插值基函数阶数为一阶时与物质点法中的传统线性插值形函数是一样的,即 样条插值基函数阶数为一阶
10、时的 样条物质点法就是物质点法,两者求解步骤基本一致 两者的不同之处在于:在物质点法中,物质点信息采用线性插值形函数映射到相应背景网格节点上,在背景网格节点上施加本质边界条件 而在 样条物质点法中,物质点信息采用 样条基函数映射到节点自由度空间的各节点自由度上,在节点自由度空间施加本质边界条件 其具体算法实现过程如下:)划分参数网格,将研究对象离散成一组物质点)初始化物质点的位置、质量、动量等物质信息)采用 样条基函数,在第 个计算时间步,将物质点上的质量和动量信息映射到节点自由度空间的节点自由度上:(,)(,),(,)(),()(,)(,),(,)()(,),()式中上标 表示计算时间步;下
11、标 示物质点的变量;表示物质点 的位置矢量;为质量;为速度矢量 )计算节点自由度上的内力(,)和外力(,):(,)(,),(,)(),()(,)(,),(,)()(,),(,)(),()式中,(,),(,)表示 样条插值基函数的梯度;为 应力张量;和 分别为体力和表面力矢量;为体积;表示面力边界厚度)计算节点自由度上的加速度和速度:(,)(,)(,)(,),()(,)(,)(,),()式中,表示加速度矢量)在节点自由度空间,施加本质边界条件)将各节点自由度的加速度和速度信息,通过 样条基函数映射回各物质点,得到第 个时间步上物质点的速度和位置:(,),(,)()(,),()(,),(,)()(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 物质 溃坝流 模拟 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。