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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 239000011159 matrix material Substances 0.000 claims abstract description 37
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 230000035945 sensitivity Effects 0.000 claims abstract description 5
- 238000012360 testing method Methods 0.000 claims description 30
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000007704 transition Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims 1
- 230000008859 change Effects 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 3
- 239000003292 glue Substances 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000009434 installation Methods 0.000 description 8
- 238000005452 bending Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000012790 adhesive layer Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling 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个时间步,重复状态预测与状态估计,最终获得推力动态过程。本发明具有较强抗干扰能力,可补偿环境温度变化引起的热变形误差,能在强干扰噪声中保持较高测量精度。此外,本发明中应变片贴在测量装置上,可在实验室中标定应变片灵敏度系数,从而消除应变片粘贴方式、胶水种类等外因影响。
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.
Description
技术领域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=εF-εM+a·(T-T0);ε2=εF+εM+a·(T-T0)ε 1 =ε F -ε M +a·(TT 0 ); ε 2 =ε F +ε M +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:
ε仪=ε1+ε2=2εF+2a·(T-T0) εmeter =ε 1 +ε 2 =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:
其中,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:
其中,ε仪(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
其中,为状态转移矩阵,wk为模拟噪声的模型误差;in, 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为测量的应变与测量的温度的变化,为观测方程,vk为测量噪声。where yk is the change in measured strain versus measured temperature, 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:
其中,为X在第k个时间步的先验估计,为X在第k-1个时间步的后验估计,为的估计误差协方差矩阵,为的估计误差协方差矩阵,Q为模型误差wk的协方差矩阵;in, is the prior estimate of X at the kth time step, is the posterior estimate of X at the k-1th time step, for The estimated error covariance matrix of , for 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
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:
R=Ck-HPk|k-1HT R=C k -HP k|k-1 H T
进一步优选地,在获取安装于船舶螺旋桨轴上的至少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
参看图5,在本实施例中,首先,需要获取安装于船舶螺旋桨轴1上的至少2片应变片4的应变信号,以及船舶螺旋桨轴1的环境温度。然后根据应变信号和环境温度构建船舶螺旋桨轴1推力与其纵向应变的第一等式。Referring to FIG. 5 , in this embodiment, first, the strain signals of at least two
当船舶螺旋桨轴1受到螺旋桨推力与弯矩作用时,两个不同应变片4的度数为:When the ship propeller shaft 1 is subjected to propeller thrust and bending moment, the degrees of two
ε1=εF-εM+a·(T-T0);ε2=εF+εM+a·(T-T0)ε 1 =ε F -ε M +a·(TT 0 ); ε 2 =ε F +ε M +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
ε仪=ε1+ε2=2εF+2a·(T-T0)。 εmeter =ε 1 +ε 2 =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:
其中,E与A分别为应变片4安装处轴段的弹性模量与截面积。Wherein, E and A are respectively the elastic modulus and cross-sectional area of the axial section where the
然后,对第一等式进行求导并按时间步离散得到第二等式,具体公式为:Then, deriving the first equation and discretizing it according to time steps to obtain the second equation, the specific formula is:
其中,ε仪(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
其中,为状态转移矩阵,wk为模拟噪声的模型误差,反映建模误差对结果影响。in, 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为测量的应变与测量的温度的变化,为观测方程,vk为测量噪声,反映测量误差对结果影响。where yk is the change in measured strain versus measured temperature, 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:
其中,为X在第k个时间步的先验估计,为X在第k-1个时间步的后验估计,为的估计误差协方差矩阵,为的估计误差协方差矩阵,Q为模型误差wk的协方差矩阵。in, is the prior estimate of X at the kth time step, is the posterior estimate of X at the k-1th time step, for The estimated error covariance matrix of , for The estimated error covariance matrix of Q is the covariance matrix of the model error w k .
接着,在第k个时间步,利用先验估计完成状态估计,具体公式为:Then, at the kth time step, use the prior estimate 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
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:
R=Ck-HPk|k-1HT R=C k -HP k|k-1 H T
图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
参看图1和图2,在本实施例中,具体地,安装组件2可拆分为第一安装件和第二安装件进行描述。第一安装件与第二安装件结构相同,均包括第一弧形件、第二弧形件和至少一连接杆。第一弧形件和第二弧形件均为半圆形的结构,连接杆的两端分别与第一弧形件一侧和第二弧形件一侧连接。当只有一根连接杆时,连接杆两端分别与第一弧形件的中间与第二弧形件的中间连接,当有两根连接杆时,连接杆两端分别将第一弧形件和第二弧形件均分为三段,连接点为三段的分割处。第一弧形件与第二弧形件对称设置,开口朝向完全相同,即第一弧形件平移后即可得到第二弧形件。Referring to FIG. 1 and FIG. 2 , in this embodiment, specifically, the mounting
至少一测试件安装于连接杆中心,可以只设置一个也可以设置多个,使用多个应变片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
第一弧形件和第二弧形件的两端均设有连接孔3,连接孔3的开设方向与连接杆相互垂直。第一安装件的第一弧形件与第二安装件的第一弧形件通过螺栓螺母穿过两者的连接孔3,从而实现螺栓连接,第二弧形件同理。从而使得第一安装件和第二安装件相互连接,从而构成环形的安装组件2,环绕船舶螺旋桨轴1贴合设置,如图3所示。Both ends of the first arc-shaped piece and the second arc-shaped piece are provided with connecting
另外,参看图4测试件为应变片4,且两个相对设置的应变片4电连接,采用半桥接线,用以消除转轴弯曲变形影响。In addition, referring to Fig. 4, the test piece is a
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式。即使对本发明作出各种变化,倘若这些变化属于本发明权利要求及其等同技术的范围之内,则仍落入在本发明的保护范围之中。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)
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)
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 |
-
2022
- 2022-11-08 CN CN202211390047.4A patent/CN115636063B/en active Active
Patent Citations (2)
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)
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 |