CN111189427B - Real-time monitoring method for bending deformation of simply supported beam - Google Patents

Real-time monitoring method for bending deformation of simply supported beam Download PDF

Info

Publication number
CN111189427B
CN111189427B CN202010090290.9A CN202010090290A CN111189427B CN 111189427 B CN111189427 B CN 111189427B CN 202010090290 A CN202010090290 A CN 202010090290A CN 111189427 B CN111189427 B CN 111189427B
Authority
CN
China
Prior art keywords
value
angle
zero
point
filtering
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
CN202010090290.9A
Other languages
Chinese (zh)
Other versions
CN111189427A (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.)
Guilin University of Technology
Original Assignee
Guilin University of Technology
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 Guilin University of Technology filed Critical Guilin University of Technology
Priority to CN202010090290.9A priority Critical patent/CN111189427B/en
Publication of CN111189427A publication Critical patent/CN111189427A/en
Application granted granted Critical
Publication of CN111189427B publication Critical patent/CN111189427B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

The invention discloses a method for monitoring bending deformation of a simply supported beam in real time. The method for establishing the real-time monitoring device for the bending deformation of the simply supported beam comprises the following steps: an embedded terminal and a computer terminal; the embedded end is fixed on the simply supported beam to be detected and used for monitoring the deformation of the simply supported beam to be detected in real time; the computer end is used for remotely monitoring real-time data, providing a data interface of software, and calculating and displaying real-time angle values, maximum deflection positions, maximum deflection values, stress points, stress sizes, real-time model animations simulating the stress deformation of the beam and the like; the method carries out filtering processing of three stages of denoising processing, zero point processing and division number processing on the collected angle data; and a deformation detection mathematical model is established; the invention can be used for constructing bridges such as slab bridges, ribbed slab bridges, box girder bridges and the like, and mechanical analysis models of container suspension arms, and can also be used for monitoring the stressed bending deformation of the solid beam arms in real time, thereby overcoming the defects of measurement methods such as strain gauges, computer vision and the like.

Description

Real-time monitoring method for bending deformation of simply supported beam
Technical Field
The invention relates to the field of deformation monitoring, in particular to a method for monitoring bending deformation of a simply supported beam in real time.
Background
The simply supported beam is a supporting structure with two end supports only providing vertical restraint and not providing corner restraint. The simply supported beam is a bridge such as a slab bridge, a ribbed slab bridge, a box girder bridge and the like, and a mechanical model of the container suspension arm.
No matter the bridge or the suspension arm is an important bearing and stressed part, and the phenomenon of bending deformation of the bridge or the suspension arm under stress is particularly prominent. And serious deformation may cause a serious safety accident. Therefore, the stress monitoring of the beam arm is very important.
The existing measurement techniques mainly include: strain gauges and computer vision. The strain gauge measuring technology is mature, the technology is high in precision and long in service life, but as an implanted measuring method, the stress cannot be judged, the bending displacement cannot be measured, and meanwhile, the technology is poor in cost and flexibility. Computer vision measurement is a new technology in these years, and the technology has high requirements on the performance of computer vision observation and installation and debugging of an image acquisition device, and is easily influenced by visual occlusion.
Based on the defects of the technology, the invention provides a device and method for monitoring stress and deformation of a simply supported beam-arm model in real time.
Disclosure of Invention
The invention aims to provide a method for monitoring bending deformation of a simply supported beam in real time.
The method comprises the following specific steps:
(1) the method for establishing the real-time monitoring device for the bending deformation of the simply supported beam comprises the following steps: an embedded terminal and a computer terminal;
an embedded port, comprising: the device comprises a wireless angle sensing module, a microcontroller module, a wireless data transmission module and a display module. The microcontroller module is electrically connected with the wireless angle sensing module, the wireless data transmission module and the display module.
The embedded end is used for acquiring sensor data, analyzing the data, filtering the data, converting the data, displaying the operation results in the form of graphs and characters, monitoring parameters in real time and alarming a threshold value, wirelessly communicating with the computer end, sending the filtered data to the computer end and receiving a control instruction of the computer end;
the embedded end is fixed on the simply supported beam to be detected and used for monitoring the deformation of the simply supported beam to be detected in real time;
the wireless angle sensing module at least comprises two triaxial angle sensors, and each angle sensor can provide instantaneous angle change of three axes in a triaxial rectangular coordinate system. The angle sensors are placed at the positions of the two supporting points of the simply supported beam to be monitored, the angle values of the two supporting points are monitored in real time and transmitted back, and a wireless data transmission function is provided.
And the microcontroller module is used for data analysis, data filtering, data conversion operation, real-time parameter monitoring and threshold value alarming, and outputting operation results in the form of graphs and characters to the display module.
And the wireless data transmission module is used for carrying out wireless communication with the wireless angle sensing module, carrying out wireless communication with a computer end and carrying out wired communication with the microcontroller module.
And the display module displays the deformation process and the result in the form of graphs and characters.
Computer end, including two parts, the first part is: a data transmission transfer device; the first part is: computer and upper computer software.
And the computer end is used for remotely monitoring the real-time data, providing a data interface of software, carrying out wireless communication with the embedded end, receiving the data processed by the embedded end and sending a control instruction to the embedded end. And the computer-side upper computer software realizes the functions of receiving data transmitted back by the serial port, sending serial port coding instructions, displaying the sampled real-time angle value, the maximum deflection position, the maximum deflection value, the stress point and the stress magnitude on a graphical interface, simulating the real-time model animation of the stress deformation of the beam and the like.
A data transfer relay for: receiving data returned by the embedded terminal in a wireless mode, sending the data to a computer by using a serial port, and processing the data by using upper computer software; and receiving a control instruction sent by the computer through the serial port, and sending the control instruction to the embedded terminal in a wireless communication mode.
(2) The embedded terminal in the step (1) provides three functional modes: a parameter setting mode, a calibration function mode and a measurement function mode.
A parameter setting mode, the mode functions including: parameter setting, threshold setting, filtering enabling, wireless data transmission enabling, wireless channel detection and angle sensor detection.
And setting parameters, namely setting the length of the beam and the bending rigidity of the beam.
Setting a threshold value, namely setting a pressure threshold value point and a maximum deflection value threshold value point, and selecting whether to start a threshold value alarm function; and starting a threshold alarm function in a measurement function mode, and if the monitored pressure value and the maximum deflection value exceed the set threshold point, sending an alarm by the device and displaying an alarm mark on a display module.
And enabling filtering, and selectively turning on or off the filtering function.
And enabling the wireless data, and selecting to turn on or turn off the function of transmitting data back to the computer terminal through the wireless mode.
And detecting a wireless channel, namely detecting whether a wireless module exists or not and whether the function can be normally operated or not.
The angle sensor detects whether there is data returned by the angle sensor.
A calibration function mode, the mode functions being: calibrating an angle zero point, correcting angle deviation and calibrating the bending rigidity of the beam to be measured.
And calibrating an angle zero point, and sampling an angle value of the simply supported beam to be tested in a static state under the static condition that the simply supported beam to be tested is not subjected to external pressure, and calibrating the angle value into an angle zero point value.
And correcting the angle deviation for calibrating the simply supported beam model. The calibration process comprises the following steps: and (3) applying downward pressure at the center of the beam by taking the two supporting points as the end points of the beam, and solving and storing the difference between the angle values of the two ends after the angle values of the two ends are stable.
A measurement function mode, the functions of which are: the angle sensor data acquisition, data analysis, data filtering, data conversion operation, operation result display in the form of graph and numerical characters, real-time monitoring of angle change and threshold value alarm, wireless communication with a computer end, transmission of filtered data with the computer end, and reception of control instructions of the computer end;
(3) and (3) carrying out deformation detection on the simply supported beam to be detected by the simply supported beam bending deformation real-time monitoring device in the step (1), and carrying out filtering processing on angle data acquired by the wireless sensing module. The filtering processing method comprises the following steps:
in the filtering process, the angle value is divided into three stages for processing, which are respectively: angle bare value, internal division number and angle view. In total, three steps can be divided: denoising processing, zero point processing and division number processing.
The angle bare value is an angle value which is sampled by the angle sensor and is not processed, namely the angle most original data;
the internal part degree is an angle calculated quantity obtained after the angle bare value is subjected to denoising treatment and zero point treatment;
the angle view is the final filtering output quantity after the division number processing is carried out on the internal division number, and is used as the angle value output by the filtering processing.
Denoising processing for filtering out fluctuation causing instability of the sampling value, the denoising processing comprising: and performing zero tracking judgment, moving average filtering and weighted filtering on the angle bare value. The denoising processing steps are as follows:
q1, inputting angle bare value;
q2, zero-tracking range judgment processing;
q3, moving average filtering;
q4, weighted filtering;
q5, processing the angle change range;
q6, outputting an angle filtered value.
And zero tracking range judgment processing for judging whether the current angle bare value is in the zero tracking range. If the zero tracking range is within the zero tracking range, zero tracking processing is performed by zero point correction, and a zero tracking output value is output. If the angle exceeds the zero tracking range, the current angle bare value is cut off.
And (5) tracking a zero range, and taking the maximum change range of the angle value acquired under the static state that the simply supported beam to be detected is not subjected to the external force.
The zero-tracking process is to increase or decrease the angle value so that the angle value approaches 0 when the angle value is within a predetermined range.
And (3) moving average filtering, wherein a moving average filtering array is used for storing zero tracking output values as a filtering operation reference quantity, the size of the array is N, and N is a natural number which is greater than or equal to 10. And (3) carrying out filtering averaging by adopting N reference quantities, and putting a new zero-tracking output value into the first element of the array every time of carrying out moving average filtering.
For the moving average filter array, two filling operation modes are provided: fast refresh fill operation, slow refresh fill operation.
And (3) a fast refreshing filling operation mode, namely, putting the first three zero-tracing output values in the array into the last three bits of the array, and filling the rest array elements with the zero-tracing output values output by the last filtering.
And (4) slowly refreshing the filling operation mode, and transmitting the zero-chasing output values in the array one by one.
And (3) moving average filtering, combining two filling operation modes, firstly running M times of slow refreshing filling operation in a filtering processing period, and then executing Q times of fast refreshing filling operation. Wherein M is a natural number greater than or equal to 2, and Q is a natural number greater than or equal to 2. After each fill, the N elements in the array are averaged and the average is output to a weighted filter.
Weighting filtering, which is used for endowing the angle bare value with different weight values in different change processes;
weighting and filtering, wherein if the simply supported beam component to be measured is in a stable state, the newly acquired angle bare value weight is weakened; and if the beam member is in an unstable state, reinforcing the newly acquired angle bare value weight.
Weighted filtering, namely judging whether the stable state of the internal division number is achieved or not;
the stable state of the internal division number is that the difference between the newly input angle bare value and the last angle visual number output by filtering does not exceed a set range.
If the stable state of the internal part number is reached, the filtering calculation formula is as follows: the filter output value (newly input angle bare value 1+ last filter output angle view 3)/4; in the above equation, the weight of the angle bare value of the new input is 25%.
If the steady state of the internal division number is not reached, the filter calculation formula is as follows: the filter output value (newly input angle bare value x 3+ last filter output angle view x 1)/4; in the above equation, the weight of the angle bare value newly input is 75%.
And the angle change range processing function is used for comparing the newly input angle bare value with the last filtering output value, judging whether the angle value is stable according to the change range interval, and determining the filtering times required when the internal division number reaches a stable state.
Zero point processing is used for correcting angle zero point deviation caused by the simply supported beam to be measured in the working process, and measuring deviation caused by the following conditions is solved: 1) the shape of the simply supported beam to be measured changes when the simply supported beam returns to the static state after being stressed and deformed every time, so that the angle zero value of the simply supported beam to be measured acquired by the angle sensor in the static state can deviate from the calibration zero value in the calibration process; 2) the mechanical loss of the simply supported beam to be tested in a stressed state for a long time causes the deviation between the actual two-end angle zero value of the simply supported beam to be tested in a static state and the calibration zero value in a calibration state to be gradually increased.
Zero point processing, comprising: zero point correction and automatic zero tracking.
Zero point correction, comprising three functions: correcting a calibration zero difference value, correcting an accumulated zero value and filtering an angle difference value at two ends.
And correcting the calibration zero point, and subtracting the calibration zero point value during calibration from the angle zero point value acquired each time.
And correcting the accumulated zero point, and calculating the zero tracking amount generated after zero tracking into the angle bare value to change the angle zero point value of the simply supported beam to be measured in a static state.
And filtering the angle difference value of the two ends, and detecting the bending deformation of the simply supported beam to be detected by calculating the relative change value of the angle values of the two ends, so that the simply supported beam to be detected is not influenced by the inclination of the simply supported beam.
And the automatic zero tracking is used for carrying out zero tracking calculation on the current angle value within the zero tracking range according to the set time frequency. If the current angle value is larger than the zero point value, the zero chasing amount is automatically increased, otherwise, the zero chasing amount is automatically decreased. When the change of the current angle value is within the set range, the change is stable.
The automatic zero-tracking process is as follows:
z1, judging whether the zero-chasing tracking time is reached, if so, entering Z2, otherwise, jumping out of the automatic zero-chasing program;
z2, judging whether the current angle value is stable, if so, entering Z3, otherwise, jumping out of the automatic zero-tracing program;
z3, judging whether the current angle value is less than or equal to zero-chasing range, if so, entering Z4, otherwise, jumping out of the automatic zero-chasing program;
z4, judging whether the current zero point is larger than the minimum difference value, if so, entering Z6, and otherwise, entering Z5;
z5, catching up the end value and completing an automatic zero catching program;
z6, judging whether the current angle value is positive, if so, entering Z8, and if not, entering Z7;
z7, the zero-chasing quantity is reduced, and the automatic zero-chasing procedure is completed;
and Z8, the zero chasing amount is automatically increased, and the automatic zero chasing program is completed.
The last value is traced for zeroing the current angle value.
After the zero point processing is executed, converting the angle value after the zero point processing into an internal division number, wherein the formula is as follows: CntIn _ Cur is Angle/32768 180 °, where Angle represents the Angle value after filtering and CntIn _ Cur represents the internal index number.
And index number processing is used for avoiding unstable measurement data caused by the phenomenon of stress deformation and jitter of the simply supported beam to be measured.
And (4) index number processing, namely performing filtering operation by judging the variation range of the index value and according to the variation range.
The differential value of the output value is processed by the fractal value, the internal division number and the division number.
Index number processing, setting 3 index number change critical values: 3C, 6C and 10C. C is the basic quantity of the set index number change, wherein C is a positive integer multiple of 0.5 degrees, the 3C critical point represents the value of 3 times the basic quantity, the 6C critical point represents the value of 6 times the basic quantity, and the 10C critical point represents the value of 10 times the basic quantity.
The division number processing program flow is as follows:
n1, judging whether the visual value is less than 3C, if so, jumping out of the index number processing program, otherwise, entering N2;
n2, judging whether the score value is less than 6C, if so, entering N2-1, otherwise, entering N3;
n2-1, accumulating the internal division number for 8 times, and obtaining the average value of 8 times as the output of division number processing;
n3, judging whether the score value is less than 10C, if so, entering N3-1, otherwise, entering N3-2;
n3-1, accumulating the internal division number for 15 times, and obtaining the average value of the 15 times as the output of division number processing;
n3-2, accumulating the internal division number for 3 times, and obtaining the average value of the 3 times as the output of division number processing;
(4) in the method for monitoring the bending deformation of the simply supported beam in real time, a mathematical model for detecting the deformation of the simply supported beam needs to be established, and the mathematical model for detecting the deformation of the simply supported beam is as follows:
key to the need for computationThe parameters are as follows: the magnitude of the vertical downward resultant force F acting on the beam surface of the simply supported beam to be measuredpResultant force FpThe beam axis position a and the maximum deflection value wmaxAnd the beam axis position d where the maximum deflection is located, and the calculated value needs to be displayed in real time through a screen and alarm is realized.
A plane rectangular coordinate system xow is constructed, the plane rectangular coordinate system xow uses the plane of the beam axis, a pivot point a of the beam as an origin O, and the horizontal extension direction of the beam as the positive half axis direction of the X axis to establish a plane rectangular coordinate system XOW, and the plane rectangular coordinate system XOW has the X axis and the W axis which are perpendicular to each other.
And the fitting calculation equation of the bending deformation shape of the simply supported beam to be detected on the stress plane is the bending line w (x) of the beam.
The resultant force received by the deformation of the simply supported beam to be measured is a downward resultant force F vertical to the beam surfacep. The two supporting points of the simply supported beam to be measured are point A and point B, the length of the beam is A, B straight line distance L between the two supporting points, the connecting line of the two supporting points A, B is the axis before beam deformation, the axis of the beam is on the X axis of the plane coordinate system XOW,
θAthe included angle between the tangent line of the deformation of the simply supported beam to be measured at the point A and the connecting line of the two points A, B is the corner of the bending line at the point A. ThetaBThe included angle between the tangent line of the deformation of the simply supported beam to be measured at the point B and the connecting line of the two points A, B is the corner of the bending line at the point A. At the same time, thetaAAnd thetaBAnd also respectively represent both end angle values output from the filtering process, and Δ θ represents θ output from the filtering processAAnd thetaBThe difference between the angles of the two ends. And L represents the length of the simply supported beam to be measured, C represents the projection position of a resultant force point on the beam surface on the x axis, D represents the projection position of the maximum displacement point of the beam surface on the x axis, a represents the distance from the projection of the resultant force point on the x axis to the starting point of the simply supported beam to be measured, and b represents the distance from the projection of the resultant force point on the x axis to the end point of the beam. b. The relational expression of a and L is b ═ L-a.
X of the flexible line w (x) is an operand. E is the modulus of elasticity coefficient of the material, I is the cross-sectional moment of inertia of the beam, and the product of the two is the bending stiffness EI.
Wherein, the angle value ratio k is calculated by formula, k is greater than 0.5, and k is less than 2:
Figure GDA0003117860110000061
calculating the beam axis position a of the resultant force:
Figure GDA0003117860110000062
magnitude of resultant force FpThe calculation formula (c) is as follows:
Figure GDA0003117860110000063
calculating the beam axis position d of the maximum deflection:
Figure GDA0003117860110000071
maximum deflection value wmaxThe calculation formula of (A) is as follows:
when a is more than or equal to b:
Figure GDA0003117860110000072
when a < b:
Figure GDA0003117860110000073
substituting the simple beam to be detected into a plane rectangular coordinate system xow to form a simple beam deformation detection mathematical model, and after determining the position a of a resultant force point vertical to the beam surface and facing downwards on the x axis, having a flexible line equation:
when a is more than or equal to b:
Figure GDA0003117860110000074
when a < b:
Figure GDA0003117860110000075
(5) and (3) displaying the condition of simulating deformation by using animation on a screen by adopting the device in the step (1) according to the mathematical model for detecting deformation of the simply supported beam in the step (4), wherein the animation drawing method comprises the following steps:
and drawing the deformation simulation animation of the simply supported beam to be detected, and observing the stress deformation condition of the beam at the embedded end and the computer end. And the real-time angle value transmitted back by the angle sensor of the deformation simulation animation is output after filtering processing, and is substituted into the deformation detection mathematical model, and finally the calculated result is converted into a drawing point, and the real-time simulation animation of the stressed deformation axis of the simply supported beam to be detected is drawn on a screen.
The animation is displayed on a screen, the screen including: display module, computer display.
The method comprises the following steps of: 1) two coordinate systems are established: the beam coordinate system xoy is coincident with a drawing coordinate system x 'o' y ', the point o of the beam coordinate system is coincident with the point o' of the drawing coordinate system, the beam coordinate system is used for representing coordinate system parameters used in beam calculation, and the drawing coordinate system is used for representing a drawing coordinate point used for drawing a graph on a screen; 2) the number of animation basic drawing points is determined. Calculating the display coordinate of each drawing point on the screen through a beam deflection line equation; 3) and performing segmented drawing on each drawing point of the simulated animation.
Displaying a simulated animation curve calculated by a flexible line equation on a screen requires converting the calculation result into coordinate points on the screen, and the conversion mainly uses the following important variables:
simulating the number DRAW _ DOP of animation curve drawing points, wherein the number is not less than 3;
drawing the x' axis length graph _ size, wherein the length is not less than 30 and is an integral multiple of the drawing point number DRAW _ DOP of the simulated animation curve;
the number of x-axis component points in the beam coordinate system beam _ x [ ], and the parenthesis is 0 to DRAW _ DOP integers of DRAW _ DOP-1; the beam x array stores the x-axis coordinate of the points uniformly distributed on the beam axis in the beam coordinate system, and the difference value of each array element and the previous element is one-half of the DRAW DOP of the beam length.
Drawing coordinate system x' axis component point number group DRAW _ x [ ], wherein the parenthesis is 0 to DRAW _ DOP integers of DRAW _ DOP-1;
drawing a y-axis scaling ratio draw _ scale _ y of the coordinate system and the beam coordinate system;
storing an array lcd _ y [ ] of coordinate values of the drawing points y' axis, wherein the number of integers DRAW _ DOP from 0 to DRAW _ DOP-1 is included in brackets;
the graph _ size variable is the x' axis length of the animation drawing area on the screen, the beam _ x [ ] arrays store the x-axis coordinates of points uniformly distributed on the beam axis in the beam coordinate system, and the difference value of each array element and the previous element is one of DRAW _ DOP of the beam length. The elements of the array beam _ x [ ] are calculated from the origin, and the corresponding formula of each element and the last element is beam _ x [ i ] ═ length/DRAW _ DOP + beam _ x [ i-1 ];
wherein, length is the length of the beam, and the value interval of i is [1, DRAW _ DOP-1 ]. In the calculation, elements in beam _ x are substituted into a DRAW _ x array one by one to store x ' axis coordinate points uniformly distributed in a drawing area, namely, the x ' axis length of the drawing area is uniformly divided into DRAW _ DOP parts, and each division point corresponds to a coordinate point on an x ' axis. The elements of the array of DRAW _ x [ ] are calculated from the origin of coordinates, and the corresponding formula of each element to the last element is DRAW _ x [ i ] ═ graph _ size/DRAW _ DOP + DRAW _ x [ i-1 ];
the draw _ scale _ y variable is a y-axis scaling ratio of the drawing coordinate system to the beam coordinate system, and is used for displaying the drawing coordinate points of the deflection of the beam on the drawing coordinate system in the beam coordinate system in a proportional mode. The formula for the calculation of draw _ scale _ y is:
draw_scale_y＝graph_size/length；
and the lcd _ y [ ] array stores coordinate values of drawing points y' in a drawing coordinate system, the lcd _ y [ ] array has DRAW _ DOP elements in total, and each element and the corresponding deflection value in the beam coordinate axis have a relationship of multiple DRAW _ scale _ y. The calculation formula of the ith element of the lcd _ y [ ] array is lcd _ y [ i ] ═ band (beam _ x [ i ]) draw _ scale _ y; bend () is the flexible line equation w (x), and bend (beam _ x [ i ]) is the i-th element flexible line equation of the beam _ x [ i ] array.
According to the application and calculation of the parameters, drawing point coordinates in a drawing coordinate system can be obtained, and a sectional line drawing operation is carried out, so that a simulated animation curve of the beam axis can be obtained, wherein the sectional line drawing mode is as follows: and connecting every two adjacent drawing points by adopting a straight line segment. The process of drawing the deformation simulation animation curve is as follows:
h1, clearing the drawing area graph;
h2, calculating coordinate values of the drawing points on the screen one by one;
h3, segmented line drawing operation.
The invention provides a device and a method for monitoring stress and deformation of a simply supported beam-arm model in real time, which can be used for constructing bridges such as slab bridges, ribbed slab bridges, box girder bridges and the like, and mechanical analysis models of container suspension arms, and can also be used for monitoring stress bending deformation of solid beam arms in real time. The defects of strain gauges, computer vision and other measurement methods are overcome.
Drawings
FIG. 1 is a block diagram of an apparatus according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of the rectangular plane coordinate system xow and parameter labeling according to the present invention.
The labels in the figure are: 1-an embedded terminal; 2-a computer terminal; 10-a microcontroller module; 11-a wireless data transmission module; 12-a display module; 13-a wireless angle sensing module; 20-a data transfer relay; 21-computer.
Detailed Description
Example (b):
the method for establishing the real-time monitoring device for the bending deformation of the simply supported beam comprises the following steps: an embedded terminal 1 and a computer terminal 2; embedded terminal 1, comprising: the device comprises a wireless angle sensing module 13, a microcontroller module 10, a wireless data transmission module 11 and a display module 12. The microcontroller module 10 is electrically connected with the wireless angle sensing module 13, the wireless data transmission module 11 and the display module 12.
The embedded terminal 1 is used for acquiring sensor data, analyzing data, filtering data, converting data, displaying operation results in the form of graphs and characters, monitoring parameters in real time and alarming threshold values, wirelessly communicating with the computer terminal 2, sending filtered data to the computer terminal 2 and receiving a control instruction of the computer terminal 2;
the embedded end 1 is fixed on the simply supported beam to be detected and used for monitoring the deformation of the simply supported beam to be detected in real time;
the wireless angle sensing module 13 at least includes two three-axis angle sensors, and each angle sensor can provide instantaneous angle changes of three axes in a three-axis rectangular coordinate system. The angle sensors are placed at the positions of the two supporting points of the simply supported beam to be monitored, the angle values of the two supporting points are monitored in real time and transmitted back, and a wireless data transmission function is provided.
And the microcontroller module 10 is used for data analysis, data filtering, data conversion operation, real-time parameter monitoring and threshold value alarming, and outputting operation results in the form of graphs and characters to the display module 12.
And the wireless data transmission module 11 is used for carrying out wireless communication with the wireless angle sensing module 13, carrying out wireless communication with the computer end 2 and carrying out wired communication with the microcontroller module 10.
The display module 12 displays the deformation process and the result in the form of graphics and characters.
The computer terminal 2 includes: a data transfer relay 20; a computer 21 and upper computer software.
And the computer end 2 is used for remotely monitoring real-time data, providing a data interface of software, carrying out wireless communication with the embedded end 1, receiving the data processed by the embedded end 1 and sending a control instruction to the embedded end 1. And the upper computer software of the computer end 2 realizes the purposes of receiving data transmitted back by the serial port, sending serial port coding instructions, displaying the sampled real-time angle value, maximum deflection position, maximum deflection value, stress point and stress size on a graphical interface and simulating the real-time model animation of the stress deformation of the beam.
A data transfer relay 20 for: 1) receiving data returned by the embedded terminal 1 in a wireless mode, sending the data to the computer 21 by using a serial port, and processing the data by using upper computer software; 2) and receiving a control instruction sent by the computer 21 through the serial port, and sending the control instruction to the embedded terminal 1 in a wireless communication mode.
The embedded terminal 1 provides three functional modes: a parameter setting mode, a calibration function mode and a measurement function mode.
A parameter setting mode, the mode functions including: parameter setting, threshold setting, filtering enabling, wireless data transmission enabling, wireless channel detection and angle sensor detection.
And setting parameters, namely setting the length of the beam and the bending rigidity of the beam.
Setting a threshold value, namely setting a pressure threshold value point and a maximum deflection value threshold value point, and selecting whether to start a threshold value alarm function; the threshold alarm function is started in the measurement function mode, and if the pressure value and the maximum deflection value are monitored to exceed the set threshold point, the device gives an alarm and displays an alarm mark on the display module 12.
And enabling filtering, and selectively turning on or off the filtering function.
The wireless data enable selects to turn on or off the function of transmitting data back to the computer terminal 2 through the wireless mode.
And detecting a wireless channel, namely detecting whether the wireless angle sensing module 13 exists or not and whether the function can be normally operated or not.
The angle sensor detects whether there is data returned by the angle sensor.
A calibration function mode, the mode functions being: calibrating an angle zero point, correcting angle deviation and calibrating the bending rigidity of the beam to be measured.
And calibrating an angle zero point, and sampling an angle value of the simply supported beam to be tested in a static state under the static condition that the simply supported beam to be tested is not subjected to external pressure, and calibrating the angle value into an angle zero point value.
And correcting the angle deviation for calibrating the simply supported beam model. The calibration process comprises the following steps: and (3) applying downward pressure at the center of the beam by taking the two supporting points as the end points of the beam, and solving and storing the difference between the angle values of the two ends after the angle values of the two ends are stable.
A measurement function mode, the functions of which are: the angle sensor data acquisition, data analysis, data filtering, data conversion operation, operation result display in the form of graph and numerical characters, real-time monitoring of angle change and threshold value alarm, wireless communication with the computer end 2, transmission of filtered data with the computer end 2, and reception of control instructions of the computer end 2.
(1) The device established by the embodiment is used for deformation detection of the simply supported beam to be detected, and filtering processing needs to be carried out on angle data acquired by the wireless sensing module. The filtering processing method comprises the following steps:
in the filtering process, the angle value is divided into three stages for processing, which are respectively: angle bare value, internal division number and angle view. In total, three steps can be divided: denoising processing, zero point processing and division number processing.
The angle bare value is an angle value which is sampled by the angle sensor and fails to be processed, namely angle most original data.
The internal part degree is an angle calculated quantity obtained after the angle bare value is subjected to denoising treatment and zero point treatment.
The angle view is the final filtering output quantity after the division number processing is carried out on the internal division number, and is used as the angle value output by the filtering processing.
Denoising processing for filtering out fluctuation causing instability of the sampling value, the denoising processing comprising: and performing zero tracking judgment, moving average filtering and weighted filtering on the angle bare value. The denoising processing steps are as follows:
q1, inputting angle bare value;
q2, zero-tracking range judgment processing;
q3, moving average filtering;
q4, weighted filtering;
q5, processing the angle change range;
q6, outputting an angle filtered value.
And zero tracking range judgment processing for judging whether the current angle bare value is in the zero tracking range. If the zero tracking range is within the zero tracking range, zero tracking processing is performed by zero point correction, and a zero tracking output value is output. If the angle exceeds the zero tracking range, the current angle bare value is cut off.
And (5) tracking a zero range, and taking the maximum change range of the angle value acquired under the static state that the simply supported beam to be detected is not subjected to the external force.
The zero-tracking process is to increase or decrease the angle value so that the angle value approaches 0 when the angle value is within a predetermined range.
And (3) moving average filtering, wherein a moving average filtering array is used for storing zero tracking output values as a filtering operation reference quantity, the size of the array is N, and N is a natural number which is greater than or equal to 10. And (3) carrying out filtering averaging by adopting N reference quantities, and putting a new zero-tracking output value into the first element of the array every time of carrying out moving average filtering.
For the moving average filter array, two filling operation modes are provided: fast refresh fill operation, slow refresh fill operation.
And (3) a fast refreshing filling operation mode, namely, putting the first three zero-tracing output values in the array into the last three bits of the array, and filling the rest array elements with the zero-tracing output values output by the last filtering.
And (4) slowly refreshing the filling operation mode, and transmitting the zero-chasing output values in the array one by one.
And (3) moving average filtering, combining two filling operation modes, firstly running M times of slow refreshing filling operation in a filtering processing period, and then executing Q times of fast refreshing filling operation. Wherein M is a natural number greater than or equal to 2, and Q is a natural number greater than or equal to 2. After each fill, the N elements in the array are averaged and the average is output to a weighted filter.
Weighting filtering, which is used for endowing the angle bare value with different weight values in different change processes;
weighting and filtering, wherein if the simply supported beam component to be measured is in a stable state, the newly acquired angle bare value weight is weakened; and if the beam member is in an unstable state, reinforcing the newly acquired angle bare value weight.
And (4) weighted filtering, namely judging whether the stable state of the internal division number is achieved or not.
The stable state of the internal division number is that the difference between the newly input angle bare value and the last angle visual number output by filtering does not exceed a set range.
If the stable state of the internal part number is reached, the filtering calculation formula is as follows: the filter output value (newly input angle bare value 1+ last filter output angle view 3)/4; in the above equation, the weight of the angle bare value of the new input is 25%.
If the steady state of the internal division number is not reached, the filter calculation formula is as follows: the filter output value (newly input angle bare value x 3+ last filter output angle view x 1)/4; in the above equation, the weight of the angle bare value newly input is 75%.
And the angle change range processing function is used for comparing the newly input angle bare value with the last filtering output value, judging whether the angle value is stable according to the change range interval, and determining the filtering times required when the internal division number reaches a stable state.
Zero point processing is used for correcting angle zero point deviation caused by the simply supported beam to be measured in the working process, and measuring deviation caused by the following conditions is solved: 1. the shape of the simply supported beam to be measured changes when the simply supported beam returns to the static state after being stressed and deformed every time, so that the angle zero value of the simply supported beam to be measured acquired by the angle sensor in the static state can deviate from the calibration zero value in the calibration process; 2. the mechanical loss of the simply supported beam to be tested in a stressed state for a long time causes the deviation between the actual two-end angle zero value of the simply supported beam to be tested in a static state and the calibration zero value in a calibration state to be gradually increased.
Zero point processing, comprising: zero point correction and automatic zero tracking.
Zero point correction, comprising three functions: correcting a calibration zero difference value, correcting an accumulated zero value and filtering an angle difference value at two ends.
And correcting the calibration zero point, and subtracting the calibration zero point value during calibration from the angle zero point value acquired each time.
And correcting the accumulated zero point, and calculating the zero tracking amount generated after zero tracking into the angle bare value to change the angle zero point value of the simply supported beam to be measured in a static state.
And filtering the angle difference value of the two ends, and detecting the bending deformation of the simply supported beam to be detected by calculating the relative change value of the angle values of the two ends, so that the simply supported beam to be detected is not influenced by the inclination of the simply supported beam.
And the automatic zero tracking is used for carrying out zero tracking calculation on the current angle value within the zero tracking range according to the set time frequency. If the current angle value is larger than the zero point value, the zero chasing amount is automatically increased, otherwise, the zero chasing amount is automatically decreased. When the change of the current angle value is within the set range, the change is stable.
The automatic zero-tracking process is as follows:
z1, judging whether the zero-chasing tracking time is reached, if so, entering Z2, otherwise, jumping out of the automatic zero-chasing program;
z2, judging whether the current angle value is stable, if so, entering Z3, otherwise, jumping out of the automatic zero-tracing program;
z3, judging whether the current angle value is less than or equal to zero-chasing range, if so, entering Z4, otherwise, jumping out of the automatic zero-chasing program;
z4, judging whether the current zero point is larger than the minimum difference value, if so, entering Z6, and otherwise, entering Z5;
z5, catching up the end value and completing an automatic zero catching program;
z6, judging whether the current angle value is positive, if so, entering Z8, and if not, entering Z7;
z7, the zero-chasing quantity is reduced, and the automatic zero-chasing procedure is completed;
and Z8, the zero chasing amount is automatically increased, and the automatic zero chasing program is completed.
The last value is traced for zeroing the current angle value.
After the zero point processing is executed, converting the angle value after the zero point processing into an internal division number, wherein the formula is as follows: CntIn _ Cur is Angle/32768 180 °, where Angle represents the Angle value after filtering and CntIn _ Cur represents the internal index number.
And index number processing is used for avoiding unstable measurement data caused by the phenomenon of stress deformation and jitter of the simply supported beam to be measured.
And (4) index number processing, namely performing filtering operation by judging the variation range of the index value and according to the variation range.
The differential value of the output value is processed by the fractal value, the internal division number and the division number.
Index number processing, setting 3 index number change critical values: 3C, 6C and 10C. C is the basic quantity of the set index number change, wherein C is a positive integer multiple of 0.5 degrees, the 3C critical point represents the value of 3 times the basic quantity, the 6C critical point represents the value of 6 times the basic quantity, and the 10C critical point represents the value of 10 times the basic quantity.
The division number processing program flow is as follows:
n1, judging whether the visual value is less than 3C, if so, jumping out of the index number processing program, otherwise, entering N2;
n2, judging whether the score value is less than 6C, if so, entering N2-1, otherwise, entering N3;
n2-1, accumulating the internal division number for 8 times, and obtaining the average value of 8 times as the output of division number processing;
n3, judging whether the score value is less than 10C, if so, entering N3-1, otherwise, entering N3-2;
n3-1, accumulating the internal division number for 15 times, and obtaining the average value of the 15 times as the output of division number processing;
n3-2, the number of internal divisions is added up 3 times, and the average value of the 3 times is obtained as the output of division number processing.
(2) In this embodiment, a mathematical model for detecting deformation of a simply supported beam is established, and the mathematical model for detecting deformation of a simply supported beam is as follows:
the key parameters that need to be calculated are: resultant force F acting on beam surface vertically and downwardspResultant force FpThe beam axis position a and the maximum deflection value wmaxAnd the beam axis position d where the maximum deflection is located, and the calculated value needs to be displayed in real time through a screen and alarm is realized.
A plane rectangular coordinate system xow is constructed, the plane rectangular coordinate system xow uses the plane of the beam axis, a pivot point a of the beam as an origin O, and the horizontal extension direction of the beam as the positive half axis direction of the X axis to establish a plane rectangular coordinate system XOW, and the plane rectangular coordinate system XOW has the X axis and the W axis which are perpendicular to each other.
And the fitting calculation equation of the bending deformation shape of the simply supported beam to be detected on the stress plane is the bending line w (x) of the beam.
The resultant force received by the deformation of the simply supported beam to be measured is a downward resultant force F vertical to the beam surfacep. The two supporting points of the simply supported beam to be measured are point A and point B, the length of the beam is A, B straight line distance L between the two supporting points, the connecting line of the two supporting points A, B is the axis before beam deformation, the axis of the beam is on the X axis of the plane coordinate system XOW,
θAthe included angle between the tangent line of the deformation of the simply supported beam to be measured at the point A and the connecting line of the two points A, B is the corner of the bending line at the point A. ThetaBThe included angle between the tangent line of the deformation of the simply supported beam to be measured at the point B and the connecting line of the two points A, B is the corner of the bending line at the point A. At the same time, thetaAAnd thetaBAnd also respectively represent both end angle values output from the filtering process, and Δ θ represents θ output from the filtering processAAnd thetaBThe difference between the angles of the two ends. And L represents the length of the simply supported beam to be measured, C represents the projection position of a resultant force point on the beam surface on the x axis, D represents the projection position of the maximum displacement point of the beam surface on the x axis, a represents the distance from the projection of the resultant force point on the x axis to the starting point of the simply supported beam to be measured, and b represents the distance from the projection of the resultant force point on the x axis to the end point of the beam. b. The relational expression of a and L is b ═ L-a.
X of the deflection line w (x) is an operation variable, E is an elastic modulus coefficient of the material, I is a section moment of inertia of the beam, and the product of the two is the bending stiffness EI.
Wherein, the angle value ratio k is calculated by formula, k is greater than 0.5, and k is less than 2:
Figure GDA0003117860110000141
calculating the beam axis position a of the resultant force:
Figure GDA0003117860110000142
magnitude of resultant force FpThe calculation formula (c) is as follows:
Figure GDA0003117860110000143
calculating the beam axis position d of the large deflection:
Figure GDA0003117860110000151
maximum deflection value wmaxThe calculation formula of (A) is as follows:
when a is more than or equal to b:
Figure GDA0003117860110000152
when a < b:
Figure GDA0003117860110000153
substituting the simple beam to be detected into a plane rectangular coordinate system xow to form a simple beam deformation detection mathematical model, and after determining the position a of a resultant force point vertical to the beam surface and facing downwards on the x axis, having a flexible line equation:
when a is more than or equal to b:
Figure GDA0003117860110000154
when a < b:
Figure GDA0003117860110000155
(3) by adopting the device, according to the mathematical model, the deformation condition is simulated by the animation displayed on the screen, and the animation drawing method comprises the following steps:
and drawing a deformation simulation animation of the simply supported beam to be detected, and observing the stress deformation condition of the beam at the embedded end 1 and the computer end 2. And the real-time angle value transmitted back by the angle sensor of the deformation simulation animation is output after filtering processing, and is substituted into the deformation detection mathematical model, and finally the calculated result is converted into a drawing point, and the real-time simulation animation of the stressed deformation axis of the simply supported beam to be detected is drawn on a screen.
The animation is displayed on a screen, the screen including: a display module 12 or a display of a computer 20.
The method comprises the following steps of: 1) two coordinate systems are established: the beam coordinate system xoy is coincident with a drawing coordinate system x 'o' y ', the point o of the beam coordinate system is coincident with the point o' of the drawing coordinate system, the beam coordinate system is used for representing coordinate system parameters used in beam calculation, and the drawing coordinate system is used for representing a drawing coordinate point used for drawing a graph on a screen; 2) the number of animation basic drawing points is determined. Calculating the display coordinate of each drawing point on the screen through a beam deflection line equation; 3) and performing segmented drawing on each drawing point of the simulated animation.
Displaying a simulated animation curve calculated by a flexible line equation on a screen requires converting the calculation result into coordinate points on the screen, and the conversion mainly uses the following important variables:
simulating the number DRAW _ DOP of animation curve drawing points, wherein the number is not less than 3;
drawing the x' axis length graph _ size, wherein the length is not less than 30 and is an integral multiple of the drawing point number DRAW _ DOP of the simulated animation curve;
the number of x-axis component points in the beam coordinate system beam _ x [ ], and the parenthesis is 0 to DRAW _ DOP integers of DRAW _ DOP-1; the beam x array stores the x-axis coordinate of the points uniformly distributed on the beam axis in the beam coordinate system, and the difference value of each array element and the previous element is one-half of the DRAW DOP of the beam length.
Drawing coordinate system x' axis component point number group DRAW _ x [ ], wherein the parenthesis is 0 to DRAW _ DOP integers of DRAW _ DOP-1;
drawing a y-axis scaling ratio draw _ scale _ y of the coordinate system and the beam coordinate system;
storing an array lcd _ y [ ] of coordinate values of the drawing points y' axis, wherein the number of integers DRAW _ DOP from 0 to DRAW _ DOP-1 is included in brackets;
the graph _ size variable is the x' axis length of the animation drawing area on the screen, the beam _ x [ ] arrays store the x-axis coordinates of points uniformly distributed on the beam axis in the beam coordinate system, and the difference value of each array element and the previous element is one of DRAW _ DOP of the beam length. The elements of the array beam _ x [ ] are calculated from the origin, and the corresponding formula of each element and the last element is beam _ x [ i ] ═ length/DRAW _ DOP + beam _ x [ i-1 ];
wherein, length is the length of the beam, and the value interval of i is [1, DRAW _ DOP-1 ]. In the calculation, elements in beam _ x are substituted into a DRAW _ x array one by one to store x ' axis coordinate points uniformly distributed in a drawing area, namely, the x ' axis length of the drawing area is uniformly divided into DRAW _ DOP parts, and each division point corresponds to a coordinate point on an x ' axis. The elements of the array of DRAW _ x [ ] are calculated from the origin of coordinates, and the corresponding formula of each element to the last element is DRAW _ x [ i ] ═ graph _ size/DRAW _ DOP + DRAW _ x [ i-1 ];
the draw _ scale _ y variable is a y-axis scaling ratio of the drawing coordinate system to the beam coordinate system, and is used for displaying the drawing coordinate points of the deflection of the beam on the drawing coordinate system in the beam coordinate system in a proportional mode. The formula for the calculation of draw _ scale _ y is:
draw_scale_y＝graph_size/length；
and the lcd _ y [ ] array stores coordinate values of drawing points y' in a drawing coordinate system, the lcd _ y [ ] array has DRAW _ DOP elements in total, and each element and the corresponding deflection value in the beam coordinate axis have a relationship of multiple DRAW _ scale _ y. The calculation formula of the ith element of the lcd _ y [ ] array is lcd _ y [ i ] ═ band (beam _ x [ i ]) draw _ scale _ y; bend () is the flexible line equation w (x), and bend (beam _ x [ i ]) is the i-th element flexible line equation of the beam _ x [ i ] array.
According to the application and calculation of the parameters, drawing point coordinates in a drawing coordinate system can be obtained, and a sectional line drawing operation is carried out, so that a simulated animation curve of the beam axis can be obtained, wherein the sectional line drawing mode is as follows: and connecting every two adjacent drawing points by adopting a straight line segment. The process of drawing the deformation simulation animation curve is as follows:
h1, clearing the drawing area graph;
h2, calculating coordinate values of the drawing points on the screen one by one;
h3, segmented line drawing operation.

Claims (1)

1. A method for monitoring bending deformation of a simply supported beam in real time is characterized by comprising the following specific steps:
(1) the method for establishing the real-time monitoring device for the bending deformation of the simply supported beam comprises the following steps: an embedded terminal and a computer terminal; an embedded port, comprising: the device comprises a wireless angle sensing module, a microcontroller module, a wireless data transmission module and a display module; the microcontroller module is electrically connected with the wireless angle sensing module, the wireless data transmission module and the display module;
the embedded end is used for acquiring sensor data, analyzing the data, filtering the data, converting the data, displaying the operation results in the form of graphs and characters, monitoring parameters in real time and alarming a threshold value, wirelessly communicating with the computer end, sending the filtered data to the computer end and receiving a control instruction of the computer end;
the embedded end is fixed on the simply supported beam to be detected and used for monitoring the deformation of the simply supported beam to be detected in real time;
the wireless angle sensing module at least comprises two triaxial angle sensors, and each angle sensor can provide instantaneous angle change of three axes in a triaxial rectangular coordinate system; the angle sensors are placed at the two fulcrum positions of the simply supported beam to be monitored, the angle values of the two fulcrums are monitored in real time and are transmitted back, and a wireless data transmission function is provided;
the microcontroller module is used for data analysis, data filtering, data conversion operation, real-time parameter monitoring and threshold value alarming, and outputting operation results in the form of graphs and characters to the display module;
the wireless data transmission module is used for carrying out wireless communication with the wireless angle sensing module, carrying out wireless communication with a computer end and carrying out wired communication with the microcontroller module;
the display module displays the deformation process and the result in the form of graphs and characters;
computer end, including two parts, the first part is: a data transmission transfer device; the first part is: computer and upper computer software;
the computer end is used for remotely monitoring real-time data, providing a data interface of software, carrying out wireless communication with the embedded end, receiving the data processed by the embedded end and sending a control instruction to the embedded end; the computer-side upper computer software realizes the purposes of receiving data transmitted back by the serial port, sending serial port coding instructions, displaying the sampled real-time angle value, maximum deflection position, maximum deflection value, stress point and stress size on a graphical interface and simulating the real-time model animation of the stress deformation of the beam;
a data transfer relay for: receiving data returned by the embedded terminal in a wireless mode, sending the data to a computer by using a serial port, and processing the data by using upper computer software; receiving a control instruction sent by a computer through a serial port, and sending the control instruction to an embedded terminal in a wireless communication mode;
(2) the embedded terminal in the step (1) provides three functional modes: a parameter setting mode, a calibration function mode and a measurement function mode;
a parameter setting mode, the mode functions including: parameter setting, threshold setting, filtering enabling, wireless data transmission enabling, wireless channel detection and angle sensor detection;
setting parameters, namely setting the length of a beam and the bending rigidity of the beam;
setting a threshold value, namely setting a pressure threshold value point and a maximum deflection value threshold value point, and selecting whether to start a threshold value alarm function; starting a threshold alarm function in a measurement function mode, and if the pressure value and the maximum deflection value are monitored to exceed a set threshold point, sending an alarm by the device and displaying an alarm mark on a display module;
enabling filtering, and selecting to turn on or off a filtering function;
enabling wireless data, and selecting to turn on or turn off the function of transmitting data back to the computer terminal in a wireless mode;
detecting a wireless channel, namely detecting whether a wireless module exists or not and whether the function can be normally operated or not;
detecting by an angle sensor, and detecting whether data returned by the angle sensor exists;
a calibration function mode, the mode functions being: calibrating an angle zero point, correcting angle deviation and calibrating the bending rigidity of the beam to be measured;
calibrating an angle zero point, and sampling an angle value of the simply supported beam to be tested in a static state under the condition that the simply supported beam to be tested is not subjected to external pressure, and calibrating the angle value as an angle zero point value;
correcting the angle deviation for calibrating the simply supported beam model; the calibration process comprises the following steps: applying downward pressure to the center of the beam by taking the two supporting points as the end points of the beam, and solving and storing the difference value between the angle values of the two ends after the angle values of the two ends are stable;
a measurement function mode, the functions of which are: the angle sensor data acquisition, data analysis, data filtering, data conversion operation, operation result display in the form of graph and numerical characters, real-time monitoring of angle change and threshold value alarm, wireless communication with a computer end, transmission of filtered data with the computer end, and reception of control instructions of the computer end;
(3) the real-time monitoring device for the bending deformation of the simply supported beam in the step (1) carries out deformation detection on the simply supported beam to be detected, and filtering processing needs to be carried out on angle data acquired by the wireless sensing module; the filtering processing method comprises the following steps:
in the filtering process, the angle value is divided into three stages for processing, which are respectively: angle bare value, internal division number and angle view number; in total, three steps can be divided: denoising, zero point processing and division number processing;
the angle bare value is an angle value which is sampled by the angle sensor and is not processed, namely the angle most original data;
the internal part degree is an angle calculated quantity obtained after the angle bare value is subjected to denoising treatment and zero point treatment;
the angle view number is the final filtering output quantity of the internal division number after division number processing and is used as the angle value output by filtering processing;
denoising processing for filtering out fluctuation causing instability of the sampling value, the denoising processing comprising: zero-tracing judgment, moving average filtering and weighted filtering are carried out on the angle bare value; the denoising processing steps are as follows:
q1, inputting angle bare value;
q2, zero-tracking range judgment processing;
q3, moving average filtering;
q4, weighted filtering;
q5, processing the angle change range;
q6, outputting an angle filtered value;
zero-tracing range judgment processing is used for judging whether the bare value of the current angle is in a zero-tracing range; if the zero tracking range is within the zero tracking range, zero tracking processing is carried out through zero point correction, and a zero tracking output value is output; if the angle exceeds the zero tracking range, the bare value of the current angle is cut off;
a zero tracking range is obtained, and the maximum change range of the angle value is obtained when the simply supported beam to be detected is in a static state without being subjected to external force;
zero-tracing processing, namely, when the angle value is in a specified range, increasing or decreasing the angle value to enable the angle value to approach 0;
the method comprises the following steps of (1) moving average filtering, wherein a moving average filtering array is used for storing zero tracking output values as a filtering operation reference quantity, the size of the array is N, and N is a natural number which is greater than or equal to 10; carrying out filtering averaging by adopting N reference quantities, and putting a new zero-tracing output value into a first element of an array every time moving average filtering is carried out;
for the moving average filter array, two filling operation modes are provided: fast refresh fill operation, slow refresh fill operation;
fast refreshing a filling operation mode, putting the first three zero-chasing output values in the array into the last three bits of the array, and filling the remaining array elements with the zero-chasing output values output by the last filtering;
slowly refreshing a filling operation mode, and transmitting zero-chasing output values in the array one by one;
the method comprises the following steps of (1) moving average filtering, combining two filling operation modes, firstly operating M times of slow refreshing filling operation in a filtering processing period, and then executing Q times of fast refreshing filling operation, wherein M is a natural number greater than or equal to 2, and Q is a natural number greater than or equal to 2; after each time of filling, averaging the N elements in the array, and outputting the average to weighted filtering;
weighting filtering, which is used for endowing the angle bare value with different weight values in different change processes;
weighting and filtering, wherein if the simply supported beam component to be measured is in a stable state, the newly acquired angle bare value weight is weakened; if the beam member is in an unstable state, reinforcing the weight of the newly acquired angle bare value;
weighted filtering, namely judging whether the stable state of the internal division number is achieved or not;
the stable state of the internal division number is that the difference between the newly input angle bare value and the last filtering output angle visual number does not exceed a set range;
if the stable state of the internal part number is reached, the filtering calculation formula is as follows: the filter output value (newly input angle bare value 1+ last filter output angle view 3)/4; in the above formula, the newly input angle bare value weight is 25%;
if the steady state of the internal division number is not reached, the filter calculation formula is as follows: the filter output value (newly input angle bare value x 3+ last filter output angle view x 1)/4; in the above formula, the weight of the newly input angle bare value is 75%;
the angle change range processing function is used for comparing a newly input angle bare value with a last filtering output value, judging whether the angle value is stable according to a change range interval, and determining the filtering times required when the internal division number reaches a stable state;
zero point processing is used for correcting angle zero point deviation caused by the simply supported beam to be measured in the working process, and measuring deviation caused by the following conditions is solved: 1) the shape of the simply supported beam to be measured changes when the simply supported beam returns to the static state after being stressed and deformed every time, so that the angle zero value of the simply supported beam to be measured acquired by the angle sensor in the static state can deviate from the calibration zero value in the calibration process; 2) the mechanical loss of the simply supported beam to be tested in a stressed state for a long time causes the deviation between the actual zero values of the two end angles of the simply supported beam to be tested in a static state and the calibration zero value in a calibration state to be gradually increased;
zero point processing, comprising: zero point correction and automatic zero tracing;
zero point correction, comprising three functions: correcting a calibration zero difference value, correcting an accumulated zero value and filtering an angle difference value at two ends;
correcting the calibration zero point, and subtracting the calibration zero point value during calibration from the angle zero point value acquired each time;
correcting the accumulated zero point, and calculating the zero tracking amount generated after zero tracking into an angle bare value to change the angle zero point value of the simply supported beam to be detected in a static state;
filtering the angle difference value of the two ends, and detecting the bending deformation of the simply supported beam to be detected by calculating the relative change value of the angle values of the two ends, so that the simply supported beam to be detected is not influenced by the inclination of the simply supported beam;
the automatic zero tracking is used for carrying out zero tracking operation on the current angle value within the zero tracking range according to the set time frequency; if the current angle value is larger than the zero value, the zero chasing amount is automatically increased, otherwise, the zero chasing amount is automatically decreased; when the change of the current angle value is within a set range, the change is stable;
the automatic zero-tracking process is as follows:
z1, judging whether the zero-chasing tracking time is reached, if so, entering Z2, otherwise, jumping out of the automatic zero-chasing program;
z2, judging whether the current angle value is stable, if so, entering Z3, otherwise, jumping out of the automatic zero-tracing program;
z3, judging whether the current angle value is less than or equal to zero-chasing range, if so, entering Z4, otherwise, jumping out of the automatic zero-chasing program;
z4, judging whether the current zero point is larger than the minimum difference value, if so, entering Z6, and otherwise, entering Z5;
z5, catching up the end value and completing an automatic zero catching program;
z6, judging whether the current angle value is positive, if so, entering Z8, and if not, entering Z7;
z7, the zero-chasing quantity is reduced, and the automatic zero-chasing procedure is completed;
z8, the zero-chasing quantity is automatically increased, and an automatic zero-chasing program is completed;
tracing the last value for returning the current angle value to zero;
after the zero point processing is executed, converting the angle value after the zero point processing into an internal division number, wherein the formula is as follows:
CntIn _ Cur ═ Angle/32768 × 180 °, where Angle represents the Angle value after filtering output, and CntIn _ Cur represents the internal division number;
the division number processing is used for avoiding unstable measurement data caused by the phenomenon of stress deformation and jitter of the simply supported beam to be measured;
the division number processing is carried out, and filtering operation is carried out according to the change range of the division value by judging the change range of the division value;
the differential value of the output value of the internal division number and the division number processing is obtained;
index number processing, setting 3 index number change critical values: 3C, 6C and 10C; c is a basic quantity of the set index number change, wherein C is a positive integer multiple of 0.5 degrees; the 3C critical point represents the value of 3 times of the basic quantity, the 6C critical point represents the value of 6 times of the basic quantity, and the 10C critical point represents the value of 10 times of the basic quantity;
the division number processing program flow is as follows:
n1, judging whether the visual value is less than 3C, if so, jumping out of the index number processing program, otherwise, entering N2;
n2, judging whether the score value is less than 6C, if so, entering N2-1, otherwise, entering N3;
n2-1, accumulating the internal division number for 8 times, and obtaining the average value of 8 times as the output of division number processing;
n3, judging whether the score value is less than 10C, if so, entering N3-1, otherwise, entering N3-2;
n3-1, accumulating the internal division number for 15 times, and obtaining the average value of the 15 times as the output of division number processing;
n3-2, accumulating the internal division number for 3 times, and obtaining the average value of the 3 times as the output of division number processing;
(4) establishing a mathematical model for detecting deformation of the simply supported beam, wherein the mathematical model for detecting deformation of the simply supported beam is as follows:
the key parameters that need to be calculated are: the magnitude of the vertical downward resultant force F acting on the beam surface of the simply supported beam to be measuredpResultant force FpThe position a of the beam axis and the maximum deflection valuewmaxAnd the beam axis position d where the maximum deflection is located, and the calculated value needs to be displayed in real time through a screen and alarm is realized;
a plane rectangular coordinate system xow is constructed, the plane rectangular coordinate system xow takes the plane of the beam axis as the origin O, takes a pivot point A of the beam as the origin O, and takes the horizontal extension direction of the beam as the positive semi-axis direction of the X axis to establish a plane rectangular coordinate system XOW, and the plane rectangular coordinate system XOW has the X axis and the W axis which are perpendicular to each other;
fitting a calculation equation of the bending deformation shape of the simply supported beam to be detected on a stress plane, wherein the equation is a bending line w (x) of the beam;
the resultant force received by the deformation of the simply supported beam to be measured is a downward resultant force F vertical to the beam surfacep(ii) a The two supporting points of the simply supported beam to be measured are point A and point B, the length of the beam is A, B straight line distance L between the two supporting points, the connecting line of the two supporting points A, B is the axis before beam deformation, the axis of the beam is on the X axis of the plane coordinate system XOW,
θAthe included angle between the tangent line of the deformation of the simply supported beam to be detected at the point A and the connecting line of the two points A, B is the corner of the bending line at the point A; thetaBThe included angle between the tangent line of the deformation of the simply supported beam to be detected at the point B and the connecting line of the two points A, B is the corner of the bending line at the point A; at the same time, thetaAAnd thetaBAnd also respectively represent both end angle values output from the filtering process, and Δ θ represents θ output from the filtering processAAnd thetaBThe difference of the angles of the two ends; l represents the length of the simply supported beam to be measured, C represents the projection position of a resultant force point on the beam surface on the x axis, D represents the projection position of the maximum displacement point of the beam surface on the x axis, a represents the distance from the projection of the resultant force point on the x axis to the starting point of the simply supported beam to be measured, and b represents the distance from the projection of the resultant force point on the x axis to the end point of the beam; b. the relational expression of a and L is b ═ L-a;
x of the flexible line w (x) is an operation variable; e is the elastic modulus coefficient of the material, I is the cross-sectional moment of inertia of the beam, and the product of the two is the bending stiffness EI;
wherein, the angle value ratio k is calculated by formula, k is greater than 0.5, and k is less than 2:
Figure FDA0003117860100000061
calculating the beam axis position a of the resultant force:
Figure FDA0003117860100000062
magnitude of resultant force FpThe calculation formula (c) is as follows:
Figure FDA0003117860100000063
calculating the beam axis position d of the maximum deflection:
Figure FDA0003117860100000064
maximum deflection value wmaxThe calculation formula of (A) is as follows:
when a is more than or equal to b:
Figure FDA0003117860100000065
when a < b:
Figure FDA0003117860100000066
substituting the simple beam to be detected into a plane rectangular coordinate system xow to form a simple beam deformation detection mathematical model, and after determining the position a of a resultant force point vertical to the beam surface and facing downwards on the x axis, having a flexible line equation:
when a is more than or equal to b:
Figure FDA0003117860100000067
when a < b:
Figure FDA0003117860100000068
(5) and (3) displaying the condition of simulating deformation by using animation on a screen by adopting the device in the step (1) according to the mathematical model for detecting deformation of the simply supported beam in the step (4), wherein the animation drawing method comprises the following steps:
drawing a deformation simulation animation of the simply supported beam to be detected, wherein the animation is used for observing the stress deformation condition of the beam at the embedded end and the computer end; the real-time angle value transmitted back by the angle sensor of the deformation simulation animation is output after filtering processing, and is substituted into the deformation detection mathematical model, and finally the calculated result is converted into a drawing point, and the real-time simulation animation of the stressed deformation axis of the simply supported beam to be detected is drawn on a screen;
the animation is displayed on a screen, the screen including: a display module and a display of a computer;
the method comprises the following steps of: 1) two coordinate systems are established: the beam coordinate system xoy is coincident with a drawing coordinate system x 'o' y ', the point o of the beam coordinate system is coincident with the point o' of the drawing coordinate system, the beam coordinate system is used for representing coordinate system parameters used in beam calculation, and the drawing coordinate system is used for representing a drawing coordinate point used for drawing a graph on a screen; 2) determining the number of basic drawing points of the animation; calculating the display coordinate of each drawing point on the screen through a beam deflection line equation; 3) each drawing point of the simulated animation is drawn in a segmented mode;
displaying a simulated animation curve calculated by a flexible line equation on a screen requires converting the calculation result into coordinate points on the screen, and the conversion mainly uses the following important variables:
simulating the number DRAW _ DOP of animation curve drawing points, wherein the number is not less than 3;
drawing the x' axis length graph _ size, wherein the length is not less than 30 and is an integral multiple of the drawing point number DRAW _ DOP of the simulated animation curve;
the number of x-axis component points in the beam coordinate system beam _ x [ ], and the parenthesis is 0 to DRAW _ DOP integers of DRAW _ DOP-1; the beam _ x [ ] arrays store the x-axis coordinates of points uniformly distributed on the beam axis in a beam coordinate system, and the difference value of each array element and the previous element is one-half of DRAW _ DOP of the beam length; drawing coordinate system x' axis component point number group DRAW _ x [ ], wherein the parenthesis is 0 to DRAW _ DOP integers of DRAW _ DOP-1;
drawing a y-axis scaling ratio draw _ scale _ y of the coordinate system and the beam coordinate system;
storing an array lcd _ y [ ] of coordinate values of the drawing points y' axis, wherein the number of integers DRAW _ DOP from 0 to DRAW _ DOP-1 is included in brackets;
the graph _ size variable is the x' axis length of an animation drawing area on the screen, the beam _ x [ ] arrays store the x-axis coordinates of points which are uniformly distributed on the beam axis in a beam coordinate system, and the difference value of each array element and the previous element is one of DRAW _ DOP of the beam length; the elements of the array beam _ x [ ] are calculated from the origin, and the corresponding formula of each element and the last element is beam _ x [ i ] ═ length/DRAW _ DOP + beam _ x [ i-1 ];
wherein, length is the length of the beam, and the value interval of i is [1, DRAW _ DOP-1 ]; in the calculation, elements in beam _ x are substituted into a DRAW _ x array one by one to store x ' axis coordinate points uniformly distributed in a drawing area, namely, the x ' axis length of the drawing area is uniformly divided into DRAW _ DOP parts, and each division point corresponds to a coordinate point on an x ' axis; the elements of the array of DRAW _ x [ ] are calculated from the origin of coordinates, and the corresponding formula of each element to the last element is DRAW _ x [ i ] ═ graph _ size/DRAW _ DOP + DRAW _ x [ i-1 ];
the draw _ scale _ y variable is a y-axis scaling ratio of the drawing coordinate system to the beam coordinate system, and is used for displaying drawing coordinate points of the deflection of the beam in the beam coordinate system on the drawing coordinate system in proportion; the formula for the calculation of draw _ scale _ y is:
draw_scale_y＝graph_size/length；
the lcd _ y [ ] array stores the coordinate value of the drawing point y' axis in the drawing coordinate system, the lcd _ y [ ] array has DRAW _ DOP elements, and each element and the corresponding deflection value in the beam coordinate axis have the relationship of multiple DRAW _ scale _ y; the calculation formula of the ith element of the lcd _ y [ ] array is lcd _ y [ i ] ═ band (beam _ x [ i ]) draw _ scale _ y; wherein, bend () is the established flexible line equation w (x) of (4), and bend (beam _ x [ i ]) is the i-th element flexible line equation of the beam _ x [ i ] array;
according to the application and calculation of the parameters, drawing point coordinates in a drawing coordinate system can be obtained, and a sectional line drawing operation is carried out, so that a simulated animation curve of the beam axis can be obtained, wherein the sectional line drawing mode is as follows: connecting every two adjacent drawing points by adopting a straight line segment; the process of drawing the deformation simulation animation curve is as follows:
h1, clearing the drawing area graph;
h2, calculating coordinate values of the drawing points on the screen one by one;
h3, segmented line drawing operation.
CN202010090290.9A 2020-02-13 2020-02-13 Real-time monitoring method for bending deformation of simply supported beam Active CN111189427B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010090290.9A CN111189427B (en) 2020-02-13 2020-02-13 Real-time monitoring method for bending deformation of simply supported beam

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010090290.9A CN111189427B (en) 2020-02-13 2020-02-13 Real-time monitoring method for bending deformation of simply supported beam

Publications (2)

Publication Number Publication Date
CN111189427A CN111189427A (en) 2020-05-22
CN111189427B true CN111189427B (en) 2021-09-24

Family

ID=70706634

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010090290.9A Active CN111189427B (en) 2020-02-13 2020-02-13 Real-time monitoring method for bending deformation of simply supported beam

Country Status (1)

Country Link
CN (1) CN111189427B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111561903B (en) * 2020-07-14 2020-10-30 北京讯腾智慧科技股份有限公司 Bridge deformation monitoring system and method
CN113654747B (en) * 2021-09-26 2024-04-16 珠海格力智能装备有限公司 Method and device for detecting joint stiffness of robot and robot
CN116026414B (en) * 2023-02-14 2023-12-19 中交第三航务工程局有限公司 Monitoring system and monitoring method for integrated bridge girder erection machine

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201731984U (en) * 2010-07-27 2011-02-02 大连理工大学 Wireless inclination based bridge structure deflection measurement system
KR101193076B1 (en) * 2011-03-08 2012-10-22 주식회사 포스코건설 System for real-time measuring penetration depth of pile using 3D photogrammetry
CN103148829B (en) * 2013-02-25 2015-11-04 桂林理工大学 Based on the structural deformation detection method of Internet of Things
CN205808651U (en) * 2016-05-09 2016-12-14 中国地震局工程力学研究所 A kind of bridge beam body deformation sensor-based system and sensor
CN107784138B (en) * 2016-08-26 2021-03-26 上海宝冶集团有限公司 Point cloud gravity deformation correction method based on structural mechanics analysis
CN106771452A (en) * 2016-11-15 2017-05-31 中国电子科技集团公司第四十研究所 A kind of waveform drawing and acquisition method based on MFC class libraries and screen resolution
CN107036581B (en) * 2017-05-16 2019-12-13 福建三鑫隆信息技术开发股份有限公司 bridge deformation online monitoring system and method based on MEMS gyroscope
CN207317725U (en) * 2017-06-20 2018-05-04 燕山大学 A kind of measuring device of beam deflection
CN110346103B (en) * 2019-06-18 2021-02-02 东南大学 Deformation analysis method of MEMS cantilever beam
CN110346104B (en) * 2019-06-18 2021-02-02 东南大学 Deformation analysis method for MEMS double-end clamped beam

Also Published As

Publication number Publication date
CN111189427A (en) 2020-05-22

Similar Documents

Publication Publication Date Title
CN111189427B (en) Real-time monitoring method for bending deformation of simply supported beam
CN100595514C (en) Cable system health monitoring method based on strain monitoring cable structure
CN101387600B (en) Cable system health monitoring method in cable structure based on mixed monitor
CN101504333B (en) Progressive health monitoring method for cable structure median cord system based on strain monitoring
US20100206040A1 (en) Method for calibrating parameter of articulated coordinate measuring apparatus
CN102221478A (en) Progressive health monitoring method for cable system based on hybrid monitoring during generalized displacement of support
CN101793621B (en) Progressive health monitoring method for cable system based on cable force monitoring during support settlement
CN111783326B (en) Wind speed and dynamic response test-based wind-induced safety assessment method for power transmission tower structure
CN111189428B (en) Real-time monitoring method for bending deformation of cantilever beam
CN101788400B (en) Progressive health monitoring method of cable system based on strain monitoring during support settlement
CN101482434B (en) Progressive health monitoring method for cable system in cable structure based on cable force monitoring
CN110160492A (en) A kind of the inclination monitoring method and device of power transmission tower
CN102323080A (en) The health monitor method of the cable system of bearing generalized displacement time space coordinate monitoring
CN102323082A (en) Health monitoring method of cable system based on strain monitoring during support generalized displacement
CN102288428A (en) Method for distinguishing slack cable based on hybrid monitoring in angular displacement of supporting seat
CN102323083A (en) Strain monitoring based method for recognizing damaged cable, loose cable and supporting seat generalized displacement
CN111896200B (en) Method for testing bending rigidity of short cable
CN102323079B (en) Angle monitoring based cable system health monitoring method applied in supporting seat generalized displacement
CN210710213U (en) Device for detecting inclination angle of elevator car based on six-axis acceleration sensor
CN102288436A (en) Method for identifying slack supporting cables based on space coordinate monitoring during angular displacement of support
CN102323093A (en) Angle monitoring based cable system progressive health monitoring method applied in supporting seat generalized displacement
CN102297771A (en) Progressive health monitoring method for cable system based on angle monitoring during supporting-base angular displacement
CN102288427A (en) Progressive health monitoring method of cable system based on strain monitoring in angular displacement of supporting seat
CN102288426A (en) Progressive health monitoring method of cable system based on spatial coordinate monitoring in angular displacement of supporting seat
CN102288434B (en) Method for identifying slack cables based on angle monitoring during angular displacement of support

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20200522

Assignee: Guilin Shenghi Tech Co.,Ltd.

Assignor: GUILIN University OF TECHNOLOGY

Contract record no.: X2022450000152

Denomination of invention: Real time monitoring method for bending deformation of simply supported beams

Granted publication date: 20210924

License type: Common License

Record date: 20221124

EE01 Entry into force of recordation of patent licensing contract