CN116956650B - 基于matlab app designer的姿控动力系统力学环境试验条件设计系统 - Google Patents
基于matlab app designer的姿控动力系统力学环境试验条件设计系统 Download PDFInfo
- Publication number
- CN116956650B CN116956650B CN202311219585.1A CN202311219585A CN116956650B CN 116956650 B CN116956650 B CN 116956650B CN 202311219585 A CN202311219585 A CN 202311219585A CN 116956650 B CN116956650 B CN 116956650B
- Authority
- CN
- China
- Prior art keywords
- module
- curve
- design
- inflection point
- psd
- 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.)
- Active
Links
- 238000013461 design Methods 0.000 title claims abstract description 94
- 238000012360 testing method Methods 0.000 title claims abstract description 70
- 238000004364 calculation method Methods 0.000 claims abstract description 59
- 230000004044 response Effects 0.000 claims abstract description 20
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 49
- 238000004140 cleaning Methods 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 10
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 230000008676 import Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 abstract description 5
- 230000001133 acceleration Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000000034 method Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012938 design process Methods 0.000 description 2
- 241000030538 Thecla Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000013102 re-test Methods 0.000 description 1
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
-
- 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)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Feedback Control In General (AREA)
Abstract
本发明提供一种基于matlab app designer的姿控动力系统力学环境试验条件设计系统，基于第一次力学环境试验以及系统试车试验得出的响应数据进行功率谱密度以及冲击响应谱分析得到每个测量点的频域响应曲线形成目标包络曲线，通过拐点计算模块的选择、曲线表格设计模块的设定经系统计算后将计算所得曲线与设计目标曲线显示在曲线显示模块。本发明在windows下，以matlab语言、matlab app designer为工具，仅需要选择对应选项以及输入关键参数即可实现对姿控动力系统力学环境试验条件的快速设计，并可视化呈现分析结果与结果对比。
Description
技术领域
本发明属于力学环境试验条件设计领域，特别涉及一种基于matlab appdesigner的姿控动力系统力学环境试验条件设计系统。
背景技术
在航天航空工程领域中，姿控动力系统的力学环境试验条件设计以及系统中搭载的各单机的力学环境试验条件设计是姿控动力系统总体设计环节中必不可少的一环，对于姿控动力系统的力学环境试验条件设计以及其搭载的各单机的力学环境试验条件设计对与姿控动力系统产品质量影响非常大，所设计得条件不论产生“过试验”还是“欠试验”都是不允许的。为了保证力学环境试验条件设计的准确性，设计人员需要进行大量的曲线数据处理以及复杂的公式计算，这不但增加了总体设计的难度，并且降低了设计效率。现有的商用软件中，针对力学环境试验条件设计的求解器少之又少，大部分基于matlab appdesigner工具实现的设计系统均未涉及航空航天姿控动力系统力学环境试验条件设计方向，并且在力学环境试验条件设计方向上，现有的设计系统及方法也较少并没有完整的、统一的路线。因此，需要一款操作简便，处理快捷，针对性强的姿控动力系统力学环境试验条件设计系统及方法来减少庞大的设计计算量，提高计算精度与效率，并且具有较强的通用性。
发明内容
本发明所解决的技术问题是提供一种基于matlab app designer的姿控动力系统力学环境试验条件设计系统，解决了人工计算需要进行大量的曲线数据处理以及复杂的公式计算难度大效率低以及目前无针对姿控动力系统力学环境试验条件的设计系统的问题。
本发明采用的技术方案是：一种基于matlab app designer的姿控动力系统力学环境试验条件设计系统，根据第一次力学环境试验以及系统试车试验得出的响应数据进行功率谱密度以及冲击响应谱分析得到每个测量点的频域响应曲线形成目标包络曲线，通过目标包络曲线对姿控动力系统力学环境试验条件设计系统进行系统参数修正用于下一次试车试验前姿控动力系统以及系统中的各单机所经过的力学环境试验，所述姿控动力系统力学环境试验条件设计系统包括PSD拐点数值计算模块、SRS拐点数值计算模块、曲线表格设计模块、曲线显示模块；
所述系统具体操作流程为：
步骤1、拐点计算模块的选择，根据发动机点火动作时序类型选取对应的PSD或SRS标签，并对应进入该标签下的拐点数值计算模块；
步骤2、正增益/负增益的选择，根据目标包络曲线选择需要进行正增益设定或负增益设定，并选择相应的正增益/负增益标签；
步骤3、根据步骤1和2的选择不同对应进入四种拐点计算模块其中之一，四种拐点计算模块分别为PSD正增益拐点数值计算、PSD负增益拐点数值计算、SRS正增益拐点数值计算和SRS负增益拐点数值计算，所述拐点计算模块包括增益设置模块、PSD或SRS上限/下限选择模块、PSD或SRS上限/下限设置模块、上限频率设置模块、下限频率设置模块、结果生成模块、结果显示模块和数值清理模块；
步骤4、进入拐点计算模块后的参数设定，
步骤4.1、在增益设置模块输入期望设计的正增益/负增益数值；
步骤4.2、在PSD上限/下限或SRS上限/下限选择模块中的下拉菜单选择数据类型Y_up或Y_down，并在PSD上限/下限或SRS上限/下限设置模块中输入对应设计PSD或设计SRS值；
步骤4.3、按照期望设计频率点分别在上限频率设置模块以及下限频率设置模块中输入对应的拐点频率f_up与f_down，在所有设计参数输入完毕后，按下结果生成模块按钮，结果会显示在下方对应的结果显示模块中；
步骤5、曲线表格设计模块设定，曲线表格设计模块包括频率设置模块和幅值设置模块，在对应步骤4计算结束后，按照所设计的各频率点频率值，在频率设置模块中分别输入对应频率点数值，并按照所涉及的PSD或SRS对应数值，在幅值设置模块中分别输入对应的幅值数值；
步骤6、曲线显示模块操作，所述曲线显示模块包括坐标区模块、绘制曲线模块和清空曲线模块，点击绘制曲线模块按钮后首先选择需要对比的目标文件PSD或SRS曲线，在曲线数据导入完成后，坐标区模块中会自动显示在曲线表格设计模块中输入的对应设计曲线以及之前所选择的对比文件曲线，点击清空曲线模块按钮后，能够清空当前显示的曲线数据，坐标区模块会显示为空白。
优选的，步骤3所述数值清理模块，点击数值清理模块按钮后，拐点计算模块中所有显示数据会清空并显示为0。
优选的，步骤3所述结果生成模块通过定义全局函数以及按钮回调函数互相调用实现，其中，全局函数通过app.xx.value获取目标参数的数值，通过switch case选择语句对app.dropdown.value选择的下拉菜单名称进行选择，当下拉菜单名称为Y_up时，计算幅值下限并赋值app.result.text从而实现结果显示模块和当前对应标签显示当前得到的结果以及结果标签名称功能；当下拉菜单名称为Y_down时，计算幅值上限并赋值app.result.text从而实现结果显示模块和当前对应标签显示当前得到的结果以及结果标签名称功能；按钮回调函数通过点击结果生成模块的结果生成按钮来调用目标全局函数用于实现结果的数值赋值。
优选的，步骤4.3所述结果显示模块前的对应标签会自动显示当前得到的结果标签名称。
优选的，步骤6所述绘制曲线模块通过定义全局函数以及按钮回调函数互相调用实现，其中全局函数通过matlab中uigetfile（）函数打开文件选择对话框来获取对比文件曲线文件格式为.txt的数据名称以及路径，并通过load（fullfile（））函数导入文件全部数据并设定第一列数据为频率，第二列数据为幅值，通过app.xx.value语句获取步骤5的曲线表格设计模块中输入的曲线设计参数，通过连续的if条件选择语句来决定曲线表格中使用的数据个数，并通过plot（）函数同时绘制对比文件曲线和设计曲线；按钮回调函数通过点击绘制曲线模块按钮来调用目标全局函数用于实现曲线的显示以及对比功能。
优选的，步骤3中所述PSD正/负增益拐点数值计算公式为
；
其中，N1为PSD目标包络曲线斜率，N1为正则为正增益，N1为负则为负增益；
Y1和Y2为PSD目标包络曲线的纵坐标值；
f1和f2为横坐标值，根据目标包络曲线拐点以及变化趋势得出，以PSD目标包络曲线的第一个峰值点的频率值取为f1，目标包络曲线的峰值最大点对应的频率值取为f2。
优选的，所述SRS正/负增益拐点数值计算公式为：
；
其中，N2为SRS目标包络曲线斜率，N2为正则为正增益，N2为负则为负增益；
G1和G2为SRS目标包络曲线纵坐标值；
f3和f4为横坐标值，根据目标包络曲线拐点以及变化趋势得出，以SRS目标包络曲线的第一个峰值点的频率值取为f3，目标包络曲线的峰值最大点对应的频率值取为f4。
本发明的有益效果是：
1、本发明在windows下，以matlab语言、matlab app designer为工具，仅需要选择对应选项以及输入关键参数即可实现对姿控动力系统力学环境试验条件的快速设计，并可视化呈现分析结果与结果对比。
2、由于现场条件的非确定性，在一些系统监测主机中无法保证软件安装的广泛性，而使用matlab app designer进行设计后的软件可以进行自定义封装并在没有安装matlab的计算机上使用，增加了设计系统软件的通用性。
附图说明
图1为基于matlab app designer的姿控动力系统力学环境试验条件设计系统流程图；
图2为基于matlab app designer的姿控动力系统力学环境试验条件设计系统计算前界面示意图；
图3为基于matlab app designer的姿控动力系统力学环境试验条件设计系统计算后界面示意图；
图4为PSD正/负增益拐点数值计算公式示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步的详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
在姿控动力系统的设计过程中，设计人员会首先接收到任务文件，其中在力学方面的要求是设计过程中的重要一环，姿控动力系统的设计要满足任务文件提出的力学环境试验要求，然而力学环境试验要求往往是经验值，由于每个系统存在差异性，任务文件中的力学环境试验条件可能存在误差，所以设计人员需要通过设计-生产-试验-再设计-再试验的闭环流程反复修正力学环境试验条件，以保证姿控动力系统的力学环境试验条件的准确性。
在姿控动力系统第一次组装完成后会进行第一次力学环境试验以及系统试车试验，试验时会在姿控动力系统的多个测量点安装传感器以测得试验的响应数据，大多数时候为加速度值，设计人员会根据试验后的响应数据进行功率谱密度（PSD- power spectraldensity）以及冲击响应谱（SRS-shock response spectrum）分析得到每个测量点的频域响应曲线，根据频域曲线的实际情况重新设计修改姿控动力系统力学环境试验条件曲线，从而保证下一次试车试验前姿控动力系统以及系统中的各单机所经过的力学环境试验是更加准确的，而不是过试验或欠试验的。
所谓的包络曲线就是在姿控动力系统进行试验后得到的时域响应曲线经设计人员分析处理后得到的频域响应曲线。
一种基于matlab app designer的姿控动力系统力学环境试验条件设计系统，包括PSD拐点数值计算模块、SRS拐点数值计算模块、曲线表格设计模块、曲线显示模块。
PSD拐点数值计算模块、SRS拐点数值计算模块为二选一的标签模式。
在PSD拐点数值计算模块中包括正增益设计模块或负增益设计模块二选一的标签模式。
选定正增益或负增益标签后显示该具体计算模块界面，如图2左上部分所示，以PSD正增益拐点数值计算模块为例，模块包括：增益设置模块、PSD上限/下限选择模块、PSD上限/下限设置模块、上限频率设置模块、下限频率设置模块、结果生成模块、结果显示模块、数值清理模块。
同理，PSD负增益拐点数值计算模块包括：增益设置模块、PSD上限/下限选择模块、PSD上限/下限设置模块、上限频率设置模块、下限频率设置模块、结果生成模块、结果显示模块、数值清理模块。
SRS正增益拐点数值计算模块包括：增益设置模块、SRS上限/下限选择模块、SRS上限/下限设置模块、上限频率设置模块、下限频率设置模块、结果生成模块、结果显示模块、数值清理模块。
SRS负增益拐点数值计算模块包括：增益设置模块、SRS上限/下限选择模块、SRS上限/下限设置模块、上限频率设置模块、下限频率设置模块、结果生成模块、结果显示模块、数值清理模块。
曲线表格设计模块包括：频率设置模块、幅值设置模块。
曲线显示模块包括：PSD曲线显示模块、SRS曲线显示模块。其中，PSD曲线显示模块包括：坐标区模块、绘制曲线模块、清空曲线模块。SRS曲线显示模块包括：坐标区模块、绘制曲线模块、清空曲线模块。
系统使用流程如图1所示。
步骤1、拐点计算模块的选择，根据发动机点火动作时序类型选取对应的PSD或SRS标签，并对应进入该标签下的拐点数值计算模块；发动机点火动作时序类型：电爆阀起爆--SRS；预点火--PSD；长程点火--PSD；脉冲点火--PSD。
步骤2、正增益/负增益的选择，根据目标包络曲线选择需要进行正增益设定或负增益设定，并选择相应的正增益/负增益标签；
步骤3、根据步骤1和2的选择不同对应进入四种拐点计算模块其中之一，四种拐点计算模块分别为PSD正增益拐点数值计算、PSD负增益拐点数值计算、SRS正增益拐点数值计算和SRS负增益拐点数值计算，下述步骤以选择PSD正增益拐点数值计算模块为例；
步骤4、进入PSD正增益拐点数值计算模块后的参数设定，如图2左上角数据界面所示。
步骤4.1、在增益设置模块输入期望设计的正增益数值，默认输入为正整数。
步骤4.2、在PSD上限/下限选择模块中的下拉菜单选择数据类型Y_up或Y_down并在PSD上限/下限设置模块中输入对应设计PSD值；一般来说，PSD上限/下限选择，就是根据目标包络曲线的第一个峰值和最大峰值点所对应的纵坐标值来向上取整确定，即下限为第一个峰值点纵坐标值向上取整；上限为最大峰值点纵坐标值向上取整。当PSD上限取值时PSD下限为对应计算所得，同理，当PSD下限取值时PSD上限为对应计算所得。PSD上限/下限选择模块中的下拉菜单选择通过matlab app designer自带的UI组件中的app.dropdown组件将下拉菜单设置为两行一列并分别命名为Y_up（g^2/hz）和Y_down（g^2/hz）。
步骤4.3、按照期望设计频率点分别在上限频率设置模块以及下限频率设置模块中输入对应的拐点频率（f_up与f_down），在所有设计参数输入完毕后，按下结果生成模块按钮，结果会显示在下方对应的结果显示模块中；
结果生成模块主要通过定义全局函数以及按钮回调函数互相调用实现。其中，全局函数通过app.xx.value获取目标参数数值，通过switch case选择语句对app.dropdown.value选择的下拉菜单名称进行选择，当下拉菜单名称为Y_up（g^2/hz）时，计算幅值下限并赋值app.result.text从而实现结果显示模块和当前对应标签显示当前得到的结果以及结果标签名称功能；当下拉菜单Y_down（g^2/hz）时，计算幅值上限并赋值app.result.text从而实现结果显示模块和当前对应标签显示当前得到的结果以及结果标签名称功能。
按钮回调函数通过鼠标按下结果生成按钮来调用目标全局函数从而实现结果的数值赋值。
点击数值清理模块按钮后，拐点计算模块中所有显示数据会清空并显示为0。数值清理模块主要通过app.xx.value语句将目标显示框内数值设为0来实现。
步骤5、曲线表格设计模块设定，在对应步骤4计算结束后，按照所设计的各频率点频率值，在频率设置模块中分别输入对应频率点数值，并按照所涉及的PSD对应数值，在幅值设置模块中分别输入对应的幅值数值。
步骤6、曲线显示模块操作，如图3所示，所述曲线显示模块包括坐标区模块、绘制曲线模块和清空曲线模块，点击绘制曲线模块按钮后首先选择需要对比的目标文件PSD曲线，默认为.txt文件格式，在曲线数据导入完成后，坐标区模块中会自动显示在曲线表格设计模块中输入的对应设计曲线以及之前所选择的对比文件曲线，其中由线段组成的阶段线为本系统计算得出的设计曲线，波浪曲线为对比文件曲线，系统设计曲线与对比文件曲线越贴近说明曲线包络的越好，即试验条件与实际产生的响应越接近。设计曲线就是为了在进行力学环境试验时所施加的试验条件（激励）尽可能接近试车时由于系统工作带来的响应值，这样可以提高姿控动力系统设计成功率。力学环境试验就是在模拟姿控动力系统实际工作中可能遇到的各种情况，如果力学环境试验条件越接近实际的工况，那么在正式试车时成功率越高。
点击清空曲线模块按钮后，能够清空当前显示的曲线数据，坐标区模块会显示为空白。
绘制曲线模块主要通过定义全局函数以及按钮回调函数互相调用实现。其中全局函数通过matlab中uigetfile（）函数打开文件选择对话框来获取对比文件曲线格式为.txt的数据名称以及路径，并通过load（fullfile（））函数导入文件全部数据并设定第一列数据为频率，第二列数据为幅值。与此同时，通过app.xx.value语句获取曲线设计表格中输入的曲线设计参数。最后通过连续的if条件选择语句来决定曲线表格中使用的数据个数，并通过plot（）函数同时绘制对比文件曲线与设计曲线。
按钮回调函数通过鼠标按下绘制曲线按钮来调用目标全局函数从而实现曲线的显示以及对比功能。
清空曲线模块通过cla（）函数清空目标区域曲线。
如图4所示，随机振动加速度功率谱密度（PSD）计算公式如下所示：
；
其中，N1为PSD目标包络曲线斜率（为已知量），N1为正则为正增益，N1为负则为负增益；
Y1和Y2为PSD目标包络曲线的纵坐标值；
f1和f2为横坐标值（为已知量），根据PSD目标包络曲线拐点以及变化趋势得出，以目标包络曲线的第一个峰值点的频率值取为f1，目标包络曲线的峰值最大点对应的频率值取为f2。一般来说，在基于随机振动分析数据进行力学环境试验条件曲线设计时，横坐标值的设定会根据目标包络曲线拐点以及变化趋势来决定，即f2、f1均为已知。当N1为正整数时（正增益）可推导Y1表达式为：
；
同理可得Y2表达式为：；
Y1和Y2其中一个为已知，根据公式计算另一个，已知的是根据目标包络曲线幅值变化趋势决定。
当N1为负整数时同理可得。
同理，冲击响应谱（SRS）加速度值的计算公式如下所示：
；
其中，N2为SRS目标包络曲线斜率，N2为正则为正增益，N2为负则为负增益；
G1和G2为SRS曲线纵坐标值；
f3和f4为横坐标值，根据目标包络曲线拐点以及变化趋势得出，通常以目标包络曲线的第一个峰值点的频率值取为f3，目标包络曲线的峰值最大点对应的频率值取为f4。一般来说，在基于冲击响应谱分析数据进行力学环境试验条件曲线设计时，横坐标值的设定会根据目标包络曲线拐点以及变化趋势来决定，即f3、f4均为已知。当N2为正整数时（正增益）可推导G1表达式为：
；
同理可得G2表达式为：
；
G1和G2其中一个为已知，根据公式计算另一个，已知的是根据目标包络曲线幅值变化趋势决定。
当N2为负整数时同理可得。
当步骤3选择结果为PSD负增益拐点数值计算流程为：在选择了PSD负增益拐点数值计算标签后，首先在增益设置模块输入期望设计的负增益数值，输入时应注意需要在输入的正整数前添加负号，随后在PSD上限/下限选择模块中的下拉菜单选择已知的数据类型（Y_up或Y_down）并在PSD上限/下限设置模块中输入对应设计PSD值，接下来需要按照期望设计频率点分别在上限频率设置模块以及下限频率设置模块中输入对应的拐点频率（f_up与f_down），在所有设计参数输入完毕后，按下结果生成模块按钮，结果会显示在下方对应显示模块中。结果显示模块前对应标签会自动显示当前得到的结果标签名称。
SRS正增益拐点数值计算：
当步骤3选择结果为SRS正增益拐点数值计算流程为：在选择了SRS正增益拐点数值计算标签后，首先在增益设置模块输入期望设计的正增益数值，默认输入为正整数，随后在SRS上限/下限选择模块中的下拉菜单选择已知的数据类型Y_up或Y_down并在SRS上限/下限设置模块中输入对应设计SRS值，接下来需要按照期望设计频率点分别在上限频率设置模块以及下限频率设置模块中输入对应的拐点频率f_up与f_down，在所有设计参数输入完毕后，按下结果生成模块按钮，结果会显示在下方对应的结果显示模块中。结果显示模块前对应标签会自动显示当前得到的结果标签名称。
当步骤3选择结果为SRS负增益拐点数值计算流程为：在选择了SRS负增益拐点数值计算标签后，首先在增益设置模块输入期望设计的负增益数值，输入时应注意需要在输入的正整数前添加负号，随后在SRS上限/下限选择模块中的下拉菜单选择已知的数据类型（Y_up或Y_down）并在SRS上限/下限设置模块中输入对应设计SRS值，接下来需要按照期望设计频率点分别在上限频率设置模块以及下限频率设置模块中输入对应的拐点频率（f_up与f_down），在所有设计参数输入完毕后，按下结果生成模块按钮，结果会显示在下方对应的结果显示模块中。结果显示模块前对应标签会自动显示当前得到的结果标签名称。
以上所述是本发明的具体实施例及所运用的技术原理，任何基于本发明技术方案基础上的任何修改、等效变换，均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于matlab app designer的姿控动力系统力学环境试验条件设计系统，其特征在于：根据第一次力学环境试验以及系统试车试验得出的响应数据进行功率谱密度以及冲击响应谱分析得到每个测量点的频域响应曲线形成目标包络曲线，通过目标包络曲线对姿控动力系统力学环境试验条件设计系统进行系统参数修正用于下一次试车试验前姿控动力系统以及系统中的各单机所经过的力学环境试验，所述姿控动力系统力学环境试验条件设计系统包括PSD拐点数值计算模块、SRS拐点数值计算模块、曲线表格设计模块、曲线显示模块；
所述系统具体操作流程为：
步骤1、拐点计算模块的选择，根据发动机点火动作时序类型选取对应的PSD或SRS标签，并对应进入该标签下的拐点数值计算模块；
步骤2、正增益/负增益的选择，根据目标包络曲线选择需要进行正增益设定或负增益设定，并选择相应的正增益/负增益标签；
步骤3、根据步骤1和2的选择不同对应进入四种拐点计算模块其中之一，四种拐点计算模块分别为PSD正增益拐点数值计算、PSD负增益拐点数值计算、SRS正增益拐点数值计算和SRS负增益拐点数值计算，所述拐点计算模块包括增益设置模块、PSD或SRS上限/下限选择模块、PSD或SRS上限/下限设置模块、上限频率设置模块、下限频率设置模块、结果生成模块、结果显示模块和数值清理模块；
步骤4、进入拐点计算模块后的参数设定，
步骤4.1、在增益设置模块输入期望设计的正增益/负增益数值；
步骤4.2、在PSD上限/下限或SRS上限/下限选择模块中的下拉菜单选择数据类型Y_up或Y_down，并在PSD上限/下限或SRS上限/下限设置模块中输入对应设计PSD或设计SRS值；
步骤4.3、按照期望设计频率点分别在上限频率设置模块以及下限频率设置模块中输入对应的拐点频率f_up与f_down，在所有设计参数输入完毕后，按下结果生成模块按钮，结果会显示在下方对应的结果显示模块中；
步骤5、曲线表格设计模块设定，曲线表格设计模块包括频率设置模块和幅值设置模块，在对应步骤4计算结束后，按照所设计的各频率点频率值，在频率设置模块中分别输入对应频率点数值，并按照所涉及的PSD或SRS对应数值，在幅值设置模块中分别输入对应的幅值数值；
步骤6、曲线显示模块操作，所述曲线显示模块包括坐标区模块、绘制曲线模块和清空曲线模块，点击绘制曲线模块按钮后首先选择需要对比的目标文件PSD或SRS曲线，在曲线数据导入完成后，坐标区模块中会自动显示在曲线表格设计模块中输入的对应设计曲线以及之前所选择的对比文件曲线，点击清空曲线模块按钮后，能够清空当前显示的曲线数据，坐标区模块会显示为空白。
2.根据权利要求1所述的基于matlab app designer的姿控动力系统力学环境试验条件设计系统，其特征在于：步骤3所述数值清理模块，点击数值清理模块按钮后，拐点计算模块中所有显示数据会清空并显示为0。
3.根据权利要求1所述的基于matlab app designer的姿控动力系统力学环境试验条件设计系统，其特征在于：步骤3所述结果生成模块通过定义全局函数以及按钮回调函数互相调用实现，其中，全局函数通过app.xx.value获取目标参数的数值，通过switch case选择语句对app.dropdown.value选择的下拉菜单名称进行选择，当下拉菜单名称为Y_up时，计算幅值下限并赋值app.result.text从而实现结果显示模块和当前对应标签显示当前得到的结果以及结果标签名称功能；当下拉菜单名称为Y_down时，计算幅值上限并赋值app.result.text从而实现结果显示模块和当前对应标签显示当前得到的结果以及结果标签名称功能；按钮回调函数通过点击结果生成模块的结果生成按钮来调用目标全局函数用于实现结果的数值赋值。
4.根据权利要求1所述的基于matlab app designer的姿控动力系统力学环境试验条件设计系统，其特征在于：步骤4.3所述结果显示模块前的对应标签会自动显示当前得到的结果标签名称。
5.根据权利要求1所述的基于matlab app designer的姿控动力系统力学环境试验条件设计系统，其特征在于：步骤6所述绘制曲线模块通过定义全局函数以及按钮回调函数互相调用实现，其中全局函数通过matlab中uigetfile（）函数打开文件选择对话框来获取对比文件曲线文件格式为.txt的数据名称以及路径，并通过load（fullfile（））函数导入文件全部数据并设定第一列数据为频率，第二列数据为幅值，通过app.xx.value语句获取步骤5的曲线表格设计模块中输入的曲线设计参数，通过连续的if条件选择语句来决定曲线表格中使用的数据个数，并通过plot（）函数同时绘制对比文件曲线和设计曲线；按钮回调函数通过点击绘制曲线模块按钮来调用目标全局函数用于实现曲线的显示以及对比功能。
6.根据权利要求1所述的基于matlab app designer的姿控动力系统力学环境试验条件设计系统，其特征在于，步骤3中所述PSD正/负增益拐点数值计算公式为：
；
其中，N1为PSD目标包络曲线斜率，N1为正则为正增益，N1为负则为负增益；
Y1和Y2为PSD目标包络曲线的纵坐标值；
f1和f2为横坐标值，根据目标包络曲线拐点以及变化趋势得出，以PSD目标包络曲线的第一个峰值点的频率值取为f1，目标包络曲线的峰值最大点对应的频率值取为f2。
7.根据权利要求1所述的基于matlab app designer的姿控动力系统力学环境试验条件设计系统，其特征在于，所述SRS正/负增益拐点数值计算公式为：
；
其中，N2为SRS目标包络曲线斜率，N2为正则为正增益，N2为负则为负增益；
G1和G2为SRS目标包络曲线纵坐标值；
f3和f4为横坐标值，根据目标包络曲线拐点以及变化趋势得出，以SRS目标包络曲线的第一个峰值点的频率值取为f3，目标包络曲线的峰值最大点对应的频率值取为f4。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311219585.1A CN116956650B (zh) | 2023-09-21 | 2023-09-21 | 基于matlab app designer的姿控动力系统力学环境试验条件设计系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311219585.1A CN116956650B (zh) | 2023-09-21 | 2023-09-21 | 基于matlab app designer的姿控动力系统力学环境试验条件设计系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116956650A CN116956650A (zh) | 2023-10-27 |
| CN116956650B true CN116956650B (zh) | 2023-11-28 |
Family
ID=88462469
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311219585.1A Active CN116956650B (zh) | 2023-09-21 | 2023-09-21 | 基于matlab app designer的姿控动力系统力学环境试验条件设计系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116956650B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998014765A1 (en) * | 1996-09-30 | 1998-04-09 | Ford Motor Company | Method to specify random vibration tests for product durability validation |
| CN104121985A (zh) * | 2013-04-29 | 2014-10-29 | 艾默生电气（美国）控股公司（智利）有限公司 | 过采样数据的选择性抽取和分析 |
| CN106468951A (zh) * | 2016-08-29 | 2017-03-01 | 华东师范大学 | 一种基于双手环传感器融合的智能遥感控制系统及其方法 |
| CN107679290A (zh) * | 2017-09-14 | 2018-02-09 | 北京空间飞行器总体设计部 | 一种多学科微振动评估优化设计方法及系统 |
| CN107688554A (zh) * | 2017-09-01 | 2018-02-13 | 南京理工大学 | 基于自适应傅里叶分解的频域辨识方法 |
| CN110276256A (zh) * | 2019-05-17 | 2019-09-24 | 中国人民解放军战略支援部队信息工程大学 | 基于调参随机共振的低信噪比通信信号调制识别方法及装置 |
| CN111964860A (zh) * | 2020-06-30 | 2020-11-20 | 中国空间技术研究院 | 一种新型冲击响应数据判读方法 |
| CN114239131A (zh) * | 2021-11-24 | 2022-03-25 | 沈阳航天新光集团有限公司 | 姿控动力系统空间飞行力学环境预示方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8542002B2 (en) * | 2009-07-31 | 2013-09-24 | Tektronix, Inc. | System and method for performing power spectral density and power level measurements |
| US9778080B2 (en) * | 2013-04-29 | 2017-10-03 | Emerson Electric (Us) Holding Corporation (Chile) Limitada | Selective decimation and analysis of oversampled data |
-
2023
- 2023-09-21 CN CN202311219585.1A patent/CN116956650B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998014765A1 (en) * | 1996-09-30 | 1998-04-09 | Ford Motor Company | Method to specify random vibration tests for product durability validation |
| CN104121985A (zh) * | 2013-04-29 | 2014-10-29 | 艾默生电气（美国）控股公司（智利）有限公司 | 过采样数据的选择性抽取和分析 |
| CN106468951A (zh) * | 2016-08-29 | 2017-03-01 | 华东师范大学 | 一种基于双手环传感器融合的智能遥感控制系统及其方法 |
| CN107688554A (zh) * | 2017-09-01 | 2018-02-13 | 南京理工大学 | 基于自适应傅里叶分解的频域辨识方法 |
| CN107679290A (zh) * | 2017-09-14 | 2018-02-09 | 北京空间飞行器总体设计部 | 一种多学科微振动评估优化设计方法及系统 |
| CN110276256A (zh) * | 2019-05-17 | 2019-09-24 | 中国人民解放军战略支援部队信息工程大学 | 基于调参随机共振的低信噪比通信信号调制识别方法及装置 |
| CN111964860A (zh) * | 2020-06-30 | 2020-11-20 | 中国空间技术研究院 | 一种新型冲击响应数据判读方法 |
| CN114239131A (zh) * | 2021-11-24 | 2022-03-25 | 沈阳航天新光集团有限公司 | 姿控动力系统空间飞行力学环境预示方法 |
Non-Patent Citations (4)
| Title |
|---|
| 基于电液伺服控制的车辆主动悬架研究;赵丁选;史舒麟;杨彬;李元鹏;王大壮;;中国公路学报(06);全文 * |
| 岩石破坏声发射时频分析算法与瞬时频率前兆研究;宫宇新;何满潮;汪政红;尹雨婷;;岩石力学与工程学报(04);全文 * |
| 星敏感器支撑结构多目标拓扑优化设计与试验;程龙;王栋;谷松;高飞;杨林;李林;;红外与激光工程(05);全文 * |
| 航天器微振动测试、隔离、抑制技术综述;谭天乐;朱春艳;朱东方;宋婷;孙宏丽;顾;杨雨;;上海航天(06);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116956650A (zh) | 2023-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Park et al. | Unstructured grid adaptation: status, potential impacts, and recommended investments towards CFD 2030 | |
| CN107391891A (zh) | 一种基于模型融合方法的大展弦比机翼优化设计方法 | |
| CN104443427A (zh) | 飞行器颤振预测系统及方法 | |
| CN111159813B (zh) | 一种基于仿真的上浪砰击下船舶流固耦合分析方法和系统 | |
| CN106650059A (zh) | 一种基于Matlab的汽车整车性能仿真与计算方法 | |
| CN102087670A (zh) | 基于多体动力学的汽车悬架与减振系统仿真试验及设计平台 | |
| EP3425466A1 (en) | Design support device, design support method and design support program | |
| CN106777820A (zh) | 一种机床基础大件动态性能优化设计方法 | |
| CN105138770A (zh) | 基于间接可靠性特征量的航天产品可靠性仿真评估方法 | |
| CN112506759A (zh) | 伺服系统控制软件的自动化测试方法、装置及存储介质 | |
| US9378315B1 (en) | Method for semiconductor process corner sweep simulation based on value selection function | |
| CN102375913B (zh) | 采用主动结构刚度法生成机翼的方法 | |
| CN116956650B (zh) | 基于matlab app designer的姿控动力系统力学环境试验条件设计系统 | |
| CN102289530A (zh) | 一种起重机通用型主梁参数化计算机建模系统 | |
| CN108228455A (zh) | 一种软件控制危险分析方法 | |
| CN103440391B (zh) | 一种基于数值选择函数的半导体工艺角扫描仿真方法 | |
| CN108090288A (zh) | 一种通过机器学习获取时序参数的方法 | |
| KR101233312B1 (ko) | 발전소 성능평가 시뮬레이션 장치 및 그 방법 | |
| CN109940462A (zh) | 铣刀切削振动变化特性的检测与高斯过程模型构建方法 | |
| CN103455663B (zh) | 新能源客车切氢切电安全系统的传感器布置方法 | |
| CN110321596A (zh) | 一种基于有限元分析的机车车辆结构仿真方法 | |
| CN115659711A (zh) | 一种仿真方法、装置、电子设备及存储介质 | |
| US8326589B2 (en) | Stable equilibrium point (SEP) calculation apparatus of power system | |
| CN102160057B (zh) | 用于粒子法的界面粒子的判定方法及装置 | |
| JP2014016830A (ja) | 電力見積支援プログラム、電力見積支援装置および電力見積支援方法 |
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 |