CN111209692A - Structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient - Google Patents

Structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient Download PDF

Info

Publication number
CN111209692A
CN111209692A CN201911258217.1A CN201911258217A CN111209692A CN 111209692 A CN111209692 A CN 111209692A CN 201911258217 A CN201911258217 A CN 201911258217A CN 111209692 A CN111209692 A CN 111209692A
Authority
CN
China
Prior art keywords
heat transfer
convective heat
temperature field
value
transfer coefficient
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
Application number
CN201911258217.1A
Other languages
Chinese (zh)
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.)
Sinomach Science And Technology Institute Co ltd
Hefei University of Technology
Hefei General Machinery Research Institute Co Ltd
Original Assignee
Sinomach Science And Technology Institute Co ltd
Hefei University of Technology
Hefei General Machinery Research Institute Co Ltd
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 Sinomach Science And Technology Institute Co ltd, Hefei University of Technology, Hefei General Machinery Research Institute Co Ltd filed Critical Sinomach Science And Technology Institute Co ltd
Priority to CN201911258217.1A priority Critical patent/CN111209692A/en
Publication of CN111209692A publication Critical patent/CN111209692A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a structural temperature field numerical calculation method based on equivalent convective heat transfer coefficients. The method comprises the following steps: s1, assigning different equivalent convection heat transfer coefficients, and performing temperature field simulation on an in-service research object by using ABAQUS finite element software to obtain a simulated temperature value of an analysis point; s2, calculating and obtaining error values of the measured temperature value and the simulated temperature value of the analysis point under different equivalent convective heat transfer coefficients according to the following formula:
Figure DDA0002310902150000011
s3, determining the error value
Figure DDA0002310902150000012
Minimum value of
Figure DDA0002310902150000013
The minimum value
Figure DDA0002310902150000014
The equivalent convective heat transfer coefficient of the corresponding assignment is optimal, and the obtained simulated temperature field is obtained. The simulated temperature field obtained by the calculation method is closer to the actual temperature field of the in-service research object, so that the reliability of the subsequent structural response analysis and the residual life prediction result of the in-service research object is improved, and the method has an important value in improving the safety of the structure in operation.

Description

Structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient
Technical Field
The invention belongs to the technical field of heat transfer, and particularly relates to a structural temperature field numerical calculation method based on equivalent convective heat transfer coefficients.
Background
The heat convection problem generally exists in the service process of industrial equipment such as aerospace, thermal machinery, energy, chemical engineering and the like, and the thermal stress caused by the large-scale change of temperature is often the main factor causing the structural failure. Therefore, the research on the convection heat exchange problem when the structure is in service has important value for ensuring the long-period safe operation of the structure. In order to obtain the structural three-dimensional space temperature field, only a test method is adopted, the workload is huge, and a large amount of manpower and material resources are consumed. In recent thirty years, with the rapid development of computer technology, numerical simulation methods have been widely used for analyzing convective heat transfer processes. In order to clarify the temperature distribution rule of the structure in the service process, a method combining tests and simulation is an effective way.
The convective heat transfer coefficient is a key parameter in convective heat transfer analysis, and plays a crucial role in the distribution of the structural three-dimensional space temperature field. Due to the complexity of the flow field when the structure is in service, the equivalent convection heat transfer coefficient is considered more reasonably. At present, scholars at home and abroad often take values by experience, so that the deviation of the temperature field numerical analysis result from the actual situation is large, the reliability of the subsequent structural response analysis and residual life prediction result is low, and finally, the safe operation in the service period of the structure is adversely affected.
Disclosure of Invention
In order to solve the technical problem, the invention provides a structural temperature field numerical calculation method based on equivalent convective heat transfer coefficients.
In order to realize the purpose of the invention, the invention adopts the following technical scheme:
a structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient comprises the following steps:
s1, assigning different equivalent convection heat transfer coefficients, and performing temperature field simulation on an in-service research object by using ABAQUS finite element software to obtain a simulated temperature value of an analysis point;
s2, calculating and obtaining error values of the measured temperature value and the simulated temperature value of the analysis point under different equivalent convective heat transfer coefficients according to the following formula:
Figure BDA0002310902130000021
in the formula:
Figure BDA0002310902130000022
the error value of the measured temperature value and the simulated temperature value is related to the equivalent convective heat transfer coefficient;
n represents the number of analysis points, i represents the ith analysis point;
m represents the number of measured temperature values of a single analysis point, and j represents the jth measured temperature value;
Figure BDA0002310902130000023
representing the equivalent convective heat transfer coefficient;
Figure BDA0002310902130000024
representing a simulated temperature value related to time and equivalent convective heat transfer coefficient;
Tij a(t) represents a measured temperature value, time dependent;
t represents time;
s3, determining the error value
Figure BDA0002310902130000025
Minimum value of
Figure BDA0002310902130000026
The minimum value
Figure BDA0002310902130000027
The equivalent convective heat transfer coefficient of the corresponding assignment is optimal, and the obtained simulated temperature field is obtained.
The further technical scheme is as follows: the temperature field simulation of the in-service research object by using ABAQUS finite element software in S1 comprises the following steps:
s1-1, establishing a numerical analysis model of an in-service research object;
s1-2, defining the technological parameters of the in-service research object in service and defining the thermophysical property data of the material used by the in-service research object;
s1-3, setting thermal boundary conditions during finite element analysis;
and S1-4, assigning an equivalent convection heat transfer coefficient, and carrying out numerical calculation to obtain a simulated temperature field.
The further technical scheme is as follows: the number of analysis points is at least 3.
The further technical scheme is as follows: and S1-4, the numerical calculation adopts transient analysis calculation, and the time step of the numerical calculation corresponds to the actual sampling time.
The further technical scheme is as follows: and when the thermal boundary condition is a dynamic thermal boundary condition, converting the interface rising speed of the time dimension into the axial height change of the in-service research object in the numerical calculation process, thereby obtaining the simulated temperature field of the in-service research object at any moment.
The invention has the beneficial effects that:
(1) the simulated temperature field obtained by the calculation method is closer to the actual temperature field of the in-service research object, so that the reliability of the subsequent structural response analysis and the residual life prediction result of the in-service research object is improved, and the method has an important value in improving the safety of the structure in operation. The method has wide applicability and is suitable for predicting or simulating the temperature field of the in-service research object in a high-temperature or low-temperature environment.
(2) The on-site temperature measuring points of the in-service research object are selected from 3 or more than 3, so that the optimized equivalent convective heat transfer coefficient is more accurate.
(3) From whether the temperature field is a function of time, thermal analysis includes steady-state and transient thermal analysis. The invention adopts transient analysis calculation related to time, and can better reveal the change rule of the temperature of the research object along with the time when the research object is actually in service.
Drawings
FIG. 1 is a schematic flow chart of the present invention.
FIG. 2 is a schematic diagram of a coke drum configuration.
FIG. 3 is a graph of measured temperatures at various points on the outer wall of the main process stage during a typical cycle of a coke drum.
FIG. 4 is a two-dimensional axisymmetric numerical analysis model of a coke drum.
FIG. 5 is a comparison graph of the measured temperature and the simulated temperature at different points.
FIG. 6 is a graph of simulated temperature profiles for various locations on the outer wall of a coke drum at the main process stage.
Detailed Description
The technical scheme of the invention is more specifically explained by combining the following embodiments:
the embodiment adopted by the invention is an in-service coke tower of a certain petrochemical plant in China, the main body part of the coke tower is formed by welding 9 cylinder sections, the top of the coke tower is a semi-elliptical head, the bottom of the coke tower is a conical head, the whole tower body is supported by a skirt, and the structural schematic diagram is shown in figure 2. The dimensional specification of the container is as follows: phi 9.8m multiplied by 36.6m, the wall thickness of the oval head and Course1-3 is 33mm, the wall thickness of Course 4-5 is 34mm, the thickness of Course 6 is 36mm, the thickness of Course 7-8 is 40mm, the thickness of Course 9 and the conical head is 44mm, and the thickness of the skirt is 38mm, wherein only the oval head and the upper segment Course1-3 contain a lining of 3 mm. Table 1 shows the operating parameters of the coke drum in service. The main body of the coke tower is made of SA387Gr11Cl2 and 410S, and the material performance parameters at different temperatures are obtained by referring to ASME boiler and pressure vessel specifications, and the specific values are shown in tables 2 and 3.
The coke tower is large-scale petrochemical equipment operated periodically and intermittently, the working period of the coke tower is 36h in the embodiment, 4 temperature measuring points are arranged on the outer wall of the coke tower, the specific positions are shown in figure 2, and the temperature sampling frequency is 3 min/time. The measured temperature profile at various points on the outer wall of the main process stage during a typical cycle of a coke drum is shown in figure 3.
According to the actual structural parameters of the coke tower, establishing a two-dimensional axisymmetric numerical analysis model shown in figure 4 by using ABAQUS finite element software, setting the technological parameters of an in-service research object in service, the thermophysical property data of the material used by the in-service research object, thermal boundary conditions and equivalent convection heat transfer coefficients in the ABAQUS finite element software, then performing numerical calculation by using transient analysis in the ABAQUS finite element software to obtain a simulated temperature field, wherein the numerical calculation time step corresponds to the actual sampling time.
Considering that the outer wall of the coke tower is provided with an insulating layer, an insulating boundary is adopted during numerical simulation, the inner wall of the coke tower is based on a third type of boundary conditions, the third type of boundary conditions comprise a static thermal boundary, a dynamic thermal boundary and a mixed thermal boundary, wherein the static thermal boundary is a thermal boundary height which does not change along with time during thermal analysis; the dynamic thermal boundary is the change of the thermal boundary height along with the time during thermal analysis; a hybrid thermal boundary is one in which a portion of the thermal boundary is time-varying and a portion of the thermal boundary is time-invariant.
For example, in the oil inlet or water inlet stage of the coke tower, because liquid below the liquid level exchanges heat with the inner wall of the tower, and the inner wall of the tower above the liquid level exchanges heat with gas in the tower, the thermal boundary formed between the liquid level and the tower wall is constantly changed along the axial direction of the tower as the liquid level rises at a constant speed, and in order to simulate the dynamic thermal boundary, the ABAQUS finite element software adopts an iterative algorithm to convert the rising speed problem of the thermal boundary interface of the time dimension into the change of the height along the axial direction of the tower for processing, so as to calculate the transient temperature field of the coke tower at any time t, the basic idea is as follows:
dividing the time interval 0-t into n identical time intervals delta t of infinite size, the starting time of each time interval being tm-1Indicating that the end time is tmIndicates, start time t 00, so tm=tm-1+ Δ t, where m ═ 1,2, 3. When the time interval Δ t → 0, the thermal boundary of the transient temperature field is assumed to be static during Δ t, i.e. the liquid and the inner wall of the column are considered to complete the heat transfer process in 2 steps during Δ t: the liquid level is raised (but the time required for the liquid level to rise is ignored), then the liquid level is kept constant, and then heat conduction is carried out. The heat transfer process above the liquid level is treated similarly.
The invention adopts the process parameters in the table 1 and continuously optimizes the equivalent convective heat transfer coefficient through the formula (1) to enable the simulated temperature field to approach the actually measured temperature field more truly, thereby obtaining the minimum value of the function E (h). The comparison curve of the temperature simulation value and the measured value of different measuring points obtained after the optimization design is shown in fig. 5. As can be seen from FIG. 5, the measured temperature distribution of each measuring point in the preheating stage is well matched with the calculation result, and although a certain error exists between the measured temperature value and the simulated temperature value of each measuring point in the oil-feeding coke-forming stage and the cooling stage, the variation trend is better in conformity. Therefore, the coke tower space three-dimensional temperature field obtained by numerical calculation of the method is closer to reality, the data reliability is high, and the heat exchange change rule of the coke tower in actual service can be truly reflected. Table 4 shows the equivalent convective heat transfer coefficients of the main process stages obtained after optimization.
FIG. 6 is a graph showing simulated temperature profiles for different locations on the outer wall of a coke drum at the main process stage of a single cycle. From the shape of the curve, the temperature changes with time in different positions have similar trends, and the temperature rises → the heat balance → the temperature falls. As can be seen from fig. 6, there is a sudden inflection point in the oil and water intake stages when the liquid level just reaches this position, because the medium temperature and equivalent convective heat transfer coefficient at this stage of the process are different from those at the previous stage. The inflection point moves to the right continuously as the liquid level rises.
TABLE 1 Coke drum production Process Main Process parameters
Stage of the process Medium Parameters of the process medium Duration of time
Steam preheating Low pressure steam 220℃ 1h
Oil gas preheating Oil gas 420℃ 3h
Oil-feeding coke Residual oil 488℃ 18h
Steam cooling Low pressure steam plus petroleum coke 220℃ 3h
Feed water cold coke Water plus Petroleum coke 60℃ 5.5h
TABLE 2 SA387Gr11Cl2 Material Performance parameters (thermophysical data)
Figure BDA0002310902130000051
Figure BDA0002310902130000061
TABLE 3410S materials Performance parameters (thermophysical data)
Figure BDA0002310902130000062
TABLE 4 main process stage equivalent convective heat transfer coefficient obtained after optimization
Figure BDA0002310902130000063

Claims (4)

1. A structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient is characterized by comprising the following steps:
s1, assigning different equivalent convection heat transfer coefficients, and performing temperature field simulation on an in-service research object by using ABAQUS finite element software to obtain a simulated temperature value of an analysis point;
s2, calculating and obtaining error values of the measured temperature value and the simulated temperature value of the analysis point under different equivalent convective heat transfer coefficients according to the following formula:
Figure FDA0002310902120000011
in the formula:
Figure FDA0002310902120000012
the error value of the measured temperature value and the simulated temperature value is related to the equivalent convective heat transfer coefficient;
n represents the number of analysis points, i represents the ith analysis point;
m represents the number of measured temperature values of a single analysis point, and j represents the jth measured temperature value;
Figure FDA0002310902120000013
representing the equivalent convective heat transfer coefficient;
Figure FDA0002310902120000014
representing a simulated temperature value related to time and equivalent convective heat transfer coefficient;
Tij a(t) represents a measured temperature value, time dependent;
t represents time;
s3, determining the error value
Figure FDA0002310902120000015
Minimum value of
Figure FDA0002310902120000016
The minimum value
Figure FDA0002310902120000017
The equivalent convective heat transfer coefficient of the corresponding assignment is optimal, and the obtained simulated temperature field is obtained.
2. The method for calculating the structural temperature field value based on the equivalent convective heat transfer coefficient as claimed in claim 1, wherein the step of simulating the temperature field of the in-service study object by using ABAQUS finite element software in S1 comprises the following steps:
s1-1, establishing a numerical analysis model of an in-service research object;
s1-2, defining the technological parameters of the in-service research object in service and defining the thermophysical property data of the material used by the in-service research object;
s1-3, setting thermal boundary conditions during finite element analysis;
and S1-4, assigning an equivalent convection heat transfer coefficient, and carrying out numerical calculation to obtain a simulated temperature field.
3. The method for calculating a structural temperature field value based on equivalent convective heat transfer coefficient of claim 1 wherein the number of analysis points is at least 3.
4. The structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient as recited in claim 2, characterized in that the numerical calculation in S1-4 adopts transient analysis calculation, and the time step of the numerical calculation corresponds to the actual sampling time.
CN201911258217.1A 2019-12-10 2019-12-10 Structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient Pending CN111209692A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911258217.1A CN111209692A (en) 2019-12-10 2019-12-10 Structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911258217.1A CN111209692A (en) 2019-12-10 2019-12-10 Structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient

Publications (1)

Publication Number Publication Date
CN111209692A true CN111209692A (en) 2020-05-29

Family

ID=70788086

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911258217.1A Pending CN111209692A (en) 2019-12-10 2019-12-10 Structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient

Country Status (1)

Country Link
CN (1) CN111209692A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112347669A (en) * 2020-10-09 2021-02-09 中国科学院国家天文台 Large-scale antenna back frame temperature measurement and real-time evaluation system and method
CN117272762A (en) * 2023-11-21 2023-12-22 中国科学院合肥物质科学研究院 Method and system for determining convective heat transfer coefficient of water-cooled magnet coil

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013091487A1 (en) * 2011-12-19 2013-06-27 上海望特能源科技有限公司 Method for monitoring inside-boiler dynamic wall temperature of power plant boiler high-temperature piping system
CN104063552A (en) * 2014-07-04 2014-09-24 上汽通用五菱汽车股份有限公司 Engine exhaust manifold thermal stress analysis and structural optimization method
CN106768909A (en) * 2016-11-15 2017-05-31 沈阳建筑大学 A kind of electro spindle temperature field prediction experimental system and method
CN107169157A (en) * 2017-04-07 2017-09-15 上海电气集团股份有限公司 A kind of structural thermal analysis finite element modeling method
CN108372208A (en) * 2018-04-20 2018-08-07 九江学院 A method of calculating continuous casting billet high-pressure water descaling convection transfer rate
CN109271688A (en) * 2018-09-03 2019-01-25 武汉理工大学 The Method for Accurate Calculation in high gear engaged transmission temperature field
CN109657368A (en) * 2018-12-24 2019-04-19 奥克斯空调股份有限公司 A kind of optimization method of air conditioner heat radiator
CN109871642A (en) * 2019-03-07 2019-06-11 上海电气集团股份有限公司 A kind of temperature field modeling method of the internal structure for steam generator

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013091487A1 (en) * 2011-12-19 2013-06-27 上海望特能源科技有限公司 Method for monitoring inside-boiler dynamic wall temperature of power plant boiler high-temperature piping system
CN104063552A (en) * 2014-07-04 2014-09-24 上汽通用五菱汽车股份有限公司 Engine exhaust manifold thermal stress analysis and structural optimization method
CN106768909A (en) * 2016-11-15 2017-05-31 沈阳建筑大学 A kind of electro spindle temperature field prediction experimental system and method
CN107169157A (en) * 2017-04-07 2017-09-15 上海电气集团股份有限公司 A kind of structural thermal analysis finite element modeling method
CN108372208A (en) * 2018-04-20 2018-08-07 九江学院 A method of calculating continuous casting billet high-pressure water descaling convection transfer rate
CN109271688A (en) * 2018-09-03 2019-01-25 武汉理工大学 The Method for Accurate Calculation in high gear engaged transmission temperature field
CN109657368A (en) * 2018-12-24 2019-04-19 奥克斯空调股份有限公司 A kind of optimization method of air conditioner heat radiator
CN109871642A (en) * 2019-03-07 2019-06-11 上海电气集团股份有限公司 A kind of temperature field modeling method of the internal structure for steam generator

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
侯文富;鲁志兵;张绍良;梁文彬;张韶伟: "操作参数对Cr-Mo钢焦炭塔瞬态温度场及热应力影响数值模拟研究", 压力容器, vol. 35, no. 1, pages 38 - 43 *
朱金花;董杰;杨小林;鲁志兵: "铬钼钢焦炭塔在线应变和温度测试", 化工机械, vol. 45, no. 4 *
王文博,赵昌昌: "统计学 经济社会统计", vol. 2005, 30 June 2005, 西安交通大学出版社, pages: 294 - 296 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112347669A (en) * 2020-10-09 2021-02-09 中国科学院国家天文台 Large-scale antenna back frame temperature measurement and real-time evaluation system and method
CN112347669B (en) * 2020-10-09 2024-04-12 中国科学院国家天文台 Large-sized antenna back frame temperature measurement and real-time evaluation system and method
CN117272762A (en) * 2023-11-21 2023-12-22 中国科学院合肥物质科学研究院 Method and system for determining convective heat transfer coefficient of water-cooled magnet coil
CN117272762B (en) * 2023-11-21 2024-02-06 中国科学院合肥物质科学研究院 Method and system for determining convective heat transfer coefficient of water-cooled magnet coil

Similar Documents

Publication Publication Date Title
CN108062427B (en) Method for reducing forging residual stress of turbine disc based on numerical calculation gradient speed control
CN111651908B (en) ANSYS-based multi-field coupling stress distribution simulation method for high-temperature carbonization furnace
CN110633520B (en) Evaluation method for flow impact of sodium-cooled fast reactor core outlet
CN108875213B (en) Reactor core thermal hydraulic multi-scale analysis method
CN111209692A (en) Structural temperature field numerical calculation method based on equivalent convective heat transfer coefficient
CN111680432B (en) Low-temperature carbonization furnace multi-coupling field stress distribution simulation method based on WORKBENCH
CN111400934B (en) Method for simulating oxygen content distribution of low-temperature carbonization furnace based on WORKBENCH
CN113569504A (en) Prediction method and prediction system for creep fatigue life of combustion chamber of aircraft engine
CN111879910A (en) Test method for optimizing forging process parameters and structure performance
CN116306301A (en) Storage tank thermal buckling failure prediction method based on neural network
CN115310311B (en) Method for calculating fatigue life of metal material under spectral load
CN112507567B (en) Method for predicting instability defect of forged microstructure of titanium alloy forging
Han et al. Numerical analysis of flow fields and temperature fields in a regenerative heating furnace for steel pipes
WO2022011724A1 (en) Workbench-based method for simulating multi-coupling field stress distribution of low-temperature carbonization furnace
Bonetti et al. A lifing method based on R5 and viscoplasticity modelling and its application to turbine rotor under thermomechanical loading
CN109959555B (en) Method for predicting induction period of central crack tensile sample by considering thickness and crack depth influence
Deng et al. Optimization of the location of the oil cooling gallery in the diesel engine piston
Yin et al. Numerical Simulation of Thermal-Mechanical Stress Field of Coke Drum Considering Deformation
CN116822299B (en) Rapid calculation method for thermal stress of aeroengine flame tube under service load course
CN113673057B (en) Design analysis method for cooling water jacket of diesel engine based on CAE simulation technology
CN113011068B (en) Three-dimensional simulation method for walking beam type plate blank heating
CN111928965B (en) Rotor temperature monitoring method and device, computer equipment and storage medium
CN115455769A (en) Rotating equipment thermal deformation finite element solving method based on measured temperature information
CN111310391B (en) Simulation method of plate-fin heat exchanger
CN109255136B (en) Method for predicting incubation period of crack containing elliptical axial inner surface in high-temperature pipeline

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