CN106156410B - Thin-wall circular tube structure axial direction Calculation of Thermal Deformation method - Google Patents
Thin-wall circular tube structure axial direction Calculation of Thermal Deformation method Download PDFInfo
- Publication number
- CN106156410B CN106156410B CN201610476246.5A CN201610476246A CN106156410B CN 106156410 B CN106156410 B CN 106156410B CN 201610476246 A CN201610476246 A CN 201610476246A CN 106156410 B CN106156410 B CN 106156410B
- Authority
- CN
- China
- Prior art keywords
- thin
- circular tube
- wall circular
- tube structure
- thermal deformation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000004044 response Effects 0.000 claims abstract description 19
- 238000004088 simulation Methods 0.000 claims abstract description 18
- 239000000835 fiber Substances 0.000 claims abstract description 17
- 239000013307 optical fiber Substances 0.000 claims description 20
- 230000005284 excitation Effects 0.000 claims description 16
- 230000000875 corresponding effect Effects 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 11
- 238000012544 monitoring process Methods 0.000 claims description 11
- 230000002596 correlated effect Effects 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 3
- 230000005669 field effect Effects 0.000 claims description 2
- 239000004744 fabric Substances 0.000 claims 3
- 230000015572 biosynthetic process Effects 0.000 abstract 1
- 238000011068 loading method Methods 0.000 abstract 1
- 230000009471 action Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000009795 derivation Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 102100038032 F-box only protein 17 Human genes 0.000 description 1
- 102100022116 F-box only protein 2 Human genes 0.000 description 1
- 102100022115 F-box only protein 27 Human genes 0.000 description 1
- 102100026080 F-box only protein 44 Human genes 0.000 description 1
- 102100024513 F-box only protein 6 Human genes 0.000 description 1
- 101000878584 Homo sapiens F-box only protein 17 Proteins 0.000 description 1
- 101000824158 Homo sapiens F-box only protein 2 Proteins 0.000 description 1
- 101000824171 Homo sapiens F-box only protein 27 Proteins 0.000 description 1
- 101000913298 Homo sapiens F-box only protein 44 Proteins 0.000 description 1
- 101001052796 Homo sapiens F-box only protein 6 Proteins 0.000 description 1
- 101100254587 Lentinula edodes RPS1 gene Proteins 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000007306 functionalization reaction Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
- Optical Transform (AREA)
Abstract
The present invention relates to thin-wall circular tube structure axial direction Calculation of Thermal Deformation methods, for two kinds of different temperature fields loading methods, the axial direction thermal deformation of thin-wall circular tube class formation is measured using fiber-optic grating sensor, this method, which is mainly comprised the steps that, carries out numerical simulation using axial thermal deformation associated arguments feature of the emulation mode to thin-wall circular tube structure, provides thin-wall circular tube structure axial direction Calculation of Thermal Deformation method according to analog result；Thermal force response signal suffered by the thin-wall circular tube structure is monitored using fibre-optical grating sensor, the response signal is substituted into Axial Thermal settlement calculation model, thin-wall circular tube structure axial direction thermal deformation is calculated.This method is by acquiring the fiber bragg grating center wavelength response signal of a small amount of discrete point, by the axial thermal deformation of relevant calculation model inference different temperatures thin-wall circular tube structure off field.
Description
Technical field
The invention belongs to the fields of monitoring structural health conditions, specifically propose a kind of based on distributed fiber grating sensing network
Different temperatures field action under thin-wall circular tube structure axial direction thermal deformation quick calculation method.
Background technique
In outer space, thermal force is to act on spaceborne Main Load.The large space structure in structure type
Often using the structure by thin-wall circular tube structure composition, temperature change occurs in Thermal Load down space structure, generates thermal change
Shape, thermal stress, it is also possible to thermal buckling and thermal vibration occur, so that spacecraft be made to fail.In order to disclose these due to thermal force and
Caused space structure thermal deformation situation it may first have to which thermal deformation of the structure under different temperatures field action is supervised in real time
It surveys, and then more accurately different Thermal Load flowering structure deformation states suffered by analysis space structure.Meanwhile analyzing structure
Deformation under different Thermal Loads will be helpful to provide corresponding size remaining when practical structures design, to guarantee
Still there is good environmental suitability after actual ambient temperature variation.Herein for thin-wall circular tube knot common in space structure
Structure has studied a kind of quick calculation method of structure axial direction thermal deformation.
Currently, structural deformation monitoring is received significant attention in aerospace field.Future space structure will be more towards more
Functionalization, multi-tasking direction are developed, this will certainly perceive planform and propose requirements at the higher level.And how quickly and easily to obtain
Deformation state information of the space structure under Thermal Load is taken, for assessment space structure military service safe condition and later maintenance
It is of great significance.Therefore, thin-walled under the present invention relates to a kind of different temperatures field action based on fiber-optic grating sensor
The quick calculation method of circular tube structure axial direction thermal deformation.
Summary of the invention
To solve above-mentioned technology, the present invention provides a kind of thin-wall circular tube structure axial direction Calculation of Thermal Deformation method, this method
It is relatively simple not need mass data sample point, process, can quickly calculate thin-wall circular tube structure heat deformable state.This method
Are as follows:
It is counted using axial thermal deformation correlated characteristic parameter of the emulation mode to thin-wall circular tube structure under Thermal Load
Value simulation, provides thin-wall circular tube structure axial direction Calculation of Thermal Deformation method according to analog result；It is sensed using distributed fiber grating
Device monitors thermal force response signal suffered by the thin-wall circular tube structure, and the response signal is substituted into axial Calculation of Thermal Deformation mould
Type calculates the axial direction thermal deformation of thin-wall circular tube structure.
This method is used for thin-wall circular tube structure axial direction thermal deformation fiber-optic monitoring under steady state excitation environment, described
Numerical simulation is carried out using axial thermal deformation correlated characteristic parameter of the emulation mode to thin-wall circular tube structure, is given according to analog result
Thin-wall circular tube structure axial direction Calculation of Thermal Deformation method out specifically:
Step 1: the thermal deformation of the thin-wall circular tube structure of steady state excitation and strain are carried out using emulation mode
Sunykatuib analysis, respectively simulation obtain thin-wall circular tube structure be under different homogeneous temperature field environment along the thermal deformation of axial path with
And strain；
Step 2: building thin-wall circular tube structure any point strain and the point in steady state excitation environment lower axis
The relational model of place thin-wall circular tube structure axial coordinate.
Thin-wall circular tube axial coordinates x is established, the thin-wall circular tube structure is under steady state excitation environment along axis
The relational model of the thermal strain of any position and thin-wall circular tube structure axial direction total deformation on the outside of the tube wall of direction are as follows:
Wherein, εr1、εr2、εr3Respectively thin-wall circular tube outer surface is along axial x=0, x=L1/ 12, x=L1It is corresponding at/6
Thermal strain, εr4For the corresponding thermal strain of thin-wall circular tube outer surface of tube wall axis midpoint, L1For thin-wall circular tube length, Δ L1
For thin-wall circular tube structure steady state excitation effect under corresponding axial deflection.
It is described to monitor thermal force response signal suffered by the thin-wall circular tube structure using fibre-optical grating sensor, it will
The response signal substitutes into Axial Thermal settlement calculation model, calculates the axial direction thermal deformation of thin-wall circular tube structure and refers specifically to:
Based on the model, thin-wall circular tube structure axial coordinates and arrangement thin-walled when without Thermal Load are established initially
Circular tube structure distributing optical fiber sensing network, by fiber arrangement in thin-wall circular tube outer surface along axial x=0, x=L1/ 12, x=
L1/ 6 and axis midpoint；Assuming that thin-wall circular tube initial temperature is T0, and record at this temperature each fiber-optic grating sensor just
Beginning central wavelength, thin-wall circular tube are placed in homogeneous temperature field, record every fiber-optic grating sensor central wavelength respectively；It will be each
Fiber-optic grating sensor center wavelength shift amount is converted into thermal strain, and substitutes into the Axial Thermal settlement calculation model, obtains steady
Thin-wall circular tube structure axial direction thermal deformation under state uniform temperature field action:
Wherein, Δ λ1、Δλ2、Δλ3、Δλ4For each fiber-optic grating sensor center wavelength shift amount.
This method is used for uneven steady-state temperature field thin-wall circular tube structure axial direction thermal deformation fiber-optic monitoring, the use
Emulation mode carries out numerical simulation to the axial thermal deformation correlated characteristic parameter of thin-wall circular tube structure, provides according to analog result thin
Wall circular tube structure Axial Thermal Method for Calculating Deformation specifically:
Step 1: using emulation mode to the thermal deformation of the thin-wall circular tube structure of uneven steady-state temperature field and strain into
Row sunykatuib analysis obtains thermal deformation and strain along thin-wall circular tube structure axial path；
Step 2: thin-wall circular tube axial coordinates x, the thin-wall circular tube knot under uneven steady-state temperature field effect are established
The structure relationship on outer surface of tube wall between any point thermal strain and the axial position coordinate of the point in the axial direction:
Wherein, coefficient a2, coefficient b2With coefficient c2, expression is as follows:
Wherein, εr1′、εr2′、εr2′Respectively thin-wall circular tube outer surface is along axial x=0, x=L2/ 2, x=L2Locate corresponding
Thermal strain, L2For the length of thin-wall circular tube.
It is described to monitor thermal force response signal suffered by the thin-wall circular tube structure using fibre-optical grating sensor, it will
The response signal substitutes into Axial Thermal settlement calculation model, calculates the axial direction thermal deformation of thin-wall circular tube structure specifically:
Based on the model, thin-wall circular tube structure axial coordinates and arrangement thin-walled when without Thermal Load are established initially
Circular tube structure distributing optical fiber sensing network, by fiber arrangement in thin-wall circular tube outer surface along axial x=0, x=L2/ 2, x=L2
Place, and optical fiber is arranged at the compensation position of each position；Assuming that thin-wall circular tube initial temperature is T0, and record each light at this temperature
Fiber grating sensor initial center wavelength, thin-wall circular tube are placed in non-uniform temperature field.Under the effect of non-uniform temperature load,
Every fiber-optic grating sensor central wavelength is recorded respectively；Each fiber-optic grating sensor center wavelength shift amount is substituted into Axial Thermal
Settlement calculation model obtains thin-wall circular tube structure axial direction thermal deformation are as follows:
Wherein,
b2=(Δ λ6-Δλ9)/1.2-a2；
Wherein, Δ λ6、Δλ7、Δλ8、Δλ9、Δλ10、Δλ11For each fiber-optic grating sensor center wavelength shift amount, Δ
L2For the length of thin-wall circular tube axial direction thermal deformation.
This method is used for thin-wall circular tube structure axial direction thermal deformation fiber-optic monitoring under uneven steady-state temperature field environment, described
Numerical simulation is carried out to the axial thermal deformation correlated characteristic parameter of thin-wall circular tube structure using emulation mode, according to analog result
Provide thin-wall circular tube structure axial direction Calculation of Thermal Deformation method specifically:
Step 1: using emulation mode to the thermal deformation of the thin-wall circular tube structure of uneven steady-state temperature field and strain into
Row sunykatuib analysis obtains thermal deformation and Temperature Distribution along axial path thin-wall circular tube structural outer surface；
Step 2: establishing thin-wall circular tube axial coordinates x, and uneven steady-state temperature field acts on lower thin-wall circular tube structure appearance
The relational expression between the temperature value of any position and the axial coordinate of the point is gone up in the axial direction in face are as follows:
Wherein, coefficient a3, coefficient b3With coefficient c3, expression is as follows:
Wherein, T1、T2、T3Respectively thin-wall circular tube outer wall is axially x=0, x=L3/ 2, x=L3The temperature value at place, L3For
The length of thin-wall circular tube.
It is described to monitor thermal force response signal suffered by thin-wall circular tube structure using fibre-optical grating sensor, it will be described
Response signal substitutes into Axial Thermal settlement calculation model, calculates the axial direction thermal deformation of thin-wall circular tube structure and refers specifically to:
Based on the model, thin-wall circular tube structure axial coordinates and arrangement thin-walled when without Thermal Load are established initially
Circular tube structure distributing optical fiber sensing network, by fiber arrangement in thin-wall circular tube outer surface along axial x=0, x=L3/ 2, x=L3
Place；Assuming that thin-wall circular tube initial temperature is T0, and record each fiber-optic grating sensor initial center wavelength at this temperature, thin-wall circular
Pipe is placed in non-uniform temperature field.Under the effect of non-uniform temperature load, every fiber-optic grating sensor center is recorded respectively
The temperature that each fiber-optic grating sensor center wavelength shift amount is converted to thin-wall circular tube is substituted into axial Calculation of Thermal Deformation mould by wavelength
Type obtains thin-wall circular tube structure axial direction thermal deformation are as follows:
Wherein,
Wherein, T1=Δ λ12/10+T0、T2=Δ λ13/10+T0、T3=Δ λ14/10+T0, Δ λ12、Δλ13、Δλ14It is each
Grating sensor center wavelength shift amount, Δ L3For the thermal deformation length of thin-wall circular tube.
Thin-wall circular tube structure axial direction thermal deformation optical fiber calculation method provided by the invention, fiber grating sensing technology have body
Product small, light weight, at low cost, being easily integrated, convenient for distributed measurement, resistance to corrosion is strong the advantages that, by more and more extensively
General concern.Based on above-mentioned analysis, the present invention proposes to use fibre-optical grating sensor institute measured data, calculates different temperature fields
Act on the thermal deformation of lower thin-wall circular tube axial direction.
This method can quickly calculate the thin-wall circular tube structure thermal deformation under two kinds of different temperatures field actions, be one
The feasible Calculation of Thermal Deformation method of kind.
Relative to other methods calculated of temperature field deformation, this method does not need mass data sample point, and operation is more
Simply.
Detailed description of the invention
Fig. 1 is the schematic diagram of thin-wall circular tube structure under steady state excitation；
Fig. 2 is thin-wall circular tube structure L under different steady state excitations1Axial elasticity strain figure when=300mm；
Fig. 3 is thin-wall circular tube structure L under different steady state excitations1Thermal deformation figure when=300mm；
Fig. 4 is the schematic diagram of thin-wall circular tube structure under uneven steady-state temperature field；
Fig. 5 is thin-wall circular tube structure L under different uneven steady-state temperature fields2Axial Thermal strain figure when=300mm；
Fig. 6 is thin-wall circular tube structure L under different uneven steady-state temperature fields2Axial Thermal deformation pattern when=300mm；
Fig. 7 is the schematic diagram of thin-wall circular tube structure under uneven steady-state temperature field；
Fig. 8 is thin-wall circular tube structure L under different uneven steady-state temperature fields3Axial temperature change curve when=300mm
Figure.
Specific embodiment
The technical solution of invention is described in detail with reference to the accompanying drawing: for the axial thermal change of thin-wall circular tube structure
Shape proposes a kind of thin-wall circular tube structure under Thermal Load, is realized using distributed fiber grating sensing network to thin-walled
The method for solving of circular tube structure axial direction thermal deformation.Specific embodiment is as follows:
1) calculation method of steady state excitation thin-wall circular tube axial direction thermal deformation, comprising the following steps:
Step 1: it to thin-wall circular tube structure axial direction thermal deformation under different steady state excitations and is answered using emulation mode
Become and carry out sunykatuib analysis, obtains the thermal deformation and strain of thin-wall circular tube structure outer wall surface in the axial direction；
Step 2: building thin-wall circular tube structure outer wall surface under steady state excitation is answered at any point in the axial direction
Become the relational model between the thin-wall circular tube structure axial coordinate of point place；
(1) as the above analysis, coordinate system is initially set up, establishes x along any axis direction of thin-wall circular tube outer wall surface
Axis.According to numerical simulation result, the thermoelastic strain at thin-wall circular tube structure both ends is distributed in conic section, and interlude tends to be steady
It is fixed, therefore strain monitoring point is set at tetra- positions point A1, A2, A3, B respectively, 4 points of A1, A2, A3, B coordinates in x-axis and
Its corresponding thermal strain is respectively A1 (0, εr1), A2 (L1/ 12, εr2), A3 (L1/ 6, εr3)、B(L1/ 2, εr4)；
(2) according to finite element stimulation as a result, thin-wall circular tube structure outer wall surface in the axial direction each point strain and its
Corresponding axial coordinate is in quadratic function relation, and interlude is constant.Therefore it may be assumed that the strain of each point on the same axis and the point
Axial coordinate between relationship are as follows:
Wherein x is the coordinate at any point in the axial direction, εrFor thermal strain corresponding with x coordinate；
In view of the symmetry of thin-wall circular tube structure itself, the deflection of thin-wall circular tube structure half length can be studied, i.e.,
Study 0≤x≤L1The deflection Δ L of/2 parts1′.The then structural strain's amount Δ L of thin-wall circular tube1=2 Δ L1′。
According to 3 thermal strains of A1, A2, A3, coefficient a can be calculated1, coefficient b1With coefficient c1It is respectively as follows:
c1=εr1
According to the strain of B point, d can be calculated1Value, expression is as follows:
d1=εr4
The then axial deflection of the front half section (OB sections) of thin-wall circular tube structure are as follows:
According to thin-wall circular tube symmetrical configuration feature, thin-wall circular tube axial direction total deformation can be deduced are as follows:
Further abbreviation is available:
Step 3: being based on above-mentioned thin-wall circular tube model, and thin-wall circular tube structure is axially sat when establishing initially without Thermal Load
Mark system and arrangement thin-wall circular tube structure distribution formula optical fiber sensing network；
(1) coordinate system is initially set up, selects any one axis of thin-wall circular tube structure outer wall surface as x-axis, according to numerical value
The thermoelastic strain of analog result, thin-wall circular tube structure both ends is distributed in conic section, and interlude tends towards stability, therefore respectively in point
Fiber-optic grating sensor FBG1, FBG2, FBG3, FBG4 are pasted at A1, A2, A3, B, coordinate defines along the x-axis direction by point A1, A2, A3
For xA1=0, xA2=L1/12、xA3=L1Abscissa is x at/6, point BB=L1/ 2, B point are thin-wall circular tube structure axis midpoint.This
The stickup direction of four fiber-optic grating sensors is located on the same axis of thin-wall circular tube structure outer wall, by optical fiber grating sensing
Device is pasted on structure outer wall, uses optical patchcord to be connected in series four fiber-optic grating sensors and is constituted distributed biography with this
Sensor network, and temperature compensation optical fiber grating FBG5 is disposed about in thin-wall circular tube.
(2) assume that thin-wall circular tube initial temperature is T0, and record each fiber-optic grating sensor initial center wave at this temperature
It is long.Thin-wall circular tube is placed in homogeneous temperature field, records the collected central wavelength lambda of every fiber-optic grating sensor respectively1、λ2、
λ3、λ4、λ5, wherein λ5For compensated optical fiber grating sensor central wavelength.Each fiber-optic grating sensor cuts its initial center wave
It is long, it obtains each fiber-optic grating sensor center wavelength shift amount and is followed successively by Δ λ1′、Δλ2′、Δλ3′、Δλ4′、Δλ5′.Here each
Fiber grating center wavelength shift amount due to caused by thermal strain is respectively as follows:
Δλ1=Δ λ1′-Δλ5′
Δλ2=Δ λ2'=Δ λ5′
Δλ3=Δ λ3′-Δλ5′
Δλ4=Δ λ4′-Δλ5′
Then strain corresponding to every fiber-optic grating sensor center wavelength shift amount is successively are as follows: εr1、εr2、εr3、εr4, i.e.,
It is respectively as follows:
εr1=Δ λ1/1.2
εr2=Δ λ2/1.2
εr3=Δ λ3/1.2
εr4=Δ λ4/1.2
According to above-mentioned derivation formula, axial deflection of the thin-wall circular tube under steady temperature environment can be obtained are as follows:
2) thin-wall circular tube Axial Thermal strain gauge algorithm packet under non-homogeneous Steady-State Thermal Field (one end applies steady temperature load)
Include following steps:
Step 1: simulation point is carried out to the thermal strain of the thin-wall circular tube structure of uneven steady-state temperature field using simulation software
Analysis, obtains the thermal strain feature along thin-wall circular tube outer wall surface axis direction；
Step 2: according to numerical simulation result, thin-wall circular tube outer wall is along axial thermal strain under non-uniform temperature field action
In specific function relationship, therefore in actual measurement, thin-wall circular tube can be estimated according to the thermal strain of several key points of structure
Axial deflection.
(1) building thin-wall circular tube structure stable state non-uniform temperature load effect lower axis on any point thermal strain with it is thin
The relational model of the axial coordinate of wall circular tube structure.X-axis is established along any axis of thin-wall circular tube outer wall surface, according to numerical simulation
As a result, exponentially curve distribution feature, fire end maximum temperature are heated with separate for the thermal strain at thin-wall circular tube structure both ends
End, body structure surface temperature gradually decay and tend towards stability, therefore strain monitoring point, point is arranged at tri- positions point A, B, C respectively
A, coordinate is defined as x along the x-axis direction by B, CA=0, xB=L2/2,xC=L2.A, 3 points of B, C coordinates in x-axis and its corresponding
Thermal strain is respectively A (0, εr1'), B (L2/ 2, εr2'), C (L2, εr3′)；
(2) by above-mentioned derivation it is found that thin-wall circular tube structural strain gradient and axis coordinate exponentially functional relation, therefore can be false
If the thin-wall circular tube outer wall in axial direction relationship on any axis between the strain and the axial coordinate of the point of each point are as follows:
According to 3 points of strain of A, B, C, coefficient a can be calculated separately out2, coefficient b2With coefficient c2, it is specific as follows:
The then thermal deformation of thin-wall circular tube structure axial direction:
It is further simplified:
Step 3: being based on above-mentioned thin-wall circular tube model, and thin-wall circular tube structure is axially sat when establishing initially without Thermal Load
Mark system and arrangement thin-wall circular tube structure distribution formula optical fiber sensing network；
(1) initially set up coordinate system, using any one axis of thin-wall circular tube structure outer wall surface as x-axis, respectively point A,
B, fiber-optic grating sensor FBG6, FBG7, FBG8 are pasted at C, and temperature is successively set near each fiber-optic grating sensor and is mended
Fiber grating FBG9, FBG10, FBG11 are repaid, this six fiber-optic grating sensors are pasted on test specimen structure table each along axial direction
Face.Three fiber-optic grating sensors and three fiber grating temperature sensors are connected in series respectively using optical patchcord, with this
Constitute distributed sensor networks.
(2) assume that thin-wall circular tube initial temperature is T0, and record each fiber-optic grating sensor initial center wave at this temperature
It is long.Thin-wall circular tube is placed in non-uniform temperature field, records the collected central wavelength lambda of every fiber-optic grating sensor respectively6、
λ7、λ8、λ9、λ10、λ11, wherein λ9、λ10、λ11For each temperature compensation optical fiber grating center sensor wavelength.Each optical fiber grating sensing
Device
Subtract its initial center wavelength, available each fiber-optic grating sensor center wavelength shift amount Δ λ6、Δλ7、Δ
λ8、Δλ9、Δλ10、Δλ11, then each fiber-optic grating sensor center wavelength shift amount as caused by thermal strain is respectively as follows:
Δλ6'=Δ λ6-Δλ9
Δλ7'=Δ λ7-Δλ10
Δλ8'=Δ λ8-Δλ11
According to fiber Bragg grating strain sensor center wavelength shift amount, the strain of available A, B, C point is respectively εr1′、
εr2′、εr3':
εr1'=(Δ λ6-Δλ9)/1.2
εr2'=(Δ λ7-Δλ10)/1.2
εr3'=(Δ λ8-Δλ11)/1.2
Axial deflection of the thin-wall circular tube under stable state non-uniform temperature environment are as follows:
Wherein,
b2=(Δ λ6-Δλ9)/1.2-a2；
3) thin-wall circular tube Axial Thermal deformation method includes under non-homogeneous Steady-State Thermal Field (one end application stationary temperature)
Following steps:
Step 1: the thin-wall circular tube structure outside wall temperature distribution of uneven steady-state temperature field is counted using emulation mode
Value simulation, obtains the variation characteristic along thin-wall circular tube outer wall surface axis direction temperature；
Step 2: according to numerical simulation result, thin-wall circular tube outer wall becomes along axial temperature under non-uniform temperature field action
Change exponentially functional relation, therefore in actual measurement, thin-wall circular can also be estimated according to the temperature value of several key nodes
The axial deflection of pipe.
(1) building thin-wall circular tube structure take up an official post in the axial direction under stable state non-uniform temperature load anticipate some temperature with it is thin
Relational model between the axial coordinate of wall circular tube structure.X-axis is established along thin-wall circular tube outer wall surface axis direction, according to numerical value
Sunykatuib analysis, thin-wall circular tube outside wall temperature add along the axial exponentially curve distribution of tube wall, fire end maximum temperature with separate
Hot end, thin-wall circular tube outside wall temperature gradually decay and tend towards stability, therefore temperature monitoring is arranged at tri- positions point A, B, C respectively
Point, coordinate is defined as x along the x-axis direction by point A, B, CA=0, xB=L3/2、xC=L3.A, 3 points of B, C coordinates in x-axis and its right
The temperature answered is respectively A (0, T1), B (L3/ 2, T2), C (L3, T3)；
(2) according to simulation result, thin-wall circular tube structure temperature gradient and axis coordinate exponentially functional relation, therefore can
Assuming that the relationship on any axis between each point temperature and the axial coordinate of the point are as follows:
According to 3 points of temperature of A, B, C, coefficient a can be calculated separately3, coefficient b3With coefficient c3, it is specific as follows:
The then thermal deformation of thin-wall circular tube structure can approximate estimation are as follows:
Wherein, α is thin-wall circular tube structural material thermal expansion coefficient, thin-wall circular tube structural material isotropic.
It is further simplified:
Step 3: it is based on above-mentioned thin-wall circular tube model, establishes thin-wall circular tube structure initially without the axial direction under Thermal Load
Coordinate system and arrangement thin-wall circular tube structure distribution formula optical fiber sensing network；
(1) coordinate system is initially set up, selects thin-wall circular tube structure axial direction as x-axis, is pasted at point A, B, C respectively
Fiber-optic grating sensor FBG12, FBG13, FBG14, three fiber-optic grating sensors are pasted on structure outer wall along axis direction.
Three fiber grating temperature sensors will be connected in series respectively using optical patchcord, distributed sensor networks are constituted with this.
(2) assume that thin-wall circular tube initial temperature is T0, and record each fiber-optic grating sensor initial center wave at this temperature
It is long.Thin-wall circular tube is placed in non-uniform temperature field, records the every collected central wavelength of fiber-optic grating sensor institute respectively
λ12、λ13、λ14, each fiber-optic grating sensor subtracts its initial center wavelength, obtains each fiber-optic grating sensor center wavelength shift
Measure Δ λ12、Δλ13、Δλ14, the temperature value T of available A, B, C point1、T2、T3It is respectively as follows:
T1=Δ λ12/10+T0
T2=Δ λ13/10+T0
T3=Δ λ14/10+T0
The then thermal deformation of thin-wall circular tube structure axial direction
Wherein,
Claims (7)
1. thin-wall circular tube structure axial direction Calculation of Thermal Deformation method, which is characterized in that using emulation mode to thin under Thermal Load
The axial thermal deformation correlated characteristic parameter of wall circular tube structure carries out numerical simulation, provides thin-wall circular tube axis of no-feathering according to analog result
To Calculation of Thermal Deformation method；The response letter of thermal force suffered by the thin-wall circular tube structure is monitored using fibre-optical grating sensor
Number, the response signal is substituted into Axial Thermal settlement calculation model, calculates the axial direction thermal deformation of thin-wall circular tube structure；
This method is used for thin-wall circular tube structure axial direction thermal deformation fiber-optic monitoring, the use under steady state excitation environment
Emulation mode carries out numerical simulation to the axial thermal deformation correlated characteristic parameter of thin-wall circular tube structure, provides according to analog result thin
Wall circular tube structure Axial Thermal Method for Calculating Deformation specifically:
Step 1: the thermal deformation and strain of the thin-wall circular tube structure of steady state excitation are simulated using emulation mode
Analysis, simulation obtains thin-wall circular tube structure and is under different homogeneous temperature field environment along the thermal deformation of axial path and answers respectively
Become；
Step 2: building thin-wall circular tube structure any point strain where with the point in steady state excitation environment lower axis
The relational model of thin-wall circular tube structure axial coordinate.
2. thin-wall circular tube structure axial direction Calculation of Thermal Deformation method according to claim 1, which is characterized in that establish thin-wall circular
Pipe axial coordinates x, the thin-wall circular tube structure any position on the outside of tube wall in the axial direction under steady state excitation environment
The relational model of the thermal strain and thin-wall circular tube structure axial direction total deformation set are as follows:
Wherein, εr1、εr2、εr3Respectively thin-wall circular tube outer surface is along axial x=0, x=L1/ 12, x=L1Corresponding heat at/6
Strain, εr4For the corresponding thermal strain of thin-wall circular tube outer surface of tube wall axis midpoint, L1For thin-wall circular tube length, Δ L1It is thin
Wall circular tube structure corresponding axial deflection under steady state excitation effect.
3. thin-wall circular tube structure axial direction Calculation of Thermal Deformation method according to claim 1, which is characterized in that described utilize is divided
Cloth fiber-optic grating sensor monitors thermal force response signal suffered by the thin-wall circular tube structure, and the response signal is substituted into axis
To Calculation of Thermal Deformation model, calculates the axial direction thermal deformation of thin-wall circular tube structure and refers specifically to:
Based on the relational model of the thin-wall circular tube structure axial coordinate, thin-wall circular tube structure when without Thermal Load is established initially
Axial coordinates and arrangement thin-wall circular tube structure distribution formula optical fiber sensing network, by fiber arrangement on thin-wall circular tube outer surface edge
Axial x=0, x=L1/ 12, x=L1/ 6 and axis midpoint；Assuming that thin-wall circular tube initial temperature is T0, and record the temperature
Under each fiber-optic grating sensor initial center wavelength, thin-wall circular tube is placed in homogeneous temperature field, records every optical fiber respectively
Grating sensor central wavelength；Each fiber-optic grating sensor center wavelength shift amount is converted into thermal strain, and substitutes into the axis
To Calculation of Thermal Deformation model, obtains steady state excitation and acts on lower thin-wall circular tube structure axial direction thermal deformation:
Wherein, Δ λ1、Δλ2、Δλ3、Δλ4For each fiber-optic grating sensor center wavelength shift amount；
L1For thin-wall circular tube length, Δ L1For the corresponding axial deformation under steady state excitation effect of thin-wall circular tube structure
Amount.
4. thin-wall circular tube structure axial direction Calculation of Thermal Deformation method according to claim 1, which is characterized in that this method is used for
To uneven steady-state temperature field thin-wall circular tube structure axial direction thermal deformation fiber-optic monitoring, the use emulation mode is to thin-wall circular tube
The axial thermal deformation correlated characteristic parameter of structure carries out numerical simulation, provides thin-wall circular tube structure axial direction thermal change according to analog result
Shape calculation method specifically:
Step 1: mould is carried out to the thermal deformation of the thin-wall circular tube structure of uneven steady-state temperature field and strain using emulation mode
Quasi- analysis, obtains thermal deformation and strain along thin-wall circular tube structure axial path；
Step 2: thin-wall circular tube axial coordinates x, the thin-wall circular tube structure edge under uneven steady-state temperature field effect are established
Relationship on axis direction outer surface of tube wall between any point thermal strain and the axial position coordinate of the point:
Wherein, coefficient a2, coefficient b2With coefficient c2, expression is as follows:
Wherein, εr1′、εr2′、εR3,Respectively thin-wall circular tube outer surface is along axial x=0, x=L2/ 2, x=L2Locate corresponding heat to answer
Become, L2For the length of thin-wall circular tube.
5. thin-wall circular tube structure axial direction Calculation of Thermal Deformation method according to claim 4, which is characterized in that described utilize is divided
Cloth fiber-optic grating sensor monitors thermal force response signal suffered by the thin-wall circular tube structure, and the response signal is substituted into axis
To Calculation of Thermal Deformation model, the axial direction thermal deformation of thin-wall circular tube structure is calculated specifically:
Based on the relational model of the thin-wall circular tube structure axial coordinate, thin-wall circular tube structure when without Thermal Load is established initially
Axial coordinates and arrangement thin-wall circular tube structure distribution formula optical fiber sensing network, by fiber arrangement on thin-wall circular tube outer surface edge
Axial x=0, x=L2/ 2, x=L2Place, and in the setting compensated optical fiber grating sensor of each position；Assuming that thin-wall circular tube is initially warm
Degree is T0, and each fiber-optic grating sensor initial center wavelength, thin-wall circular tube are placed on non-uniform temperature field to record at this temperature
In；Under the effect of non-uniform temperature load, every fiber-optic grating sensor central wavelength is recorded respectively；By each optical fiber grating sensing
Device center wavelength shift amount substitutes into Axial Thermal settlement calculation model, obtains thin-wall circular tube structure axial direction thermal deformation are as follows:
Wherein,
b2=(Δ λ6-Δλ9)/1.2-a2；
Wherein, Δ λ6、Δλ7、Δλ8、Δλ9、Δλ10、Δλ11For each fiber-optic grating sensor center wavelength shift amount, Δ L2For
The length of thin-wall circular tube axial direction thermal deformation.
6. thin-wall circular tube structure axial direction Calculation of Thermal Deformation method according to claim 1, which is characterized in that this method is used for
To thin-wall circular tube structure axial direction thermal deformation fiber-optic monitoring under uneven steady-state temperature field environment, the use emulation mode is to thin
The axial thermal deformation correlated characteristic parameter of wall circular tube structure carries out numerical simulation, provides thin-wall circular tube axis of no-feathering according to analog result
To Calculation of Thermal Deformation method specifically:
Step 1: mould is carried out to the thermal deformation of the thin-wall circular tube structure of uneven steady-state temperature field and strain using emulation mode
Quasi- analysis, obtains thermal deformation and Temperature Distribution along axial path thin-wall circular tube structural outer surface；
Step 2: establishing thin-wall circular tube axial coordinates x, and uneven steady-state temperature field acts on lower thin-wall circular tube structural outer surface edge
Relational expression on axis direction between the temperature value of any position and the axial coordinate of the point are as follows:
Wherein, coefficient a3, coefficient b3With coefficient c3, expression is as follows:
Wherein, T1、T2、T3Respectively thin-wall circular tube outer wall is axially x=0, x=L3/ 2, x=L3The temperature value at place, L3For thin-walled
The length of round tube.
7. thin-wall circular tube structure axial direction Calculation of Thermal Deformation method according to claim 6, which is characterized in that described utilize is divided
Cloth fiber-optic grating sensor monitors thermal force response signal suffered by thin-wall circular tube structure, and the response signal is substituted into Axial Thermal
Settlement calculation model calculates the axial direction thermal deformation of thin-wall circular tube structure and refers specifically to:
Based on the relational model of the thin-wall circular tube structure axial coordinate, thin-wall circular tube structure when without Thermal Load is established initially
Axial coordinates and arrangement thin-wall circular tube structure distribution formula optical fiber sensing network, by fiber arrangement on thin-wall circular tube outer surface edge
Axial x=0, x=L3/ 2, x=L3Place；Assuming that thin-wall circular tube initial temperature is T0, and record each optical fiber grating sensing at this temperature
Device initial center wavelength, thin-wall circular tube are placed in non-uniform temperature field；Under the effect of non-uniform temperature load, record is every respectively
Each fiber-optic grating sensor center wavelength shift amount is converted to the temperature of thin-wall circular tube by root fiber-optic grating sensor central wavelength
Axial Thermal settlement calculation model is substituted into, thin-wall circular tube structure axial direction thermal deformation is obtained are as follows:
Wherein,
Wherein, T1=Δ λ12/10+T0、T2=Δ λ13/10+T0、T3=Δ λ14/10+T0, Δ λ12、Δλ13、Δλ14For each grating
Center sensor wavelength shift, Δ L3For the thermal deformation length of thin-wall circular tube；
α is thin-wall circular tube structural material thermal expansion coefficient.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610476246.5A CN106156410B (en) | 2016-06-24 | 2016-06-24 | Thin-wall circular tube structure axial direction Calculation of Thermal Deformation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610476246.5A CN106156410B (en) | 2016-06-24 | 2016-06-24 | Thin-wall circular tube structure axial direction Calculation of Thermal Deformation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106156410A CN106156410A (en) | 2016-11-23 |
| CN106156410B true CN106156410B (en) | 2019-08-16 |
Family
ID=57349272
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610476246.5A Active CN106156410B (en) | 2016-06-24 | 2016-06-24 | Thin-wall circular tube structure axial direction Calculation of Thermal Deformation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106156410B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108519061B (en) * | 2018-03-29 | 2020-10-30 | 中原工学院 | Method and device for measuring deformation strain gradient of component |
| CN114370826B (en) * | 2021-12-21 | 2022-09-20 | 大连理工大学 | Measurement reference deviation correction method based on thermal expansion coefficient reverse calculation |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1384341A (en) * | 2002-06-14 | 2002-12-11 | 清华大学 | Optical-fiber grating sensor detecting pressure temperature simultaneously |
| WO2012038784A1 (en) * | 2010-09-20 | 2012-03-29 | Indian Institute Of Science | Strain and temperature discrimination using fiber bragg gratings in a cross-wire configuration |
| CN103217454A (en) * | 2013-03-22 | 2013-07-24 | 南京航空航天大学 | Fiber bragg grating measurement system and method for cylindrical structure thermal diffusivity |
| CN104019739A (en) * | 2014-05-08 | 2014-09-03 | 南京航空航天大学 | Plate structure impact load positioning method based on distributed fiber bragg grating sensing network |
| CN104567705A (en) * | 2014-12-19 | 2015-04-29 | 南京航空航天大学 | Strain and temperature aliasing signal decoupling method for optical fiber grating sensor under dynamic load |
| CN104834760A (en) * | 2014-02-12 | 2015-08-12 | 上海通用汽车有限公司 | Processing method of extracting deformation of each order of engine cylinder bore |
| CN105277130A (en) * | 2015-10-26 | 2016-01-27 | 南京航空航天大学 | Folded and unfolded structure fold area axial deformation optical fiber fast calculation method |
-
2016
- 2016-06-24 CN CN201610476246.5A patent/CN106156410B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1384341A (en) * | 2002-06-14 | 2002-12-11 | 清华大学 | Optical-fiber grating sensor detecting pressure temperature simultaneously |
| WO2012038784A1 (en) * | 2010-09-20 | 2012-03-29 | Indian Institute Of Science | Strain and temperature discrimination using fiber bragg gratings in a cross-wire configuration |
| CN103217454A (en) * | 2013-03-22 | 2013-07-24 | 南京航空航天大学 | Fiber bragg grating measurement system and method for cylindrical structure thermal diffusivity |
| CN104834760A (en) * | 2014-02-12 | 2015-08-12 | 上海通用汽车有限公司 | Processing method of extracting deformation of each order of engine cylinder bore |
| CN104019739A (en) * | 2014-05-08 | 2014-09-03 | 南京航空航天大学 | Plate structure impact load positioning method based on distributed fiber bragg grating sensing network |
| CN104567705A (en) * | 2014-12-19 | 2015-04-29 | 南京航空航天大学 | Strain and temperature aliasing signal decoupling method for optical fiber grating sensor under dynamic load |
| CN105277130A (en) * | 2015-10-26 | 2016-01-27 | 南京航空航天大学 | Folded and unfolded structure fold area axial deformation optical fiber fast calculation method |
Non-Patent Citations (5)
| Title |
|---|
| "光纤光栅传感网络的优化研究";王桂娜 等;《激光与红外》;20150131;第45卷(第1期);第66-69页 |
| "基于光纤传感技术的结构热属性与应变场监测研究";周雅斌;《中国优秀硕士学位论文全文数据库信息科技辑》;20150115;第I140-305页 |
| "机械零件热变形计算及分析";李倩龙;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20140815(第08期);第B022-346页 |
| "稳态均匀温度场中圆柱轴热变形研究";罗哉 等;《农业机械学报》;20070430;第38卷(第4期);第127-130页 |
| "非均匀温度场中厚壁圆筒变形研究";罗哉;《黑龙江科技学院学报》;20041130;第14卷(第6期);第331-333页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106156410A (en) | 2016-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Peng et al. | High precision strain monitoring for lithium ion batteries based on fiber Bragg grating sensors | |
| US9835506B2 (en) | Self-calibrated flow meter | |
| CN103453833B (en) | A kind of long gauge-length carbon fiber strain sensing devices and method of testing thereof | |
| CN106156410B (en) | Thin-wall circular tube structure axial direction Calculation of Thermal Deformation method | |
| EP3121586B1 (en) | Fatigue tester | |
| Najafi et al. | Real-time heat flux measurement using directional flame thermometer | |
| Nie et al. | Dew point measurement using dual quartz crystal resonator sensor | |
| Song et al. | Dynamic displacement prediction of beam structures using fiber bragg grating sensors | |
| CN105486242A (en) | Temperature compensation method for fiber grating strain sensor with package | |
| CN103217454B (en) | Fiber bragg grating measurement method for cylindrical structure thermal diffusivity | |
| CN106442727B (en) | A kind of method and system recognizing hard coat material mechanics parameters | |
| CN113483723B (en) | Passive excitation-based bridge structure strain monitoring system online calibration method | |
| Guzmán et al. | Determination of the constitutive relation parameters of a metallic material by measurement of temperature increment in compressive dynamic tests | |
| CN108953004B (en) | A kind of high-thrust rocket lateral force test predictor method | |
| CN102183540A (en) | Stress-free strain related analysis method in concrete stress strain monitoring | |
| Kalita et al. | Design and uncertainty evaluation of a strain measurement system | |
| CN108318529A (en) | For the temperature-compensation method of voltage detecting, electric field fingerprint detection method and system | |
| White et al. | The hysteresis characteristics of some industrial PRTs | |
| Zhang et al. | Development of strain measurement in superconducting magnet through fiber Bragg grating | |
| Rabin et al. | Infrared temperature sensing of mechanically loaded specimens: thermal analysis | |
| Najafi | Real-time heat flux estimation using filter based solutions for inverse heat conduction problems | |
| Cheng et al. | Calibrating static measurement data from distributed fiber optics by the integration of limited FBG sensors based on the extended kernel regression method | |
| CN105426660B (en) | A kind of Forecasting Methodology of machine components table planar residual stress field | |
| Chudzik | Measuring system with a dual needle probe for testing the parameters of heat-insulating materials | |
| CN114111579B (en) | Impact positioning method and system based on fiber bragg grating sensor under low sampling rate |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |