GB/T29081-2012
航天器模态计算方法
Modalanalysisforspacecraft
- 中国标准分类号（CCS）V70
- 国际标准分类号（ICS）49.020
- 实施日期2013-07-01
- 文件格式PDF
- 文本页数14页
- 文件大小386.97KB
以图片形式预览航天器模态计算方法
航天器模态计算方法
国家标准 GB/29081一2012 航天器模态计算方法 Modalanalysisforspaceraft 2012-12-31发布 2013-07-01实施 国家质量监督检监检疫总局 发布 国家标准花管理委员会国家标准
GB/T29081一2012 前 言 本标准按照GB/T1.1一2009给出的规则起草
本标准由航天科技集团公司提出
本标准由全国宇航技术及其应用标准化技术委员会(SAC/TC425)归口
本标雅起草单位;空间技术研究院总体部 本标准主要起草人:柴洪友、钱志英
GB/T29081一2012 引 言 本标准属于航天国家标准体系
航天国家标准体系适用于航天领域国家标准的制修订和 管理,覆盖航天管理、航天技术,航天应用与服务三大领域,是指导航天器和运载火箭项目管理、工程研 制、航天发射服务、卫星在轨应用等活动的依据
随着我国航天事业的迅猛发展,长寿命、高可靠性、功能复杂的航天器不断涌现
在航天器结构设 计和验证过程中,对航天器模态的要求(固有频率的要求)是大多数航天器结构设计的重要要求
航天 器结构模态分析是航天器结构分析的重要内容之一,在航天器结构研制过程中具有重要的作用
本标准为获取合理.准确、可靠的航天器模态计算结果提供依据和指导,对航天器模态分析的模型 一的规范,目的是保证工程研制过程中航天器模态分析 建立、模型检查以及模态特性分析等方面进行统一 的准确性
航天器模态计算在航天器研制过程中具有重要作用,如对避免频率羁合引起的共振,航天器 构型方案选择、结构故障诊断、仪器设备的动力学环境预示等都有着重要的影响
当前国内外航天事业 迅猛发展,各类型号的航天器都在紧密的研制开发过程中,因而有必要将航天器模态计算的过程进行统 -规定,确保计算,分析过程的规范化和结果的有效性
GB/T29081一2012 航天器模态计算方法 范围 本标准规定了航天器的模态计算方法
本标准适用于航天器的模态计算和航天器界面广义刚度矩阵和广义质量矩阵的计算
舱段、组件 的模态计算可参照使用
术语和定义 下列术语和定义适用于本文件
静定界面statieallydetermimateinterfaee 在结构运动时可视为刚体的界面
注:结构通过界面与其他结构相连,对于静定界面只需通过某个参考点的自由度就能表示出所有界面自由度的运动 2 v 结构理想化strueturalidealizatiom 把实际结构转化成一种能够表现结构的几何特征和连接关系,并能有效反映结构的力学特性和进 行力学分析的形式的过程
注;如将具有一定几何特征的构件抽象为点、线、面或体,进而简化为结构力学中的质点、弹簧、杆、梁、板、壳、体等
航天器模态计算步骤 航天器模态计算一般步骤如图1所示 模态计算相关 设计参数 试脸数据 参数获取 结构理想化 有限元建树 限元模型检 查 否通过 模态计算 图1航天器模态计算步骤
GB/T29081一2012 航天器模态计算相关参数获取 航天器模态计算结果主要取决于以下物理参数 a)航天器系统的质量和惯量分布 b航天器结构的刚度特性; 航天器结构的边界条件
其中,刚度特性和边界条件的确定除考虑设计参数外,还应参考试验结果不断修正
航天器有限元模型的建模 5.1结构理想化 5.1.1概述 由于航天器系统的质量分布相对于结构刚度更容易准确地反映在模型中,因此在建模时,应当充分 关注结构刚度的模拟
对于影响结构刚度较大的因索,不应过分简化 5.1.2有限单元类型的确定 单元类型根据结构的几何特征,受力情况及分析软件确定,主要应用的单元类型包括质点单元、弹 黄元.杆单元、梁单元.板壳单元.,体单元等
在选择单元时.以下内容将有助于建立合理的计算模型 a 衔架结构采用杆单元,构架结构采用梁单元,板和壳结构采用板壳元 b)当连接结构刚度很大而无必要细化时,采用刚性单元模拟; c)接头和特殊连接处用弹簧元或多点约束模拟; D 包带连接解锁机构用弹簧元模拟; 胶接,铆接和焊接一般理想化为刚性连接 连接形式要视具体情况而定,有时理想化为刚性连接(限制6个自由度),有时理想化为限制 5 个自由度的连接或限制3个平移自由度的钞接 较大的设备用质点单元模拟,并用多点约束或刚性单元连接到结构板的安装位置上,有较多固 定点的设备作为均布的非结构质量定义于适当区域的结构板单元上; h 某些仪器参与结构承力时应适当模拟仪器结构刚度; iD 液体质量可用等效质量模拟; j 单元之间的位移应当协调,包括平动和转动 k)对于不易进行简化或者进行了简化但难以确定参数的局部结构,可进行局部试验,根据试验数 据进行建模
5.1.3 位移边界条件的确定 根据具体情况采用固支边界条件,简支边界条件或弹性边界条件
为了解简化程度对分析结果的影响,可以采用不同类型的边界条件进行分析,比较后最终确定合适 的边界条件 5.2坐标系选择 坐标系选择应遵从以下原则: 航天器有限元模型的整体坐标系一般应采用航天器总体定义的机械总体坐标系; a
GB/T29081一2012 b根据航天器具体的结构形式建立局部坐标系以方便计算
单位制选择 采用千克、米、秒作为基本单位
5 单元网格划分 按照结构理想化选定的单元,将结构离散为有限个单元的集合体
单元网格划分应考虑以下因素: 单元划分的疏密应能够有效描述所关心频率对应的结构振型 a b)单元的几何形状应满足一定的要求,在划分三角形单元时,三条边应尽量相等，一般不应出现 钝角三角形;在划分四边形单元时,应尽量接近于正方形,通常情况下,四边形的内角不小于 45°,不大于135"; 对于几儿何特性(如尺寸和截面等)或物理特性(如弹性系数等)有突变之处,在划分单元时应作 为单元的边线,必要时还应对这些部位的单元进行细化处理 单元形状比例应符合软件使用说明书的要求
5 单元特性确定 5.5.1 物理特性确定 定义材料的密度、弹性性能(或刚度性能),非结构质量等物理特性参数
5.5.2几何特性确定 根据所选单元定义截面几何特性参数
在应用杆单元时,需要定义截面面积(对于抗扭杆还应定义抗扭截面系数 在应用梁单元时,需要定义梁截面的参考坐标系和截面形状,大小,也可直接输人截面面积、抗弯和 抗扭截面系数、剪切修正系数等
在应用板单元时,普通金属薄板只需定义厚度;对于夹层板或层合复合材料板可按铺层进行定义， 也可以利用试验或分析结果直接定义整个夹层板或层合复合材料板的刚度系数
当板单元或梁单元的节点与自身中性面不重合时,需要定义中性面的偏置
当采用梁单元来模拟加强筋时,需要定义加强筋中性面的偏置位置
航天器有限元模型检查 6.1几何检查 利用建模软件的前处理功能检查有限元模型 几何形状检查 b)单元雅克比acobi)行列式检查; c自由边检查和单元重叠检查
质量特性检查 计算有限元模型的质量质心位置和转动惯量,检查模型的质量特性(主要是质量、质心)是否符合 设计值或测试值
一般情况下,模型的质量误差应小于0.3%
当航天器结构本体尺寸不大于1m时， 质心位置误差应小于3mm;当航天器结构本体尺寸大于1m时,质心位置误差应小于0.3%
为了掌握更准确的质厨分布,通常要对模型各结构部件、大型天线、太阳翼相机、舱段等分别进行
GB/T29081一2012 质量特性检查
3 协调性检查 o 下列情况可能造成计算中断,需调整节点单元或约束条件予以解决: a)单元间位移不协调; b) 单元比例失调; e)不符合单元形成规则 d)结构刚度矩阵病态或奇异; 节点约束不符合所选用的单元类型等
a 刚体模态检查 在无位移约束边界条件下(去掉所有位移边界条件),计算得到前六阶模态应满足固有频率均小于 静定结构最小固有频率的0.1%,并且六个模态的有效质量比之和在六个方向上均大于99.9%
否则 应查找原因,修改模型后再调试,直到满足要求
6.5力平衡检查 给定边界约束条件,分别在三个方向上施加1片的加速度载荷,确认支反力仅产生在被约束的节点 上,每方向的支反力相加后应该等于模型的总重量,否则检查结构间是否应用了不适当的弹簧元或刚性 元等 航天器模态计算内容 7.1模态计算 7.1.1固有频率及振型 7.1.1.1计算公式 以节点位移为未知变量,将结构离散成具有n个自由度的系统,则结构的自由振动方程见公式(1). 结构的特征方程见公式(2)
模态计算就是通过结构建模计算得到刚度矩阵K和质量矩阵M,然后求 解特征方程(2)得到和与之对应的,并分析必的特征
Mr十Kr 0 式中 M 质量矩阵,nXn阶; 加速度向量,n×1阶 . K 刚度矩阵,n×n阶 -位移向量,n×1阶
K一oM)e=0 式中 o 特征值,工程上称为固有圆频率; -特征向量,n×1阶,工程上称为振型 9 解特征方程(2)的广义特征值,得到第阶固有圆频率及其对应的振型9.(k=1,2,,m;mn)
为正则化模态(质量归一) 第走阶固有额率按公式(3)计算
GB/T29081一2012 3 f= 2x 式中: 第k阶固有频率; 第阶固有圆频率
o 7.1.1.2计算要求 计算航天器前m个模态的固有频率和振型,m一般为100Hz之内的所有模态阶数
如果每个方 向上模态有效质量的总和都大于80%,可认为主要模态均已算出,否则,需追加模态计算阶数
对于装有大型部件(如大型天线、太阳翼、相机、贮箱等)或由多舱段组成的航天器,可以将航天器本 体,大型部件,舱段等分别作为子结构,应用模态综合法计算模态
推荐采用固定界面模态综合法进行 模态计算
固定界面模态综合法的相关内容见附录A
7.1.2有效模态质量比 7.1.2.1计算公式 第人阶模态的有效质量按公式(4)计算
M,=4从 式中: M -第人阶模态的有效质量矩阵,6×6阶 其中4.(1×6阶)按公式(5)计算 (5 4 =T'Me 式中: T 结构刚体运动变换矩阵,n×6阶
-第阶模态向量,n×1阶
P 结构内部自由度刚体质量矩阵按公式(6)计算
M,
=下IMr,-习Me 式中: -结构内部自由度刚体总质量矩阵.6X6阶 M
第人阶模态第1个自由度的有效质量比劝按公式(7)计算
=1,2,,6 功=m;/m(k=1l,2,,m ;l 式中: n -第人阶模态在第1个自由度上的模态有效质量相对刚体质量贡献的百分比 第阶模态有效质量矩阵第！个自由度的有效质量,是有效质量矩阵的第！个对角线元素 m 刚体总质量矩阵第1个自由度的质量,是刚体总质量矩阵中第！个对角线元索
n 7.1.2.2计算要求 当需要定性判断一个模态的运动特征时,应开展模态的有效质量比的计算
如果某阶模态在某方 向上的有效模态质量比大于10%,则可定性认为此阶模态是主要模态
模态的有效质量计算还可以作为模态截断的一种判定方法
当所得到的前阶模态的模态有效 质量之和占总质量的比例很大时(在航天器模态计算中一般为80%),就可以认为主要模态已包含在前 m阶模态中
GB/T29081一2012 7.1.3模态动能比 7.1.3.1计算公式 第k个正则化模态9的总动能T按公式(8)计算
几-喜M,,-叫 该模态中第1个自由度的位移为少A,其对应的广义质量为m,则在该模态中,该自由度的模态动能 a按公式(9)计算
=一m;Y 在计算时,可以把节点分别归并于若干个部件(或部位),对各部件进行模态动能计算,得到部件的 模态动能
部件的模态动能就是该部件所包含的节点模态动能的和 在实际计算中,通常按公式(10)计算模态动能比6 (10 -一-n学 7.1.3.2计算要求 当需要判定结构的哪些节点在模态运动中起主导作用时,需要开展模态动能比的计算
7.1.4模态势能比 7.1.4.1计算公式 第k个正则化模态9的总势能V按公式(11)计算
V .(11l -Kp;-" 若某个单元在总体坐标系中的单元刚度矩阵为k",该模态中该单元各节点的位移列阵为外,则在 该模态中,该单元的模态势能砍按公式(12)计算
(12 欢=,'k好 在计算时,可以把单元分别归并为若干个部件(或部位).对各部件进行模态势能计算
部件的模态 势能是该部件所包含单元模态势能的和
在实际计算中,通常按公式(13)计算模态势能比
U -丝 (13 =岗 o 7.1.4.2计算要求 当需要判定在某阶模态下结构的哪些部位在模态运动中变形更严重时,应开展模态势能比的计算
7.2广义刚度矩阵和广义质量矩阵计算 7.2.1计算公式 对无阻尼系统,结构的振动方程按照内部自由度r,(n个自由度)和界面自由度r,(b个自由度)可 写为公式(14)
mnm ]]-谈 (14 nl
GB/T29081一2012 式中: 内部自由度对应的质量矩阵 m 内部自由度与边界自由度合的质量矩阵,其中m,是m,的转置; m;,m 边界自由度对应的质量矩阵 m 内部自由度对应的刚度矩阵 k， 内部自由度与边界自由度稠合的刚度矩阵,其中A;是友,的转置 k,k 边界自由度对应的刚度矩阵
e f 界面力
航天器界面一般视为静定界面,由静定界面的定义,所有界面自由度r,可以用一个参考点的自由 度r,来表示,见公式(15)
r=T,r (15 式中 界面自由度的刚体位移变换矩阵,"=1,2,6,b×6阶
公式(14)可写成公式(16)
m1 Imn 16 ]] Lmm 其中m n、从、k.、,友，和/,分别按公式(17)~公式(21)计算
17 nm
=T,m,;T 一m, =mnT 18 71 A，=TT,k,T 19 20 k
=k,=k,T f，=TT 21 式中 各分量的下脚标,r=1,2,,6
广义出口刚度矩阵关和广义出口质量矩阵"分别按公式(22)和公式(23)计算 [2? 22 部 23 5 " 4M
式中: -各分量的下脚标，=1,2,m 其中Q'、)，和,分别按公式(24)公式(26)计算
Q=diag(of“ w” 25 M
=p'm,T,十 9 n
26 T,m,T十T,m，十m.T Mm=m十 式中: 模态向量矩阵,n×阶; Pa 内部自由度相对于界面参考点的刚体位移变换矩阵,n×6阶 7.2.2计算要求 在计算航天器界面的广义刚度矩阵和广义质量矩阵时,应选择航天器与运载火箭连接面的中心作
GB/T29081一2012 为界面参考点
航天器的模态自由度至少取100H2之内的所有模态阶数
在提供航天器界面的广义刚度矩阵和广义质量矩阵用于航天器与运载火箭耦合载荷分析时,通常 需要提供航天器刚体模态检查结果
GB/T29081一2012 附 录A 规范性附录 固定界面模态综合法 模态综合法将整体结构划分为若干个子结构,子结构之间的相互连接处形成界面
固定界面模态 综合法(Craig-IBampton法)采用固定界面的正则化振动模态(所有界面自由度固定时的正则振动模态 和对于界面坐标的约束模态(界面上某一个自由度单位位移,其他界面自由度固定的静平衡模式)构成 极态失量矩阵来描述结构变形 第一步,将整体结构系统)分割成若干子结构,使子结构之间的连接界面完全固定,进行子结构模 态计算及第一次自由度变换 对无阻尼系统,系统自由振动时,非界面力为零,子结构的振动方程按照内部自由度r，和界面自由 度r,可写为公式(A.1)
Imnim l]" (A.1 [m
m 对于固定界面,即界面自由度为零(工=0),于是得到公式(A.2) m;.r，十人r，=0 A.2 ]"为截断模态个数
对于截断模态个 由方程(A.2)可以得到正则化模态内=[内声 -般要求子结构截断模态的最高频率大于整体最高分析频率 数， 与公式(A.1)相对应的子结构静力平衡方程见公式(A.3)
人" A.3 IC- Lkk 其中r,可表示为公式(A.4)
.r，=中j.x 其中中,可表示为公式(A.5). -k
k， A.5 做自由度变换,把子结构自由度r变换到模态自由度q上,见公式(A.6) p 1g = A.6 lxr 式中 -对应振动模态的模态自由度; g -对应约束模态的模态自由度.g=r,,即约束模态自由度就是界面自由度; -从物理坐标到模态坐标的变换矩阵
按公式(A.6)将子结构运动方程(A.1)变换到模态自由度q上,见公式(A.7)
丽十9十局=了7 式中: 丽 与模态坐标对应的质量矩阵向量， -与模态坐标对应的刚度矩阵向量; -与模态坐标对应的力矩阵向量
其中丽、瓦和了可按公式(A.8)一公式(A.16)计算
GB/T29081一2012 =重Tmnp A.8 m1一 =p'ke A.9 =Tf A.10 m A.11 m一 m [万 A.12) A.13 m,=m=中,m,中,十中,m A.14 m,=m,中,m.中，十pm;十m,中 飞
=diag(oo A.15 A.16 人;=k》十k,? 第二步,建立系统方程
不失一般性,考虑有两个子结构的系统
建立两个子结构A和B在模态 自由度下的振动方程,见公式(A.17)
7m1A 9gA i-"
A.17) m lo 式中: -与A子结构相关的质量矩阵 mA m 与B子结构相关的质量矩阵; 与A子结构相关的模态自由度 q -与B子结构相关的模态自由度; 人
与A子结构相关的刚度矩阵; 万 与B子结构相关的刚度矩阵; 了 与A子结构相关的力向量; 了 -与B子结构相关的力向量
在界面上,由界面位移协调条件4 选择系统的广义自由度为Q=[g ]' gB qA=qB=qj
代人公式A.17),得到系统无阻尼自由振动方程,见公式(A.18) nmQ十AQ=0 A.18 其中相关系数m、和可分别按公式(A.19)和公式(A.20)计算
mgA (A.19 n I mB LmnAmnBmA十m)B o o (A.20 k= o 万;八十飞m 式中: 与A子结构相关的质量矩阵向量分块,=1,2,,m; mn八 与B子结构相关的质量矩阵向量分块，=1.2, ,n; mB 10
GB/T29081一2012 分A--与A子结构相关的刚度矩阵向量分块，=1,2,,m; 飞8与B子结构相关的刚度矩阵向量分块，！=1,2,,m
第三步,解公式(A.18)的特征值,可得到模态自由度下的振型,再按公式(A.6)计算结构系统自由 度下的振型
航天器模态计算方法GB/T29081-2012
航天器的模态计算是指通过对其结构进行分析和计算，得出其不同振动模式及其频率、振型等参数，用于评估航天器在运行过程中的动态响应和稳定性。
GB/T29081-2012标准是我国航天器模态计算方面的规范，主要涉及到航天器结构动力学分析的基本原理、模态计算方法、模态试验方法等方面。
在进行航天器模态计算时，需要对其结构进行建模，并采用相应的数值计算方法进行求解。常见的数值计算方法包括有限元法、边界元法、有限差分法等。
除了数值计算方法外，还需要注意到一些影响航天器模态计算精度的因素，例如材料的非线性特性、结构的非均匀性等。针对这些因素，需要进行相应的修正和补偿。
在航天器结构的设计与分析中，模态计算是一项重要的工作。通过合理的模态计算方法和技术手段，可以有效提高航天器的设计质量和运行稳定性，从而保证航天任务的顺利实施。
