CN107391786A - 适用于卫星有限元热变形分析的温度场快速赋值方法 - Google Patents

适用于卫星有限元热变形分析的温度场快速赋值方法 Download PDF

Info

Publication number
CN107391786A
CN107391786A CN201710419042.2A CN201710419042A CN107391786A CN 107391786 A CN107391786 A CN 107391786A CN 201710419042 A CN201710419042 A CN 201710419042A CN 107391786 A CN107391786 A CN 107391786A
Authority
CN
China
Prior art keywords
temperature field
thermal deformation
satellite
finite element
deformation analysis
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710419042.2A
Other languages
English (en)
Inventor
苏若斌
孔祥森
周静
陶炯鸣
赵川
李�昊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Satellite Engineering
Original Assignee
Shanghai Institute of Satellite Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Satellite Engineering filed Critical Shanghai Institute of Satellite Engineering
Priority to CN201710419042.2A priority Critical patent/CN107391786A/zh
Publication of CN107391786A publication Critical patent/CN107391786A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种适用于卫星有限元热变形分析的温度场快速赋值方法，其包括以下步骤：步骤一，建立需要进行热变形分析的卫星有限元模型，包括网格划分、多点约束定义、材料定义，导出模型文件；步骤二，对温度场测试数据进行坐标转换，使测试数据坐标与有限元模型坐标一致；步骤三，将模型文件导入NX软件，根据热变形分析方案对模型进行分组设置等。本发明能够在短时间内完成结构热变形分析，进而指导设计生产，缩短卫星研制周期。

Description

适用于卫星有限元热变形分析的温度场快速赋值方法
技术领域
本发明涉及一种赋值方法，特别是涉及一种适用于卫星有限元热变形分析的温度场快速赋值方法。
背景技术
适用于卫星有限元热变形分析的温度场快速赋值方法直接关乎到卫星结构设计的合理性、有效载荷安装面的精度稳定性，以及为卫星在轨热变形有效控制提供数据参考，进而影响卫星技术指标能否实现。
传统的结构有限元温度场赋值方法为，根据结构温度实测数据，选中对应结构有限元模型的一组节点，对这些节点进行温度赋值，再选中另一组节点，赋对应温度值，如此一组组选中分别手动赋值，直至对所需计算的模型全部加载上温度场的数值。由于有限元模型节点多，采用此方法工作量大，耗时长，并且容易存在疏漏或加载节点存在位置偏差，导致模型温度设置不精准，难以分析出卫星结构的热变形情况，进而无法指导卫星在轨热变形的控制。
发明内容
本发明所要解决的技术问题是提供一种适用于卫星有限元热变形分析的温度场快速赋值方法，其能够在短时间内完成结构热变形分析，进而指导设计生产，缩短卫星研制周期。
本发明是通过下述技术方案来解决上述技术问题的：一种适用于卫星有限元热变形分析的温度场快速赋值方法，其包括以下步骤：
步骤一，建立需要进行热变形分析的卫星有限元模型，包括网格划分、多点约束定义、材料定义，导出模型文件；
步骤二，对温度场测试数据进行坐标转换，使测试数据坐标与有限元模型坐标一致；
步骤三，将模型文件导入NX(代号)软件，根据热变形分析方案对模型进行分组设置；
步骤四，对NX模型文件中的分组有限元分别进行温度场加载；
步骤五，将NX模型文件保存并导出有限元格式文件；
步骤六，在MSC Patran(工业领域最著名的有限元前、后处理器，是一个开放式、多功能的三维MCAE软件包，具有集工程设计、工程分析和结果评估功能于一体的、交互图形界面的CAE集成环境)中导入温度场模型文件用于热变形分析可获得分析数据。
优选地，所述步骤一是利用MSC Patran有限元软件建立需要进行热变形分析的卫星有限元模型。
优选地，所述步骤一中导出模型文件为为bdf(基于ASCLL的在windows系统和mac系统通用的字体格式)格式的Nastran(大型应用有限元程序)模型文件，获得包含卫星结构所有热相关特征信息。
优选地，所述步骤二是利用MATLAB(商业数学软件)编写的坐标转换程序对卫星温度场实测数据进行坐标转换。
优选地，所述步骤二中的坐标转换的文件均为excel(电子表格软件)编制的csv(Comma-Separated Values，逗号分隔值)格式文件，包含空间坐标信息和温度值。
优选地，所述步骤四、步骤五和步骤六都是联合NX软件和MSC Patran软件进行有限元模型温度场赋值，实现Nastran格式有限元的快速赋值。
优选地，所述步骤四采用坐标转换后的温度场数据对有限元模型进行温度加载，NX软件自动将坐标信息与对应位置的有限元节点进行匹配，使温度值数据加载到对应的有限元节点上。
优选地，所述步骤五保存并导出的文件为可进行MSC Nastran分析计算的bdf格式有限元模型文件，包含了温度场特征信息。
优选地，所述步骤六利用MSC Patran或Nastran软件导入温度场模型后可直接计算出卫星结构的热变形结果。
本发明的积极进步效果在于：本发明能够在进行卫星热变形试验后利用实测数值快速对MSC Patran/Nastran有限元模型进行温度场赋值，从而能在短时间内完成结构热变形分析，进而指导设计生产，缩短卫星研制周期。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合附图给出本发明较佳实施例，以详细说明本发明的技术方案。
如图1所示，本发明适用于卫星有限元热变形分析的温度场快速赋值方法包括以下步骤：
步骤一，建立需要进行热变形分析的卫星有限元模型，包括网格划分、多点约束定义、材料定义，导出模型文件；
步骤二，对温度场测试数据进行坐标转换，使测试数据坐标与有限元模型坐标一致；
步骤三，将模型文件导入NX软件，根据热变形分析方案对模型进行分组设置；
步骤四，对NX模型文件中的分组有限元分别进行温度场加载；
步骤五，将NX模型文件保存并导出有限元格式文件；
步骤六，在MSC Patran中导入温度场模型文件用于热变形分析可获得分析数据。
所述步骤一是利用MSC Patran有限元软件建立需要进行热变形分析的卫星有限元模型。
所述步骤一中导出模型文件为为bdf格式的Nastran模型文件，获得包含卫星结构所有热相关特征信息，特别是热膨胀系数等。
所述步骤二是利用MATLAB编写的坐标转换程序对卫星温度场实测数据进行坐标转换。
所述步骤二中的坐标转换的文件均为excel编制的csv格式文件，包含空间坐标信息和温度值。
所述步骤二经坐标转换后的温度场的坐标信息与有限元模型的坐标信息一致。
所述步骤四、步骤五和步骤六都是联合NX软件和MSC Patran软件进行有限元模型温度场赋值，实现Nastran格式有限元的快速赋值。
所述步骤四采用坐标转换后的温度场数据对有限元模型进行温度加载，NX软件自动将坐标信息与对应位置的有限元节点进行匹配，使温度值数据加载到对应的有限元节点上。
所述步骤五保存并导出的文件为可进行MSC Nastran分析计算的bdf格式有限元模型文件，包含了温度场特征信息。
所述步骤六利用MSC Patran或Nastran软件导入温度场模型后可直接计算出卫星结构的热变形结果，具体地，可进行静力分析，获得温度场作用下的应力应变数据。
以上所述的具体实施例，对本发明的解决的技术问题、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种适用于卫星有限元热变形分析的温度场快速赋值方法，其特征在于，其包括以下步骤：
步骤一，建立需要进行热变形分析的卫星有限元模型，包括网格划分、多点约束定义、材料定义，导出模型文件；
步骤二，对温度场测试数据进行坐标转换，使测试数据坐标与有限元模型坐标一致；
步骤三，将模型文件导入NX软件，根据热变形分析方案对模型进行分组设置；
步骤四，对NX模型文件中的分组有限元分别进行温度场加载；
步骤五，将NX模型文件保存并导出有限元格式文件；
步骤六，在MSC Patran中导入温度场模型文件用于热变形分析可获得分析数据。
2.如权利要求1所述的适用于卫星有限元热变形分析的温度场快速赋值方法，其特征在于，所述步骤一是利用MSC Patran有限元软件建立需要进行热变形分析的卫星有限元模型。
3.如权利要求2所述的适用于卫星有限元热变形分析的温度场快速赋值方法，其特征在于，所述步骤一中导出模型文件为为bdf格式的Nastran模型文件，获得包含卫星结构所有热相关特征信息。
4.如权利要求1所述的适用于卫星有限元热变形分析的温度场快速赋值方法，其特征在于，所述步骤二是利用MATLAB编写的坐标转换程序对卫星温度场实测数据进行坐标转换。
5.如权利要求4所述的适用于卫星有限元热变形分析的温度场快速赋值方法，其特征在于，所述步骤二中的坐标转换的文件均为excel编制的csv格式文件，包含空间坐标信息和温度值。
6.如权利要求1所述的适用于卫星有限元热变形分析的温度场快速赋值方法，其特征在于，所述步骤四、步骤五和步骤六都是联合NX软件和MSCPatran软件进行有限元模型温度场赋值，实现Nastran格式有限元的快速赋值。
7.如权利要求6所述的适用于卫星有限元热变形分析的温度场快速赋值方法，其特征在于，所述步骤四采用坐标转换后的温度场数据对有限元模型进行温度加载，NX软件自动将坐标信息与对应位置的有限元节点进行匹配，使温度值数据加载到对应的有限元节点上。
8.如权利要求6所述的适用于卫星有限元热变形分析的温度场快速赋值方法，其特征在于，所述步骤五保存并导出的文件为可进行MSC Nastran分析计算的bdf格式有限元模型文件，包含了温度场特征信息。
9.如权利要求6所述的适用于卫星有限元热变形分析的温度场快速赋值方法，其特征在于，所述步骤六利用MSC Patran或Nastran软件导入温度场模型后可直接计算出卫星结构的热变形结果。
CN201710419042.2A 2017-06-06 2017-06-06 适用于卫星有限元热变形分析的温度场快速赋值方法 Pending CN107391786A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710419042.2A CN107391786A (zh) 2017-06-06 2017-06-06 适用于卫星有限元热变形分析的温度场快速赋值方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710419042.2A CN107391786A (zh) 2017-06-06 2017-06-06 适用于卫星有限元热变形分析的温度场快速赋值方法

Publications (1)

Publication Number Publication Date
CN107391786A true CN107391786A (zh) 2017-11-24

Family

ID=60332130

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710419042.2A Pending CN107391786A (zh) 2017-06-06 2017-06-06 适用于卫星有限元热变形分析的温度场快速赋值方法

Country Status (1)

Country Link
CN (1) CN107391786A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113283018A (zh) * 2021-07-21 2021-08-20 星河动力(北京)空间科技有限公司 运载火箭的温度场信息的重构方法、装置及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101727527A (zh) * 2009-12-28 2010-06-09 中国农业大学 动力学分析软件与有限元分析软件间数据自动转换的方法
CN102682171A (zh) * 2012-05-15 2012-09-19 中国电子科技集团公司第五十四研究所 一种对高精度复合材料天线反射体单元模具热变形进行补偿的方法
CN103177153A (zh) * 2013-01-25 2013-06-26 北京空间飞行器总体设计部 一种有限元温度场结果的映射方法
CN103279586A (zh) * 2013-01-25 2013-09-04 北京空间飞行器总体设计部 一种有限元温度场的二次映射方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101727527A (zh) * 2009-12-28 2010-06-09 中国农业大学 动力学分析软件与有限元分析软件间数据自动转换的方法
CN102682171A (zh) * 2012-05-15 2012-09-19 中国电子科技集团公司第五十四研究所 一种对高精度复合材料天线反射体单元模具热变形进行补偿的方法
CN103177153A (zh) * 2013-01-25 2013-06-26 北京空间飞行器总体设计部 一种有限元温度场结果的映射方法
CN103279586A (zh) * 2013-01-25 2013-09-04 北京空间飞行器总体设计部 一种有限元温度场的二次映射方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113283018A (zh) * 2021-07-21 2021-08-20 星河动力(北京)空间科技有限公司 运载火箭的温度场信息的重构方法、装置及设备

Similar Documents

Publication Publication Date Title
CN103593518B (zh) 一种基于模态试验数据的飞行器模型修正系统
CN110516325A (zh) 一种cae自动化仿真分析方法和系统
US11281824B2 (en) Authoring loading and boundary conditions for simulation scenarios
CN106682301A (zh) 一种面向多点激励振动试验的系统模型修正方法
CN103207926A (zh) 使用显式有限元分析、采用子循环和质量缩放的组合技术数值模拟产品的结构性能
CN106897527A (zh) 一种车辆悬架台架耐久载荷分析方法及装置
CN104298823B (zh) 一种高、低温天平分析方法及分析系统
CN112414668B (zh) 一种风洞试验数据静气弹修正方法、装置、设备及介质
CN110032756A (zh) 基于流函数分数坐标系变换的流动边界层数值分析方法
JP2019507405A (ja) 統合システム用の共同シミュレーションを構成する方法
CN115438598A (zh) 基于一般时间根方尺度的雷诺应力湍流模型的数值方法
CN109711002A (zh) 新车型载荷分解方法及装置
CN112861374B (zh) 基于预控制器的多物理耦合仿真处理方法、装置和设备
CN107391786A (zh) 适用于卫星有限元热变形分析的温度场快速赋值方法
CN105404751A (zh) 基于热-力-电磁场网络统一的实现方法
CN107153755A (zh) 一种页岩气井数值模拟的求解方法
US10599879B2 (en) Optimal pressure-projection method for incompressible transient and steady-state navier-stokes equations
Mani et al. Geometry optimization in three-dimensional unsteady flow problems using the discrete adjoint
CN107545097A (zh) 一种高效铝合金榫卯连接结构的cae模拟优化分析方法
US8180598B2 (en) System and method for performing thermal analysis on a building through universal meshing
CN115293069B (zh) 一种用于飞行器外流场仿真控制参数智能优化的系统
CN115146402A (zh) 高温合金的氧化过程模拟方法及装置
Anand et al. Process automation of simulation using Toolkit/Tool Command Language (TK/TCL) scripting
Liou et al. A novel approach of three-dimensional hybrid grid methodology: Part 2. Flow solution
CN111199085A (zh) 抽象表达模型及其状态的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20171124