JP2006247684A - Method, apparatus and program for simulating die-casting, and recording medium for recording the program - Google Patents

Method, apparatus and program for simulating die-casting, and recording medium for recording the program Download PDF

Info

Publication number
JP2006247684A
JP2006247684A JP2005065666A JP2005065666A JP2006247684A JP 2006247684 A JP2006247684 A JP 2006247684A JP 2005065666 A JP2005065666 A JP 2005065666A JP 2005065666 A JP2005065666 A JP 2005065666A JP 2006247684 A JP2006247684 A JP 2006247684A
Authority
JP
Japan
Prior art keywords
solid phase
mold
filling
die casting
cavity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2005065666A
Other languages
Japanese (ja)
Other versions
JP4373353B2 (en
Inventor
Yasushi Iwata
靖 岩田
Yoshio Sugiyama
義雄 杉山
Jushin Tou
ジュシン トウ
Hiroaki Iwabori
弘昭 岩堀
Masafumi Nishida
雅文 西田
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.)
Toyota Motor Corp
Toyota Central R&D Labs Inc
Original Assignee
Toyota Motor Corp
Toyota Central R&D Labs Inc
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 Toyota Motor Corp, Toyota Central R&D Labs Inc filed Critical Toyota Motor Corp
Priority to JP2005065666A priority Critical patent/JP4373353B2/en
Publication of JP2006247684A publication Critical patent/JP2006247684A/en
Application granted granted Critical
Publication of JP4373353B2 publication Critical patent/JP4373353B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Molds, Cores, And Manufacturing Methods Thereof (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a method for simulating a die-casting with which rippled surface defect developed caused by solidification of molten material on the way of filling-up, can further surely be presumed. <P>SOLUTION: This method is performed through the following processes, that is; a pre-treating process S100 containing a factor-making step S101 for dividing the shape of a mold into fine factors by positioning on the coordinate axes, and a factor-defining step S102 for defining the mold factors and cavity factors; a filling-up heat conduction analyzing process S200 containing a filling-up heat conduction analyzing step S201 for performing a filling-up range analysis to respective cavity factors and also, calculating the temperature in the molten material filling-up factors and a solid-phase ratio calculating step S202 for calculating the solid-phase ratio of the molten metal material; and an evaluating process S300 containing a solid-phase thickness calculating step S301 for calculating the solid-phase thickness as the thickness in the cavity factors which continuously becomes the solid-phase 1 from the inner wall surface in the mold factors at the completing time of the filling-up, and a surface-defect developing evaluating step S302 for presuming the development of the surface defect from the solid-phase thickness. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

本発明は、ダイカスト鋳造(以下、単にダイカストという)により得られた鋳物に生じうる表面欠陥を精度よく予測できるダイカストシミュレーションに関する。   The present invention relates to a die casting simulation capable of accurately predicting surface defects that may occur in a casting obtained by die casting (hereinafter simply referred to as die casting).

例えば、アルミニウム合金等よりなる鋳物(成形品)を製造する場合には、鋳型内に溶融した金属材料(以下、適宜、溶融金属材料又は単に溶融材料という)を加圧充填する方法であるダイカストが多用されている。
ダイカストにより得られる鋳物における課題としては、ひけ巣欠陥、空気等の巻き込み欠陥、湯回り不良欠陥、表面欠陥等がある。このうち、特に大きな問題となる欠陥の1つとして表面欠陥がある。
For example, in the case of producing a casting (molded product) made of an aluminum alloy or the like, there is a die casting that is a method of pressurizing and filling a molten metal material (hereinafter, appropriately referred to as a molten metal material or simply a molten material) into a mold. It is used a lot.
Problems in castings obtained by die casting include sinkhole defects, entrainment defects such as air, defective hot water defects, surface defects, and the like. Among these, there is a surface defect as one of the particularly serious defects.

ダイカストでは20〜80MPaと大きな鋳造圧力をかけて鋳物を作製するが、鋳造条件によっては、鋳物表面各部にしわ状の形態を示す表面欠陥(湯じわ)が発生する。このような表面欠陥をなくすためには、成形品(鋳物)の形状、成形方案(ランナ、ゲート、オーバーフロー)、射出条件(低速速度、切り替えタイミング、高速速度等)、型温度制御等を適正化する必要がある。   In die casting, a casting is produced by applying a large casting pressure of 20 to 80 MPa. However, depending on casting conditions, surface defects (wrinkles) showing wrinkle-like shapes are generated in each part of the casting surface. In order to eliminate such surface defects, the shape of the molded product (casting), molding method (runner, gate, overflow), injection conditions (low speed, switching timing, high speed, etc.), mold temperature control, etc. are optimized. There is a need to.

しかしながらダイカストで製造される成形品は、通常3次元的に複雑な形状を有すると共に肉薄であるので、溶融材料の流れ、凝固現象は非常に複雑かつ短時間の現象である。そのため、凝固現象等を理論的に必ずしも正確に解明できず、それゆえ適正条件を見出すことは容易ではない。また、実験的に表面欠陥を系統的に解析することは容易ではなく、試行錯誤を繰り返すのが現状である。   However, since a molded article manufactured by die casting usually has a three-dimensionally complicated shape and is thin, the flow of the molten material and the solidification phenomenon are extremely complicated and short-time phenomena. For this reason, the coagulation phenomenon or the like cannot theoretically be accurately clarified, and therefore it is not easy to find an appropriate condition. Moreover, it is not easy to systematically analyze surface defects experimentally, and the current situation is that trial and error are repeated.

ところで近年、コンピュータの計算能力の向上に伴い、ダイカストにおける溶融材料を型内に充填するときの溶融材料の挙動について、コンピュータを利用して解析する手法、いわばダイカストシミュレーションの適用範囲が広がっている。
これまでのダイカストシミュレーションは、溶融材料の流れ及び流れに伴う凝固の挙動について理解を深めることを目的としており、適正なダイカスト鋳造条件を探索する有用な手段として期待されている。
By the way, in recent years, with the improvement of computer computing power, the range of application of a technique for analyzing the behavior of a molten material when filling the mold with the molten material in a die using a computer, that is, a so-called die casting simulation has been expanded.
The conventional die casting simulation aims to deepen the understanding of the flow of the molten material and the solidification behavior accompanying the flow, and is expected as a useful means for searching for appropriate die casting conditions.

従来のダイカストシミュレーション方法では、溶融材料の流れと凝固の同時解析を行うものが多いが、充填終了後の凝固挙動を解析するものも開発されている。例えば、特許文献1に開示された金属溶湯の流動凝固解析方法は、凝固解析において、型を微小要素に分割し、その微小要素毎に熱伝導及び溶質移動を解析している。   Many conventional die casting simulation methods perform simultaneous analysis of the flow of molten material and solidification, but those that analyze solidification behavior after filling have also been developed. For example, the method for analyzing solidification of molten metal disclosed in Patent Document 1 divides a mold into microelements in the solidification analysis, and analyzes heat conduction and solute movement for each microelement.

また、溶湯充填時の内部欠陥についても開発が行われており、例えば特許文献2に開示された数値解析を利用した鋳造欠陥の予測方法は、型を微小要素に分割して、その微小要素毎に対して所定時間毎に溶湯温度、溶湯圧力及びガス圧力を求め、溶湯圧力よりもガス圧力が高い部位に空孔(ひけ巣欠陥)が発生することを予測する。
しかしながら、従来のダイカストシミュレーションにおける充填、凝固時の解析は、内部欠陥を予測することが主体となっており、表面欠陥の発生を子測する手法は未だ開発されていない。
Further, internal defects at the time of filling a molten metal have also been developed. For example, a casting defect prediction method using numerical analysis disclosed in Patent Document 2 divides a mold into minute elements, and each minute element is divided. Then, the molten metal temperature, molten metal pressure and gas pressure are obtained every predetermined time, and it is predicted that vacancies (shrinkage defects) will be generated at a site where the gas pressure is higher than the molten metal pressure.
However, the analysis at the time of filling and solidification in the conventional die casting simulation is mainly for predicting internal defects, and a method for measuring the occurrence of surface defects has not been developed yet.

特開平10−34320号公報JP-A-10-34320 特開平8−257741号公報JP-A-8-257741

本発明は、かかる従来の問題点に鑑みてなされたもので、溶融材料が充填途中で凝固することに起因して発生する湯じわ状の表面欠陥について、より正確に推測できるダイカストシミュレーション方法、その装置及びそのプログラム、並びに当該プログラムを記録した記録媒体を提供しようとするものである。   The present invention has been made in view of such conventional problems, and a die casting simulation method that can more accurately estimate the surface defects of the wrinkle generated due to solidification of the molten material during filling, The present invention intends to provide the apparatus, the program thereof, and a recording medium on which the program is recorded.

第1の発明は、鋳型内に溶融金属材料を加圧注入して凝固させることにより所望形状の鋳物を得るダイカストを行う際に、少なくとも上記鋳物に生じる表面欠陥の発生を予測するダイカストシミュレーション方法であって、
上記溶融金属材料の成形に用いる上記鋳型の形状を座標系上に位置づけ、該座標系の空間を複数の微小要素に分割する要素作成ステップと、
上記微小要素のそれぞれについて、上記鋳型の型領域に位置する場合には型要素と、上記鋳型のキャビティ領域に位置する場合にはキャビティ要素と定義する要素定義ステップと、を含む前処理工程と、
上記キャビティ要素のそれぞれについて、上記溶融金属材料の充填状態を経時的に解析する充填領域解析を行うと共に、上記溶融金属材料が充填された上記キャビティ要素である溶融材料充填要素について、該溶融材料充填要素同士の間、及び該溶融材料充填要素と上記型要素との間の伝熱をそれぞれの熱伝達係数に基づいて経時的に解析して、上記溶融材料充填要素の温度を経時的に算出する充填伝熱解析ステップと、
算出された上記溶融材料充填要素の温度に応じて、該溶融材料充填要素内の上記溶融金属材料の固相率を経時的に算出する固相率算出ステップと、を含む充填伝熱解析工程と、
充填完了時における、上記型要素の内壁面から連続的に固相率1となる上記キャビティ要素の厚さである固体相厚さを算出する固体相厚さ算出ステップと、
該固体相厚さから表面欠陥の発生を予測する表面欠陥発生評価ステップと、を含む評価工程とを有することを特徴とするダイカストシミュレーション方法にある(請求項1)。
The first invention is a die casting simulation method for predicting at least the occurrence of surface defects in the casting when performing die casting to obtain a casting having a desired shape by pressurizing and solidifying a molten metal material into a mold. There,
An element creating step of positioning the shape of the mold used for forming the molten metal material on a coordinate system and dividing the space of the coordinate system into a plurality of microelements;
For each of the microelements, a preprocessing step including a mold element when positioned in the mold region of the mold and an element defining step that defines a cavity element when positioned in the cavity region of the mold;
For each of the cavity elements, a filling region analysis for analyzing the filling state of the molten metal material over time is performed, and the molten material filling element that is the cavity element filled with the molten metal material is filled with the molten material. The heat transfer between the elements and between the molten material filling element and the mold element is analyzed over time based on the respective heat transfer coefficients, and the temperature of the molten material filling element is calculated over time. Filling heat transfer analysis step;
A solid phase ratio calculating step for calculating, over time, a solid phase ratio of the molten metal material in the molten material filling element according to the calculated temperature of the molten material filling element, ,
A solid phase thickness calculating step of calculating a solid phase thickness which is a thickness of the cavity element having a solid phase ratio of 1 continuously from an inner wall surface of the mold element when filling is completed;
A die casting simulation method comprising: an evaluation step including a surface defect occurrence evaluation step for predicting the occurrence of surface defects from the solid phase thickness (Claim 1).

本発明のダイカストシミュレーション方法では、少なくとも、上記前処理工程と、上記充填伝熱解析工程とを行い、かつ、上記固体相厚さ算出ステップと上記表面欠陥発生評価ステップとを含む上記評価工程とを行う。これによって、従来解析できなかった表面欠陥の発生を予測することが可能となる。   In the die casting simulation method of the present invention, at least the pretreatment step and the filling heat transfer analysis step, and the evaluation step including the solid phase thickness calculation step and the surface defect occurrence evaluation step Do. This makes it possible to predict the occurrence of surface defects that could not be analyzed conventionally.

すなわち、上記課題を解決する目的で本発明者が鋭意研究を行った結果、従来のダイカストシミュレーション方法では、充填途中の凝固状態の解析は行われるが、鋳型に接する表面の固相の厚さの算出は行われず、そのために、最終的に鋳造圧力で鋳物が型に密着することを阻害する固相の影響を充分に把握することができないことを見出した。また、表面欠陥の発生が、この固体相厚さが厚い部分に集中することも見出したのである。固相が厚いと、鋳造圧力で固相を機械的に変形させ鋳型と密着させて急冷することが困難となり、湯じわが発生する。本発明では、この点を評価しうる上記評価工程を積極的に追加したのである。   That is, as a result of intensive research conducted by the present inventors for the purpose of solving the above-mentioned problems, in the conventional die casting simulation method, the solidification state during filling is analyzed, but the thickness of the solid phase on the surface in contact with the mold is It was found that the calculation was not performed, and therefore, it was not possible to fully grasp the influence of the solid phase that ultimately prevented the casting from adhering to the mold at the casting pressure. It has also been found that the occurrence of surface defects concentrates on the portion where the solid phase thickness is thick. If the solid phase is thick, it becomes difficult to rapidly cool the solid phase by mechanically deforming it with the casting pressure and bringing it into close contact with the mold, thereby generating water wrinkles. In this invention, the said evaluation process which can evaluate this point was added positively.

具体的には、上記充填伝熱解析工程の固相率算出ステップによって、充填に伴う各要素の固相率を算出する。そして、上記評価工程の固体相厚さ算出ステップにおいて、充填完了時における、上記型要素の内壁面から連続的に固相率1となる上記キャビティ要素の厚さである固体相厚さを算出する。つまり、充填直後における、キャビティ内の溶融金属材料表面に固化形成された表皮部分の厚さを算出する。そして、この固体相厚さを用いて上記表面欠陥発生評価ステップを行うことによって、従来行えなかった判断が可能となる。   Specifically, the solid phase ratio of each element accompanying the filling is calculated by the solid phase ratio calculating step of the filling heat transfer analysis step. Then, in the solid phase thickness calculation step of the evaluation step, the solid phase thickness, which is the thickness of the cavity element having a solid phase ratio of 1 continuously from the inner wall surface of the mold element when filling is completed, is calculated. . That is, the thickness of the skin portion solidified on the surface of the molten metal material in the cavity immediately after filling is calculated. Then, by performing the surface defect occurrence evaluation step using the solid phase thickness, it becomes possible to make a determination that could not be made conventionally.

それ故、本発明のダイカストシミュレーション方法によれば、溶融材料が充填途中で凝固することに起因して発生する湯じわ状の表面欠陥について、より正確に推測できる。   Therefore, according to the die-casting simulation method of the present invention, it is possible to more accurately estimate the surface defects that are caused by the melting of the molten material during solidification.

第2の発明は、上記第1の発明のダイカストシミュレーション方法をコンピュータに実行させるためのダイカストシミュレーションプログラムにある(請求項5)。
第3の発明は、上記第2の発明のダイカストシミュレーションプログラムを記録したコンピュータ読み取り可能な記憶媒体にある(請求項6)。
第4の発明は、上記第1の発明のダイカストシミュレーション方法を実行するよう構成されたコンピュータを有することを特徴とするダイカストシミュレーション装置にある(請求項7)。
A second invention is a die casting simulation program for causing a computer to execute the die casting simulation method of the first invention.
A third invention resides in a computer-readable storage medium storing the die casting simulation program of the second invention.
According to a fourth aspect of the present invention, there is provided a die casting simulation apparatus comprising a computer configured to execute the die casting simulation method according to the first aspect of the present invention.

これら第2〜第4の発明は、いずれも、これらを用いて上記ダイカストシミュレーションを実施することができ、上述したように、従来解析できなかった表面欠陥の発生を予測することができる。   In any of these second to fourth inventions, the above-described die casting simulation can be performed using them, and as described above, it is possible to predict the occurrence of surface defects that could not be analyzed conventionally.

第1の発明においては、上記のごとく、少なくとも、上記要素作成ステップと上記要素定義ステップとを含む前処理工程と、上記充填伝熱解析ステップと上記固相率算出ステップとを含む充填伝熱解析工程と、上記固体相厚さ算出ステップと上記表面欠陥発生評価ステップとを含む評価工程とを実施する。以下に、これらについてさらに説明する。   In the first invention, as described above, a packed heat transfer analysis including at least a pretreatment process including the element creating step and the element defining step, the filling heat transfer analyzing step, and the solid phase ratio calculating step. A process and an evaluation process including the solid phase thickness calculation step and the surface defect occurrence evaluation step are performed. These will be further described below.

<前処理工程>
前処理工程は(1)要素作成ステップと(2)要素定義ステップとを含み、鋳型のモデルデータを作成して、後述する充填伝熱解析工程を行う準備をする工程である。
(1)要素作成ステップ
要素作成ステップは、本ダイカストシミュレーション方法の対象である鋳型(型)を座標系上に位置づけ、その座標系上の空間を多面体からなる複数の微小要素に分割するステップである。すなわち、座標系上の空間を解析用の微小要素に細分化するステップである。
座標系は、任意のものを選択することが可能である。この座標系上の空間には必要に応じた大きさで微小要素が形成される。
<Pretreatment process>
The pre-processing step includes (1) an element creation step and (2) an element definition step, and is a step of preparing model data for a mold and preparing for a filling heat transfer analysis step described later.
(1) Element creation step The element creation step is a step in which a template (mold) that is the object of the present die casting simulation method is positioned on a coordinate system, and a space on the coordinate system is divided into a plurality of minute elements made of a polyhedron. . That is, this is a step of subdividing the space on the coordinate system into microelements for analysis.
Any coordinate system can be selected. In this space on the coordinate system, microelements are formed with a size as required.

微小要素に分割する方法としては有限差分法で採用されるような直交6面体の微小要素で分割する方法、有限要素法のように要素の形状を鋳造型の形状に応じた多面体として比較的自由に変更できる方法等がある。有限差分法は微小要素への分割が容易であり、かつ解析が数学的に簡潔であるという利点がある。   As a method of dividing into minute elements, a method of dividing with small elements of orthogonal hexahedron as used in the finite difference method, and a relatively polyhedral shape according to the shape of the casting mold as in the finite element method. There are methods that can be changed. The finite difference method is advantageous in that it can be easily divided into small elements and the analysis is mathematically simple.

なお、座標系空間のすべてに微小要素を規定する必要はなく、必要な部分(溶融材料が注入されるキャビティ領域とその周りに接する型領域等のように後述する充填伝熱解析工程で必要な部分)を最小限含むような範囲で規定すれば充分である。但し、鋳型温度等をより正確に解析するために、型領域のすべてを含むように微小要素を作成することが好ましい。そして作成する微小要素の大きさはできるだけ小さい方が解析の精度が向上できるが、より多くの解析時間が必要となる。
また、鋳型の構造を充分に再現できるような微小要素の大きさを採用することが好ましい。従って、微小要素の大きさは要求される精度やシミュレーションの原理的な制約、解析時間等から適正に決定できる。なお、微小要素の大きさはすべての部分について同じ大きさとする必要はなく、解析部位によって大きさを変更することができる。例えば、ダイカスト品の肉薄な部分では、局所的に微小要素の大きさを小さく設定し、解析精度を向上することが好ましい。
In addition, it is not necessary to define minute elements in the entire coordinate system space, and it is necessary in a filling heat transfer analysis process described later such as a necessary part (such as a cavity region into which molten material is injected and a mold region in contact therewith). It is sufficient to define it within a range that includes a minimum of (part). However, in order to analyze the mold temperature and the like more accurately, it is preferable to create a microelement so as to include the entire mold region. The accuracy of analysis can be improved if the size of the minute element to be created is as small as possible, but more analysis time is required.
In addition, it is preferable to employ a size of minute elements that can sufficiently reproduce the structure of the mold. Therefore, the size of the minute element can be appropriately determined from the required accuracy, the fundamental constraint of simulation, the analysis time, and the like. Note that the size of the microelements does not have to be the same for all portions, and the size can be changed depending on the analysis site. For example, in a thin part of a die-cast product, it is preferable to locally set a small element size to improve analysis accuracy.

ところで、型を座標系上に位置づけるためには、鋳型の形状がCADデータ型式等の数値データに変換されている必要がある。鋳型の形状を数値データに変換する方法としては、特に限定されず、例えば、最初から鋳型の形状をCADにより設計したり、試作品の形状を3次元スキャナ等の何らかの方法で数値化しても良い。ここでCADにより型の数値データを作成した場合には、CAD等により作成された型のデータを読み込み、型の外形データを抽出する必要がある。その方法については公知の方法が使用できる。また、本方法においてCADデータをそのまま使用できるようにしても良い。なお、本ステップにおいて、鋳型に代えてダイカストされるダイカスト品(鋳物)を座標系上に位置づけることもできる。   By the way, in order to position the mold on the coordinate system, the shape of the mold needs to be converted into numerical data such as a CAD data model. The method for converting the shape of the mold into numerical data is not particularly limited. For example, the shape of the mold may be designed from the beginning by CAD, or the shape of the prototype may be digitized by some method such as a three-dimensional scanner. . Here, when the numerical value data of the type is created by CAD, it is necessary to read the type data created by CAD or the like and extract the outer shape data of the type. A known method can be used for the method. Further, the CAD data may be used as it is in this method. In this step, a die-cast product (cast) that is die-cast instead of the mold can be positioned on the coordinate system.

(2)要素定義ステップ
要素定義ステップは、前述の要素作成ステップにおいて規定された微小要素のそれぞれについて、鋳型の型領域に位置する場合には型要素と定義し、鋳型のキャビティ領域に位置する場合にはキャビティ要素と定義してモデルデータとするステップである。すなわち、後述の充填伝熱解析工程用に各微小要素の属性を定義し、座標系上に型の形状を微小要素により構築するステップである。
なお、本ステップは、前述の要素作成ステップにおいて微小要素が規定された後に行われるステップであるが、すべての微小要素が規定された後に行う必要は必ずしもなく、1以上の微小要素が規定される毎に本ステップを行い、その後に要素作成ステップを再度行うことを繰り返すこともできる。
(2) Element definition step In the element definition step, each of the microelements defined in the element creation step is defined as a mold element when positioned in the mold area of the mold and positioned in the cavity area of the mold. Is a step of defining model data by defining a cavity element. That is, it is a step of defining the attributes of each microelement for a later-described filling heat transfer analysis step and constructing the shape of the mold on the coordinate system with the microelements.
Note that this step is a step that is performed after the microelements are defined in the above-described element creation step. However, this step is not necessarily performed after all the microelements are defined, and one or more microelements are defined. It is also possible to repeat this step every time and then repeat the element creation step.

ここで、型の「型領域」とは型自身を形成する領域であって、溶融材料が流れない部分であり、型の「キャビティ領域」とは溶融材料が流れ最終的に成形品(鋳物)が形成される部分である領域をそれぞれ意味する。
具体的に各微小要素を型要素とキャビティ要素とに定義する方法としては特に限定されず、公知の方法が採用可能である。以下に図1、図2を参照しながら一例を説明する。
Here, the “mold area” of the mold is an area where the mold itself is formed and the molten material does not flow. The “cavity area” of the mold is the flow of the molten material and finally the molded product (casting). Each means a region that is a portion where is formed.
Specifically, a method for defining each microelement as a mold element and a cavity element is not particularly limited, and a known method can be adopted. An example will be described below with reference to FIGS.

図1には、鋳型の形状及び微小要素の一部を拡大して示す。また、同図は記載及び説明の便宜上2次元上にて鋳型及び微小要素を示し、以下の説明も2次元の図に基づいて行うが、その本質は3次元のものと異なるところはない。
図1に示すように、座標として直交座標を採用し、その座標系上に正方形の微小要素1(形状は特に正方形に限定されるものではない。また、3次元上に適用する場合には直方体・立方体その他任意形状の多面体が要素の形として例示できる。以下同じ。)が規則的に規定されている。また、座標上には、型のモデルデータの境界線Kが位置づけられている。
FIG. 1 shows an enlarged view of the shape of the mold and a part of the microelements. In addition, this figure shows a mold and minute elements on a two-dimensional view for convenience of description and explanation, and the following explanation is also based on a two-dimensional figure, but the essence is not different from a three-dimensional one.
As shown in FIG. 1, orthogonal coordinates are adopted as coordinates, and a square minute element 1 on the coordinate system (the shape is not particularly limited to a square. In addition, when applied to three dimensions, a rectangular parallelepiped.・ Cubes and other polyhedrons of arbitrary shape can be exemplified as element shapes, and so on. Further, the boundary line K of the model data of the mold is positioned on the coordinates.

図1、図2に示すごとく、各微小要素1それぞれの重心2の位置が、型の型領域S1(斜線部分)に存在する場合にはその微小要素1を型要素と定義し、上記境界線Kに囲まれた領域であるキャビティ領域S2に存在する場合にはその微小要素をキャビティ要素と定義する。各微小要素1を型要素及びキャビティ要素に定義した状態を図2に示す。図2では型領域S1に存在する重心2を白丸で、キャビティ領域S2に存在する重心2を黒丸で表す。なお、型領域S1及びキャビティ領域S2のいずれにも該当しない微小要素の扱いは特に限定しないが、計算上の負荷とならないように規定することが好ましい。   As shown in FIGS. 1 and 2, when the position of the center of gravity 2 of each microelement 1 exists in the mold area S1 (shaded portion) of the mold, the microelement 1 is defined as a mold element, and the boundary line If it exists in the cavity region S2, which is a region surrounded by K, the microelement is defined as a cavity element. FIG. 2 shows a state in which each microelement 1 is defined as a mold element and a cavity element. In FIG. 2, the center of gravity 2 existing in the mold region S1 is represented by a white circle, and the center of gravity 2 present in the cavity region S2 is represented by a black circle. In addition, although handling of the microelement which does not correspond to any of mold area | region S1 and cavity area | region S2 is not specifically limited, It is preferable to prescribe | regulate so that it may not become a computational load.

<充填伝熱解析工程>
充填伝熱解析工程は、(1)充填伝熱解析ステップと(2)固相率算出ステップとを行う。
充填伝熱解析工程においては、基準となる時間を設定することが好ましい。つまり、解析の進行と共に、その基準となる時間を微少量ずつ進めていき、その基準となる時間に基づいて各ステップを適用することで、鋳型内の溶融材料の挙動を解析できる。従って、充填伝熱解析工程においては、実時間の進行に関わらず充填伝熱解析工程で設定された任意の時間経過に基づいて解析を行う。そして、各ステップは充填伝熱解析工程内に設定された時間経過に対して、常に同じ頻度で解析を行うことは必須ではない。例えば、これらのステップを行う時間間隔はすべて同じ間隔を採用することもできるし、異なる間隔を採用することもできる。
<Filling heat transfer analysis process>
The filling heat transfer analysis step performs (1) a filling heat transfer analysis step and (2) a solid phase ratio calculation step.
In the filling heat transfer analysis step, it is preferable to set a reference time. That is, as the analysis progresses, the reference time is advanced little by little, and the behavior of the molten material in the mold can be analyzed by applying each step based on the reference time. Therefore, in the filling heat transfer analysis step, the analysis is performed based on an arbitrary time set in the filling heat transfer analysis step regardless of the progress of real time. And it is not indispensable that each step always analyzes at the same frequency with respect to the time passage set in the filling heat transfer analysis process. For example, the time intervals for performing these steps can all be the same interval, or can be different intervals.

(1)充填伝熱解析ステップ
充填伝熱解析ステップでは、充填領域解析を行うと共に、溶融材料充填要素の温度を経時的に算出する。
まず、充填領域解析では、キャビティ要素のそれぞれについて、溶融金属材料の充填状態の解析を経時的に行う。すなわち、鋳型内における注入された溶融金属材料の物理的挙動を、微小要素毎に、かつ微小時間毎に解析する。そして、溶融金属材料が充填された微小要素は溶融材料充填要素として扱う。
基本的な溶湯の充填解析の方法については、特に限定されるものではない。例えば、VOF(Volume of Fluid)、SOLA、FAN及びそれらの改良された方法等の公知技術・慣用技術等を適用することができる。
(1) Filling heat transfer analysis step In the filling heat transfer analysis step, the filling region analysis is performed and the temperature of the molten material filling element is calculated over time.
First, in the filling region analysis, the filling state of the molten metal material is analyzed over time for each of the cavity elements. That is, the physical behavior of the injected molten metal material in the mold is analyzed for every minute element and every minute time. And the microelement with which the molten metal material was filled is handled as a molten material filling element.
There is no particular limitation on the basic method for analysis of molten metal filling. For example, known techniques and conventional techniques such as VOF (Volume of Fluid), SOLA, FAN, and improved methods thereof can be applied.

また、上記溶融材料充填要素について、それぞれの溶融材料充填要素同士の間、及び該溶融材料充填要素と上記型要素との間の伝熱をそれぞれの熱伝達係数に基づいて経時的に解析して、それぞれの上記溶融材料充填要素の温度を経時的に算出する。
伝熱の解析は、溶融材料充填要素同士の間でのほか、溶融材料充填要素と鋳造型最表面(すなわちキャビティ要素に接する型要素)との間で経時的に行い、各微小要素の温度を算出する。
Further, for the molten material filling element, the heat transfer between the molten material filling elements and between the molten material filling element and the mold element is analyzed over time based on the respective heat transfer coefficients. The temperature of each of the molten material filling elements is calculated over time.
The heat transfer analysis is performed not only between the molten material filling elements but also over time between the molten material filling element and the outermost surface of the casting mold (that is, the mold element in contact with the cavity element). calculate.

また、充填伝熱解析ステップでは、計算が発散せずかつ許容時間内で計算が終了するように設定された時間間隔で、各要素間の伝熱を各モデルに設定された伝熱係数に基づいて計算する。充填伝熱解析ステップで行う伝熱解析方法は特に限定されるものではないが、例えば、熱移流、潜熱を考慮した非定常熱伝導解析に差分法とADI法とを併用する等の計算方法を用いてそれぞれの要素について熱の伝導を計算することができる。   In the heat transfer analysis step, the heat transfer between each element is based on the heat transfer coefficient set in each model at a time interval set so that the calculation does not diverge and the calculation is completed within the allowable time. To calculate. The heat transfer analysis method performed in the filling heat transfer analysis step is not particularly limited. For example, a calculation method such as using the differential method and the ADI method in combination with unsteady heat conduction analysis in consideration of heat advection and latent heat. Can be used to calculate the heat conduction for each element.

(2)固相率算出ステップ
固相率算出ステップは、充填伝熱解析ステップで算出された溶融材料充填要素の温度に基づいてその溶融材料充填要素に充填された溶融材料の固相率を算出する。固相率は、その溶融材料充填要素内の溶融金属材料のすべてが固化した状態を1、すべてが溶融状態にある場合を0として、その中間の状態は0〜1の間の数値で表す。固相率の算出は状態図等により求めたり、シャイルの式等の理論式又は近似式により算出することが可能である。ここで、上記固相率は、固相率そのものだけでなく、これに関連のあるパラメータを含む概念であり、固相率と関連のあるパラメータである液相率や温度等を用いて計算を行っても良い。液相率は1から固相率を減じた値である。固相率に関連するパラメータとして温度を用いる場合には、すべて液相になる温度以上及びすべて固相になる温度以下はそれぞれ同一に扱う。
(2) Solid phase ratio calculation step The solid phase ratio calculation step calculates the solid phase ratio of the molten material filled in the molten material filling element based on the temperature of the molten material filling element calculated in the filling heat transfer analysis step. To do. The solid phase ratio is represented by a numerical value between 0 and 1 when the state in which all of the molten metal material in the molten material filling element is solidified is 1, and when the state is all in the molten state, and the intermediate state is 0. The solid phase ratio can be calculated from a phase diagram or the like, or can be calculated by a theoretical formula or approximate formula such as a Seil formula. Here, the solid phase ratio is a concept that includes not only the solid phase ratio itself but also parameters related to the solid phase ratio, and is calculated using the liquid phase ratio, temperature, etc., which are parameters related to the solid phase ratio. You can go. The liquid phase ratio is a value obtained by subtracting the solid phase ratio from 1. When temperature is used as a parameter related to the solid phase ratio, the temperature above the temperature at which all become liquid phase and the temperature below at which all becomes solid phase are treated in the same manner.

<評価工程>
評価工程では、(1)固体相厚さ算出ステップと、(2)表面欠陥発生評価ステップとを行う。
(1)固体相厚さ算出ステップ
固体相厚さ算出ステップは、充填完了時における、上記型要素の内壁面から連続的に固相率1となる上記キャビティ要素の厚さである固体相厚さを算出するステップである。固体相厚さは、上述したごとく、充填直後におけるキャビティ内の溶融金属材料表面に固化形成された表皮部分の厚さを示すものである。このステップは、充填完了直後に行う。
<Evaluation process>
In the evaluation process, (1) a solid phase thickness calculation step and (2) a surface defect occurrence evaluation step are performed.
(1) Solid phase thickness calculation step The solid phase thickness calculation step is a solid phase thickness that is the thickness of the cavity element that has a solid phase ratio of 1 continuously from the inner wall surface of the mold element when filling is completed. Is a step of calculating. As described above, the solid phase thickness indicates the thickness of the skin portion solidified on the surface of the molten metal material in the cavity immediately after filling. This step is performed immediately after completion of filling.

通常、溶融金属材料は、鋳型に接する壁面から熱を奪われ凝固するために、鋳型の内壁面(キャビティ面)に隣接するキャビティ要素から固相率は高くなる。固相率が1となる要素は鋳型の内壁面側に位置する。そして、固相率1となる要素の内側には固相率1未満となる要素が存在することとなる。
鋳型(型要素)の内壁面から固相率<1となる要素までの距離を求めることにより、固相率1となる固体相厚さを求めることができる。そして、内壁面から連続的に固相率1となる場合の最終要素に隣接する固相率<1となる要素の固体相厚さは、要素体積と固相率、隣接部の面積から求めることができる。
そして、前記の連続的に固相率1となる要素の厚さの和と、その最終要素に隣接する固相率<1となる要素内の固体相厚さの和が求める固体相厚さとなる。また、内壁面と隣接する要素が固相率<1のとき、固体相厚さはその要素体積と固相率、隣接部の面積により求まる値となる。
Usually, since the molten metal material is deprived of heat from the wall surface in contact with the mold and solidified, the solid phase ratio is increased from the cavity element adjacent to the inner wall surface (cavity surface) of the mold. The element having a solid phase ratio of 1 is located on the inner wall surface side of the mold. An element having a solid phase ratio of less than 1 exists inside the element having a solid phase ratio of 1.
By obtaining the distance from the inner wall surface of the mold (mold element) to the element having the solid phase ratio <1, the solid phase thickness at which the solid phase ratio is 1 can be obtained. The solid phase thickness of the element having a solid phase ratio <1 adjacent to the final element when the solid phase ratio is continuously 1 from the inner wall surface is obtained from the element volume, the solid phase ratio, and the area of the adjacent portion. Can do.
The sum of the thicknesses of the elements having a solid phase ratio of 1 and the sum of the solid phase thicknesses in the elements having a solid phase ratio of <1 adjacent to the final element is obtained as the solid phase thickness. . When the element adjacent to the inner wall surface has a solid phase ratio <1, the solid phase thickness is a value determined by the element volume, the solid phase ratio, and the area of the adjacent part.

上記固体相厚さを求める算出式としては、例えば次のものがある。   Examples of the calculation formula for obtaining the solid phase thickness include the following.

Figure 2006247684
Figure 2006247684

ここで、L:固体相厚さ、n:内壁面から連続的に固相率=1となる要素の数、Vn+1:壁面からn+1番目の要素体積、fsn+1:壁面からn+1番目の要素の固相率、Sn+1:壁面からn+1番目の要素のn番目要素に接する部位の面積、Lm:壁面からm番目の要素の長さ。 Where L: solid phase thickness, n: number of elements with solid phase ratio = 1 continuously from the inner wall surface, V n + 1 : n + 1-th element volume from wall surface, fs n + 1 : n + 1 from wall surface The solid phase ratio of the n th element, S n + 1 : the area of the n + 1 th element from the wall surface in contact with the n th element, L m : the length of the m th element from the wall surface.

(2)表面欠陥発生評価ステップ
表面欠陥発生評価ステップは、求めた固体相厚さから表面欠陥の発生を予測するステップである。
具体的には、求めた固体相厚さが予め定めた基準厚さよりも厚いか否かを、各部分において行い、厚い部分に表面欠陥が発生すると予測することができる。また、この基準を定めるに当たっては、鋳造圧力等の実際の鋳造条件によって基準値を変えることにより、より精度の高い評価を行うことができる。
(2) Surface defect occurrence evaluation step The surface defect occurrence evaluation step is a step of predicting the occurrence of surface defects from the obtained solid phase thickness.
Specifically, whether or not the obtained solid phase thickness is thicker than a predetermined reference thickness is performed in each part, and it can be predicted that surface defects will occur in the thick part. Further, in determining this standard, more accurate evaluation can be performed by changing the standard value according to actual casting conditions such as casting pressure.

この表面欠陥発生評価ステップは、上記固体相厚さが0.4mm以上の場所に表面欠陥が発生すると予測することが好ましい(請求項2)。後述する実施例に示すように、実験の結果では、少なくとも本発明のダイカストシミュレーション方法において0.4mm以上の場所に表面欠陥が発生すると予測することにより高い精度を得ることができる。   In this surface defect occurrence evaluation step, it is preferable to predict that a surface defect will occur at a location where the solid phase thickness is 0.4 mm or more. As shown in the examples described later, it is possible to obtain high accuracy by predicting that a surface defect is generated at a location of 0.4 mm or more in at least the die casting simulation method of the present invention.

また、上記表面欠陥発生評価ステップは、上記固体相厚さが0.1×e0.0297×Pmm(ここで、Pは鋳造圧力(MPa))以上の場所に表面欠陥が発生すると予測することがより好ましい(請求項3)。すなわち、通常の場合には、上記のごとく、1つの固定した基準値で評価可能であるが、例えば、鋳造圧力の条件で実際の表面欠陥の発生状態が変わる場合には、評価基準にその想定している鋳造圧力を反映させることが好ましい。その一例が上記の算出式である。このように、鋳造圧力によって異なる基準を設定することにより、より実際に近い解析を行うことが可能となる。なお、上記固体相厚さ算出ステップまでの工程では、解析に鋳造圧力を用いる必要は特にない。 Further, the surface defect generation evaluation step may predict that a surface defect is generated at a location where the solid phase thickness is 0.1 × e 0.0297 × P mm (where P is a casting pressure (MPa)) or more. More preferred (Claim 3). That is, in the normal case, as described above, it is possible to evaluate with one fixed reference value. However, for example, when the actual surface defect occurrence state changes under casting pressure conditions, the evaluation standard assumes that It is preferable to reflect the casting pressure currently being performed. One example is the above calculation formula. In this way, by setting different standards depending on the casting pressure, it is possible to perform an analysis that is closer to actuality. In the process up to the solid phase thickness calculation step, it is not particularly necessary to use casting pressure for analysis.

また、本発明のダイカストシミュレーション方法は、解析対象となる上記溶融金属材料は、凝固時に表皮側から凝固相が形成される表皮生成型凝固を示す合金であることが最適である(請求項3)。ダイカスト用の合金の中には、むしろ上記のような表皮生成型凝固を示す合金よりも、内部から凝固しやすいタイプの方が多く、その場合には、表面欠陥が起こりにくい。一方、上記表皮生成型凝固を示す合金は、ダイカストにおいて表面欠陥を発生しやすいので、本発明の方法の利用が特に有効である。   In the die casting simulation method of the present invention, it is optimal that the molten metal material to be analyzed is an alloy exhibiting skin-generating solidification in which a solidified phase is formed from the skin side during solidification (Claim 3). . Among die casting alloys, there are more types that are more likely to solidify from the inside than alloys that exhibit skin-generating solidification as described above, in which case surface defects are less likely to occur. On the other hand, the alloy showing the skin formation type solidification is likely to cause surface defects in die casting, and therefore, the use of the method of the present invention is particularly effective.

具体的には、上記溶融金属材料は、Si含有量が8〜14mass%であるAl−Si系アルミニウム合金であることが好ましい(請求項4)。上記特定のSi含有量を有するAl−Si系アルミニウム合金は、特に、ダイカストによって表面欠陥の発生率が高い。そのため、上記本発明のダイカストシミュレーション方法の活用が非常に有効である。   Specifically, the molten metal material is preferably an Al—Si based aluminum alloy having a Si content of 8 to 14 mass%. The Al—Si-based aluminum alloy having the specific Si content particularly has a high incidence of surface defects due to die casting. Therefore, the use of the above-described die casting simulation method of the present invention is very effective.

(実施例1)
本発明の実施例に係るダイカストシミュレーション方法及びその装置につき、図3〜図5を用いて説明する。
本例のダイカストシミュレーション装置1は、図4に示すごとく、各解析ステップを実行する手段を備えたコンピュータ10を有する装置である。コンピュータ10には、キーボード、記憶媒体読み取り装置等の入力装置11と、ディスプレイ、プリンタ等の出力装置12が接続されている。
Example 1
A die casting simulation method and apparatus according to an embodiment of the present invention will be described with reference to FIGS.
The die casting simulation apparatus 1 of this example is an apparatus having a computer 10 provided with means for executing each analysis step as shown in FIG. An input device 11 such as a keyboard and a storage medium reading device and an output device 12 such as a display and a printer are connected to the computer 10.

コンピュータ10が実行可能な手段は、後述するダイカストシミュレーション方法における各ステップを実行する手段であり、具体的には、要素作成ステップ実行手段21と、要素定義ステップ実行手段22と、充填伝熱解析ステップ実行手段23と、固相率算出ステップ実行手段24と、固体相厚さ算出ステップ実行手段25と、表面欠陥発生評価ステップ実行手段26とを有している。なお、これらの実行手段は、後述するダイカストシミュレーション方法をコンピュータに実行させるためのダイカストシミュレーションプログラムとして作成したものを読み込んで実行可能としたものである。   The means that can be executed by the computer 10 is means for executing each step in a die casting simulation method to be described later. Specifically, the element creation step execution means 21, the element definition step execution means 22, and the filling heat transfer analysis step. An execution unit 23, a solid phase ratio calculation step execution unit 24, a solid phase thickness calculation step execution unit 25, and a surface defect occurrence evaluation step execution unit 26 are included. These execution means read and execute what is created as a die casting simulation program for causing a computer to execute a die casting simulation method to be described later.

本例のダイカストシミュレーション方法は、図3に示すごとく、要素作成ステップS101と要素定義ステップS102とを含む前処理工程S100と、充填伝熱解析ステップS201と固相率算出ステップS202とを含む充填伝熱解析工程S200と、固体相厚さ算出ステップS301と表面欠陥発生評価ステップS302とを含む評価工程S300を行う方法である。
以下、これをさらに詳説する。
As shown in FIG. 3, the die casting simulation method of the present example includes a pretreatment process S100 including an element creation step S101 and an element definition step S102, a filling heat transfer analysis step S201, and a solid phase ratio calculation step S202. This is a method of performing an evaluation step S300 including a thermal analysis step S200, a solid phase thickness calculation step S301, and a surface defect occurrence evaluation step S302.
This will be described in further detail below.

本例では、鋳型のモデルデータをCADで作成し、そのモデルデータを用いてダイカストシミュレーションを行う。
(1)前処理工程S100
座標系として直交座標系を採用した。鋳造型の形状はCADデータとして作成される。説明を簡略化するためにx、yで表す2次元での説明を図5に示す。以下に行う2次元での説明は単純にそのまま3次元に拡張可能である。
まず、CADデータを2次元の座標系上に配置する。そして、座標系をx、y軸方向で微小要素に分割する(要素作成ステップS101)。
In this example, model data of a mold is created by CAD, and die casting simulation is performed using the model data.
(1) Pretreatment process S100
An orthogonal coordinate system was adopted as the coordinate system. The shape of the casting mold is created as CAD data. In order to simplify the description, a two-dimensional description represented by x and y is shown in FIG. The following description in two dimensions can be simply extended to three dimensions.
First, CAD data is arranged on a two-dimensional coordinate system. Then, the coordinate system is divided into minute elements in the x and y axis directions (element creation step S101).

次に、この微小要素1の重心位置がCADデータの鋳造型領域S1内に位置する要素を型要素1a、キャビティ領域S2内に位置する要素をキャビティ要素1bと定義する(要素定義ステップS102)。また、キャビティ要素1b間、型要素1a間及びキャビティ要素1bと型要素1aとの間のそれぞれの熱伝達係数は予め適正値を設定する。   Next, an element in which the position of the center of gravity of the minute element 1 is located in the casting mold area S1 of the CAD data is defined as a mold element 1a, and an element located in the cavity area S2 is defined as a cavity element 1b (element definition step S102). Further, the heat transfer coefficients between the cavity elements 1b, between the mold elements 1a, and between the cavity element 1b and the mold element 1a are set to appropriate values in advance.

なお、図3に示すごとく、本例の鋳造型においては、図面右側に溶湯が注入されるゲートGが、左側にキャビティ領域S2内の空気を排出するためのベントホールVがそれぞれ配置されている。また、上述したごとく、四角はそれぞれ微小要素1を示し、微小要素1の境界線のうちの太線Kは鋳型のCADデータを表す。   As shown in FIG. 3, in the casting mold of this example, the gate G into which the molten metal is injected is arranged on the right side of the drawing, and the vent hole V for discharging the air in the cavity region S2 is arranged on the left side. . As described above, each square represents the microelement 1, and the bold line K of the boundary line of the microelement 1 represents the CAD data of the template.

(2)充填伝熱解析工程S200
充填伝熱解析工程S200では、充填領域解析を行うと共に溶融材料充填要素の温度を経時的に算出する充填伝熱解析ステップS201と、溶融材料充填要素内の上記溶融金属材料の固相率を経時的に算出する固相率算出ステップS202とを行う。各解析は、シミュレーション内の時間の進行に従って、それぞれに設定された微小時間間隔毎に行う。
(2) Filling heat transfer analysis step S200
In the filling heat transfer analysis step S200, a filling heat transfer analysis step S201 for performing the filling region analysis and calculating the temperature of the molten material filling element over time, and the solid phase ratio of the molten metal material in the molten material filling element over time are performed. The solid phase ratio calculating step S202 is calculated. Each analysis is performed at each minute time interval set according to the progress of time in the simulation.

充填領域解析では、キャビティ要素について溶融金属の充填割合を所定時間間隔で順次計算していく。充填領域解析は、有限差分法によって解析を進行させ、射出開始から所定時間後に、いくつかのキャビティ要素に溶融金属材料が充填され、更に時間が進行すると、溶融金属材料が充填されたキャビティ要素が増加していく。そして、これと共に行う伝熱の解析では、非定常熱伝導計算法により溶湯が充填されたキャビティ各微小要素(溶融材料充填要素)の温度を解析する(充填伝熱解析ステップS201)。   In the filling region analysis, the filling ratio of the molten metal for the cavity element is sequentially calculated at predetermined time intervals. The filling region analysis proceeds by the finite difference method. After a predetermined time from the start of injection, several cavity elements are filled with the molten metal material, and when the time further advances, the cavity element filled with the molten metal material is changed. It will increase. In the heat transfer analysis performed together with this, the temperature of each microelement (molten material filling element) filled with the molten metal is analyzed by an unsteady heat conduction calculation method (filling heat transfer analysis step S201).

また、固相率算出ステップS202では、算出された温度に基づき、上記各溶融材料充填要素について、固相率を算出する。固相率は、各溶融材料充填要素の温度をシャイルの式にあてはめて算出する。この固相率算出ステップS202は、所定時間間隔として上記充填伝熱解析ステップS201と同じ微小時間間隔毎に行う。この充填伝熱解析ステップS201と固相率算出ステップS202は、充填完了するまで繰り返し行う。そして、ステップS209において充填完了したことが確認された後、評価工程S300に進む。   In the solid phase ratio calculating step S202, the solid phase ratio is calculated for each of the molten material filling elements based on the calculated temperature. The solid phase ratio is calculated by applying the temperature of each molten material filling element to the Seil's equation. This solid phase ratio calculation step S202 is performed at the same minute time interval as the filling heat transfer analysis step S201 as a predetermined time interval. The filling heat transfer analysis step S201 and the solid phase ratio calculation step S202 are repeated until filling is completed. Then, after confirming that the filling is completed in step S209, the process proceeds to the evaluation step S300.

(3)評価工程S300
評価工程S300では、まず固体相厚さ算出ステップS301を行う。このステップは、上記の充填伝熱解析工程S200の結果を踏まえ、キャビティへの溶融金属の充填が完了した直後に行う。
(3) Evaluation process S300
In the evaluation step S300, first, a solid phase thickness calculation step S301 is performed. This step is performed immediately after the filling of the molten metal into the cavity is completed based on the result of the above-described filling heat transfer analysis step S200.

固体相厚さ算出ステップS301では、まず、型要素と隣接するキャビティ要素の中で、固相率が0を越えた要素のある領域を検出する。
次いで、検出された要素の固相率fsが、0<fs<1の場合には、その要素をCs要素とし、その要素内における固体相厚さLcのみが、その部分での固体相厚さとなる。このときの算出式としては、固体相厚さLc=(要素体積×固相率)/(型要素と接した面の面積)を採用することができる。
In the solid phase thickness calculation step S301, first, a region where an element having a solid phase ratio exceeding 0 is detected in the cavity element adjacent to the mold element.
Next, when the solid phase ratio fs of the detected element is 0 <fs <1, the element is set as a Cs element, and only the solid phase thickness Lc in the element is equal to the solid phase thickness in the portion. Become. As a calculation formula at this time, the solid phase thickness Lc = (element volume × solid phase ratio) / (area of the surface in contact with the mold element) can be adopted.

一方、上記の型要素と隣接するキャビティ要素のうち固相率が0を越えたと検出された領域に存在する要素のうち、固相率が1のものは、準型要素として認定し、その要素の固体相厚さLjは、Lj=(要素体積)/(型要素と接した面の面積)という算出式で算出する。
次いで、上記準型要素に対して厚み方向に隣接するキャビティ要素の固相率が1の場合には、この要素も準型要素として、その固体相厚さLjは、Lj=(要素体積)/(上記の準型要素と接した面の面積)という算出式で算出する。以下、厚み方向に接する要素の固相率が1である限り、同様の算出式で各要素の厚さLjを算出する。
また、上記の準型要素に対して厚み方向に隣接するキャビティ要素が上述したCs要素である場合には、その固体相厚さLcを、Lc=(要素体積×固相率)/(上記の準型要素と接した面の面積)という算出式で算出する。
On the other hand, among the elements in the cavity element adjacent to the mold element, those having a solid phase ratio of 1 in the region where the solid phase ratio is detected to exceed 0 are recognized as quasi-type elements, and the elements The solid phase thickness Lj is calculated by the following formula: Lj = (element volume) / (area of the surface in contact with the mold element).
Next, when the solid phase ratio of the cavity element adjacent to the quasi-type element in the thickness direction is 1, this element is also a quasi-type element, and the solid phase thickness Lj is Lj = (element volume) / (Area of the surface in contact with the above-mentioned quasi-type element) Hereinafter, as long as the solid phase ratio of the element in contact with the thickness direction is 1, the thickness Lj of each element is calculated with the same calculation formula.
When the cavity element adjacent in the thickness direction to the quasi-type element is the Cs element described above, the solid phase thickness Lc is expressed as Lc = (element volume × solid phase ratio) / (the above The area of the surface in contact with the quasi-type element is calculated using a calculation formula.

そして、型要素と厚み方向に順次接する要素のうち、すべての準型要素の固体相厚さLjの総和(ΣLj)と、Cs要素の固体相厚さLcとの和、つまりΣLj+Lcが、求める固体相厚さLとなる。
なお、固相率が0の要素の固体相厚さLは当然0である。
The sum of the solid phase thickness Lj of all quasi-type elements (ΣLj) and the solid phase thickness Lc of the Cs element, that is, ΣLj + Lc, among the elements that are sequentially in contact with the mold element in the thickness direction, that is, ΣLj + Lc The phase thickness is L.
The solid phase thickness L of the element having a solid phase ratio of 0 is naturally 0.

次に、表面欠陥発生評価ステップS302を行う。このステップでは、型要素に隣接するキャビティ要素毎に求められた固体相厚さLから評価を行う。
例えば、固体相厚さが0.4mmを越えた場合に、表面欠陥発生要素とする。なお、こ0.4mmという基準は、実際に行われる鋳造条件や溶融金属材料の材質等に応じて変更することも可能である。
いずれにしても、型要素に隣接するすべてのキャビティ要素について、その厚み方向の合計の固体相厚さLが基準値を超える場合には、そのキャビティ要素を特定して他と区別することにより、表面欠陥が発生する分布を、コンピュター上に表現することができる。具体的には、得られる鋳物(成形品)の画像に対し、表面欠陥が発生する部分を所定のマークで塗りつぶす等の手法で描画的に出力するすることができる。また、画像ではなく、座標データのみで出力することも可能である。
Next, a surface defect occurrence evaluation step S302 is performed. In this step, evaluation is performed from the solid phase thickness L obtained for each cavity element adjacent to the mold element.
For example, when the solid phase thickness exceeds 0.4 mm, the surface defect generating element is used. The standard of 0.4 mm can be changed according to the actual casting conditions and the material of the molten metal material.
In any case, for all the cavity elements adjacent to the mold element, if the total solid phase thickness L in the thickness direction exceeds the reference value, the cavity element is identified and distinguished from the others, The distribution in which surface defects occur can be expressed on a computer. Specifically, it is possible to output the image of a casting (molded product) to be drawn by a technique such as painting a portion where a surface defect occurs with a predetermined mark. It is also possible to output only coordinate data, not an image.

(実施例2)
本例は、実施例1の方法を用いてシミュレーションを行った結果についての具体例を示す。
まず、本例では、上記シミュレーションの有効性を評価すべく、実際にダイカスト用鋳型を作製し、図6に示す形状の鋳物8を実際に作製し、表面欠陥を観察した。鋳物8は平板状であって、これには、溶湯の注入加圧部81及びランナ部82、ゲート部83がつながっている。
(Example 2)
This example shows a specific example of the result of simulation using the method of the first embodiment.
First, in this example, in order to evaluate the effectiveness of the simulation, a die casting mold was actually manufactured, a casting 8 having the shape shown in FIG. 6 was actually manufactured, and surface defects were observed. The casting 8 has a flat plate shape, to which a molten metal injection pressurizing portion 81, a runner portion 82, and a gate portion 83 are connected.

鋳造圧力としては、30MPaと65MPaの2種類の条件を用いた。材質は、ADC12合金とした。また、その他の主な条件としては、射出温度:640℃、低速速度:0.3m/s、高速速度:1.5m/sとした。上記の低速速度及び高速速度は、溶融金属材料を射出するのに用いる加圧用プランジャ(図示略)の前進スピードを意味する。本例では、その速度を、最初は上記低速速度の条件とし、途中で上記高速速度の条件に切り替える二段射出方法を採用した。   As casting pressure, two kinds of conditions of 30 MPa and 65 MPa were used. The material was ADC12 alloy. Other main conditions were injection temperature: 640 ° C., low speed: 0.3 m / s, and high speed: 1.5 m / s. The low speed and the high speed mean the forward speed of a pressurizing plunger (not shown) used for injecting the molten metal material. In this example, a two-stage injection method is adopted in which the speed is initially set to the condition of the low speed and then switched to the condition of the high speed in the middle.

そして、得られた鋳物8の表面を観察して、表面欠陥(湯じわ)発生部9を特定して、その結果を図7(a)(b)に示した。同図の(a)は、鋳造圧力30MPaにより製造したものであり、同図(b)は、鋳造圧力65MPaにより製造したものである。また、これらの図には、観察された表面欠陥(湯じわ)発生部9を示す。   And the surface of the obtained casting 8 was observed, the surface defect (hot water wrinkle) generating part 9 was specified, and the result was shown to Fig.7 (a) (b). (A) of the figure is manufactured with a casting pressure of 30 MPa, and (b) of the figure is manufactured with a casting pressure of 65 MPa. Moreover, these figures show the observed surface defect (water bath) generating portion 9.

次に、実施例1に示したダイカストシミュレーション方法を実施した。本例では、図6に示す平板状の鋳物8を作ることを想定し、また、溶湯の注入加圧部81及びランナ部82、ゲート部83も鋳物の一部としたモデルを想定する。そして、これらを囲う部分を鋳型部分とし、その形状をCADデータとして作成し、上述した前処理工程S100を行った。   Next, the die casting simulation method shown in Example 1 was performed. In this example, it is assumed that the flat casting 8 shown in FIG. 6 is made, and a model in which the molten metal injection pressurizing portion 81, the runner portion 82, and the gate portion 83 are also part of the casting is assumed. And the part which encloses these was made into the casting_mold | template part, the shape was created as CAD data, and pre-processing process S100 mentioned above was performed.

また、上記充填伝熱解析工程S200を行った後の評価工程S300の表面欠陥発生評価ステップS302においては、まずは、基準値を0.4mmとして、固体相厚さLが0.4mm以上の場所に表面欠陥が発生すると予測することとした。この場合、シミュレーションでは、鋳造圧力については全く考慮していない。
その結果を図8に示す。同図に示すごとく、シミュレーション結果の出力は、鋳物8の画像の表面に、表面欠陥が発生しうる部分5を示した。
図8と前述した図7との比較から知られるように、今回のシミュレーション結果は、図7(b)に示す鋳造圧力65MPaの場合とよく一致していることが分かる。一方、図7(a)の鋳造圧力30MPaの場合と比べると、若干表面欠陥が発生しうる部分5が少ない結果となった。
In addition, in the surface defect occurrence evaluation step S302 of the evaluation step S300 after the filling heat transfer analysis step S200 is performed, first, the reference value is set to 0.4 mm, and the solid phase thickness L is set to 0.4 mm or more. It was predicted that surface defects would occur. In this case, the simulation does not consider the casting pressure at all.
The result is shown in FIG. As shown in the figure, the simulation result output shows a portion 5 where surface defects may occur on the surface of the image of the casting 8.
As can be seen from the comparison between FIG. 8 and FIG. 7 described above, it can be seen that the simulation results this time are in good agreement with the casting pressure of 65 MPa shown in FIG. On the other hand, as compared with the case of the casting pressure of 30 MPa in FIG.

次に、本例では、再度ダイカストシミュレーションを実施した。評価工程S300の表面欠陥発生評価ステップS302においては、その判断の基準値を、0.1×e0.0297×Pmmとした。そして、鋳造圧力Pの値としては、30MPaを用いた。
そのシミュレーション結果を図9に示す。同図は、図8の場合と同様に、鋳物8の画像の表面に表面欠陥が発生しうる部分5を示したものである。
図9と前述した図7との比較から知られるように、今回のシミュレーション結果は、図7(a)に示す鋳造圧力30MPaの場合とよく一致していることが分かる。
Next, in this example, the die casting simulation was performed again. In the surface defect occurrence evaluation step S302 of the evaluation step S300, the reference value for the determination was set to 0.1 × e 0.0297 × P mm. And as a value of casting pressure P, 30 MPa was used.
The simulation result is shown in FIG. This figure shows a portion 5 where surface defects may occur on the surface of the image of the casting 8 as in the case of FIG.
As can be seen from the comparison between FIG. 9 and FIG. 7 described above, it can be seen that the simulation result of this time is in good agreement with the case of the casting pressure of 30 MPa shown in FIG.

本発明における、要素作成方法の一例を示す説明図。Explanatory drawing which shows an example of the element creation method in this invention. 本発明における、要素定義方法の一例示す説明図。Explanatory drawing which shows an example of the element definition method in this invention. 実施例1における、ダイカストシミュレーション方法を示すフロー図。FIG. 3 is a flowchart showing a die casting simulation method in the first embodiment. 実施例1における、ダイカストシミュレーション装置の構成を示す説明図。BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is an explanatory diagram illustrating a configuration of a die casting simulation apparatus in Embodiment 1. 実施例1における、要素作成結果を示す説明図。Explanatory drawing which shows the element creation result in Example 1. FIG. 実施例2における、実際に作製した鋳物を示す斜視図。The perspective view which shows the casting actually produced in Example 2. FIG. 実施例2における、(a)鋳造圧力30MPaの条件で作製した鋳物における表面欠陥の発生分布を示す説明図、(b)鋳造圧力65MPaの条件で作製した鋳物における表面欠陥の発生分布を示す説明図。In Example 2, (a) Explanatory view showing the generation distribution of surface defects in a casting produced under a casting pressure of 30 MPa, (b) Explanatory view showing the occurrence distribution of surface defects in a casting produced under a casting pressure of 65 MPa. . 実施例2における、第1のシミュレーション結果を示す説明図。FIG. 9 is an explanatory diagram showing a first simulation result in Example 2. 実施例2における、第2のシミュレーション結果を示す説明図。Explanatory drawing which shows the 2nd simulation result in Example 2. FIG.

符号の説明Explanation of symbols

1 ダイカストシミュレーション装置
10 コンピュータ
21 要素作成ステップ実行手段
22 要素定義ステップ実行手段
23 充填伝熱解析ステップ実行手段
24 固相率算出ステップ実行手段
25 固相厚さ算出ステップと実行手段
26 表面欠陥発生評価ステップ実行手段
S100 前処理工程
S101 要素作成ステップ
S102 要素定義ステップ
S200 充填伝熱解析工程
S201 充填伝熱解析ステップ
S202 固相率算出ステップ
S300 評価工程
S301 固体厚さ算出ステップ
S302 表面欠陥発生評価ステップ
DESCRIPTION OF SYMBOLS 1 Die-casting simulation apparatus 10 Computer 21 Element creation step execution means 22 Element definition step execution means 23 Filling heat transfer analysis step execution means 24 Solid phase ratio calculation step execution means 25 Solid phase thickness calculation step and execution means 26 Surface defect generation evaluation step Execution means S100 Preprocessing step S101 Element creation step S102 Element definition step S200 Filling heat transfer analysis step S201 Filling heat transfer analysis step S202 Solid phase ratio calculation step S300 Evaluation step S301 Solid thickness calculation step S302 Surface defect occurrence evaluation step

Claims (8)

鋳型内に溶融金属材料を加圧注入して凝固させることにより所望形状の鋳物を得るダイカストを行う際に、少なくとも上記鋳物に生じる表面欠陥の発生を予測するダイカストシミュレーション方法であって、
上記溶融金属材料の成形に用いる上記鋳型の形状を座標系上に位置づけ、該座標系の空間を複数の微小要素に分割する要素作成ステップと、
上記微小要素のそれぞれについて、上記鋳型の型領域に位置する場合には型要素と、上記鋳型のキャビティ領域に位置する場合にはキャビティ要素と定義する要素定義ステップと、を含む前処理工程と、
上記キャビティ要素のそれぞれについて、上記溶融金属材料の充填状態を経時的に解析する充填領域解析を行うと共に、上記溶融金属材料が充填された上記キャビティ要素である溶融材料充填要素について、該溶融材料充填要素同士の間、及び該溶融材料充填要素と上記型要素との間の伝熱をそれぞれの熱伝達係数に基づいて経時的に解析して、上記溶融材料充填要素の温度を経時的に算出する充填伝熱解析ステップと、
算出された上記溶融材料充填要素の温度に応じて、該溶融材料充填要素内の上記溶融金属材料の固相率を経時的に算出する固相率算出ステップと、を含む充填伝熱解析工程と、
充填完了時における、上記型要素の内壁面から連続的に固相率1となる上記キャビティ要素の厚さである固体相厚さを算出する固体相厚さ算出ステップと、
該固体相厚さから表面欠陥の発生を予測する表面欠陥発生評価ステップと、を含む評価工程とを有することを特徴とするダイカストシミュレーション方法。
A die casting simulation method for predicting at least the occurrence of surface defects occurring in the casting when performing die casting to obtain a casting having a desired shape by pressurizing and solidifying a molten metal material into a mold,
An element creating step of positioning the shape of the mold used for forming the molten metal material on a coordinate system and dividing the space of the coordinate system into a plurality of microelements;
For each of the microelements, a preprocessing step including a mold element when positioned in the mold region of the mold and an element defining step that defines a cavity element when positioned in the cavity region of the mold;
For each of the cavity elements, a filling region analysis for analyzing the filling state of the molten metal material over time is performed, and the molten material filling element that is the cavity element filled with the molten metal material is filled with the molten material. The heat transfer between the elements and between the molten material filling element and the mold element is analyzed over time based on the respective heat transfer coefficients, and the temperature of the molten material filling element is calculated over time. Filling heat transfer analysis step;
A solid phase ratio calculating step for calculating, over time, a solid phase ratio of the molten metal material in the molten material filling element according to the calculated temperature of the molten material filling element, ,
A solid phase thickness calculating step of calculating a solid phase thickness that is a thickness of the cavity element that is continuously at a solid phase ratio of 1 from an inner wall surface of the mold element when filling is completed;
A die casting simulation method comprising: an evaluation step including a surface defect occurrence evaluation step for predicting the occurrence of surface defects from the solid phase thickness.
請求項1において、上記表面欠陥発生評価ステップは、上記固体相厚さが0.4mm以上の場所に表面欠陥が発生すると予測することを特徴とするダイカストシミュレーション方法。   2. The die casting simulation method according to claim 1, wherein the surface defect occurrence evaluation step predicts that a surface defect is generated at a location where the solid phase thickness is 0.4 mm or more. 請求項1において、上記表面欠陥発生評価ステップは、上記固体相厚さが0.1×e0.0297×Pmm(ここで、Pは鋳造圧力(MPa))以上の場所に表面欠陥が発生すると予測することを特徴とするダイカストシミュレーション方法。 3. The surface defect occurrence evaluation step according to claim 1, wherein the surface defect occurrence evaluation step predicts that a surface defect occurs at a location where the solid phase thickness is 0.1 × e 0.0297 × P mm (where P is a casting pressure (MPa)) or more. A die casting simulation method characterized by: 請求項1〜3のいずれか1項において、解析対象となる上記溶融金属材料は、凝固時に表皮側から凝固相が形成される表皮生成型凝固を示す合金であることを特徴とするダイカストシミュレーション方法。   The die casting simulation method according to any one of claims 1 to 3, wherein the molten metal material to be analyzed is an alloy exhibiting skin-generating solidification in which a solidified phase is formed from the skin side during solidification. . 請求項4において、上記溶融金属材料は、Si含有量が8〜14mass%であるAl−Si系アルミニウム合金であることを特徴とするダイカストシミュレーション方法。   5. The die casting simulation method according to claim 4, wherein the molten metal material is an Al—Si based aluminum alloy having a Si content of 8 to 14 mass%. 請求項1〜5のいずれか1項に記載のダイカストシミュレーション方法をコンピュータに実行させるためのダイカストシミュレーションプログラム。   A die casting simulation program for causing a computer to execute the die casting simulation method according to any one of claims 1 to 5. 請求項6に記載のダイカストシミュレーションプログラムを記録したコンピュータ読み取り可能な記憶媒体。   A computer-readable storage medium on which the die casting simulation program according to claim 6 is recorded. 請求項1〜5のいずれか1項に記載のダイカストシミュレーション方法を実行するよう構成されたコンピュータを有することを特徴とするダイカストシミュレーション装置。   A die-casting simulation apparatus comprising a computer configured to execute the die-casting simulation method according to claim 1.
JP2005065666A 2005-03-09 2005-03-09 DIE-CAST SIMULATION METHOD, ITS DEVICE, ITS PROGRAM, AND RECORDING MEDIUM CONTAINING THE PROGRAM Expired - Fee Related JP4373353B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005065666A JP4373353B2 (en) 2005-03-09 2005-03-09 DIE-CAST SIMULATION METHOD, ITS DEVICE, ITS PROGRAM, AND RECORDING MEDIUM CONTAINING THE PROGRAM

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005065666A JP4373353B2 (en) 2005-03-09 2005-03-09 DIE-CAST SIMULATION METHOD, ITS DEVICE, ITS PROGRAM, AND RECORDING MEDIUM CONTAINING THE PROGRAM

Publications (2)

Publication Number Publication Date
JP2006247684A true JP2006247684A (en) 2006-09-21
JP4373353B2 JP4373353B2 (en) 2009-11-25

Family

ID=37088698

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005065666A Expired - Fee Related JP4373353B2 (en) 2005-03-09 2005-03-09 DIE-CAST SIMULATION METHOD, ITS DEVICE, ITS PROGRAM, AND RECORDING MEDIUM CONTAINING THE PROGRAM

Country Status (1)

Country Link
JP (1) JP4373353B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010110792A (en) * 2008-11-06 2010-05-20 Toyota Central R&D Labs Inc Solidification analysis method for alloy molten metal and solidification analysis program thereof

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010110792A (en) * 2008-11-06 2010-05-20 Toyota Central R&D Labs Inc Solidification analysis method for alloy molten metal and solidification analysis program thereof

Also Published As

Publication number Publication date
JP4373353B2 (en) 2009-11-25

Similar Documents

Publication Publication Date Title
Ravi et al. Casting solidification analysis by modulus vector method
CN102395972B (en) Simulation of ejection after mold filling
Humphreys et al. Modelling and validation: Casting of Al and TiAl alloys in gravity and centrifugal casting processes
Yan et al. Numerical simulation of AZ91D alloy automobile plug in pressure die casting process
KR101312474B1 (en) Molten alloy solidification analyzing method and computer-readable storage medium with molten alloy solidification analyzing program for performing the same
Danylchenko Comparative analysis of computer systems for casting processes simulation
Hu et al. Flow modeling in high-pressure die-casting processes using SPH model
JP4582067B2 (en) Shape prediction apparatus, shape prediction method, and computer program therefor
JP4373353B2 (en) DIE-CAST SIMULATION METHOD, ITS DEVICE, ITS PROGRAM, AND RECORDING MEDIUM CONTAINING THE PROGRAM
Mitchell et al. Verified reduction of a model for a continuous casting process
JP4032848B2 (en) Molding simulation method, molding simulation apparatus, molding simulation program, and computer-readable recording medium recording the molding simulation program
Hernandez-Ortega et al. An experimental and numerical study of flow patterns and air entrapment phenomena during the filling of a vertical die cavity
Babaei et al. Improved advection algorithm of computational modeling of free surface flow using structured grids
TKADLEČKOVÁ et al. Testing of numerical model settings for simulation of steel ingot casting and solidification
JP4032755B2 (en) Molding simulation method, molding simulation apparatus, molding simulation program, and computer-readable recording medium recording the molding simulation program
Egole et al. Evaluation of the effect of cooling conditions and heat transfer coefficients on solidification of Al-Cu binary alloy in static casting process
JP2003334650A (en) Method and apparatus for simulating mold, mold simulating program and computer-readable recording medium having mold simulation program recorded thereon
Li Numerical simulation of flow and temperature evolution during the initial phase of steady-state solidification
Kimatsuka et al. Mold filling simulation for predicting gas porosity
Kuo et al. Development of an interactive simulation system for die cavity filling and its application to the operation of a low-pressure casting process
Catalina et al. Simplified pressure model for quantitative shrinkage porosity prediction in steel castings
JP2006159221A (en) Prediction method for occurrence of micro-porosity, prediction program for occurrence of micro-porosity and recording medium therefor
JPH0415761A (en) Analysis of temperature of thermal cycle structure and design device for metallic mold device system
JP4145078B2 (en) Casting simulation method
JP2003170269A (en) Simulated casting method, simulated casting device, simulated casting program and recording medium from which recorded simulated casting program can be read by computer

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070702

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20090813

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: 20090901

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090903

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120911

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313532

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120911

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120911

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130911

Year of fee payment: 4

LAPS Cancellation because of no payment of annual fees