Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP5994086B2 - Deformation calculation method for brick structure - Google Patents
[go: Go Back, main page]

JP5994086B2 - Deformation calculation method for brick structure - Google Patents

Deformation calculation method for brick structure Download PDF

Info

Publication number
JP5994086B2
JP5994086B2 JP2013082919A JP2013082919A JP5994086B2 JP 5994086 B2 JP5994086 B2 JP 5994086B2 JP 2013082919 A JP2013082919 A JP 2013082919A JP 2013082919 A JP2013082919 A JP 2013082919A JP 5994086 B2 JP5994086 B2 JP 5994086B2
Authority
JP
Japan
Prior art keywords
brick
brick structure
coke oven
bricks
deformation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2013082919A
Other languages
Japanese (ja)
Other versions
JP2014205741A (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.)
JFE Steel Corp
Nagoya City University
Original Assignee
JFE Steel Corp
Nagoya City University
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 JFE Steel Corp, Nagoya City University filed Critical JFE Steel Corp
Priority to JP2013082919A priority Critical patent/JP5994086B2/en
Publication of JP2014205741A publication Critical patent/JP2014205741A/en
Application granted granted Critical
Publication of JP5994086B2 publication Critical patent/JP5994086B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Furnace Housings, Linings, Walls, And Ceilings (AREA)

Description

本発明は、コークス炉などのようなレンガ構造物の変形量を構造解析により算出する方法に関するものである。   The present invention relates to a method for calculating a deformation amount of a brick structure such as a coke oven by structural analysis.

レンガ構造物の一例として、コークス炉(室炉式コークス炉)が挙げられる。図1はコークス炉の概略斜視図、図2は炉体内部のレンガ構造の概略斜視図である。図1に示すように、コークス炉は、石炭をコークスに乾留するための多数の炭化室20と、石炭加熱用の燃料ガスを燃焼するための多数の燃焼室21とを交互に炉幅方向に配置し構成されている。図2に示すように、炭化室20は、燃焼室21に対して炉壁22で仕切られた空間であり、燃焼室21は燃焼ガスが通るフリュー23を持つレンガと目地の構造体である。そして、燃焼室21と炭化室20の炉体レンガは金物(図示せず)によって表面を覆われ、金物をバックステー24が支持し、これを炉頂部や炉下部に配置したテンションロッド(図示せず)で締め付けて炉体を保持している。乾留したコークスは押出し機25によって炉外に押し出される。   An example of the brick structure is a coke oven (a chamber type coke oven). FIG. 1 is a schematic perspective view of a coke oven, and FIG. 2 is a schematic perspective view of a brick structure inside the furnace body. As shown in FIG. 1, the coke oven has a number of carbonization chambers 20 for carbonizing coal into coke and a number of combustion chambers 21 for burning coal heating fuel gas alternately in the furnace width direction. Arranged and configured. As shown in FIG. 2, the carbonization chamber 20 is a space partitioned by a furnace wall 22 with respect to the combustion chamber 21, and the combustion chamber 21 is a brick-joint structure having a flue 23 through which combustion gas passes. The furnace body bricks of the combustion chamber 21 and the carbonization chamber 20 are covered with hardware (not shown), and the hardware is supported by the backstay 24, which is arranged on the top of the furnace and the lower part of the furnace (not shown). To hold the furnace body. The carbonized coke is extruded out of the furnace by an extruder 25.

レンガ構造物であるコークス炉は、長年の操業中に押出し負荷や熱負荷が掛かり、次第に目地やレンガに亀裂が発生する。テンションロッドが切れて炉締めが解放された場合にも、レンガ構造が緩んで目地やレンガに亀裂が発生し、損傷が進むと炉壁全体の変形へと進んでいく。この変形が進行するとレンガを積替える大規模な補修を行う必要が生じるため、コークス炉の変形程度の把握は極めて重要である。
従来、コークス炉の変形程度や損傷程度を把握するために、以下のような方法が提案されている。
特許文献1には、炭化室の任意高さにおける長さ方向の炉壁間距離を測定し、設計時の炉壁間距離と比較して、カーボン付着量やレンガ欠損量を判断する方法が提案されている。
Coke ovens, which are brick structures, are subject to extrusion load and heat load during many years of operation, and cracks in joints and bricks gradually develop. Even when the tension rod breaks and the furnace tightening is released, the brick structure is loosened, cracks occur in the joints and bricks, and when damage progresses, the entire furnace wall is deformed. As this deformation progresses, it will be necessary to carry out a large-scale repair of transshipment of bricks, so it is extremely important to understand the degree of deformation of the coke oven.
Conventionally, in order to grasp the degree of deformation and damage of a coke oven, the following methods have been proposed.
Patent Document 1 proposes a method for measuring the distance between furnace walls in the longitudinal direction at an arbitrary height of the carbonization chamber and comparing the distance between furnace walls at the time of design to determine the amount of carbon adhesion and the amount of brick defects. Has been.

また、特許文献2には、電磁波距離計による炭化室壁面までの距離計測結果を、押出しラム移動軌跡を用いて補正して、炭化室壁面の絶対形状を算出し、壁面補修終了後、操業毎に壁面の絶対形状を計測し、補修終了直後の壁面絶対形状の変化から壁面へのカーボン付着状態あるいは壁面の損耗状態を定量的に把握する方法が提案されている。
また、特許文献3には、コークス炉の燃焼室に観察ランスを挿入し、燃焼室の壁面を観察ランスにより撮像して欠陥部画像を得る方法が提案されている。
In Patent Document 2, the distance measurement result to the wall surface of the carbonization chamber by the electromagnetic distance meter is corrected using the extrusion ram movement trajectory to calculate the absolute shape of the wall surface of the carbonization chamber. A method has been proposed in which the absolute shape of the wall surface is measured and the carbon adhering state or the worn state of the wall surface is quantitatively determined from the change in the absolute wall shape immediately after the repair is completed.
Patent Document 3 proposes a method in which an observation lance is inserted into a combustion chamber of a coke oven and a wall surface of the combustion chamber is imaged with the observation lance to obtain a defect portion image.

特開2007−332382号公報JP 2007-332382 A 特開2006−36958号公報JP 2006-36958 A 特開2002−226862号公報JP 2002-226862 A

上記のように、コークス炉の変形程度や損傷程度を把握する方法はいくつか提案されており、その手法は、主に距離計による測定若しくは観察装置による撮像である。しかしながら、特許文献1の方法では、コークス炉が老朽化すると設計時と基準位置がずれている場合があり、設計時の距離を基準とすることは必ずしも正確ではない。また、特許文献2の方法では、壁面補修後に毎度壁面の絶対形状を測定する必要があり、労力のかかる方法である。また、特許文献3の方法では、撮像からの判断になり、定量的な把握が困難な場合がある。   As described above, several methods for grasping the degree of deformation and damage of the coke oven have been proposed, and the method is mainly measurement with a distance meter or imaging with an observation device. However, in the method of Patent Document 1, when the coke oven is aged, the reference position may be shifted from that at the time of design, and it is not always accurate to use the distance at the time of design as a reference. Moreover, in the method of patent document 2, it is necessary to measure the absolute shape of the wall surface every time after repairing the wall surface, which is a laborious method. Further, in the method of Patent Document 3, it is determined from imaging, and it may be difficult to quantitatively grasp.

したがって本発明の目的は、以上のような従来技術の課題を解決し、レンガ構造物の変形量を少ない労力で精度よく算出することができる算出方法を提供することにある。   Accordingly, an object of the present invention is to solve the above-described problems of the prior art and to provide a calculation method capable of accurately calculating the deformation amount of a brick structure with little effort.

本発明者らは、目地を間に含む2つ以上のレンガからなり、そのレンガ及び目地が、解析対象となるレンガ構造物のレンガ及び目地と同等の構成を有するレンガ構成体について、予めその機械物性値を求め、その機械物性値を利用して構造解析を行うことにより、コークス炉などのレンガ構造物の変形量を精度よく算出できることを見出した。
本発明は、このような知見に基づきなされたもので、その要旨は以下のとおりである。
The inventors of the present invention have previously made a machine for a brick structure having two or more bricks including joints, the bricks and joints having a structure equivalent to the bricks and joints of the brick structure to be analyzed. It has been found that the amount of deformation of a brick structure such as a coke oven can be accurately calculated by obtaining the physical property value and performing structural analysis using the mechanical property value.
The present invention has been made based on such findings, and the gist thereof is as follows.

[1]レンガ構造物を形状モデル化し、該レンガ構造物の機械物性値及び境界条件を与えて有限要素法による構造解析を行い、レンガ構造物の変形量を算出する方法であって、レンガ構造物がコークス炉であり、該レンガ構造物の機械物性値として、目地(b)を間に含む2つ以上のレンガ(a)からなり、該レンガ(a)及び目地(b)がレンガ構造物のレンガ及び目地と同等の構成を有するレンガ構成体(A)を用いた試験により予め求めた値を用いることを特徴とするレンガ構造物の変形量算出方法。
[2]上記[1]の変形量算出方法において、構造解析において算出される主応力が引張強度を超えると主応力方向と直交する方向にひび割れが発生すると判定し、ひび割れと直交する方向には軟化させる引張挙動の構成則を適用することを特徴とするレンガ構造物の変形量算出方法。
[1] A method of calculating a deformation amount of a brick structure by modeling a shape of the brick structure, performing a structural analysis by a finite element method by giving mechanical property values and boundary conditions of the brick structure, The object is a coke oven, and the mechanical properties of the brick structure are two or more bricks (a) including joints (b). The bricks (a) and joints (b) are brick structures. The deformation | transformation amount calculation method of the brick structure characterized by using the value calculated | required previously by the test using the brick structure (A) which has a structure equivalent to a brick and a joint.
[2] In the deformation amount calculation method of [1] above, if the principal stress calculated in the structural analysis exceeds the tensile strength, it is determined that a crack will occur in the direction perpendicular to the principal stress direction, and in the direction perpendicular to the crack, A method for calculating a deformation amount of a brick structure, wherein a constitutive law of tensile behavior to be softened is applied .

本発明によれば、目地(b)を間に含む2つ以上のレンガ(a)からなり、該レンガ(a)及び目地(b)が、解析対象となるレンガ構造物のレンガ及び目地と同等の構成を有するレンガ構成体(A)について、予めその機械物性値を求め、その機械物性値を利用して構造解析を行うことにより、コークス炉などのレンガ構造物の変形量を少ない労力で精度よく算出することができる。さらに、構造解析において、レンガ構造物の主応力が引張強度を超えた場合には軟化する特性を含むようにすることにより、レンガ構造物の変形量をより精度よく算出することができる。   According to the present invention, the brick (a) is composed of two or more bricks including the joint (b), and the brick (a) and joint (b) are equivalent to the brick and joint of the brick structure to be analyzed. For the brick structure (A) having the following structure, the mechanical property value is obtained in advance, and the mechanical property value is used to perform the structural analysis, thereby reducing the amount of deformation of the brick structure such as a coke oven with less effort. It can be calculated well. Furthermore, in the structural analysis, the amount of deformation of the brick structure can be calculated more accurately by including a characteristic that softens when the main stress of the brick structure exceeds the tensile strength.

室炉式コークス炉の概略斜視図Schematic perspective view of chamber furnace coke oven 室炉式コークス炉の炉体内部におけるレンガ構造の概略斜視図Schematic perspective view of the brick structure inside the furnace body of a chamber-type coke oven 本発明において、機械物性値を求めるために用いるレンガ構成体の一例を示す説明図In this invention, explanatory drawing which shows an example of the brick structure used in order to obtain | require a mechanical property value 本発明で用いる圧縮挙動における構成則の一例を示す図面Drawing which shows an example of the constitutive law in the compression behavior used by this invention 本発明で用いる引張挙動における構成則の一例を示す図面Drawing which shows an example of the constitutive law in the tensile behavior used by this invention 本発明のレンガ構造物の変形量算出方法のフローチャートThe flowchart of the deformation amount calculation method of the brick structure of the present invention 本発明法の実施に供される演算システム(変形量算出装置)の一例を示す説明図Explanatory drawing which shows an example of the arithmetic system (deformation amount calculation apparatus) which is provided for implementation of the method of the present invention 実施例におけるコークス炉解析モデルの図面Drawing of coke oven analysis model in the embodiment

以下、レンガ構造物がコークス炉である場合を例に説明する。
本発明者らは、コークス炉は高温で、しかもセンサーなどを挿入したり、設置したりできる空間が限られているため、コークス炉の変形量の把握には構造解析技術を適用することが適当であると考えた。すなわち、この方法は、コークス炉のなかで変形量を算出したい領域(レンガ構造を含む炉体部分)を形状モデル化し、その領域の機械物性値及び境界条件などを与えて有限要素法による構造解析を行い、コークス炉の当該特定領域の変形量を算出するものである。しかしながら、コークス炉は主にレンガと目地を組み合わせた構造体であり、その材質の非一様性から、目地切れやレンガの亀裂などの局所的な損傷挙動が見られ、コークス炉の損傷の予測は困難である。コークス炉のひとつの窯を考えた場合であっても、レンガと目地の挙動を個別に取り扱って評価することはほとんど不可能である。
Hereinafter, a case where the brick structure is a coke oven will be described as an example.
Since the coke oven is hot and the space in which sensors can be inserted and installed is limited, it is appropriate to apply a structural analysis technique to grasp the amount of deformation of the coke oven. I thought. In other words, in this method, the area of the coke oven where the amount of deformation is to be calculated (furnace body part including the brick structure) is modeled, and the mechanical properties and boundary conditions of that area are given to analyze the structure using the finite element method. To calculate the deformation amount of the specific region of the coke oven. However, the coke oven is a structure that mainly combines bricks and joints, and due to the non-uniformity of the material, local damage behavior such as joint breaks and cracks in bricks can be seen, and prediction of damage to the coke ovens It is difficult. Even when considering a single coke oven, it is almost impossible to handle and evaluate the behavior of bricks and joints separately.

そこで、本発明では、構造解析によりコークス炉の変形量を算出するに当たり、コークス炉の機械物性値として、目地を間に含む2つ以上のレンガからなるレンガ構成体を用いて予め求めた値を利用することを考えた。レンガ構造物の変形の度合いはレンガの大きさ、材質、ダボの有無、目地の厚さなどに大きく依存するので、2つ以上のレンガと目地が一体となったレンガ構造体で機械物性を評価し、その機械物性を構造解析に利用することが、コークス炉などのレンガ構造物の変形量を精度よく算出するのに有効であると考えられる。図3は、このように予め機械物性を求めるレンガ構成体Aの一例を示している。このレンガ構成体Aは、目地bを間に含む2つのレンガaからなり、各レンガaはダボc(及びダボ穴)を有している。レンガa及び目地bは、解析対象となるコークス炉(レンガ構造物)で使用されているレンガ及び目地と同等の構成(レンガのサイズ・材質・ダボの有無及び構成、目地の厚さなど)を有する。ここで、レンガa及び目地bが、解析対象となるコークス炉(レンガ構造物)で使用されているレンガ及び目地と同等の構成を有するとは、レンガ及び目地の材質が解析対象となる領域のレンガ構造物と同じであり、且つレンガの寸法(ダボの有無及び構成を含む)及び目地の厚みが当該領域のレンガ構造物を代表する寸法及び厚み(つまり、解析対象となる領域のレンガ構造物において最も多く使われているレンガの寸法と目地の厚み)であることを意味する。1つのコークス炉に使用されるレンガの形状や寸法は全て同じである訳ではなく、また、目地の厚みも部分的に変更されているのが一般的だからである。
レンガ構成体Aは、目地bを間に含む3つ以上のレンガaで構成してもよい。ただし、技術的・経済的に容易に試験を実施できる「レンガ構成体A」を用いるという観点からして、レンガ構成体Aを構成するレンガaの数は多くても10個程度とするのが適当であり、普通には2〜3個程度である。
Therefore, in the present invention, when calculating the deformation amount of the coke oven by structural analysis, as a mechanical property value of the coke oven, a value obtained in advance using a brick structure composed of two or more bricks including joints between them is used. I thought about using it. The degree of deformation of a brick structure greatly depends on the size, material, presence of dowels, joint thickness, etc., so the mechanical properties are evaluated with a brick structure in which two or more bricks and joints are integrated. However, it is considered that using the mechanical properties for structural analysis is effective in accurately calculating the deformation amount of a brick structure such as a coke oven. FIG. 3 shows an example of the brick structure A for which the mechanical properties are obtained in advance as described above. The brick structure A is composed of two bricks a including a joint b, and each brick a has a dowel c (and a dowel hole). Brick a and joint b have the same structure as the brick and joint used in the coke oven (brick structure) to be analyzed (brick size, material, presence and configuration of dowels, joint thickness, etc.) Have. Here, the bricks a and joints b have the same configuration as the bricks and joints used in the coke oven (brick structure) to be analyzed. It is the same as the brick structure, and the size and thickness of the brick (including the presence / absence and configuration of the dowel) and the joint thickness are representative of the brick structure in the region (that is, the brick structure in the region to be analyzed) The most commonly used brick size and joint thickness). This is because the shapes and dimensions of the bricks used in one coke oven are not all the same, and the joint thickness is generally partially changed.
The brick structure A may be composed of three or more bricks a including the joint b therebetween. However, from the viewpoint of using a “brick structure A” that can be easily tested technically and economically, the number of bricks a constituting the brick structure A is about 10 at most. Appropriate, usually about 2-3.

このレンガ構成体Aで求められる機械物性値には、例えば、密度、弾性係数、ポアソン比、圧縮強度、引張強度などが含まれる。また、対象がコークス炉などのように高温で使用され、温度変化が大きいレンガ構造物の場合には、熱膨張率や各物性の温度依存性を考慮することが望ましく、例えば、幾つかの代表温度での機械物性を含めることが好ましい。必要であれば、温度物性(熱伝導率や比熱など)を与えて温度解析を実施することが望ましい。これらの機械物性値は、主に、レンガ構成体Aについて圧縮試験、引張試験などの必要な試験を実施することで求められる。   The mechanical property values required for the brick structure A include, for example, density, elastic modulus, Poisson's ratio, compressive strength, tensile strength, and the like. In addition, when the object is a brick structure that is used at a high temperature such as a coke oven and has a large temperature change, it is desirable to consider the coefficient of thermal expansion and the temperature dependence of each property. It is preferable to include mechanical properties at temperature. If necessary, it is desirable to perform temperature analysis by giving temperature properties (thermal conductivity, specific heat, etc.). These mechanical property values are mainly determined by performing necessary tests such as a compression test and a tensile test on the brick component A.

さらに、上記の機械物性値から圧縮挙動および引張挙動における構成則を決定し、これを構造解析に利用することが望ましい。構成則の一例を図4及び図5に示す。図4は圧縮挙動における構成則を、図5は引張挙動における構成則を、それぞれ示している。この例では、機械物性値として、弾性係数Eo、圧縮強度fc、引張強度ftに加えて、極限圧縮ひずみεu、極限引張ひずみεmを用い、応力−ひずみ関係を決定する。ここで、極限圧縮ひずみεuは、最大圧縮応力値に対応するひずみであり、その数値は、構成則のモデルによって異なるが、例えば図4のモデルの場合、経験的に0.001〜0.05の範囲で決められる。また、極限引張ひずみεmは、破壊時に対応する引張ひずみであり、その数値は、構成則のモデルによって異なるが、例えば図5のモデルの場合、経験的に0.1〜0.5の範囲で決められる。
圧縮挙動は、図4のように、作用する圧縮応力σが大きくなればひずみεも大きくなるように仮定し、極限に達すれば材料が圧壊する。また、引張挙動は、図5のように、引張強度ftに達するまでは引張応力σが大きくなればひずみεも大きくなる。引張強度ftに達すると材料がひび割れ、引張が作用してもひずみεだけが大きくなり、極限に達すると破断する。
引張挙動における構成則では、引張強度ftを超えた場合には、ひび割れが発生し、軟化する特性が考慮されていることが望ましい。
Furthermore, it is desirable to determine the constitutive law in the compression behavior and the tensile behavior from the above mechanical property values, and use this for the structural analysis. An example of the constitutive law is shown in FIGS. FIG. 4 shows a constitutive law for compression behavior, and FIG. 5 shows a constitutive law for tensile behavior. In this example, the stress-strain relationship is determined by using the ultimate compressive strain εu and the ultimate tensile strain εm in addition to the elastic modulus Eo, the compressive strength fc, and the tensile strength ft as mechanical property values. Here, the ultimate compressive strain εu is a strain corresponding to the maximum compressive stress value, and the numerical value varies depending on the model of the constitutive law. For example, in the case of the model of FIG. It is decided in the range. Further, the ultimate tensile strain εm is a tensile strain corresponding to at the time of fracture, and its numerical value varies depending on the model of the constitutive law. For example, in the case of the model of FIG. It is decided.
As shown in FIG. 4, it is assumed that the compressive behavior is such that the strain ε increases as the compressive stress σ acting increases, and the material collapses when reaching the limit. In addition, as shown in FIG. 5, the tensile behavior increases as the tensile stress σ increases until the tensile strength ft is reached, as shown in FIG. When the tensile strength ft is reached, the material cracks, and even if tension is applied, only the strain ε increases, and when the limit is reached, the material breaks.
According to the constitutive law in the tensile behavior, it is desirable to take into account the property of cracking and softening when the tensile strength exceeds ft.

本発明では、コークス炉(レンガ構造物)を形状モデル化し、有限要素法による構造解析を行い、コークス炉の変形量を算出するものであるが、この構造解析は、上述したレンガ構成体Aにおいて予め求めた機械物性値(さらに好ましくは、この機械物性値から決定された構成則)を利用して行われる。
図6に、本発明法のフローチャートを示す。本発明法では、まず、対象となるコークス炉のなかで解析によって変形量を算出したい領域(レンガ構造を含む炉体部分)を選び、その領域を形状モデル化する。解析の対象となる領域とは、例えば、1つの燃焼室や天井部分などである。この形状モデル化では、例えば、適切なプリプロセッサーを使用してその形状を作成し、計算メッシュを作成する。
In the present invention, the shape of the coke oven (brick structure) is modeled, the structural analysis is performed by the finite element method, and the amount of deformation of the coke oven is calculated. This is performed using a mechanical property value obtained in advance (more preferably, a constitutive law determined from this mechanical property value).
FIG. 6 shows a flowchart of the method of the present invention. In the method of the present invention, first, a region (furnace body portion including a brick structure) for which the amount of deformation is to be calculated by analysis is selected from the target coke oven, and the region is modeled. The region to be analyzed is, for example, one combustion chamber or a ceiling portion. In this shape modeling, for example, the shape is created using an appropriate preprocessor, and a calculation mesh is created.

次に、形状モデルに含まれる炉体部分(レンガ構造部)および付帯構造物(レンガ構造を拘束するバックステーや金物など)の機械物性値、さらに好ましくは機械物性値に基づいて決定された圧縮挙動及び引張挙動の構成則を入力する。炉体部分(レンガ構造部)の機械物性値と圧縮挙動及び引張挙動の構成則については、上述したとおりである。
また、コークス炉は高温で使用されるため熱膨張率、熱伝導率などの温度物性も入力する。
Next, compression determined based on the mechanical property values of the furnace body part (brick structure part) and incidental structures (backstays and hardware that restrain the brick structure) included in the shape model, more preferably based on the mechanical property values Enter constitutive laws for behavior and tensile behavior. The constitutive rules of the mechanical property value, compression behavior, and tensile behavior of the furnace body portion (brick structure portion) are as described above.
Also, since coke ovens are used at high temperatures, temperature properties such as thermal expansion coefficient and thermal conductivity are also input.

さらに、解析に必要な初期条件、境界条件として、温度条件、拘束条件、外力を入力する。温度条件は定常操業時には等温で扱ってよい場合もあるが、例えば、(i)補修時に一時操業を停止した場合には、窯の温度が数百℃降下した後に操業温度1200℃程度まで上昇すること、(ii)温度分布が一様でなく、部分的に温度が低い(例えば端部のみが温度が低い)場合があること、などが考えられ、温度条件は重要になる。また、拘束条件とは、解析対象となる炉体部分(レンガ構造部)をその位置に留め置いている付帯構造物や他の炉体部分などによる拘束条件のことである。特に、上述したようにコークス炉はレンガ構造が金物やバックステーによって拘束されており、どのような力で締め付けが行われているか(例えば、何らかの理由で締め付け力が緩むような場合を含む)が計算には必要になる。また、石炭の乾留中には炉壁に石炭膨張圧が作用し、コークス押出し時には炉壁に押出し負荷がかかり、これらの作用も外力として考慮する。また、以上のような温度条件や拘束条件、外力は、窯に石炭を装入した直後か否か、乾留の末期か否か、コークス押出し時か否か、補修時か否か、隣窯が空窯か否か、といった操業条件に応じて変わってくるので、それらの違いに応じて温度条件や拘束条件、外力を決める。   Furthermore, temperature conditions, constraint conditions, and external forces are input as initial conditions and boundary conditions necessary for the analysis. Although the temperature condition may be handled isothermally during steady operation, for example, (i) When temporary operation is stopped during repair, the temperature of the kiln rises to about 1200 ° C after the temperature drops by several hundred degrees (Ii) the temperature distribution is not uniform and the temperature may be partially low (for example, only the edge is low in temperature). In addition, the constraint condition is a constraint condition by an incidental structure or another furnace body part that holds the furnace body part (brick structure part) to be analyzed at the position. In particular, as described above, in the coke oven, the brick structure is constrained by hardware or a backstay, and what kind of force is used for tightening (including the case where the tightening force is loosened for some reason). Necessary for calculation. In addition, during the carbonization of coal, the coal expansion pressure acts on the furnace wall, and during the coke extrusion, an extrusion load is applied to the furnace wall, and these actions are also considered as external forces. In addition, the temperature conditions, restraint conditions, and external force as described above are determined immediately after charging coal in the kiln, whether it is at the end of dry distillation, whether it is during coke extrusion, whether it is repaired, Since it changes depending on the operating conditions such as whether or not it is an empty kiln, temperature conditions, restraint conditions, and external force are determined according to these differences.

以上の各項目を入力した上で、有限要素法によるコークス炉の構造解析を実施する。この有限要素法による構造解析は、上記のような各項目の入力に基づき、コンピュータ上で汎用的な構造解析用プログラムを用いて実行される。この構造解析では、ひび割れの発生を判定し、このひび割れによる部材の軟化を考慮して変形量を算出することが好ましい。すなわち、算出される主応力が引張強度を超えると主応力方向と直交する方向にひび割れが発生すると判定し、ひび割れ発生まで或いはひび割れと平行する方向には弾性体として扱うが、ひび割れと直交する方向には軟化させる引張挙動の構成則を適用する。ひび割れの発生により変形量は著しく大きくなるため、ひび割れを判定し、ひび割れ後の軟化特性を考慮することは重要である。   After inputting the above items, the structural analysis of the coke oven by the finite element method is performed. The structural analysis by the finite element method is executed on a computer using a general-purpose structural analysis program based on the input of each item as described above. In this structural analysis, it is preferable to determine the occurrence of a crack and calculate the deformation amount in consideration of the softening of the member due to the crack. That is, if the calculated principal stress exceeds the tensile strength, it is determined that a crack will occur in the direction perpendicular to the principal stress direction, and the crack is generated as an elastic body until the occurrence of the crack or in the direction parallel to the crack, but the direction perpendicular to the crack The constitutive law of the tensile behavior to be softened is applied. Since the amount of deformation is significantly increased due to the occurrence of cracks, it is important to determine cracks and consider the softening properties after cracks.

繰り返し計算にて解が得られた後、コークス炉の変形を出力する。出力方法は任意であり、例えば、コークス炉の形状モデルを変形させた可視化表示であっても、コークス炉形状モデル上に色で区別した等変形線図であっても、X−Yグラフのような出力であってもよい。また、解析で得られる応力などの他のパラメータを表示してもよい。適切なポストプロセッサーを用いて実施してもよい。   After the solution is obtained by repeated calculations, the deformation of the coke oven is output. The output method is arbitrary. For example, even a visualization display obtained by deforming a coke oven shape model, or an iso-deformation diagram distinguished by color on a coke oven shape model, as in an XY graph Output may be sufficient. Moreover, you may display other parameters, such as the stress obtained by analysis. It may be implemented using a suitable post processor.

本発明に実施に供される演算システム(算出装置)の一例を図7に示す。この演算システム(算出装置)は、物性値(機械物性、温度物性)及び初期条件・境界条件(温度条件、拘束条件、外力)を入力する入力部1と、入力された物性値(機械物性、温度物性)及び初期条件・境界条件(温度条件、拘束条件、外力)に基づいて有限要素法による構造解析を実行し、変形量などを演算する演算処理部2と、この演算処理部2で得られた変形量などを出力する出力部3を備えている。
本発明は、コークス炉に限らず、加熱炉や転炉などをはじめとするあらゆる種類のレンガ構造物を対象とすることができる。
An example of an arithmetic system (calculation device) provided for implementation in the present invention is shown in FIG. This computing system (calculation device) includes an input unit 1 for inputting physical property values (mechanical properties, temperature physical properties) and initial conditions / boundary conditions (temperature conditions, constraint conditions, external force), and input physical property values (mechanical physical properties, The calculation processing unit 2 performs structural analysis by the finite element method based on the temperature physical properties) and the initial conditions / boundary conditions (temperature conditions, constraint conditions, external force), and calculates the deformation amount. An output unit 3 is provided for outputting the deformation amount and the like.
The present invention is not limited to a coke oven, but can be applied to all types of brick structures including a heating furnace and a converter.

老朽化した室炉式コークス炉において、保守点検のために、ある窯を空窯にして炉蓋が開放された。操業影響にて炉蓋開放時間が長くなり、通常温度が約1200℃のところ約800℃程度まで温度低下が見られた。再度昇温して約1200℃にしたところ、炉壁レンガが部分的に約35mm張り出した。張り出し箇所は奥行きが窯口から2m、高さ方向が炉底から5mの付近であった。張り出し部分の窯口炉壁レンガを回収して調べたところ、本来100mmのはずの厚みが45mmしかなかった。   In an aging furnace-type coke oven, a certain kiln was made empty and the furnace lid was opened for maintenance. Due to operational effects, the furnace lid opening time became longer, and when the normal temperature was about 1200 ° C., the temperature dropped to about 800 ° C. When the temperature was raised again to about 1200 ° C., the furnace wall brick partially protruded about 35 mm. The overhanging portion had a depth of 2 m from the kiln entrance and a height direction of about 5 m from the furnace bottom. When the overhanging furnace wall bricks of the overhanging part were collected and examined, the original thickness that should have been 100 mm was only 45 mm.

本発明の効果を確認するために、燃焼室1つを解析領域として形状モデル化した。この形状モデル化では、所定のプリプロセッサーを使用してその形状を作成し、計算メッシュを作成した。その形状を図8に示す。炉壁で厚みが減少していた部分(図8参照)は、厚み45mmとした。図6のフローに従い、物性値および境界条件等を入力し、有限要素法による構造解析を実施した。付帯構造物(金物)関して使用した物性値、ケース1(比較例)及びケース2(本発明例)で使用した物性値を表1に示す。また、境界条件を図8に示す。天井部と底面はY方向への変位をゼロ(Uy=0)として拘束した。端部下端では炉長方向には伸びない(Ux=0)とした。温度については、事前に温度の伝熱解析を行った。外気温度30℃と仮定し、炉頂、炉端部で熱伝達係数50W/mKとした。炉底部は断熱とした。燃焼室温度が炉蓋開放時800℃から通常操業時1200℃になったときの伝熱解析結果から熱膨張量を計算し、熱応力として構造解析で考慮した。 In order to confirm the effect of the present invention, one combustion chamber was modeled as an analysis region. In this shape modeling, the shape was created using a predetermined preprocessor, and a calculation mesh was created. The shape is shown in FIG. The portion where the thickness was reduced on the furnace wall (see FIG. 8) was 45 mm in thickness. According to the flow of FIG. 6, physical property values, boundary conditions and the like were input, and structural analysis was performed by a finite element method. Table 1 shows the physical property values used for the incidental structure (metal) and the physical property values used in Case 1 (Comparative Example) and Case 2 (Invention Example). The boundary conditions are shown in FIG. The ceiling and bottom surfaces were constrained with zero displacement (Yy = 0) in the Y direction. The lower end of the end portion was not extended in the furnace length direction (Ux = 0). Regarding temperature, heat transfer analysis of temperature was performed in advance. Assuming an outside air temperature of 30 ° C., the heat transfer coefficient was set to 50 W / m 2 K at the top and end of the furnace. The bottom of the furnace was insulated. The amount of thermal expansion was calculated from the heat transfer analysis results when the combustion chamber temperature was 800 ° C. when the furnace lid was opened and 1200 ° C. during normal operation, and was taken into account in the structural analysis as thermal stress.

ケース1(比較例)では、引張強度と極限引張ひずみは目地の物性値を利用し、それ以外はレンガの物性値を利用し、構造解析を実施した。
ケース2(本発明例)では、図3に示すような目地bを間に含む2つのレンガaからなるレンガ構成体Aにより表1に示す機械物性値を求めた。レンガ構成体Aのレンガa及び目地bは、コークス炉のレンガ及び目地と同等の構成を有するものであり、レンガaの寸法は、長さL=200mm、幅W=100mm、高さH=130mm、目地bは厚みT=4mmであった。目地部にはダボcがあり、凸部はR12mmであった。また、レンガおよび目地(モルタル)の材質は珪石質であった。
ケース1,2ともに、圧縮挙動および引張挙動における構成則は図4、図5に示す形態とした。
In Case 1 (Comparative Example), the tensile strength and the ultimate tensile strain utilized the physical property values of the joints, and other than that, the physical property values of the bricks were utilized for structural analysis.
In case 2 (example of the present invention), the mechanical property values shown in Table 1 were obtained from a brick structure A composed of two bricks a including a joint b as shown in FIG. The brick a and the joint b of the brick structure A have the same configuration as the brick and joint of the coke oven, and the dimensions of the brick a are as follows: length L = 200 mm, width W = 100 mm, height H = 130 mm. The joint b had a thickness T = 4 mm. The joint has a dowel c, and the convex portion was R12 mm. The material of bricks and joints (mortar) was siliceous.
In both cases 1 and 2, the constitutive rules for the compression behavior and the tensile behavior are shown in FIGS.

ケース1(比較例)の場合には7mmの張り出しが計算されたのに対して、ケース2(本発明例)では33mmの張り出しが計算された。ケース2が実際のコークス炉の変形を精度よく再現できていることが判った。
このように、目地bを間に含む2つ以上のレンガaからなるレンガ構成体A(レンガa及び目地bはレンガ構造物のレンガ及び目地と同じ構成を有する)の機械物性を利用して構造解析を行うことで、コークス炉の変形を精度よく把握できることが判った。
In case 1 (comparative example), a 7 mm overhang was calculated, whereas in case 2 (invention example), a 33 mm overhang was calculated. It was found that Case 2 was able to accurately reproduce the actual deformation of the coke oven.
As described above, the structure is made by utilizing the mechanical properties of the brick structure A composed of two or more bricks a including the joints b (the bricks a and joints b have the same configuration as the bricks and joints of the brick structure). The analysis revealed that the deformation of the coke oven can be grasped with high accuracy.

Figure 0005994086
Figure 0005994086

A レンガ構成体
a レンガ
b 目地
c ダボ
1 入力部
2 演算処理部
3 出力部
A brick structure a brick b joint c dowel 1 input unit 2 arithmetic processing unit 3 output unit

Claims (2)

レンガ構造物を形状モデル化し、該レンガ構造物の機械物性値及び境界条件を与えて有限要素法による構造解析を行い、レンガ構造物の変形量を算出する方法であって、
レンガ構造物がコークス炉であり、該レンガ構造物の機械物性値として、目地(b)を間に含む2つ以上のレンガ(a)からなり、該レンガ(a)及び目地(b)がレンガ構造物のレンガ及び目地と同等の構成を有するレンガ構成体(A)を用いた試験により予め求めた値を用いることを特徴とするレンガ構造物の変形量算出方法。
The brick structure is modeled, the mechanical property value of the brick structure and boundary conditions are given, the structure analysis is performed by the finite element method, and the deformation amount of the brick structure is calculated,
The brick structure is a coke oven, and the mechanical properties of the brick structure are two or more bricks (a) including joints (b). The bricks (a) and joints (b) are bricks. A method for calculating a deformation amount of a brick structure, wherein a value obtained in advance by a test using a brick structure (A) having a structure equivalent to a brick and joints of the structure is used.
構造解析において算出される主応力が引張強度を超えると主応力方向と直交する方向にひび割れが発生すると判定し、ひび割れと直交する方向には軟化させる引張挙動の構成則を適用することを特徴とする請求項1に記載のレンガ構造物の変形量算出方法。   When the principal stress calculated in the structural analysis exceeds the tensile strength, it is judged that cracks will occur in the direction perpendicular to the principal stress direction, and the constitutive law of tensile behavior that softens in the direction perpendicular to the crack is applied. The deformation | transformation amount calculation method of the brick structure of Claim 1 to do.
JP2013082919A 2013-04-11 2013-04-11 Deformation calculation method for brick structure Active JP5994086B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013082919A JP5994086B2 (en) 2013-04-11 2013-04-11 Deformation calculation method for brick structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013082919A JP5994086B2 (en) 2013-04-11 2013-04-11 Deformation calculation method for brick structure

Publications (2)

Publication Number Publication Date
JP2014205741A JP2014205741A (en) 2014-10-30
JP5994086B2 true JP5994086B2 (en) 2016-09-21

Family

ID=52119578

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013082919A Active JP5994086B2 (en) 2013-04-11 2013-04-11 Deformation calculation method for brick structure

Country Status (1)

Country Link
JP (1) JP5994086B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6350487B2 (en) * 2015-11-06 2018-07-04 Jfeスチール株式会社 Mortar application method, mortar application nozzle, and arm type robot for mortar application

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3723618B2 (en) * 1996-01-09 2005-12-07 ソニー株式会社 Test apparatus and method
JP2001323278A (en) * 2000-05-16 2001-11-22 Nkk Corp Simulation method for coke oven wall deformation in coke oven, coke oven operation method, repair method and design method

Also Published As

Publication number Publication date
JP2014205741A (en) 2014-10-30

Similar Documents

Publication Publication Date Title
Dai et al. Determination of the fracture behaviour of MgO-refractories using multi-cycle wedge splitting test and digital image correlation
Destrebecq et al. Analysis of cracks and deformations in a full scale reinforced concrete beam using a digital image correlation technique
JP6601762B2 (en) Steel heat treatment simulation method and steel heat treatment simulation program
Zhang et al. Creep crack growth behavior analysis of the 9Cr-1Mo steel by a modified creep-damage model
CN110188451A (en) A kind of analysis method of the residual stress of polyvinyl piping materials welding point
Li et al. Unified viscoplastic constitutive model under axial-torsional thermo-mechanical cyclic loading
CN111044351B (en) A method for predicting creep deformation of welded joints based on DIC technology
CN103323351A (en) Cantilever bending load metal material fatigue damage test measurement method
CN108519437A (en) A Multiple Regression Prediction Model of Uniaxial Compressive Strength of Coal Samples and Its Establishment Method
Remadi et al. Prediction of fatigue crack growth life under variable-amplitude loading using finite element analysis
CN117954022A (en) A method for estimating fatigue crack growth rate of metal materials under elastic-plastic conditions
CN104330300B (en) Superhigh temperature ceramic material heat-damage stiffness of coupling indirect measurement method
Brochen et al. Improved thermal stress resistance parameters considering temperature gradients for bricks in refractory linings
JP5994086B2 (en) Deformation calculation method for brick structure
KR101533542B1 (en) Method for tgmf life predicting of thermal barrier coating
JP6339969B2 (en) Deformation resistance identification method for thin inspection materials
Zheng et al. A novel fatigue assessment approach by Direct Steady Cycle Analysis (DSCA) considering the temperature-dependent strain hardening effect
Zheng et al. Axial compression behavior of high-strength fire-resistant steel columns in fire
CN105606255B (en) The prediction technique of metal blank simple tension process temperature variation
Seitl et al. Effect of rivet holes on calibration curves for edge cracks under various loading types in steel bridge structure
JP5983951B2 (en) Blast furnace stave design method
Li et al. Numerical simulation of ratcheting and fatigue behaviour of mitred pipe bends under in-plane bending and internal pressure
Gasser et al. Thermomechanical behaviour analysis and simulation of steel/refractory composite linings
CN118673739A (en) A method for evaluating cracks in arrest steel based on finite element simulation
CN204514837U (en) For the portable test unit of thermal expansion coefficient of concrete test

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20141029

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150821

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150901

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20151029

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160405

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160531

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160621

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160719

R150 Certificate of patent or registration of utility model

Ref document number: 5994086

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250