理解和表征模型的不确定性
“不确定性量化模块”可用于理解模型不确定性的影响——研究关注量与模型输入变化的相关性。其中提供的通用接口可用于筛选、灵敏度分析、不确定性传播和可靠性分析。
“不确定性量化模块”可以有效地测试模型假设的有效性,可靠地简化模型,帮助您了解关注量的关键输入,探索关注量的概率分布以及确认设计的可靠性。模型的准确性得到保障以及对关注量的深入理解可以帮助您降低生产、开发和制造成本。
“不确定性量化模块”可与整个 COMSOL 产品库中的产品相结合,用于分析电磁、结构、声学、流体流动、传热和化学工程仿真的不确定性。不仅如此,您还可以将其与“CAD 导入模块”、“设计模块”或任意 CAD LiveLink™ 产品结合使用。
输入参数和关注量
在运行不确定性量化研究时,您可以根据 COMSOL Multiphysics® 模型的解来定义一组关注量。这样一来,关注量就成了输入参数的函数。
在结构分析的情况下,关注量可以是最大位移、应力或偏转角。对于传热或 CFD 分析,关注量可以是最高温度、总热损耗或总流体流率。对于电磁仿真,则可以是电阻、电容或电感。由于“不确定性量化模块”适用于使用 COMSOL Multiphysics® 软件计算的任何物理场模型,以及各种求解场量的数学表达式,因此关注量的选择是无穷无尽的。
任何不确定的模型输入,无论是物理场设置、几何尺寸、材料属性还是离散化设置,都可以作为输入参数,并且任何模型输出都可以用来定义关注量。
筛选
筛选,MOAT 研究类型实现了一种轻量级的全局筛选方法,可以对每个输入参数的重要性进行定性度量。该方法纯粹基于样本,使用 Morris one-at-a-time(MOAT)方法,所需的 COMSOL 模型计算量相对较小。当输入参数的数量太大而无法进行计算成本较高的不确定性量化研究时,这使其成为一种理想的方法。
对于每个关注量,MOAT 方法计算每个输入参数的 MOAT 平均值和标准差,并将它们显示在 MOAT 散点图中。MOAT 平均值和标准差的排序给出了输入参数的相对重要性。前者越高意味着该参数显著影响关注量,后者越高则意味着该参数会产生影响,并且它要么与其他参数有很强的相互作用,要么具有非线性影响,或两者兼而有之。
灵敏度分析
灵敏度分析 研究类型用于计算关注量对输入参数的敏感性,包含两种方法:Sobol 法和相关法。
Sobol 法用于分析整个输入参数分布,并将每个关注量的方差分解为输入参数及其相互作用的贡献之和。
Sobol 法可以计算每个输入参数的 Sobol 指数。一阶 Sobol 指数显示由各个输入参数的方差产生的关注量的方差。总 Sobol 指数显示由每个输入参数的方差及其与其他输入参数的相互作用而产生的关注量的方差。每个关注量和所有参数的 Sobol 指数都显示在专用的 Sobol 图中,其中直方图按总 Sobol 指数排序。关注量对具有最高总 Sobol 指数的输入参数最敏感。输入参数的总 Sobol 指数和一阶 Sobol 指数之间的差异可以衡量该输入与其他输入之间相互作用的效果。
与筛选方法相比,灵敏度分析用于定量分析关注量的不确定性如何分配给不同的输入参数。这种方法需要更多的计算资源,原因是精确的 Sobol 指数计算依赖于高质量的代理模型。
相关法可以计算每个输入参数与关注量之间的线性和单调关系。对于基于相关法的灵敏度分析,可以计算四种类型的相关性:双变量相关、等级双变量相关、偏相关或等级偏相关。
不确定性传播
不确定性传播 研究类型用于分析输入参数的不确定性如何通过计算其概率密度函数(PDF)传播到每个关注量。大多数通过 COMSOL Multiphysics® 模型计算来将输入参数与关注量进行映射的底层物理应用通常是无法进行解析计算的。
因此,需要进行蒙特卡罗分析来近似得到 PDF。与 Sobol 法类似,使用代理模型可以显著降低蒙特卡罗分析的计算成本。对于每个关注量,都可以执行核密度估计(KDE),并为其生成可视化效果图,作为 PDF 的近似值。此外,置信区间表可以基于此分析给出每个关注量的平均值、标准差、最小值、最大值以及对应于 90%、95% 和 99% 置信水平的下限值和上限值。
可靠性分析
与用于探索关注量的总体不确定性的其他不确定性量化研究类型相比,可靠性分析,EGRA 方法可以解决更直接的问题。给定一个标称设计和一些特定的不确定输入,设计失败的概率是多少?失败的原因可以是设计的完全失效,也可以用质量准则来表述。
为了确保可靠性,传统的建模仿真方法是使用安全裕度和最坏情况。通过适当的可靠性分析,可以计算实际概率,从而避免高估和低估结果。从每个关注量的不确定性传播的置信区间表中可以得到粗略估计。但是通过可靠性分析,您可以根据关注量和相应阈值的组合来定义更复杂的可靠性准则。用于可靠性分析研究的高效全局可靠性分析(EGRA)方法可以有效地将计算资源引导至区分设计失败和成功的极限状态。
代理模型和响应面
使用 Sobol 法计算的灵敏度分析、不确定性传播和可靠性分析都依赖于精确的蒙特卡罗型分析,这通常需要大量的计算才能达到良好的准确性。对于 COMSOL Multiphysics® 模型计算可能需要大量资源且不确定性量化分析涉及多个参数的现实问题,仅使用 COMSOL Multiphysics® 模型计算进行的蒙特卡罗分析在计算上并不可行。“不确定性量化模块”的一个关键特点是,它能够针对特定的 UQ 分析训练和使用代理模型(也称为元模型),从而节省计算资源。
代理模型是一种紧凑的数学模型,用于表示和计算由输入参数定义的关注域中的关注量。这种模型完全独立于基础 COMSOL Multiphysics® 模型,并且在经过适当的训练后,可用于代替 COMSOL Multiphysics® 模型来预测其他输入参数值的关注量的值,而不是求解的值。构造代理模型的过程通常是自适应的,并且代理模型可以高精度地近似原始模型。用户定义的容差可以提高代理模型的准确性。更高级别的准确性需要额外的 COMSOL Multiphysics® 模型计算。
一旦建立代理模型,您就可以进行独立验证以进一步测试其有效性,并可以快速计算整个输入参数空间的响应面数据。然后,可以将响应面进行可视化,其中每次绘制一个关注量与两个输入参数的关系。
逆不确定性量化
当一些输入参数具有未知的概率分布时,使用逆不确定性量化(逆UQ ),称为校准参数。使用逆UQ,实验数据可以反向传播,以深入了解这些校准参数的统计特性。为了应用逆UQ,在进行分析之前,需要每个校准参数的先验概率分布。
实验数据通常可用于感兴趣的量和实验中使用的参数。还有一些无法直接测量的校准输入参数。例如,考虑一个实验,我们想要校准一个机械零件的杨氏模量。我们应该做一个实验,测量作为指定材料位移函数的拉伸应力。然后,应建立逆UQ研究,以使用实验数据和杨氏模量的先验知识来校准概率分布,该概率分布将最好地再现来自实验的拉伸应力的测量值。逆UQ可应用于各种基于物理的模型,包括与结构力学、流体流动、声学、热传递、电磁学和化学工程相关的模型。
为了使校准参数的后验概率分布的计算可行,代理模型与马尔可夫链蒙特卡罗(MCMC)方法一起使用。在计算之后,校准的输入参数的联合和边缘概率分布可以被可视化。此外,生成置信区间表,其提供诸如平均值的信息;标准偏差;最小值和最大值;并且对于每个校准的输入参数,下限值和上限值对应于90%、95%和99%的置信水平。