专利名称：一种数据提取的热分析建模方法
技术领域：
本发明属于热分析有限元建模领域，具体涉及一种数据提取的热分析建模方法。
背景技术：
有限元方法是在理论研究、工程应用中对特定对象进行热分析的常用方法，由于它的通用性和有效性，受到工程技术界的高度重视，广泛应用于国防、建筑等各个领域。自有限元方法出现以来，在国内外的理论研究和工程应用上，特定对象的热分析广泛采用此方法。同时，出现了大量的商用有限元软件，由于有限元软件建模、划分网格功能强大，解算能力强，可以快速地对分析各种对象、各种条件进行热分析，利用有限元软件进行热分析广泛应用于各个领域。然而，对于针对某些情况，例如当加载的边界条件较为复杂，随外界多种因素不断变化；以及考虑热分析模型与其他模型联合仿真时，采用有限元软件进行分析受到了一系列的限制，手动编写有限元程序仍然是十分必要的，例如在分析空间热环境对航天器的柔性部件的影响上，国内外均采用编写有限元程序进行分析。但是编写有限元程序具有以下缺点(a)对于复杂的模型建模、划分网格困难；(b)建模工作量大； (c)有限元编程过程复杂。由此可见，无论是利用有限元软件进行热分析还是编写有限元程序进行热分析，都存在一系列的限制，有必要探索新的热分析方法以解决现有方法的不足。
发明内容
针对现有技术中的问题，本发明提出一种数据提取的热分析建模方法，使在科学理论研究和实际工程应用中，更方便、更快捷地建立热分析模型。本发明基于有限元软件提取的数据进行热分析建模。由于某些情况下利用有限元软件进行分析存在诸多不便，因此需要通过有限元编程建立进行建模，其一般形式的有限元方程组如公式(1)所示。Cf+ KT = Fq(1)公式中C、Kc和F,分别代表组装的整体热容矩阵、热传导矩阵和热流载荷矩阵，T 为组装后的节点温度列阵。然而有限元建模过程复杂，工作量大，容易出错。观察公式(1)中的热分析有限元方程组，热传导矩阵、热容矩阵都是常系数矩阵，不随时间而变化，这些矩阵计算量较大，若从有限元软件中提取矩阵，可以节省计算量，减小建模出错的概率；对于方程中随时间变化的系数，可从有限元软件中提取模型的几何信息和网格信息，经过相关处理，基于有限元原理积分求得。有限元方程中的积分和微分方程求解分别采用高斯积分和开源数据库提供的函数提供。本发明提出一种数据提取的热分析建模方法，具体包括以下几个过程步骤一有限元软件中建模、划分网格。在有限元软件中建立要分析的有限元模型，有限元模型方程组的基本形式为公式 (1) ^^ +夂7 =巧，其中C代表组装的整体热容矩阵，K。为热传导矩阵，Ftl为热流载荷列阵，T为组装后的节点温度列阵，F,为随时间变化的系数。设置热分析单元类型，热分析单元类型为正六面体。设置密度、热容系数和热传导系数等材料参数，以及模型的几何尺寸，并进行网格划分，网格划分规则要求划分出的单元均为规则的正六面体。步骤二 基于有限元软件的信息计算随时间变化的系数。在有限元方程组一般形式的公式(1)中随时间变化的系数都需要通过有限元原理进行积分计算，例如热流载荷列阵Ftl为随时间变化的系数，这些系数往往是对模型的外表面或者某些表面积分。不随时间变化的系数可以从有限元软件进行提取。首先在有限元软件中输入单元信息文件和节点信息文件提取命令，提取模型的单元信息文件和节点信息文件。热流载荷列阵F,是对模型外表面进行积分计算得到的，因此需要筛选出模型的外表面，其他随时间变化的系数也一般是对模型的某些表面进行积分计算，在筛选出外表面后，进一步得到其这些系数的积分表面，通过有限元原理进行积分计算。单元信息文件中记录了每个单元的所有节点。由于所选择的单元为正六面体，可以列出每个单元所包含的6个表面，这些表面可由表面的顶点节点号来表示。拆分模型中所有单元得到模型的所有表面后，需要筛选出外表面。判定一个面是不是外表面的原则如下若一个面被两个单元所共有，则为内表面，若一个面仅属于一个单元，则为外表面。遍历模型的所有表面，筛选出顶点节点号不完全与其他面相同的表面，即得到模型的所有外表面。节点信息文件中记录了每个节点的三维坐标。得知三维坐标后，根据筛选出的外表面，可对方程公式(1)中的F,系数进行积分计算。步骤三基于有限元软件提取不随时间变化的常系数矩阵。从有限元软件中提取热分析有限元方程组中不随时间变化的常系数矩阵(方程 (1)中热容矩阵C，热传导矩阵K。为常系数矩阵)的相关文件(例如矩阵文件和节点顺序文件等)，解析文件，提取得到常系数矩阵。提取常系数矩阵的方法共有两种，分别为提取整体矩阵文件进行处理方法和提取单元矩阵文件进行处理方法。提取整体矩阵文件进行处理方法，在有限元软件中输入提取整体矩阵文件命令以及提取节点顺序文件命令，得到哈维尔-波音格式的整体矩阵文件和矩阵节点顺序文件。通过解析哈维尔-波音格式文件得到其代表的大型稀疏矩阵。解析出的大型稀疏矩阵是以特定的点顺序进行排列，对于有限元方程组的建立，通常需要把该矩阵按照节点的升序来排列。在矩阵节点顺序文件中所引矩阵的节点顺序，得到节点升序排列的整体矩阵即得到常系数矩阵。提取单元矩阵文件进行处理方法，在有限元软件中提取单元矩阵文件。单元矩阵文件中记录了每个单元的单元矩阵。提取单元信息文件，根据其记录的每个单元的节点号， 寻找单元矩阵中的元素在整体矩阵中的位置，进行有限元组集，生成节点升序排列的整体矩阵即得到常系数矩阵。步骤四配置求解环境。配置有限元方程组的解算环境，以完成有限元方程组的求解。通过编写微分方程组求解器或者利用开源的解算器所提供的微分方程组进行求解。步骤五设置初始条件及边界条件。设置要分析的有限元模型各节点的初始温度等初始参数，设置外界温度、外界热流量或热辐射边界条件。步骤六解算模型。解算热分析有限元方程组，求解模型中所有节点随时间变化的温度值。本发明具有以下优点(1)本发明提出一种数据提取的热分析建模方法，能够利用有限元软件强大的建模、划分网格功能，解决了复杂模型建模、划分网格困难的问题。(2)本发明提出一种数据提取的热分析建模方法能够减小有限元建模工作量；(3)本发明提出一种数据提取的热分析建模方法可便捷地加载复杂边界条件；(4)本发明提出一种数据提取的热分析建模方法支持与其他模型联合仿真。
图1是本发明提出的一种数据提取的热分析建模方法流程图；图2是本发明提出的提取整体矩阵文件处理得到整体矩阵流程图；图3是本发明提出的提取单元矩阵文件处理流程图；图4 实施例中某节点温度随时间变化的曲线。
具体实施例方式下面将结合附图和实例对本发明作进一步的详细说明。本发明提出一种数据提取的热分析建模方法，如图1所示，具体包括以下几个过程步骤一有限元软件中建模、划分网格。在有限元软件中建立要分析的有限元模型，(所述的有限元软件为ansys软件)， 有限元模型方程组的基本形式为公式(1) Cf=巧，其中C代表组装的整体热容矩阵， Kc为热传导矩阵，Ftl为热流载荷列阵，T为组装后的节点温度列阵，Ftl为随时间变化的系数。 设置热分析单元类型，热分析单元类型为正六面体。设置密度、热容系数和热传导系数等材料参数，以及模型的几何尺寸，并进行网格划分，网格划分规则要求划分出的单元均为规则的正六面体。步骤二 基于有限元软件的信息计算随时间变化的系数。在有限元方程组一般形式的公式(1)中随时间变化的系数都需要通过有限元原理进行积分计算，例如热流载荷列阵Fq为随时间变化的系数，这些系数往往是对模型的外表面或者某些表面积分。不随时间变化的系数可以从有限元软件进行提取。计算随时间变化的系数共需要三个步骤(1)从有限元软件中提取单元信息文件和节点信息文件；(2)根据单元信息文件筛选出模型的外表面；(3)对随时间变化的系数进行有限元积分计算。首先在有限元软件中输入单元信息文件和节点信息文件提取命令，提取模型的单元信息文件和节点信息文件。热流载荷列阵F,是对模型外表面进行积分计算得到的，因此需要筛选出模型的外表面，其他随时间变化的系数也一般是对模型的某些表面进行积分计算，在筛选出外表面后，进一步得到其这些系数的积分表面，通过有限元原理进行积分计算。单元信息文件中记录了每个单元的所有节点。由于所选择的单元为正六面体，可以列
6出每个单元所包含的6个表面，这些表面可由表面的顶点节点号来表示。拆分模型中所有单元得到模型的所有表面后，需要筛选出外表面。判定一个面是不是外表面的原则如下若一个面被两个单元所共有，则为内表面，若一个面仅属于一个单元，则为外表面。遍历模型的所有表面，筛选出顶点节点号不完全与其他面相同的表面，即得到模型的所有外表面。节点信息文件中记录了每个节点的三维坐标。得知三维坐标后，根据筛选出的外表面，可对方程公式(1)中的F,系数进行积分计算。步骤三基于有限元软件提取不随时间变化的常系数矩阵。从有限元软件中提取热分析有限元方程组中不随时间变化的常系数矩阵(方程 (1)中热容矩阵C，热传导矩阵K。为常系数矩阵)的相关文件(例如矩阵文件和节点顺序文件等)，解析文件，提取得到常系数矩阵。提取常系数矩阵的方法共有两种，分别为提取整体矩阵文件进行处理方法和提取单元矩阵文件进行处理方法。提取整体矩阵文件进行处理方法如图2所示，在有限元软件中输入提取整体矩阵文件命令以及提取节点顺序文件命令，得到哈维尔-波音格式的整体矩阵文件和矩阵节点顺序文件。矩阵哈维尔-波音格式文件的基本格式是前面有4或5行描述数据，其后为单位列矩阵元素值，矩阵哈维尔-波音格式文件的格式为第1行为文件头的字符型解释，如刚度矩阵或质量矩阵等标题。第2行表示该文件的总行数(不包括文件头即第1行)、矩阵列指针的总行数、 矩阵行索引的总行数、矩阵元素数值的总行数和右边项总行数。第3行表示矩阵类型、矩阵行数、矩阵列数、矩阵行索引数(对组装后的矩阵，该值等于矩阵行索引数)和单元元素数(对于组装后的矩阵此值为0)。第4行表示列指针格式、行索引格式、系数矩阵数值格式和右边项数值格式。第5行表示右边项格式、应用高斯起始矢量、求解矢量；两个整数分别表示右边项列数和行索引数。第6行后表示矩阵元素值(单列)，记载着矩阵中的非零元素。通过解析哈维尔-波音格式文件得到其代表的大型稀疏矩阵。解析出的大型稀疏矩阵是以特定的点顺序进行排列，对于有限元方程组的建立，通常需要把该矩阵按照节点的升序来排列。在矩阵节点顺序文件中所引矩阵的节点顺序，得到节点升序排列的整体矩阵即得到常系数矩阵。提取单元矩阵文件进行处理方法如图3所示，在有限元软件中提取单元矩阵文件。单元矩阵文件中记录了每个单元的单元矩阵。提取单元信息文件，根据其记录的每个单元的节点号，寻找单元矩阵中的元素在整体矩阵中的位置，进行有限元组集，生成节点升序排列的整体矩阵即得到常系数矩阵。步骤四配置求解环境。配置有限元方程组(即步骤一种的有限元模型的有限元方程组)的解算环境，以完成有限元方程组的求解。通过编写微分方程组求解器或者利用开源的解算器所提供的微分方程组进行求解。步骤五设置初始条件及边界条件。设置要分析的有限元模型各节点的初始温度等初始参数，设置外界温度、外界热流量或热辐射边界条件。步骤六解算模型。解算热分析有限元方程组，求解模型中所有节点随时间变化的温度值。实施例下面针对空间热环境对航天器柔件部件的热分析的实施例说明本发明。 航天器柔性部件在空间受到太阳辐射、地球反射太阳辐射以及地球红外热流的影响，同时柔性部件会向太空辐射大量的热量。由于空间热环境恶劣，导致柔性部件温度可能会变化剧烈，从而可能对航天器造成各种影响，而这种分析一般采用编写有限元程序进行建模，本实施例选取航天器太阳能帆板的简化模型——三维平板作为研究对象，由于加入了辐射边界条件，热分析的有限元方程组变为CT+ (Kc+Kr)T=Fq(2)公式中C，K。，&和F,分别代表组装的整体热容矩阵、热传导矩阵、热辐射矩阵和热流载荷矩阵，T为组装后的节点温度列阵。选取ANSYS软件为数据提取的有限元软件进行数据提取的热分析建模方法，包括以下几个步骤步骤一有限元软件中建模、划分网格。在ANSYS有限元软件中进行建模，热分析单元选择三维单元类型为S0LID70，所选单元为正六面体。材料参数设置如表1所示。设置三维平板的几何尺寸，长宽高分别为細、 Im和0. 15m。随后对模型进行网格划分，划分中的单元为规则的正六面体。表1平板材料参数
权利要求
1.一种数据提取的热分析建模方法，其特征在于具体包括以下几个过程步骤一有限元软件中建模、划分网格在有限元软件中建立要分析的有限元模型的有限元方程组形式，设置热分析单元类型，设置密度、热容系数和热传导系数，以及模型的几何尺寸，并进行网格划分；(1)步骤二基于有限元软件的信息计算随时间变化的系数从有限元软件中提取单元信息文件和节点信息文件；(2)根据单元信息文件筛选出模型的外表面；(3)对随时间变化的系数进行有限元积分计算；步骤三基于有限元软件提取不随时间变化的常系数矩阵从有限元软件中提取热分析有限元方程组中不随时间变化的常系数矩阵的相关文件， 解析文件，提取得到常系数矩阵；步骤四配置求解环境配置有限元方程组的解算环境，通过编写微分方程组求解器或者利用开源的解算器所提供的微分方程组进行求解；步骤五设置初始条件及边界条件设置要分析的有限元模型各节点的初始参数，设置外界温度、外界热流量或热辐射边界条件；步骤六解算模型解算热分析有限元方程组，求解模型中所有节点随时间变化的温度值。
2.根据权利要求1所述的一种数据提取的热分析建模方法，其特征在于所述的步骤一中有限元模型的方程组基本形式为0 + =巧，其中C代表组装的整体热容矩阵，K。为热传导矩阵，Fq为热流载荷矩阵，T为组装后的节点温度列阵，Fq为随时间变化的系数。
3.根据权利要求1所述的一种数据提取的热分析建模方法，其特征在于所述的步骤一中热分析单元类型为正六面体。
4.根据权利要求1所述的一种数据提取的热分析建模方法，其特征在于所述的步骤一中网格划分规则要求划分出的单元均为规则的正六面体。
5.根据权利要求1所述的一种数据提取的热分析建模方法，其特征在于所述的步骤三中不随时间变化的常系数矩阵的相关文件包括矩阵文件和节点顺序文件。
6.根据权利要求1所述的一种数据提取的热分析建模方法，其特征在于所述的步骤三中提取常系数矩阵的方法包括提取整体矩阵文件进行处理方法和提取单元矩阵文件进行处理方法。
7.根据权利要求6所述的一种数据提取的热分析建模方法，其特征在于所述的提取整体矩阵文件进行处理方法具体为在有限元软件中输入提取整体矩阵文件命令以及提取节点顺序文件命令，得到哈维尔-波音格式的整体矩阵和矩阵节点顺序文件，通过解析哈维尔-波音格式文件得到大型稀疏矩阵，根据节点顺序调整该矩阵顺序，得到节点升序排列的常系数矩阵。
8.根据权利要求6所述的一种数据提取的热分析建模方法，其特征在于所述的提取单元矩阵文件进行处理方法具体为在有限元软件中提取单元矩阵文件，单元矩阵文件中记录了每个单元的单元矩阵；提取单元信息文件，根据其记录的每个单元的节点号，寻找单元矩阵中的元素在整体矩阵中的位置，进行有限元组集，生成节点升序排列的整体矩阵，得到常系数矩阵。
全文摘要
本发明提出一种数据提取的热分析建模方法，具体包括步骤一有限元软件中建模、划分网格；步骤二基于有限元软件的信息计算随时间变化的系数；步骤三基于有限元软件提取不随时间变化的常系数矩阵；步骤四配置求解环境；步骤五设置初始条件及边界条件；步骤六解算模型。本发明的建模方法，能够利用有限元软件强大的建模、划分网格功能，解决了复杂模型建模、划分网格困难的问题，能够减小有限元建模工作量，并可便捷地加载复杂边界条件。
文档编号G06F17/50GK102156775SQ20111008514
公开日2011年8月17日 申请日期2011年4月6日 优先权日2011年4月6日
发明者王江云 申请人:北京航空航天大学