[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN115636063A - High-precision real-time identification method and device for thrust of ship propeller shaft - Google Patents

High-precision real-time identification method and device for thrust of ship propeller shaft Download PDF

Info

Publication number
CN115636063A
CN115636063A CN202211390047.4A CN202211390047A CN115636063A CN 115636063 A CN115636063 A CN 115636063A CN 202211390047 A CN202211390047 A CN 202211390047A CN 115636063 A CN115636063 A CN 115636063A
Authority
CN
China
Prior art keywords
equation
ship
propeller shaft
thrust
strain
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.)
Granted
Application number
CN202211390047.4A
Other languages
Chinese (zh)
Other versions
CN115636063B (en
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.)
Shanghai Jiao Tong University
Original Assignee
Shanghai Jiao Tong 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 Shanghai Jiao Tong University filed Critical Shanghai Jiao Tong University
Priority to CN202211390047.4A priority Critical patent/CN115636063B/en
Publication of CN115636063A publication Critical patent/CN115636063A/en
Application granted granted Critical
Publication of CN115636063B publication Critical patent/CN115636063B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

本发明公开了一种船舶螺旋桨轴推力高精度实时辨识方法与装置,包括如下步骤:加工测量装置并安装在转轴上;获取测量装置纵向应变及环境温度,构建推力、温度与应变的第一等式;对第一等式求导并按时间离散得到第二等式,基于第二等式得到系统矩阵方程与观测方程;基于矩阵方程完成第k个时间步推力状态预测,基于观测方程完成推力状态估计;将时间递推至第k+1个时间步,重复状态预测与状态估计,最终获得推力动态过程。本发明具有较强抗干扰能力,可补偿环境温度变化引起的热变形误差,能在强干扰噪声中保持较高测量精度。此外,本发明中应变片贴在测量装置上,可在实验室中标定应变片灵敏度系数,从而消除应变片粘贴方式、胶水种类等外因影响。

Figure 202211390047

The invention discloses a method and device for high-precision real-time identification of propeller shaft thrust of a ship, comprising the following steps: processing a measuring device and installing it on a rotating shaft; Formula; deriving the first equation and discretizing it according to time to obtain the second equation, based on the second equation, the system matrix equation and observation equation are obtained; based on the matrix equation, the thrust state prediction of the kth time step is completed, and the thrust state is completed based on the observation equation State estimation: recursively advance the time to the k+1th time step, repeat the state prediction and state estimation, and finally obtain the thrust dynamic process. The invention has strong anti-interference ability, can compensate thermal deformation error caused by ambient temperature change, and can maintain high measurement accuracy in strong interference noise. In addition, in the present invention, the strain gauge is pasted on the measuring device, and the sensitivity coefficient of the strain gauge can be calibrated in the laboratory, thereby eliminating the influence of external factors such as the sticking method of the strain gauge and the type of glue.

Figure 202211390047

Description

一种船舶螺旋桨轴推力高精度实时辨识方法及装置A high-precision real-time identification method and device for ship propeller shaft thrust

技术领域technical field

本发明属于船舶动力领域,尤其涉及一种船舶螺旋桨轴推力高精度实时辨识方法及装置。The invention belongs to the field of ship power, in particular to a method and device for high-precision real-time identification of ship propeller shaft thrust.

背景技术Background technique

在线、实时、准确监测船舶推进系统在各个转速下的推力,可为船舶快速性预报、航行状态监测、船-机-桨优化匹配等提供准确的数据需求,具有重要工程意义。然而,当前面临的技术瓶颈主要是在实船干扰环境下难以准确、实时监测螺旋桨推力,主要原因是推力在轴系上产生的微弱形变信号极易被轮机舱的温度热变形、机械噪声、电磁噪声、轴系振动等噪声淹没,难以从测量的形变信号中辨识微弱推力信号。Online, real-time, and accurate monitoring of the thrust of the ship's propulsion system at various speeds can provide accurate data requirements for rapid ship prediction, navigation status monitoring, ship-engine-propeller optimization matching, etc., which is of great engineering significance. However, the current technical bottleneck is that it is difficult to monitor the propeller thrust accurately and in real time under the interference environment of the real ship. Noise such as noise and shafting vibration is submerged, and it is difficult to identify weak thrust signals from the measured deformation signals.

当前已有一些专利提出推力测试相关专利,包括方法类和装置类。At present, some patents have proposed patents related to thrust testing, including methods and devices.

方法类:一种水轮发电机转轴应力状态分析方法。该专利包括如下几个步骤:S1、根据测试转轴,建立模型有限元计算模型;S2、通过理论计算验证有限元模型正确性;S3、在轴上布置应变片,测试转轴轴向等应力;S4、在有限有模型中施加边界条件,使计算应力与测试应力相等;S5、通过分析有限元计算模型的计算结果和现场应力测试结果,建立转轴最大应力值与现场应力测试值之间的关系,实现由现场应力测试结果直接评估转轴寿命和安全状态的目的。Method category: a method for analyzing the stress state of the shaft of a hydroelectric generator. The patent includes the following steps: S1. Establish a finite element calculation model based on the test shaft; S2. Verify the correctness of the finite element model through theoretical calculations; S3. Arrange strain gauges on the shaft to test the axial stress of the shaft; S4 1. Apply boundary conditions in the finite model to make the calculated stress equal to the test stress; S5. By analyzing the calculation results of the finite element calculation model and the on-site stress test results, establish the relationship between the maximum stress value of the rotating shaft and the on-site stress test value, The purpose of directly evaluating the life and safety state of the rotating shaft from the results of the field stress test is realized.

装置类:水下吊舱或舵桨推进装置模型驱动及推力、扭矩测试系统。该专利包括驱动模块、测量模块和传输模块。驱动模块采用直流电机,电机旋转轴连接测量模块,测量模块另一端连接螺旋桨。测量模块上的应变片可测量推力产生的应变。传输模块通过非接触耦合方式传输信号及应变片的输入电压。Device category: underwater pod or rudder propulsion device model drive and thrust and torque test system. The patent includes a drive module, a measurement module and a transmission module. The driving module adopts a DC motor, the rotating shaft of the motor is connected to the measuring module, and the other end of the measuring module is connected to the propeller. Strain gauges on the measurement module measure the strain induced by the thrust. The transmission module transmits the signal and the input voltage of the strain gauge through non-contact coupling.

商用车推力杆轴向力测试装置。该专利包括推力杆、应变片、数据采集仪。其中,推力杆与应变片固定安装,应变片通过粘合层贴合固定安装在平行于推力杆轴线的推力杆表面;应变片与数据采集仪电连接;数据采集仪采集应变片的数据。Commercial vehicle thrust rod axial force test device. The patent includes thrust rods, strain gauges, and data acquisition instruments. Among them, the thrust rod and the strain gauge are fixedly installed, and the strain gauge is fixedly installed on the surface of the thrust rod parallel to the axis of the thrust rod through the adhesive layer; the strain gauge is electrically connected to the data acquisition instrument; the data acquisition instrument collects the data of the strain gauge.

然而,上述专利中均具有抗干扰能力低的主要缺陷,在温度变化、强干扰噪声环境中测试精度低的技术问题。However, the above-mentioned patents all have the main defect of low anti-interference ability, and the technical problem of low test accuracy in the environment of temperature changes and strong interference noise.

发明内容Contents of the invention

本发明的技术目的是提供一种船舶螺旋桨轴推力高精度实时辨识方法及装置,以解决环境温度变化及强干扰噪声造成的测量精度低的技术问题。The technical purpose of the present invention is to provide a high-precision real-time identification method and device for ship propeller shaft thrust to solve the technical problem of low measurement accuracy caused by environmental temperature changes and strong interference noise.

为解决上述问题,本发明的技术方案为:In order to solve the above problems, the technical solution of the present invention is:

一种船舶螺旋桨轴推力高精度实时辨识方法,包括如下步骤:A high-precision real-time identification method for ship propeller shaft thrust, comprising the following steps:

获取安装于船舶螺旋桨轴上的至少2片应变片的应变信号,以及船舶螺旋桨轴的环境温度;Obtain the strain signals of at least two strain gauges installed on the propeller shaft of the ship, and the ambient temperature of the propeller shaft of the ship;

根据应变信号和环境温度构建船舶螺旋桨轴推力与其纵向应变的第一等式;对第一等式进行求导并按时间步离散得到第二等式;基于第二等式表示成矩阵方程和观测方程;Construct the first equation of the propeller shaft thrust and its longitudinal strain according to the strain signal and the ambient temperature; derivate the first equation and discretize it according to time steps to obtain the second equation; based on the second equation, express it as a matrix equation and observation equation;

基于矩阵方程和观测方程完成状态预测,得到第k个时间步相对应的先验估计;基于先验估计完成第k个时间步下的状态估计,得到第k个时间步相对应的估计变量;将时间递推至k+1个时间步,重复状态预测与状态估计,以此得到新的估计变量,基于新的估计变量与第一等式计算得到船舶螺旋桨轴的轴系推力。The state prediction is completed based on the matrix equation and the observation equation, and the prior estimation corresponding to the k-th time step is obtained; the state estimation at the k-th time step is completed based on the prior estimation, and the estimated variable corresponding to the k-th time step is obtained; Recurse the time to k+1 time steps, repeat the state prediction and state estimation, so as to obtain new estimated variables, and calculate the shafting thrust of the ship propeller shaft based on the new estimated variables and the first equation.

其中,根据应变信号和环境温度构建船舶螺旋桨轴推力与其纵向应变的第一等式具体步骤为Among them, according to the strain signal and the ambient temperature, the specific steps for constructing the first equation of the ship propeller shaft thrust and its longitudinal strain are as follows

分别读取不同应变片的应变信号,其数度满足以下公式:Read the strain signals of different strain gauges respectively, and the magnitudes satisfy the following formula:

ε1=εFM+a·(T-T0);ε2=εFM+a·(T-T0)ε 1FM +a·(TT 0 ); ε 2FM +a·(TT 0 )

其中,ε1、ε2分别为应变片R1、R2产生的应变,εF为船舶螺旋桨轴推力所产生的纵向应变,εM为弯矩产生的应变,a为轴系热膨胀系数,T0为环境初始温度,T为环境温度;Among them, ε 1 and ε 2 are the strains produced by the strain gauges R 1 and R 2 respectively, ε F is the longitudinal strain produced by the shaft thrust of the propeller of the ship, ε M is the strain produced by the bending moment, a is the coefficient of thermal expansion of the shaft system, T 0 is the initial temperature of the environment, T is the ambient temperature;

分别与不同应变片连接的应变仪的读数ε的公式为:The formula of the reading ε meter of the strain gauge connected with different strain gauges is:

ε=ε12=2εF+2a·(T-T0) εmeter =ε 12 =2ε F +2a·(TT 0 )

船舶螺旋桨轴推力FT与纵向应变εF关系满足以下第一等式:The relationship between ship propeller shaft thrust F T and longitudinal strain ε F satisfies the following first equation:

Figure BDA0003931644910000031
Figure BDA0003931644910000031

其中,E与A分别为应变片安装处轴段的弹性模量与截面积。Among them, E and A are the elastic modulus and cross-sectional area of the shaft section where the strain gauge is installed, respectively.

其中,对第一等式进行求导并按时间步离散得到第二等式,具体公式为:Among them, deriving the first equation and discretizing it according to time steps to obtain the second equation, the specific formula is:

Figure BDA0003931644910000032
Figure BDA0003931644910000032

其中,ε(k+1)为第k+1个时间步下的应变仪读数,ε(k)为第k个时间步下的应变仪读数,ΔFT(k)为第k个时间步下的推力变化值,ΔT(k)为第k个时间步下环境温度变化值。Among them, εmeter (k+1) is the strain gauge reading at the k+1th time step, εmeter (k) is the strain gauge reading at the kth time step, and ΔF T (k) is the kth time The change value of the thrust at the next step, ΔT(k) is the change value of the ambient temperature at the kth time step.

其中,基于第二等式表示成矩阵方程具体为Among them, based on the second equation, it is expressed as a matrix equation, specifically as

令Xk+1=[ε(k+1)ΔFT(k)ΔT(k+1)]T,则矩阵方程为Let X k+1 =[ε(k+1)ΔF T (k)ΔT(k+1)] T , then the matrix equation is

Xk+1=Φ·Xk+wk X k+1 =Φ·X k +w k

其中,

Figure BDA0003931644910000033
为状态转移矩阵,wk为模拟噪声的模型误差;in,
Figure BDA0003931644910000033
is the state transition matrix, w k is the model error of the simulated noise;

基于第二等式表示成观测方程具体为Based on the second equation, the observation equation is expressed as

yk=H·Xk+vk y k =H·X k +v k

其中,yk为测量的应变与测量的温度的变化,

Figure BDA0003931644910000034
为观测方程,vk为测量噪声。where yk is the change in measured strain versus measured temperature,
Figure BDA0003931644910000034
is the observation equation, and v k is the measurement noise.

其中,在第k个时间步,得到第k个时间步相对应的先验估计具体为:Among them, at the kth time step, the prior estimate corresponding to the kth time step is specifically:

Figure BDA0003931644910000041
Figure BDA0003931644910000041

Figure BDA0003931644910000042
Figure BDA0003931644910000042

其中,

Figure BDA0003931644910000043
为X在第k个时间步的先验估计,
Figure BDA0003931644910000044
为X在第k-1个时间步的后验估计,
Figure BDA0003931644910000045
Figure BDA0003931644910000046
的估计误差协方差矩阵,
Figure BDA0003931644910000047
Figure BDA0003931644910000048
的估计误差协方差矩阵,Q为模型误差wk的协方差矩阵;in,
Figure BDA0003931644910000043
is the prior estimate of X at the kth time step,
Figure BDA0003931644910000044
is the posterior estimate of X at the k-1th time step,
Figure BDA0003931644910000045
for
Figure BDA0003931644910000046
The estimated error covariance matrix of ,
Figure BDA0003931644910000047
for
Figure BDA0003931644910000048
The estimated error covariance matrix of , Q is the covariance matrix of the model error w k ;

在第k个时间步,得到第k个时间步相对应的估计变量具体为:At the kth time step, the estimated variables corresponding to the kth time step are specifically:

Kk=Pk|k-1HT(HPk|k-1HT+R)-1 K k =P k|k-1 H T (HP k|k-1 H T +R) -1

Figure BDA0003931644910000049
Figure BDA0003931644910000049

Figure BDA00039316449100000410
Figure BDA00039316449100000410

Pk|k=(1-KkH)Pk|k-1 P k|k =(1-K k H)P k|k-1

其中,Kk为增益矩阵,Zk为测量信号与先验估计的偏差序列,R为测量噪声vk的协方差矩阵;Among them, K k is the gain matrix, Z k is the deviation sequence between the measurement signal and the prior estimate, and R is the covariance matrix of the measurement noise v k ;

协方差矩阵Q与协方差矩阵R的公式分别为:The formulas of covariance matrix Q and covariance matrix R are respectively:

Figure BDA00039316449100000411
Figure BDA00039316449100000411

R=Ck-HPk|k-1HT R=C k -HP k|k-1 H T

Figure BDA00039316449100000412
Figure BDA00039316449100000412

Figure BDA00039316449100000413
Figure BDA00039316449100000413

进一步优选地,在获取安装于船舶螺旋桨轴上的至少2片应变片的应变信号以及船舶螺旋桨轴的环境温度之前还包括如下步骤Further preferably, the following steps are also included before obtaining the strain signals of at least 2 strain gauges installed on the ship's propeller shaft and the ambient temperature of the ship's propeller shaft

在船舶螺旋桨轴选取测点位置,并测量测点处的轴段直径以获得截面积A;Select the position of the measuring point on the propeller shaft of the ship, and measure the diameter of the shaft section at the measuring point to obtain the cross-sectional area A;

在测点位置布置安装环以粘贴2片应变片,并校核应变片的灵敏度系数;Arrange the mounting ring at the measuring point to paste 2 strain gauges, and check the sensitivity coefficient of the strain gauges;

将两片应变片采用半桥接线方式连接。The two strain gauges are connected by a half-bridge connection.

一种船舶螺旋桨轴推力高精度实时辨识装置,配置有上述的船舶螺旋桨轴推力高精度实时辨识方法,A high-precision real-time identification device for ship propeller shaft thrust, configured with the above-mentioned high-precision real-time identification method for ship propeller shaft thrust,

包括安装组件、测试件;安装组件安装于船舶螺旋桨轴上测试点处;测试件贴设于安装组件上;Including installation components and test pieces; the installation components are installed at the test points on the propeller shaft of the ship; the test pieces are attached to the installation components;

安装组件包括第一安装件和第二安装件;The mounting assembly includes a first mounting part and a second mounting part;

第一安装件与第二安装件结构相同,均包括第一弧形件、第二弧形件和至少一连接杆;The first mounting part and the second mounting part have the same structure, and both include a first arc-shaped part, a second arc-shaped part and at least one connecting rod;

第一弧形件和第二弧形件均为半圆形,连接杆的两端分别与第一弧形件一侧和第二弧形件一侧连接,且第一弧形件与第二弧形件对称设置;至少一测试件安装于连接杆上;Both the first arc and the second arc are semicircular, and the two ends of the connecting rod are respectively connected to one side of the first arc and one side of the second arc, and the first arc and the second The arc parts are arranged symmetrically; at least one test piece is installed on the connecting rod;

第一弧形件和第二弧形件的两端设有连接孔,连接孔的开设方向与连接杆相互垂直,第一安装件与第二安装件通过连接孔与螺栓螺母配合,实现螺栓连接。Both ends of the first arc-shaped piece and the second arc-shaped piece are provided with connecting holes, and the opening direction of the connecting holes is perpendicular to the connecting rod. The first mounting part and the second mounting part cooperate with bolts and nuts through the connecting holes to realize bolt connection .

具体地,第一安装件与第二安装件均设有两个连接杆;Specifically, both the first mounting part and the second mounting part are provided with two connecting rods;

且两个连接杆的两端均分别与第一弧形件一侧和第二弧形件一侧连接,连个连接杆的连接位置将第一弧形件或第二弧形件均分为三段弧形;每个连接杆上均设有至少一测试件。And both ends of the two connecting rods are respectively connected to one side of the first arc-shaped piece and one side of the second arc-shaped piece, and the connecting positions of the two connecting rods divide the first arc-shaped piece or the second arc-shaped piece Three arcs; each connecting rod is provided with at least one test piece.

其中,测试件为应变片,且两个相对设置的应变片电连接,采用半桥接线。Wherein, the test piece is a strain gauge, and two oppositely arranged strain gauges are electrically connected by using a half-bridge connection.

本发明由于采用以上技术方案,使其与现有技术相比具有以下的优点和积极效果:Compared with the prior art, the present invention has the following advantages and positive effects due to the adoption of the above technical scheme:

本发明一方面采用预估-校正的加权递推最小二乘估计从带噪声干扰的应变信号中辨识螺旋桨推力,具有较强的抗干扰能力,能在强干扰噪声环境中仍能保持较高的测量精度。另一方面,递推方程中建立了环境温度与热变形间函数关系,使得变量递推过程中可补偿温度引起的热变形误差。最后,本发明的应变片贴在安装组件上,可在实验室中标定粘贴后的应变片灵敏度系数并校核,从而消除应变片粘贴方式、胶水种类等外因的影响。两片应变片采用半桥接线,可消除转轴弯曲变形影响,应变片通过滑环或无线传输等方式直接将数据传输至电脑中进行后续处理。On the one hand, the present invention adopts the weighted recursive least square estimation of estimation-correction to identify the propeller thrust from the strain signal with noise interference, which has strong anti-interference ability and can still maintain a high level in the environment of strong interference noise. measurement accuracy. On the other hand, the functional relationship between ambient temperature and thermal deformation is established in the recurrence equation, so that the thermal deformation error caused by temperature can be compensated during the variable recurrence process. Finally, the strain gauge of the present invention is pasted on the installation assembly, and the sensitivity coefficient of the pasted strain gauge can be calibrated and checked in the laboratory, thereby eliminating the influence of external factors such as the pasting method of the strain gauge and the type of glue. The two strain gauges adopt half-bridge connection, which can eliminate the influence of the bending deformation of the shaft. The strain gauges directly transmit the data to the computer through slip rings or wireless transmission for subsequent processing.

附图说明Description of drawings

通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiment. The drawings are only for the purpose of illustrating a preferred embodiment and are not to be considered as limiting the invention.

图1为本发明的安装组件的结构图;Fig. 1 is the structural diagram of mounting assembly of the present invention;

图2为本发明的安装组件的另一实施方式的结构图;Fig. 2 is a structural diagram of another embodiment of the mounting assembly of the present invention;

图3为本发明的船舶螺旋桨轴推力高精度实时辨识装置的结构图;Fig. 3 is the structural diagram of the ship propeller shaft thrust high-precision real-time identification device of the present invention;

图4为本发明的应变片半桥接线示意图;Fig. 4 is the schematic diagram of strain gauge half-bridge wiring of the present invention;

图5为本发明的一种船舶螺旋桨轴推力高精度实时辨识方法;Fig. 5 is a kind of ship propeller shaft thrust high-precision real-time identification method of the present invention;

图6为本发明的应变与推力辨识效果比较图。Fig. 6 is a comparative diagram of strain and thrust identification effects of the present invention.

附图标记说明Explanation of reference signs

1:船舶螺旋桨轴;2:安装组件;3:连接孔;4:应变片。1: Ship propeller shaft; 2: Mounting component; 3: Connection hole; 4: Strain gauge.

具体实施方式Detailed ways

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对照附图说明本发明的具体实施方式。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图,并获得其他的实施方式。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the specific implementation manners of the present invention will be described below with reference to the accompanying drawings. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention, and those skilled in the art can obtain other accompanying drawings based on these drawings and obtain other implementations.

为使图面简洁,各图中只示意性地表示出了与本发明相关的部分,它们并不代表其作为产品的实际结构。另外,以使图面简洁便于理解,在有些图中具有相同结构或功能的部件,仅示意性地绘示了其中的一个,或仅标出了其中的一个。在本文中,“一个”不仅表示“仅此一个”,也可以表示“多于一个”的情形。In order to make the drawing concise, each drawing only schematically shows the parts related to the present invention, and they do not represent the actual structure of the product. In addition, to make the drawings concise and easy to understand, in some drawings, only one of the components having the same structure or function is schematically shown, or only one of them is marked. Herein, "a" not only means "only one", but also means "more than one".

以下结合附图和具体实施例对本发明提出的一种船舶螺旋桨轴1推力高精度实时辨识方法及装置作进一步详细说明。根据下面说明和权利要求书,本发明的优点和特征将更清楚。The method and device for high-precision real-time identification of the propeller shaft 1 thrust of a ship proposed by the present invention will be further described in detail below in conjunction with the accompanying drawings and specific embodiments. Advantages and features of the present invention will be apparent from the following description and claims.

实施例1Example 1

参看图5,本实施例提供一种船舶螺旋桨轴1推力高精度实时辨识方法。在实施本实施例之前需要在船舶螺旋桨轴1选取一测点位置,通常将测点位置选在船舶螺旋桨轴1靠近螺旋桨处,并测量测点处的轴段直径以获得截面积A。Referring to FIG. 5 , this embodiment provides a high-precision real-time identification method for the propeller shaft 1 thrust of a ship. Before implementing this embodiment, it is necessary to select a measuring point on the propeller shaft 1 of the ship. Usually, the measuring point is selected near the propeller shaft 1 of the ship, and the diameter of the shaft section at the measuring point is measured to obtain the cross-sectional area A.

根据船舶螺旋桨轴1的轴段直径在测点位置安装一安装环,如图1、图2所示。参看图2至图4,安装环上粘贴有至少2片应变片4,且2片应变片4分别安装于安装环的两侧,即相对设置。安装完成后,在实验室中校核应变片4的灵敏度系数(校核时安装环还未安装在船舶螺旋桨轴1上)。安装环利用螺栓-螺母固定在轴上,如图3所示。对安装环上的两片应变片4采用半桥接线方式连接,以此消除转轴弯曲造成的弯曲变形影响,至此完成准备阶段。According to the diameter of the shaft section of the propeller shaft 1 of the ship, install a mounting ring at the position of the measuring point, as shown in Fig. 1 and Fig. 2 . Referring to FIGS. 2 to 4 , at least two strain gauges 4 are pasted on the mounting ring, and the two strain gauges 4 are installed on both sides of the mounting ring respectively, that is, oppositely arranged. After the installation is completed, check the sensitivity coefficient of the strain gauge 4 in the laboratory (the installation ring has not been installed on the propeller shaft 1 of the ship when checking). The mounting ring is secured to the shaft with bolt-nuts, as shown in Figure 3. The two strain gauges 4 on the mounting ring are connected by a half-bridge connection, so as to eliminate the influence of bending deformation caused by the bending of the rotating shaft, and the preparatory stage is completed so far.

参看图5,在本实施例中,首先,需要获取安装于船舶螺旋桨轴1上的至少2片应变片4的应变信号,以及船舶螺旋桨轴1的环境温度。然后根据应变信号和环境温度构建船舶螺旋桨轴1推力与其纵向应变的第一等式。Referring to FIG. 5 , in this embodiment, first, the strain signals of at least two strain gauges 4 installed on the propeller shaft 1 of the ship and the ambient temperature of the propeller shaft 1 of the ship need to be acquired. Then the first equation of the thrust of the propeller shaft 1 of the ship and its longitudinal strain is constructed according to the strain signal and the ambient temperature.

当船舶螺旋桨轴1受到螺旋桨推力与弯矩作用时,两个不同应变片4的度数为:When the ship propeller shaft 1 is subjected to propeller thrust and bending moment, the degrees of two different strain gauges 4 are:

ε1=εFM+a·(T-T0);ε2=εFM+a·(T-T0)ε 1FM +a·(TT 0 ); ε 2FM +a·(TT 0 )

其中,ε1、ε2分别为应变片4R1、R2产生的应变,εF为船舶螺旋桨轴1推力所产生的纵向应变,εM为弯矩产生的应变,a为轴系热膨胀系数,T0为环境初始温度,T为环境温度。Among them, ε 1 and ε 2 are the strains produced by the strain gauges 4R 1 and R 2 respectively, ε F is the longitudinal strain produced by the thrust of the propeller shaft 1 of the ship, ε M is the strain produced by the bending moment, a is the thermal expansion coefficient of the shaft system, T 0 is the initial temperature of the environment, and T is the ambient temperature.

应变仪分别与不同应变片4电连接,其读数ε的满足:The strain gauges are electrically connected to different strain gauges 4 respectively, and the readings of the ε gauges meet:

ε=ε12=2εF+2a·(T-T0)。 εmeter12 =2ε F +2a·(TT 0 ).

从而得到船舶螺旋桨轴1推力FT与纵向应变εF之间的关系满足以下第一等式:Thus, the relationship between the thrust FT of the propeller shaft 1 of the ship and the longitudinal strain ε F satisfies the following first equation:

Figure BDA0003931644910000081
Figure BDA0003931644910000081

其中,E与A分别为应变片4安装处轴段的弹性模量与截面积。Wherein, E and A are respectively the elastic modulus and cross-sectional area of the axial section where the strain gauge 4 is installed.

然后,对第一等式进行求导并按时间步离散得到第二等式,具体公式为:Then, deriving the first equation and discretizing it according to time steps to obtain the second equation, the specific formula is:

Figure BDA0003931644910000082
Figure BDA0003931644910000082

其中,ε(k+1)为第k+1个时间步下的应变仪读数,ε(k)为第k个时间步下的应变仪读数,ΔFT(k)为第k个时间步下的推力变化值,ΔT(k)为第k个时间步下环境温度变化值。Among them, εmeter (k+1) is the strain gauge reading at the k+1th time step, εmeter (k) is the strain gauge reading at the kth time step, and ΔF T (k) is the kth time The change value of the thrust at the next step, ΔT(k) is the change value of the ambient temperature at the kth time step.

进而,基于第二等式表示成矩阵方程和观测方程。Furthermore, it is expressed as a matrix equation and an observation equation based on the second equation.

令Xk+1=[ε(k+1)ΔFT(k)ΔT(k+1)]T,则可表示成如下矩阵方程:Let X k+1 =[ε(k+1)ΔF T (k)ΔT(k+1)] T , it can be expressed as the following matrix equation:

Xk+1=Φ·Xk+wk X k+1 =Φ·X k +w k

其中,

Figure BDA0003931644910000083
为状态转移矩阵,wk为模拟噪声的模型误差,反映建模误差对结果影响。in,
Figure BDA0003931644910000083
is the state transition matrix, w k is the model error of the simulated noise, reflecting the influence of the modeling error on the result.

观测方程可表示为The observation equation can be expressed as

yk=H·Xk+vk y k =H·X k +v k

其中,yk为测量的应变与测量的温度的变化,

Figure BDA0003931644910000084
为观测方程,vk为测量噪声,反映测量误差对结果影响。where yk is the change in measured strain versus measured temperature,
Figure BDA0003931644910000084
is the observation equation, and v k is the measurement noise, which reflects the influence of measurement errors on the results.

上述两个方程分别组成预估-校正加权递推最小二乘估计的状态方程和观测方程。递推计算过程如下:The above two equations compose the state equation and the observation equation of the predictor-correction weighted recursive least squares estimation respectively. The recursive calculation process is as follows:

在第k个时间步,完成状态一步预测,得到第k个时间步相对应的先验估计,具体公式为:At the kth time step, the one-step prediction of the state is completed, and the prior estimate corresponding to the kth time step is obtained. The specific formula is:

Figure BDA0003931644910000085
Figure BDA0003931644910000085

Figure BDA0003931644910000086
Figure BDA0003931644910000086

其中,

Figure BDA0003931644910000087
为X在第k个时间步的先验估计,
Figure BDA0003931644910000088
为X在第k-1个时间步的后验估计,
Figure BDA0003931644910000089
Figure BDA00039316449100000810
的估计误差协方差矩阵,
Figure BDA00039316449100000811
Figure BDA0003931644910000091
的估计误差协方差矩阵,Q为模型误差wk的协方差矩阵。in,
Figure BDA0003931644910000087
is the prior estimate of X at the kth time step,
Figure BDA0003931644910000088
is the posterior estimate of X at the k-1th time step,
Figure BDA0003931644910000089
for
Figure BDA00039316449100000810
The estimated error covariance matrix of ,
Figure BDA00039316449100000811
for
Figure BDA0003931644910000091
The estimated error covariance matrix of Q is the covariance matrix of the model error w k .

接着,在第k个时间步,利用先验估计

Figure BDA0003931644910000092
完成状态估计,具体公式为:Then, at the kth time step, use the prior estimate
Figure BDA0003931644910000092
Complete the state estimation, the specific formula is:

Kk=Pk|k-1HT(HPk|k-1HT+R)-1 K k =P k|k-1 H T (HP k|k-1 H T +R) -1

Figure BDA0003931644910000093
Figure BDA0003931644910000093

Figure BDA0003931644910000094
Figure BDA0003931644910000094

Pk|k=(1-KkH)Pk|k-1 P k|k =(1-K k H)P k|k-1

其中,Kk为增益矩阵,Zk为测量信号与先验估计的偏差序列,R为测量噪声vk的协方差矩阵。Among them, K k is the gain matrix, Z k is the deviation sequence between the measurement signal and the prior estimate, and R is the covariance matrix of the measurement noise v k .

将时间递推至k+1个时间步,重复状态预测与状态估计,以此得到新的估计变量,基于新的估计变量与第一等式计算得到船舶螺旋桨轴1的轴系推力。在递推计算中,需估计建模误差引起的协方差矩阵,以及测量噪声引起的协方差矩阵。然而,实际应用中,这两个矩阵不易确定,需根据轴系及现场环境调整。本实施例给出一种可确定协方差矩阵Q与协方差矩阵R的公式,具体为:Recurse the time to k+1 time steps, repeat the state prediction and state estimation, so as to obtain new estimated variables, and calculate the shafting thrust of ship propeller shaft 1 based on the new estimated variables and the first equation. In the recursive calculation, the covariance matrix caused by modeling error and the covariance matrix caused by measurement noise need to be estimated. However, in practical applications, these two matrices are not easy to determine and need to be adjusted according to the shafting and site environment. This embodiment provides a formula that can determine the covariance matrix Q and the covariance matrix R, specifically:

Figure BDA0003931644910000095
Figure BDA0003931644910000095

R=Ck-HPk|k-1HT R=C k -HP k|k-1 H T

Figure BDA0003931644910000096
Figure BDA0003931644910000096

Figure BDA0003931644910000097
Figure BDA0003931644910000097

图6为本方法在信噪比(SNR)低至20dB时的应变与推力辨识效果,推力辨识最大相对误差仅为4.8%Figure 6 shows the strain and thrust identification effect of this method when the signal-to-noise ratio (SNR) is as low as 20dB, and the maximum relative error of thrust identification is only 4.8%

实施例2Example 2

参看图1至图4,本实施例提供一种船舶螺旋桨轴推力高精度实时辨识装置,配置有实施例1的船舶螺旋桨轴推力高精度实时辨识方法。其包括安装组件2和测试件。安装组件2用于安装于在船舶螺旋桨轴1上测试点处,测试件则贴设于安装组件2上。Referring to FIGS. 1 to 4 , this embodiment provides a device for high-precision real-time identification of ship propeller shaft thrust, which is equipped with the high-precision real-time identification method for ship propeller shaft thrust of Embodiment 1. It includes the mounting assembly 2 and the test piece. The installation assembly 2 is used to be installed at the test point on the propeller shaft 1 of the ship, and the test piece is attached to the installation assembly 2 .

参看图1和图2,在本实施例中,具体地,安装组件2可拆分为第一安装件和第二安装件进行描述。第一安装件与第二安装件结构相同,均包括第一弧形件、第二弧形件和至少一连接杆。第一弧形件和第二弧形件均为半圆形的结构,连接杆的两端分别与第一弧形件一侧和第二弧形件一侧连接。当只有一根连接杆时,连接杆两端分别与第一弧形件的中间与第二弧形件的中间连接,当有两根连接杆时,连接杆两端分别将第一弧形件和第二弧形件均分为三段,连接点为三段的分割处。第一弧形件与第二弧形件对称设置,开口朝向完全相同,即第一弧形件平移后即可得到第二弧形件。Referring to FIG. 1 and FIG. 2 , in this embodiment, specifically, the mounting assembly 2 can be divided into a first mounting part and a second mounting part for description. The first mounting part and the second mounting part have the same structure, and both include a first arc-shaped part, a second arc-shaped part and at least one connecting rod. Both the first arc-shaped piece and the second arc-shaped piece have a semicircular structure, and the two ends of the connecting rod are respectively connected to one side of the first arc-shaped piece and one side of the second arc-shaped piece. When there is only one connecting rod, the two ends of the connecting rod are respectively connected with the middle of the first arc-shaped piece and the middle of the second arc-shaped piece; when there are two connecting rods, the two ends of the connecting rod connect the first arc-shaped piece and the second arc are divided into three sections, and the connection point is the division of the three sections. The first arc-shaped part and the second arc-shaped part are arranged symmetrically, and the openings have the same orientation, that is, the second arc-shaped part can be obtained after the first arc-shaped part is translated.

至少一测试件安装于连接杆中心,可以只设置一个也可以设置多个,使用多个应变片4可消除随机误差影响。同样的,当连接杆设置为一根时,连接杆上相对应地设有一个或多个测试件;当连接杆设置多根时,每根连接杆上也相应地设有一个或多个测试件。At least one test piece is installed in the center of the connecting rod, and there may be only one or more than one test piece, and the use of multiple strain gauges 4 can eliminate the influence of random errors. Similarly, when there is one connecting rod, one or more test pieces are correspondingly provided on the connecting rod; when multiple connecting rods are provided, one or more test pieces are correspondingly provided on each connecting rod. pieces.

第一弧形件和第二弧形件的两端均设有连接孔3,连接孔3的开设方向与连接杆相互垂直。第一安装件的第一弧形件与第二安装件的第一弧形件通过螺栓螺母穿过两者的连接孔3,从而实现螺栓连接,第二弧形件同理。从而使得第一安装件和第二安装件相互连接,从而构成环形的安装组件2,环绕船舶螺旋桨轴1贴合设置,如图3所示。Both ends of the first arc-shaped piece and the second arc-shaped piece are provided with connecting holes 3, and the opening direction of the connecting holes 3 is perpendicular to the connecting rod. The first arc-shaped part of the first mounting part and the first arc-shaped part of the second mounting part pass through the connection holes 3 of the two through bolts and nuts, so as to realize bolt connection, and the same is true for the second arc-shaped part. In this way, the first mounting part and the second mounting part are connected to each other, thereby forming a ring-shaped mounting assembly 2, which is fitted around the propeller shaft 1 of the ship, as shown in FIG. 3 .

另外,参看图4测试件为应变片4,且两个相对设置的应变片4电连接,采用半桥接线,用以消除转轴弯曲变形影响。In addition, referring to Fig. 4, the test piece is a strain gauge 4, and two opposite strain gauges 4 are electrically connected, and a half-bridge line is used to eliminate the influence of the bending deformation of the rotating shaft.

上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式。即使对本发明作出各种变化,倘若这些变化属于本发明权利要求及其等同技术的范围之内,则仍落入在本发明的保护范围之中。The embodiments of the present invention have been described in detail above with reference to the accompanying drawings, but the present invention is not limited to the above embodiments. Even if various changes are made to the present invention, if these changes fall within the scope of the claims of the present invention and equivalent technologies, they still fall within the protection scope of the present invention.

Claims (9)

1. A high-precision real-time identification method for thrust of a ship propeller shaft is characterized by comprising the following steps:
acquiring strain signals of at least 2 strain gauges arranged on a propeller shaft of a ship and the ambient temperature of the propeller shaft of the ship;
constructing a first equation of the thrust of the propeller shaft of the ship and the longitudinal strain of the propeller shaft of the ship according to the strain signal and the environment temperature; carrying out derivation on the first equation and dispersing the first equation according to time steps to obtain a second equation; expressing into a matrix equation and an observation equation based on the second equation;
state prediction is completed based on the matrix equation and the observation equation, and prior estimation corresponding to the kth time step is obtained; completing state estimation at the kth time step based on the prior estimation to obtain an estimation variable corresponding to the kth time step; and (4) recursion is carried out on the time to k +1 time steps, state prediction and state estimation are repeated, so that a new estimation variable is obtained, and shafting thrust of the ship propeller shaft is obtained through calculation based on the new estimation variable and the first equation.
2. The method for identifying the thrust force of the ship propeller shaft in high precision in real time according to claim 1, wherein the specific steps of constructing the first equation of the thrust force of the ship propeller shaft and the longitudinal strain of the ship propeller shaft according to the strain signal and the environment temperature comprise
Respectively reading the strain signals of different strain gauges, wherein the degree of the strain signals satisfies the following formula:
ε 1 =ε FM +a·(T-T 0 );ε 2 =ε FM +a·(T-T 0 )
wherein epsilon 1 、ε 2 Are respectively strain gauges R 1 、R 2 Strain generated, epsilon F Longitudinal strain, epsilon, generated by thrust of propeller shaft of ship M A is the coefficient of thermal expansion of the shafting, T 0 Is an ambient initial temperature, T is the ambient temperature;
readings epsilon of strain gauges respectively connected with different strain gauges Instrument for measuring the shape of a human body The formula of (1) is:
ε instrument for measuring the shape of a human body =ε 12 =2ε F +2a·(T-T 0 )
Propeller shaft thrust F of ship T With longitudinal strain epsilon F The relationship satisfies the first equation:
Figure FDA0003931644900000011
wherein E and A are respectively the elastic modulus and the sectional area of the shaft section at the mounting position of the strain gauge.
3. The method for identifying the thrust of the propeller shaft of the ship in real time with high precision as claimed in claim 2, wherein the first equation is derived and dispersed in time steps to obtain a second equation, and the specific formula is as follows:
Figure FDA0003931644900000021
wherein epsilon Instrument for measuring the shape of a human body (k + 1) is the strain gauge reading at time step k +1, ε Instrument for measuring the shape of a human body (k) For strain gauge readings at the kth time step, Δ F T (k) The thrust variation value at the kth time step is shown, and Δ T (k) is the ambient temperature variation value at the kth time step.
4. The method for high-precision real-time identification of thrust of a propeller shaft of a marine vessel as claimed in claim 3, wherein said expression into a matrix equation based on said second equation is specifically
Let X k+1 =[ε Instrument for measuring the shape of a human body (k+1)ΔF T (k) ΔT(k+1)] T Then the matrix equation is
X k+1 =Φ·X k +w k
Wherein,
Figure FDA0003931644900000022
is a state transition matrix, wk isSimulating a model error of the noise;
expressing an observation equation based on the second equation specifically is
y k =H·X k +v k
Wherein, y k For changes in measured strain and measured temperature,
Figure FDA0003931644900000023
to observe the equation, v k To measure noise.
5. The method for identifying the thrust of the propeller shaft of the ship in high precision in real time according to claim 4, wherein at the kth time step, the obtaining of the prior estimate corresponding to the kth time step specifically comprises:
Figure FDA0003931644900000024
Figure FDA0003931644900000025
wherein,
Figure FDA0003931644900000026
for the a priori estimate of X at the kth time step,
Figure FDA0003931644900000027
for the a posteriori estimate of X at the k-1 time step,
Figure FDA0003931644900000028
is composed of
Figure FDA0003931644900000029
The covariance matrix of the estimated error of (a),
Figure FDA00039316449000000210
is composed of
Figure FDA0003931644900000031
Q is the model error w k The covariance matrix of (a);
at the kth time step, the obtaining of the estimated variable corresponding to the kth time step specifically includes:
K k =P k|k-1 H T (HP k|k-1 H T +R) -1
Figure FDA0003931644900000032
Figure FDA0003931644900000033
P k|k =(1-K k H)P k|k-1
wherein, K k Is a gain matrix, Z k For a sequence of deviations of the measurement signal from said a priori estimate, R is the measurement noise v k The covariance matrix of (a);
the covariance matrix Q and the covariance matrix R have the respective formulae:
Figure FDA0003931644900000034
R=C k -HP k|k-1 H T
Figure FDA0003931644900000035
Figure FDA0003931644900000036
6. the method for identifying the thrust of the propeller shaft of the ship in real time according to any one of claims 1 to 5, wherein the method further comprises the following steps before acquiring the strain signals of at least 2 strain gauges installed on the propeller shaft of the ship and the ambient temperature of the propeller shaft of the ship
Selecting a measuring point position on a propeller shaft of the ship, and measuring the diameter of a shaft section at the measuring point position to obtain the sectional area A;
arranging a mounting ring at the measuring point position to paste 2 strain gauges and checking the sensitivity coefficient of the strain gauge;
and connecting the two strain gauges by adopting a half-bridge connection line mode.
7. A high-precision real-time identification device for the thrust of a ship propeller shaft, which is configured with the high-precision real-time identification method for the thrust of the ship propeller shaft as claimed in claims 1 to 6,
the device comprises a mounting assembly and a test piece; the mounting assembly is mounted at a test point on a propeller shaft of the ship; the test piece is attached to the mounting assembly;
the mounting assembly comprises a first mount and a second mount;
the first mounting part and the second mounting part are same in structure and respectively comprise a first arc-shaped part, a second arc-shaped part and at least one connecting rod;
the first arc-shaped part and the second arc-shaped part are both semicircular, two ends of the connecting rod are respectively connected with one side of the first arc-shaped part and one side of the second arc-shaped part, and the first arc-shaped part and the second arc-shaped part are symmetrically arranged; at least one test piece is arranged on the connecting rod;
the two ends of the first arc-shaped part and the second arc-shaped part are provided with connecting holes, the opening directions of the connecting holes are perpendicular to the connecting rod, and the first mounting part and the second mounting part are matched with a bolt and a nut through the connecting holes to realize bolt connection.
8. The device for identifying the thrust of the propeller shaft of the ship in real time according to claim 7, wherein the first mounting piece and the second mounting piece are provided with two connecting rods;
two ends of the two connecting rods are respectively connected with one side of the first arc-shaped part and one side of the second arc-shaped part, and the connecting positions of the two connecting rods divide the first arc-shaped part or the second arc-shaped part into three sections of arcs; each connecting rod is provided with at least one test piece.
9. The device for identifying the thrust of the propeller shaft of the ship in real time according to claim 7 or 8, wherein the test piece is a strain gauge, and the two strain gauges arranged oppositely are electrically connected and adopt a half-bridge connection.
CN202211390047.4A 2022-11-08 2022-11-08 A method and device for high-precision real-time identification of ship propeller shaft thrust Active CN115636063B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211390047.4A CN115636063B (en) 2022-11-08 2022-11-08 A method and device for high-precision real-time identification of ship propeller shaft thrust

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211390047.4A CN115636063B (en) 2022-11-08 2022-11-08 A method and device for high-precision real-time identification of ship propeller shaft thrust

Publications (2)

Publication Number Publication Date
CN115636063A true CN115636063A (en) 2023-01-24
CN115636063B CN115636063B (en) 2025-02-25

Family

ID=84948602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211390047.4A Active CN115636063B (en) 2022-11-08 2022-11-08 A method and device for high-precision real-time identification of ship propeller shaft thrust

Country Status (1)

Country Link
CN (1) CN115636063B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017089644A1 (en) * 2015-11-26 2017-06-01 Wärtsilä Finland Oy Marine vessel performance monitoring
KR20220001280A (en) * 2020-06-29 2022-01-05 한국해양대학교 산학협력단 System and method for prediction diagnosis of ship propeller and shaft condition using digital twin

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017089644A1 (en) * 2015-11-26 2017-06-01 Wärtsilä Finland Oy Marine vessel performance monitoring
KR20220001280A (en) * 2020-06-29 2022-01-05 한국해양대학교 산학협력단 System and method for prediction diagnosis of ship propeller and shaft condition using digital twin

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
徐鹏: "船舶推进轴系弯—纵耦合效应的非线性振动特性", 中国舰船研究, no. 2019, 31 October 2019 (2019-10-31), pages 49 - 57 *
李惊东;陈文俊;苍晓羽;: "一种舰船推力测量装置的结构和电磁兼容性设计", 中国舰船研究, no. 01, 21 January 2014 (2014-01-21), pages 125 - 130 *
邹冬林: "桨-轴系统双向流固耦合动力学分析及其纵向激励力特性研究", 博士电子期刊, no. 2019, 15 August 2019 (2019-08-15), pages 1 - 185 *
邹冬林;张建波;塔娜;饶柱石;: "船舶推进轴系纵横耦合非线性动力学分析――叶频激励下横向主共振响应", 船舶力学, no. 02, 15 February 2017 (2017-02-15), pages 77 - 86 *

Also Published As

Publication number Publication date
CN115636063B (en) 2025-02-25

Similar Documents

Publication Publication Date Title
US11860066B2 (en) Dynamic strain field measuring method and system for rotor blade based on blade tip timing
US4550617A (en) Multi axis force and moments transducer
GB2592775A (en) Rotating blade non-contact dynamic strain field measurement method and system
CN112033672B (en) Calibration device and calibration method for static and dynamic load identification of ship radial bearing
US5922967A (en) Method and apparatus for estimating loads imposed on structural body
CN212363509U (en) Fastener with stress sensing function and film transduction sensing system
CN112326081A (en) Method and device for monitoring torque of aero-engine
CN103728136B (en) Bush(ing) bearing oil film dynamic stiffness on-line testing method
CN115077901B (en) Bearing dynamic load indirect measurement method considering basic excitation
CN115636063A (en) High-precision real-time identification method and device for thrust of ship propeller shaft
CN115560898A (en) Method and device for measuring axial thrust of ship propeller
CN112857632B (en) Axial force measuring equipment
CN117269326B (en) An ultrasonic stress measurement method with temperature detection
CN109357808B (en) Test platform for testing excitation force of underwater impeller
CN208860596U (en) Aluminum enclosure bearing block Thermal-mechanical Coupling strain measurement system
CN109612608B (en) Installation method of strain bridge for rotating shaft and underwater impeller excitation force measuring platform
CN115452214A (en) Harmonic reducer flexible gear stress-strain measurement and method
Islam et al. Method of torque measurement based on volumetric strain
Chahmi Study and realization of a torque measurement sensor based on strain gauges dedicated to bending and torsion
CN114184352B (en) Nacelle power meter model test balance device and use method
CN115200859B (en) Torque meter and use method thereof
KR102608312B1 (en) Thrust and torque measurement equipment using data post-processing for ship model tests
CN113607299B (en) Structural health monitoring system and method for pressure container
Séchaud Development of a new rig with enhanced specifications for wind tunnel tests of “isolated” CROR propulsion systems
CN112781763B (en) Ship shaft power measurement method based on Monte Carlo Kalman filtering

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