CN109974892A - Based on optical fiber temperature-measurement the temperature field forming method and system complementary with Digital Measurement of Temperature - Google Patents

Based on optical fiber temperature-measurement the temperature field forming method and system complementary with Digital Measurement of Temperature Download PDF

Info

Publication number
CN109974892A
CN109974892A CN201910229044.4A CN201910229044A CN109974892A CN 109974892 A CN109974892 A CN 109974892A CN 201910229044 A CN201910229044 A CN 201910229044A CN 109974892 A CN109974892 A CN 109974892A
Authority
CN
China
Prior art keywords
temperature
optical fiber
measurement
dam
measurement point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910229044.4A
Other languages
Chinese (zh)
Other versions
CN109974892B (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.)
Tsinghua University
Shanghai Investigation Design and Research Institute Co Ltd SIDRI
China Three Gorges Projects Development Co Ltd CTG
Original Assignee
Tsinghua University
Shanghai Investigation Design and Research Institute Co Ltd SIDRI
China Three Gorges Projects Development Co Ltd CTG
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 Tsinghua University, Shanghai Investigation Design and Research Institute Co Ltd SIDRI, China Three Gorges Projects Development Co Ltd CTG filed Critical Tsinghua University
Priority to CN201910229044.4A priority Critical patent/CN109974892B/en
Publication of CN109974892A publication Critical patent/CN109974892A/en
Application granted granted Critical
Publication of CN109974892B publication Critical patent/CN109974892B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/32Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K13/00Thermometers specially adapted for specific purposes

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)

Abstract

The temperature field forming method and system complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement that the invention discloses a kind of.Wherein, this method comprises: the space coordinate of optical fiber temperature-measurement point is obtained according to the whole coordinate of concrete dam and the buried line of optical fiber, wherein, concrete dam is embedded with a plurality of optical fiber, and optical fiber temperature-measurement point is on optical fiber for measuring the monitoring point of concrete dam temperature;The temperature value that the space coordinate that will acquire is measured with optical fiber temperature-measurement point is set as mapping one by one;According to the temperature value that optical fiber temperature-measurement point measures, it is distributed using the two-dimension temperature of plane where three dimensional finite element method optical fiber, is further distributed according to the two-dimension temperature of two neighboring optical fiber, calculates the distribution of three-dimensional temperature between two neighboring optical fiber plane.The present invention solves can not be according to the bulk temperature field distribution of optical fiber temperature-measurement degree data acquisition and dam coordinate phase mapping the technical issues of.

Description

Based on optical fiber temperature-measurement the temperature field forming method and system complementary with Digital Measurement of Temperature
Technical field
The present invention relates to construction engineering technical fields, in particular to complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement Temperature field forming method and system.
Background technique
Mass concrete is widely used in dam, Longspan Bridge and the isostructural main force part of skyscraper.Greatly Volume concrete one-time-concreting amount is big, and volume of concrete is thick, and thermal coefficient is relatively low, is poor conductor of heat, therefore it was hardened The heat overwhelming majority occurred in journey cannot dissipate, and be contained after inside concrete, concreting, due to the work of the heat of hydration With internal temperature increases, so as to cause concrete temperature raising, volume expansion.In During Initial Stage Construction, due to the production of hydration heat of cement Raw, internal temperature is more raised than external temperature fast, and volume of concrete expansion is big, to generate tensile stress on the surface of structure.? During post drop, due to the constraint by basis, periphery concrete and concrete internal temperature difference, and can be mixed Solidifying soil is internal to generate tensile stress.Since concrete is fragile material, tensile strength is low, when the drawing that concrete temperature variation generates is answered When power is greater than tensile strength, concrete just will appear crack, seriously endanger engineering safety.The large volume coagulation as caused by temperature The development in soil structure crack and crack is increasingly subject to the attention of the engineering circles personage such as building and water conservancy.
Large volume concrete structural temperature field is the important topic in design and construction.For concrete dam, big The construction period on dam will generally take measures to control the temperature of concrete, to reduce temperature gradient in dam body, prevent on Concrete cracking, it is ensured that dam safety.Therefore, the temperature field of acquisition dam entirety promptly and accurately, to be dam concrete Temperature control provide foundation, the Security Construction of concrete dam is had a very important significance.
Large volume Study on Temperature Field method is generally divided into two classes: numerical approximation, analytical Calculation.Numerical approximation is former according to calculating Reason can be divided into calculus of finite differences, boundary element method, FInite Element.(1) calculus of finite differences: suitable for some boundaries more rule, simply ask Topic.(2) boundary element method: along computational domain boundary discrete method, reduce the dimension of problem, the given data needed is few.(3) finite element The temperature field solved under certain boundary value condition: being converted into the functional problem that solves by method, first computational domain it is discrete be limited Unit uses certain interpolating function in unit, establishes the system of linear equations of unit junction temperature, then solving equations find out node Temperature value.
The Wilson's of nineteen sixty-eight California, USA university are that US Army engineer team develops analog large volume coagulation The two-dimensional finite metaprogram DOT-DICE in soil structure stage construction temperature field;Nineteen eighty-two, the engineer of US Army engineering division S.B.Tatro and E.K.Schrader has made modification to DOT-DICE program;1992, Barret (P.K.Barrett) etc., Three dimensional temperature Stress calculation software ANACAP is described, creativeness is cracking model to be introduced into the analysis of temperature stress In;1973, self-editing the having of first, China concrete temperature Creep Stress of Chinese Water Conservancy water power academy of sciences Zhu Baifang academician Metaprogram is limited, China's simulation calculation of mass concrete temperature stress for the first time in history is realized;Wu Feng et al. secondary development ANSYS analyzes the temperature crack control of concrete structure;Zhang Ziming et al. calculates temperature using crack band model and splits Seam;Liu Xinghong et al. has carried out simulation calculation to the thermal cracking process of concrete using non-mesh method;Chen Hui et al. is utilized 3D solid unit in software ANSYS ANSYS has carried out thermosetting coupling analysis to beams of concrete;NAKAMURE Et al. propose the boundary element method of thermal parameters inverting.
The monitoring of traditional dam body temperature generally uses point type thermometer, has the disadvantage that (1) thermometer is only surveyed Measure the temperature of a point;(2) tight to operating environment requirements, poor anti jamming capability, installation complex jamming construction;(3) thermometer is believed Breath amount is very little, is difficult to grasp the variation of entire dam body internal temperature field.
Chinese patent 201810308038.3 discloses a kind of concrete dam distributed temperature measuring optical fiber bifilar " zigzag " and leads to storehouse Burial and construction method, this method is by designing bifilar embedded, " zigzag " laying and the embedded distributed temperature measuring promoted of odd-even alternation Optical fiber buries technique and method, it is ensured that the relatively high-survival rate of distributed temperature measuring optical fiber, is more preferably applied more monitoring data amounts Work technique and measurement and computational accuracy.But the advantage for how combining optical fiber temperature-measurement do not suggested that, is obtained corresponding with dam body coordinate Bulk temperature field distribution.
For above-mentioned problem, currently no effective solution has been proposed.
Summary of the invention
The embodiment of the invention provides the temperature field forming method and system complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement, with The technical issues of at least solution can not be according to the bulk temperature field distribution of optical fiber temperature-measurement degree data acquisition and dam coordinate phase mapping.
According to an aspect of an embodiment of the present invention, a kind of temperature complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement is provided Field forming method, comprising: according to the whole coordinate of concrete dam and the buried line of optical fiber, obtain the sky of optical fiber temperature-measurement point Between coordinate, wherein the concrete dam is embedded with a plurality of optical fiber, and the optical fiber temperature-measurement point is on the optical fiber for surveying Measure the monitoring point of the concrete dam temperature;The temperature value that the space coordinate and the optical fiber temperature-measurement point that will acquire measure is arranged To map one by one;According to the temperature value that the optical fiber temperature-measurement point measures, optical fiber where calculating the optical fiber using finite difference calculus The two-dimension temperature of plane is distributed;It is distributed according to the two-dimension temperature of two neighboring optical fiber plane, calculates the two neighboring light Distribution of three-dimensional temperature between fine plane.
Further, according to the whole coordinate of concrete dam and the buried line of optical fiber, optical fiber temperature-measurement point is obtained Space coordinate includes: to obtain optical fiber inbuilt route in the concrete dam;Multiple monitorings are taken in the route Point is used as the optical fiber temperature-measurement point;It is sat according to the space that the whole coordinate of the concrete dam obtains the optical fiber temperature-measurement point Mark.
Further, the temperature value measured according to the optical fiber temperature-measurement point calculates the optical fiber institute using finite difference calculus It include: that the plane where the optical fiber is divided into multiple triangular units in the two-dimension temperature distribution of optical fiber plane, wherein described The node of triangular unit is the optical fiber temperature-measurement point;The two dimension in the triangular unit is calculated using the finite difference calculus Temperature Distribution.
Further, it is distributed according to the two-dimension temperature of two neighboring optical fiber plane, calculates the two neighboring optical fiber Distribution of three-dimensional temperature between plane includes: acquisition curve side element, the curve side element by the two neighboring optical fiber plane, The upstream dam facing of the concrete dam, the downstream dam facing of the concrete dam and the concrete dam are contacted with air Two faces surround;Utilize the distribution of three-dimensional temperature of curve side element described in three dimensional finite element method.
Further, described before the distribution of three-dimensional temperature using curve side element described in three dimensional finite element method Method further include: the curve side element is converted to by regular hexahedron unit using coordinate transformation method.
According to another aspect of an embodiment of the present invention, a kind of temperature complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement is additionally provided It spends field and forms system, comprising: acquiring unit, for obtaining according to the whole coordinate of concrete dam and the buried line of optical fiber Take the space coordinate of optical fiber temperature-measurement point, wherein the concrete dam is embedded with a plurality of optical fiber, and the optical fiber temperature-measurement point is For measuring the monitoring point of the concrete dam temperature on the optical fiber;Setting unit, space coordinate for will acquire with The temperature value that the optical fiber temperature-measurement point measures is set as mapping one by one;First computing unit, for according to the optical fiber temperature-measurement point The temperature value measured, the two-dimension temperature of optical fiber plane is distributed where calculating the optical fiber using finite difference calculus;Second calculates list Member calculates between the two neighboring optical fiber plane for being distributed according to the two-dimension temperature of two neighboring optical fiber plane Distribution of three-dimensional temperature.
Further, the acquiring unit includes: the first acquisition module, for obtaining the fibre in the concrete dam In inbuilt route;Selecting module, for selecting multiple monitoring points as the optical fiber temperature-measurement point in the route;Second obtains Modulus block obtains the space coordinate of the optical fiber temperature-measurement point for the whole coordinate according to the concrete dam.
Further, first computing unit includes: division module, for the plane where the optical fiber to be divided into Multiple triangular units, wherein the node of the triangular unit is the optical fiber temperature-measurement point;First computing module, for utilizing institute The two-dimension temperature that finite difference calculus calculates in the triangular unit is stated to be distributed.
Further, first computing unit includes: that third obtains module, for obtaining curve side element, the curl Unit by the two neighboring optical fiber plane, the upstream dam facing of the concrete dam, the concrete dam downstream dam facing And two faces that the concrete dam is contacted with air surround;Second computing module, based on using three-dimensional finite element method Calculate the distribution of three-dimensional temperature of the curve side element.
Further, the system also includes conversion modules, for being converted the curve side element using coordinate transformation method Be positive hexahedral element.
For the prior art in temperature control, what is generallyd use is all the mode of Level temperature difference control, is guaranteed any in plane The temperature difference of point-to-point transmission is in safe range.Such control mode, control range is limited to plane, and precision is not high, can not The temperature change of completely true reflection mass concrete.The temperature complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement provided by the invention Field forming method and system are spent, using the three-dimensional temperature field reconstructed, arbitrary neighborhood optical fiber can be buried mixed between plane Solidifying soil carries out temperature difference control.In conjunction with intelligent water communication temperature control system, it is ensured that in three-dimensional space, the temperature difference dimension between any two points It holds in safe range.Original temperature is controlled, three-dimensional space is extended to by two-dimensional surface, improves temperature controlled peace Quan Xing, science are avoided since the excessive caused temperature stress of local temperature difference is excessive, are avoided concrete and crack occur, mention The high safety of engineering.
In embodiments of the present invention, it using according to the whole coordinate of concrete dam and the buried line of optical fiber, obtains The space coordinate of optical fiber temperature-measurement point, wherein the concrete dam is embedded with a plurality of optical fiber, and the optical fiber temperature-measurement point is institute It states on optical fiber for measuring the monitoring point of the concrete dam temperature;The space coordinate that will acquire and the optical fiber temperature-measurement point are surveyed The temperature value obtained is set as mapping one by one;According to the temperature value that the optical fiber temperature-measurement point measures, institute is calculated using finite difference calculus The two-dimension temperature distribution of optical fiber plane where stating optical fiber;It is distributed, is calculated according to the two-dimension temperature of two neighboring optical fiber plane The mode of distribution of three-dimensional temperature between the two neighboring optical fiber plane, measures concrete dam by distribution type fiber-optic Monitoring point temperature, and finite difference calculus is combined, achieve the purpose that monitor dam temperature field by optical fiber temperature-measurement, to realize The technical effect of the corresponding bulk temperature field distribution of dam body coordinate is obtained using optical fiber temperature-measurement, and then solving can not be according to light The technical issues of bulk temperature field distribution of fine testing temperature data acquisition and dam coordinate phase mapping.
Detailed description of the invention
The drawings described herein are used to provide a further understanding of the present invention, constitutes part of this application, this hair Bright illustrative embodiments and their description are used to explain the present invention, and are not constituted improper limitations of the present invention.In the accompanying drawings:
Fig. 1 is a kind of temperature field forming method complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to an embodiment of the present invention Flow chart;
Fig. 2 (a) is the schematic diagram according to an embodiment of the present invention that optical fiber temperature-measurement point is arranged in fibre circuit;
Fig. 2 (b) is the schematic diagram according to an embodiment of the present invention that optical fiber temperature-measurement point is arranged in fibre circuit;
Fig. 3 is the schematic diagram of triangular unit according to an embodiment of the present invention;
Fig. 4 is the schematic diagram of curve side element according to an embodiment of the present invention;
Fig. 5 is the schematic diagram of regular hexahedron unit of the curve side element according to an embodiment of the present invention after coordinate transform;
Fig. 6 is that a kind of temperature field complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to an embodiment of the present invention forms system Schematic diagram.
Specific embodiment
In order to enable those skilled in the art to better understand the solution of the present invention, below in conjunction in the embodiment of the present invention Attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is only The embodiment of a part of the invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people The model that the present invention protects all should belong in member's every other embodiment obtained without making creative work It encloses.
It should be noted that description and claims of this specification and term " first " in above-mentioned attached drawing, " Two " etc. be to be used to distinguish similar objects, without being used to describe a particular order or precedence order.It should be understood that using in this way Data be interchangeable under appropriate circumstances, so as to the embodiment of the present invention described herein can in addition to illustrating herein or Sequence other than those of description is implemented.In addition, term " includes " and " having " and their any deformation, it is intended that cover Cover it is non-exclusive include, for example, the process, method, system, product or equipment for containing a series of steps or units are not necessarily limited to Step or unit those of is clearly listed, but may include be not clearly listed or for these process, methods, product Or other step or units that equipment is intrinsic.
Firstly, several nouns involved by the embodiment of the present application are introduced:
1, finite difference calculus: finite difference method is exactly a kind of numerical solution, its basic thought is that first problem is determined Adopted domain carries out mesh generation and difference quotient in Definite problem is changed into difference by Numerical Differentiation Formulae appropriate then on mesh point Quotient so that problem is discrete turns to difference scheme former, and then finds out numerical solution.
2, temperature field: the collection of temperature is collectively referred to as temperature field on each point in material system.It is time and space coordinate Function reflects distribution of the temperature on room and time.
3, upstream dam facing: the side of dam high water level is upstream, and the upstream face of dam upstream side is the upstream dam of dam Face.
4, FInite Element: being a kind of high-effect, common numerical computation method, and principle is that continuous solution domain is discrete For the assembly of one group of unit, with the approximate function assumed in each unit come fragment indicate to solve it is to be asked unknown on domain Field function, approximate function are usually expressed by unknown field function and its derivative in the numerical interpolation function of each node of unit.To A continuous infinite degrees of freedom problem is set to become discrete finite degrees of freedom problem.
According to embodiments of the present invention, a kind of embodiment of the method for concrete dam reconstructing temperature field method is provided, is needed Illustrate, step shown in the flowchart of the accompanying drawings can be in a computer system such as a set of computer executable instructions It executes, although also, logical order is shown in flow charts, and it in some cases, can be to be different from herein suitable Sequence executes shown or described step.
Fig. 1 is a kind of temperature field forming method complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to an embodiment of the present invention Flow chart, as shown in Figure 1, this method comprises the following steps:
Step S102 obtains optical fiber temperature-measurement point according to the whole coordinate of concrete dam and the buried line of optical fiber Space coordinate, wherein concrete dam is embedded with a plurality of optical fiber, and optical fiber temperature-measurement point is on optical fiber for measuring concrete dam temperature The monitoring point of degree;
Step S104, the temperature value that the space coordinate that will acquire is measured with optical fiber temperature-measurement point are set as mapping one by one;One It, can be by writing program in a optional embodiment, the temperature data for measuring each optical fiber temperature-measurement point and each optical fiber temperature-measurement point Space coordinate map one by one;
Step S106, according to the temperature value that optical fiber temperature-measurement point measures, optical fiber is flat where calculating optical fiber using finite difference calculus The two-dimension temperature in face is distributed;
Step S108 is distributed according to the two-dimension temperature of two neighboring optical fiber plane, calculates between two neighboring optical fiber plane Distribution of three-dimensional temperature.
Optionally, after obtaining distribution of three-dimensional temperature, this method further include: using thermometer measure temperature with again The three-dimensional temperature field of structure compares and analyzes.
In embodiments of the present invention, the temperature value that optical fiber temperature-measurement point measures is mapped one by one with actual spatial value, Obtain the point temperature in mass concrete.Line temperature is formed along the embedded path of optical fiber using a temperature.Utilize line temperature Degree is buried plane to optical fiber and is divided, using line temperature as boundary value.The face that optical fiber buries plane is obtained using finite difference Temperature.Finally, reconstructing with FInite Element as boundary value and obtaining three-dimensional temperature the face temperature that adjacent fiber buries plane. It is reconstructed in this way, completing optical fiber temperature-measurement by putting to line, by line to face, by the three-dimensional temperature field in face to body.
In embodiments of the present invention, it using according to the whole coordinate of concrete dam and the buried line of optical fiber, obtains The space coordinate of optical fiber temperature-measurement point, wherein concrete dam is embedded with a plurality of optical fiber, and optical fiber temperature-measurement point is on optical fiber for measuring The monitoring point of concrete dam temperature;The temperature value that the space coordinate that will acquire is measured with optical fiber temperature-measurement point is set as reflecting one by one It penetrates;According to the temperature value that optical fiber temperature-measurement point measures, the two-dimension temperature point of optical fiber plane where calculating optical fiber using finite difference calculus Cloth;It is distributed according to the two-dimension temperature of two neighboring optical fiber plane, calculates the three-dimensional temperature field point between two neighboring optical fiber plane The mode of cloth is measured the monitoring point temperature of concrete dam by distribution type fiber-optic, and combines finite difference calculus, has been reached and has been passed through Optical fiber temperature-measurement monitors the purpose in dam temperature field, obtains the corresponding whole temperature of dam body coordinate using optical fiber temperature-measurement to realize The technical effect of field distribution is spent, and then solving can not be according to the entirety of optical fiber temperature-measurement degree data acquisition and dam coordinate phase mapping The technical issues of thermo parameters method.
Optionally, according to the whole coordinate of concrete dam and the buried line of optical fiber, the sky of optical fiber temperature-measurement point is obtained Between coordinate include: obtain optical fiber inbuilt route in concrete dam;Take multiple monitoring points as optical fiber temperature-measurement in the line Point;The space coordinate of optical fiber temperature-measurement point is obtained according to the whole coordinate of concrete dam.
In an optional embodiment, the buried line of optical fiber is recorded in CAD design figure, according to optical fiber in coagulation The CAD design figure and actual monitoring demand of buried line, take multiple monitoring points as light in CAD design figure in native dam Fine temperature measuring point, and according to the whole from mark of dam, extract the space coordinate of each optical fiber temperature-measurement point.Optionally, it can incite somebody to action Optical fiber buried line in CAD utilizes " spacing equal part " function, makes to generate a point every fixed range on route, by the point As optical fiber temperature-measurement point.Such as: a point can be generated every 1m and extracts the space coordinate of each point, such as in CAD design figure Shown in Fig. 2 (a), Fig. 2 (b).Optionally, during optical fiber is embedded, using the linear heat detector of distribution type fiber-optic (DTS), It monitors on fibre circuit, the temperature data at optical fiber temperature-measurement point at the position 1m.
Since optical fiber runs through entire volume of concrete, the point temperature data obtained by optical fiber temperature-measurement point will have quantity The features such as more, regularity of distribution, precision are high, timing updates, and the thus obtained line temperature for burying path along optical fiber is true and reliable. In the inbuilt plane of optical fiber, optical fiber is buried into path as boundary condition, using the method for finite difference to the two-dimentional warm of plane Degree field is reconstructed.Due to optical fiber bury the temperature in path it is known that so, reconstruction plane is all First Boundary Condition, and side Boundary's shape levels off to quadrangle, and boundary condition is simple, and the two-dimensional temperature field precision that finite difference obtains is high.Utilize adjacent two dimension Temperature field selects hexahedral element to divide concrete, with increasing for element number, precision using FInite Element It correspondinglys increase.Accordingly, the three-dimensional temperature field of the concrete obtained is true and reliable, and more traditional reconstructing temperature field method is more smart Really, and then the Thermal Stress Field of acquisition is also more accurate.
In an alternative embodiment, optical fiber is spaced at regular intervals, acquires a temperature data, by distributed light Temperature data along fibre is shown in mass concrete by true spatial position, and then is calculated at regular intervals Dam thermo parameters method, calculated thermo parameters method situation obtain the Temperature Distribution feelings of different monoliths, different moments Condition.Time interval can self-setting as needed, as Fig. 2 (a) show the temperature of the 15:00 collecting fiber on the 22nd of May in 2018 Data, Fig. 2 (b) show the temperature data of the May in 2018 of 15:05 collecting fiber on the 22nd.
Optical fiber acquires a temperature data at regular intervals, carries out a data more according to freshly harvested temperature data Newly, and then update reconstruct three-dimensional temperature field, the temperature field guaranteed be in the construction process it is continuous at any time, by continuous Temperature field, it can be seen that the heat dispersal situations of concrete, temperature change, to grasp more true concrete three-dimensional temperature field Variation.On the one hand, it is changed with time situation according to three-dimensional temperature field, the Thermal Stress Field changed over time can be reconstructed, And then analyze and obtain the transmission process of concrete internal temperature, grasp true temperature stress distribution, more accurate, safe pair The exothermic process of concrete is monitored.On the other hand, temperature field is updated with the update of fiber data, it can be seen that is poured The heat release situation of concrete in the process, and according to the variation in temperature field, water flowing is controlled, while also can control the pouring of concrete, The variation in temperature field mutually coordinated with temperature change by be allowed pouring progress, guarantee that concrete pours as one of the condition that control is poured Filling is not in the excessively high situation of temperature, and more scientific and reasonable instructs site operation, makes the pouring progress of concrete more equal Weighing apparatus.
Optionally, the temperature value measured according to optical fiber temperature-measurement point, optical fiber plane where calculating optical fiber using finite difference calculus Two-dimension temperature distribution include: that the plane where optical fiber is divided into multiple triangular units, wherein the node of triangular unit be light Fine temperature measuring point;The two-dimension temperature in triangular unit is calculated using finite difference calculus to be distributed.
In an optional embodiment, according to the coordinate value that step S102 is got, plane where optical fiber is regarded as It is two dimensions and unstable temperature field, calculates temperature value using the method for finite difference, and then obtain two-dimension temperature distribution.Optionally, Difference unit selects triangular element, i.e., the plane where optical fiber is divided into many triangular elements, wherein each triangle list All there are three nodes for member.Not such as Fig. 3 institute, coordinate system x, the y table not plane where optical fiber, the node of triangular unit e are i, j and m, For any 3 optical fiber temperature-measurements point on optical fiber, coordinate is respectively (xi, yi), (xj, yj) and (xm, ym).If in triangular unit e Temperature be coordinate linear function, i.e.,
Tj e(x, y)=γ12xj3yj,
3 points of i, j and m of coordinate is substituted into above formula to obtain:
Ti e(x, y)=γ1+y2xi3yi
Tj e(x, y)=γ12xj3yj
Tj e(x, y)=γ12xj3yj
Simultaneous solution above formula can find out coefficient gamma1, γ2, γ3, and then it is available,
Te(x, y)=NiTi+NjTj+NmTm, wherein Ni、Nj、NmFor shape function,
ai=xiym-xmyi, bi=yi-ym, ci=xm-xjOptionally, according to the two-dimension temperature of two neighboring optical fiber plane point Cloth, calculating the distribution of three-dimensional temperature between two neighboring optical fiber plane includes: acquisition curve side element, and curve side element is by adjacent two A optical fiber plane, the upstream dam facing of concrete dam, the downstream dam facing of concrete dam and concrete dam are contacted with air Two faces surround;Utilize the distribution of three-dimensional temperature of three dimensional finite element method curve side element.
According to heat transfer theory, three dimensions temperature field T (x, y, z) should meet following partial differential equation and it is corresponding just Beginning condition and boundary condition:
Partial differential equation:
Primary condition T |τ=0=T0(x, y, z)
Boundary condition: T=Tb, wherein TbFor the temperature value on six faces of curve side element, wherein the six of curve side element Face includes two adjacent optical fiber planes, upstream dam facing, two faces that downstream dam facing and concrete dam are contacted with air.
The temperature of any point can use shape function N in curve side elementiWith unit junction temperature interpolation representation are as follows:
And then it is available,
For the three-dimensional space domain R of curve side element, obtained by weighted residual method
And then using gal the Liao Dynasty golden method weighting function W in spatial domainiEqual to shape function Ni, it is available:
It is rightIntegration by parts is carried out to obtain:
By formulaWrite as the form of matrix:
(i= 1,2 ..., 8) it, sums to all units, and is included in boundary condition:It enables
Then obtain:Discretization is carried out in time-domain to it, using linear interpolation function, In 0≤τ of time-domain≤Δ τ, junction temperature { T } can be indicated are as follows:
Wherein, N0(τ) N1(τ) is the shape function in time-domain,
Due toSo the time-derivative of junction temperature are as follows:
Initial node temperature { T }0It is known, when that ask is τ=Δ τ junction temperature { T }1, by the golden formula of gal the Liao Dynasty, Take time-domain weight function W1(τ)=N1(τ), by:
It may further obtain:
Time τ is integrated, and abbreviation obtains:
Equally { P } is indicated are as follows:
{P}0{ P }1Respectively indicate τ=0 and τ=Δ τ moment { P } value, then:
ByAvailable solution is non-steady The equation for determining temperature field is as follows:
The equation is solved, the distribution of three-dimensional temperature of available curve side element.
Optionally, before the distribution of three-dimensional temperature using three dimensional finite element method curve side element, method further include: Curve side element is converted to by regular hexahedron unit using coordinate transformation method.
Any curve side element as shown in Figure 4, grid dividing are not limited by boundary shape, and cell size can not phase Deng being that a kind of precision is high and widely used unit.However directly carry out unit analysis to it and be difficult, this is because it Irregular, the ununified shape of geometry, press different formulas calculating one by one to each unit, it is difficult due to its heavy workload To carry out.
In an alternative embodiment, using the method for coordinate transform, any curve side element in xyz coordinate system is turned The regular hexahedron unit being changed in another coordinate system xhz, optionally, can this regular hexahedron unit be referred to as basic unit or mother Any curve side element is considered as the image of basic unit, referred to as actual cell or subelement by unit.8 nodes as shown in Figure 5 are just Hexahedral element is female unit, and establishing origin can be obtained in the local coordinate system (ξ, η, ζ) of unit centroid by coordinate transform 8 node isoparametric elements of space, coordinate conversion relation are as follows:
The displacement function of unit are as follows:
Wherein, μi, vi, ωiAnd xi, yi, ziThe respectively actual displacement of node i and seat Scale value.
Element displacement function is indicated with matrix form are as follows:
Wherein, { δi}=[μi vi ωi]T(i=1,2 ..., 8) is displacement of joint array;{δ}e=[{ δ1} {δ2} … {δ8}]TFor the node position of entire unit Move array.The unified shape function expression formula of 8 nodes are as follows:
Wherein, ξ1、η1、ζ1It is node i in part Coordinate in coordinate system (ξ, η, ζ).
Derivative of the shape function to local coordinate are as follows:
The geometric equation of space problem are as follows:Wherein, [B], { δ }eRespectively unit Geometric matrix and unit displacement of joint matrix, the wherein geometric matrix [B] of unit are as follows:
It can be obtained by compound function derivation rule:
Matrix is cited as the three-dimensional Jacobian matrix of coordinate transform:
Its inverse matrix [J]-1It can press Formula obtains:
According to another aspect of an embodiment of the present invention, a kind of temperature complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement is additionally provided It spends field and forms system, as shown in Figure 6, comprising: acquiring unit 601, for the whole coordinate and optical fiber according to concrete dam Buried line, obtain optical fiber temperature-measurement point space coordinate, wherein concrete dam is embedded with a plurality of optical fiber, optical fiber temperature-measurement point For the monitoring point for being used to measure concrete dam temperature on optical fiber;Setting unit 602, space coordinate and optical fiber for will acquire The temperature value that temperature measuring point measures is set as mapping one by one;First computing unit 603, the temperature for being measured according to optical fiber temperature-measurement point Value, the two-dimension temperature of optical fiber plane is distributed where calculating optical fiber using finite difference calculus;Second computing unit 604 is used for basis The two-dimension temperature of two neighboring optical fiber plane is distributed, and calculates the distribution of three-dimensional temperature between two neighboring optical fiber plane.
In embodiments of the present invention, using acquiring unit, for according to the whole coordinate of concrete dam and optical fiber Buried line obtains the space coordinate of optical fiber temperature-measurement point, wherein concrete dam is embedded with a plurality of optical fiber, and optical fiber temperature-measurement point is For measuring the monitoring point of concrete dam temperature on optical fiber;Setting unit, space coordinate and optical fiber temperature-measurement for will acquire The temperature value that point measures is set as mapping one by one;First computing unit, the temperature value for being measured according to optical fiber temperature-measurement point utilize The two-dimension temperature distribution of optical fiber plane where finite difference calculus calculates optical fiber;Second computing unit, for according to two neighboring light The two-dimension temperature of fine plane is distributed, and is calculated the mode of the distribution of three-dimensional temperature between two neighboring optical fiber plane, is passed through distribution Formula optical fiber measures the monitoring point temperature of concrete dam, and combines finite difference calculus, has reached and has monitored dam by optical fiber temperature-measurement The purpose in temperature field is imitated to realize using the technology that optical fiber temperature-measurement obtains the corresponding bulk temperature field distribution of dam body coordinate Fruit, and then solving can not be according to the technology of optical fiber temperature-measurement degree data acquisition and the bulk temperature field distribution of dam coordinate phase mapping Problem.
Optionally, acquiring unit includes: the first acquisition module, for obtaining fine route inbuilt in concrete dam; Selecting module, for selecting multiple monitoring points as optical fiber temperature-measurement point in the line;Second obtains module, for according to concrete The whole coordinate of dam obtains the space coordinate of optical fiber temperature-measurement point.
In an optional embodiment, the buried line of optical fiber is recorded in CAD design figure, according to optical fiber in coagulation The CAD design figure and actual monitoring demand of buried line, take multiple monitoring points as light in CAD design figure in native dam Fine temperature measuring point, and according to the whole from mark of dam, extract the space coordinate of each optical fiber temperature-measurement point.Optionally, it can incite somebody to action Optical fiber buried line in CAD utilizes " spacing equal part " function, makes to generate a point every fixed range on route, by the point As optical fiber temperature-measurement point.Such as: a point can be generated every 1m and extracts the space coordinate of each point, such as in CAD design figure Shown in Fig. 2 (a) and Fig. 2 (b).Optionally, during optical fiber is embedded, using the linear heat detector of distribution type fiber-optic (DTS), It monitors on fibre circuit, the temperature data at optical fiber temperature-measurement point at the position 1m.In an alternative embodiment, optical fiber It is spaced at regular intervals, acquires a temperature data, by the temperature data along distribution type fiber-optic, by true spatial position It is shown in mass concrete, and then calculates dam thermo parameters method at regular intervals, calculated temperature Field distribution situation obtains the profiling temperatures of different monoliths, different moments.Time interval can self-setting as needed, such as Fig. 2 (a) show the temperature data of the 15:00 collecting fiber on the 22nd of May in 2018, and Fig. 2 (b) show on May 22nd, 2018 The temperature data of 15:05 collecting fiber.
Optionally, the first computing unit includes: division module, for the plane where optical fiber to be divided into multiple triangle lists Member, wherein the node of triangular unit is optical fiber temperature-measurement point;First computing module, for calculating triangle list using finite difference calculus Two-dimension temperature distribution in member.
In an optional embodiment, according to the coordinate value that step S102 is got, plane where optical fiber is regarded as It is two dimensions and unstable temperature field, calculates temperature value using the method for finite difference, and then obtain two-dimension temperature distribution.Optionally, Difference unit selects triangular element, i.e., the plane where optical fiber is divided into many triangular elements, wherein each triangle list All there are three nodes for member.As shown in figure 3, coordinate system x, y indicate the plane where optical fiber, the node of triangular unit e is i, j and m, For any 3 optical fiber temperature-measurements point on optical fiber, coordinate is respectively (xi, yi), (xj, yj) and (xm, ym).If in triangular unit e Temperature be coordinate linear function, i.e.,
Tj e(x, y)=γ1+y2xj3yj,
3 points of i, j and m of coordinate is substituted into above formula to obtain:
Ti e(x, y)=γ12xi3yi
Tj e(x, y)=γ12xj3yj
Tj e(x, y)=γ12xj3yj
Simultaneous solution above formula can find out coefficient gamma1, γ2, γ3, and then it is available,
Te(x, y)=NiTi+NjTj+NmTm, wherein Ni、Nj、NmFor shape function,
ai=xiym-xmyi, bi=yi-ym, ci=xm-xj,
A is the area of triangular unit e, likewise, also available aj, bj, cj, am, bm, cm, can also further obtain
According to shape function NiThe characteristics of, there is also following conditions:
Therefore, for triangular unit e, have
And then obtain,
Solve system of equation:
Wherein,
For all units related with node i and,
It is τ=τ that it, which is solved,n+1When junction temperature { Tn+1}。
Optionally, the first computing unit includes: that third obtains module, and for obtaining curve side element, curve side element is by adjacent Two optical fiber planes, the upstream dam facing of concrete dam, the downstream dam facing of concrete dam and concrete dam and air connect Two faces of touching surround;Second computing module, for the distribution of three-dimensional temperature using three dimensional finite element method curve side element.
According to heat transfer theory, three dimensions temperature field T (x, y, z) should meet following partial differential equation and it is corresponding just Beginning condition and boundary condition:
Partial differential equation:
Primary condition: T |τ=0=T0(x, y, z)
Boundary condition: T=Tb, wherein TbFor the temperature value on six faces of curve side element, wherein the six of curve side element Face includes two adjacent optical fiber planes, upstream dam facing, two faces that downstream dam facing and concrete dam are contacted with air.
The temperature of any point can use shape function N in curve side elementiWith unit junction temperature interpolation representation are as follows:
And then it is available,
For the three-dimensional space domain R of curve side element, obtained by weighted residual method
And then using gal the Liao Dynasty golden method weighting function W in spatial domainiEqual to shape function Ni, it is available:
It is rightIntegration by parts is carried out to obtain:
By formulaWrite as the form of matrix:
(i=1, 2 ..., 8) it, sums to all units, and is included in boundary condition:It enables
Then obtain:Discretization is carried out in time-domain to it, using linear interpolation function, In 0≤τ of time-domain≤Δ τ, junction temperature { T } can be indicated are as follows:
Wherein, N0(τ) N1(τ) is the shape function in time-domain,
Due toSo the time-derivative of junction temperature are as follows:
Initial node temperature { T }0It is known, when that ask is τ=Δ τ junction temperature { T }1, by the golden formula of gal the Liao Dynasty, Take time-domain weight function W1(τ)=N1(τ), by:
It may further obtain:
Time τ τ is integrated, and abbreviation obtains:
Equally { P } is indicated are as follows:
{P}0{ P }1Respectively indicate τ=0 and τ=Δ τ moment { P } value, then:
ByAvailable solution is non-steady The equation for determining temperature field is as follows:
Wherein,
{T}0={ T (τ0), { T }1={ T (τ0+Δτ)}
{P}0={ P (τ0), { P }1={ P (τ0+Δτ)}
The equation is solved, the distribution of three-dimensional temperature of available curve side element.
Optionally, system further include: conversion module, for curve side element to be converted to regular hexahedron using coordinate transformation method Unit.
Any curve side element as shown in Figure 4, grid dividing are not limited by boundary shape, and cell size can not phase Deng being that a kind of precision is high and widely used unit.However directly carry out unit analysis to it and be difficult, this is because it Irregular, the ununified shape of geometry, press different formulas calculating one by one to each unit, it is difficult due to its heavy workload To carry out.
In an alternative embodiment, using the method for coordinate transform, any curve side element in xyz coordinate system is turned The regular hexahedron unit being changed in another coordinate system xhz, optionally, can this regular hexahedron unit be referred to as basic unit or mother Any curve side element is considered as the image of basic unit, referred to as actual cell or subelement by unit.8 nodes as shown in Figure 5 are just Hexahedral element is female unit, and establishing origin can be obtained in the local coordinate system (ξ, η, ζ) of unit centroid by coordinate transform 8 node isoparametric elements of space, coordinate conversion relation are as follows:
The displacement function of unit are as follows:
Wherein, μi, vi, ωiAnd xi, yi, ziThe respectively actual displacement of node i and seat Scale value.
Element displacement function is indicated with matrix form are as follows: Wherein, { δi}=[μi vi ωi]T(i=1,2 ..., 8) is displacement of joint array;{δ}e= [{δ1} {δ2} … {δ8}]TFor the displacement of joint array of entire unit.The unified shape function expression formula of 8 nodes are as follows:
Wherein, ξ1、η1、ζ1It is node i in part Coordinate in coordinate system (ξ, η, ζ).
Derivative of the shape function to local coordinate are as follows:
The geometric equation of space problem are as follows:Wherein, [B], { δ }eRespectively unit Geometric matrix and unit displacement of joint matrix, the wherein geometric matrix [B] of unit are as follows:
It can be obtained by compound function derivation rule:
Matrix [J] is the three-dimensional Jacobian matrix of coordinate transform:
Its inverse matrix [J]-1It can press Formula obtains:
The serial number of the above embodiments of the invention is only for description, does not represent the advantages or disadvantages of the embodiments.
In the above embodiment of the invention, it all emphasizes particularly on different fields to the description of each embodiment, does not have in some embodiment The part of detailed description, reference can be made to the related descriptions of other embodiments.
In several embodiments provided herein, it should be understood that disclosed technology contents can pass through others Mode is realized.Wherein, system embodiment described above is only schematical, such as the division of the unit, Ke Yiwei A kind of logical function partition, there may be another division manner in actual implementation, for example, multiple units or components can combine or Person is desirably integrated into another system, or some features can be ignored or not executed.Another point, shown or discussed is mutual Between coupling, direct-coupling or communication connection can be through some interfaces, the INDIRECT COUPLING or communication link of unit or module It connects, can be electrical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple On unit.It can some or all of the units may be selected to achieve the purpose of the solution of this embodiment according to the actual needs.
It, can also be in addition, the functional units in various embodiments of the present invention may be integrated into one processing unit It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list Member both can take the form of hardware realization, can also realize in the form of software functional units.
If the integrated unit is realized in the form of SFU software functional unit and sells or use as independent product When, it can store in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially The all or part of the part that contributes to existing technology or the technical solution can be in the form of software products in other words It embodies, which is stored in a storage medium, including some instructions are used so that a computer Equipment (can for personal computer, server or network equipment etc.) execute each embodiment the method for the present invention whole or Part steps.And storage medium above-mentioned includes: that USB flash disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited Reservoir (RAM, Random Access Memory), mobile hard disk, magnetic or disk etc. be various to can store program code Medium.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered It is considered as protection scope of the present invention.

Claims (10)

1. a kind of temperature field forming method complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement characterized by comprising
According to the whole coordinate of concrete dam and the buried line of optical fiber, the space coordinate of optical fiber temperature-measurement point is obtained, wherein The concrete dam is embedded with a plurality of optical fiber, and the optical fiber temperature-measurement point is on the optical fiber for measuring the concrete The monitoring point of dam temperature;
The temperature value that the space coordinate that will acquire is measured with the optical fiber temperature-measurement point is set as mapping one by one;
According to the temperature value that the optical fiber temperature-measurement point measures, the two of optical fiber plane where calculating the optical fiber using finite difference calculus Tie up Temperature Distribution;
It is distributed according to the two-dimension temperature of two neighboring optical fiber plane, calculates the three-dimensional between the two neighboring optical fiber plane Thermo parameters method.
2. the temperature field forming method complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to claim 1, feature exist According to the whole coordinate of concrete dam and the buried line of optical fiber, the space coordinate for obtaining optical fiber temperature-measurement point includes:
Obtain the optical fiber inbuilt route in the concrete dam;
Take multiple monitoring points as the optical fiber temperature-measurement point in the route;
The space coordinate of the optical fiber temperature-measurement point is obtained according to the whole coordinate of the concrete dam.
3. the temperature field forming method complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to claim 1, feature exist In, according to the temperature value that the optical fiber temperature-measurement point measures, the two of optical fiber plane where calculating the optical fiber using finite difference calculus Tieing up Temperature Distribution includes:
Plane where the optical fiber is divided into multiple triangular units, wherein the node of the triangular unit is the optical fiber Temperature measuring point;
The two-dimension temperature in the triangular unit is calculated using the finite difference calculus to be distributed.
4. the temperature field forming method complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to claim 1, feature exist In being distributed according to the two-dimension temperature of two neighboring optical fiber plane, calculate the three-dimensional between the two neighboring optical fiber plane Thermo parameters method includes:
Obtain curve side element, the curve side element by the two neighboring optical fiber plane, the concrete dam upstream dam facing, Two faces that the downstream dam facing and the concrete dam of the concrete dam are contacted with air surround;
Utilize the distribution of three-dimensional temperature of curve side element described in three dimensional finite element method.
5. the temperature field forming method complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to claim 4, feature exist In, before the distribution of three-dimensional temperature using curve side element described in three dimensional finite element method, the method also includes:
The curve side element is converted to by regular hexahedron unit using coordinate transformation method.
6. a kind of temperature field complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement forms system characterized by comprising
Acquiring unit, for obtaining optical fiber temperature-measurement point according to the whole coordinate of concrete dam and the buried line of optical fiber Space coordinate, wherein the concrete dam is embedded with a plurality of optical fiber, and the optical fiber temperature-measurement point is to be used on the optical fiber Measure the monitoring point of the concrete dam temperature;
Setting unit, the temperature value that the space coordinate for will acquire is measured with the optical fiber temperature-measurement point are set as mapping one by one;
First computing unit, the temperature value for being measured according to the optical fiber temperature-measurement point calculate the light using finite difference calculus The two-dimension temperature distribution of optical fiber plane where fine;
Second computing unit calculates described two neighboring for being distributed according to the two-dimension temperature of two neighboring optical fiber plane Distribution of three-dimensional temperature between optical fiber plane.
7. the temperature field complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to claim 6 forms system, feature exists In the acquiring unit includes:
First obtains module, for obtaining the optical fiber inbuilt route in the concrete dam;
Selecting module, for selecting multiple monitoring points as the optical fiber temperature-measurement point in the route;
Second obtains module, and the space for obtaining the optical fiber temperature-measurement point for the whole coordinate according to the concrete dam is sat Mark.
8. the temperature field complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to claim 6 forms system, feature exists In first computing unit includes:
Division module, for the plane where the optical fiber to be divided into multiple triangular units, wherein the knot of the triangular unit Point is the optical fiber temperature-measurement point;
First computing module is distributed for calculating the two-dimension temperature in the triangular unit using the finite difference calculus.
9. the temperature field complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to claim 6 forms system, feature exists In first computing unit includes:
Third obtains module, and for obtaining curve side element, the curve side element is by the two neighboring optical fiber plane, the coagulation Two faces that the upstream dam facing of native dam, the downstream dam facing of the concrete dam and the concrete dam are contacted with air It surrounds;
Second computing module, for the distribution of three-dimensional temperature using curve side element described in three dimensional finite element method.
10. the temperature field complementary with Digital Measurement of Temperature based on optical fiber temperature-measurement according to claim 9 forms system, feature exists In, the system also includes:
Conversion module, for the curve side element to be converted to regular hexahedron unit using coordinate transformation method.
CN201910229044.4A 2019-03-25 2019-03-25 Temperature field forming method and system based on complementation of optical fiber temperature measurement and digital temperature measurement Active CN109974892B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910229044.4A CN109974892B (en) 2019-03-25 2019-03-25 Temperature field forming method and system based on complementation of optical fiber temperature measurement and digital temperature measurement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910229044.4A CN109974892B (en) 2019-03-25 2019-03-25 Temperature field forming method and system based on complementation of optical fiber temperature measurement and digital temperature measurement

Publications (2)

Publication Number Publication Date
CN109974892A true CN109974892A (en) 2019-07-05
CN109974892B CN109974892B (en) 2020-07-03

Family

ID=67080351

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910229044.4A Active CN109974892B (en) 2019-03-25 2019-03-25 Temperature field forming method and system based on complementation of optical fiber temperature measurement and digital temperature measurement

Country Status (1)

Country Link
CN (1) CN109974892B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110440950A (en) * 2019-07-29 2019-11-12 湖北省水利水电规划勘测设计院 Mass concrete temperature monitoring visualization system and method based on optical fiber temperature-measurement
CN111693175A (en) * 2020-07-22 2020-09-22 广州大学 Method, system, device and medium for measuring space temperature field based on fiber bragg grating string
CN113447159A (en) * 2021-08-30 2021-09-28 中国空气动力研究与发展中心设备设计与测试技术研究所 Wind tunnel cold leakage monitoring method based on global normal distribution and local outlier factor
CN113609198A (en) * 2021-07-26 2021-11-05 三峡大学 Method for reconstructing concrete temperature field in pouring bin based on optical fiber measured temperature data
CN115311805A (en) * 2022-08-18 2022-11-08 重庆大学 Fire-fighting safety system integrating multi-parameter sensing and signal processing method
CN117057239A (en) * 2023-08-14 2023-11-14 上海频准激光科技有限公司 Optical fiber temperature measuring point positioning system based on laser

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104298884A (en) * 2014-10-17 2015-01-21 武汉科技大学 Finite element and finite difference coupling method for fast calculating rolled piece section temperature
CN104809334A (en) * 2015-04-08 2015-07-29 河海大学 Calculating method of large-volume concrete cooling temperature field
CN106021720A (en) * 2016-05-19 2016-10-12 河海大学 Method for simulating large-volume concrete water pipe cooling temperature field
CN106500867A (en) * 2016-10-26 2017-03-15 上海市政工程设计研究总院(集团)有限公司 The installation method of dykes and dams monitoring internal temperature system
US20170199093A1 (en) * 2015-01-08 2017-07-13 Nec Laboratories America, Inc. Multi-parameter Sensing based on Few-mode Fiber Bragg Gratings using Femtosecond IR Laser
KR20170092067A (en) * 2016-02-02 2017-08-10 서울과학기술대학교 산학협력단 Finite element analysis method for optimally designing dam on wafer
CN108612058A (en) * 2018-04-08 2018-10-02 三峡大学 A kind of concrete dam distributed temperature measuring optical fiber bifilar " zigzag " leads to storehouse burial and construction method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104298884A (en) * 2014-10-17 2015-01-21 武汉科技大学 Finite element and finite difference coupling method for fast calculating rolled piece section temperature
US20170199093A1 (en) * 2015-01-08 2017-07-13 Nec Laboratories America, Inc. Multi-parameter Sensing based on Few-mode Fiber Bragg Gratings using Femtosecond IR Laser
CN104809334A (en) * 2015-04-08 2015-07-29 河海大学 Calculating method of large-volume concrete cooling temperature field
KR20170092067A (en) * 2016-02-02 2017-08-10 서울과학기술대학교 산학협력단 Finite element analysis method for optimally designing dam on wafer
CN106021720A (en) * 2016-05-19 2016-10-12 河海大学 Method for simulating large-volume concrete water pipe cooling temperature field
CN106500867A (en) * 2016-10-26 2017-03-15 上海市政工程设计研究总院(集团)有限公司 The installation method of dykes and dams monitoring internal temperature system
CN108612058A (en) * 2018-04-08 2018-10-02 三峡大学 A kind of concrete dam distributed temperature measuring optical fiber bifilar " zigzag " leads to storehouse burial and construction method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
丁勇 等: "空间结构瞬态温度场的Fourier-有限元分析", 《清华大学学报(自然科学版)》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110440950A (en) * 2019-07-29 2019-11-12 湖北省水利水电规划勘测设计院 Mass concrete temperature monitoring visualization system and method based on optical fiber temperature-measurement
CN111693175A (en) * 2020-07-22 2020-09-22 广州大学 Method, system, device and medium for measuring space temperature field based on fiber bragg grating string
CN113609198A (en) * 2021-07-26 2021-11-05 三峡大学 Method for reconstructing concrete temperature field in pouring bin based on optical fiber measured temperature data
CN113609198B (en) * 2021-07-26 2023-08-25 三峡大学 Method for reconstructing concrete temperature field in pouring warehouse based on optical fiber measured temperature data
CN113447159A (en) * 2021-08-30 2021-09-28 中国空气动力研究与发展中心设备设计与测试技术研究所 Wind tunnel cold leakage monitoring method based on global normal distribution and local outlier factor
CN113447159B (en) * 2021-08-30 2021-11-09 中国空气动力研究与发展中心设备设计与测试技术研究所 Wind tunnel cold leakage monitoring method based on global normal distribution and local outlier factor
CN115311805A (en) * 2022-08-18 2022-11-08 重庆大学 Fire-fighting safety system integrating multi-parameter sensing and signal processing method
CN115311805B (en) * 2022-08-18 2024-02-27 重庆大学 Fire safety system integrating multi-parameter sensing and signal processing method
CN117057239A (en) * 2023-08-14 2023-11-14 上海频准激光科技有限公司 Optical fiber temperature measuring point positioning system based on laser
CN117057239B (en) * 2023-08-14 2024-04-12 上海频准激光科技有限公司 Optical fiber temperature measuring point positioning system based on laser

Also Published As

Publication number Publication date
CN109974892B (en) 2020-07-03

Similar Documents

Publication Publication Date Title
CN109974892A (en) Based on optical fiber temperature-measurement the temperature field forming method and system complementary with Digital Measurement of Temperature
CN106021755B (en) Nuclear island raft foundation of nuclear power station mass concrete temperature stress simulating analysis
CN106339820B (en) Real-time dynamic multidimensional analysis method and system for hydraulic building based on virtual reality
Morgan et al. The computation of three-dimensional flows using unstructured grids
CN106021720B (en) A kind of analogy method of Massive Volume Concrete Temperature Field with Water Pipe Cooling
CN106294963A (en) Direct-buried cable carrying current calculation method
CN102819647B (en) A kind of heterogeneous material random microscopic structure finite element modeling method
Teng et al. Temperature and displacement monitoring to steel roof construction of Shenzhen Bay Stadium
JP2011514587A (en) Computer-implemented method and corresponding computer program product for modeling a nuclear reactor core
CN106570252A (en) Probability technology based thermal protection system design method
CN109101675A (en) A kind of simulation solid material thermal fracture method
CN104424380A (en) Mechanical strain gauge simulation
CN103245691A (en) Method for measuring thermal resistance of three-core cable packing layer
Doukas et al. Identification of a polymer growth process with an equilibrium multicritical collapse phase transition: The meeting point of swollen, collapsed, and crystalline polymers
CN109959467A (en) A kind of universe dam determines Wen Fangfa
CN107871045A (en) A kind of transmission building element based on BIM quantifies extracting method
CN109559377A (en) A method of it is cut using longitude and latitude polygon and generates three-dimensional map model
CN110440950A (en) Mass concrete temperature monitoring visualization system and method based on optical fiber temperature-measurement
CN111159936A (en) Cable joint thermal field calculation method based on generalized time domain finite difference
CN109598059A (en) A kind of thermal protection system optimum design method and designing system based on agent model
CN108009562A (en) A kind of hydrographic water resource feature space variability knows method for distinguishing
Ciancio et al. A method for the calculation of inter-element stresses in 3D
CN106156537A (en) The Force Feedback Model modeling method of radial direction base mesh free soft tissue data based on wheat quart algorithm
Si et al. Development of simulation program for temperature field of mass concrete structures
Bialecki et al. Reconstruction of time-dependent boundary heat flux by a BEM-based inverse algorithm

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