CN115165298B - 一种旋转轴天平实时动态载荷监控系数的监测方法 - Google Patents
一种旋转轴天平实时动态载荷监控系数的监测方法 Download PDFInfo
- Publication number
- CN115165298B CN115165298B CN202211098375.7A CN202211098375A CN115165298B CN 115165298 B CN115165298 B CN 115165298B CN 202211098375 A CN202211098375 A CN 202211098375A CN 115165298 B CN115165298 B CN 115165298B
- Authority
- CN
- China
- Prior art keywords
- stress
- load
- balance
- real
- rotary shaft
- 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
- 238000012544 monitoring process Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000004590 computer program Methods 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 5
- 238000005096 rolling process Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000005452 bending Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
- G01M9/062—Wind tunnel balances; Holding devices combined with measuring arrangements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L19/00—Details of, or accessories for, apparatus for measuring steady or quasi-steady pressure of a fluent medium insofar as such details or accessories are not special to particular types of pressure gauges
- G01L19/08—Means for indicating or recording, e.g. for remote indication
- G01L19/12—Alarms or signals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Fluid Mechanics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明是一种旋转轴天平实时动态载荷监控系数的监测方法。本发明涉及旋转轴天平载荷监控技术领域,本发明根据旋转轴天平的应力状态,确定载荷作用时的应力分量和MISES应力;根据旋转天平的应力分量,测量旋转天平中心的受载,确定载荷点所受的应力状态;确定六种基元载荷分别单独作用时,对应的6个应力分量,计算线弹性情况下旋转轴天平上某个指定节点在任意载荷下的MISES应力及应力分量;根据节点的载荷系数得到实时状态下该点的实时应力值,将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时,则报警,对旋转轴天平进行载荷预警。本发明可实现判断天平实时应力状态,极大减小了旋转轴天平在试验过程中发生事故的机率。
Description
技术领域
本发明涉及旋转轴天平载荷监控技术领域,是一种旋转轴天平实时动态载荷监控系数的监测方法。
背景技术
旋转轴天平是复杂、精密的结构装配体。在设计过程中,对气动载荷的预估偏差会对天平与风洞试验带来风险挑战。气动载荷预估偏大会导致天平灵敏度下降,气动载荷预估偏小导致天平设计载荷无法满足试验需求,天平超载,有损坏危险。这种现象暴露出了一个问题,对旋转轴天平的全面的实时载荷监控不到位。通过对旋转轴天平进行载荷监控可以判断天平实时应力状态,一旦发生应力超过预警值得情况,可以立即采取必要措施,能够有效避免安全事故的发生。
发明内容
本发明为克服现有技术的不足,本发明通过研究旋转轴天平实时动态载荷监控系数,确定旋转轴天平载荷监控系数。应用载荷系数计算方法得出监测应力值,衡量旋转轴天平在不同载荷条件下天平是否超载。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
一种旋转轴天平实时动态载荷监控系数的监测方法,所述方法包括以下步骤:
步骤1:根据旋转轴天平的应力状态,确定载荷作用时的应力分量和MISES应力;
步骤2:根据旋转天平的应力分量,测量旋转天平中心的受载,确定载荷点所受的应力状态;
步骤3:根据步骤2确定六种基元载荷分别单独作用时,对应的6个应力分量,计算线弹性情况下旋转轴天平上某个指定节点在任意载荷下的MISES应力及应力分量;
步骤4:根据节点的载荷系数得到实时状态下该点的实时应力值,将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时,则报警,对旋转轴天平进行载荷预警。
优选地,所述步骤1具体为:
根据叠加原理可知,确定载荷作用时的应力分量和MISES应力,应力分量和MISES应力通过下式确定:
其中,为x方向应力分量,、为不同载荷单独作用下的x方向应力分量,为y方向应力分量,、为不同载荷单独作用下y方向应力分量,为剪切应力分量,、为不同载荷单独作用下剪切应力分量;通过、、、、、来计算载荷作用时的应力分量和MISES应力。
优选地,所述步骤2具体为:
测出天平中心受载通过下式表示:
其中,x为天平受轴向力,y为天平受升力,z为天平受侧向力,为x方向受力系数,为y方向受力系数,为z方向受力系数,L为天平受滚转力矩,为滚转方向载荷系数,M为天平受偏航力矩,为偏航方向载荷系数,N'为天平受俯仰力矩,为俯仰方向载荷系数,N·m为扭矩单位牛米,N为单位牛顿;
其中,为受轴向力x方向产生的应力分量,为受升力x方向产生的应力分量,为受侧向力x方向产生的应力分量,为受滚转力矩x方向产生的应力分量,为受偏航力矩x方向产生的应力分量,为受俯仰力矩x方向产生的应力分量,、、、、、分别为受六元载荷下y方向产生的应力分量,、、、、、分别为受六元载荷下z方向产生的应力分量,、、、、、分别为受六元载荷下xy方向的剪切应力分量,、、、、、分别为六元载荷下yz方向的剪切应力分量,、、、、、为分别受六元载荷下xz方向的剪切应力分量,为天平受到的实际应力。
优选地,所述步骤3具体为:
此时的应力分量为:
同理可得:
优选地,所述步骤4具体为:
旋转轴天平在使用时,通过设计载荷得到各个需要注意位置节点的应力最大值即应力阈值,通过加载基元载荷,根据旋转轴天平的实时输出,根据天平所受的基元载荷,,,,,,通过得到节点的载荷系数,实时状态下该点的实时应力值,将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时报警,即对旋转轴天平进行载荷预警。
优选地,当计算X=2000N时某点的应力,只需X=1N时6个应力分量乘以2000即得X=2000N时的六个应力分量;
当计算Y=200N,L=10N*M时某点MISES应力,只需Y=1N时六个应力分量乘以200,对应加上L=1N*M时的六个应力分量乘以10的值,就得到最终的应力分量,最后计算得到MISES应力值。
一种旋转轴天平实时动态载荷监控系数的监测系统,所述系统包括:
应力确定模块,所述应力确定模块根据旋转轴天平的应力状态,确定载荷作用时的应力分量和MISES应力;
应力状态确定模块,所述应力状态确定模块根据旋转天平的应力分量,测量旋转天平中心的受载,确定载荷点所受的应力状态;
天平应力分量确定模块,所述天平应力分量确定模块确定六种基元载荷分别单独作用时,对应的6个应力分量,计算线弹性情况下旋转轴天平上某个指定节点在任意载荷下的MISES应力及应力分量;
预警模块,所述预警模块根据节点的载荷系数得到实时状态下该点的实时应力值,将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时,则报警,对旋转轴天平进行载荷预警。
一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行,以用于实现一种旋转轴天平实时动态载荷监控系数的监测方法。
一种计算机设备,包括存储器和处理器,所述存储器中存储有计算机程序,当所述处理器运行所述存储器存储的计算机程序时,所述处理器执行一种旋转轴天平实时动态载荷监控系数的监测方法。
本发明具有以下有益效果:
本发明在各种不同的载荷条件下,使用基于应力分量的载荷系数计算方法得出的MISES应力(监测应力值)都与有限元计算得出的MISES应力值相等,这种方法可以实时监控旋转轴天平体的应力情况。本发明使用基于应力分量的载荷系数计算方法可以对旋转轴天平进行实时载荷监控,判断天平实时应力状态,极大减小了旋转轴天平在试验过程中发生事故的机率。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为单独受载荷A作用下应力状态图;
图2为单独受载荷B作用下应力状态图;
图3为受载荷A与B同时作用下的应力状态图;
图4为梁根部取点图;
图5为外环与梁连接部分取点图;
图6为梁中部取点图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
此外,下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
以下结合具体实施例,对本发明进行了详细说明。
具体实施例一:
根据图1至图6所示,本发明为解决上述技术问题采取的具体优化技术方案是:本发明涉及一种旋转轴天平实时动态载荷监控系数的监测方法。
一种旋转轴天平实时动态载荷监控系数的监测方法,所述方法包括以下步骤:
步骤1:根据旋转轴天平的应力状态,确定载荷作用时的应力分量和MISES应力;
步骤2:根据旋转天平的应力分量,测量旋转天平中心的受载,确定载荷点所受的应力状态;
步骤3:根据步骤2确定六种基元载荷分别单独作用时,对应的6个应力分量,计算线弹性情况下旋转轴天平上某个指定节点在任意载荷下的MISES应力及应力分量;
步骤4:根据节点的载荷系数得到实时状态下该点的实时应力值,将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时,则报警,对旋转轴天平进行载荷预警。
具体实施例二:
本申请实施例二与实施例一的区别仅在于:
所述步骤1具体为:
根据叠加原理可知,确定载荷作用时的应力分量和MISES应力,应力分量和MISES应力通过下式确定:
其中,为x方向应力分量,、为不同载荷单独作用下的x方向应力分量,为y方向应力分量,、为不同载荷单独作用下y方向应力分量,为剪切应力分量,、为不同载荷单独作用下剪切应力分量;通过、、、、、来计算载荷作用时的应力分量和MISES应力。
具体实施例三:
本申请实施例三与实施例二的区别仅在于:
所述步骤2具体为:
测出天平中心受载通过下式表示:
其中,x为天平受轴向力,y为天平受升力,z为天平受侧向力,为x方向受力系数,为y方向受力系数,为z方向受力系数,L为天平受滚转力矩,为滚转方向载荷系数,M为天平受偏航力矩,为偏航方向载荷系数,N'为天平受俯仰力矩,为俯仰方向载荷系数,N·m为扭矩单位牛米,N为单位牛顿;
其中,为受轴向力x方向产生的应力分量,为受升力x方向产生的应力分量,为受侧向力x方向产生的应力分量,为受滚转力矩x方向产生的应力分量,为受偏航力矩x方向产生的应力分量,为受俯仰力矩x方向产生的应力分量,、、、、、分别为受六元载荷下y方向产生的应力分量,、、、、、分别为受六元载荷下z方向产生的应力分量,、、、、、分别为受六元载荷下xy方向的剪切应力分量,、、、、、分别为六元载荷下yz方向的剪切应力分量,、、、、、为分别受六元载荷下xz方向的剪切应力分量,为天平受到的实际应力。
具体实施例四:
本申请实施例四与实施例三的区别仅在于:
所述步骤3具体为:
此时的应力分量为:
同理可得:
具体实施例五:
本申请实施例五与实施例四的区别仅在于:
具体实施例六:
本申请实施例六与实施例五的区别仅在于:
所述步骤4具体为:
旋转轴天平在使用时,通过设计载荷得到各个需要注意位置节点的应力最大值即应力阈值,通过加载基元载荷,根据旋转轴天平的实时输出,根据天平所受的基元载荷,,,,,,通过得到节点的载荷系数,实时状态下该点的实时应力值,将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时报警,即对旋转轴天平进行载荷预警。
具体实施例七:
本申请实施例七与实施例六的区别仅在于:
当计算X=2000N时某点的应力,只需X=1N时6个应力分量乘以2000即得X=2000N时的六个应力分量;
当计算Y=200N,L=10N*M时某点MISES应力,只需Y=1N时六个应力分量乘以200,对应加上L=1N*M时的六个应力分量乘以10的值,就得到最终的应力分量,最后计算得到MISES应力值。
具体实施例八:
本申请实施例八与实施例七的区别仅在于:
本发明提供一种旋转轴天平实时动态载荷监控系数的监测系统,所述系统包括:
应力确定模块,所述应力确定模块根据旋转轴天平的应力状态,确定载荷作用时的应力分量和MISES应力;
应力状态确定模块,所述应力状态确定模块根据旋转天平的应力分量,测量旋转天平中心的受载,确定载荷点所受的应力状态;
天平应力分量确定模块,所述天平应力分量确定模块确定六种基元载荷分别单独作用时,对应的6个应力分量,计算线弹性情况下旋转轴天平上某个指定节点在任意载荷下的MISES应力及应力分量;
预警模块,所述预警模块根据节点的载荷系数得到实时状态下该点的实时应力值,将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时,则报警,对旋转轴天平进行载荷预警。
具体实施例九:
本申请实施例九与实施例八的区别仅在于:
本发明提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行,以用于实现一种旋转轴天平实时动态载荷监控系数的监测方法。
具体实施例十:
本申请实施例十与实施例九的区别仅在于:
本发明提供一种计算机设备,包括存储器和处理器,所述存储器中存储有计算机程序,当所述处理器运行所述存储器存储的计算机程序时,所述处理器执行一种旋转轴天平实时动态载荷监控系数的监测方法。
具体实施例十一:
本申请实施例十一与实施例十的区别仅在于:
以某旋转轴天平为例,在做比较分析时所用节点来自天平的三个区域,如图4 -图6所示。节点号N1,N2,N3,N4取自天平梁根部,节点号N5,N6,N7取自梁与外环连接部,节点号N8,N9,N10取自应变梁中部。
以二维平面为例,假设此二维平面单独受到某载荷A的作用,那么其某点W应力状态如图1所示;假设此二维平面单独受到某载荷B的作用,那么其某点W应力状态如图2所示;最后,假设此二维平面受到某载荷A与B的共同作用,那么其某点W应力状态如图3所示。
因为始终处于线弹性,所以根据叠加原理可知,有:
现测出天平中心受载为:
只要得到六种基元载荷分别单独作用时,对应的6个应力分量(6×6=36,一个点36个值),就可以准确计算线弹性情况下旋转轴天平上某个指定节点在任意载荷下的MISES应力及应力分量。
这种方法确保了监测应力值替代MISES应力来评估旋转轴天平安全性能的可靠性和安全性。
需注意的是,一个点虽有36个系数,但是并不是在计算时都会用到。比如计算X=2000N时某点的应力,只需X=1N时6个应力分量乘以2000即得X=2000N时的六个应力分量。又比如计算Y=200N,L=10N*M时某点MISES应力,只需Y=1N时六个应力分量乘以200,对应加上L=1N*M时的六个应力分量乘以10的值,就可以得到最终的应力分量,最后计算得到MISES应力值。
因此,对于一些特殊情况,基于应力分量的载荷系数计算方法是需要简化的。在6种载荷,,,,,作用下,应力分量及MISES应力如式(1)-(10)所示。假设,6种载荷等比例扩大到K倍,即6种载荷为,,,,,。
此时的应力分量为:
同理可得:
旋转轴天平在使用时,可以通过天平的设计载荷得到各个需要注意位置节点的应力最大值即应力阈值。通过加载基元载荷,可以得到,须注意位置节点的载荷系数。在试验过程中,根据旋转轴天平的实时输出,通过天平公式可以算出天平所受的六元力,,,,,,通过需要注意位置的节点的载荷系数可以得到,实时状态下该点的实时应力值。将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时,就会报警,即可对旋转轴天平进行载荷预警。
将对两种方法计算得来的监测应力值与有限元计算得来的MISES应力值作对比,以验证其理论正确性。
在对应力分量载荷进行监测时,对天平依次单独施加基元载荷(共6种),提取36个系数(应力分量)。格式如表1所示。
表1 各节点应力分量(MPa)系数表
当所选定监测节点的载荷监测系数全部确定后,则考虑进行简单工况下(如:单分量加载)的对比验证。
首先考虑单个载荷加载情况,即,,,,,只有一项存在。若天平单独受力Y向500N的载荷,那么监测应力值应该由单独受Y向1N力的情况下的6个应力分量系数乘以500,再做MISES等效得出,然后将算得的MISES应力(监测应力值)与天平单独受Y向500N力的情况下的MISES应力值(有限元计算直接得出)进行对比,如果两者计算结果一致则说明本项目所提出的全域载荷监测方法在单分量加载情况下是行得通的。
其次可考虑等比例加载情况。考虑多分量载荷同时作用的情况,但是加载为等比例加载。
外载荷:法兰盘上下表面载荷(等效到几何中心点)为X轴方向200N,Y轴方向200N,Z轴方向200N,绕X轴200N·m扭矩,绕Y轴200N·m弯矩,绕Z轴200N·m弯矩。
这里依旧将算得的MISES应力值(监测应力值)与有限元计算得出的MISES应力值进行对比,如果两者计算结果一致则说明本项目所提出的全域载荷监测方法在等比例加载情况下是行得通的。
最后考虑非等比例加载情况,即对天平施加六分量载荷,但各分量载荷大小方向均不一致。例如:在天平校心分别施加X轴方向200N,Y轴方向300N,Z轴方向400N,绕X轴转动100N·m,绕Y轴转动50N·m,绕Z轴转动60N·m的六分量载荷,使用载荷分量方法计算监测应力值并与有限元法直接计算得到的MISES应力值进行对比,两者计算结果一致,这种方法可以很稳定的衡量天平是否超载。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、 “示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或 者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表 述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或 N 个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下, 本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特 征进行结合和组合。 此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性 或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“N个”的含义是至少两个,例如 两个,三个等,除非另有明确具体的限定。 流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个 或更N个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本发明的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本发明的 实施例所属技术领域的技术人员所理解。 在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实 现逻辑功能的可执行指令的定序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设 备而使用。就本说明书而言,"计算机可读介质"可以是任何可以包含、存储、通信、传播 或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下:具有一个或N个布线的电连接部(电子装置),便携式计算机盘盒(磁装置),随机存取存储器(RAM),只读存储器(ROM),可擦除可编辑只读存储器(EPROM 或闪速存储器),光纤装置,以及便携式光盘只读存储器(CDROM)。另外,计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质,因为可以例如通过对纸或其他介质进行光学扫描,接着进 行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序,然后将其存 储在计算机存储器中。 应当理解,本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实 施方式中,N 个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或 固件来实现。如,如果用硬件来实现和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离 散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场 可编程门阵列(FPGA)等。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可 以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。 此外,在本发明各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。
以上所述仅是一种旋转轴天平实时动态载荷监控系数的监测方法的优选实施方式,一种旋转轴天平实时动态载荷监控系数的监测方法的保护范围并不仅局限于上述实施例,凡属于该思路下的技术方案均属于本发明的保护范围。应当指出,对于本领域的技术人员来说,在不脱离本发明原理前提下的若干改进和变化,这些改进和变化也应视为本发明的保护范围。
Claims (7)
1.一种旋转轴天平实时动态载荷监控系数的监测方法,其特征是:所述方法包括以下步骤:
步骤1:根据旋转轴天平的应力状态,确定载荷作用时的应力分量和MISES应力;
步骤2:根据旋转天平的应力分量,测量旋转天平中心的受载,确定载荷点所受的应力状态;
步骤3:根据步骤2确定六种基元载荷分别单独作用时,对应的6个应力分量,计算线弹性情况下旋转轴天平上某个指定节点在任意载荷下的MISES应力及应力分量;
步骤4:根据节点的载荷系数得到实时状态下该点的实时应力值,将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时,则报警,对旋转轴天平进行载荷预警;
所述步骤2具体为:
测出天平中心受载通过下式表示:
x=Fx×1N y=Fy×1N z=Fz×1N
L=Myz×1N·m M=Mxz×1N·m N'=Myx×1N·m
其中,x为天平受轴向力,y为天平受升力,z为天平受侧向力,Fx为x方向受力系数,Fy为y方向受力系数,Fz为z方向受力系数,L为天平受滚转力矩,Myz为滚转方向载荷系数,M为天平受偏航力矩,Mxz为偏航方向载荷系数,N'为天平受俯仰力矩,Myx为俯仰方向载荷系数,N·m为扭矩单位牛米,N为单位牛顿;
确定载荷点所受的应力状态,包括6个应力分量,监测应力值σA及MISES应力,通过下式表示:
σx=σx1Fx+σx2Fy+σx3Fz+σx4Myz+σx5Mxz+σx6Mxy (4)
σy=σy1Fx+σy2Fy+σy3Fz+σy4Myz+σy5Mxz+σy6Mxy (5)
σz=σz1Fx+σz2Fy+σz3Fz+σz4Myz+σz5Mxz+σz6Mxy (6)
τxy=τxy1Fx+τxy2Fy+τxy3Fz+τxy4Myz+τxy5Mxz+τxy6Mxy (7)
τxz=τxz1Fx+τxz2Fy+τxz3Fz+τxz4Myz+τxz5Mxz+τxz6Mxy (8)
τyz=τyz1Fx+τyz2Fy+τyz3Fz+τyz4Myz+τyz5Mxz+τyz6Mxy (9)
其中,σx1为受轴向力x方向产生的应力分量,σx2为受升力x方向产生的应力分量,σx3为受侧向力x方向产生的应力分量,σx4为受滚转力矩x方向产生的应力分量,σx5为受偏航力矩x方向产生的应力分量,σx6为受俯仰力矩x方向产生的应力分量,σy1、σy2、σy3、σy4、σy5、σy6分别为受六元载荷下y方向产生的应力分量,;σz1、σz2、σz3、σz4、σz5、σz6分别为受六元载荷下z方向产生的应力分量,τxy1、τxy2、τxy3、τxy4、τxy5、τxy6分别为受六元载荷下xy方向的剪切应力分量,τyz1、τyz2、τyz3、τyz4、τyz5、τyz6分别为六元载荷下yz方向的剪切应力分量,τxz1、τxz2、τxz3、τxz4、τxz5、τxz6为分别受六元载荷下xz方向的剪切应力分量,σmises为天平受到的实际应力;
所述步骤3具体为:
基于应力分量的载荷系数,在6种载荷Fx,Fy,Fz,Myz,Mxz,Mxy作用下,应力分量及MISES应力如式(1)-(10)所示,6种载荷等比例扩大到K倍,即6种载荷为KFx,KFy,KFz,KMyz,KMxz,KMxy;
此时的应力分量为:
σx'=σx1KFx+σx2KFy+σx3KFz+σx4KMyz+σx5KMxz+σx6KMxy=Kσx (11)
σy'=Kσy (12)
σz'=Kσz (13)
τxy'=Kτxy (14)
τxz'=Kτxz (15)
τyz'=Kτyz (16)
监测应力值σA'及MISES应力通过下式确定:
2.根据权利要求1所述的一种旋转轴天平实时动态载荷监控系数的监测方法,其特征是:所述步骤1具体为:
根据叠加原理可知,确定载荷作用时的应力分量和MISES应力,应力分量和MISES应力通过下式确定:
σx=σ'x+σ″x (1)
σy=σ'y+σ″y (2)
τxy=τ'xy+τ″xy (3)
其中,σx为x方向应力分量,σ'x、σ″x为不同载荷单独作用下的x方向应力分量,σy为y方向应力分量,σ'y、σ″y为不同载荷单独作用下y方向应力分量,τxy为剪切应力分量,τ'xy、τ″xy为不同载荷单独作用下剪切应力分量;
通过σ'x、σ″x、σ'y、σ″y、τ'xy、τ″xy来计算载荷作用时的应力分量和MISES应力。
3.根据权利要求1所述的一种旋转轴天平实时动态载荷监控系数的监测方法,其特征是:
所述步骤4具体为:
旋转轴天平在使用时,通过设计载荷得到各个需要注意位置节点的应力最大值即应力阈值,通过加载基元载荷,根据旋转轴天平的实时输出,根据天平所受的基元载荷Fx,Fy,Fz,Myz,Mxz,Mxy,通过得到节点的载荷系数,实时状态下该点的实时应力值,将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时报警,即对旋转轴天平进行载荷预警。
4.根据权利要求3所述的一种旋转轴天平实时动态载荷监控系数的监测方法,其特征是:
当计算X=2000N时某点的应力,只需X=1N时6个应力分量乘以2000即得X=2000N时的六个应力分量;
当计算Y=200N,L=10N*M时某点MISES应力,只需Y=1N时六个应力分量乘以200,对应加上L=1N*M时的六个应力分量乘以10的值,就得到最终的应力分量,最后计算得到MISES应力值。
5.一种旋转轴天平实时动态载荷监控系数的监测系统,其特征是:用于实现权利要求1-4任意一项权利要求所述的监测方法,所述系统包括:
应力确定模块,所述应力确定模块根据旋转轴天平的应力状态,确定载荷作用时的应力分量和MISES应力;
应力状态确定模块,所述应力状态确定模块根据旋转天平的应力分量,测量旋转天平中心的受载,确定载荷点所受的应力状态;
天平应力分量确定模块,所述天平应力分量确定模块确定六种基元载荷分别单独作用时,对应的6个应力分量,计算线弹性情况下旋转轴天平上某个指定节点在任意载荷下的MISES应力及应力分量;
预警模块,所述预警模块根据节点的载荷系数得到实时状态下该点的实时应力值,将实时应力值与应力阈值进行对比,当实时应力值超过应力阈值时,则报警,对旋转轴天平进行载荷预警。
6.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行,以用于实现如权利要求1-4任意一项权利要求所述的一种旋转轴天平实时动态载荷监控系数的监测方法。
7.一种计算机设备,其特征在于,包括存储器和处理器,所述存储器中存储有计算机程序,当所述处理器运行所述存储器存储的计算机程序时,所述处理器执行根据权利要求1-4任意一项权利要求所述的一种旋转轴天平实时动态载荷监控系数的监测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211098375.7A CN115165298B (zh) | 2022-09-09 | 2022-09-09 | 一种旋转轴天平实时动态载荷监控系数的监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211098375.7A CN115165298B (zh) | 2022-09-09 | 2022-09-09 | 一种旋转轴天平实时动态载荷监控系数的监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115165298A CN115165298A (zh) | 2022-10-11 |
CN115165298B true CN115165298B (zh) | 2022-12-23 |
Family
ID=83482443
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211098375.7A Active CN115165298B (zh) | 2022-09-09 | 2022-09-09 | 一种旋转轴天平实时动态载荷监控系数的监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115165298B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115752999B (zh) * | 2022-11-11 | 2024-02-13 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种低速风洞天平检测加载头 |
CN115931280B (zh) * | 2023-03-09 | 2023-05-09 | 中国空气动力研究与发展中心低速空气动力研究所 | 铰链力矩风洞试验天平动态载荷实时监视预警方法及系统 |
CN115993229B (zh) * | 2023-03-24 | 2023-05-16 | 中国航空工业集团公司哈尔滨空气动力研究所 | 测量飞机起降过程中非定常气动力系数的风洞试验方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101826131A (zh) * | 2009-03-04 | 2010-09-08 | 中国核电工程有限公司 | 一种大批量标准部件在不同载荷下的应力计算评定方法 |
CN104483044A (zh) * | 2014-12-19 | 2015-04-01 | 深圳供电局有限公司 | 一种瓷柱绝缘子应力检测报警装置及方法 |
CN111780921A (zh) * | 2020-08-07 | 2020-10-16 | 中国科学院武汉岩土力学研究所 | 一种光纤光栅三维应力监测传感器的标定方法 |
CN112800633A (zh) * | 2021-04-06 | 2021-05-14 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种多元校准风洞天平数据的处理方法 |
CN114199436A (zh) * | 2021-12-10 | 2022-03-18 | 哈尔滨工程大学 | 一种基于传感器数据的船舶监测结构整体应力状态快速获取方法 |
CN114444324A (zh) * | 2022-02-15 | 2022-05-06 | 贵州贵飞飞机设计研究院有限公司 | 一种通过多区域分别受载叠加后评估结构件上结构应力关注点处应力大小的方法 |
CN115014694A (zh) * | 2021-11-23 | 2022-09-06 | 中国直升机设计研究所 | 一种旋转天平的数据采集处理方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3217492B2 (ja) * | 1992-09-30 | 2001-10-09 | 大阪瓦斯株式会社 | 応力測定装置 |
KR100500999B1 (ko) * | 2003-08-07 | 2005-07-18 | 삼성중공업 주식회사 | 선체응력 감시장치에서의 변형률 성분 분리법 |
CN110435812A (zh) * | 2012-05-30 | 2019-11-12 | 赛创尼克株式会社 | 通过对海洋结构物的实时测量监视的控制方法 |
-
2022
- 2022-09-09 CN CN202211098375.7A patent/CN115165298B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101826131A (zh) * | 2009-03-04 | 2010-09-08 | 中国核电工程有限公司 | 一种大批量标准部件在不同载荷下的应力计算评定方法 |
CN104483044A (zh) * | 2014-12-19 | 2015-04-01 | 深圳供电局有限公司 | 一种瓷柱绝缘子应力检测报警装置及方法 |
CN111780921A (zh) * | 2020-08-07 | 2020-10-16 | 中国科学院武汉岩土力学研究所 | 一种光纤光栅三维应力监测传感器的标定方法 |
CN112800633A (zh) * | 2021-04-06 | 2021-05-14 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种多元校准风洞天平数据的处理方法 |
CN115014694A (zh) * | 2021-11-23 | 2022-09-06 | 中国直升机设计研究所 | 一种旋转天平的数据采集处理方法 |
CN114199436A (zh) * | 2021-12-10 | 2022-03-18 | 哈尔滨工程大学 | 一种基于传感器数据的船舶监测结构整体应力状态快速获取方法 |
CN114444324A (zh) * | 2022-02-15 | 2022-05-06 | 贵州贵飞飞机设计研究院有限公司 | 一种通过多区域分别受载叠加后评估结构件上结构应力关注点处应力大小的方法 |
Non-Patent Citations (1)
Title |
---|
提高风洞实验数据精度的方法;恽起麟 等;《气动实验与测量控制》;19940930;第8卷(第3期);正文第81页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115165298A (zh) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115165298B (zh) | 一种旋转轴天平实时动态载荷监控系数的监测方法 | |
CN101051055B (zh) | 倾斜角运算方法以及倾斜角运算装置 | |
JP7184698B2 (ja) | 力覚センサ | |
US20120016623A1 (en) | Use of Multiple Internal Sensors for Measurement Validation | |
CN104884926A (zh) | 轴承装置振动分析方法、轴承装置振动分析装置和滚动轴承状态监视装置 | |
CN106644499B (zh) | 电动汽车动力总成强度测试方法和系统 | |
CN110220722A (zh) | 一种控制臂载荷测试系统和标定实验系统 | |
US11391617B2 (en) | Weight-bearing measurement device and method and weight-bearing equipment | |
CN110807285B (zh) | 一种工程运输车车架疲劳分析方法、系统及装置 | |
JP2013032916A (ja) | 多軸力検出器 | |
CN111985134B (zh) | 一种起重机回转支承状态确定方法、装置及相关设备 | |
CN110807283B (zh) | 基于有限元计算模型的多载荷应变测量点获取方法及系统 | |
Gobbi et al. | Race motorcycle smart wheel | |
CN110807284A (zh) | 基于有限元和应变测量的载荷谱识别方法、系统及装置 | |
CN114132324B (zh) | 整车质量估算方法、装置、设备及存储介质 | |
CN113656994B (zh) | 一种汽车悬置系统的悬置力获取方法及装置 | |
dos Santos et al. | The use of dynamic strain sensors and measurements on the ground vibration testing of an f-16 aircraft | |
Gobbi et al. | 6-Axis measuring wheels for trucks or heavy vehicles | |
CN110686631A (zh) | 一种测量t形截面钢压杆初始弯扭缺陷的方法 | |
CN112525310A (zh) | 轮荷测试系统,方法及装置 | |
JP3256709B2 (ja) | 車輪作用力測定装置 | |
EP4089551A1 (en) | Computer implemented method of determining a transfer function of a module or a component and generating such component | |
Yadav et al. | Design and Fabrication of low capacity Torque Transducer | |
US20240192081A1 (en) | Transducer geometry for amplification of axial load response and suppression of non-axial response | |
CN118549083A (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 |