JP4715462B2 - Molding analysis method and molding analysis apparatus - Google Patents

Molding analysis method and molding analysis apparatus Download PDF

Info

Publication number
JP4715462B2
JP4715462B2 JP2005335645A JP2005335645A JP4715462B2 JP 4715462 B2 JP4715462 B2 JP 4715462B2 JP 2005335645 A JP2005335645 A JP 2005335645A JP 2005335645 A JP2005335645 A JP 2005335645A JP 4715462 B2 JP4715462 B2 JP 4715462B2
Authority
JP
Japan
Prior art keywords
analysis
initial
predetermined time
coordinate information
point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2005335645A
Other languages
Japanese (ja)
Other versions
JP2007141046A (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.)
Toyo Seikan Kaisha Ltd
Original Assignee
Toyo Seikan Kaisha Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toyo Seikan Kaisha Ltd filed Critical Toyo Seikan Kaisha Ltd
Priority to JP2005335645A priority Critical patent/JP4715462B2/en
Publication of JP2007141046A publication Critical patent/JP2007141046A/en
Application granted granted Critical
Publication of JP4715462B2 publication Critical patent/JP4715462B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、多層成形解析に関し、特に、多層圧縮成形過程における中間層の位置を算出する成形解析手法、及び成形解析装置に関する。   The present invention relates to a multilayer molding analysis, and more particularly to a molding analysis method and a molding analysis apparatus for calculating the position of an intermediate layer in a multilayer compression molding process.

圧縮成形とは、溶融した樹脂を型に落とし込み、型締力にて所望の形状に成形する成形方法である。圧縮成形は、古くはガラスのコップや熱硬化性樹脂などに用いられてきた。
近年、熱可塑性樹脂への適用が増えており、従来、射出成形にて成形していた製品を圧縮成形で成形するケースもある。
Compression molding is a molding method in which molten resin is dropped into a mold and molded into a desired shape with a clamping force. Compression molding has long been used for glass cups and thermosetting resins.
In recent years, application to thermoplastic resins has increased, and there are cases where products that have been conventionally molded by injection molding are molded by compression molding.

圧縮成形は、射出成形に比べて低温成形が可能であるので、熱履歴による樹脂の劣化を抑制できる特長がある。また、多層化については、射出成形は層構成に制限があるのに比べ、圧縮成形は自由度が高く、さまざまな樹脂の多層化が可能である。さらに、生産性においては、射出成形は、多数個取りの金型で大量生産すると、広口成形品の場合は同じ金型面積に対する金型取数が少なくなるので生産効率が低下するが、圧縮成形は、ロータリー式の金型数を増やすのみで対応でき、この効率低下がないので相対的に有利となる。   Since compression molding can be performed at a lower temperature than injection molding, it has a feature that it can suppress deterioration of the resin due to thermal history. As for multilayering, compression molding has a higher degree of freedom than injection molding has a limited layer structure, and multilayering of various resins is possible. Furthermore, in terms of productivity, injection molding is mass-produced with multi-cavity molds, but in the case of wide-mouth molded products, the number of molds for the same mold area is reduced, so the production efficiency decreases. Can be handled only by increasing the number of rotary molds, and since this efficiency is not reduced, it is relatively advantageous.

ボトルやジャーなどの飲料容器のキャップは、熱可塑性樹脂を圧縮成形にて生産している例の一つである。
このような食品等の包装容器の分野では、内容品の酸化防止や防湿等のために、ガス透過性の少ない樹脂を使用して製品を多層化することが多い。
Caps for beverage containers such as bottles and jars are one example of producing thermoplastic resin by compression molding.
In the field of such packaging containers for foods and the like, products are often multilayered using a resin with low gas permeability in order to prevent oxidation and moisture prevention of the contents.

多層圧縮成形では、こういった製品性能を確保するため、層分布や厚みを規定の形態に成形することが求められる。特に、食品等の包装容器として用いる場合には、外層は衛生性の面から内容物に対して影響の少ない樹脂で完全に覆われ、酸素バリアなどの機能性樹脂は中間層として包含されつつも均一に広がることが求められる。   In multilayer compression molding, in order to ensure such product performance, it is required to mold the layer distribution and thickness into a prescribed form. In particular, when used as a packaging container for food, etc., the outer layer is completely covered with a resin that has little influence on the contents from the viewpoint of hygiene, while a functional resin such as an oxygen barrier is included as an intermediate layer. It is required to spread evenly.

中間層が最適に分布した所望の成形品を得るためには、圧縮成形される前の樹脂において、適切な形状に中間層が配置されることが重要である。また、圧縮成形される前の層全体の初期形状における中間層の最適な配置を求めるためには、圧縮工程中の樹脂挙動を把握する必要がある。
この樹脂挙動を実験的に求めることは極めて煩雑であり、現実的に行うことは困難である。
このため、数値解析によって、圧縮工程中の樹脂の挙動を可視化する手法の開発が求められていた。
In order to obtain a desired molded article in which the intermediate layer is optimally distributed, it is important that the intermediate layer is arranged in an appropriate shape in the resin before being subjected to compression molding. In addition, in order to obtain the optimal arrangement of the intermediate layer in the initial shape of the entire layer before compression molding, it is necessary to grasp the resin behavior during the compression process.
Obtaining this resin behavior experimentally is very complicated and difficult to carry out realistically.
For this reason, the development of a method for visualizing the behavior of the resin during the compression process by numerical analysis has been demanded.

このような樹脂挙動を可視化する手法としては、従来、アダプティブメッシュ法が用いられている(例えば、非特許文献1参照)。
アダプティブメッシュ法とは、成形に応じて変形する領域のメッシュを繰り返し修正しながら解析を進める手法である。
As a method for visualizing such resin behavior, an adaptive mesh method has been conventionally used (see, for example, Non-Patent Document 1).
The adaptive mesh method is a method of performing analysis while repeatedly correcting a mesh in a region that is deformed according to molding.

すなわち、このアダプティブメッシュ法では、図11に示すように、樹脂の圧縮成形により、メッシュも樹脂全体の変形に応じて変形する。しかし、樹脂の形状変化が大きい場合には、メッシュの形状が、それ以上計算を続行することができない程度にまで変形してしまう。
そこで、図12に示すように、メッシュの潰れが大きくなり、それ以上計算の続行が不可能になるまでに、メッシュの再設定(リメッシング)を繰り返す。
これによって、圧縮成形によって樹脂の形状が大きく変化する場合でも、成形が完了するまでの樹脂挙動を逐次追跡することを可能にしている。
That is, in this adaptive mesh method, as shown in FIG. 11, the mesh is also deformed in accordance with the deformation of the entire resin by compression molding of the resin. However, when the shape change of the resin is large, the shape of the mesh is deformed to such an extent that the calculation cannot be continued any more.
Therefore, as shown in FIG. 12, mesh resetting (remeshing) is repeated until the collapse of the mesh becomes large and further calculation cannot be continued.
Thus, even when the shape of the resin changes greatly due to compression molding, it is possible to sequentially track the resin behavior until the molding is completed.

このようなアダプティブメッシュ法を多層圧縮成形に適用する場合は、それぞれの樹脂層に対応した領域ごとにリメッシングを行うことになる。リメッシングの際には界面はそのまま保持されるので、層間の境界を明確に得ることができる。   When such an adaptive mesh method is applied to multilayer compression molding, remeshing is performed for each region corresponding to each resin layer. Since the interface is maintained as it is during remeshing, the boundary between the layers can be clearly obtained.

深堀 穂高:成形加工シンポジア‘02要旨集,283(2002)Hodaka Fukahori: Molding Symposia '02 Abstract, 283 (2002)

しかしながら、このアダプティブメッシュ法には、中間層の量が少ない樹脂構成に適用した場合、成形により中間層が薄くなりすぎると、その層にメッシュを作成することが不可能になり、計算を続行することができなくなってしまうという問題があった。
このような問題に対し、より細かなメッシュを作成することで対処することも考えられる。しかしながら、メッシュが細かくなればなるほど、計算時間が長くなることから、中間層の厚みによっては、現実的に用いることは困難である。
したがって、近年増えつつある機能性樹脂などの薄い中間層については、圧縮成形における樹脂挙動を解析することが、大変難しいという問題があった。
このため、このような中間層を備えた樹脂の多層圧縮成形において、樹脂の挙動を適切に解析可能な手法が求められている。
However, when this adaptive mesh method is applied to a resin structure having a small amount of intermediate layer, if the intermediate layer becomes too thin due to molding, it becomes impossible to create a mesh in that layer, and the calculation is continued. There was a problem that it was impossible.
It may be possible to cope with such a problem by creating a finer mesh. However, the finer the mesh, the longer the calculation time. Therefore, depending on the thickness of the intermediate layer, it is difficult to use it practically.
Therefore, for thin intermediate layers such as functional resins, which have been increasing in recent years, it has been very difficult to analyze the resin behavior in compression molding.
For this reason, in the multilayer compression molding of resin provided with such an intermediate layer, a method capable of appropriately analyzing the behavior of the resin is required.

そこで、本発明の発明者は鋭意検討した結果、アダプティブメッシュ法による多層圧縮成形解析の欠点を補うべく、点追跡法及び大域的追跡法を開発した。
点追跡法とは、圧縮成形工程の解析をアダプティブメッシュ法を用いて単層で実施し、この解析結果を用いて、層全体の初期形状上に点を配置し、その点を追跡することにより成形後の中間層の位置を求める手法である。
また、大域的追跡法とは、中間層に相当する部分だけでなく、層全体に点を配置して点追跡を行い、成形の前後における樹脂の移動の対応関係を求めることで、成形後の所望の中間層に位置する点群のそれぞれの初期座標値を、初期の樹脂形状に重ねることにより、初期の中間層形状のあるべき配置を求める手法である。
Therefore, as a result of intensive studies, the inventors of the present invention have developed a point tracking method and a global tracking method to compensate for the disadvantages of multilayer compression molding analysis by the adaptive mesh method.
In the point tracking method, the compression molding process is analyzed in a single layer using the adaptive mesh method, and points are placed on the initial shape of the entire layer using this analysis result, and the points are tracked. This is a method for obtaining the position of the intermediate layer after molding.
In addition, the global tracking method is not only the portion corresponding to the intermediate layer, but also by placing points on the entire layer to perform point tracking, and by determining the correspondence of resin movement before and after molding, In this method, the initial coordinate value of each point group located in a desired intermediate layer is overlapped with the initial resin shape to obtain the desired arrangement of the initial intermediate layer shape.

すなわち、本発明は、多層圧縮成形において、成形後の中間層の位置を求めることと、成形後の中間層の分布が最適となる、初期の樹脂形状における中間層分布を生成することの可能な成形解析手法、及び成形解析装置の提供を目的とする。   That is, according to the present invention, in the multilayer compression molding, the position of the intermediate layer after molding can be obtained, and the intermediate layer distribution in the initial resin shape can be generated in which the distribution of the intermediate layer after molding is optimal. An object is to provide a molding analysis method and a molding analysis apparatus.

上記目的を達成するため、本発明の成形解析手法は、成形解析装置に、成形過程における層分布を算出させる成形解析手法であって、成形解析装置における層形状変化解析手段が、層を所定の解析条件にもとづき解析して、層に割り振られたメッシュの節点の座標情報と該節点の速度ベクトルを含む物理量とからなる解析結果を所定の時間ステップごとに算出し、成形解析装置における点追跡手段が、所定の時間ステップごとに、解析結果を読み込んで、この解析結果と層内に想定された一又は二以上の任意の点の座標情報とにもとづき各点における速度ベクトルを算出するとともに、所定の時間ステップ後の各点の座標情報を算出し、かつ、解析結果の読み込み、速度ベクトルの算出、及び座標情報の算出を所定時間が経過するまで繰り返し、層内に想定された任意の点の所定時間後の座標情報を算出する手法としてある。   In order to achieve the above object, the molding analysis method of the present invention is a molding analysis method for causing a molding analysis device to calculate a layer distribution in a molding process, wherein the layer shape change analysis means in the molding analysis device defines a predetermined layer. Analyzing based on the analysis conditions, calculating the analysis result consisting of the coordinate information of the nodes of the mesh assigned to the layer and the physical quantity including the velocity vector of the nodes for each predetermined time step, and the point tracking means in the molding analyzer However, at every predetermined time step, the analysis result is read, the velocity vector at each point is calculated based on the analysis result and the coordinate information of one or more arbitrary points assumed in the layer. The coordinate information of each point after the time step is calculated, and the analysis result reading, the velocity vector calculation, and the coordinate information calculation are repeated until a predetermined time elapses. Is a method of calculating the coordinate information of the predetermined time after the arbitrary point that was assumed in the layer.

成形解析手法をこのような手法にすれば、層全体の初期形状における中間層の位置に点を配置して、成形過程におけるその点の移動を追跡することができる。
このため、あらゆる初期中間層形状について、成形による中間層分布の変化を計算することができ、従来は困難であった少量の中間層の形状変化を適切に把握することが可能となる。
また、本発明によれば、層形状変化の解析を一度行うのみで、様々な形状の中間層形状について点追跡を行うことができる。
このため、従来のアダプティブメッシュ法のみによる成形解析手法に比べて、より効率的な成形解析を行うことが可能となる。
If the forming analysis method is such a method, a point can be arranged at the position of the intermediate layer in the initial shape of the entire layer, and the movement of the point in the forming process can be traced.
For this reason, for any initial intermediate layer shape, the change in the intermediate layer distribution due to molding can be calculated, and a small amount of intermediate layer shape change, which has been difficult in the past, can be properly grasped.
In addition, according to the present invention, it is possible to perform point tracking for various shapes of the intermediate layer by performing a single layer shape change analysis.
For this reason, it becomes possible to perform more efficient forming analysis as compared with the forming analysis method based only on the conventional adaptive mesh method.

また、本発明の成形解析手法は、成形解析装置における大域的点追跡手段が、所定の時間ステップごとに、解析結果を読み込んで、解析結果における最初の時間ステップの全ての節点の座標情報及び初期座標値を大域的点追跡結果記憶手段に記憶させ、当該節点を点追跡手段により所定の時間ステップの間点追跡して、所定の時間ステップ後の各節点の座標情報を算出し、当該所定の時間ステップ後の各節点の初期座標値分布を生成して、かつ、解析結果の読み込み、当該解析結果にもとづく所定の時間ステップ後の各節点の座標情報及び初期座標値の記憶、所定の時間ステップ後の各節点に対する点追跡によるさらに所定の時間ステップ後の各節点の座標情報の算出、及び当該所定の時間ステップ後の各節点の初期座標値分布の生成を所定時間が経過するまで繰り返し、解析結果における最初の時間ステップの全ての節点の所定時間後の座標情報と初期座標値分布を算出する手法としてある。   Further, according to the molding analysis method of the present invention, the global point tracking means in the molding analysis apparatus reads the analysis result for each predetermined time step, and the coordinate information and initial values of all the nodes in the first time step in the analysis result. The coordinate value is stored in the global point tracking result storage unit, the node is tracked for a predetermined time step by the point tracking unit, the coordinate information of each node after the predetermined time step is calculated, Generate initial coordinate value distribution of each node after time step, read analysis result, store coordinate information and initial coordinate value of each node after predetermined time step based on analysis result, predetermined time step Predetermined calculation of coordinate information of each node after a predetermined time step by point tracking for each subsequent node and generation of initial coordinate value distribution of each node after the predetermined time step Repeat to between elapses, there as a method of calculating the coordinate information and the initial coordinate value distribution after a predetermined time of all the nodes in the first time step in the analysis results.

成形解析手法をこのような手法にすれば、層に割り振られたメッシュの全ての節点について、初期座標値を保有させつつ点追跡を行うことができる。
このため、成形後の層全体の最終形状における追跡された点の初期座標値分布を把握することが可能となる。
If the forming analysis method is such a method, point tracking can be performed while retaining the initial coordinate values for all the nodes of the mesh assigned to the layer.
For this reason, it is possible to grasp the initial coordinate value distribution of the tracked points in the final shape of the entire layer after molding.

また、本発明の成形解析手法は、成形解析装置における初期中間層形状生成手段が、大域的点追跡手段により算出された全ての節点の所定時間後の座標情報により構成される層全体の最終形状上に特定された任意の中間層形状における多数の点の座標情報にもとづいて、当該点に対応する初期座標値を大域的点追跡結果記憶手段に記憶された座標情報及び初期座標値にもとづき算出し、これらの初期座標値の点を層全体の初期形状情報に重ね合わせて、中間層形状に対応する初期形状情報を生成する手法としてある。   Further, the molding analysis method of the present invention is such that the initial intermediate layer shape generation means in the molding analysis apparatus has a final shape of the entire layer constituted by coordinate information after a predetermined time of all nodes calculated by the global point tracking means. Based on the coordinate information of a large number of points in the arbitrary intermediate layer shape specified above, the initial coordinate value corresponding to the point is calculated based on the coordinate information and the initial coordinate value stored in the global point tracking result storage means. Then, these initial coordinate value points are superimposed on the initial layer shape information of the entire layer to generate initial shape information corresponding to the intermediate layer shape.

成形解析手法をこのような手法にすれば、成形後の任意の形状の中間層に対応する、層全体の初期形状における中間層の初期形状情報を生成することができる。
このため、少ない量の中間層についても、成形後に最適な中間層を得ることが可能な、層全体の初期形状における理想的な中間層の分布を把握することが可能となる。
If the forming analysis method is such a method, it is possible to generate the initial shape information of the intermediate layer in the initial shape of the entire layer corresponding to the intermediate layer having an arbitrary shape after forming.
For this reason, even with a small amount of the intermediate layer, it is possible to obtain an ideal intermediate layer distribution in the initial shape of the entire layer, which can obtain an optimal intermediate layer after molding.

また、本発明の成形解析手法は、点追跡手段が、層内に想定された一又は二以上の任意の点の座標情報が解析結果における節点の座標情報に一致しない場合に、解析結果における、任意の点の周囲の節点に対応する速度ベクトルを用いて、補間計算により任意の点の速度ベクトルを算出する手法としてある。   Further, in the molding analysis method of the present invention, the point tracking means, when the coordinate information of one or more arbitrary points assumed in the layer does not match the coordinate information of the node in the analysis result, in the analysis result, This is a technique for calculating a velocity vector at an arbitrary point by interpolation calculation using velocity vectors corresponding to nodes around the arbitrary point.

成形解析手法をこのような手法にすれば、層形状変化の解析によって、メッシュの節点が割り振られていない部分における点についても、補間計算を行うことによって、その点の速度ベクトルなどの物理量を算出することができる。
このため、このような点についても適切に点追跡を行うことが可能となる。
If the forming analysis method is such a method, the physical quantity such as the velocity vector of the point is calculated by performing the interpolation calculation for the point in the part where the mesh node is not allocated by analyzing the layer shape change. can do.
For this reason, it becomes possible to perform point tracking appropriately also about such a point.

また、本発明の成形解析装置は、成形過程における層分布を算出する成形解析装置であって、層を所定の解析条件にもとづき解析し、層に割り振られたメッシュの節点の座標情報と該節点の速度ベクトルを含む物理量とからなる解析結果を所定の時間ステップごとに算出する層形状変化解析手段と、所定の時間ステップごとに、解析結果を読み込んで、この解析結果と層内に想定された一又は二以上の任意の点の座標情報とにもとづき各点における速度ベクトルを算出するとともに、所定の時間ステップ後の各点の座標情報を算出し、かつ、解析結果の読み込み、速度ベクトルの算出、及び座標情報の算出を所定時間が経過するまで繰り返し、層内に想定された任意の点の所定時間後の座標情報を算出する点追跡手段と、を備えた構成としてある。   The molding analysis apparatus of the present invention is a molding analysis apparatus that calculates a layer distribution in a molding process, analyzes a layer based on a predetermined analysis condition, and coordinates information of the nodes of the mesh allocated to the layer and the nodes. The layer shape change analysis means for calculating the analysis result consisting of the physical quantity including the velocity vector at every predetermined time step, and reading the analysis result at every predetermined time step, and this analysis result and the layer assumed Calculate the velocity vector at each point based on the coordinate information of one or more arbitrary points, calculate the coordinate information of each point after a predetermined time step, read the analysis result, and calculate the velocity vector And a point tracking unit that repeats the calculation of the coordinate information until a predetermined time elapses, and calculates coordinate information after a predetermined time of an arbitrary point assumed in the layer. That.

成形解析装置をこのような構成にすれば、層全体の初期形状における任意の点について、成形過程における移動を追跡することができるため、少量の中間層の形状変化を把握することが可能となる。
また、層形状変化の解析を一度行うのみで、様々な形状の中間層形状について点追跡を行うことができるため、従来の成形解析手法に比べて、より効率的な成形解析を行うことが可能となる。
With such a configuration of the molding analysis apparatus, it is possible to track the movement in the molding process at any point in the initial shape of the entire layer, and therefore it is possible to grasp the shape change of a small amount of the intermediate layer. .
In addition, it is possible to perform more efficient forming analysis than conventional forming analysis methods because it is possible to perform point tracking for various shapes of intermediate layer shapes by performing layer shape change analysis only once. It becomes.

また、本発明の成形解析装置は、所定の時間ステップごとに、解析結果を読み込んで、解析結果における最初の時間ステップの全ての節点の座標情報及び初期座標値を大域的点追跡結果記憶手段に記憶させ、当該節点を点追跡手段により所定の時間ステップの間点追跡して、所定の時間ステップ後の各節点の座標情報を算出し、当該所定の時間ステップ後の各節点の初期座標値分布を生成して、かつ、解析結果の読み込み、当該解析結果にもとづく所定の時間ステップ後の各節点の座標情報及び初期座標値の記憶、所定の時間ステップ後の各節点に対する点追跡によるさらに所定の時間ステップ後の各節点の座標情報の算出、及び当該所定の時間ステップ後の各節点の初期座標値分布の生成を所定時間が経過するまで繰り返し、解析結果における最初の時間ステップの全ての節点の所定時間後の座標情報と初期座標値分布を算出する大域的点追跡手段を備えた構成としてある。   Further, the molding analysis apparatus of the present invention reads the analysis result at every predetermined time step, and stores the coordinate information and initial coordinate values of all the nodes in the first time step in the analysis result in the global point tracking result storage means. Storing the node, tracking the node for a predetermined time step by the point tracking means, calculating coordinate information of each node after the predetermined time step, and initial coordinate value distribution of each node after the predetermined time step And by reading the analysis result, storing the coordinate information and initial coordinate value of each node after a predetermined time step based on the analysis result, and further tracking the point for each node after the predetermined time step. The calculation of the coordinate information of each node after the time step and the generation of the initial coordinate value distribution of each node after the predetermined time step are repeated until the predetermined time elapses. A configuration equipped with first global point tracking means for calculating the coordinate information and the initial coordinate value distribution after a predetermined time of all the nodes in time step that.

成形解析装置をこのような構成にすれば、層に割り振られたメッシュの全ての節点につき初期座標値を保有させて点追跡を行うことができ、成形後の層全体の最終形状における追跡された点の初期座標値分布を把握することが可能となる。   If the forming analysis apparatus is configured in this way, it is possible to carry out point tracking by retaining the initial coordinate values for all the nodes of the mesh assigned to the layer, and it is tracked in the final shape of the entire layer after forming. It becomes possible to grasp the initial coordinate value distribution of the points.

また、本発明の成形解析装置は、大域的点追跡手段により算出された全ての節点の所定時間後の座標情報により構成される層全体の最終形状上に特定された任意の中間層形状における多数の点の座標情報にもとづいて、当該点に対応する初期座標値を大域的点追跡結果記憶手段に記憶された座標情報及び初期座標値にもとづき算出し、これらの初期座標値の点を層全体の初期形状情報に重ね合わせて、中間層形状に対応する初期形状情報を生成する初期中間層形状生成手段を備えた構成としてある。   In addition, the molding analysis apparatus of the present invention has a large number of arbitrary intermediate layer shapes specified on the final shape of the entire layer constituted by coordinate information after a predetermined time of all nodes calculated by the global point tracking means. Based on the coordinate information of the point, the initial coordinate value corresponding to the point is calculated based on the coordinate information and the initial coordinate value stored in the global point tracking result storage means, and the point of these initial coordinate values is calculated for the entire layer. The initial intermediate layer shape generating means for generating initial shape information corresponding to the intermediate layer shape is superimposed on the initial shape information.

成形解析装置をこのような構成にすれば、成形後に最適な中間層を得ることが可能な、層全体の初期形状における中間層の初期形状情報を生成することができ、理想的な中間層の初期形状を把握することが可能となる。   By configuring the molding analysis device in this way, it is possible to generate the initial shape information of the intermediate layer in the initial shape of the entire layer, which can obtain the optimal intermediate layer after molding, and the ideal intermediate layer It becomes possible to grasp the initial shape.

また、本発明の成形解析装置は、点追跡手段が、層内に想定された一又は二以上の任意の点の座標情報が解析結果における節点の座標情報に一致しない場合に、解析結果における、任意の点の周囲の節点に対応する速度ベクトルを用いて、補間計算により任意の点の速度ベクトルを算出する構成としてある。   Further, in the molding analysis apparatus of the present invention, when the point tracking means, the coordinate information of one or more arbitrary points assumed in the layer does not match the coordinate information of the nodes in the analysis result, in the analysis result, A velocity vector at an arbitrary point is calculated by interpolation calculation using velocity vectors corresponding to nodes around the arbitrary point.

成形解析装置をこのような構成にすれば、層形状変化の解析で、メッシュの節点が割り振られていない部分についても、補間計算により速度ベクトルを算出することができ、適切に点追跡を行うことが可能となる。   If the forming analysis device is configured in this way, the velocity vector can be calculated by interpolation calculation even for the part where the mesh nodes are not allocated in the analysis of the layer shape change, and the point tracking should be performed appropriately. Is possible.

本発明によれば、圧縮成形のような大変形の成形工程における多層の樹脂挙動を可視化することができる。また、層全体に配置した多数の点を追跡することにより、成形後に所望の層配置になるような初期の層構成を予測することができる。
このため、中間層の厚みが小さい場合でも、解析を行うことができ、層全体の初期形状における理想的な中間層分布を得ることが可能となる。また、層全体の形状変化の解析を一度行うのみで、複数の中間層形状について解析を行うことが可能となる。
したがって、本発明によれば、多層圧縮成形等を用いた製品開発を効率化することが可能となる。
According to the present invention, it is possible to visualize the multilayer resin behavior in a large deformation molding process such as compression molding. In addition, by tracking a large number of points arranged in the entire layer, it is possible to predict an initial layer configuration that achieves a desired layer arrangement after molding.
For this reason, even when the thickness of the intermediate layer is small, analysis can be performed, and an ideal intermediate layer distribution in the initial shape of the entire layer can be obtained. In addition, it is possible to analyze a plurality of intermediate layer shapes by performing analysis of the shape change of the entire layer once.
Therefore, according to the present invention, product development using multilayer compression molding or the like can be made efficient.

以下、本発明に係る圧縮成形解析装置の好ましい実施形態について、図面を参照しつつ説明する。
なお、以下の実施形態に示す本発明の圧縮成形解析装置は、プログラムに制御されたコンピュータにより動作するようになっている。コンピュータのCPUは、プログラムにもとづいてコンピュータの各構成要素に指令を送り、圧縮成形解析装置の動作に必要となる所定の処理、例えば、層形状変化の解析処理、追跡点の座標取得処理、追跡点の速度ベクトルの算出処理、初期座標分布の生成処理等を行わせる。このように、本発明の圧縮成形解析装置における各処理,動作は、プログラムとコンピュータとが協働した具体的手段により実現できるものである。
Hereinafter, a preferred embodiment of a compression molding analysis apparatus according to the present invention will be described with reference to the drawings.
In addition, the compression molding analysis apparatus of this invention shown to the following embodiment operate | moves with the computer controlled by the program. The CPU of the computer sends a command to each component of the computer based on the program, and performs predetermined processing necessary for the operation of the compression molding analysis apparatus, for example, analysis processing of layer shape change, coordinate acquisition processing of tracking points, tracking A point velocity vector calculation process, an initial coordinate distribution generation process, and the like are performed. Thus, each process and operation in the compression molding analysis apparatus of the present invention can be realized by specific means in which the program and the computer cooperate.

プログラムは予めROM,RAM等の記録媒体に格納され、コンピュータに実装された記録媒体から当該コンピュータにプログラムを読み込ませて実行されるが、例えば通信回線を介してコンピュータに読み込ませることもできる。
また、プログラムを格納する記録媒体は、例えば半導体メモリ,磁気ディスク,光ディスク、その他任意のコンピュータで読取り可能な任意の記録手段により構成できる。
The program is stored in advance in a recording medium such as a ROM and a RAM, and is executed by causing the computer to read the program from a recording medium mounted on the computer. For example, the program may be read by the computer via a communication line.
Further, the recording medium for storing the program can be constituted by, for example, a semiconductor memory, a magnetic disk, an optical disk, or any other recording means that can be read by any computer.

[第一実施形態]
まず、本発明の第一実施形態の構成について、図1を参照して説明する。同図は、本実施形態の圧縮成形解析装置の構成を示すブロック図である。
図1に示すように、本実施形態の圧縮成形解析装置10は、モデルファイル記憶部10a,解析条件入力部10b,層形状変化解析部10c,解析結果記憶部10d,中間層形状情報入力部10e,点追跡部10f,点追跡結果記憶部10gを備えている。
[First embodiment]
First, the configuration of the first embodiment of the present invention will be described with reference to FIG. This figure is a block diagram showing the configuration of the compression molding analysis apparatus of the present embodiment.
As shown in FIG. 1, the compression molding analysis apparatus 10 of this embodiment includes a model file storage unit 10a, an analysis condition input unit 10b, a layer shape change analysis unit 10c, an analysis result storage unit 10d, and an intermediate layer shape information input unit 10e. , A point tracking unit 10f and a point tracking result storage unit 10g.

モデルファイル記憶部10aは、層形状変化の解析処理における解析対象データとしてのモデルファイルを記憶する。
このモデルファイルは、金型の形状情報、初期状態における樹脂領域の形状(層全体の初期形状)情報等を保有する。
また、モデルファイルは、CADソフトウェアやメッシュ生成ソフトウェアなどを用いて、流体解析ソフトウェアの実行により形成される層形状変化解析部10cが読み込めるような状態に作成されている。
なお、形状情報とは、多数の点の座標情報により構成される。本実施形態では、座標情報として二次元の情報を想定しているが、座標情報を三次元の情報とすることも可能である。
The model file storage unit 10a stores a model file as analysis target data in the layer shape change analysis process.
This model file holds information on the shape of the mold, information on the shape of the resin region in the initial state (initial shape of the entire layer), and the like.
The model file is created in a state that can be read by the layer shape change analysis unit 10c formed by the execution of the fluid analysis software using CAD software, mesh generation software, or the like.
The shape information is constituted by coordinate information of a large number of points. In the present embodiment, two-dimensional information is assumed as the coordinate information, but the coordinate information may be three-dimensional information.

解析条件入力部10bは、圧縮成形工程における層形状変化を解析するための条件を圧縮成形解析装置10に入力し、層形状変化解析部10cによりこの解析条件を使用可能にするため、圧縮成形解析装置10に記憶させる。
この解析条件としては、金型及び樹脂の物性、成形速度、金型と樹脂間の熱、すべり境界条件、メッシュ修正の条件、解析の時間ステップ等を挙げることができる。
The analysis condition input unit 10b inputs a condition for analyzing the layer shape change in the compression molding process to the compression molding analysis apparatus 10, and the layer shape change analysis unit 10c enables the use of this analysis condition. It is stored in the device 10.
Examples of the analysis conditions include physical properties of the mold and the resin, molding speed, heat between the mold and the resin, slip boundary conditions, mesh correction conditions, analysis time step, and the like.

層形状変化解析部10cは、層の形状変化を解析するものであり、所定の流体解析ソフトウェアを実行することで、圧縮成形解析装置10におけるCPU及びメモリ上に構成される。
この流体解析ソフトウェアとしては、例えばFLUENT社から販売された「POLYFLOW」を用いることができる。
層形状変化解析部10cは、モデルファイル記憶部10aから流体解析のためのモデルファイルを読み込むとともに、解析条件入力部10bから入力された解析条件を読み込んで、解析処理を実行する。
The layer shape change analysis unit 10c analyzes a layer shape change, and is configured on a CPU and a memory in the compression molding analysis device 10 by executing predetermined fluid analysis software.
As this fluid analysis software, for example, “POLYFLOW” sold by FLUENT can be used.
The layer shape change analysis unit 10c reads a model file for fluid analysis from the model file storage unit 10a, and also reads analysis conditions input from the analysis condition input unit 10b, and executes analysis processing.

この層形状変化解析部10cによる解析処理は、本発明の発明者により提案された上記非特許文献1に記載の「圧縮成形工程の数値解析」に示された手順により実行することができる。
すなわち、層形状変化解析部10cは、アダプティブメッシュ法を用いた解析手法により、層全体の形状変化を解析条件にもとづき計算する。
The analysis processing by the layer shape change analysis unit 10c can be executed by the procedure shown in “Numerical analysis of compression molding process” described in Non-Patent Document 1 proposed by the inventor of the present invention.
That is, the layer shape change analysis unit 10c calculates the shape change of the entire layer based on the analysis conditions by an analysis method using the adaptive mesh method.

このとき、層形状変化解析部10cは、成形により潰れたメッシュが計算不可能になる前に逐次修正(リメッシング)して、その状態から新たに形状変化の計算を続けることを最終形状まで繰り返す。
これによって、層全体の初期形状から最終形状までの形状変化を解析する。
At this time, the layer shape change analysis unit 10c repeats correction (remeshing) before the mesh crushed by molding becomes impossible to calculate, and continues to calculate the shape change from that state until the final shape.
Thereby, the shape change from the initial shape to the final shape of the entire layer is analyzed.

また、この際、形状により、必要に応じて、適宜メッシュサイズや解析の時間ステップを変更することができる。
なお、層の形状変化が小さく、メッシュを修正する必要がない場合、層形状変化解析部10cは、アダプティブメッシュ法を用いることなく、層形状変化の解析を行うことができる。
At this time, depending on the shape, the mesh size and the time step of analysis can be changed as needed.
When the layer shape change is small and the mesh need not be corrected, the layer shape change analysis unit 10c can analyze the layer shape change without using the adaptive mesh method.

解析結果記憶部10dは、層形状変化解析部10cにより算出された解析結果を記憶する。
この解析結果には、層全体を初期形状から最終形状まで成形する過程における所定の時間ステップごとに、メッシュにおける全ての節点の座標情報及び物理量が含まれている。
物理量は、上記解析条件として、算出対象とするものを設定することができる。例えば、物理量に、各節点における速度ベクトル,圧力,温度,応力等を含めることができる。
The analysis result storage unit 10d stores the analysis result calculated by the layer shape change analysis unit 10c.
This analysis result includes coordinate information and physical quantities of all nodes in the mesh for each predetermined time step in the process of forming the entire layer from the initial shape to the final shape.
As the physical quantity, a physical quantity to be calculated can be set as the analysis condition. For example, the physical quantity can include a velocity vector, pressure, temperature, stress and the like at each node.

中間層形状情報入力部10eは、点追跡の対象とする中間層形状に対応する点の座標情報を圧縮成形解析装置10に入力し、点追跡部10fで使用可能にするために、圧縮成形解析装置10に記憶させる。
この中間層形状に対応する点としては、その中間層形状の領域の輪郭上及び内部に、中間層形状を判別することができる程度のものを特定する。なお、この特定については、コンピュータにより自動的に行う他、人手により行うこともできる。
The intermediate layer shape information input unit 10e inputs the coordinate information of the point corresponding to the intermediate layer shape that is the target of point tracking to the compression molding analysis apparatus 10 so that it can be used by the point tracking unit 10f. It is stored in the device 10.
As the points corresponding to the intermediate layer shape, a point that can determine the intermediate layer shape is specified on the contour and inside the region of the intermediate layer shape. Note that this specification can be performed manually as well as automatically by a computer.

点追跡部10fは、層全体を初期形状から最終形状まで成形する過程における所定の時間ステップごとに、中間層形状に対応する点の座標の変化を追跡するものであり、所定の流体解析ソフトを実行することで、圧縮成形解析装置10におけるCPU及びメモリ上に構成される。
この流体解析ソフトウェアとしては、例えば株式会社ヴァイナスにより販売された流体解析用ポストプロセッサソフトウェア「FIELDVIEW」(開発元:米国Inteligent Light社)を用いることができる。
The point tracking unit 10f tracks the change in the coordinates of the point corresponding to the intermediate layer shape at every predetermined time step in the process of forming the entire layer from the initial shape to the final shape. By executing, it is configured on the CPU and memory in the compression molding analysis apparatus 10.
As the fluid analysis software, for example, post-processor software “FIELDVIEW” (developed by Intelligent Light Corporation of the United States) sold by Vinus Corporation can be used.

点追跡部10fは、解析結果記憶部10dから層形状変化の解析結果を読み込むとともに、中間層形状情報入力部10eから入力された中間層形状に対応する点の座標情報を読み込んで、点追跡処理を実行する。
この点追跡処理において、点追跡部10fは、点追跡の対象となる点(以下、追跡点と称する場合がある。)の座標における速度ベクトルを、層形状変化の解析結果にもとづき算出する。そして、この速度ベクトルと点の座標情報にもとづいて、次の時間ステップにおける点の座標情報を算出する。
The point tracking unit 10f reads the analysis result of the layer shape change from the analysis result storage unit 10d, and also reads the coordinate information of the point corresponding to the intermediate layer shape input from the intermediate layer shape information input unit 10e, and performs point tracking processing Execute.
In this point tracking process, the point tracking unit 10f calculates a velocity vector at the coordinates of a point to be point-tracked (hereinafter sometimes referred to as a tracking point) based on the analysis result of the layer shape change. Then, based on the velocity vector and the point coordinate information, the point coordinate information in the next time step is calculated.

このとき、点追跡部10fは、追跡点が解析結果における節点に一致する場合は、その節点に対応づけて解析結果記憶部10dに記憶されている速度ベクトルを用いて、これをそのまま追跡点の速度ベクトルとして算出する。
また、追跡点が解析結果における節点に一致しない場合は、その節点の周囲を取り囲む節点に対応づけて解析結果記憶部10dに記憶されている速度ベクトルを用いて、後述する補間計算により速度ベクトルを算出する。
点追跡部10fは、これらの計算処理を、全ての追跡対象とする点について、解析終了時間まで、所定の時間ステップごとに実行する。
これによって、点追跡部10fは、初期中間層形状に対応する各点が、解析終了時間にどの座標に移動するかを計算でき、成形後の中間層分布を算出することができる。
At this time, if the tracking point matches the node in the analysis result, the point tracking unit 10f uses the velocity vector stored in the analysis result storage unit 10d in association with the node, and uses this as it is for the tracking point. Calculated as a velocity vector.
If the tracking point does not coincide with the node in the analysis result, the velocity vector is calculated by interpolation calculation described later using the velocity vector stored in the analysis result storage unit 10d in association with the node surrounding the node. calculate.
The point tracking unit 10f executes these calculation processes for every point to be tracked for every predetermined time step until the analysis end time.
Thereby, the point tracking unit 10f can calculate to which coordinate each point corresponding to the initial intermediate layer shape moves at the analysis end time, and can calculate the intermediate layer distribution after molding.

図2は、このような点追跡の様子を示す図である。本実施形態の圧縮成形解析装置10による点追跡法によれば、同図に示すように、初期形状(a)上に、中間層の形状に応じた点を追跡対象として配置し(b)、この点を上記点追跡処理により追跡することで、成形後の中間層分布(c)を算出することができる。   FIG. 2 is a diagram showing the state of such point tracking. According to the point tracking method by the compression molding analysis apparatus 10 of the present embodiment, as shown in the figure, a point corresponding to the shape of the intermediate layer is placed on the initial shape (a) as a tracking target (b), By tracking this point by the point tracking process, the intermediate layer distribution (c) after molding can be calculated.

点追跡結果記憶部10gは、点追跡部10fにより算出された追跡結果を記憶する。
この追跡結果には、所定の時間ステップごとに、メッシュにおける全ての節点の座標情報及び物理量が含まれているとともに、追跡対象とした全ての点の座標情報が含まれている。
The point tracking result storage unit 10g stores the tracking result calculated by the point tracking unit 10f.
This tracking result includes coordinate information and physical quantities of all nodes in the mesh for every predetermined time step, and also includes coordinate information of all points to be tracked.

次に、本実施形態の圧縮成形解析装置10における処理手順について、図3及び図4を参照して説明する。図3は、本実施形態の圧縮成形解析装置10により実行される層形状変化の解析処理を示すフローチャートであり、図4は、同装置により実行される点追跡処理を示すフローチャートである。なお、図3及び図4において、点線で表示されている部分は、人により行われる部分である。以下の実施形態においても同様の形式で表示している。   Next, a processing procedure in the compression molding analysis apparatus 10 of the present embodiment will be described with reference to FIGS. FIG. 3 is a flowchart showing a layer shape change analysis process executed by the compression molding analysis apparatus 10 of the present embodiment, and FIG. 4 is a flowchart showing a point tracking process executed by the apparatus. In FIGS. 3 and 4, a portion displayed by a dotted line is a portion performed by a person. In the following embodiments, the same format is used.

まず、図3を参照して、層形状変化の解析処理について説明する。
はじめに、層形状変化の解析処理における解析対象データとしてのモデルファイルを作成する(ステップ10)。
このモデルファイルは、金型形状情報、及び層全体の初期形状情報を保有する。
そして、圧縮成形解析装置10におけるモデルファイル記憶部10aに、作成したモデルファイルを記憶させる(ステップ11)。
First, with reference to FIG. 3, the analysis process of a layer shape change is demonstrated.
First, a model file is created as analysis target data in the layer shape change analysis process (step 10).
This model file holds mold shape information and initial shape information of the entire layer.
Then, the created model file is stored in the model file storage unit 10a in the compression molding analysis apparatus 10 (step 11).

次に、層形状変化の解析条件を設定する(ステップ12)。この解析条件は、上述したように、金型及び樹脂の物性,成形速度,解析の時間ステップ等とすることができる。その他、用いる流体解析ソフトウェアにより、例えば各時間ステップにおける全ての節点の速度ベクトルを解析結果として出力するか否か等の解析条件を設定する。   Next, analysis conditions for layer shape change are set (step 12). As described above, the analysis conditions may be physical properties of the mold and resin, molding speed, time step of analysis, and the like. In addition, the analysis conditions such as whether or not to output the velocity vectors of all the nodes at each time step as analysis results are set by the fluid analysis software to be used.

圧縮成形解析装置10における解析条件入力部10bは、この設定された解析条件を圧縮成形解析装置10に入力して記憶させ、層形状変化解析部10cにより利用可能とする(ステップ13)。
そして、層形状変化解析部10cは、モデルファイル記憶部10aに記憶されたモデルファイルを用いて、解析条件入力部10bにより入力された解析条件にもとづき層形状変化の解析処理を実行する(ステップ14)。
最後に、層形状変化解析部10cは、解析結果を解析結果記憶部10dに記憶させる(ステップ15)。
The analysis condition input unit 10b in the compression molding analysis apparatus 10 inputs and stores the set analysis condition in the compression molding analysis apparatus 10 and makes it available to the layer shape change analysis unit 10c (step 13).
Then, the layer shape change analysis unit 10c uses the model file stored in the model file storage unit 10a to execute a layer shape change analysis process based on the analysis conditions input by the analysis condition input unit 10b (step 14). ).
Finally, the layer shape change analysis unit 10c stores the analysis result in the analysis result storage unit 10d (step 15).

次に、図4を参照して、点追跡処理について説明する。
はじめに、解析したい初期中間層形状を任意に決定する(ステップ30)。これは、例えば圧縮成形解析装置10に接続されたディスプレイ等に層全体の初期形状を表示して、この初期形状の内部に任意の中間層形状を画面上で描画して入力することなどにより行うことができる。
Next, the point tracking process will be described with reference to FIG.
First, an initial intermediate layer shape to be analyzed is arbitrarily determined (step 30). For example, the initial shape of the entire layer is displayed on a display or the like connected to the compression molding analysis apparatus 10, and an arbitrary intermediate layer shape is drawn and input on the screen inside the initial shape. be able to.

次に、初期中間層の形状が分かるように、その輪郭上及び/又は内部に点を作成して、これを圧縮成形解析装置10に入力する(ステップ31)。これは、ステップ30で圧縮成形解析装置10に入力された初期中間層形状情報にもとづいて、圧縮成形解析装置10における中間層形状情報入力部10eにより自動的に行うことができる。
次に、圧縮成形解析装置10における点追跡部10fは、解析結果記憶部10dから最初の時間ステップ(層形状変化を開始する時間)における解析結果を読み込む(ステップ32)。以降、当該ステップ32を繰り返し実行する際には、時間ステップごとに解析結果記憶部10dに記憶された解析結果を順次読込む。
Next, a point is created on and / or inside the contour so that the shape of the initial intermediate layer can be understood, and this is input to the compression molding analyzer 10 (step 31). This can be automatically performed by the intermediate layer shape information input unit 10e in the compression molding analyzer 10 based on the initial intermediate layer shape information input to the compression molding analyzer 10 in step 30.
Next, the point tracking unit 10f in the compression molding analysis apparatus 10 reads the analysis result at the first time step (time to start the layer shape change) from the analysis result storage unit 10d (step 32). Thereafter, when the step 32 is repeatedly executed, the analysis results stored in the analysis result storage unit 10d are sequentially read for each time step.

次に、圧縮成形解析装置10における点追跡部10fは、追跡対象の点(追跡点)を決定する(ステップ33)。この追跡点としては、ステップ31において入力された全ての点とすることができる。また、これらの点のうち、一部のもののみを追跡点として決定することもできる。   Next, the point tracking unit 10f in the compression molding analyzer 10 determines a tracking target point (tracking point) (step 33). The tracking points can be all points input in step 31. Also, only some of these points can be determined as tracking points.

なお、ステップ32の処理とステップ33の処理の順番を入れ替えることも可能である。すなわち、まずステップ31において入力された点から追跡点を決定した後に、その決定した追跡点についてのみ解析結果を読み込むようにすることができる。このとき、追跡点に一致する節点については、その節点の解析結果を読み込むことができる。また、節点に一致しない追跡点については、その追跡点の周囲の節点の解析結果を読み込むことができる。
このようにすれば、実際に追跡する点についてのみ解析結果を読み込むことができるため、圧縮成形解析装置10の処理負荷を軽減することができる。このため、メモリに余裕がない場合や、ステップ31において入力された点の数が多い場合などには、このようにすることも好ましい。
Note that the order of the process of step 32 and the process of step 33 can be switched. That is, first, after the tracking point is determined from the point input in step 31, the analysis result can be read only for the determined tracking point. At this time, for a node that matches the tracking point, the analysis result of the node can be read. For tracking points that do not match the nodes, the analysis results of the nodes around the tracking points can be read.
In this way, the analysis result can be read only for the points that are actually tracked, so that the processing load of the compression molding analysis apparatus 10 can be reduced. For this reason, it is also preferable to do this when there is not enough memory or when the number of points input in step 31 is large.

そして、点追跡部10fは、追跡対象の点の座標情報を取得する(ステップ34)。この座標情報は、ステップ31において入力された点の座標情報を用いることができる。
次に、点追跡部10fは、追跡点における速度ベクトルを算出する(ステップ35)。この速度ベクトルは、追跡点が、解析結果に含まれる節点のうちの一つに一致する場合には、その一致する節点の物理量として解析結果に保有されている速度ベクトルをそのまま用いることができる。また、追跡点が節点に一致しない場合には、その追跡点の周囲の節点の速度ベクトルから補間計算により算出する。
And the point tracking part 10f acquires the coordinate information of the point of tracking object (step 34). As this coordinate information, the coordinate information of the point input in step 31 can be used.
Next, the point tracking unit 10f calculates a velocity vector at the tracking point (step 35). When the tracking point matches one of the nodes included in the analysis result, the velocity vector held in the analysis result can be used as the physical quantity of the matching node as it is. Further, when the tracking point does not coincide with the node, it is calculated by interpolation calculation from the velocity vectors of the nodes around the tracking point.

この補間計算の方法としては、従来公知の種々の方法を用いることができるが、例えば、形状関数を用いたアイソパラメトリック要素を使用した有限要素法等を用いることができる。
アイソパラメトリック要素を使用した有限要素法の具体例としては、例えば、4節点からなる二次元要素の内部の物理量を4節点の物理量から形状関数により内挿して求める方法を挙げることができる。勿論、より多くの節点を用いて補間計算を行うことも可能である。
As the interpolation calculation method, various conventionally known methods can be used. For example, a finite element method using an isoparametric element using a shape function can be used.
As a specific example of the finite element method using an isoparametric element, for example, a method in which a physical quantity inside a two-dimensional element having four nodes is interpolated from a physical quantity of four nodes by a shape function can be cited. Of course, it is also possible to perform interpolation calculation using more nodes.

次に、点追跡部10fは、追跡点の座標情報及び速度ベクトルとから次の時間ステップにおける点の位置の座標情報を算出し、これを時間ステップごとに点追跡結果記憶部10gに記憶させる(ステップ36)。
そして、全ての追跡点について、ステップ34〜ステップ36の処理を繰り返す(ステップ37のNo)。
Next, the point tracking unit 10f calculates the coordinate information of the position of the point in the next time step from the coordinate information and the velocity vector of the tracking point, and stores this in the point tracking result storage unit 10g for each time step ( Step 36).
And the process of step 34-step 36 is repeated about all the tracking points (No of step 37).

全ての追跡点について、次の時間ステップにおける点の位置の座標情報の算出と、その座標情報の点追跡結果記憶部10gへの記憶が終了すると(ステップ37のYes)、次に、解析終了時間になったかどうかを判定する(ステップ38)。
そして、解析終了時間に達していなければ(ステップ38のNo)、次の時間ステップにステップ32からステップ37の処理を繰り返し、解析終了時間に達していれば(ステップ38のYes)、追跡点処理を終了する。
When the calculation of the coordinate information of the position of the point in the next time step and the storage of the coordinate information in the point tracking result storage unit 10g for all the tracking points is completed (Yes in step 37), the analysis end time is then reached. Is determined (step 38).
If the analysis end time has not been reached (No in step 38), the processing from step 32 to step 37 is repeated for the next time step. If the analysis end time has been reached (Yes in step 38), the tracking point processing is performed. Exit.

以上説明したように、本実施形態の圧縮成形解析装置10によれば、層全体の初期形状における中間層の位置に点を配置して、これを追跡することができる。
このため、あらゆる初期中間層形状を追跡することができ、従来は困難であった少量の中間層の形状変化を追跡することが可能となる。
また、本実施形態の圧縮成形解析装置10によれば、層形状変化の解析を一度行うのみで、様々な形状の中間層形状について点追跡を行うことができる。
このため、従来のアダプティブメッシュ法のみによる成形解析方法に比べて、より効率的な成形解析を行うことが可能となる。
As described above, according to the compression molding analysis apparatus 10 of this embodiment, it is possible to place a point at the position of the intermediate layer in the initial shape of the entire layer and track this.
Therefore, any initial intermediate layer shape can be tracked, and a small amount of intermediate layer shape change, which has been difficult in the past, can be tracked.
In addition, according to the compression molding analysis apparatus 10 of the present embodiment, it is possible to perform point tracking for various shapes of the intermediate layer shape by performing the analysis of the layer shape change only once.
For this reason, it becomes possible to perform more efficient forming analysis as compared with the forming analysis method based only on the conventional adaptive mesh method.

[第二実施形態]
次に、本発明の第二実施形態について、図5〜図7を参照して説明する。図5は、本実施形態の圧縮成形解析装置の構成を示すブロック図である。図6は、本実施形態の大域的点追跡法における点追跡の様子を示す図である。図7は、本実施形態の大域的点追跡法における初期形状上の中間層分布を算出する様子を示す図である。
[Second Embodiment]
Next, a second embodiment of the present invention will be described with reference to FIGS. FIG. 5 is a block diagram showing the configuration of the compression molding analysis apparatus of the present embodiment. FIG. 6 is a diagram showing a state of point tracking in the global point tracking method of the present embodiment. FIG. 7 is a diagram showing how the intermediate layer distribution on the initial shape is calculated in the global point tracking method of the present embodiment.

本実施形態は、第一実施形態における点追跡法をさらに発展させたものであり、層全体の形状変化の解析にあたり、全ての節点について、その初期座標値を保有させたまま、点追跡を実行し、最終形状における初期座標値分布を得ることで、成形後の中間層の分布が最適となるような初期形状における中間層分布を算出することの可能な大域的点追跡法に関するものである。   This embodiment is a further development of the point tracking method in the first embodiment. In analyzing the shape change of the entire layer, point tracking is executed while retaining the initial coordinate values of all nodes. In addition, the present invention relates to a global point tracking method capable of calculating the intermediate layer distribution in the initial shape such that the distribution of the intermediate layer after molding is optimal by obtaining the initial coordinate value distribution in the final shape.

本実施形態の圧縮成形解析装置10は、図5に示すように、第一実施形態における構成に加え、さらに、大域的点追跡部10h,大域的点追跡結果記憶部10i,成形後中間層形状情報入力部10j,初期中間層形状情報生成部10k,初期中間層形状情報記憶部10lを備えている。   As shown in FIG. 5, the compression molding analysis apparatus 10 of the present embodiment has a global point tracking unit 10 h, a global point tracking result storage unit 10 i, a post-molding intermediate layer shape, in addition to the configuration in the first embodiment. An information input unit 10j, an initial intermediate layer shape information generation unit 10k, and an initial intermediate layer shape information storage unit 10l are provided.

大域的点追跡部10hは、大域的点追跡処理を実行するものであり、所定の流体解析ソフトウェアを実行することで、圧縮成形解析装置10におけるCPU及びメモリ上に構成される。
この流体解析ソフトウェアとしては、層形状変化解析部10cと同様に、例えばFLUENT社から販売された「POLYFLOW」を用いることができる。
The global point tracking unit 10h executes a global point tracking process, and is configured on a CPU and a memory in the compression molding analysis apparatus 10 by executing predetermined fluid analysis software.
As this fluid analysis software, for example, “POLYFLOW” sold by FLUENT can be used as in the layer shape change analysis unit 10c.

大域的点追跡部10hは、所定の時間ステップごとに、解析結果記憶部10dから解析結果を読み込み、すべての節点の初期座標値を、時間ステップ及び節点ごとに大域的点追跡結果記憶部10iに記憶させる。
なお、節点の初期座標値としては、変形に際してアダプティブメッシュ法が用いられなかった時間ステップにおいては、移動前の初期座標値と同じものを用いることができる。一方、アダプティブメッシュ法が用いられた時間ステップにおいては、補間計算によって、リメッシングにより新たに設定された節点の初期座標値を算出する。したがって、初期座標値を大域的点追跡結果記憶部10iに一度記憶させた後は、この初期座標値をそのまま用いるか、あるいは更新して用いることができる。
The global point tracking unit 10h reads the analysis result from the analysis result storage unit 10d at every predetermined time step, and stores the initial coordinate values of all the nodes in the global point tracking result storage unit 10i for each time step and node. Remember.
As the initial coordinate values of the nodes, the same initial coordinate values as before the movement can be used in the time step when the adaptive mesh method is not used in the deformation. On the other hand, in a time step in which the adaptive mesh method is used, initial coordinate values of nodes newly set by remeshing are calculated by interpolation calculation. Therefore, after the initial coordinate value is stored once in the global point tracking result storage unit 10i, the initial coordinate value can be used as it is or can be updated and used.

また、大域的点追跡部10hは、層全体の初期形状におけるすべての節点について、点追跡部10fにより、所定の時間ステップ分ずつ点追跡を行わせ、その追跡結果を大域的点追跡結果記憶部10iに記憶させる。この追跡結果には、時間ステップごとに、追跡対象とした全ての点の移動後の座標情報と速度ベクトルなどの物理量が含まれている。
なお、層形状変化の解析結果に、座標情報と物理量が記憶されていない部分における点(節点でない点)については、補間計算を行うことで速度ベクトルを算出し、点追跡することができる。この補間計算の方法については、第一実施形態の処理手順のステップ35において上述した方法と同様の方法とすることができる。
The global point tracking unit 10h causes the point tracking unit 10f to perform point tracking for each predetermined time step for all nodes in the initial shape of the entire layer, and the tracking result is stored in the global point tracking result storage unit. 10i is stored. This tracking result includes, for each time step, physical information such as coordinate information after moving all the points to be tracked and a velocity vector.
It should be noted that with respect to a point (a point that is not a node) in a portion where coordinate information and a physical quantity are not stored in the analysis result of the layer shape change, a velocity vector can be calculated and point tracking can be performed by performing interpolation calculation. The interpolation calculation method can be the same method as described above in step 35 of the processing procedure of the first embodiment.

そして、大域的点追跡部10hは、所定の時間ステップごとに、追跡した点の初期座標値分布を生成する。
この際、追跡点のすべてが節点に一致する場合は、その節点の初期座標値分布がそのまま追跡点の初期座標値分布として生成される。また、追跡点の一部もしくはすべてが節点に一致しない場合は、その一致しない追跡点について各追跡点の周囲の節点の初期座標値にもとづき補間計算により、それぞれの追跡点の初期座標値を算出して、初期座標値分布が生成される。
これによって、大域的点追跡部10hは、解析終了時間まで上記処理を繰り返すことで、層全体の最終形状における初期座標値分布を生成し、これを大域的点追跡結果記憶部10iに記憶させることができる。
Then, the global point tracking unit 10h generates an initial coordinate value distribution of the tracked points every predetermined time step.
At this time, if all of the tracking points coincide with the nodes, the initial coordinate value distribution of the nodes is generated as it is as the initial coordinate value distribution of the tracking points. If some or all of the tracking points do not match the nodes, the initial coordinate value of each tracking point is calculated by interpolation calculation based on the initial coordinate values of the nodes around each tracking point. Thus, an initial coordinate value distribution is generated.
Thereby, the global point tracking unit 10h generates the initial coordinate value distribution in the final shape of the entire layer by repeating the above processing until the analysis end time, and stores this in the global point tracking result storage unit 10i. Can do.

図6は、このような大域的点追跡法における点追跡の様子を示す図である。本実施形態の圧縮成形解析装置10による大域的点追跡法によれば、同図に示すように、初期形状(a)における層全体に点を配置し(b)、層を圧縮する工程において各点を点追跡して(c)、各点の初期座標値を大域的点追跡結果記憶部10iに記憶させておくことで、成形後の追跡点の初期座標値分布を取得することが可能となる。   FIG. 6 is a diagram showing a state of point tracking in such a global point tracking method. According to the global point tracking method by the compression molding analysis apparatus 10 of the present embodiment, as shown in FIG. By tracking the points (c) and storing the initial coordinate values of each point in the global point tracking result storage unit 10i, it is possible to obtain the initial coordinate value distribution of the tracking points after molding. Become.

大域的点追跡結果記憶部10iは、各追跡点の初期座標値と、点追跡部10fにより計算された追跡結果を記憶する。
この追跡結果には、層全体を初期形状から最終形状まで成形する過程における、時間ステップごとの、層全体の初期形状の節点から追跡したすべての点の座標情報、及び各点における速度ベクトル等の物理量が含まれる。
The global point tracking result storage unit 10i stores the initial coordinate value of each tracking point and the tracking result calculated by the point tracking unit 10f.
The tracking results include the coordinate information of all points tracked from the nodes of the initial shape of the entire layer and the velocity vector at each point in the process of forming the entire layer from the initial shape to the final shape. Includes physical quantities.

成形後中間層形状情報入力部10jは、層全体の最終形状上に決定された任意の領域上の点の座標情報により構成される成形後中間層形状情報を圧縮成形解析装置10に入力し、初期中間層形状情報生成部10kで使用するために、圧縮成形解析装置10に記憶させる。層全体の最終形状とは、層形状変化の解析により解析結果記憶部10dに記憶された層全体の最終形状情報を意味する。
任意の領域として、圧縮成形により得ることを希望する理想的な中間層の形状を決定することで、初期中間層形状情報生成部10kにより、その理想的な中間層を成形するために必要な、層全体の初期形状における中間層の形状を計算することが可能となる。
The post-molding intermediate layer shape information input unit 10j inputs post-molding intermediate layer shape information constituted by coordinate information of points on an arbitrary region determined on the final shape of the entire layer to the compression molding analysis apparatus 10, For use in the initial intermediate layer shape information generation unit 10k, the compression molding analysis apparatus 10 stores the information. The final shape of the entire layer means the final shape information of the entire layer stored in the analysis result storage unit 10d by analyzing the layer shape change.
By determining the ideal intermediate layer shape desired to be obtained by compression molding as an arbitrary region, the initial intermediate layer shape information generation unit 10k is necessary to form the ideal intermediate layer. It is possible to calculate the shape of the intermediate layer in the initial shape of the entire layer.

この任意の領域上の点としては、その任意の領域の輪郭上及び内部に、任意の領域を判別することができる程度のものを特定する。その特定方法としては、例えば領域に含まれる全ての節点を特定する方法とすることができる。また、領域の全周をn分割した点などを特定する方法などとすることもできる。なお、この特定については、コンピュータにより自動的に行う他、人手により行うこともできる。   As the points on the arbitrary region, those on the outline and inside of the arbitrary region that can discriminate the arbitrary region are specified. As the specifying method, for example, it is possible to specify all the nodes included in the region. Also, a method of identifying points obtained by dividing the entire circumference of the region into n may be used. Note that this specification can be performed manually as well as automatically by a computer.

初期中間層形状情報生成部10kは、層全体の最終形状上に決定された任意の領域上の点の初期座標値にもとづいて、当該任意の領域に対応する、層全体の初期形状における初期中間層形状情報を生成する。
図7は、この初期中間層形状情報生成部10kによって、層全体の初期形状上に初期中間層分布を描いた様子を示す図である。
同図に示すように、層全体の最終形状上に理想とする中間層の形状を決定して、この形状に対応する点を特定し、特定した点の初期座標値を取得して、この初期座標値の点を層全体の初期形状情報に重ねて配置することで、成形後に理想的な中間層を構成する初期中間層形状情報を生成することができる。
The initial intermediate layer shape information generation unit 10k generates an initial intermediate layer in the initial shape of the entire layer corresponding to the arbitrary region based on the initial coordinate value of the point on the arbitrary region determined on the final shape of the entire layer. Layer shape information is generated.
FIG. 7 is a diagram illustrating a state in which the initial intermediate layer distribution is drawn on the initial shape of the entire layer by the initial intermediate layer shape information generation unit 10k.
As shown in the figure, the ideal shape of the intermediate layer is determined on the final shape of the entire layer, the point corresponding to this shape is identified, the initial coordinate value of the identified point is obtained, and this initial value is obtained. By arranging the coordinate value points so as to overlap the initial shape information of the entire layer, it is possible to generate initial intermediate layer shape information constituting an ideal intermediate layer after molding.

初期中間層形状情報記憶部10lは、初期中間層形状情報生成部10kにより生成された初期中間層形状情報を記憶する。この初期中間層形状情報は、層全体の初期形状を表す形状情報、及びこの初期形状における理想の中間層を表す点の座標情報(初期座標値)を含むものである。   The initial intermediate layer shape information storage unit 101 stores the initial intermediate layer shape information generated by the initial intermediate layer shape information generation unit 10k. The initial intermediate layer shape information includes shape information indicating the initial shape of the entire layer, and coordinate information (initial coordinate values) of points indicating the ideal intermediate layer in the initial shape.

次に、本実施形態の圧縮成形解析装置10における処理手順について、図8及び図9を参照して説明する。図8は、本実施形態の圧縮成形解析装置により実行される大域的点追跡処理を示すフローチャートであり、図9は、同装置により実行される初期中間層分布生成処理を示すフローチャートである。   Next, a processing procedure in the compression molding analysis apparatus 10 of the present embodiment will be described with reference to FIGS. FIG. 8 is a flowchart showing a global point tracking process executed by the compression molding analysis apparatus of the present embodiment, and FIG. 9 is a flowchart showing an initial intermediate layer distribution generation process executed by the apparatus.

まず、本実施形態の大域的点追跡処理を実行するにあたっては、最初に図3を用いて上述した層形状変化の解析処理を実行する。そして、次に、図8に示す大域的点追跡処理を実行する。
図8に示すように、圧縮成形解析装置10における大域的点追跡部10hは、解析結果記憶部10dから層形状変化の解析結果を読み込む(ステップ50)。
次に、大域的点追跡部10hは、読み込んだ解析結果におけるすべての節点の初期座標値を大域的点追跡結果記憶部10iに記憶させる(ステップ51)。
First, when executing the global point tracking process of the present embodiment, the layer shape change analysis process described above with reference to FIG. 3 is first executed. Then, the global point tracking process shown in FIG. 8 is executed.
As shown in FIG. 8, the global point tracking unit 10h in the compression molding analysis apparatus 10 reads the analysis result of the layer shape change from the analysis result storage unit 10d (step 50).
Next, the global point tracking unit 10h stores the initial coordinate values of all nodes in the read analysis result in the global point tracking result storage unit 10i (step 51).

このとき、大域的点追跡部10hは、時間ステップごとにこれらの処理を繰り返し実行する。
すなわち、第一回目の処理において、大域的点追跡部10hは、層全体の初期形状におけるすべての節点の座標値、及びそれぞれの節点における物理量を読み込む。
また、第二回目以降の処理において、大域的点追跡部10hは、次の時間ステップの形状におけるすべての節点の座標情報、及びそれぞれの節点における物理量を読み込む。
At this time, the global point tracking unit 10h repeatedly executes these processes for each time step.
That is, in the first process, the global point tracking unit 10h reads the coordinate values of all the nodes in the initial shape of the entire layer, and the physical quantities at the respective nodes.
In the second and subsequent processes, the global point tracking unit 10h reads the coordinate information of all the nodes in the next time step shape and the physical quantities at the respective nodes.

さらに、大域的点追跡部10hは、第一回目の処理において、層全体の初期形状における各節点の座標情報を初期座標値として大域的点追跡結果記憶部10iに記憶させる。
また、第二回目以降の処理において、大域的点追跡部10hは、前回の処理で算出した次の時間ステップの追跡点の座標情報(ステップ52−4において後述する)と、解析結果記憶部10dから読み込んだ、次の時間ステップの層全体の形状におけるすべての節点の座標情報とを比較する。
Further, in the first processing, the global point tracking unit 10h stores the coordinate information of each node in the initial shape of the entire layer in the global point tracking result storage unit 10i as an initial coordinate value.
In the second and subsequent processes, the global point tracking unit 10h includes the coordinate information of the tracking point of the next time step calculated in the previous process (described later in Step 52-4), and the analysis result storage unit 10d. Compare the coordinate information of all nodes in the shape of the whole layer of the next time step read from.

そして、これらの節点が一致する場合には、前回記憶した初期座標値分布を、そのまま次の時間ステップの層全体の形状における節点の初期座標値分布として、大域的点追跡結果記憶部10iに記憶させる。
また、これらの節点が一致しない場合、すなわち、層形状変化の解析において、リメッシングが行われた時間ステップでは、補間計算により、前回記憶した初期座標値分布にもとづいて、次の時間ステップの層全体の形状における節点の初期座標値分布を生成する。
この初期座標値分布の補間計算方法は、ステップ35の説明において上述した方法と同様の方法とすることができる。
If these nodes coincide, the initial coordinate value distribution stored last time is stored in the global point tracking result storage unit 10i as the initial coordinate value distribution of the nodes in the shape of the entire layer at the next time step as it is. Let
Also, if these nodes do not match, that is, in the time step when remeshing is performed in the analysis of the layer shape change, the entire layer of the next time step is calculated based on the previously stored initial coordinate value distribution by interpolation calculation. An initial coordinate value distribution of nodes in the shape of is generated.
The interpolation calculation method of the initial coordinate value distribution can be the same method as described above in the description of step 35.

次に、大域的点追跡部10hは、点追跡部10fに、第一実施形態における点追跡と同様の処理を実行させる(ステップ52)。
まず、点追跡部10fは、追跡点を決定する(ステップ52−1)。本実施形態の場合、この追跡点は、第一回目の処理(最初の時間ステップにおける処理)においては、ステップ50で読み込まれた全ての節点とすることができる。第二回目以降の処理(第二の時間ステップ以降における処理)においては、前回追跡した点の移動後の点となる。
そして、点追跡部10fは、追跡点の座標情報を取得する(ステップ52−2)。この座標情報は、第一回目の処理においては、ステップ50で読み込まれた節点の座標情報を用いることができる。第二回目以降の処理においては、前回追跡した点の移動後の点の座標情報を用いることができる。
Next, the global point tracking unit 10h causes the point tracking unit 10f to execute the same processing as the point tracking in the first embodiment (step 52).
First, the point tracking unit 10f determines a tracking point (step 52-1). In the case of the present embodiment, this tracking point can be all nodes read in step 50 in the first processing (processing in the first time step). In the second and subsequent processes (processes after the second time step), the point after the last tracked point is moved.
And the point tracking part 10f acquires the coordinate information of a tracking point (step 52-2). As this coordinate information, the coordinate information of the node read in step 50 can be used in the first processing. In the second and subsequent processing, the coordinate information of the point after the movement of the point tracked last time can be used.

次に、点追跡部10fは、追跡点における速度ベクトルを算出する(ステップ52−3)。この速度ベクトルは、第一回目の処理においては、ステップ50で読み込まれた節点の物理量に含まれるものを用いることができる。第二回目以降の処理においては、追跡点が、解析結果に含まれる節点のうちの一つに一致する場合には、その一致する節点の物理量として解析結果に保有されている速度ベクトルをそのまま用いることができる。また、追跡点が節点に一致しない場合には、その追跡点の周囲の節点の速度ベクトルから補間計算により算出する。
この補間計算の方法としては、ステップ35の説明において上述した方法と同様の方法とすることができる。
Next, the point tracking unit 10f calculates a velocity vector at the tracking point (step 52-3). As the velocity vector, a vector included in the physical quantity of the node read in step 50 can be used in the first processing. In the second and subsequent processing, when the tracking point matches one of the nodes included in the analysis result, the velocity vector held in the analysis result is used as it is as the physical quantity of the matching node. be able to. Further, when the tracking point does not coincide with the node, it is calculated by interpolation calculation from the velocity vectors of the nodes around the tracking point.
As the interpolation calculation method, the same method as described above in the description of step 35 can be used.

次に、点追跡部10fは、追跡点の座標情報と速度ベクトルとから次の時間ステップにおける位置の座標を算出して、大域的点追跡部10hに受け渡し、大域的点追跡部10hは、この座標情報を時間ステップごとに、大域的点追跡結果記憶部10iに記憶させる(ステップ52−4)。
そして、全ての追跡点について、ステップ52−2〜ステップ52−4の処理を繰り返す(ステップ52−5のNo)。
Next, the point tracking unit 10f calculates the coordinates of the position at the next time step from the coordinate information of the tracking point and the velocity vector, and passes the coordinates to the global point tracking unit 10h. The global point tracking unit 10h The coordinate information is stored in the global point tracking result storage unit 10i for each time step (step 52-4).
And the process of step 52-2-step 52-4 is repeated about all the tracking points (No of step 52-5).

全ての追跡点について、次の時間ステップにおける位置の座標の算出と、その座標情報の大域的点追跡結果記憶部10iへの記憶が終了すると(ステップ52−5のYes)、大域的点追跡部10hは、時間ステップごとに、これらの追跡点の初期座標値分布を生成する(ステップ53)。   When calculation of the coordinates of the position in the next time step and storage of the coordinate information in the global point tracking result storage unit 10i for all the tracking points is completed (Yes in step 52-5), the global point tracking unit 10h generates an initial coordinate value distribution of these tracking points for each time step (step 53).

このとき、各追跡点がステップ50で読み込まれた解析結果における節点に一致する場合は、ステップ51で記憶した各節点の初期座標値をそのまま用いて、追跡点の初期座標値分布を生成する。また、各追跡点がステップ50で読み込まれた解析結果における節点に一致しない場合は、ステップ51で記憶した各節点の初期座標値にもとづいて、補間計算により各追跡点の初期座標値を算出し、算出した初期座標値を用いて、追跡点の初期座標値分布を生成する。
なお、初期座標値分布とは、各追跡点の座標情報、及びそれぞれの初期座標値からなる情報を意味する。
At this time, if each tracking point matches the node in the analysis result read in step 50, the initial coordinate value distribution of the tracking point is generated using the initial coordinate value of each node stored in step 51 as it is. If each tracking point does not coincide with the node in the analysis result read in step 50, the initial coordinate value of each tracking point is calculated by interpolation based on the initial coordinate value of each node stored in step 51. The initial coordinate value distribution of the tracking points is generated using the calculated initial coordinate values.
Note that the initial coordinate value distribution means coordinate information of each tracking point and information including each initial coordinate value.

次に、大域的点追跡部10hは、解析終了時間になったかどうかを判定する(ステップ54)。
そして、解析終了時間に達していなければ(ステップ54のNo)、次の時間ステップにステップ50からステップ53の処理を繰り返し、解析終了時間に達していれば(ステップ54のYes)、大域的追跡点処理を終了する。
Next, the global point tracking unit 10h determines whether or not the analysis end time has come (step 54).
If the analysis end time has not been reached (No in step 54), the processing from step 50 to step 53 is repeated in the next time step. If the analysis end time has been reached (Yes in step 54), global tracking is performed. End point processing.

次に、図9を参照して、初期中間層分布生成処理について説明する。
まず、初期中間層分布生成処理に先立って、解析結果記憶部10dに記憶されている層全体の最終形状上に、任意の領域(以下、最終中間層領域と称する場合がある。)を決定する(ステップ70)。この層全体の最終形状は、解析結果記憶部10dに記憶されている層全体の最終形状の形状情報により把握することができる。
最終中間層領域としては、希望する理想的な成形後の中間層の領域を決定することが好ましい。そうすれば、初期中間層分布生成処理により、理想的な成形後の中間層を得るために必要な、初期形状における中間層分布を算出することができる。この最終中間層領域の決定は、人により行われる。
Next, the initial intermediate layer distribution generation process will be described with reference to FIG.
First, prior to the initial intermediate layer distribution generation process, an arbitrary region (hereinafter sometimes referred to as a final intermediate layer region) is determined on the final shape of the entire layer stored in the analysis result storage unit 10d. (Step 70). The final shape of the entire layer can be grasped from the shape information of the final shape of the entire layer stored in the analysis result storage unit 10d.
As the final intermediate layer region, it is preferable to determine a desired ideal intermediate layer region after molding. Then, the intermediate layer distribution in the initial shape necessary for obtaining the ideal intermediate layer after molding can be calculated by the initial intermediate layer distribution generation process. This final intermediate layer region is determined by a person.

次に、圧縮成形解析装置10における成形後中間層形状情報入力部10jは、最終中間層領域上に、その領域を判別することができる程度に特定された点の座標情報を入力する(ステップ71)。
そして、圧縮成形解析装置10における初期中間層形状情報生成部10kは、最終中間層領域上の点の初期座標値を大域的点追跡結果記憶部10iに記憶された情報にもとづき算出する(ステップ72)。
Next, the post-molding intermediate layer shape information input unit 10j in the compression molding analysis apparatus 10 inputs the coordinate information of the points specified to such an extent that the region can be discriminated on the final intermediate layer region (step 71). ).
Then, the initial intermediate layer shape information generation unit 10k in the compression molding analysis apparatus 10 calculates the initial coordinate values of the points on the final intermediate layer region based on the information stored in the global point tracking result storage unit 10i (step 72). ).

このとき、最終中間層領域上の点の座標情報が、大域的点追跡結果記憶部10iに記憶されている場合には、初期中間層形状情報生成部10kは、その最終中間層領域上の点に対応する初期座標値を大域的点追跡結果記憶部10iから読み込むことで、初期座標値を取得することができる。
また、最終中間層領域上の点の座標情報が、大域的点追跡結果記憶部10iに記憶されていない場合には、初期中間層形状情報生成部10kは、その最終中間層領域上の点を取り囲む周囲の点として、大域的点追跡結果記憶部10iに記憶されている点に対応する初期座標値を大域的点追跡結果記憶部10iから読み込んで、上述したような補間計算を行うことにより、最終中間層領域上の点の初期座標値を算出することができる。
At this time, when the coordinate information of the point on the final intermediate layer area is stored in the global point tracking result storage unit 10i, the initial intermediate layer shape information generation unit 10k By reading the initial coordinate values corresponding to the global point tracking result storage unit 10i, the initial coordinate values can be acquired.
When the coordinate information of the point on the final intermediate layer area is not stored in the global point tracking result storage unit 10i, the initial intermediate layer shape information generation unit 10k selects the point on the final intermediate layer area. By reading the initial coordinate values corresponding to the points stored in the global point tracking result storage unit 10i as the surrounding points from the global point tracking result storage unit 10i and performing the interpolation calculation as described above, An initial coordinate value of a point on the final intermediate layer region can be calculated.

次に、初期中間層形状情報生成部10kは、解析結果記憶部10dから層全体の初期形状情報を読み込んで、この初期形状情報に最終中間層領域に対応する点の初期座標値の点を重ねて配置する(ステップ73)。これによって、最終中間層領域に対応する中間層の初期形状情報を生成することができる。
そして、初期中間層形状情報生成部10kは、得られた中間層の初期形状情報を初期中間層形状情報記憶部10lに記憶させる(ステップ74)。
Next, the initial intermediate layer shape information generation unit 10k reads the initial shape information of the entire layer from the analysis result storage unit 10d, and superimposes the initial coordinate value points of the points corresponding to the final intermediate layer region on this initial shape information. (Step 73). Thereby, the initial shape information of the intermediate layer corresponding to the final intermediate layer region can be generated.
Then, the initial intermediate layer shape information generation unit 10k stores the obtained initial intermediate layer shape information in the initial intermediate layer shape information storage unit 10l (step 74).

以上説明したように、本実施形態の圧縮成形解析装置によれば、成形後の任意の形状の中間層に対応する、層全体の初期形状における中間層の初期形状情報を生成することができ、中間層分布を算出することができる。
このため、少ない量の中間層についても、成形後に最適な中間層を得ることが可能な、層全体の初期形状における理想的な中間層の分布を把握することが可能となる。
As described above, according to the compression molding analysis apparatus of the present embodiment, it is possible to generate the initial shape information of the intermediate layer in the initial shape of the entire layer corresponding to the intermediate layer of any shape after molding, An intermediate layer distribution can be calculated.
For this reason, even with a small amount of the intermediate layer, it is possible to obtain an ideal intermediate layer distribution in the initial shape of the entire layer, which can obtain an optimal intermediate layer after molding.

[第三実施形態]
次に、本発明の第三実施形態について、図10を参照して説明する。同図は、本実施形態の圧縮成形解析装置により実行される物性の違いを考慮した層形状変化の解析処理を示すフローチャートである。
本実施形態は、第一実施形態における点追跡法をさらに発展させたものであり、複数の材料からなる層全体の形状変化の解析にあたり、材料ごとの物性の違いを考慮して、形状変化の解析を行うものである。その他の点については、第一実施形態と同様である。
[Third embodiment]
Next, a third embodiment of the present invention will be described with reference to FIG. This figure is a flowchart showing a layer shape change analysis process taking into account the difference in physical properties executed by the compression molding analysis apparatus of the present embodiment.
This embodiment is a further development of the point tracking method in the first embodiment.In analyzing the shape change of the entire layer composed of a plurality of materials, the shape change is considered in consideration of the difference in physical properties for each material. Analyze. About another point, it is the same as that of 1st embodiment.

まず、層形状変化の解析処理における解析対象データとしてのモデルファイルを作成する(ステップ90)。
このモデルファイルは、第一実施形態におけるものと同様に、金型形状、及び初期樹脂領域形状の情報を保有する。
そして、圧縮成形解析装置10におけるモデルファイル記憶部10aに、作成したモデルファイルを記憶させる(ステップ91)。
First, a model file is created as analysis target data in the layer shape change analysis process (step 90).
This model file holds information on the mold shape and the initial resin region shape as in the first embodiment.
Then, the created model file is stored in the model file storage unit 10a in the compression molding analysis apparatus 10 (step 91).

次に、層を構成する材料ごとに、その材料を配置する領域を決定する(ステップ92)。これは、中間層の初期形状を決定することに相当する。
そして、各領域に仮想スカラー量を設定する(ステップ93)。具体的には、ある領域に含まれる節点とそうでない節点にそれぞれ別の値のスカラー量(例えば、中間層は1、外層は2など)を割り当てる。
このスカラー量は、層全体が初期形状から最終形状に至るまでの全過程において、時間ステップごとに、節点に対応づけて解析結果記憶部10dに記憶される。
これによって、解析結果記憶部10dに記憶された時間ステップごとの各節点のスカラー量にもとづいて、時間ステップごとのスカラー量分布を得ることが可能になる。
Next, for each material constituting the layer, a region in which the material is arranged is determined (step 92). This corresponds to determining the initial shape of the intermediate layer.
Then, a virtual scalar quantity is set for each area (step 93). Specifically, different scalar amounts (for example, 1 for the intermediate layer, 2 for the outer layer, etc.) are assigned to nodes included in a certain region and nodes that are not.
This scalar amount is stored in the analysis result storage unit 10d in association with the node for each time step in the entire process from the initial shape to the final shape of the entire layer.
Accordingly, it is possible to obtain a scalar amount distribution for each time step based on the scalar amount of each node for each time step stored in the analysis result storage unit 10d.

次に、仮想スカラー量と物性との関連づけを行う(ステップ94)。例えば、物理量としての粘度について、スカラー量1の部分は粘度3000、スカラー量2の部分は粘度6000というように関連づける。
また、その他の解析条件、例えば金型及び樹脂の物性,成形速度,解析の時間ステップ等を設定する(ステップ95)。
そして、圧縮成形解析装置10における解析条件入力部10bは、ステップ92〜ステップ95により設定された解析条件を圧縮成形解析装置10に入力して記憶させ、層形状変化解析部10cにより利用可能とする(ステップ96)。
Next, the virtual scalar quantity and physical properties are associated (step 94). For example, regarding the viscosity as a physical quantity, the portion of the scalar amount 1 is related to the viscosity of 3000, and the portion of the scalar amount 2 is related to the viscosity of 6000.
Further, other analysis conditions, such as physical properties of the mold and resin, molding speed, time step of analysis, etc. are set (step 95).
Then, the analysis condition input unit 10b in the compression molding analysis apparatus 10 inputs and stores the analysis conditions set in steps 92 to 95 in the compression molding analysis apparatus 10 and makes them usable by the layer shape change analysis unit 10c. (Step 96).

次に、層形状変化解析部10cは、第一実施形態における場合と同様に、モデルファイル記憶部10aに記憶されたモデルファイルを用いて、解析条件入力部10bにより入力された解析条件にもとづき層形状変化の解析処理を実行する(ステップ97)。
そして、層形状変化解析部10cは、解析結果を解析結果記憶部10dに記憶させる(ステップ98)。
Next, as in the first embodiment, the layer shape change analysis unit 10c uses the model file stored in the model file storage unit 10a to create a layer based on the analysis conditions input by the analysis condition input unit 10b. A shape change analysis process is executed (step 97).
Then, the layer shape change analysis unit 10c stores the analysis result in the analysis result storage unit 10d (step 98).

以上説明したように、本実施形態の圧縮成形解析装置によれば、物性の違いを考慮して、層形状変化の解析を行うことができ、点追跡を行うことが可能となる。   As described above, according to the compression molding analysis apparatus of the present embodiment, the layer shape change can be analyzed in consideration of the difference in physical properties, and the point tracking can be performed.

本発明は、以上の実施形態に限定されるものではなく、本発明の範囲内において、種々の変更実施が可能であることは言うまでもない。
例えば、上記実施形態では、圧縮成形解析について説明したが、本発明は射出成形解析に適用することも可能である。
また、上記実施形態では、層の形状変化を二次元でシミュレーションする場合について説明したが、同様の方法を三次元に適用することが可能である。
It goes without saying that the present invention is not limited to the above embodiment, and that various modifications can be made within the scope of the present invention.
For example, although the compression molding analysis has been described in the above embodiment, the present invention can be applied to the injection molding analysis.
Moreover, although the said embodiment demonstrated the case where the shape change of a layer was simulated in two dimensions, it is possible to apply the same method to three dimensions.

さらに、上記実施形態では、層形状変化解析部10cによって、まず層全体の形状変化を初期形状から最終形状まで解析した後に、点追跡又は大域的点追跡を行っているが、層全体を初期形状から最終形状まで成形する過程における解析の時間ステップごとに、逐次層形状変化解析と点追跡又は大域的点追跡を行うようにすることも可能である。
また、上記第二実施形態と第三実施形態とを組み合わせて、層を構成する材料の物性の違いを考慮して、大域的点追跡を行うようにするなど適宜変更することが可能である。
Furthermore, in the above embodiment, the layer shape change analysis unit 10c first analyzes the shape change of the entire layer from the initial shape to the final shape, and then performs point tracking or global point tracking. It is also possible to perform sequential layer shape change analysis and point tracking or global point tracking at each time step of analysis in the process from forming to final shape.
In addition, the second embodiment and the third embodiment can be combined and changed appropriately so as to perform global point tracking in consideration of the difference in physical properties of the material constituting the layer.

本発明は、包装容器などの成形を行うにあたり、成形後に所望の層配置になるような初期の層構成の予測に利用することができる。特に、大変形を伴う多層の樹脂の成形解析に好適に利用することが可能である。   INDUSTRIAL APPLICABILITY The present invention can be used for predicting an initial layer configuration so as to obtain a desired layer arrangement after molding when molding a packaging container or the like. In particular, it can be suitably used for molding analysis of multi-layered resins with large deformation.

本発明の第一実施形態の圧縮成形解析装置の構成を示すブロック図である。It is a block diagram which shows the structure of the compression molding analyzer of 1st embodiment of this invention. 本発明の第一実施形態の点追跡法における点追跡の様子を示す図である。It is a figure which shows the mode of the point tracking in the point tracking method of 1st embodiment of this invention. 本発明の第一実施形態の圧縮成形解析装置により実行される層形状変化の解析処理を示すフローチャートである。It is a flowchart which shows the analysis process of the layer shape change performed by the compression molding analyzer of 1st embodiment of this invention. 本発明の第一実施形態の圧縮成形解析装置により実行される点追跡処理を示すフローチャートである。It is a flowchart which shows the point tracking process performed by the compression molding analyzer of 1st embodiment of this invention. 本発明の第二実施形態の圧縮成形解析装置の構成を示すブロック図である。It is a block diagram which shows the structure of the compression molding analyzer of 2nd embodiment of this invention. 本発明の第二実施形態の大域的点追跡法における点追跡の様子を示す図である。It is a figure which shows the mode of the point tracking in the global point tracking method of 2nd embodiment of this invention. 本発明の第二実施形態の大域的点追跡法における初期形状上の中間層分布を算出する様子を示す図である。It is a figure which shows a mode that the intermediate layer distribution on the initial shape is calculated in the global point tracking method of 2nd embodiment of this invention. 本発明の第二実施形態の圧縮成形解析装置における大域的点追跡処理を示すフローチャートである。It is a flowchart which shows the global point tracking process in the compression molding analyzer of 2nd embodiment of this invention. 本発明の第二実施形態の圧縮成形解析装置における初期中間層分布生成処理を示すフローチャートである。It is a flowchart which shows the initial stage intermediate layer distribution production | generation process in the compression molding analyzer of 2nd embodiment of this invention. 本発明の第三実施形態の圧縮成形解析装置における物性の違いを考慮した層形状変化の解析処理を示すフローチャートである。It is a flowchart which shows the analysis process of the layer shape change in consideration of the difference in the physical property in the compression molding analyzer of 3rd embodiment of this invention. 従来のアダプティブメッシュ法による層形状変化の解析の様子(1)を示す図である。It is a figure which shows the mode (1) of the analysis of the layer shape change by the conventional adaptive mesh method. 従来のアダプティブメッシュ法による層形状変化の解析の様子(2)を示す図である。It is a figure which shows the mode (2) of the analysis of the layer shape change by the conventional adaptive mesh method.

符号の説明Explanation of symbols

10 圧縮成形解析装置
10a モデルファイル記憶部
10b 解析条件入力部
10c 層形状変化解析部
10d 解析結果記憶部
10e 中間層形状情報入力部
10f 点追跡部
10g 点追跡結果記憶部
10h 大域的点追跡部
10i 大域的点追跡結果記憶部
10j 成形後中間層形状情報入力部
10k 初期中間層形状情報生成部
10l 初期中間層形状情報記憶部
DESCRIPTION OF SYMBOLS 10 Compression molding analyzer 10a Model file memory | storage part 10b Analysis condition input part 10c Layer shape change analysis part 10d Analysis result memory | storage part 10e Intermediate | middle layer shape information input part 10f Point tracking part 10g Point tracking result memory | storage part 10h Global point tracking part 10i Global point tracking result storage unit 10j Post-molding intermediate layer shape information input unit 10k Initial intermediate layer shape information generation unit 10l Initial intermediate layer shape information storage unit

Claims (8)

成形解析装置に、成形過程における層分布を算出させる成形解析手法であって、
前記成形解析装置における層形状変化解析手段が、層を所定の解析条件にもとづき解析して、前記層に割り振られたメッシュの節点の座標情報と該節点の速度ベクトルを含む物理量とからなる解析結果を所定の時間ステップごとに算出し、
前記成形解析装置における点追跡手段が、前記所定の時間ステップごとに、前記解析結果を読み込んで、この解析結果と前記層内に想定された一又は二以上の任意の点の座標情報とにもとづき各点における速度ベクトルを算出するとともに、前記所定の時間ステップ後の各点の座標情報を算出し、かつ、前記解析結果の読み込み、前記速度ベクトルの算出、及び前記座標情報の算出を所定時間が経過するまで繰り返し、前記層内に想定された任意の点の所定時間後の座標情報を算出する
ことを特徴とする成形解析手法。
A molding analysis method for causing a molding analysis device to calculate a layer distribution in the molding process,
The layer shape change analysis means in the forming analysis device analyzes the layer based on a predetermined analysis condition, and an analysis result comprising the coordinate information of the mesh node allocated to the layer and a physical quantity including the velocity vector of the node For each given time step,
The point tracking means in the molding analysis apparatus reads the analysis result at each predetermined time step, and based on the analysis result and the coordinate information of one or more arbitrary points assumed in the layer. The velocity vector at each point is calculated, the coordinate information of each point after the predetermined time step is calculated, the analysis result is read, the velocity vector is calculated, and the coordinate information is calculated for a predetermined time. It repeats until it passes, and calculates the coordinate information after the predetermined time of the arbitrary points assumed in the said layer. The shaping | molding analysis method characterized by the above-mentioned.
前記成形解析装置における大域的点追跡手段が、前記所定の時間ステップごとに、前記解析結果を読み込んで、前記解析結果における最初の時間ステップの全ての節点の座標情報及び初期座標値を大域的点追跡結果記憶手段に記憶させ、当該節点を前記点追跡手段により前記所定の時間ステップの間点追跡して、前記所定の時間ステップ後の各節点の座標情報を算出し、当該所定の時間ステップ後の各節点の初期座標値分布を生成して、かつ、前記解析結果の読み込み、当該解析結果にもとづく前記所定の時間ステップ後の各節点の座標情報及び初期座標値の記憶、前記所定の時間ステップ後の各節点に対する点追跡によるさらに前記所定の時間ステップ後の各節点の座標情報の算出、及び当該所定の時間ステップ後の各節点の初期座標値分布の生成を所定時間が経過するまで繰り返し、前記解析結果における最初の時間ステップの全ての節点の所定時間後の座標情報と初期座標値分布を算出する
ことを特徴とする請求項1記載の成形解析手法。
The global point tracking means in the molding analysis apparatus reads the analysis result at each predetermined time step, and sets the coordinate information and initial coordinate values of all the nodes in the first time step in the analysis result as global points. The tracking result is stored in the tracking result storage unit, the node is tracked during the predetermined time step by the point tracking unit, the coordinate information of each node after the predetermined time step is calculated, and after the predetermined time step Generating an initial coordinate value distribution of each node and reading the analysis result, storing the coordinate information and initial coordinate value of each node after the predetermined time step based on the analysis result, the predetermined time step Further calculation of the coordinate information of each node after the predetermined time step by point tracking for each subsequent node, and the initial coordinate value of each node after the predetermined time step 2. The molding analysis according to claim 1, wherein generation of coordinates is repeated until a predetermined time elapses, and coordinate information and initial coordinate value distribution after a predetermined time of all nodes in the first time step in the analysis result are calculated. Technique.
前記成形解析装置における初期中間層形状生成手段が、前記大域的点追跡手段により算出された前記全ての節点の所定時間後の座標情報により構成される前記層全体の最終形状上に特定された任意の中間層形状における多数の点の座標情報にもとづいて、当該点に対応する初期座標値を前記大域的点追跡結果記憶手段に記憶された座標情報及び初期座標値にもとづき算出し、これらの初期座標値の点を前記層全体の初期形状情報に重ね合わせて、前記中間層形状に対応する初期形状情報を生成する
ことを特徴とする請求項2記載の成形解析手法。
Arbitrary intermediate layer shape generating means in the molding analysis device is specified on the final shape of the entire layer composed of coordinate information after a predetermined time of all the nodes calculated by the global point tracking means Based on the coordinate information of a large number of points in the intermediate layer shape, initial coordinate values corresponding to the points are calculated based on the coordinate information and initial coordinate values stored in the global point tracking result storage means, and these initial values are calculated. The molding analysis method according to claim 2, wherein the initial shape information corresponding to the intermediate layer shape is generated by superimposing the coordinate value points on the initial shape information of the entire layer.
前記点追跡手段が、
前記層内に想定された一又は二以上の任意の点の座標情報が前記解析結果における節点の座標情報に一致しない場合に、前記解析結果における、前記任意の点の周囲の節点に対応する速度ベクトルを用いて、補間計算により前記任意の点の速度ベクトルを算出する
ことを特徴とする請求項1〜3のいずれかに記載の成形解析手法。
The point tracking means comprises:
When the coordinate information of one or more arbitrary points assumed in the layer does not match the coordinate information of the nodes in the analysis result, the speed corresponding to the nodes around the arbitrary point in the analysis result The molding analysis method according to claim 1, wherein a velocity vector of the arbitrary point is calculated by interpolation calculation using a vector.
成形過程における層分布を算出する成形解析装置であって、
層を所定の解析条件にもとづき解析し、前記層に割り振られたメッシュの節点の座標情報と該節点の速度ベクトルを含む物理量とからなる解析結果を所定の時間ステップごとに算出する層形状変化解析手段と、
前記所定の時間ステップごとに、前記解析結果を読み込んで、この解析結果と前記層内に想定された一又は二以上の任意の点の座標情報とにもとづき各点における速度ベクトルを算出するとともに、前記所定の時間ステップ後の各点の座標情報を算出し、かつ、前記解析結果の読み込み、前記速度ベクトルの算出、及び前記座標情報の算出を所定時間が経過するまで繰り返し、前記層内に想定された任意の点の所定時間後の座標情報を算出する点追跡手段と、を備えた
ことを特徴とする成形解析装置。
A molding analysis device for calculating a layer distribution in a molding process,
Layer shape change analysis that analyzes a layer based on a predetermined analysis condition and calculates an analysis result consisting of the coordinate information of the nodes of the mesh assigned to the layer and a physical quantity including the velocity vector of the node for each predetermined time step Means,
For each predetermined time step, read the analysis result, calculate the velocity vector at each point based on the analysis result and the coordinate information of one or more arbitrary points assumed in the layer, Calculate the coordinate information of each point after the predetermined time step, and repeat the reading of the analysis result, the calculation of the velocity vector, and the calculation of the coordinate information until a predetermined time elapses, and assume in the layer And a point tracking means for calculating coordinate information of the given arbitrary point after a predetermined time.
前記所定の時間ステップごとに、前記解析結果を読み込んで、前記解析結果における最初の時間ステップの全ての節点の座標情報及び初期座標値を大域的点追跡結果記憶手段に記憶させ、当該節点を前記点追跡手段により前記所定の時間ステップの間点追跡して、前記所定の時間ステップ後の各節点の座標情報を算出し、当該所定の時間ステップ後の各節点の初期座標値分布を生成して、かつ、前記解析結果の読み込み、当該解析結果にもとづく前記所定の時間ステップ後の各節点の座標情報及び初期座標値の記憶、前記所定の時間ステップ後の各節点に対する点追跡によるさらに前記所定の時間ステップ後の各節点の座標情報の算出、及び当該所定の時間ステップ後の各節点の初期座標値分布の生成を所定時間が経過するまで繰り返し、前記解析結果における最初の時間ステップの全ての節点の所定時間後の座標情報と初期座標値分布を算出する大域的点追跡手段を備えた
ことを特徴とする請求項5記載の成形解析装置。
For each predetermined time step, the analysis result is read, and the coordinate information and initial coordinate values of all nodes in the first time step in the analysis result are stored in the global point tracking result storage means, and the node is The point tracking means performs point tracking during the predetermined time step, calculates coordinate information of each node after the predetermined time step, and generates initial coordinate value distribution of each node after the predetermined time step. And reading the analysis result, storing the coordinate information and initial coordinate value of each node after the predetermined time step based on the analysis result, and further tracking the point for each node after the predetermined time step. Repeat the calculation of the coordinate information of each node after the time step and the generation of the initial coordinate value distribution of each node after the predetermined time step until the predetermined time has passed. The first time all nodes of a predetermined time after the forming analysis apparatus according to claim 5, further comprising a global point tracking means for calculating the coordinate information and the initial coordinate value distribution step in the analysis results.
前記大域的点追跡手段により算出された前記全ての節点の所定時間後の座標情報により構成される前記層全体の最終形状上に特定された任意の中間層形状における多数の点の座標情報にもとづいて、当該点に対応する初期座標値を前記大域的点追跡結果記憶手段に記憶された座標情報及び初期座標値にもとづき算出し、これらの初期座標値の点を前記層全体の初期形状情報に重ね合わせて、前記中間層形状に対応する初期形状情報を生成する初期中間層形状生成手段を備えた
ことを特徴とする請求項6記載の成形解析装置。
Based on the coordinate information of a large number of points in any intermediate layer shape specified on the final shape of the entire layer composed of coordinate information after a predetermined time of all the nodes calculated by the global point tracking means The initial coordinate value corresponding to the point is calculated based on the coordinate information and the initial coordinate value stored in the global point tracking result storage means, and the point of these initial coordinate values is used as the initial shape information of the entire layer. The molding analysis apparatus according to claim 6, further comprising initial intermediate layer shape generation means for generating initial shape information corresponding to the intermediate layer shape in a superimposed manner.
前記点追跡手段が、
前記層内に想定された一又は二以上の任意の点の座標情報が前記解析結果における節点の座標情報に一致しない場合に、前記解析結果における、前記任意の点の周囲の節点に対応する速度ベクトルを用いて、補間計算により前記任意の点の速度ベクトルを算出する
ことを特徴とする請求項5〜7のいずれかに記載の成形解析装置。
The point tracking means comprises:
When the coordinate information of one or more arbitrary points assumed in the layer does not match the coordinate information of the nodes in the analysis result, the speed corresponding to the nodes around the arbitrary point in the analysis result The molding analysis apparatus according to claim 5, wherein a velocity vector of the arbitrary point is calculated by interpolation calculation using a vector.
JP2005335645A 2005-11-21 2005-11-21 Molding analysis method and molding analysis apparatus Active JP4715462B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005335645A JP4715462B2 (en) 2005-11-21 2005-11-21 Molding analysis method and molding analysis apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005335645A JP4715462B2 (en) 2005-11-21 2005-11-21 Molding analysis method and molding analysis apparatus

Publications (2)

Publication Number Publication Date
JP2007141046A JP2007141046A (en) 2007-06-07
JP4715462B2 true JP4715462B2 (en) 2011-07-06

Family

ID=38203810

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005335645A Active JP4715462B2 (en) 2005-11-21 2005-11-21 Molding analysis method and molding analysis apparatus

Country Status (1)

Country Link
JP (1) JP4715462B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6420881B1 (en) * 2017-10-31 2018-11-07 株式会社Jsol Compression molding analysis system, compression molding analysis method, and compression molding analysis program
CN109859250B (en) * 2018-11-20 2023-08-18 北京悦图遥感科技发展有限公司 Aviation infrared video multi-target detection and tracking method and device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6234282A (en) * 1985-08-08 1987-02-14 Hitachi Ltd Molding process simulation system
JPH071560A (en) * 1993-06-18 1995-01-06 Sekisui Chem Co Ltd Layer distribution control condition-setting device in molding process of multilayered resin film
JPH07137108A (en) * 1993-11-16 1995-05-30 Sekisui Chem Co Ltd Method and device for setting compressing time in injection compression molding
JP2001277320A (en) * 2000-03-29 2001-10-09 Polyplastics Co Molding testing system and method for testing molding
JP2004288121A (en) * 2003-03-25 2004-10-14 Toray Ind Inc Analytical model generating method, generator, program and storage medium
JP2004318863A (en) * 2003-04-02 2004-11-11 Toray Ind Inc Design support method of molded article and equipment therefor

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6234282A (en) * 1985-08-08 1987-02-14 Hitachi Ltd Molding process simulation system
JPH071560A (en) * 1993-06-18 1995-01-06 Sekisui Chem Co Ltd Layer distribution control condition-setting device in molding process of multilayered resin film
JPH07137108A (en) * 1993-11-16 1995-05-30 Sekisui Chem Co Ltd Method and device for setting compressing time in injection compression molding
JP2001277320A (en) * 2000-03-29 2001-10-09 Polyplastics Co Molding testing system and method for testing molding
JP2004288121A (en) * 2003-03-25 2004-10-14 Toray Ind Inc Analytical model generating method, generator, program and storage medium
JP2004318863A (en) * 2003-04-02 2004-11-11 Toray Ind Inc Design support method of molded article and equipment therefor

Also Published As

Publication number Publication date
JP2007141046A (en) 2007-06-07

Similar Documents

Publication Publication Date Title
US11256832B2 (en) Replica selection
JP7164295B2 (en) Orienting real objects for 3D printing
JP7343275B2 (en) Part design using topology optimization
JP6850116B2 (en) Notch detection
US11604902B2 (en) Modeling using a weak type definition
JP7456767B2 (en) Mechanical part design
JP2023071722A (en) Data set for learning function using image as input
CN105761303A (en) Creation Of Bounding Boxes On 3d Modeled Assembly
JP2018106713A (en) B-rep of result of two-axis 3d printing process
CN107428086A (en) Object data for increasing material manufacturing represents
JP6557029B2 (en) Designing industrial products using geometry connected by geometric constraints
US20100156936A1 (en) Deformation method of analysis model and computer
JP4715462B2 (en) Molding analysis method and molding analysis apparatus
JP7477279B2 (en) 3D design of B-Rep skin
WO2006027882A1 (en) Automatic construction system for three-dimensional model
JP6867797B2 (en) Exterior design of composite parts
US20180169955A1 (en) Method and equipment for generating a numerical representation of a three-dimensional object, said numerical representation being suited to be used for making said three-dimensional object through stereolithography
JP6949415B2 (en) Mold CAD model data creation device and mold CAD model data creation method
US8712737B1 (en) Use of physical deformation during scanning of an object to generate views of the object
EP3674931A1 (en) Flexible modeling using a weak type definition
JP6972647B2 (en) 3D shape data editing device and 3D shape data editing program
KR20140147761A (en) Designing a folded sheet object
JP2021131852A (en) Simulation for structure of machine component
JP5889077B2 (en) Molded product shrinkage deformation prediction apparatus, molded product shrinkage deformation prediction method, and molded product shrinkage deformation prediction program
JP5899004B2 (en) Molded product warpage deformation prediction device, molded product warpage deformation prediction method, and molded product warpage deformation prediction program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20081023

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20110214

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110314

R150 Certificate of patent or registration of utility model

Ref document number: 4715462

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20140408

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20140408

Year of fee payment: 3

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

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

Free format text: PAYMENT UNTIL: 20140408

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

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

Free format text: PAYMENT UNTIL: 20140408

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350