一种求解不确定热传导问题的随机正交展开方法技术

本发明专利技术公开了一种求解不确定热传导问题的随机正交展开方法，步骤如下：引入随机变量对热传导问题中的不确定参数进行定量化表示；结合随机变量建立热传导问题的随机微分控制方程；根据随机变量分布类型选用正交多项式基底函数，将随机温度响应进行正交展开；给定每个随机变量的配点数量，利用张量积法则构造整个不确定空间的配点集合；计算所有配点处的温度响应，利用矩阵的广义逆求得温度响应正交展开式中的各项系数；根据基底函数的正交关系计算随机温度响应的均值和标准差。本发明专利技术可系统化解决含有随机不确定参数的热传导问题，进一步提高了随机不确定分析方法的计算精度，这是一般商用软件所不能实现的。

【技术实现步骤摘要】
一种求解不确定热传导问题的随机正交展开方法
本专利技术属于机械工程领域，具体涉及一种求解不确定热传导问题的随机正交展开方法。
技术介绍
在对工程系统进行分析和设计时，要求能正确地计算系统响应，以便保证工程系统满足某些指标的要求。但是，由于人们认识客观世界水平和手段的限制，实际工程中经常存在着与材料属性、外力载荷、初始条件、边界约束以及加工装配有关的不确定性。对复杂系统而言，即使很小的不确定因素，通过各子系统之间的传播和扩散，也可能会对最终的响应输出产生明显的扰动。因此，研究这些不确定性对系统响应的影响有着广泛的工程背景和重要的学术价值。热分析在工程中普遍存在，尤其是在航空航天领域，其核心问题是先要确定结构的温度场。传热模型输入参数的不确定性必然会导致结构温度场的不确定性特征。传统的安全因子法仅仅根据工程师的经验，在计算过程中对传热模型进行不同程度的修正，以完善不确定性带来的不足。安全因子的大小是根据经验粗略制定的，这使得计算结果非常粗糙，不能满足精细化要求。而随机模型是将系统中的不确定性看作随机变量或随机过程，进而利用概率论和统计方法研究不确定性传播规律。随机建模及数值计算方法具有成熟的理论基础，在不确定分析领域发挥了重要作用。目前，将随机分析理论与有限元计算方法相结合衍生出来的随机有限元法在不确定结构的静、动力特性分析方面已经取得了不少研究成果，但在传热领域的应用还十分有限。另外，传统的随机模拟方法，通过数字模拟和统计分析来求取系统响应的概率特征。尽管操作简单，但其计算精度依赖于大量的抽样实验，因此很难应用于复杂的工程系统。随机摄动方法计算量小，但由于忽略了部分高阶项，计算精度难以满足工程需求。因此，如何建立准确高效的随机分析方法对不确定传热问题进行数值求解，是目前学术领域的一个研究热点，对于弥补现有传热数值计算方法的不足，具有重要的工程应用价值。
技术实现思路
本专利技术所要解决的技术问题为：克服现有技术在热传导问题求解中存在的不足，基于正交多项式展开理论和配点分析技术，提出了一种有效预测结构温度场概率特征的数值计算方法，可系统化解决含有随机不确定参数的热传导问题，在保证计算效率的同时，进一步提高了随机分析方法的计算精度。本专利技术为解决上述技术问题采用的技术方案为：一种求解不确定热传导问题的随机正交展开方法，包括以下步骤：步骤一：引入随机变量对热传导问题中的不确定参数进行定量化表示；步骤二：结合步骤一中引入的随机变量，建立热传导问题的随机微分控制方程；步骤三：根据步骤一中随机变量的分布类型选用正交多项式基底函数，将步骤二随机微分控制方程中涉及的温度响应进行正交展开，得到随机温度响应的正交展开式；步骤四：给定每个随机变量所对应的配点数量，利用张量积法则构造整个不确定空间的配点集合；步骤五：利用现有软件或程序计算步骤四配点集合中所有配点处的温度响应，建立关于步骤三随机温度响应正交展开式中各项系数的线性方程组，利用矩阵的广义逆对此线性方程组进行求解，得到各项系数的一组值；步骤六：将步骤五中得到的各项系数的一组值代回到步骤三随机温度响应的正交展开式中，根据基底函数的正交关系，计算随机温度响应的均值和标准差。其中，所述步骤三中利用正交多项式对随机微分控制方程中涉及的温度响应进行正交展开时，正交多项式的类型和截断阶数并不是固定不变的，根据随机变量分布类型和逼近精度要求进行选取，例如高斯分布的随机变量对应埃尔米特正交多项式，均匀分布的随机变量对应勒让德正交多项式，另外正交多项式截断阶数越高，逼近精度就越高。其中，所述步骤四中配点集合的建立并不是固定不变的，根据计算耗费和计算精度的要求来选取每个随机变量所对应的配点数量，配点数量越多，计算精度就越高，而计算耗费就越大。上述各步骤具体包括以下过程：步骤一：引入n个随机变量ξ1,ξ2,...,ξn对热传导问题中的不确定参数进行定量化表示，并将其统一记为向量的形式ξ＝(ξ1,ξ2,...,ξn)。步骤二：结合步骤一中引入的随机变量，建立热传导问题的随机微分控制方程：其中x为物理坐标，T为温度响应，k为材料热传导系数，f表示系统的热源强度。步骤三：根据步骤一中随机变量的分布类型选用合适的正交多项式基底函数，将步骤二随机微分控制方程中涉及的温度响应T(x；ξ)进行正交展开，得到随机温度响应的正交展开式：其中Φi(ξ)为事先选定的正交多项式基底函数，Ti(x)为对应的各项系数，i＝(i1,i2,...,in)且满足|i|＝i1+i2+...+in，N为正交多项式的截断阶数。上述正交多项式中展开项的个数可用随机变量个数n和截断阶数N计算n为随机变量的个数。步骤四：给定每个随机变量所对应的配点数量，利用张量积法则构造整个不确定空间的配点集合。首先，对于随机变量ξi而言，确定其分布区间其中ξi和表示此分布区间的下界和上界。其次，给定配点数量mi，则在区间中各个配点的具体位置为：其中和称作区间的中点和半径。然后，用点集表示随机变量ξi在分布区间内所有配点组成的集合，那么对于n个随机变量组成的整个不确定空间而言，直接利用张量积法则可得配点集合Θ：而配点集合Θ中的配点总数M为：在此基础上，将配点集合Θ改写为的形式，用来表示整个不确定空间中所有的配点其中上标node为配点符号。步骤五：利用现有软件或程序计算步骤四配点集合中所有配点处的温度响应，建立关于步骤三随机温度响应正交展开式中各项系数的线性方程组，利用矩阵的广义逆对此线性方程组进行求解，得到各项系数的一组值。首先，步骤二中的随机微分控制方程在配点处可表示为：其次，利用现有软件或程序对上述方程进行求解，得到所有配点处的温度响应然后，基于步骤三中随机温度响应的正交展开式，建立关于各项系数Ti(x)的线性方程组：紧接着，利用矩阵的广义逆对此线性方程组进行求解，得到各项系数Ti(x)的一组值。步骤六：将步骤五中得到的各项系数Ti(x)的一组值代回到步骤三随机温度响应的正交展开式中，根据基底函数的正交关系，最终可以得到随机温度响应T(x；ξ)的均值E[T(x；ξ)]和标准差σ[T(x；ξ)]：E[T(x；ξ)]≈E[TN(x；ξ)]＝T0(x)其中γi为表征基底函数正交关系的归一化因子。本专利技术与现有技术相比的优点在于：(1)与传统的传热数值计算方法相比，所提出的随机正交展开方法充分计及传热模型的不确定因素，计算结果对温度场分析具有更重要的指导意义。(2)利用随机正交展开方法对温度响应进行近似表示，可有效提高逼近精度。同时，利用基底函数的正交关系，可快速得到随机温度响应均值、标准差等概率特征。(3)基于配点理论对随机温度响应正交展开式的各项系数进行求解，可以充分利用原有确定性模型的计算程序，操作简单，实施方便。附图说明图1为本专利技术的一种求解不确定热传导问题的随机正交展开方法流程图；图2为本专利技术的二维平板传热结构模型示意图。具体实施方式下面结合附图和实施例对本专利技术做进一步说明。本专利技术适用于含有随机不确定参数的热传导问题的温度场预测。本专利技术实施方式以某二维平板传热结构模型为例，具体说明所述的一种求解不确定热传导问题的随机正交展开方法。另外，此二维平板传热结构模型的温度响应随机正交展开方法可以推广到其他含有随机不确定参数的热传导问题温度场预测中。一种求本文档来自技高网...

【技术保护点】
一种求解不确定热传导问题的随机正交展开方法，其特征在于包括以下步骤：步骤一：引入随机变量对热传导问题中的不确定参数进行定量化表示；步骤二：结合步骤一中引入的随机变量，建立热传导问题的随机微分控制方程；步骤三：根据步骤一中随机变量的分布类型选用正交多项式基底函数，将步骤二随机微分控制方程中涉及的温度响应进行正交展开，得到随机温度响应的正交展开式；步骤四：给定每个随机变量所对应的配点数量，利用张量积法则构造整个不确定空间的配点集合；步骤五：计算步骤四配点集合中所有配点处的温度响应，建立关于步骤三随机温度响应正交展开式中各项系数的线性方程组，利用矩阵的广义逆对此线性方程组进行求解，得到各项系数的一组值；步骤六：将步骤五中得到的各项系数的一组值代回到步骤三随机温度响应的正交展开式中，根据基底函数的正交关系，计算随机温度响应的均值和标准差。

【技术特征摘要】
1.一种求解不确定热传导问题的随机正交展开方法，其特征在于包括以下步骤：步骤一：引入n个随机变量ξ1,ξ2,...,ξn对热传导问题中的不确定参数进行定量化表示，并将其统一记为向量的形式ξ＝(ξ1,ξ2,...,ξn)，步骤二：结合步骤一中引入的随机变量，建立热传导问题的随机微分控制方程：其中x为物理坐标，T为温度响应，k为材料热传导系数，f表示系统的热源强度；步骤三：根据步骤一中随机变量的分布类型选用合适的正交多项式基底函数，将步骤二随机微分控制方程中涉及的温度响应T(x；ξ)进行正交展开，得到随机温度响应的正交展开式：其中Φi(ξ)为事先选定的正交多项式基底函数，Ti(x)为对应的各项系数，i＝(i1,i2,...,in)且满足|i|＝i1+i2+...+in，N为正交多项式的截断阶数；上述正交多项式中展开项的个数可用随机变量个数n和截断阶数N计算n为随机变量的个数；步骤四：给定每个随机变量所对应的配点数量，利用张量积法则构造整个不确定空间的配点集合，首先，对于随机变量ξi而言，确定其分布区间其中ξi和表示此分布区间的下界和上界；其次，给定配点数量mi，则在区间中各个配点的具体位置为：其中和称作区间的中点和半径；然后，用点集表示随机变量ξi在分布区间内所有配点组成的集合，那么对于n个随机变量组成的整个不确定空间而言，直接利用张量积法则可得配点集合Θ：而配点集合Θ中的配点总数M为：在此基础上，将配点集合Θ改写为的形式，用来表示整个不确定空间中所有的配点其中上标node为配点符号；步骤五：计算步骤四配点集合中所有配点处的温度响应，建立关于步骤三随机温度响应正交展开式中各项...

【专利技术属性】
技术研发人员：邱志平，王冲，王晓军，许孟辉，李云龙，陈贤佳，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11