CN107844618B - 用于测量推力和冲量的扭摆系统的设计方法 - Google Patents

用于测量推力和冲量的扭摆系统的设计方法 Download PDF

Info

Publication number
CN107844618B
CN107844618B CN201611048882.4A CN201611048882A CN107844618B CN 107844618 B CN107844618 B CN 107844618B CN 201611048882 A CN201611048882 A CN 201611048882A CN 107844618 B CN107844618 B CN 107844618B
Authority
CN
China
Prior art keywords
measurement
design
measuring
thrust
torsional pendulum
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
Application number
CN201611048882.4A
Other languages
English (en)
Other versions
CN107844618A (zh
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.)
Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Original Assignee
Peoples Liberation Army Strategic Support Force Aerospace Engineering University
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 Peoples Liberation Army Strategic Support Force Aerospace Engineering University filed Critical Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Priority to CN201611048882.4A priority Critical patent/CN107844618B/zh
Publication of CN107844618A publication Critical patent/CN107844618A/zh
Application granted granted Critical
Publication of CN107844618B publication Critical patent/CN107844618B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

本发明提供了一种用于测量推力和冲量的扭摆系统的设计方法。主要包括如下三个步骤:明确扭摆系统的设计条件,关键系统参数设计,实现扭摆系统的设计方案。首先明确测量系统的设计条件,包括推力和冲量的测量量程与误差、推力作用时间、测量噪声误差、位移传感器量程与误差等;其次根据关键系统参数的设计方法,确定扭转刚度系数、力臂、测量臂、转动惯量等关键系统参数;最后,以关键系统参数为依据,提出扭摆测量系统的设计方案,实施并实现扭摆测量系统的设计方案。校核关键系统参数,反复迭代,直到满足设计条件要求,扭摆测量系统设计完成。

Description

用于测量推力和冲量的扭摆系统的设计方法
技术领域
本发明属于空间推进技术领域,是一种测量推力和冲量的扭摆装置的实现方法。
背景技术
激光微推进以激光为动力,烧蚀工作介质产生冲量,实现对空间飞行器的姿态和轨道控制。激光微推进的推力和冲量测量对于研究激光微推进技术是至关重要的。激光微推进所产生的推力范围通常在10-6~10-3N量级之间,冲量范围通常在10-7~10-4N·s量级之间,测量要求精度高、响应快、量程大,测量装置需要进行特殊的设计以满足要求。
在推力测量中,如果推力作用时间与测量系统周期比较足够大,就可以根据推力作用下测量系统响应,分析和计算推力与时间关系,这种情况就是推力测量的问题;如果推力作用时间与测量系统周期比较很小,此时关注的是推力产生的冲量大小,就可以根据推力作用下测量系统响应,分析和计算推力的冲量,这种情况就是冲量测量的问题。
扭摆测量系统既可用于推力的测量,又可用于冲量的测量,在两种情况下,扭摆测量系统的设计方法不同。在推力和冲量的扭摆测量系统设计过程中,迫切需要推力和冲量的扭摆测量系统设计方法,以保证设计质量和实现设计指标要求。
推力和冲量测量的主要区别在于:推力测量采用持续力作用下系统响应,分析和计算推力;冲量测量采用脉冲力作用下系统响应,分析和计算冲量。推力测量和冲量测量的两种情况下,尽管扭摆测量系统的设计方法不同,但是扭摆测量系统的组成和功能相似。
测量系统的结构布局如图2所示。扭摆测量系统组成如下:
(1)扭摆横梁、扭摆转轴、枢轴、配重、阻尼器;
(2)位移传感器、标定力施力装置;
(3)水平测量平面的调平装置;
(4)真空舱的舱外控制装置。
发明内容
本发明提供了一种用于测量推力和冲量的扭摆系统的设计方法。主要包括如下三个步骤:明确扭摆系统的设计条件,关键系统参数设计,实现扭摆系统的设计方案。首先明确测量系统的设计条件,包括推力量程与误差、推力作用时间、测量噪声误差、位移传感器量程与误差等;其次根据关键系统参数的设计方法,确定扭转刚度系数、力臂、测量臂、转动惯量等关键系统参数;最后,以关键系统参数为依据,提出扭摆测量系统的设计方案,实施并实现扭摆测量系统的设计方案。校核关键系统参数,反复迭代,直到满足设计条件要求,扭摆测量系统设计完成。
本发明提供了一种用于测量推力和冲量的扭摆系统的设计方法,其特征在于包括如下步骤:
(1)明确扭摆系统的设计条件
推力量程为[fd,fu](其中,fd<fu),推力测量的相对误差为εf,冲量测量量程为[Sd,Su](其中,Sd<Su),冲量测量的相对误差为εS,在该量程范围内,测量误差满足指定要求;推力作用时间为T0,为推力器指定的工作时间;测量噪声的误差限为Δθ,环境的外力激励和位移激励、位移传感器、推力加载等干扰,将造成测量噪声,采用隔振平台和隔振物体抑制测量噪声。位移传感器的量程为[hsd,hsu](其中,hsd<hsu),在该量程范围内,位移测量的误差限为Δhs
(2)关键系统参数设计
首先,根据公式
Figure GDA0002743381010000031
Cs=5~10,确定测量臂Ls,式中Cs越大对位移传感器误差限要求越严格,一般可取Cs=5。其次,根据公式
Figure GDA0002743381010000032
确定力臂Lf与扭转刚度系数k,式中当阻尼比ζ无法确定时,可取Cf=2(此时,ζ=0)。最后,根据公式
Figure GDA0002743381010000033
确定转动惯量J,式中当阻尼比ζ无法确定时,可取Cζ=1(此时,ζ=0)。
(3)实现扭摆系统的设计方案
确定关键系统参数后,以关键系统参数为依据,进一步确定扭摆横梁、扭摆转轴、枢轴、配重、阻尼器、位移传感器、标定力施力装置、水平测量平面的调平装置、真空舱的舱外控制装置等,提出扭摆测量系统的设计方案。
实施并实现扭摆测量系统的设计方案。通过调节扭转刚度系数与转动惯量乘积,将阻尼比控制在0<ζ≤0.4范围内。
实现扭摆测量系统的设计方案后,采用阶跃标定力确定阻尼比,核实阻尼比为0<ζ≤0.4范围内的某个值。在阻尼比已知条件下,根据已知设计条件,采用步骤(2)方法,由公式
Figure GDA0002743381010000034
Cζ=1-ζ2,重新校核关键系统参数,反复迭代,直到满足已知设计条件要求,扭摆测量系统的设计完成。
本发明的扭摆测量系统设计方法的特点为:
(1)提出了测量系统的设计条件。提出了以推力量程与误差、推力作用时间、测量噪声误差、位移传感器量程与误差等作为初始设计条件的方法,这些设计条件便于工程中提炼和量化。解决了初始设计条件不合理,造成的设计条件过少不能完成参数设计、设计条件过多参数设计冲突的问题。例如,推力作用时间条件,为转动惯量设计提供了依据。
(2)提出了关键系统参数设计方法。紧紧围绕设计条件,通过研究推力与系统响应之间关系,建立了推力量程、系统响应量程、位移传感器量程之间的不等式关系,建立了推力误差、测量噪声误差、传感器误差之间的不等式关系,揭示了扭转刚度系数、力臂、测量臂、转动惯量等关键系统参数之间内在制约关系,提出了关键系统参数设计方法。解决了以满量程的最大误差表示推力误差,造成在量程低端相对误差大的问题,使得推力量程范围内都满足相对误差要求。
(3)提出了扭摆测量系统设计方法。首先,采用提出的关键系统参数设计方法,确定扭转刚度系数、力臂、测量臂、转动惯量等关键系统参数;其次,以关键系统参数为依据,提出扭摆测量系统的设计方案;最后,实施并实现扭摆测量系统的设计方案,校核关键系统参数,反复迭代,直到满足设计条件要求。
附图说明
图1为本发明的用于测量推力和冲量的扭摆系统的设计流程图;
图2为扭摆测量系统的结构布局图;
图3为位移传感器测量扭转角示意图;
图4为系数Kθ随着扭转角θ的变化图。
具体实施方式
现结合附图和实例对本发明扭摆系统的实现方法做进一步详细描述。1.用于测量推力的扭摆系统设计方法
推力的扭摆测量系统设计过程,是根据推力量程与相对误差、推力作用时间、测量噪声、位移传感器量程与测量误差等已知设计条件,提出扭转刚度系数、力臂、测量臂、横梁转动惯量等关键系统参数的设计方案,实施并实现系统设计,最后校核是否达到推力量程与相对误差的设计要求。
1.1已知设计条件
(1)测量系统的推力量程、测量误差
推力量程为[fd,fu](fd<fu),推力测量的相对误差为εf,在该量程范围内,推力的测量误差满足指定要求。
(2)推力作用时间
推力作用时间为T0,为推力器指定的工作时间。
(3)测量噪声
测量噪声的误差限为Δθ,环境的外力激励和位移激励、位移传感器、推力加载等干扰,将造成测量噪声,采用隔振平台和隔振物体抑制测量噪声。
(4)位移传感器的量程、误差限
位移传感器的量程为[hsd,hsu](hsd<hsu),在该量程范围内,位移测量的误差限为Δhs
1.2系统参数的设计依据
(1)位移传感器测量原理和小扭转角条件
如图3所示,扭摆横梁在推力作用下在水平面内扭转振动,位移传感器的探测方向垂直扭摆横梁对称面,位移传感器的测量臂为Ls,测量位移为hs,扭摆横梁的扭转角为
Figure GDA0002743381010000051
位移传感器的测量误差限为Δhs,由于位移传感器测量误差造成的扭转角测量误差为
Figure GDA0002743381010000061
由于位移传感器测量误差造成的扭转角测量的相对误差为
Figure GDA0002743381010000062
Figure GDA0002743381010000063
如图4所示,为系数Kθ随着扭转角θ的变化,当θ<<5°时0.995<Kθ≤1,所以在小扭转角条件下(θ≤5°),由于位移传感器测量误差造成的扭转角测量的相对误差,可近似为
Figure GDA0002743381010000064
扭转角为
Figure GDA0002743381010000065
工程中为了计算方便,采用近似公式
Figure GDA0002743381010000066
具体计算表明,当扭转角θ≤5°时,相对误差小于0.255%,当扭转角θ≤3°时,相对误差小于0.1%,因此小扭转角条件下可采用近似计算公式。
(2)推力测量原理和最大扭转角
在推力f(t)作用下,理想的系统响应为
Figure GDA0002743381010000067
式中,Lf为力臂,J为横梁转动惯量,ζ为阻尼比,ωn为固有振动频率,
Figure GDA0002743381010000068
为振动频率。实际系统响应为
Figure GDA0002743381010000071
式中,Δθ(t)表示测量噪声。
推力测量原理:根据实际系统响应测量值[ti,Θ(ti)](i=0,1,2,…),分析和计算推力f(ti)(i=0,1,2,…)。扭转角为
Figure GDA0002743381010000072
即可通过调节扭转刚度系数与转动惯量的乘积,使得扭转角满足小扭转角条件。
在恒定力f0作用下,扭摆的实际系统响应为
Figure GDA0002743381010000073
Figure GDA0002743381010000074
式中,
Figure GDA0002743381010000077
为扭转刚度系数,Δθ(t)~N(0,σ2)为环境、位移传感器和推力加载等产生的测量噪声,采用零均值正态分布表示,其误差限可取Δθ=3σ。
极值点对应扭转角为
Figure GDA0002743381010000075
最大扭转角为
Figure GDA0002743381010000076
设测量系统的推力量程为[fd,fu](fd<fu),推力量程对应的最大扭转角为
Figure GDA0002743381010000081
其中
Figure GDA0002743381010000082
当阻尼比为零ζ=0时,Cf=2。
在系统参数设计阶段,阻尼比为未知,令阻尼比为零ζ=0,取Cf=2进行设计,测量系统设计并实现后,根据实际阻尼比进行系统参数校核。
1.3系统参数的设计方法
(1)扭转刚度系数与力臂的设计
设恒定力作用下,推力的相对误差为εf,由于恒定力作用下推力的相对误差与扭转角最大值的相对误差相等,令
Figure GDA0002743381010000083
只要满足
Figure GDA0002743381010000084
式中,n≥1为裕度系数,Δθ为扭转角误差限。
因此,要求扭转刚度系数与力臂满足
Figure GDA0002743381010000085
又由于最大扭转角为
Figure GDA0002743381010000086
应在位移传感器量程之内
Figure GDA0002743381010000087
还要求扭转刚度系数与力臂满足
Figure GDA0002743381010000088
即有
Figure GDA0002743381010000091
为了保证上述不等式成立,要求
Figure GDA0002743381010000092
(2)测量臂的设计
设位移传感器的误差限为Δhs,位移传感器误差造成的扭转角误差为Δθs,在小扭转角条件下(θ≤5°),有
Figure GDA0002743381010000093
要求位移传感器能够辨识环境造成的测量噪声,取
Figure GDA0002743381010000094
式中,可取Cs=5~10。
因此,要求测量臂满足
Figure GDA0002743381010000095
又由于
Figure GDA0002743381010000096
Figure GDA0002743381010000097
还要求测量臂满足
Figure GDA0002743381010000098
又由于位移传感器的量程为[hsd,hsu](hsd<hsu),小扭转角条件(θ≤5°)可写作为
Figure GDA0002743381010000099
还要求测量臂满足
Figure GDA00027433810100000910
即有
Figure GDA0002743381010000101
为了保证上述不等式成立,要求
Figure GDA0002743381010000102
(3)转动惯量的设计
设已知推力的作用时间T0,推力作用时间T0与测量系统周期Td比值为
Figure GDA0002743381010000108
为了能够得到推力随着时间变化关系曲线,要求
Figure GDA0002743381010000109
Figure GDA0002743381010000103
由于
Figure GDA0002743381010000104
Figure GDA0002743381010000105
要求转动惯量满足
Figure GDA0002743381010000106
其中
Cζ=1-ζ2
当阻尼比为零ζ=0时,Cζ=1。
为了标定系统参数方便(在阶跃力作用下,系统响应出现明显突出的若干个极值点),一般要求阻尼比0<ζ≤0.4,因此有0.84<Cζ≤1。由于阻尼比
Figure GDA0002743381010000107
通过调节扭转刚度系数与转动惯量乘积,将阻尼比控制在0<ζ≤0.4范围内。
1.4系统参数的设计与校核
系统参数的设计与校核的具体步骤为:
(1)根据系统参数的设计方法,确定扭转刚度系数、力臂、测量臂、转动惯量等关键系统参数。
首先,根据以下公式确定测量臂
Figure GDA0002743381010000111
式中,Cs=5~10,Cs越大对位移传感器误差限要求越严格,一般可取Cs=5。
其次,根据以下公式确定力臂与扭转刚度系数
Figure GDA0002743381010000112
式中,当阻尼比无法确定时,可取Cf=2(ζ=0)。
最后,根据以下公式确定转动惯量
Figure GDA0002743381010000113
式中,当阻尼比无法确定时,可取Cζ=1(ζ=0)。
(2)确定关键系统参数后,以关键系统参数为依据,进一步确定扭摆横梁、扭摆转轴、枢轴、配重、阻尼器、位移传感器、标定力施力装置、水平测量平面的调平装置、真空舱的舱外控制装置等,提出扭摆测量系统的设计方案。
(3)实施并实现扭摆测量系统的设计方案。通过调节扭转刚度系数与转动惯量乘积,将阻尼比控制在0<ζ≤0.4范围内。
(4)校核关键系统参数。
首先,实现扭摆测量系统的设计方案后,采用阶跃标定力确定阻尼比,核实阻尼比为0<ζ≤0.4范围内的某个值。
其次,在阻尼比已知条件下,根据已知设计条件,采用步骤2方法,由
Figure GDA0002743381010000114
重新校核关键系统参数,反复迭代,直到满足已知设计条件要求。最后,根据校核后的关键系统参数,对推力测量系统设计进行调整,完成推力测量系统的设计与实现。
2.冲量的扭摆测量系统设计方法
冲量的扭摆测量系统设计过程,是根据冲量量程与相对误差、测量噪声、位移传感器量程与测量误差等已知设计条件,提出扭转刚度系数、力臂、测量臂、横梁转动惯量等关键系统参数的设计方案,实施并实现系统设计,最后校核是否达到冲量量程与相对误差的设计要求。
2.1已知设计条件
(1)测量系统的冲量量程、测量误差
冲量量程为[Sd,Su](Sd<Su),冲量测量的相对误差为εs,在该量程范围内,冲量的测量误差满足指定要求。
(2)推力作用时间
推力作用时间T0与测量系统周期Td比值为
Figure GDA0002743381010000121
为了满足冲量瞬间作用条件,要求
Figure GDA0002743381010000122
3测量噪声
测量噪声的误差限为Δθ,环境的外力激励和位移激励、位移传感器、冲量加载等干扰,将造成测量噪声,采用隔振平台和隔振物体抑制测量噪声。
4位移传感器的量程、误差限
位移传感器的量程为[hsd,hsu](hsd<hsu),在该量程范围内,位移测量的误差限为0<ζ≤0.3。
2.2系统参数的设计依据
(1)位移传感器测量原理和小扭转角条件
与推力测量情况相同,见1.2中相关内容。
2冲量测量原理和最大扭转角
在推力f(t)作用下,理想的系统响应为
Figure GDA0002743381010000136
式中,Lf为力臂,J为横梁转动惯量,ζ为阻尼比,ωn为固有振动频率,
Figure GDA0002743381010000131
为振动频率。
脉冲力是指作用时间很短的瞬间作用推力,冲量为S的脉冲力为f(τ)=Sδ(τ),此时,理想系统响应为
Figure GDA0002743381010000132
实际系统响应为
Figure GDA0002743381010000133
式中,Δθ(t)~N(0,σ2)为环境、位移传感器和冲量加载等产生的测量噪声,采用零均值正态分布表示,其误差限可取Δθ=3σ。
冲量测量原理:根据实际系统响应测量值[ti,Θ(ti)](i=0,1,2,…),分析和计算冲量S。
扭转角为
Figure GDA0002743381010000134
即可通过调节扭转刚度系数与转动惯量的乘积,使得扭转角满足小扭转角条件。
在瞬间作用冲量S作用下,最大扭转角为
Figure GDA0002743381010000135
式中,ωdtM1=π/2。
设测量系统的冲量量程为[Sd,Su](Sd<Su),冲量量程对应的最大扭转角为
Figure GDA0002743381010000141
其中
Figure GDA0002743381010000142
当阻尼比为零ζ=0时,CS=1。
在系统参数设计阶段,阻尼比为未知,令阻尼比为零ζ=0,取CS=1进行设计,测量系统设计并实现后,根据实际阻尼比进行系统参数校核。
2.3系统参数的设计方法
(1)扭转刚度系数、转动惯量、力臂的设计
设冲量测量的相对误差为εS,由于冲量测量的相对误差与扭转角测量的相对误差相等,令
Figure GDA0002743381010000143
只要
Figure GDA0002743381010000144
式中,n≥1为裕度系数。
因此,要求扭转刚度系数、转动惯量、力臂满足
Figure GDA0002743381010000145
又由于最大扭转角为
Figure GDA0002743381010000146
应在位移传感器量程之内
Figure GDA0002743381010000147
还要求扭转刚度系数、转动惯量、力臂满足
Figure GDA0002743381010000148
即有
Figure GDA0002743381010000151
设已知冲量的作用时间T0,冲量作用时间T0与测量系统周期Td比值为
Figure GDA0002743381010000159
为了满足冲量瞬间作用条件,要求
Figure GDA00027433810100001510
具体计算表明当阻尼比0<ζ≤0.3时,只要
Figure GDA00027433810100001511
冲量瞬间作用模型误差不大于1%。因此,设
Figure GDA0002743381010000152
由于
Figure GDA0002743381010000153
Figure GDA0002743381010000154
要求扭转刚度系数、转动惯量满足
Figure GDA0002743381010000155
其中
Cζ=1-ζ2
当阻尼比为零ζ=0时,Cζ=1。
为了标定系统参数方便(在阶跃力作用下,系统响应出现明显突出的若干个极值点),以及冲量瞬间作用下模型误差不大于1%的条件,一般要求阻尼比0<ζ≤0.3,因此有0.91≤Cζ<1。由于阻尼比
Figure GDA0002743381010000156
通过调节扭转刚度系数与转动惯量乘积,将阻尼比控制在0<ζ≤0.3范围内。
因此,要求转动惯量满足、扭转刚度系数、力臂满足
Figure GDA0002743381010000157
为了保证上述不等式成立,要求
Figure GDA0002743381010000158
2测量臂的设计
设位移传感器的误差限为Δhs,位移传感器误差造成的扭转角误差为ΔθS,在小扭转角条件下(θ≤5°),有
Figure GDA0002743381010000161
要求位移传感器能够辨识环境造成的测量噪声,取
Figure GDA0002743381010000162
式中,可取Cs=5~10。
因此,要求测量臂满足
Figure GDA0002743381010000163
又由于
Figure GDA0002743381010000168
Figure GDA0002743381010000169
还要求测量臂满足
Figure GDA0002743381010000164
又由于位移传感器的量程为[hsd,hsu](hsd<hsu),小扭转角条件(θ≤5°)可写作为
Figure GDA0002743381010000165
还要求测量臂满足
Figure GDA00027433810100001610
即有
Figure GDA0002743381010000166
为了保证上述不等式成立,要求
Figure GDA0002743381010000167
2.4系统参数的设计与校核
系统参数的设计与校核的具体步骤为:
(1)根据2.3系统参数的设计方法,确定扭转刚度系数、力臂、测量臂、转动惯量等关键系统参数。
首先,根据以下公式确定测量臂
Figure GDA0002743381010000171
式中,Cs=5~10,Cs越大对位移传感器误差限要求越严格,一般可取Cs=5。
其次,在枢轴货架产品中,选择枢轴类型,确定扭转刚度系数的选择范围[kd,ku](kd<ku)。
最后,选择不同力臂,根据以下公式,通过反复迭代计算,确定扭转刚度系数、力臂、转动惯量
Figure GDA0002743381010000172
式中,当阻尼比为零ζ=0时,Cζ=1和CS=1。
(2)确定关键系统参数后,以关键系统参数为依据,进一步确定扭摆横梁、扭摆转轴、枢轴、配重、阻尼器、位移传感器、标定力施力装置、水平测量平面的调平装置、真空舱的舱外控制装置等,提出扭摆测量系统的设计方案。
(3)实施并实现扭摆测量系统的设计方案。通过调节扭转刚度系数与转动惯量乘积,将阻尼比控制在0<ζ≤0.3范围内。
(4)校核关键系统参数。
首先,实现扭摆测量系统的设计方案后,采用阶跃标定力确定阻尼比,核实阻尼比为0<ζ≤0.3范围内的某个值。
其次,在阻尼比已知条件下,根据已知设计条件,采用步骤2.2方法,由
Figure GDA0002743381010000181
重新校核关键系统参数,反复迭代,直到满足已知设计条件要求。
最后,根据校核后的关键系统参数,对冲量测量系统设计进行调整,完成冲量测量系统的设计与实现。

Claims (1)

1.一种用于测量推力和冲量的扭摆系统的设计方法,其特征在于包括如下步骤:
(1)明确扭摆系统的设计条件
推力测量量程为[fd,fu],fd<fu且单位为N,推力测量的相对误差为εf,冲量测量量程为[Sd,Su],Sd<Su且单位为N·s,冲量测量的相对误差为εS,在该量程范围内,测量误差满足指定要求;推力作用时间为T0,为推力器指定的工作时间;测量噪声的误差限为Δθ,环境的外力激励和位移激励、位移传感器、推力加载将造成测量噪声,采用隔振平台和隔振物体抑制测量噪声;位移传感器的量程为[hsd,hsu],hsd<hsu且单位为mm,在该量程范围内,位移测量的误差限为Δhs
(2)关键系统参数设计
首先,根据公式
Figure FDA0002768056990000011
确定测量臂Ls,式中,Cs为测量噪声误差限的放大系数,εf为裕度系数;其次,根据公式
Figure FDA0002768056990000012
确定力臂Lf与扭转刚度系数k,其中,
Figure FDA0002768056990000013
当阻尼比ζ无法确定时,取ζ=0、Cf=2;最后,根据公式
Figure FDA0002768056990000014
确定转动惯量J,式中当阻尼比ζ无法确定时,取ζ=0、Cζ=1;
(3)实现扭摆系统的设计方案
确定关键系统参数后,以关键系统参数为依据,进一步确定扭摆横梁、扭摆转轴、枢轴、配重、阻尼器、位移传感器、标定力施力装置、水平测量平面的调平装置、真空舱的舱外控制装置,提出扭摆测量系统的设计方案;
实施并实现扭摆测量系统的设计方案;通过调节扭转刚度系数与转动惯量乘积,将推力测量的阻尼比控制在0<ζ≤0.4范围内,将冲量测量的阻尼比控制在0<ζ≤0.3范围内;
实现扭摆测量系统的设计方案后,采用阶跃标定力确定阻尼比,核实阻尼比为0<ζ≤0.4或0<ζ≤0.3范围内的某个值;在阻尼比已知条件下,根据已知设计条件,采用步骤(2)方法,由公式
Figure FDA0002768056990000021
Cζ=1-ζ2,重新校核关键系统参数,反复迭代,直到满足已知设计条件要求,扭摆测量系统的设计完成。
CN201611048882.4A 2016-11-22 2016-11-22 用于测量推力和冲量的扭摆系统的设计方法 Active CN107844618B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611048882.4A CN107844618B (zh) 2016-11-22 2016-11-22 用于测量推力和冲量的扭摆系统的设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611048882.4A CN107844618B (zh) 2016-11-22 2016-11-22 用于测量推力和冲量的扭摆系统的设计方法

Publications (2)

Publication Number Publication Date
CN107844618A CN107844618A (zh) 2018-03-27
CN107844618B true CN107844618B (zh) 2021-01-01

Family

ID=61661046

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611048882.4A Active CN107844618B (zh) 2016-11-22 2016-11-22 用于测量推力和冲量的扭摆系统的设计方法

Country Status (1)

Country Link
CN (1) CN107844618B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108680303B (zh) * 2018-07-06 2019-07-30 中国人民解放军战略支援部队航天工程大学 用于扭摆测力的系统响应误差计算方法、扭摆测力系统
CN109813496A (zh) * 2019-03-26 2019-05-28 北京卫星环境工程研究所 平面气浮式转动惯量测量台校准方法
CN112284589B (zh) * 2020-09-02 2021-12-07 上海新力动力设备研究所 一种对称钟摆式微推力测量装置
CN113447180B (zh) * 2021-05-10 2022-03-22 中山大学 一种基于差分电容的扭摆测试装置、扭摆解耦读取方法
CN113392482B (zh) * 2021-06-29 2023-03-10 中国科学院长春光学精密机械与物理研究所 变扭转刚度装置及扭转刚度修正方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6032904A (en) * 1998-02-23 2000-03-07 Space Systems/Loral, Inc. Multiple usage thruster mounting configuration
CN102175376A (zh) * 2011-01-27 2011-09-07 哈尔滨工业大学 多光束激光外差测量微冲量的装置及方法
CN102252794A (zh) * 2011-05-31 2011-11-23 哈尔滨工业大学 基于多光束激光外差法和扭摆法测量微冲量的方法
CN105784237A (zh) * 2016-05-13 2016-07-20 中国科学院力学研究所 一种微推力测试系统及方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6032904A (en) * 1998-02-23 2000-03-07 Space Systems/Loral, Inc. Multiple usage thruster mounting configuration
CN102175376A (zh) * 2011-01-27 2011-09-07 哈尔滨工业大学 多光束激光外差测量微冲量的装置及方法
CN102252794A (zh) * 2011-05-31 2011-11-23 哈尔滨工业大学 基于多光束激光外差法和扭摆法测量微冲量的方法
CN105784237A (zh) * 2016-05-13 2016-07-20 中国科学院力学研究所 一种微推力测试系统及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《一种用于微小推力冲量测量的扭摆系统参数标定方法》;金星 等;《推进技术》;20151030;第36卷(第10期);全文 *

Also Published As

Publication number Publication date
CN107844618A (zh) 2018-03-27

Similar Documents

Publication Publication Date Title
CN107844618B (zh) 用于测量推力和冲量的扭摆系统的设计方法
Cobleigh et al. Flush airdata sensing (FADS) system calibration procedures and results for blunt forebodies
Hamel et al. Advances in rotorcraft system identification
US5214596A (en) System for determining the airspeed of helicopters
US9909954B2 (en) System and method to quantify viscous damping steering feel of a vehicle equipped with an electric power steering system
Wang et al. Force measurement using strain-gauge balance in a shock tunnel with long test duration
CN111220347B (zh) 一种飞行器气动协调修正方法
US8332081B2 (en) Methods and systems for reducing the phenomenon of structural coupling in the control system of an in-flight refuelling boom
JP2018020759A (ja) 航空機の空力角度のための適合フィルタリングシステム
US6711476B2 (en) Method and computer program product for estimating at least one state of a dynamic system
CN113848963B (zh) 一种飞行控制系统的控制律参数设计方法
CN108469734A (zh) 考虑状态约束的电机伺服系统自抗扰控制方法
CN105527455B (zh) 用于自动地估计飞行器的至少一个速度的方法及装置
Roizner et al. Safe flutter tests using parametric flutter margins
CN112068444B (zh) 一种采用非线性自适应滑模的飞行器攻角控制方法
CN109960831B (zh) 用于扭摆系统的微推力平滑降噪优化还原方法
CN115793449B (zh) 直升机飞行姿态控制器设计方法及装置、设备、存储介质
CN107843389B (zh) 用于扭摆系统的冲量测量误差评定方法
CN112432746B (zh) 一种确定直升机桨叶挥摆载荷耦合程度的方法
CN112124620A (zh) 应用于飞行器的伺服负载力矩计算方法、设备及存储介质
CN113325709B (zh) 转台变负载情况下的自适应控制方法
CN111442786B (zh) 一种飞行器陀螺仪的零点漂移偏差和姿态估计方法
Jain et al. A novel methodology to update table in air data system of a high performance fighter aircraft
Pieniążek et al. Integrated measurement system for UAV
CN111373877B (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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 101416 3380 mailbox No. 86, equipment college, Bayi Road, Huairou District, Beijing, No. 86

Applicant after: China People's Liberation Army Strategic Support Unit Space Engineering University

Address before: 101416 3380 mailbox No. 86, equipment college, Bayi Road, Huairou District, Beijing, No. 86

Applicant before: PLA 'S Equipment College

GR01 Patent grant
GR01 Patent grant