一种基于物质点法的多尺度泥石流现象模拟方法技术

本发明专利技术公开了一种基于物质点法的多尺度泥石流现象模拟方法，涉及计算机物理动画技术领域，该方法包括以下步骤：S1：基于自适应物质点方法模拟泥石流液相浆体的多尺度特征；S2：基于自适应物质点有限元方法模拟泥石流固相物体的多尺度特征；S3：基于物质点法模拟多尺度液固耦合；S4：对泥石流动画进行多尺度真实感渲染。本发明专利技术解决了在泥石流动画模拟过程中，如何构造描述泥石流灾害场景多种尺度特征的有效物理模型，并建立不同模型间一致性约束以实现多尺度数值求解和渲染的问题。本发明专利技术能够真实感模拟泥石流现象多尺度特征，生成具有照片级逼真度的泥石流动画，为防灾救灾、数字娱乐、影视特效等应用提供泥石流现象的影像数据。

A method of multi-scale debris flow simulation based on material point method

The invention discloses a multi-scale debris flow phenomenon simulation method based on material point method, which relates to the technical field of computer physical animation. The method comprises the following steps: S1: simulating the multi-scale characteristics of debris flow liquid slurry based on adaptive material point method; S2: simulating debris flow based on adaptive material point finite element method Multi-scale characteristics of solid objects; S3: simulation of multi-scale liquid-solid coupling based on material point method; S4: multi-scale realistic rendering of debris flow painting. The invention solves the problem of how to construct an effective physical model describing multi-scale characteristics of debris flow disaster scene in the process of debris flow simulation, and how to establish consistency constraints between different models to realize multi-scale numerical solution and rendering. The invention can simulate the multi-scale characteristics of debris flow phenomena realistically, generate debris flow pictures with photo-level fidelity, and provide image data of debris flow phenomena for disaster prevention and relief, digital entertainment, film and television effects and other applications.

【技术实现步骤摘要】
一种基于物质点法的多尺度泥石流现象模拟方法
本专利技术涉及计算机物理动画
，具体涉及一种基于物质点法的多尺度泥石流现象模拟方法。
技术介绍
泥石流作为自然环境中普遍存在的流体类地质灾害，给人类生命财产与经济建设带来了巨大损失。长期以来，人们研究泥石流的运动规律以期对其进行预防和治理。但泥石流灾害的不可预知性和强大破坏性，使人们不仅难以通过预测泥石流发生的时间、地点来安排观测和拍摄，而且难以承受其瞬时产生的破坏力带来的危险和不确定性，这造成了泥石流灾害真实影像数据的缺失。然而，现实生活中防灾救灾、影视特效以及数字娱乐等领域都对泥石流灾害的真实感图像和动画具有迫切的需求。因此，对泥石流灾害场景的真实感建模和仿真具有重大现实意义。近年来，为了获得真实物体和现象的照片级真实感模拟效果，人们综合连续介质力学、数值计算以及计算机图形学的真实感绘制算法，提出了基于物理的模拟技术，真实地模拟了多种复杂自然现象，如：波涛汹涌的流水，袅袅上升的烟雾，熊熊燃烧的火焰，固体的变形和断裂以及固流交互等。然而，对于复杂泥石流灾害场景的真实感模拟，计算机图形学领域中却极少有人涉及，少数学者利用传统的基于物理的模拟技术对泥石流模拟进行了初步探索，但主要关注液相浆体的流动及其与环境的简单交互，不能够完整而逼真地重现泥石流灾害场景以满足实际的应用需求。造成泥石流灾害场景真实感模拟困难的最主要原因在于难于实现其多尺度特征的统一物理建模。泥石流灾害场景的多尺度特征具体体现在：既有大尺度的液相浆体、固相物体的宏观运动和双向交互，又有小尺度的非均质组成特性、漩涡细节以及固体断裂和裂纹扩展。传统的物理模拟方法中，欧拉方法适合模拟大尺度的特大变形和不可压缩流体场景，却难以跟踪流体界面和小尺度漩涡细节特征；拉格朗日方法避免了非线性对流项且方便跟踪材料的界面，尤其是其中的无网格粒子方法非常适于模拟小尺度细节特征，但在模拟大尺度流体场景时计算过度复杂，并且本质边界条件施加繁琐。由此可见，欧拉方法和拉格朗日法具有各自的优势，但也都存在各自的不足，难以利用单一的方法对泥石流场景的多尺度特征进行模拟。目前，有许多研究工作对不同物理方法进行耦合形成的混合方法能充分吸收各自的优势，克服各自的缺点，可解决只用纯粹的单一方法所解决不了的问题，但存在的关键问题是不同模型交界面的无缝处理极其复杂，难以实现具有多尺度物理特征的泥石流现象的真实感模拟。
技术实现思路
为了克服现有技术中存在的问题，本专利技术提供一种基于物质点法的多尺度泥石流现象模拟方法，能够对泥石流场景的多种尺度特征构造具有一致性约束的有效物理模型，基于物质点法对物理模型进行多尺度数值求解，同时实现泥石流场景的多尺度渲染。为实现上述目的，本专利技术所述的基于物质点法的多尺度泥石流现象模拟方法包括包括以下步骤：S1：基于自适应物质点方法模拟泥石流液相浆体的多尺度特征；S2：基于自适应物质点有限元方法模拟泥石流固相物体的多尺度特征；S3：基于物质点法模拟多尺度液固耦合；S4：对泥石流动画进行多尺度真实感渲染。所述步骤S1进一步包括：确定泥石流浆体成分比重；构造描述泥石流液相浆体非均质特征的近似物理模型；建立描述场景复杂度的自适应采样判断准则来控制物质点的分裂和合并。对于步骤S1所述的泥石流液相浆体非均质特征的近似物理模型的构造，在任一物质点p上通过分别求解弹性固体和不可压缩牛顿流体的本构方程并进行融合来计算应力σp，本构方程为：其中，固体应力部分σs为位移向量u的函数采用非线性弹性材料的本构方程进行计算；流体应力部分σf为速度向量v的函数采用不可压缩流体本构方程求解；对于关键的随机函数模型η(x)∈[0,1]，根据实验流体力学中的测量数据建立统计模型来建立，从而准确表示泥石流浆体中固液两相的比重。所述步骤S2进一步包括：初始时采用高效、精确的非线性有限元法离散固相物体并模拟其宏观运动，在模拟过程中综合等效塑性应变和单元特征长度比值来计算有限单元的畸变程度，然后把畸变程度超过用户设定阈值的有限单元标记为潜在断裂区域并转化为物质点，最后在物质点区域稳定模拟断裂现象。所述步骤S3中对液相浆体和固相物体的多尺度耦合模拟进一步包括：利用背景网格多重速度场设计物质点法的双向交互计算模型以模拟宏观尺度运动，基于边界层理论建立物质点法的湍流模型以模拟小尺度漩涡细节。优选地，所述步骤S3中对液相浆体和固相物体的多尺度耦合模拟进一步包括：S3.1：基于背景网格多重速度场进行物质点法的双向交互计算，模拟宏观尺度上的运动和形变；对于液相浆体和固相物体物质点区域的交互，利用一种基于局部多重背景网格的交互算法来计算物质点法表示的两个不同物体之间的双向交互；对于液相浆体和固相物体有限元区域的交互，利用基于同一全局背景网格的交互算法来计算分别由物质点法和有限元法表示的两个不同物体之间的双向交互；S3.2：在物质点法流体中对边界层理论漩涡模型进行近似模拟小尺度漩涡细节；假定边界层由到固体表面距离小于设定阈值的液相物质点组成；计算固体表面结点i与相邻液相物质点的加权平均相对速度并根据其与法线ni的几何关系检测边界层脱落点SP；计算边界层物质点和脱落点的位置关系以确定实际脱落的边界层物质点BLPj，并设计随机函数赋予其涡度ωj和湍流动能kj；把远离固体的BLPj转化为涡度粒子VPj而形成涡度场；在物质点法背景网格上结合Biot-Savart定律和k-ε湍流模型计算涡度场引起的扰动速度场uv，并将其叠加到积分动量方程得到的平均速度场U上的方式来合成漩涡细节。所述步骤S4包括：首先设计一种高效的基于八叉树的自适应表面重建方法以捕获泥石流表面细节特征的几何结构，然后设计一种点和面的光线跟踪混合绘制方法实现多尺度绘制方法。所述的高效的基于八叉树的自适应表面重建方法包括检测出位于泥石流表面的物质点，只在泥石流表面附近建立隐式距离场，内存消耗和计算复杂度只取决于泥石流表面积。步骤S4所述的对泥石流动画进行多尺度真实感渲染需要建立一种高效的基于八叉树的自适应流体表面重建方法，包括：首先，对于任一物质点i，本专利技术计算其再归一化矩阵Bi：其中，Vj为物质点体积，W为再归一化高斯核函数，h＝4d0，d0为物质点初始间距，▽为梯度算子,为张量积算子。然后，计算矩阵Bi的最小特征值并根据Bi和实验统计阈值λthreshold的大小关系判定粒子i是否属于泥石流表面；最后，根据泥石流表面的物质点创建基于八叉树的自适应距离场，通过在流体表面附近建立隐式距离场使内存消耗和计算复杂度只取决于流体表面积。所述的点和面的光线跟踪混合绘制方法包括对隐式流体表面和表示飞溅液滴、泡沫的溅射粒子采用光线跟踪算法在屏幕空间进行直接绘制，根据光在气-液交界面传播的菲涅耳效应计算光线的反射和折射。步骤S4所述的对泥石流动画进行多尺度真实感渲染还需要实现表示泥石流的点和面的光线跟踪混合绘制，包括：建立隐式流体表面和溅射物质点的高效的八叉树组织结构；对于光线和溅射粒子的求交计算，把光线假想成具有一定半径的光柱体，求交计算转化为寻找光柱体内是否包含溅射粒子，如果检测到交点，则搜索交点附近的溅射粒子，计算真实交点及其属性；对于光线和隐式流体表面的求交，采用三线性插值计算交点及其属性；根据光线与场景的求交结本文档来自技高网...

【技术保护点】
1.一种基于物质点法的多尺度泥石流现象模拟方法，其特征在于，包括以下步骤：S1：基于自适应物质点方法模拟泥石流液相浆体的多尺度特征；S2：基于自适应物质点有限元方法模拟泥石流固相物体的多尺度特征；S3：基于物质点法模拟多尺度液固耦合；S4：对泥石流动画进行多尺度真实感渲染。

【技术特征摘要】
1.一种基于物质点法的多尺度泥石流现象模拟方法，其特征在于，包括以下步骤：S1：基于自适应物质点方法模拟泥石流液相浆体的多尺度特征；S2：基于自适应物质点有限元方法模拟泥石流固相物体的多尺度特征；S3：基于物质点法模拟多尺度液固耦合；S4：对泥石流动画进行多尺度真实感渲染。2.如权利要求1所述的基于物质点法的多尺度泥石流现象模拟方法，其特征在于，所述步骤S1进一步包括：确定泥石流浆体成分比重；构造描述泥石流液相浆体非均质特征的近似物理模型；建立描述场景复杂度的自适应采样判断准则来控制物质点的分裂和合并。3.如权利要求2所述的基于物质点法的多尺度泥石流现象模拟方法，其特征在于，对于步骤S1所述的泥石流液相浆体非均质特征的近似物理模型的构造，在任一物质点p上通过分别求解弹性固体和不可压缩牛顿流体的本构方程并进行融合来计算应力σp，本构方程为：其中，固体应力部分σs为位移向量u的函数采用非线性弹性材料的本构方程进行计算；流体应力部分σf为速度向量v的函数采用不可压缩流体本构方程求解；对于关键的随机函数模型η(x)∈[0,1]，根据实验流体力学中的测量数据建立统计模型来建立，从而准确表示泥石流浆体中固液两相的比重。4.如权利要求1所述的基于物质点法的多尺度泥石流现象模拟方法，其特征在于，所述步骤S2进一步包括：初始时采用高效、精确的非线性有限元法离散固相物体并模拟其宏观运动，在模拟过程中综合等效塑性应变和单元特征长度比值来计算有限单元的畸变程度，然后把畸变程度超过用户设定阈值的有限单元标记为潜在断裂区域并转化为物质点，最后在物质点区域稳定模拟断裂现象。5.如权利要求1所述的基于物质点法的多尺度泥石流现象模拟方法，其特征在于，所述步骤S3中对液相浆体和固相物体的多尺度耦合模拟进一步包括：利用背景网格多重速度场设计物质点法的双向交互计算模型以模拟宏观尺度运动，基于边界层理论建立物质点法的湍流模型以模拟小尺度漩涡细节。6.如权利要求1所述的基于物质点法的多尺度泥石流现象模拟方法，其特征在于，所述步骤S3中对液相浆体和固相物体的多尺度耦合模拟进一步包括：S3.1：基于背景网格多重速度场进行物质点法的双向交互计算，模拟宏观尺度上的运动和形变；对于液相浆体和固相物体物质点区域的交互，利用一种基于局部多重背景网格的交互算法来计算物质点法表示的两个不同物体之间的双向交互；对于液相浆体和固相物...

【专利技术属性】
技术研发人员：邵绪强，徐伟峰，王洪涛，
申请(专利权)人：华北电力大学保定，
类型：发明
国别省市：河北,13