形态系数的计算

分类: 热解析 | 集成版 2026-04-06
CAE visualization for view factor theory - technical simulation diagram
形态系数的计算

形态系数计算的理论基础

概要

🙋

老师！今天是形态系数计算的话题是吧？那是什么东西呢？


🎓

面间的几何学放射交换系数。利用相反法则和总和法则进行计算。复杂形状采用Monte Carlo法。




支配方程




$$ F_{12} = \frac{1}{A_1}\int_{A_1}\int_{A_2}\frac{\cos\theta_1\cos\theta_2}{\pi r^2}dA_2\,dA_1 $$
$$ \sum_{j=1}^{N} F_{ij} = 1 $$



🙋

前辈说"形态系数计算的记述一定要好好做"的意思，现在理解了。


离散化手法

🙋

这个方程式，在计算机上实际怎么解呢？


🎓

用有限元法（FEM）进行空间离散化。组装要素刚性矩阵，构建全局刚性方程。


🎓

进行弱形式（变分形式）的转换，使用试验函数和形状函数，用Galerkin法进行定式化。要素类型的选择（低次要素 vs. 高次要素、完全积分 vs. 低减积分）直接关系到解的精度与计算成本的权衡。




矩阵求解算法

🙋

矩阵求解算法，具体是怎么回事呢？


🎓

直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程。大规模问题采用有前处理的迭代法效果显著。



求解法分类内存使用量适用规模
LU分解直接法O(n²)小～中规模
Cholesky分解直接法（对称正定值）O(n²)小～中规模
PCG法迭代法O(n)大规模
GMRES法迭代法O(n·m)大规模·非对称
AMG前处理前处理O(n)超大规模
🙋

也就是说，在有限元法的地方偷工减料的话，后来会吃苦头吧。铭记于心！


商用工具的实现

🙋

那么做形态系数计算的话，用什么样的软件呢？


工具名开发方/现属主要文件格式
Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
COMSOL MultiphysicsCOMSOL AB.mph
Ansys Mechanical (旧ANSYS Structural)Ansys Inc..cdb, .rst, .db, .ans, .mac
Abaqus FEA (SIMULIA)Dassault Systèmes SIMULIA.inp, .odb, .cae, .sta, .msg

厂商系谱与产品整合的经过

🙋

各个软件的来历，结果多少有点戏剧性呢？



Ansys Fluent

🙋

接下来是Ansys Fluent的话题吧。是什么内容呢？


🎓

由Fluent Inc.开发。2006年被Ansys收购。非结构网格基础的通用CFD求解器。

现属: Ansys Inc.



Simcenter STAR-CCM+

🙋

接下来是Simcenter STAR的话题吧。是什么内容呢？


🎓

由CD-adapco开发。2016年被西门子收购，整合到Simcenter品牌。多面体网格为特色。

现属: Siemens Digital Industries Software


🙋

听到这里，终于明白了开发为什么那么重要！



COMSOL Multiphysics

🙋

"COMSOL Multiphysics"请教一下！


🎓

1986年在瑞典成立。作为与MATLAB联动的FEMLAB开始，后改名为COMSOL。擅长多物理场。

现属: COMSOL AB


🙋

哦～，开发的故事，超有意思！请告诉我更多。


文件格式与互操作性

🙋

在不同软件之间交换数据时有什么注意事项吗？


格式扩展名种类概要
STEP.stp/.step中立CADISO 10303准拟的3D CAD数据交换格式。支持形状+PMI。
CGNS.cgnsCFD数据CFD General Notation System。CFD结果的标准交换格式。
VTK.vtk/.vtu可视化Visualization Toolkit格式。用于ParaView等。
🎓

在不同求解器间转换模型时，需要注意要素类型的对应关系、材料模型的兼容性、荷载·边界条件的表示差异。特别是高次要素或特殊要素（内聚要素、用户定义要素等）往往无法在求解器间直接转换。


🙋

原来如此……格式看似简单，实际上深度非常大呢。


实务上的注意事项

🙋

教科书里没有的"现场知识"之类的东西有吗？


🎓

网格收敛性确认、边界条件妥当性验证、材料参数敏感性分析非常重要。


🎓
  • 网格依赖性验证：用3水平以上的网格密度确认收敛性
  • 边界条件妥当性：设置物理上有意义的约束条件
  • 结果验证：与理论解、实验数据、已知基准问题的比较


    • 🙋

    形态系数计算的整体情况掌握了！明天开始在实务中加以意识。


    🎓

    嗯，不错！实际动手实践是最好的学习。有不懂的随时问我。


    Coffee Break 闲谈

    形态系数的定义与物理意义

    形态系数（View Factor）Fij表示从面i放射的扩散能中到达面j的比率。也称角系数、configuration factor，是0≦Fij≦1的无次元数。相互关系法则Ai·Fij＝Aj·Fji（Hottel-Sparrow，1960年代）和总和法则ΣjFij＝1是形态系数的根本约束，CAE代码是否满足这些条件可以用来检验结果的物理妥当性。仅由几何性质决定，与材料特性无关是其特征。

    形态系数计算的数值计算手法

    数值手法的详细内容

    🙋

    具体用什么样的算法来解形态系数计算呢？




    离散化的定式化



    🎓

    使用形状函数 $N_i$ 来近似未知量：



    $$ u^h(\mathbf{x}) = \sum_{i=1}^{n} N_i(\mathbf{x}) \, u_i $$




    🎓

    用式子表示如下。


    $$ K_e = \int_{\Omega_e} B^T \, D \, B \, d\Omega \approx \sum_{g=1}^{n_g} w_g \, B^T(\xi_g) \, D \, B(\xi_g) \, |J(\xi_g)| $$

    基础方程式的离散形式


    🎓

    用式子表示如下。


    $$ F_{12} = \frac{1}{A_1}\int_{A_1}\int_{A_2}\frac{\cos\theta_1\cos\theta_2}{\pi r^2}dA_2\,dA_1 $$
    $$ \sum_{j=1}^{N} F_{ij} = 1 $$

    🙋

    嗯～，光看式子不容易理解…… 表示什么呢？


    🎓

    将连续体的支配方程离散化，得到下列代数方程组：



    $$ [K]\{u\} = \{F\} $$


    🎓

    这里 $[K]$ 是全局刚性矩阵（或同等的系统矩阵），$\{u\}$ 是未知节点变量向量，$\{F\}$ 是外力向量。


    🙋

    啊，原来如此！连续体的支配方程……原来是这么个机制啊。


    要素技术

    🙋

    "要素技术"听过，但可能理解得不够深……


    要素类型次数节点数(3D)精度计算成本
    四面体1次线性4低（剪切锁定）低
    四面体2次二次10高中
    六面体1次线性8中中
    六面体2次二次20非常高高
    棱柱线性/二次6/15中～高中

    积分方案

    🙋

    积分方案，具体是怎么回事？


    🎓
    • 完全积分：所有项进行精确积分。刚性过大评估的倾向（锁定）
    • 低减积分：削减积分点数。计算效率提高，但有沙漏模式发生的风险
    • 选择性低减积分 (B-bar法)：体积项和偏差项分别积分。回避锁定

      • 🙋

      听到这里，终于明白了要素类型为什么那么重要！


      收敛性与稳定性

      🙋

      不收敛的话，首先要检查什么？


      🎓
      • h-细化：细分网格（减小要素尺寸 h）来提高精度
      • p-细化：提高要素的多项式次数来提高精度
      • hp-细化：同时优化 h 和 p

        • 🎓

        收敛速度：二次要素为 $O(h^2)$ 的阶数误差减少（光滑解的情况）


        🙋

        原来……细分网格看似简单，实际上深度非常大呢。


        求解器设置的建议

        🙋

        具体用什么样的算法来解形态系数计算呢？


        参数推荐值备注
        迭代法收敛判定$10^{-6}$残差范数基准
        前处理手法ILU(0) or AMG根据问题规模
        最大迭代次数1000不收敛时需重新设定
        内存模式In-core尽可能

        线性要素 vs 2次要素

        热传导分析中线性要素往往可以获得充分精度。温度梯度急剧的区域（热冲击等）推荐2次要素。

        热流束的评估

        从要素内的温度梯度计算。需要平滑化处理，如同节点应力。

        对流-扩散问题

        当Peclet数高（对流支配）时需要风上稳定化（SUPG等）。纯热传导问题无需。

        非定常分析的时间步

        热扩散的特征时间 $\tau = L^2 / \alpha$（$\alpha$: 热扩散率）相比要充分小。对急剧温度变化采用自动时间步控制有效。

        非线性收敛

        温度依赖物性导致的非线性往往较温和，Picard迭代（直接代换法）多数情况下就足够。放射的强非线性情况推荐Newton法。

        定常分析的判定

        全节点温度变化在阈值以下（$|\Delta T| / T_{max} < 10^{-5}$等）时判定为收敛。

        形态系数计算的实务应用

        实践指南

        🙋

        老师，请教一下"实践指南"！


        🎓

        讲解形态系数计算的实务性分析流程和注意事项。


        🙋

        老师的说明好懂！形态系数计算实务的困惑消散了。


        分析流程

        🙋

        从第一步开始请教！首先要做什么？


        🎓

        1. 前处理 (Pre-processing)

        • CAD数据的导入和形状简化
        • 材料特性的定义
        • 网格生成（要素类型·尺寸的决定）
        • 边界条件和荷载条件的设置

        🎓

        2. 求解 (Solving)

        • 求解器设置（解法、收敛基准、输出控制）
        • 作业投入和计算执行
        • 收敛监视

        🎓

        3. 后处理 (Post-processing)

        • 结果可视化（位移、应力、其他物理量）
        • 结果验证和妥当性确认
        • 报告制作


        网格生成的最佳实践

        🙋

        网格的好坏怎么判断？



        要素品质指标

        🙋

        "要素品质指标"请教一下！


        指标理想值允许范围影响
        纵横比1.0< 5.0精度下降
        Jacobian比1.0> 0.3要素退化
        翘曲0°< 15°精度下降
        歪度0°< 45°收敛性恶化
        锥形比0< 0.5精度下降

        网格密度的决定

        🙋

        网格密度的决定，具体是怎么回事？


        🎓
        • 应力集中部：配置最少3层以上的要素
        • 应力梯度大的区域：要素尺寸缩小到周围的1/3～1/5
        • 荷载施加点附近：局部细化
        • 远方区域：粗网格以确保计算效率


          • 边界条件设置指南

          🙋

          边界条件，听说这儿要是搞错了全部完蛋……


          🎓
          • 注意过约束：刚体运动的约束只需6自由度
          • 活用对称条件：削减计算规模
          • 荷载的等价分配：集中荷载 vs. 分布荷载的选择

            • 🙋

            啊，原来如此！过约束注意……原来是这么个机制啊。


            商用工具的实现步骤

            🙋

            各种软件，各有特点吧？具体分别介绍一下特征！


            工具名开发方/现属主要文件格式
            Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
            Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
            COMSOL MultiphysicsCOMSOL AB.mph
            Ansys Mechanical (旧ANSYS Structural)Ansys Inc..cdb, .rst, .db, .ans, .mac
            Abaqus FEA (SIMULIA)Dassault Systèmes SIMULIA.inp, .odb, .cae, .sta, .msg

            Ansys Fluent

            🙋

            接下来是Ansys Fluent的话题吧。是什么内容呢？


            🎓

            由Fluent Inc.开发。2006年被Ansys收购。非结构网格基础的通用CFD求解器。

            现属: Ansys Inc.



            Simcenter STAR-CCM+

            🙋

            接下来是Simcenter STAR的话题吧。是什么内容呢？


            🎓

            由CD-adapco开发。2016年被西门子收购，整合到Simcenter品牌。多面体网格为特色。

            现属: Siemens Digital Industries Software


            🙋

            听到这里，终于明白了开发为什么那么重要！


            常见失败与对策

            🙋

            初学者容易犯什么样的失败？希望事前知道！


            症状原因对策
            计算不收敛网格品质不良、不适切的边界条件网格改善、拘束条件重新检视
            应力异常大应力特异点、网格依赖特异点回避、局部网格细化
            位移非现实材料常数误差、单位系不一致输入数据确认
            计算时间过长不必要的细化、低效的解法网格最优化、并行计算

            品质保证检查清单

            🙋

            教科书里没有的"现场知识"之类的东西有吗？


            🎓
            • 网格收敛性用3水平以上确认了吗
            • 力的平衡（反力合计）验证了吗
            • 结果在物理上的妥当范围内吗确认了吗
            • 与已知理论解或基准问题比较了吗


              • 🙋

              形态系数计算的整体情况掌握了！明天开始在实务中加以意识。


              🎓

              嗯，不错！实际动手实践是最好的学习。有不懂的随时问我。


              Coffee Break 闲谈

              原子炉格纳容器的热设计

              加压水型原子炉（PWR）的格纳容器内部热设计中，炉心·蒸汽发生器·格纳容器壁间的形态系数是衰变热除去（余热除去）设计的根据。三菱重工制APWR的初期安全解析（2010年代）中构建了面数3000超的3D Enclosure，进行LOCA（冷却材丧失事故）后的自然对流＋放射连成解析，评价了容器最高温度。形态系数的误差直接联系到事故时温度预测，因此算法的验证采用上述Hamilton & Morgan解析解的照合是强制的。

              形态系数计算的软件比较

              商用工具比较

              🙋

              各种软件，各有特点吧？具体分别介绍一下特征！


              🎓

              对形态系数计算进行支持的主要商用CAE工具的功能比较，以及各产品的历史背景进行详细记述。


              🙋

              老师的说明好懂！形态系数计算支持的困惑消散了。


              支持工具列表

              🙋

              那么做形态系数计算的话，用什么样的软件呢？


              工具名开发方/现属主要文件格式
              Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
              Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
              COMSOL MultiphysicsCOMSOL AB.mph
              Ansys Mechanical (旧ANSYS Structural)Ansys Inc..cdb, .rst, .db, .ans, .mac
              Abaqus FEA (SIMULIA)Dassault Systèmes SIMULIA.inp, .odb, .cae, .sta, .msg

              Ansys Fluent

              🙋

              接下来是Ansys Fluent的话题吧。是什么内容呢？


              🎓

              由Fluent Inc.开发。2006年被Ansys收购。非结构网格基础的通用CFD求解器。

              现属: Ansys Inc.



              Simcenter STAR-CCM+

              🙋

              接下来是Simcenter STAR的话题吧。是什么内容呢？


              🎓

              由CD-adapco开发。2016年被西门子收购，整合到Simcenter品牌。多面体网格为特色。

              现属: Siemens Digital Industries Software


              🙋

              听到这里，终于明白了开发为什么那么重要！



              COMSOL Multiphysics

              🙋

              "COMSOL Multiphysics"请教一下！


              🎓

              1986年在瑞典成立。作为与MATLAB联动的FEMLAB开始，后改名为COMSOL。擅长多物理场。

              现属: COMSOL AB



              Ansys Mechanical (旧ANSYS Structural)

              🙋

              "Ansys Mechanical"请教一下！


              🎓

              1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（Ansys Parametric Design Language）。

              现属: Ansys Inc.


              🙋

              啊，原来如此！开发……原来是这么个机制啊。


              功能比较矩阵

              🙋

              预算和时间都有限，最划算的是哪个？


              功能FluentStar-CCM+COMSOLAnsys Mechanical
              基本功能○○○○
              高级功能○○○△
              自动化/脚本○○○○
              并行计算○○○○
              GPU支持△△△○

              转换时的风险

              🙋

              转换时的风险，具体是怎么回事？


              🎓
              • 要素类型的不兼容：求解器固有要素无法用中立格式表示
              • 材料模型的差异：虽然同名但内部实装不同的情况
              • 边界条件的重新定义：大多数情况下需要手动重新设置
              • 结果数据的比较：输出变量的定义（节点值 vs. 要素值、积分点值）有差异

                • 🙋

                啊，原来如此！不同工具间的模……原来是这么个机制啊。


                许可证形式

                🙋

                "许可证形式"听过，但可能理解得不够深……


                工具许可证特色
                商用FEA节点锁定/浮动高额但附带官方支持
                OpenFOAMGPL免费但支持需付费
                COMSOL节点锁定/浮动按模块购买
                Code_AsterGPLEDF开发的开源求解器

                选择指南

                🙋

                最后到底选哪个，判断标准教我一下？


                🎓

                形态系数计算的工具选择时需要考虑以下方面：


                🎓
                • 分析规模：针对数万～数亿DOF的可扩展性
                • 物理模型：必要的构成则·要素类型的支持情况
                • 工作流：与CAD的联动、自动化的容易程度
                • 成本：初期投资 + 年度保修 + 教育成本
                • 支持：技术支持的质量和响应速度


                  • 🙋

                  形态系数计算的整体情况掌握了！明天开始在实务中加以意识。


                  🎓

                  嗯，不错！实际动手实践是最好的学习。有不懂的随时问我。


                  Coffee Break 闲谈

                  按求解器分的视角因子计算法

                  形态系数计算算法因求解器而异，各具特色。ANSYS Mechanical 2024R1用Hemicube法高速，Abaqus/Standard R2024用Ray Shooting对复杂形状强。MSC Nastran SOL153是解析积分（仅简单形状）和数值积分的混合。OpenFOAM的viewFactors是Hemi-Cube法的开源实装，定制改造自由度高。选择基准是面数·形状复杂度·精度要求的3轴，是实务上的常识。

                  形态系数计算的先端研究

                  先端话题与研究动向

                  🙋

                  形态系数计算的领域，将来如何演进呢？


                  🎓

                  见一下形态系数计算领域的最新研究动向和先进手法。



                  最新的数值手法

                  🙋

                  接下来是最新数值手法的话题吧。是什么内容呢？



                  🙋

                  嗯～，光看式子不容易理解…… 表示什么呢？


                  🎓
                  • 等几何分析 (IGA)：直接使用NURBS基函数，实现CAD-CAE间的无缝联动
                  • 粒子法 (SPH, MPM)：网格自由手法追踪大变形·破坏
                  • 位相场法 (Phase-Field)：通过界面的隐式表示进行复杂界面追踪
                  • 机器学习支援：代理模型、物理信息神经网络 (PINN)


                    • 高性能计算 (HPC) 的对应


                    并行化手法概要适用求解器
                    MPI (领域分割)分布式内存型。大规模问题的标准所有主要求解器
                    OpenMP共享内存型。节点内并行许多求解器
                    GPU (CUDA/OpenCL)GPGPU活用。特别对显式法有效LS-DYNA, Fluent等
                    混合 MPI+OpenMP节点间+节点内并行大规模HPC环境

                    形态系数计算的故障排除

                    故障排除




                    常见错误与对策

                    🙋

                    老师也在形态系数计算上熬过夜调试吗？（笑）



                    1. 收敛失败

                    🙋

                    收敛失败，具体是怎么回事？


                    🎓

                    症状：求解器在指定迭代次数内不收敛，异常终止


                    🎓

                    可能的原因：

                    • 网格品质不足（过度扭曲的要素）
                    • 材料参数设置不当
                    • 不适切的初始条件
                    • 非线性性过强（荷载步不足）

                    🎓

                    对策：

                    • 进行网格品质检查（纵横比、Jacobian）
                    • 确认材料参数的单位系
                    • 荷载分割为多个步骤（增加子步数）
                    • 放宽收敛判定基准（但要注意精度）

                    🙋

                    也就是说，在收敛失败的地方偷工减料的话，后来会吃苦头吧。铭记于心！



                    2. 非物理的结果