CN112231861B - 一种抑制调姿共振的集群式控制力矩陀螺隔振方法 - Google Patents
一种抑制调姿共振的集群式控制力矩陀螺隔振方法 Download PDFInfo
- Publication number
- CN112231861B CN112231861B CN202011111653.9A CN202011111653A CN112231861B CN 112231861 B CN112231861 B CN 112231861B CN 202011111653 A CN202011111653 A CN 202011111653A CN 112231861 B CN112231861 B CN 112231861B
- Authority
- CN
- China
- Prior art keywords
- model
- spring
- control moment
- damper
- vibration
- 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.)
- Expired - Fee Related
Links
- 238000002955 isolation Methods 0.000 title claims abstract description 43
- 230000002401 inhibitory effect Effects 0.000 title claims abstract description 10
- 238000006073 displacement reaction Methods 0.000 claims abstract description 40
- 238000013016 damping Methods 0.000 claims description 20
- 238000000034 method Methods 0.000 claims description 20
- 239000007787 solid Substances 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 13
- 238000005381 potential energy Methods 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 8
- 239000000463 material Substances 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 8
- 150000001875 compounds Chemical class 0.000 claims description 3
- 230000000149 penetrating effect Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 description 18
- 238000004458 analytical method Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 8
- 230000004044 response Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- 238000009434 installation Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000011967 cystometrography Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000003380 propellant Substances 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
- F16F15/00—Suppression of vibrations in systems; Means or arrangements for avoiding or reducing out-of-balance forces, e.g. due to motion
- F16F15/02—Suppression of vibrations of non-rotating, e.g. reciprocating systems; Suppression of vibrations of rotating systems by use of members not moving with the rotating systems
- F16F15/04—Suppression of vibrations of non-rotating, e.g. reciprocating systems; Suppression of vibrations of rotating systems by use of members not moving with the rotating systems using elastic means
-
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mechanical Engineering (AREA)
- Computational Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Acoustics & Sound (AREA)
- Vibration Prevention Devices (AREA)
Abstract
一种抑制调姿共振的集群式控制力矩陀螺隔振方法,涉及隔振领域。本发明是为了解决目前控制力矩陀螺上控制电机产生的振动位移与装有隔振器的控制力矩陀螺某一阶振型会发生耦合,导致控制力矩陀螺共振的问题。构建集群式控制力矩陀螺实体模型和动力学模型;根据该两个模型得到固有频率,根据固有频率得到隔振器下支点空间坐标理论值;纠正隔振器下支点空间坐标理论值,得到最终的隔振器下支点空间坐标;根据最终的隔振器下支点空间坐标得到多个弹簧和阻尼器的结构参数,并设计弹簧模型和阻尼器模型;按照弹簧模型、阻尼器模型和最终阻尼器下支点空间坐标调整隔振器结构和4个隔振器下支点与固定平面的连接位置。它用于确定隔振器安装位置。
Description
技术领域
本发明涉及星载控制力矩陀螺的隔振设计,特别是一种抑制调姿共振的集群式控制力矩陀螺隔振设计。属于隔振领域。
背景技术
控制力矩陀螺(Control Moment Gyroscopic,简称CMG)作为卫星姿态控制中的重要机构,由于其不消耗推进剂、调控精度高、相应速度快等特点,使其广泛服务于大型高性能卫星平台。然而,CMG主要是依靠其内部动量轮的旋转以及控制电机驱动的往复运动部分来完成规定任务的,动量轮和控制电机都会产生较大程度的振动,这对卫星的调姿精度会有不利影响。针对这一现象,通过设计合适的隔振器,使得隔振器的起始隔振频率尽可能低于CMG的振动频率,以实现良好的隔振。与目前普遍采用柔性铰的隔振器设计不同,考虑到CMG在上升段的过载以及柔性铰设计/加工过程中便利性的不足,可采用固接的方式将CMG与隔振器相连,既能保证隔振效果,又能满足隔振器对过载条件的需求。
然而,通过目前对单一CMG隔振器的设计和试验发现,CMG控制电机产生的振动位移与装有隔振器的CMG某一阶振型会发生耦合,导致CMG的共振。
发明内容
本发明是为了解决目前控制力矩陀螺上控制电机产生的振动位移与装有隔振器的控制力矩陀螺某一阶振型会发生耦合,导致控制力矩陀螺共振的问题。现提供一种抑制调姿共振的集群式控制力矩陀螺隔振方法。
一种抑制调姿共振的集群式控制力矩陀螺隔振方法,所述方法包括以下步骤:
步骤1、支架安装在固定平面上,将多个控制力矩陀螺安装在支架上,根据多个控制力矩陀螺与支架的连接位置,构建集群式控制力矩陀螺实体模型,将4个隔振器安装在支架和固定平面之间,根据隔振器上支点和下支点分别与支架和固定平面的连接位置,构建动力学模型,所述每个隔振器由一个弹簧和沿弹簧伸缩方向穿在弹簧内部的阻尼器组成;
步骤2、根据集群式控制力矩陀螺实体模型和动力学模型得到固有频率,根据固有频率得到隔振器下支点空间坐标理论值;
步骤3、根据集群式控制力矩陀螺实体模型和动力学模型建立含有集群式控制力矩陀螺和隔振器的有限元模型,将步骤2中得到的隔振器下支点空间坐标理论值添加到所述有限元模型中,使用有限元软件对隔振器下支点空间坐标理论值进行解析,得到隔振器的固有频率,判断隔振器的固有频率是否处于预设固有频率范围,如果是,则根据固有频率得到最终的隔振器下支点空间坐标,如果否,则调整隔振器下支点空间坐标,直到计算的固定频率达到预设固有频率范围,根据最终生成的固有频率得到最终的隔振器下支点空间坐标;
根据最终的隔振器下支点空间坐标得到多个弹簧和阻尼器的结构参数,根据每个所述结构参数设计弹簧模型和阻尼器模型;
步骤4、从得到的多个弹簧模型和阻尼器模型任选出一个弹簧模型和阻尼器模型,按照选出的弹簧模型、阻尼器模型和步骤3中得到的最终的阻尼器下支点空间坐标调整步骤1中隔振器的结构和4个隔振器下支点与固定平面的连接位置,建立隔振器与集群式控制力矩陀螺的最终连接,实现对集群式控制力矩陀螺的隔振。
优选地,步骤1中,动力学模型包括集群式控制力矩陀螺的动能和隔振器的弹性势能,
集群式控制力矩陀螺的动能E表示为;
式中,Ix、Iy和Iz表示集群式控制力矩陀螺所处的平台质心且平行于坐标系轴的转动惯量,所述平台包括集群式控制力矩陀螺和支架,m表示平台的质量,和是集群式控制力矩陀螺质心的三个平动位移xc、yc和zc的导数,是平台质心位置转动位移的导数,
其中,集群式控制力矩陀螺质心的三个平动位移xc、yc和zc表示为:
式中,x、y和z是集群式控制力矩陀螺在全局坐标系原点位置的三个平动位移,x0、y0和z0是集群式控制力矩陀螺质心在全局坐标系下的坐标,α、β和γ是平台质心位置转动位移。
优选地,单个隔振器的弹性势能Epi为:
式中,v为设定的隔振器上弹簧与集群式控制力矩陀螺连接点的位移向量,e为弹簧的轴向向量,vge为弹簧在轴向的位移,v-(vge)e为弹簧在轴向的位移与轴向垂直方向的位移,ka和kr分别为设定的弹簧轴向和径向刚度,e=ei,i=1,2,3,4,
四个隔振器上支点弹簧的轴向向量ei,i=1,2,3,4,分别表示为:
式中,θ1、θ2分别表示隔振器与全局坐标系x轴和z轴方向的夹角,
v由四个上支点的位移分量xi、yi、zi,i=1,2,3,4表示,xi、yi、zi,i=1,2,3,4,表示为:
式中,a、b分别表示隔振器上支点距全局坐标系x轴和y轴的垂线距离;
假设四个弹簧的刚度相同,则总的弹性势能Ep为:
式中,vi和ei分别表示第i个位移向量和第i个隔振器上支点弹簧的轴向向量;
将公式4和公式5代入公式6,得到:
优选地,步骤3中,固有频率表示为:
|K-ω2M|=0 公式8,
式中,ω是减振系统的固有频率,K是刚度矩阵,M是质量矩阵,
优选地,步骤3中,建立含有集群式控制力矩陀螺的有限元模型的具体过程为:
根据集群式控制力矩陀螺模型中支架空间梁的分布情况,在有限元模型中的有限元软件中选择梁单元并赋予梁单元的截面属性,构建支架模型,根据集群式控制力矩陀螺模型中多个控制力矩陀螺的运动特征,在有限元软件中分别用点单元描述运动特征中多个控制力矩陀螺与支架模型连接的固定部分和多个控制力矩陀螺转动部分,并定义固定部分和转动部分的连接属性,完成含有集群式控制力矩陀螺的有限元模型的建立。
优选地,步骤3中,根据每个所述结构参数设计一个弹簧模型和阻尼器模型的具体过程为:
将每个弹簧和阻尼器结构参数输入到CAD软件中,由CAD软件建立3D模型,将建立的3D模型再导入到有限元软件中计算此时弹簧和阻尼器的刚度,并与各自预设的刚度进行比较,如果弹簧和阻尼器的刚度与各自预设的刚度均相同,则将建立的3D模型作为弹簧模型和阻尼器模型,如果弹簧和阻尼器的刚度与各自预设的刚度有一个不相同,则修改CAD软件中的3D模型,并再次使用有限元软件计算修正后的3D模型中的弹簧和阻尼器的刚度,重复在CAD中修改3D模型和使用有限元软件计算修正后的3D模型中的弹簧和阻尼器的刚度的过程,直到修正后的最终3D模型中的弹簧和阻尼器的刚度与各自预设的刚度均相同,则将修正后的最终3D模型作为弹簧模型和阻尼器模型。
优选地,所述弹簧结构参数包括弹簧长度、槽数和厚度,阻尼器结构参数包括阻尼系数、长度和阻尼材料耗散系数。
优选地,所述步骤3还包括:
对得到的弹簧模型和阻尼器模型进行验证。
优选地,对步骤3得到的弹簧模型和阻尼器模型进行验证,具体过程为:
有限元软件从步骤3中得到的多个弹簧模型和多个阻尼器模型中任意挑选一个弹簧模型和一个阻尼器模型与支架建立组合模型,直到完成所有弹簧模型和阻尼器与支架的组合情况,得到多个组合模型,每个组合模型输出一个固体频率,将固有频率满足预设固有频率范围的弹簧模型和阻尼器模型保留,从保留的弹簧模型和阻尼器模型中选择一个弹簧模型和阻尼器模型。
本发明的有益效果为:
本申请采用集群式控制力矩陀螺的方案,即多个控制力矩陀螺连接在支架上,通过在支架上合适的位置安装隔振器,就可以避免控制电机产生的振动位移与隔振系统振型之间的耦合,保证隔振性能的实现。
基于此思想,本申请先是建立由集群式控制力矩陀螺实体模型和动力学模型组成的实体结构,根据该两个模型得到固有频率,根据固有频率得到隔振器下支点空间坐标理论值,即得到隔振器下支点空间坐标理论值就知道隔振器的长度;再采用有限元模型验证该理论值是否正确,如果不正确纠正该理论值,最终得到满足固定频率的隔振器下支点空间坐标,根据满足固定频率的隔振器下支点空间坐标得到多个弹簧和阻尼器结构参数(几何参数:弹簧长度、弹簧槽数、弹簧厚度、阻尼器阻尼系数、阻尼器长度等),通过CAD软件和有限元软件结合判断每个弹簧和阻尼器结构参数组成的模型安装在支架上的刚度,最终选出符合刚度要求的弹簧和阻尼器模型,此时选出的参数可能有多个模型满足要求,后期工作人员根据喜好或者手边的材料选择适合自己的能够制作出来的一个弹簧和阻尼器模型,按照此弹簧和阻尼器模型的结构参数进行制作隔振器实体结构,将做好的隔振器实体结构安装在支架上。此时的隔振器位置不会与集群式控制力矩陀螺中的控制电机发生耦合,从而保证隔振性能。
有限元软件的作用是在软件中验证弹簧和阻尼器模型的刚度,避免在将实体隔振器与支架连接后,隔振器发生变形等刚度不够的问题。经过有限元软件验证后得到的弹簧模型和阻尼器模型,可以按照此模型建立实体隔振器结构。
结构设计人员可以基于此方法进行星载控制力矩陀螺的隔振设计,即在合适位置安装隔振器,以便对卫星产品进行基于振动控制的优化,极大地提高了工作效率。
附图说明
图1为一种抑制调姿共振的集群式控制力矩陀螺隔振方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是发明一部分实时例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
下面结合附图和具体实施例对本发明作进一步说明。
实施例1:
本发明提供的一种抑制调姿共振的集群式控制力矩陀螺隔振方法,如图1所示,其特征在于,所述方法包括以下步骤:
步骤1、支架安装在固定平面上,将多个控制力矩陀螺安装在支架上,根据多个控制力矩陀螺与支架的连接位置,构建集群式控制力矩陀螺实体模型,将4个隔振器安装在支架和固定平面之间,根据隔振器上支点和下支点分别与支架和固定平面的连接位置,构建动力学模型,所述每个隔振器由一个弹簧和沿弹簧伸缩方向穿在弹簧内部的阻尼器组成;
步骤2、根据集群式控制力矩陀螺实体模型和动力学模型得到固有频率,根据固有频率得到隔振器下支点空间坐标理论值;
步骤3、根据集群式控制力矩陀螺实体模型和动力学模型建立含有集群式控制力矩陀螺和隔振器的有限元模型,将步骤2中得到的隔振器下支点空间坐标理论值添加到所述有限元模型中,使用有限元软件对隔振器下支点空间坐标理论值进行解析,得到隔振器的固有频率,判断隔振器的固有频率是否处于预设固有频率范围,如果是,则根据固有频率得到最终的隔振器下支点空间坐标,如果否,则调整隔振器下支点空间坐标,直到计算的固定频率达到预设固有频率范围,根据最终生成的固有频率得到最终的隔振器下支点空间坐标;
根据最终的隔振器下支点空间坐标得到多个弹簧和阻尼器的结构参数,根据每个所述结构参数设计弹簧模型和阻尼器模型;
步骤4、从得到的多个弹簧模型和阻尼器模型任选出一个弹簧模型和阻尼器模型,按照选出的弹簧模型、阻尼器模型和步骤3中得到的最终的阻尼器下支点空间坐标调整步骤1中隔振器的结构和4个隔振器下支点与固定平面的连接位置,建立隔振器与集群式控制力矩陀螺的最终连接,实现对集群式控制力矩陀螺的隔振。
具体地,为了避免控制力矩陀螺中控制电机产生的振动位移,选择将一定数量控制力矩陀螺安装在支架上,形成集群式控制力矩陀螺。而后,根据隔振器固接在集群式控制力矩陀螺上的安装方式,建立其进行动力学模型。
对于步骤3,分别完成集群式CMG动力学模型和有限元模型的建立后,根据起始隔振频率(第一阶固有频率)分布区间的约束条件,优化其动力学模型中的设计参数(如隔振器空间坐标、刚度、阻尼等),而后在有限元模型的对应位置建立隔振器的单元并赋予优化后的参数特性,进一步完善集群式CMG的有限元模型。该步骤将集群式CMG的隔振系统动力学参数与有限元模型中的参数连接起来,建立了有限元模型与理论模型间的传递关系,借助于理论模型的初步优化,能大大降低使用有限元模型优化的次数,提升计算效率。
通过对集群式控制力矩陀螺起始隔振频率(固有频率)的需求,优化其动力学模型中的参数,确定隔振器的空间坐标和参数,进而指导阻尼器参数的设计。而后建立有限元模型与动力学模型之间的数据传递关系,补齐有限元模型中有关隔振器和阻尼器的空间坐标和参数特征。
在完成所有设计后,将集群式CMG的实体模型以及隔振器的实体模型导入到有限元软件中,定义它们的材料属性、接触关系、位移约束等,进行整体性仿真。首先根据模态分析的结果验证集群式CMG的固有频率是否满足需求,接着再将激励施加再CMG的相应位置,进行瞬态响应分析,以验证隔振前后输出的力和位移的改善程度。该步骤通过模态分析和瞬态响应分析,分别从频域和时域两种途径验证了本方案设计的可行性,在实际工程中时域仿真的价值在于能很好对CMG工作时振动环境进行预示,而通过瞬态响应分析,能为隔振系统的定型提供合理有效的依据。
本申请还可以输出有限元软件的所有参数,具体为:设计完成后,由有限元软件中输出模型文件,通过建立星载CMG产品设计参数与动力学模型之间的接口,解析并读取隔振器的设计参数,并按要求输出固有频率、隔振效果等计算结果及图形文件。
本申请提供的方法实现了动力学建模和有限元分析用于星载CMG的隔振设计,显著降低有限元分析的计算量,大大地提高了技术人员的工作效率。
步骤3中有限元软件采用网格划分及求解设置的方法更新隔振器下支点空间坐标,根据步骤3隔振器弹簧和阻尼器的空间坐标,更新有限元模型中相应节点的坐标,并重新生成集群式CMG有限元模型文件;此后,对坐标更新后的有限元模型重新划分网格,并把原始有限元模型的参数传递到新的有限元模型,其中与隔振器弹簧和阻尼器相关的参数需要用更新后的参数输入;最后,指定有限元软件的求解类型,并结合所使用计算机处理器的配置信息,驱动有限元软件求解器中相应模块进行计算,计算得到带隔振系统的集群式CMG的固有频率。至此,可以得到隔振器/阻尼器在更新参数后集群式CMG的振动特性,若该振动特性满足工程需求,则导出目前状态隔振器/阻尼器的参数,若该振动特性不满足工程需求,则还需要重复上述步骤,对集群式CMG隔振系统进行新一轮的更新,直至计算结果满足需求。通过该步骤可以验证隔振器/阻尼器参数的有效性,而有效的参数在后续瞬态分析和实体仿真中至关重要,甚至是最终实验验证的基础。
本发明的一种较优实施例中,步骤1中,动力学模型包括集群式控制力矩陀螺的动能和隔振器的弹性势能,
集群式控制力矩陀螺的动能E表示为;
式中,Ix、Iy和Iz表示集群式控制力矩陀螺所处的平台质心且平行于坐标系轴的转动惯量,所述平台包括集群式控制力矩陀螺和支架,m表示平台的质量,和是集群式控制力矩陀螺质心的三个平动位移xc、yc和zc的导数,是平台质心位置转动位移的导数,
其中,集群式控制力矩陀螺质心的三个平动位移xc、yc和zc表示为:
式中,x、y和z是集群式控制力矩陀螺在全局坐标系原点位置的三个平动位移,x0、y0和z0是集群式控制力矩陀螺质心在全局坐标系下的坐标,α、β和γ是平台质心位置转动位移。
具体地,使用动力学建模获取固有频率的方式,能先大致确定最优参数的大致分布域,避免了直接使用有限元计算而带来的大规模协同优化计算,有助于加快设计进度。此外,与已有方法采用柔性铰连接的不同,本方案采用隔振器固接的安装方式,使得隔振器的安装对冲击的抗性更强、可靠性更高。而刚体的转动描述在所有平行的坐标系是相同的,所以坐标系的转动分量大小等于质心的转动分量大小,于是集群式CMG的动能为公式2。
本发明的一种较优实施例中,单个隔振器的弹性势能Epi为:
式中,v为设定的隔振器上弹簧与集群式控制力矩陀螺连接点的位移向量,e为弹簧的轴向向量,vge为弹簧在轴向的位移,v-(vge)e为弹簧在轴向的位移与轴向垂直方向的位移,ka和kr分别为设定的弹簧轴向和径向刚度,e=ei,i=1,2,3,4,
四个隔振器上支点弹簧的轴向向量ei,i=1,2,3,4,分别表示为:
式中,θ1、θ2分别表示隔振器与全局坐标系x轴和z轴方向的夹角,
v由四个上支点的位移分量xi、yi、zi,i=1,2,3,4表示,xi、yi、zi,i=1,2,3,4,表示为:
式中,a、b分别表示隔振器上支点距全局坐标系x轴和y轴的垂线距离;
假设四个弹簧的刚度相同,则总的弹性势能Ep为:
式中,vi和ei分别表示第i个位移向量和第i个隔振器上支点弹簧的轴向向量;
将公式4和公式5代入公式6,得到:
具体地,为了计算系统的弹性势能,需要计算四个隔振器上端点的位移,即四个上支点的位移分量。由于集群式CMG的刚度要远大于隔振器,所以可根据刚体的运动学关系,可以得到四个上支点的位移分量xi、yi、zi。
本发明的一种较优实施例中,步骤3中,固有频率表示为:
|K-ω2M|=0 公式8,
式中,ω是减振系统的固有频率,K是刚度矩阵,M是质量矩阵,
具体地,质量矩阵是根据第二类拉格朗日方程简化得到的。
刚度矩阵是根据第二类拉格朗日方程简化得到的。
构建频率行列式,并令频率行列式等于0可以得到频率方程,如公式8所示,求解后即可得到减振系统的六阶固有频率,通过调整减振系统动力学模型中各个参数的取值,以使减振系统的固有频率满足需求。
本发明的一种较优实施例中,步骤3中,建立含有集群式控制力矩陀螺的有限元模型的具体过程为:
根据集群式控制力矩陀螺模型中支架空间梁的分布情况,在有限元模型中的有限元软件中选择梁单元并赋予梁单元的截面属性,构建支架模型,根据集群式控制力矩陀螺模型中多个控制力矩陀螺的运动特征,在有限元软件中分别用点单元描述运动特征中多个控制力矩陀螺与支架模型连接的固定部分和多个控制力矩陀螺转动部分,并定义固定部分和转动部分的连接属性,完成含有集群式控制力矩陀螺的有限元模型的建立。
具体地,根据控制力矩陀螺在空间上的分布位置及刚性支架的设计,建立集群式控制力矩陀螺的实体模型。此外,为了减少有限元分析的计算量,选择点单元对控制力矩陀螺的组件进行建模、选择梁单元对刚性支架进行建模,与实体模型中的参数进行对比,修正并验证有限元模型的合理性。
本发明的一种较优实施例中,步骤3中,根据每个所述结构参数设计一个弹簧模型和阻尼器模型的具体过程为:
将每个弹簧和阻尼器结构参数输入到CAD软件中,由CAD软件建立3D模型,将建立的3D模型再导入到有限元软件中计算此时弹簧和阻尼器的刚度,并与各自预设的刚度进行比较,如果弹簧和阻尼器的刚度与各自预设的刚度均相同,则将建立的3D模型作为弹簧模型和阻尼器模型,如果弹簧和阻尼器的刚度与各自预设的刚度有一个不相同,则修改CAD软件中的3D模型,并再次使用有限元软件计算修正后的3D模型中的弹簧和阻尼器的刚度,重复在CAD中修改3D模型和使用有限元软件计算修正后的3D模型中的弹簧和阻尼器的刚度的过程,直到修正后的最终3D模型中的弹簧和阻尼器的刚度与各自预设的刚度均相同,则将修正后的最终3D模型作为弹簧模型和阻尼器模型。
具体地,如果弹簧和阻尼器的刚度与各自预设的刚度均相同,则说明隔振器刚度或阻尼满足需求,将建立的3D模型作为弹簧模型和阻尼器模型。得到弹簧模型和阻尼器模型后,人为的选出一个弹簧模型和阻尼器模型,按照此模型分别选用横槽弹簧和圆柱形阻尼材料来组成隔振器,与支架连接,实现对应的隔振/阻尼性能。
具体地,根据多个CMG在空间上的分布位置及安装方式,设计一个刚性支架将它们连接成一个整体,建立集群式CMG的实体模型;而后,根据刚性支架中空间梁的分布情况,在有限元软件中选择梁单元并赋予相应的截面属性,来建立刚性支架的有限元模型,为了进一步减少有限元分析的计算量,还需要将各个CMG根据运动特征将其固定部分和转动部分分别用点单元描述,并定义连接属性;最后,对比有限元模型与实体模型之间质量、惯量、质心等信息的吻合度,以修正简化模型的合理性。该步骤的意义在于通过合理有效的途径简化有限元模型中单元类型的数量,从而只计算我们需要的参数,这能让集群式CMG有限元模型的计算效率更高,更利于后续的参数优化。
本发明的一种较优实施例中,所述弹簧结构参数包括弹簧长度、槽数和厚度,阻尼器结构参数包括阻尼系数、长度和阻尼材料耗散系数。
具体地,待确定隔振器刚度、阻尼、几何等参数后,首先对弹簧进行设计,根据横槽弹簧的长度、槽数、厚度等,对弹簧的三向刚度(一个轴向和两个径向)进行设计、仿真和优化,以此来修正弹簧的刚度参数使之与设计值相符。而对于阻尼器,采用圆柱形橡胶材料进行制造,考虑到圆柱形橡胶阻尼器的阻尼系数与长度、直径、材料耗散系数的比例关系,据此确定橡胶阻尼器的形状参数。该步骤的意义在于,通过选择合适构型的隔振弹簧,不仅能降低弹簧的设计和加工难度,减少加工成本,还有利于后续标准化作业,实现隔振弹簧的大批量生产。此外,阻尼器的构型选择也涉及到隔振阻尼参数表征的完好性,通过设计适合的橡胶阻尼器,能极大降低共振区域的响应幅值,对振动传递的衰减起着至关重要的作用。
本发明的一种较优实施例中,在步骤3和步骤4之间,还包括:
对步骤3得到的弹簧模型和阻尼器模型进行进一步验证。
本发明的一种较优实施例中,对步骤3得到的弹簧模型和阻尼器模型进行进一步验证,具体过程为:
有限元软件从步骤3中得到的多个弹簧模型和多个阻尼器模型中任意挑选一个弹簧模型和一个阻尼器模型与支架建立组合模型,直到完成所有弹簧模型和阻尼器与支架的组合情况,得到多个组合模型,每个组合模型输出一个固体频率,将固有频率满足预设固有频率范围的弹簧模型和阻尼器模型保留,从保留的弹簧模型和阻尼器模型中选择一个弹簧模型和阻尼器模型。
具体地,这步骤中比较固有频率的目的是为了看一下采用实体模型计算得到的固有频率是否仍然保持一致,虽然我们知道这肯定一致,但还是要重新计算一下进一步验证下。
Claims (8)
1.一种抑制调姿共振的集群式控制力矩陀螺隔振方法,其特征在于,所述方法包括以下步骤:
步骤1、支架安装在固定平面上,将多个控制力矩陀螺安装在支架上,根据多个控制力矩陀螺与支架的连接位置,构建集群式控制力矩陀螺实体模型,将4个隔振器安装在支架和固定平面之间,根据隔振器上支点和下支点分别与支架和固定平面的连接位置,构建动力学模型,每个隔振器由一个弹簧和沿弹簧伸缩方向穿在弹簧内部的阻尼器组成;
步骤2、根据集群式控制力矩陀螺实体模型和动力学模型得到固有频率,根据固有频率得到隔振器下支点空间坐标理论值;
步骤3、根据集群式控制力矩陀螺实体模型和动力学模型建立含有集群式控制力矩陀螺和隔振器的有限元模型,将步骤2中得到的隔振器下支点空间坐标理论值添加到所述有限元模型中,使用有限元软件对隔振器下支点空间坐标理论值进行解析,得到隔振器的固有频率,判断隔振器的固有频率是否处于预设固有频率范围,如果是,则根据固有频率得到最终的隔振器下支点空间坐标,如果否,则调整隔振器下支点空间坐标,直到计算的固定频率达到预设固有频率范围,根据最终生成的固有频率得到最终的隔振器下支点空间坐标;
根据最终的隔振器下支点空间坐标得到多个弹簧和阻尼器的结构参数,根据每个所述结构参数设计弹簧模型和阻尼器模型;
步骤4、从得到的多个弹簧模型和阻尼器模型任选出一个弹簧模型和阻尼器模型,按照选出的弹簧模型、阻尼器模型和步骤3中得到的最终的阻尼器下支点空间坐标调整步骤1中隔振器的结构和4个隔振器下支点与固定平面的连接位置,建立隔振器与集群式控制力矩陀螺的最终连接,实现对集群式控制力矩陀螺的隔振;
步骤3中,根据每个所述结构参数设计一个弹簧模型和阻尼器模型的具体过程为:
将每个弹簧和阻尼器结构参数输入到CAD软件中,由CAD软件建立3D模型,将建立的3D模型再导入到有限元软件中计算此时弹簧和阻尼器的刚度,并与各自预设的刚度进行比较,如果弹簧和阻尼器的刚度与各自预设的刚度均相同,则将建立的3D模型作为弹簧模型和阻尼器模型,如果弹簧和阻尼器的刚度与各自预设的刚度有一个不相同,则修改CAD软件中的3D模型,并再次使用有限元软件计算修正后的3D模型中的弹簧和阻尼器的刚度,重复在CAD中修改3D模型和使用有限元软件计算修正后的3D模型中的弹簧和阻尼器的刚度的过程,直到修正后的最终3D模型中的弹簧和阻尼器的刚度与各自预设的刚度均相同,则将修正后的最终3D模型作为弹簧模型和阻尼器模型;弹簧采用横槽弹簧,阻尼器采用橡胶阻尼器。
2.根据权利要求1所述一种抑制调姿共振的集群式控制力矩陀螺隔振方法,其特征在于,步骤1中,动力学模型包括集群式控制力矩陀螺的动能和隔振器的弹性势能,
集群式控制力矩陀螺的动能E表示为;
式中,Ix、Iy和Iz表示集群式控制力矩陀螺所处的平台质心且平行于坐标系轴的转动惯量,所述平台包括集群式控制力矩陀螺和支架,m表示平台的质量,和是集群式控制力矩陀螺质心的三个平动位移xc、yc和zc的导数,是平台质心位置转动位移的导数,
其中,集群式控制力矩陀螺质心的三个平动位移xc、yc和zc表示为:
式中,x、y和z是集群式控制力矩陀螺在全局坐标系原点位置的三个平动位移,x0、y0和z0是集群式控制力矩陀螺质心在全局坐标系下的坐标,α、β和γ是平台质心位置转动位移。
3.根据权利要求2所述一种抑制调姿共振的集群式控制力矩陀螺隔振方法,其特征在于,单个隔振器的弹性势能Epi为:
式中,v为设定的隔振器上弹簧与集群式控制力矩陀螺连接点的位移向量,e为弹簧的轴向向量,vge为弹簧在轴向的位移,v-(vge)e为弹簧在轴向的位移与轴向垂直方向的位移,ka和kr分别为设定的弹簧轴向和径向刚度,e=ei,i=1,2,3,4,
四个隔振器上支点弹簧的轴向向量ei,i=1,2,3,4,分别表示为:
式中,θ1、θ2分别表示隔振器与全局坐标系x轴和z轴方向的夹角,
v由四个上支点的位移分量xi、yi、zi,i=1,2,3,4表示,xi、yi、zi,i=1,2,3,4,表示为:
式中,a、b分别表示隔振器上支点距全局坐标系x轴和y轴的垂线距离;
假设四个弹簧的刚度相同,则总的弹性势能Ep为:
式中,vi和ei分别表示第i个位移向量和第i个隔振器上支点弹簧的轴向向量;
将公式4和公式5代入公式6,得到:
5.根据权利要求1所述一种抑制调姿共振的集群式控制力矩陀螺隔振方法,其特征在于,步骤3中,建立含有集群式控制力矩陀螺的有限元模型的具体过程为:
根据集群式控制力矩陀螺模型中支架空间梁的分布情况,在有限元模型中的有限元软件中选择梁单元并赋予梁单元的截面属性,构建支架模型,根据集群式控制力矩陀螺模型中多个控制力矩陀螺的运动特征,在有限元软件中分别用点单元描述运动特征中多个控制力矩陀螺与支架模型连接的固定部分和多个控制力矩陀螺转动部分,并定义固定部分和转动部分的连接属性,完成含有集群式控制力矩陀螺的有限元模型的建立。
6.根据权利要求1所述一种抑制调姿共振的集群式控制力矩陀螺隔振方法,其特征在于,所述弹簧结构参数包括弹簧长度、槽数和厚度,阻尼器结构参数包括阻尼系数、长度和阻尼材料耗散系数。
7.根据权利要求1所述一种抑制调姿共振的集群式控制力矩陀螺隔振方法,其特征在于,所述步骤3还包括:
对得到的弹簧模型和阻尼器模型进行验证。
8.根据权利要求7所述一种抑制调姿共振的集群式控制力矩陀螺隔振方法,其特征在于,对步骤3得到的弹簧模型和阻尼器模型进行验证,具体过程为:
有限元软件从步骤3中得到的多个弹簧模型和多个阻尼器模型中任意挑选一个弹簧模型和一个阻尼器模型与支架建立组合模型,直到完成所有弹簧模型和阻尼器与支架的组合情况,得到多个组合模型,每个组合模型输出一个固体频率,将固有频率满足预设固有频率范围的弹簧模型和阻尼器模型保留,从保留的弹簧模型和阻尼器模型中选择一个弹簧模型和阻尼器模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011111653.9A CN112231861B (zh) | 2020-10-16 | 2020-10-16 | 一种抑制调姿共振的集群式控制力矩陀螺隔振方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011111653.9A CN112231861B (zh) | 2020-10-16 | 2020-10-16 | 一种抑制调姿共振的集群式控制力矩陀螺隔振方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112231861A CN112231861A (zh) | 2021-01-15 |
CN112231861B true CN112231861B (zh) | 2021-07-06 |
Family
ID=74117770
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011111653.9A Expired - Fee Related CN112231861B (zh) | 2020-10-16 | 2020-10-16 | 一种抑制调姿共振的集群式控制力矩陀螺隔振方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112231861B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114838083B (zh) * | 2022-07-04 | 2022-11-15 | 中建安装集团有限公司 | 一种基于LoRa通信的分布式可变阻尼复合减振系统及减振方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109226753A (zh) * | 2018-09-20 | 2019-01-18 | 北京理工大学 | 基于3d打印技术制备钨颗粒增强金属基复合材料的方法 |
CN111323010A (zh) * | 2020-04-21 | 2020-06-23 | 桂林电子科技大学 | 一种室内外智能导航系统 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104166758B (zh) * | 2014-08-07 | 2017-03-22 | 东北大学 | 一种转子‑叶片耦合系统固有频率的确定方法 |
US9963978B2 (en) * | 2015-06-09 | 2018-05-08 | Ebert Composites Corporation | 3D thermoplastic composite pultrusion system and method |
CN106005484B (zh) * | 2016-05-18 | 2018-05-04 | 北京空间飞行器总体设计部 | 一种并联桁架式控制力矩陀螺群减隔振装置 |
CN106742066B (zh) * | 2016-12-29 | 2019-05-24 | 中国空间技术研究院 | 一种轻巧型星上执行机构隔振装置 |
JP2019028192A (ja) * | 2017-07-28 | 2019-02-21 | キヤノン株式会社 | 撮像装置 |
CN107559371B (zh) * | 2017-09-07 | 2020-09-18 | 北京控制工程研究所 | 一种基于磁流变技术的半主动隔振器 |
US11299981B2 (en) * | 2017-12-04 | 2022-04-12 | Hrl Laboratories, Llc | Continuous trajectory calculation for directional drilling |
CN109189086A (zh) * | 2018-07-27 | 2019-01-11 | 中国人民解放军战略支援部队航天工程大学 | 一种基于磁悬浮控制敏感陀螺的航天器姿态与振动一体化控制方法 |
-
2020
- 2020-10-16 CN CN202011111653.9A patent/CN112231861B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109226753A (zh) * | 2018-09-20 | 2019-01-18 | 北京理工大学 | 基于3d打印技术制备钨颗粒增强金属基复合材料的方法 |
CN111323010A (zh) * | 2020-04-21 | 2020-06-23 | 桂林电子科技大学 | 一种室内外智能导航系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112231861A (zh) | 2021-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104850697B (zh) | 基于ansys和adams的大型天线动力学建模方法 | |
CN109255188B (zh) | 基于有限元的六轴工业机器人动态性能优化方法 | |
CN111310380B (zh) | 一种电动车动力总成悬置橡胶衬套结构的设计开发方法 | |
CN108897973B (zh) | 一种弹簧-变截面盘-叶片系统的动力学建模方法 | |
CN110188417B (zh) | 一种基于多级超单元的复杂结构有限元模型修正方法 | |
CN108920740A (zh) | 悬置内骨架模态优化方法及系统 | |
US20160139577A1 (en) | Method for determining a movement limit | |
CN112231861B (zh) | 一种抑制调姿共振的集群式控制力矩陀螺隔振方法 | |
CN108416083B (zh) | 一种高耸电视塔结构二维动力模型分析方法及系统 | |
JPH0689322A (ja) | 最適設計システム | |
Patil et al. | A scalable time-parallel solution of periodic rotor dynamics in X3D | |
US20220035352A1 (en) | Virtualized Cable Modeling For Manufacturing Resource Simulation | |
CN102945002A (zh) | 基于非线性数学模型的通用型无人机仿真方法及系统 | |
CN111008417B (zh) | 一种大长细比连续结构的风致振动分析方法 | |
Wu | A robust optimization for the frequency and decoupling ratio of a powertrain mounting system based on interval analysis | |
CN111125829A (zh) | 一种优化全动平尾静气动弹性和颤振的方法 | |
JP2003156413A (ja) | 車体振動のシミュレーション方法 | |
CN111651923B (zh) | 一种多自由度被动隔振系统优化设计方法 | |
CN111814277B (zh) | 一种隔振平台设备及其动力学模型构建方法 | |
CN114818118A (zh) | 空气弹簧系统结构件设计方法、装置、设备及存储介质 | |
JP2002149713A (ja) | 6自由度シミュレーション方法 | |
TWI399496B (zh) | Two-way coupled tuned mass damper design method, computer program products and bi-directional coupled tuned mass damper | |
CN109726454B (zh) | 管路系统的流固耦合建模方法及装置 | |
CN103324032A (zh) | 一种多面体主基板及其制造方法和加工方法 | |
CN105259907A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210706 |