CN110442967A - Back Analysis of Concrete Thermal Parameters method - Google Patents
Back Analysis of Concrete Thermal Parameters method Download PDFInfo
- Publication number
- CN110442967A CN110442967A CN201910716033.9A CN201910716033A CN110442967A CN 110442967 A CN110442967 A CN 110442967A CN 201910716033 A CN201910716033 A CN 201910716033A CN 110442967 A CN110442967 A CN 110442967A
- Authority
- CN
- China
- Prior art keywords
- concrete
- parameters
- temperature rise
- days
- parameter
- 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.)
- Pending
Links
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
Back Analysis of Concrete Thermal Parameters method belongs to the technical field of the concrete structure monitoring and construction in hydraulic and hydroelectric engineering, and in particular to a kind of mass concrete thermal parameters inverse model utilizes the method for the model back analysis concrete thermal parameters.This method is made of parameters sensitivity analysis, improved response phase method and optimizing algorithm.The present invention provides a kind of mass concrete thermal parameters inverse model and utilize its inverse analysis method, by it is quick, accurately calculate two parameter values of final adiabatic temperature rise and heat of hydration half number of days, it will appreciate that concrete internal temperature field changing rule, to ensure construction quality.
Description
Technical field
The invention belongs to the technical fields of concrete structure monitoring and construction in hydraulic and hydroelectric engineering, and in particular to a kind of
Mass concrete thermal parameters inverse model utilizes the method for the model back analysis concrete thermal parameters.
Background technique
In current Practical Project, concrete is most widely used as a kind of most common construction material.Especially
Nowadays such as high concrete gravity dam, induced joint, Offshore Concrete Platforms, over strait big of the large volume concrete structural under complex environment
The appearance of bridge etc., so that people are increasingly urgent to the performance understanding of concrete especially mass concrete.Large volume coagulation
Soil is often relatively fewer in unit volume concrete cement dosage, and the heat of hydration is lower, but the requirement of its speed of perfusion is very high, in order to
Achieve the purpose that its continuous placing rapid construction, generally use concreting without longitudinal joint, pour is larger, radiates slower, generally requires several
Natural cooling in 10 years can be only achieved equilibrium temperature.However higher temperature field can cause biggish temperature stress, this is to substantially
Product concrete construction it is very unfavorable, therefore understand mass concrete inside temperature field and its changing rule it is quite important.
Presently described concrete internal temperature field changing rule generallys use the description of adiabatic temperature rise formula, adiabatic temperature rise formula
Mainly use exponential type and hyperbolic-type, both models require to know final adiabatic temperature rise and heat of hydration half number of days this two
A important parameter.The two parameters can be obtained by laboratory test, but in Practical Project, true execution conditions are very multiple
Miscellaneous, laboratory test is difficult to simulate true construction environment, and the parameter that laboratory test value obtains can not fit like a glove and construct
The concrete temperature field changing rule of completion.
Meanwhile on the other hand, with the development of sensor technology and computer technology, temperature of the people for inside concrete
Degree field and stress field monitoring become more convenient, provide new approach to understand the basic performance inside mass concrete.
Back analysis is a kind of effective ways that the practical thermal parameters of concrete are obtained by field monitoring data.Under normal conditions
The method of back analysis can be divided into two class of direct method and indirect method.Parameter inversion problems in engineering, will frequently with direct method
Material Parameters Inversion problem is converted to the minimum problem of Mathematical Planning, objective function is constructed according to inverting purpose, by linear
Law of planning, Direct Iterative Method or matrix inversion method seek its minimum, to obtain best parameter group.It is anti-in the above parameter
It drills in method, requires the solution for carrying out a large amount of forward modelling.But when inverted parameters are more, computational efficiency just becomes
It is time-consuming also longer to be low, therefore urgently need a kind of computational efficiency high, the accurate inverse analysis method of inverted parameters.
Summary of the invention
The object of the present invention is to provide a kind of mass concrete thermal parameters inverse analysis methods, and this method is by parameter sensitivity
Property analysis, improved response phase method and optimizing algorithm constitute.
Method of the invention uses ANSYS system, is a large-scale general finite element of ANSYS company, U.S. exploitation
Analyze software.Solid70 unit is the included unit of ANSYS program in software, and solid70 is for calculating temperature field.In ANSYS
It is middle to create three-dimensional finite element model with solid70 unit, calculate temperature field.
The Morris method that parameters sensitivity analysis uses is (see document " mining safety and environmental protection " 2018,45 (01): 102-
106, Liu Bo, Xiao Hongfei, Yang Ya just write " the deformation of the surrounding rock in tunnel overall situation sensitive analog analysis based on Morris method ") it is to pass through
It constructs several groups " orthogonal test ", " basic to influence " of the primary each parameter of value wheel stream calculation for only changing a parameter, thus
Can the influence to mode input parameter to output data assess, and obtain the comparison and parameter phase of parameter global sensitivity
Closing property and nonlinear qualitative description.Specific calculating process is:
If model output function y=f (x1,x2,x3…xn), containing n parameter, combine parameters to be obeyed first
The variation range of each parameter is mapped in closed interval [0,1] by probability distribution, and by the horizontal p of preset sampling by its from
Dispersion, each parameter can only be fromMiddle value, allows all random samplings in these sample points of n parameter
Once, and all only one parameters in two groups of adjacent sampling tests is allowed to change, and all occur Δ variable quantity, in this way according to
It is secondary to be sampled test, one group of sampling test is once completed until this n parameter changes in turn；It is calculated after the completion of test by following formula
The basic influence of each parameter:
In formula: EEjThe as basic influence of parameter x, then their mean μ and variances sigma are calculated, finally carry out sensibility point
Analysis.The value of μ is bigger, and influence of the corresponding parameter to computation model is more obvious, and the value of variances sigma is bigger, then illustrates that the parameter exists
Influence model output when and other parameters interaction it is bigger, that is, it is easier influenced by other parameters, in other words should
The influence that parameter exports model is nonlinear.
The optimizing algorithm uses genetic algorithm, encodes to the input parameter of response phase method, and in specific sections
Random seed is substituted into response phase method and is calculated, finds the individual for being best suitable for Response Face Function by interior generation random seed.
Back Analysis of Concrete Thermal Parameters method of the present invention, which is characterized in that method includes the following steps:
S1: by the feature of concrete dam, three-dimensional finite element model is created using solid70 unit in ANSYS；
S2: global sensitivity analysis is carried out to concrete thermal parameters using Morris method, is determined to inverted parameters and mould
Type output；
Choose final adiabatic temperature rise θ0It is used as with heat of hydration half number of days two parameters of n to inverted parameters, dam body inside is mixed
The maximum temperature for coagulating soil is exported as model；
S3: observed temperature data is extracted from the thermometer or other sensors for be embedded in inside concrete；
S4: solving Response Face Function,
Function is specifically:
In formula: N is final temperature rise number of parameters in reflection adiabatic temperature rise function；
M is the number of parameters of heat of hydration half number of days；
xiAnd yjFor the parameter to inverting, the final adiabatic temperature rise θ of selection is respectively corresponded0With the heat of hydration half number of days n two
Parameter；
S5: the actual measurement highest adiabatic temperature rise that will be obtained in step S3It substitutes into response surface equation, utilizes genetic algorithm
The nonlinear equation is solved, optimal solution is sought obtaining, finds the parameter xi and yj for being best suitable for response surface equation, that is, obtains final insulation temperature
Rise θ0With the optimal value of heat of hydration half number of days n.
The feature of the concrete dam includes concrete dam figure, concrete material subregion, and pours situation；Its
In, in the three-dimensional finite element model of creation, the not respective material parameter in same district is assigned according to dam body materials partitioning scenario and is used to have
Limit member calculates.
The Response Face Function solution is by randomly selecting 5 groups of final adiabatic temperature rise θ0With heat of hydration half number of days n
Two parameters substitute into the finite element model established in step S1, calculate separately the highest thermal insulation temperature rise under different parameters, and
Using the highest thermal insulation temperature rise asConstruct five independent systems of linear equations, simultaneous solution response surface coefficient a, bi、
ci、djAnd fj。
It, can the present invention provides a kind of mass concrete thermal parameters inverse model and using its inverse analysis method
Quickly and accurately inverting obtains the thermal parameters of concrete, and having based on true concrete structure can be used in calculating process
Limit meta-model is small using genetic algorithm Searching efficiency height, error without modeling again.It is calculated finally absolutely by quick, accurate
Two parameter values of hot temperature rise and heat of hydration half number of days, will appreciate that concrete internal temperature field changing rule, to ensure engineering
Quality.
Detailed description of the invention
Fig. 1 is a kind of concrete calorifics Parameter Inversion Model flow chart；
Fig. 2 is that Huang steps on power station Compacted Concrete Gravity Dam Section three-dimensional finite element model.
Specific embodiment
Understand for the ease of those of ordinary skill in the art and implement the present invention, with reference to the accompanying drawing to of the invention specific
Implementation is described in further detail.
Embodiment 1: the present invention proposes a kind of concrete thermal parameters so that Huang steps on power station Compacted Concrete Gravity Dam Section as an example
Inverse analysis method, comprising the following steps:
S1: by the feature of concrete dam, three-dimensional finite element model is created using solid70 unit in ANSYS；
The feature of concrete dam includes concrete dam figure, concrete material subregion, and pours situation；
Wherein, in the three-dimensional finite element model of creation, the not respective material in same district is assigned according to dam body materials partitioning scenario
Parameter is used for FEM calculation, the model built such as attached drawing 2.
S2: carrying out global sensitivity analysis to concrete thermal parameters using Morris method, is chosen by analysis final exhausted
Hot temperature rise θ0It is used as with heat of hydration half number of days two parameters of n to inverted parameters, the maximum temperature conduct of dam body inner concrete
Model output；
It treats inverted parameters and carries out sensitivity analysis, calculate maximum temperature under different input Parameter Conditions influences substantially
Mean μ and variances sigma；Calculated result mean μ is larger, and variances sigma is smaller, illustrates adiabatic temperature rise parameter to the output shadow of computation model
Sound is more obvious, and is affected by other parameters smaller, and stability is preferable during parametric inversion.
S3: Huang is extracted from the thermometer or other sensors for being embedded in inside concrete and steps on the dam concrete gravity dam 14#
Section observed temperature data.
S4: solving Response Face Function, which uses based on the mathematical model of Adiabatic temperature rise of concrete, comprehensively consider
The functional form of final adiabatic temperature rise and heat of hydration half number of days, constructs the response surface side of following form in Adiabatic temperature rise of concrete
Journey:
In formula: N is final temperature rise number of parameters in reflection adiabatic temperature rise function；
M is the number of parameters of heat of hydration half number of days；
xiAnd yjFor the parameter to inverting, the final adiabatic temperature rise θ of selection is respectively corresponded0With the heat of hydration half number of days n two
Parameter；
Randomly select 5 groups of final adiabatic temperature rise θ0With heat of hydration half number of days two parameters of n, substitute into S1 in establish it is limited
Meta-model, calculates separately the highest thermal insulation temperature rise under different parameters, and using the highest thermal insulation temperature rise asStructure
Make five independent systems of linear equations, simultaneous solution response surface coefficient a, bi、ci、djAnd fj；
Adiabatic temperature rise θ0Computer can be used to generate random number with two parameters of heat of hydration half number of days n, it can also be taking human as
At will choose；Usually based on test value or empirical value, choose in the reasonable scope, usually take laboratory experiment value 80%~
As selection numerical value between 120%.
S5: the actual measurement highest adiabatic temperature rise that will be obtained in S3It substitutes into response surface equation, is solved using genetic algorithm
The nonlinear equation seeks obtaining optimal solution, finds the parameter x for being best suitable for response surface equationiAnd yj, that is, obtain final adiabatic temperature rise θ0
With the optimal value of heat of hydration half number of days n.
Claims (3)
1. Back Analysis of Concrete Thermal Parameters method, which is characterized in that method includes the following steps:
S1: by the feature of concrete dam, three-dimensional finite element model is created using solid70 unit in ANSYS；
S2: global sensitivity analysis is carried out to concrete thermal parameters using Morris method, determination is defeated to inverted parameters and model
Out；
Choose final adiabatic temperature rise θ0It is used as with heat of hydration half number of days two parameters of n to inverted parameters, dam body inner concrete
Maximum temperature is exported as model；
S3: observed temperature data is extracted from the thermometer or other sensors for be embedded in inside concrete；
S4: solving Response Face Function,
Function is specifically:
In formula: N is final temperature rise number of parameters in reflection adiabatic temperature rise function；
M is the number of parameters of heat of hydration half number of days；
xiAnd yjFor the parameter to inverting, the final adiabatic temperature rise θ of selection is respectively corresponded0Join with the heat of hydration half number of days n two
Number；
S5: the actual measurement highest adiabatic temperature rise that will be obtained in step S3It substitutes into response surface equation, is solved using genetic algorithm
The nonlinear equation seeks obtaining optimal solution, finds the parameter xi and yj for being best suitable for response surface equation, that is, obtains final adiabatic temperature rise θ0
With the optimal value of heat of hydration half number of days n.
2. Back Analysis of Concrete Thermal Parameters method as described in claim 1, which is characterized in that the spy of the concrete dam
Sign includes concrete dam figure, concrete material subregion, and pours situation；Wherein, it in the three-dimensional finite element model of creation, presses
The not respective material parameter in same district, which is assigned, according to dam body materials partitioning scenario is used for FEM calculation.
3. Back Analysis of Concrete Thermal Parameters method as described in claim 1, which is characterized in that the Response Face Function is asked
Solution is by randomly selecting 5 groups of final adiabatic temperature rise θ0With heat of hydration half number of days two parameters of n, substitute into S1 in establish it is limited
Meta-model, calculates separately the highest thermal insulation temperature rise under different parameters, and using the highest thermal insulation temperature rise asStructure
Make five independent systems of linear equations, simultaneous solution response surface coefficient a, bi、ci、djAnd fj。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910716033.9A CN110442967A (en) | 2019-08-05 | 2019-08-05 | Back Analysis of Concrete Thermal Parameters method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910716033.9A CN110442967A (en) | 2019-08-05 | 2019-08-05 | Back Analysis of Concrete Thermal Parameters method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110442967A true CN110442967A (en) | 2019-11-12 |
Family
ID=68433118
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910716033.9A Pending CN110442967A (en) | 2019-08-05 | 2019-08-05 | Back Analysis of Concrete Thermal Parameters method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110442967A (en) |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62103557A (en) * | 1985-07-11 | 1987-05-14 | Sumitomo Cement Co Ltd | Adiabatic temperature rise tester for concrete and mortar |
| JP2008241520A (en) * | 2007-03-28 | 2008-10-09 | Tokyo Institute Of Technology | Adiabatic calorimeter and quality control method of cement and concrete using it |
| CN102721480A (en) * | 2012-06-27 | 2012-10-10 | 清华大学 | Method for calculating equivalent temperature field of large-size concrete based on cooling water monitoring |
| CN102979307A (en) * | 2012-12-12 | 2013-03-20 | 新疆生产建设兵团金来建设工程技术研发有限责任公司 | Temperature-controlled crack prevention construction method for concrete structure |
| CN103605841A (en) * | 2013-11-07 | 2014-02-26 | 河海大学 | Method for building numerical simulation model of hydraulic engineering asphalt concrete |
| CN103726434A (en) * | 2014-01-07 | 2014-04-16 | 东南大学 | Method for establishing temperature field model of steel box beam bridge road system under high-temperature asphalt concrete paving |
| CN103837533A (en) * | 2014-01-16 | 2014-06-04 | 河海大学 | Method for concrete temperature monitoring and simulation back analysis based on thermal imager |
| CN103914594A (en) * | 2014-03-26 | 2014-07-09 | 河海大学 | Concrete thermodynamic parameter intelligent recognition method based on support vector machine |
| CN105787174A (en) * | 2016-02-25 | 2016-07-20 | 武汉大学 | High-rockfill-dam transient-rheological-parameter inversion method based on response surface method |
| CN107301282A (en) * | 2017-06-12 | 2017-10-27 | 天津大学 | The concrete dam mechanics parameter inversion method of time series data is monitored based on multi-source |
| US20170370050A1 (en) * | 2016-06-23 | 2017-12-28 | University Of Florida Research Foundation, Inc. | Waste to energy ash and engineered aggregate in road construction |
| CN108133111A (en) * | 2017-12-29 | 2018-06-08 | 中铁十二局集团有限公司 | A kind of Study on Temperature Field method and temprature control method based on mass concrete |
| CN108981966A (en) * | 2018-07-02 | 2018-12-11 | 雷元新 | A kind of mass concrete temperature gradient Analysis of Limit Value method and device |
| CN109800473A (en) * | 2018-12-26 | 2019-05-24 | 武汉大学 | Rock mechanical parameters inversion method based on differential evolution method |
-
2019
- 2019-08-05 CN CN201910716033.9A patent/CN110442967A/en active Pending
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62103557A (en) * | 1985-07-11 | 1987-05-14 | Sumitomo Cement Co Ltd | Adiabatic temperature rise tester for concrete and mortar |
| JP2008241520A (en) * | 2007-03-28 | 2008-10-09 | Tokyo Institute Of Technology | Adiabatic calorimeter and quality control method of cement and concrete using it |
| CN102721480A (en) * | 2012-06-27 | 2012-10-10 | 清华大学 | Method for calculating equivalent temperature field of large-size concrete based on cooling water monitoring |
| CN102979307A (en) * | 2012-12-12 | 2013-03-20 | 新疆生产建设兵团金来建设工程技术研发有限责任公司 | Temperature-controlled crack prevention construction method for concrete structure |
| CN103605841A (en) * | 2013-11-07 | 2014-02-26 | 河海大学 | Method for building numerical simulation model of hydraulic engineering asphalt concrete |
| CN103726434A (en) * | 2014-01-07 | 2014-04-16 | 东南大学 | Method for establishing temperature field model of steel box beam bridge road system under high-temperature asphalt concrete paving |
| CN103837533A (en) * | 2014-01-16 | 2014-06-04 | 河海大学 | Method for concrete temperature monitoring and simulation back analysis based on thermal imager |
| CN103914594A (en) * | 2014-03-26 | 2014-07-09 | 河海大学 | Concrete thermodynamic parameter intelligent recognition method based on support vector machine |
| CN105787174A (en) * | 2016-02-25 | 2016-07-20 | 武汉大学 | High-rockfill-dam transient-rheological-parameter inversion method based on response surface method |
| US20170370050A1 (en) * | 2016-06-23 | 2017-12-28 | University Of Florida Research Foundation, Inc. | Waste to energy ash and engineered aggregate in road construction |
| CN107301282A (en) * | 2017-06-12 | 2017-10-27 | 天津大学 | The concrete dam mechanics parameter inversion method of time series data is monitored based on multi-source |
| CN108133111A (en) * | 2017-12-29 | 2018-06-08 | 中铁十二局集团有限公司 | A kind of Study on Temperature Field method and temprature control method based on mass concrete |
| CN108981966A (en) * | 2018-07-02 | 2018-12-11 | 雷元新 | A kind of mass concrete temperature gradient Analysis of Limit Value method and device |
| CN109800473A (en) * | 2018-12-26 | 2019-05-24 | 武汉大学 | Rock mechanical parameters inversion method based on differential evolution method |
Non-Patent Citations (1)
| Title |
|---|
| 徐江宇: "RCC重力坝固结灌浆期间仓面裂缝成因探讨", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技II辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | Numerical simulation of a flow-like landslide using the particle finite element method | |
| Bauer et al. | Transient 3D analysis of borehole heat exchanger modeling | |
| Khoei et al. | A mesh-independent finite element formulation for modeling crack growth in saturated porous media based on an enriched-FEM technique | |
| Raymond et al. | Numerical analysis of thermal response tests with a groundwater flow and heat transfer model | |
| CN102979307B (en) | A kind of Temperature-controllcrack crack prevention construction method for concrete structure | |
| CN109974892B (en) | Temperature field forming method and system based on complementation of optical fiber temperature measurement and digital temperature measurement | |
| CN109408971A (en) | A kind of early concrete condition Changeement method under novel effect of multiple parameters | |
| Bazazzadeh et al. | Simulation of chemo-thermo-mechanical problems in cement-based materials with Peridynamics | |
| CN111210519A (en) | High core rock-fill dam deformation parameter inversion method | |
| Abeka et al. | Thermal effect of mass concrete structures in the tropics: Experimental, modelling and parametric studies | |
| Jin et al. | Characterisation of temperature-dependent heat conduction in heterogeneous concrete | |
| Yu et al. | A geomechanics-coupled embedded discrete fracture model and its application in geothermal reservoir simulation | |
| Qiang et al. | A p-version embedded model for simulation of concrete temperature fields with cooling pipes | |
| Wang et al. | Novel interval theory‐based parameter identification method for engineering heat transfer systems with epistemic uncertainty | |
| Bi et al. | An analytical model for the thermal conductivity of soils during a freezing process | |
| Yu et al. | Inverse analysis of thermal conductivities in transient non-homogeneous and non-linear heat conductions using BEM based on complex variable differentiation method | |
| CN115828704A (en) | Rapid prediction method for underground water pollution | |
| Xiong et al. | Modification of thermo-elasto-viscoplastic model for soft rock and its application to THM analysis of heating tests | |
| CN108088870A (en) | The method for building up of Nonlinear Heat Conduction model test similarity criterion | |
| Yu et al. | An improved form of SPH method for simulating the thermo-mechanical-damage coupling problems and its applications | |
| CN100377152C (en) | Method for confirming stress intensity factor distribution on member crack tip | |
| Shi et al. | Thermal stress analyses and reinforcement design of massive RC structures | |
| Carka et al. | The Dirichlet-to-Neumann map for two-dimensional crack problems | |
| Mohammadiun et al. | Real-time evaluation of severe heat load over moving interface of decomposing composites | |
| Gadi et al. | Modeling soil-crack–water–atmospheric interactions: a novel root water uptake approach to simulate the evaporation through cracked soil and experimental validation |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191112 |