[go: up one dir, main page]

JP7645207B2 - Method, program, recording medium, and apparatus for analyzing resin molded products - Google Patents

Method, program, recording medium, and apparatus for analyzing resin molded products Download PDF

Info

Publication number
JP7645207B2
JP7645207B2 JP2022006020A JP2022006020A JP7645207B2 JP 7645207 B2 JP7645207 B2 JP 7645207B2 JP 2022006020 A JP2022006020 A JP 2022006020A JP 2022006020 A JP2022006020 A JP 2022006020A JP 7645207 B2 JP7645207 B2 JP 7645207B2
Authority
JP
Japan
Prior art keywords
orientation
molded product
nodes
resin molded
shrinkage
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
JP2022006020A
Other languages
Japanese (ja)
Other versions
JP2023104808A (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.)
Mazda Motor Corp
DaikyoNishikawa Corp
Toray Engineering Co Ltd
Original Assignee
Mazda Motor Corp
DaikyoNishikawa Corp
Toray Engineering Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mazda Motor Corp, DaikyoNishikawa Corp, Toray Engineering Co Ltd filed Critical Mazda Motor Corp
Priority to JP2022006020A priority Critical patent/JP7645207B2/en
Publication of JP2023104808A publication Critical patent/JP2023104808A/en
Application granted granted Critical
Publication of JP7645207B2 publication Critical patent/JP7645207B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Injection Moulding Of Plastics Or The Like (AREA)
  • Moulds For Moulding Plastics Or The Like (AREA)

Description

この発明は、繊維を含む樹脂材料によって成形される樹脂成形品の収縮量を解析する樹脂成形品解析方法、プログラム、記録媒体、および、樹脂成形品解析装置に関する。 This invention relates to a resin molded product analysis method, program, recording medium, and resin molded product analysis device for analyzing the amount of shrinkage of a resin molded product made from a resin material containing fibers.

従来、繊維強化樹脂を用いた射出成形品の形状予測を行う形状予測装置が知られている(たとえば、特許文献1参照)。 Conventionally, a shape prediction device that predicts the shape of an injection molded product using fiber-reinforced resin is known (see, for example, Patent Document 1).

上記特許文献1には、樹脂に混入された繊維の繊維配向を解析することによって、繊維強化樹脂を用いた射出成形品の形状を予測する形状予測装置が開示されている。この形状予測装置では、入力された射出成形品の形状に対応する樹脂流動解析モデルが、有限個の樹脂流動解析要素に分割される。そして、この形状予測装置は、分割された樹脂流動解析要素ごとに樹脂流動解析を実行することによって、樹脂内に混入された繊維の繊維配向データを取得する。そして、この形状予測装置は、樹脂流動解析要素よりも大きい構造解析用要素ごとの繊維配向データを取得するために、複数の樹脂流動解析要素の繊維配向データの合成ベクトルを求めることによって繊維配向データを変換する。上記特許文献1に記載されている形状予測装置は、合成された繊維配向データを用いて構造解析を実行することによって、射出成形品の形状予測を行う。 The above-mentioned Patent Document 1 discloses a shape prediction device that predicts the shape of an injection molded product using fiber-reinforced resin by analyzing the fiber orientation of the fibers mixed in the resin. In this shape prediction device, a resin flow analysis model corresponding to the shape of an input injection molded product is divided into a finite number of resin flow analysis elements. Then, this shape prediction device acquires fiber orientation data of the fibers mixed in the resin by performing resin flow analysis for each divided resin flow analysis element. Then, in order to acquire fiber orientation data for each structural analysis element that is larger than the resin flow analysis element, this shape prediction device converts the fiber orientation data by determining a composite vector of the fiber orientation data of multiple resin flow analysis elements. The shape prediction device described in the above-mentioned Patent Document 1 predicts the shape of an injection molded product by performing structural analysis using the composite fiber orientation data.

特開2008-120089号公報JP 2008-120089 A

ここで、上記特許文献1には明記されていないが、上記特許文献1に記載の形状予測装置のように、比較的細かい要素に分割した状態で解析された繊維配向を合成することによって射出成形品(樹脂成形品)の形状予測を行う場合に、樹脂成形品の厚み方向に沿って比較的細かい要素(仮想節点)に分割して繊維配向の分布を取得するとともに、樹脂成形品の厚み方向における繊維配向の分布を合成することによって処理負担を軽減しながら収縮量を解析することが考えられる。しかしながら、樹脂成形品は、金型内において冷却されることによって厚み方向の表層側から順に固化していくため、厚み方向における中心部分が冷却されて固化する場合には、すでに先に固化した表層側の樹脂からの影響を受けながら固化する。そのため、樹脂成形品の厚み方向における繊維配向の分布を合成することによって、合成された繊維配向に基づいて樹脂成形品の収縮量を解析する場合には、厚み方向の表層側から順に固化することによる影響が考慮されていないため、繊維配向に基づいて解析された収縮量が実際の収縮量とは異なるという問題点がある。 Here, although it is not specified in the above Patent Document 1, when predicting the shape of an injection molded product (resin molded product) by synthesizing fiber orientations analyzed in a state where the product is divided into relatively fine elements, as in the shape prediction device described in the above Patent Document 1, it is possible to obtain the distribution of fiber orientation by dividing the resin molded product into relatively fine elements (virtual nodes) along the thickness direction, and to analyze the shrinkage amount while reducing the processing load by synthesizing the distribution of fiber orientation in the thickness direction of the resin molded product. However, since the resin molded product solidifies from the surface layer side in the thickness direction as it is cooled in the mold, when the center part in the thickness direction is cooled and solidified, it solidifies while being influenced by the resin on the surface side that has already solidified. Therefore, when analyzing the shrinkage amount of a resin molded product based on the synthesized fiber orientation by synthesizing the distribution of fiber orientation in the thickness direction of the resin molded product, there is a problem that the shrinkage amount analyzed based on the fiber orientation differs from the actual shrinkage amount because the influence of solidification from the surface layer side in the thickness direction is not taken into account.

この発明は、上記のような課題を解決するためになされたものであり、この発明の1つの目的は、樹脂成形品の厚み方向における繊維配向の分布を合成することによって、合成された繊維配向に基づいて樹脂成形品の収縮量を解析する場合に、樹脂成形品の収縮量を正確に解析することが可能な樹脂成形品解析方法、プログラム、記録媒体、および、樹脂成形品解析装置を提供することである。 This invention has been made to solve the above problems, and one object of the invention is to provide a method, program, recording medium, and device for analyzing a resin molded product that can accurately analyze the amount of shrinkage of a resin molded product when analyzing the amount of shrinkage of the resin molded product based on the synthesized fiber orientation by synthesizing the distribution of fiber orientation in the thickness direction of the resin molded product.

上記目的を達成するために、この発明の第1の局面による樹脂成形品解析方法は、繊維を含む樹脂材料によって成形される樹脂成形品の収縮量を解析する樹脂成形品解析方法であって、樹脂成形品を複数の微小要素に分割するステップと、複数の微小要素の各々、もしくは、各微小要素の頂点を構成する複数の節点の各々において、樹脂成形品の厚み方向に沿って複数の仮想節点を設定するステップと、各微小要素もしくは各節点の複数の仮想節点の各々において、繊維の配向を解析することによって解析配向を取得するステップと、各微小要素もしくは各節点の複数の仮想節点の各々において、樹脂成形品の表層側から順に解析配向を補正する計算を行うことによって、収縮量を解析するための補正配向を算出するステップと、各微小要素もしくは各節点の複数の仮想節点の各々において算出された補正配向に基づいて、複数の微小要素もしくは複数の節点の各々の収縮量を解析するとともに、解析された複数の微小要素もしくは複数の節点の各々における収縮量に基づいて、樹脂成形品の収縮量を解析するステップと、を備え、補正配向を算出するステップは、各微小要素もしくは各節点の複数の仮想節点の各々において、自身よりも表層側の仮想節点における補正結果に基づいて解析配向を補正する計算を行うことによって、補正配向を算出するステップを含む。 In order to achieve the above object, a resin molded product analysis method according to a first aspect of the present invention is a resin molded product analysis method for analyzing the amount of shrinkage of a resin molded product molded from a resin material containing fibers, comprising the steps of: dividing the resin molded product into a plurality of microelements; setting a plurality of virtual nodes along the thickness direction of the resin molded product for each of the plurality of microelements or for each of the plurality of nodes constituting the vertices of each microelement; obtaining an analytical orientation by analyzing the orientation of fibers for each of the plurality of virtual nodes of each microelement or for each of the nodes; and obtaining an analytical orientation for each of the plurality of virtual nodes of each microelement or for each of the nodes in order from the surface side of the resin molded product. The method includes the steps of: calculating a correction orientation for analyzing the amount of shrinkage by performing a calculation to correct the amount of shrinkage; analyzing the amount of shrinkage of each of the multiple microelements or multiple nodes based on the correction orientation calculated for each of the multiple virtual nodes of each microelement or each node, and analyzing the amount of shrinkage of the resin molded product based on the amount of shrinkage at each of the analyzed multiple microelements or multiple nodes; the step of calculating the correction orientation includes the step of calculating the correction orientation by performing a calculation to correct the analysis orientation at each of the multiple virtual nodes of each microelement or each node based on the correction result at the virtual node on the surface side of the microelement or each node.

この第1の局面による樹脂成形品解析方法は、上記のように、各微小要素もしくは各節点の複数の仮想節点の各々において、樹脂成形品の表層側から順に解析配向を補正する計算を行うことによって、収縮量を解析するための補正配向を算出する。そして、各微小要素もしくは各節点の複数の仮想節点の各々において、自身よりも表層側の仮想節点における補正結果に基づいて解析配向を補正する計算を行うことによって、補正配向を算出する。これにより、各微小要素もしくは各節点の複数の仮想節点の各々において解析された繊維配向である解析配向に対して、樹脂成形品の表層側から順に、自身よりも表層側の仮想節点における補正結果に基づいて補正する計算を行うため、表層側から順に樹脂材料が固化することによる影響を考慮して、解析された繊維配向(解析配向)を補正することができる。そのため、表層側から順に樹脂材料が固化することによる影響を考慮して補正された繊維配向の分布を合成することによって、樹脂成形品の収縮量を解析することができる。その結果、樹脂成形品の厚み方向における繊維配向の分布を合成することによって、合成された繊維配向に基づいて樹脂成形品の収縮量を解析する場合に、樹脂成形品の収縮量を正確に解析することができる。 In the resin molded product analysis method according to the first aspect, as described above, a correction orientation for analyzing the amount of shrinkage is calculated by performing a calculation to correct the analytical orientation in each of the multiple virtual nodes of each microelement or each node, starting from the surface layer side of the resin molded product. Then, a correction orientation is calculated by performing a calculation to correct the analytical orientation based on the correction result in the virtual node on the surface side of each microelement or each node, starting from the surface side of the resin molded product. As a result, a correction calculation is performed based on the correction result in the virtual node on the surface side of each microelement or each node, starting from the surface side of the resin molded product, for the analytical orientation, which is the fiber orientation analyzed in each of the multiple virtual nodes of each microelement or each node, starting from the surface side of the resin molded product, to correct the analyzed fiber orientation (analysis orientation) taking into account the effect of the resin material solidifying in sequence from the surface side. Therefore, the amount of shrinkage of the resin molded product can be analyzed by synthesizing the distribution of the fiber orientation corrected taking into account the effect of the resin material solidifying in sequence from the surface side. As a result, by synthesizing the distribution of fiber orientation in the thickness direction of a resin molded product, the amount of shrinkage of the resin molded product can be accurately analyzed when analyzing the amount of shrinkage of the resin molded product based on the synthesized fiber orientation.

上記第1の局面による樹脂成形品解析方法において、好ましくは、解析配向を取得するステップは、複数の仮想節点のうちの一の仮想節点である対象仮想節点において、対象仮想節点における解析配向である対象解析配向を取得するステップを含み、補正配向を算出するステップは、複数の仮想節点のうちの対象仮想節点の表層側に隣接する隣接仮想節点において算出された補正配向である隣接補正配向に基づいて、対象解析配向を補正する計算を行うことによって、対象仮想節点における補正配向である対象補正配向を算出するステップと、対象補正配向を算出するステップを、複数の仮想節点の各々において、樹脂成形品の表層側から順に繰り返すステップとを含む。このように構成すれば、生成される樹脂成形品の厚み方向において設定された複数の仮想節点の各々において、表層側に隣り合う仮想節点における補正結果を用いて、1層ずつ順に表層側から内側に向かって補正配向を算出することができる。そのため、複数の仮想節点の各々において、自身よりも表層側の全ての仮想節点における補正結果を用いる場合に比べて、補正配向の算出に用いる仮想節点の数を少なくすることができるので、収縮量を解析するための解析配向に対する補正を容易に行うことができる。その結果、表層側に隣り合う仮想節点における補正結果を用いることによって、合成された補正後の繊維配向に基づいて樹脂成形品の収縮量を容易にかつ正確に解析することができる。 In the resin molded product analysis method according to the first aspect, preferably, the step of acquiring the analysis orientation includes a step of acquiring a target analysis orientation, which is an analysis orientation at a target virtual node, which is one of the multiple virtual nodes, and the step of calculating the correction orientation includes a step of calculating a target correction orientation, which is a correction orientation at the target virtual node, by performing a calculation to correct the target analysis orientation based on an adjacent correction orientation, which is a correction orientation calculated at an adjacent virtual node adjacent to the surface side of the target virtual node among the multiple virtual nodes, and a step of repeating the step of calculating the target correction orientation in order from the surface side of the resin molded product for each of the multiple virtual nodes. With this configuration, it is possible to calculate the correction orientation from the surface side to the inside, layer by layer, using the correction result at the virtual node adjacent to the surface side at each of the multiple virtual nodes set in the thickness direction of the generated resin molded product. Therefore, for each of the multiple virtual nodes, the number of virtual nodes used to calculate the corrected orientation can be reduced compared to when the correction results for all virtual nodes on the surface side of the virtual node itself are used, making it easier to correct the analytical orientation for analyzing the amount of shrinkage. As a result, by using the correction results for the virtual nodes adjacent to the surface side, the amount of shrinkage of the resin molded product can be easily and accurately analyzed based on the composite corrected fiber orientation.

この場合、好ましくは、対象補正配向を算出するステップは、隣接補正配向と隣接仮想節点における樹脂材料の弾性率とに基づいて、対象解析配向を補正する計算を行うことによって、対象補正配向を算出するステップを含む。ここで、表層側に隣接する層が先に冷却されて固化した場合には、すでに固化した隣接する層との弾性率の差異に起因して、実際の収縮量と解析された繊維配向に基づいて算出された収縮量とに差が生じると考えられる。これを考慮して、本発明では、隣接補正配向と隣接仮想節点における樹脂材料の弾性率とに基づいて、対象解析配向を補正する計算を行うことによって、対象補正配向を算出する。このように構成すれば、表層側に隣り合う仮想節点における補正配向に加えて、表層側に隣り合う仮想節点における弾性率をも考慮にいれて、解析配向を補正することができる。そのため、収縮量を解析するために解析された繊維配向である解析配向をより精度よく補正することができる。その結果、補正された繊維配向である補正配向に基づいて収縮量を解析する場合に、収縮量をより精度よく正確に解析することができる。 In this case, the step of calculating the target correction orientation preferably includes a step of calculating the target correction orientation by performing a calculation to correct the target analysis orientation based on the adjacent correction orientation and the elastic modulus of the resin material at the adjacent virtual node. Here, if the layer adjacent to the surface layer side is cooled and solidified first, it is considered that a difference occurs between the actual shrinkage amount and the shrinkage amount calculated based on the analyzed fiber orientation due to the difference in elastic modulus with the adjacent layer that has already solidified. In consideration of this, in the present invention, the target correction orientation is calculated by performing a calculation to correct the target analysis orientation based on the adjacent correction orientation and the elastic modulus of the resin material at the adjacent virtual node. With this configuration, in addition to the correction orientation at the virtual node adjacent to the surface layer side, the elastic modulus at the virtual node adjacent to the surface layer side can also be taken into consideration to correct the analysis orientation. Therefore, the analysis orientation, which is the fiber orientation analyzed to analyze the shrinkage amount, can be corrected more accurately. As a result, when the shrinkage amount is analyzed based on the correction orientation, which is the corrected fiber orientation, the shrinkage amount can be analyzed more accurately and precisely.

上記隣接補正配向と隣接仮想節点における樹脂材料の弾性率とに基づいて、対象解析配向を補正する計算を行うことによって、対象補正配向を算出する樹脂成形品解析方法において、好ましくは、対象補正配向を算出するステップは、隣接仮想節点における樹脂材料の弾性率と、対象仮想節点と隣接仮想節点との間の温度差とに基づいて、補正係数を算出するステップと、隣接補正配向と補正係数とに基づいて、対象解析配向を補正する計算を行うステップとを含む。このように構成すれば、樹脂材料の弾性率は温度に依存して変化するため、表層側に隣接する仮想節点との温度差と隣接する仮想節点における弾性率とに基づいて算出された補正係数を用いることによって、すでに固化した隣接する層との弾性率の差異に起因する影響を考慮しながら、収縮量を解析するために解析配向をより一層精度よく補正することができる。その結果、補正された繊維配向である補正配向を合成することによって収縮量を解析する場合に、収縮量をより一層精度よく正確に解析することができる。 In the resin molding analysis method for calculating the target correction orientation by performing a calculation to correct the target analysis orientation based on the adjacent correction orientation and the elastic modulus of the resin material at the adjacent virtual node, preferably, the step of calculating the target correction orientation includes a step of calculating a correction coefficient based on the elastic modulus of the resin material at the adjacent virtual node and the temperature difference between the target virtual node and the adjacent virtual node, and a step of performing a calculation to correct the target analysis orientation based on the adjacent correction orientation and the correction coefficient. With this configuration, since the elastic modulus of the resin material changes depending on the temperature, the analytical orientation can be corrected more accurately to analyze the amount of shrinkage while considering the influence caused by the difference in the elastic modulus with the adjacent layer that has already solidified by using a correction coefficient calculated based on the temperature difference with the virtual node adjacent to the surface layer side and the elastic modulus at the adjacent virtual node. As a result, when the amount of shrinkage is analyzed by synthesizing the correction orientation, which is the corrected fiber orientation, the amount of shrinkage can be analyzed more accurately and precisely.

この場合、好ましくは、対象解析配向を補正する計算を行うステップは、対象解析配向と、補正係数によって重み付けされた隣接補正配向との平均を計算することによって、対象解析配向を補正する計算を行うステップを含む。このように構成すれば、補正係数によって重み付けされた隣接補正配向との平均を計算することによって、解析配向を補正する計算をより容易に行うことができる。その結果、補正された解析配向である補正配向に基づいて、収縮量をより容易に解析することができる。 In this case, preferably, the step of performing a calculation to correct the target analysis orientation includes a step of performing a calculation to correct the target analysis orientation by calculating an average of the target analysis orientation and an adjacent correction orientation weighted by a correction coefficient. With this configuration, the calculation to correct the analysis orientation can be more easily performed by calculating an average of the target analysis orientation and an adjacent correction orientation weighted by a correction coefficient. As a result, the amount of shrinkage can be more easily analyzed based on the correction orientation, which is the corrected analysis orientation.

上記第1の局面による樹脂成形品解析方法において、好ましくは、収縮量を解析するステップは、各微小要素もしくは各節点の複数の仮想節点の各々において算出された補正配向の平均に基づいて、複数の微小要素もしくは複数の節点の各々の収縮量を解析するとともに、解析された複数の微小要素もしくは複数の節点の各々における収縮量に基づいて、樹脂成形品の収縮量を解析するステップを含む。このように構成すれば、仮想節点の各々において算出された補正配向の平均を取得することによって、厚み方向において合成された繊維配向を容易に取得することができる。そのため、補正配向に基づいて、各微小要素もしくは各節点ごとに合成された補正配向を容易に取得することができるので、各微小要素もしくは各節点ごとの収縮量を容易に解析することができる。その結果、表層側から順に樹脂材料が固化することによる影響を考慮して、樹脂成形品の収縮量を容易に正確に解析することができる。 In the resin molded product analysis method according to the first aspect, the step of analyzing the shrinkage preferably includes a step of analyzing the shrinkage of each of the multiple microelements or multiple nodes based on the average of the corrected orientation calculated at each of the multiple virtual nodes of each microelement or each node, and analyzing the shrinkage of the resin molded product based on the shrinkage of each of the analyzed multiple microelements or multiple nodes. With this configuration, the fiber orientation synthesized in the thickness direction can be easily obtained by obtaining the average of the corrected orientation calculated at each virtual node. Therefore, the corrected orientation synthesized for each microelement or each node can be easily obtained based on the corrected orientation, so that the shrinkage of each microelement or each node can be easily analyzed. As a result, the shrinkage of the resin molded product can be easily and accurately analyzed, taking into account the effect of the resin material solidifying in order from the surface side.

上記第1の局面による樹脂成形品解析方法において、好ましくは、複数の仮想節点を設定するステップは、複数の微小要素もしくは複数の節点の各々において、樹脂成形品の厚み方向に沿って偶数個の複数の仮想節点を設定するステップを含む。このように構成すれば、偶数個の仮想節点が設定されるため、厚み方向における一方側からの計算と他方側からの計算とを同じ回数ずつ行うことによって、複数の仮想節点の全体において補正配向の算出を表層側から順に行うことができる。その結果、奇数個の仮想節点が設定されている場合とは異なり、余りがでることなく全ての仮想節点において一方側と他方側との両方から順に同様の補正の計算を行うことができるので、全ての仮想節点において補正の計算をより容易に行うことができる。 In the resin molded product analysis method according to the first aspect, the step of setting a plurality of virtual nodes preferably includes a step of setting an even number of virtual nodes along the thickness direction of the resin molded product for each of a plurality of infinitesimal elements or a plurality of nodes. With this configuration, an even number of virtual nodes are set, and by performing calculations from one side and the other side in the thickness direction the same number of times, it is possible to calculate the correction orientation for all of the plurality of virtual nodes in order from the surface layer side. As a result, unlike the case where an odd number of virtual nodes are set, similar correction calculations can be performed in order from both one side and the other side for all virtual nodes without any remainders, making it easier to perform correction calculations for all virtual nodes.

上記第1の局面による樹脂成形品解析方法において、好ましくは、収縮量を解析するステップは、各微小要素もしくは各節点の複数の仮想節点の各々において算出された補正配向に基づいて、複数の微小要素もしくは複数の節点の各々において、樹脂材料の流動方向と、樹脂成形品の表面に沿う面内において流動方向と直交する直交方向との収縮量を解析するとともに、解析された複数の微小要素もしくは複数の節点の各々における流動方向と直交方向との収縮量に基づいて、樹脂成形品の流動方向と直交方向との収縮量を解析するステップを含む。このように構成すれば、樹脂材料の流動方向と直交方向との2つの方向における繊維配向を補正する計算を行うため、流動方向および直交方向に厚み方向をも加えた3つの方向における繊維配向の補正を行う場合に比べて、補正の計算に要する処理負担を軽減することができる。その結果、算出された補正配向に基づいて流動方向と直交方向との2つの方向における収縮量の解析を行うことによって、樹脂成形品における収縮量を解析する場合における処理負担を軽減することができる。 In the resin molded product analysis method according to the first aspect, the step of analyzing the shrinkage amount preferably includes a step of analyzing the shrinkage amount in the flow direction of the resin material and the orthogonal direction orthogonal to the flow direction in a plane along the surface of the resin molded product at each of the multiple microelements or multiple nodes based on the correction orientation calculated at each of the multiple virtual nodes of each microelement or each node, and analyzing the shrinkage amount in the flow direction and the orthogonal direction of the resin molded product based on the shrinkage amount in the flow direction and the orthogonal direction at each of the analyzed multiple microelements or multiple nodes. With this configuration, a calculation is performed to correct the fiber orientation in two directions, the flow direction and the orthogonal direction of the resin material, so the processing load required for the correction calculation can be reduced compared to the case where the fiber orientation is corrected in three directions, including the flow direction, the orthogonal direction, and the thickness direction. As a result, the processing load when analyzing the shrinkage amount in the resin molded product can be reduced by analyzing the shrinkage amount in two directions, the flow direction and the orthogonal direction, based on the calculated correction orientation.

この発明の第2の局面によるプログラムは、第1の局面による樹脂成形品解析方法をコンピュータに実行させる。 A program according to a second aspect of the present invention causes a computer to execute the resin molded product analysis method according to the first aspect.

この発明の第2の局面によるプログラムでは、上記第1の局面による樹脂成形品解析方法をコンピュータに実行させることにより、樹脂成形品の厚み方向における繊維配向の分布を合成することによって、合成された繊維配向に基づいて樹脂成形品の収縮量を解析する場合に、樹脂成形品の収縮量を正確に解析することができる。 In the program according to the second aspect of the present invention, the resin molded product analysis method according to the first aspect is executed on a computer, whereby the distribution of fiber orientation in the thickness direction of the resin molded product is synthesized, and when the shrinkage amount of the resin molded product is analyzed based on the synthesized fiber orientation, the shrinkage amount of the resin molded product can be accurately analyzed.

この発明の第3の局面による記憶媒体は、第2の局面によるプログラムが記録され、コンピュータにより読み取り可能である。 A storage medium according to the third aspect of the present invention has a program according to the second aspect recorded thereon and is computer-readable.

この発明の第3の局面による記憶媒体は、上記第2の局面によるプログラムが記録され、コンピュータにより読み取り可能であることにより、樹脂成形品の厚み方向における繊維配向の分布を合成することによって、合成された繊維配向に基づいて樹脂成形品の収縮量を解析する場合に、樹脂成形品の収縮量を正確に解析することができる。 The storage medium according to the third aspect of the present invention has the program according to the second aspect recorded thereon and is computer-readable, so that when the distribution of fiber orientation in the thickness direction of the resin molded product is synthesized and the amount of shrinkage of the resin molded product is analyzed based on the synthesized fiber orientation, the amount of shrinkage of the resin molded product can be accurately analyzed.

この発明の第4の局面による樹脂成形品解析装置は、繊維を含む樹脂材料によって成形される樹脂成形品の収縮量を解析する樹脂成形品解析装置であって、樹脂成形品を複数の微小要素に分割する微小要素分割部と、複数の微小要素の各々、もしくは、各微小要素の頂点を構成する複数の節点の各々において、樹脂成形品の厚み方向に沿って複数の仮想節点を設定する仮想節点設定部と、各微小要素もしくは各節点の複数の仮想節点の各々において、繊維の配向を解析することによって解析配向を取得する解析配向取得部と、各微小要素もしくは各節点の複数の仮想節点の各々において、樹脂成形品の表層側から順に解析配向を補正する計算を行うことによって、収縮量を解析するための補正配向を算出する補正配向算出部と、各微小要素もしくは各節点の複数の仮想節点の各々において算出された補正配向に基づいて、複数の微小要素もしくは複数の節点の各々の収縮量を解析するとともに、解析された複数の微小要素もしくは複数の節点の各々における収縮量に基づいて、樹脂成形品の収縮量を解析する収縮量解析部と、を備え、補正配向算出部は、各微小要素もしくは各節点の複数の仮想節点の各々において、自身よりも表層側の仮想節点における補正結果に基づいて解析配向を補正する計算を行うことによって、補正配向を算出するように構成されている。 The resin molded product analysis device according to a fourth aspect of the present invention is a resin molded product analysis device that analyzes the amount of shrinkage of a resin molded product molded from a resin material containing fibers, and includes a microelement division unit that divides the resin molded product into a plurality of microelements, a virtual node setting unit that sets a plurality of virtual nodes along the thickness direction of the resin molded product in each of the plurality of microelements or in each of the plurality of nodes constituting the vertices of each microelement, an analytical orientation acquisition unit that acquires an analytical orientation by analyzing the orientation of the fibers in each of the plurality of virtual nodes of each microelement or each node, and a correction unit that corrects the analytical orientation in each of the plurality of virtual nodes of each microelement or each node, in order from the surface side of the resin molded product. The apparatus includes a correction orientation calculation unit that calculates a correction orientation for analyzing the amount of shrinkage by performing a calculation to calculate a correction orientation for analyzing the amount of shrinkage, and a shrinkage amount analysis unit that analyzes the amount of shrinkage of each of the multiple microelements or multiple nodes based on the correction orientation calculated at each of the multiple virtual nodes of each microelement or each node, and analyzes the amount of shrinkage of the resin molded product based on the shrinkage amount at each of the analyzed multiple microelements or multiple nodes. The correction orientation calculation unit is configured to calculate the correction orientation by performing a calculation to correct the analysis orientation at each of the multiple virtual nodes of each microelement or each node based on the correction result at the virtual node on the surface side of the microelement or each node.

この発明の第4の局面による樹脂成形品解析装置は、上記のように、各微小要素もしくは各節点の複数の仮想節点の各々において、樹脂成形品の表層側から順に解析配向を補正する計算を行うことによって、収縮量を解析するための補正配向を算出する。そして、各微小要素もしくは各節点の複数の仮想節点の各々において、自身よりも表層側の仮想節点における補正結果に基づいて解析配向を補正する計算を行うことによって、補正配向を算出する。これにより、各微小要素もしくは各節点の複数の仮想節点の各々において解析された繊維配向である解析配向に対して、樹脂成形品の表層側から順に、自身よりも表層側の仮想節点における補正結果に基づいて補正する計算を行うため、表層側から順に樹脂材料が固化することによる影響を考慮して、解析された繊維配向(解析配向)を補正することができる。そのため、表層側から順に樹脂材料が固化することによる影響を考慮して補正された繊維配向の分布を合成することによって、樹脂成形品の収縮量を解析することができる。その結果、樹脂成形品の厚み方向における繊維配向の分布を合成することによって、合成された繊維配向に基づいて樹脂成形品の収縮量を解析する場合に、樹脂成形品の収縮量を正確に解析することが可能な樹脂成形品解析装置を提供することができる。 The resin molded product analysis device according to the fourth aspect of the present invention, as described above, calculates the corrected orientation for analyzing the shrinkage amount by performing a calculation to correct the analytical orientation in each of the multiple virtual nodes of each microelement or each node, starting from the surface side of the resin molded product. Then, in each of the multiple virtual nodes of each microelement or each node, a calculation is performed to correct the analytical orientation based on the correction result in the virtual node on the surface side of the microelement or each node, and the corrected orientation is calculated. As a result, the analytical orientation, which is the fiber orientation analyzed in each of the multiple virtual nodes of each microelement or each node, is corrected based on the correction result in the virtual node on the surface side of the microelement or each node, starting from the surface side of the resin molded product, so that the analyzed fiber orientation (analysis orientation) can be corrected taking into account the effect of the resin material solidifying in sequence from the surface side. Therefore, the shrinkage amount of the resin molded product can be analyzed by synthesizing the distribution of the fiber orientation corrected taking into account the effect of the resin material solidifying in sequence from the surface side. As a result, it is possible to provide a resin molded product analysis device that can accurately analyze the amount of shrinkage of a resin molded product by synthesizing the distribution of fiber orientation in the thickness direction of the resin molded product and analyzing the amount of shrinkage of the resin molded product based on the synthesized fiber orientation.

本発明によれば、上記のように、樹脂成形品の厚み方向における繊維配向の分布を合成することによって、合成された繊維配向に基づいて樹脂成形品の収縮量を解析する場合に、樹脂成形品の収縮量を正確に解析することができる。 According to the present invention, as described above, by synthesizing the distribution of fiber orientation in the thickness direction of a resin molded product, the amount of shrinkage of the resin molded product can be accurately analyzed when analyzing the amount of shrinkage of the resin molded product based on the synthesized fiber orientation.

本発明の一実施形態における解析装置の構成を示したブロック図である。1 is a block diagram showing a configuration of an analysis device according to an embodiment of the present invention. 樹脂成形品の微小要素の分割を示した模式図である。FIG. 2 is a schematic diagram showing division of minute elements of a resin molded product. 1つの微小要素における仮想節点の設定を説明するための模式図である。FIG. 13 is a schematic diagram for explaining the setting of a virtual node in one infinitesimal element. 樹脂成形品における繊維配向を示した模式図である。FIG. 2 is a schematic diagram showing fiber orientation in a resin molded product. 解析配向と補正配向とを説明するための図である。FIG. 13 is a diagram for explaining an analytical orientation and a corrective orientation. 補正配向の算出を説明するための図である。FIG. 13 is a diagram for explaining calculation of a corrected orientation. 一実施形態による樹脂成形品解析方法を説明するためのフローチャート図である。FIG. 1 is a flowchart illustrating a resin molded product analysis method according to an embodiment. 一実施形態による解析配向の補正方法の詳細を説明するためのフローチャート図である。FIG. 11 is a flowchart for explaining details of a method for correcting analytical orientation according to one embodiment. 本発明の変形例による1つの節点における仮想節点の設定を説明するための模式図である。FIG. 13 is a schematic diagram for explaining setting of a virtual node at one node according to a modified example of the present invention.

以下、本発明を具体化した実施形態を図面に基づいて説明する。 The following describes an embodiment of the present invention with reference to the drawings.

(解析装置の構成)
図1~図6を参照して、本実施形態による解析装置100の構成について説明する。なお、解析装置100は、特許請求の範囲における「樹脂成形品解析装置」および「コンピュータ」の一例である。
(Configuration of the analysis device)
The configuration of an analysis device 100 according to the present embodiment will be described with reference to Figures 1 to 6. The analysis device 100 is an example of the "resin molding analysis device" and "computer" in the claims.

解析装置100は、繊維(強化繊維)を含む樹脂材料によって成形される樹脂成形品1(図2参照)の収縮量を解析することによって、樹脂成形品1の形状(そり変形)の解析を行うことが可能であるように構成されている。樹脂成形品1を成形する樹脂材料は、たとえば、PP樹脂およびABS樹脂などの熱可塑性樹脂に対して、ガラス繊維および炭素繊維などの繊維を複合させた繊維強化樹脂である。そして、樹脂成形品1は、金型内に樹脂材料を充填させる射出成形により成形される。なお、樹脂成形品1は、金属部品などと一体成形されたインサート成形品などを含んでいてもよい。 The analysis device 100 is configured to be able to analyze the shape (warpage) of a resin molded product 1 (see FIG. 2) by analyzing the amount of shrinkage of the resin molded product 1 (see FIG. 2) that is molded from a resin material containing fibers (reinforced fibers). The resin material that molds the resin molded product 1 is, for example, a fiber-reinforced resin that is a composite of thermoplastic resin such as PP resin and ABS resin with fibers such as glass fiber and carbon fiber. The resin molded product 1 is molded by injection molding, in which the resin material is filled into a mold. The resin molded product 1 may also include an insert molded product that is molded integrally with a metal part or the like.

図1に示すように、解析装置100は、CPU(Central Processing Unit)などからなる1または複数のプロセッサ10と、ROM(Read Only Memory)、RAM(Random Access Memory)および記憶装置などを含んだ記憶部20とを備える。記憶部20は、たとえば、ハードディスクドライブや半導体記憶装置などである。解析装置100は、樹脂成形品1の収縮量を解析する作業者に用いられるPC(パーソナルコンピュータ)である。なお、解析装置100にプログラム21を実行させることにより行われる処理の一部または全部が、専用の演算回路等のハードウェアによって行われてもよい。 As shown in FIG. 1, the analysis device 100 includes one or more processors 10, such as a CPU (Central Processing Unit), and a storage unit 20 including a ROM (Read Only Memory), a RAM (Random Access Memory), and a storage device. The storage unit 20 is, for example, a hard disk drive or a semiconductor storage device. The analysis device 100 is a personal computer (PC) used by an operator who analyzes the amount of shrinkage of the resin molded product 1. Note that some or all of the processing performed by the analysis device 100 executing the program 21 may be performed by hardware such as a dedicated arithmetic circuit.

解析装置100は、記憶部20に記憶されたプログラム21をプロセッサ10に実行させることにより、樹脂成形品1の収縮量の解析を行う。プログラム21は、樹脂材料の成形をシミュレートすることによって解析するためのCAE(Computer Aided Engineering)ソフトウェアである。また、プログラム21は、記録媒体101から読み出される他、インターネットなどのネットワークやLAN(Local Area Network)などの伝送経路を介して外部サーバなどから提供されてもよい。記録媒体101は、光学ディスク、磁気ディスク、不揮発性半導体メモリなどのコンピュータ読み取り可能な記録媒体であり、プログラム21が記録されている。 The analysis device 100 analyzes the amount of shrinkage of the resin molded product 1 by having the processor 10 execute the program 21 stored in the storage unit 20. The program 21 is CAE (Computer Aided Engineering) software for performing analysis by simulating the molding of a resin material. The program 21 may be read from the recording medium 101, or may be provided from an external server via a transmission path such as a network such as the Internet or a LAN (Local Area Network). The recording medium 101 is a computer-readable recording medium such as an optical disk, a magnetic disk, or a non-volatile semiconductor memory, and the program 21 is recorded therein.

記憶部20には、プログラム21の他、収縮量の解析を行うために用いられる各種の解析用データ22が記憶される。解析用データ22は、たとえば、樹脂材料の粘度、弾性率、熱伝導率などの熱物性、および、状態線図などの材料物性データを含む。また、解析用データ22は、収縮量の解析を行う際の各種パラメータを含む。解析用データ22の各種パラメータは、たとえば、図示しない操作部に対する作業者の入力操作に基づいて取得される境界条件を含む。境界条件は、ゲート位置、射出速度、樹脂温度、および、金型温度などである。 In addition to the program 21, the memory unit 20 stores various analysis data 22 used to analyze the amount of shrinkage. The analysis data 22 includes, for example, thermophysical properties of the resin material, such as viscosity, elastic modulus, and thermal conductivity, and material property data such as a state diagram. The analysis data 22 also includes various parameters used when analyzing the amount of shrinkage. The various parameters of the analysis data 22 include, for example, boundary conditions acquired based on input operations by an operator to an operation unit (not shown). The boundary conditions include the gate position, injection speed, resin temperature, and mold temperature.

また、解析装置100は、液晶表示装置などの表示部30と、記録媒体101からプログラム21や各種データを読み取るための読取部40とを備えている。読取部40は、記録媒体101の種類に応じたリーダ装置などである。また、記憶部20に記憶される解析用データ22は、作業者が作成した記録媒体101から読取部40により読み出されたり、作業者が外部サーバなどに作成しておいて、伝送経路を介して外部サーバから取得されたりしてもよい。 The analysis device 100 also includes a display unit 30 such as a liquid crystal display device, and a reading unit 40 for reading the program 21 and various data from the recording medium 101. The reading unit 40 is a reader device or the like according to the type of recording medium 101. The analysis data 22 stored in the memory unit 20 may be read by the reading unit 40 from the recording medium 101 created by the worker, or may be created by the worker on an external server and acquired from the external server via a transmission path.

(プロセッサの構成)
本実施形態では、プロセッサ10は、微小要素分割部11、仮想節点設定部12、解析配向取得部13、補正配向算出部14、および、収縮量解析部15を含んでいる。具体的には、ハードウェアとしてのプロセッサ10は、ソフトウェア(プログラム21)の機能ブロックとして、微小要素分割部11、仮想節点設定部12、解析配向取得部13、補正配向算出部14、および、収縮量解析部15を含むように構成されている。
(Processor Configuration)
In this embodiment, the processor 10 includes an infinitesimal element division unit 11, a virtual node setting unit 12, an analytical orientation acquisition unit 13, a correction orientation calculation unit 14, and a shrinkage amount analysis unit 15. Specifically, the processor 10 as hardware is configured to include the infinitesimal element division unit 11, the virtual node setting unit 12, the analytical orientation acquisition unit 13, the correction orientation calculation unit 14, and the shrinkage amount analysis unit 15 as functional blocks of the software (program 21).

図2に示すように、本実施形態では、プロセッサ10の微小要素分割部11は、樹脂成形品1を複数の微小要素に分割する。具体的には、解析装置100は、樹脂成形品1の解析を行う作業者による入力操作などに基づいて、樹脂成形品1の3次元CAD情報などを取得することによって、樹脂成形品1の形状データを取得する。そして、微小要素分割部11は、樹脂成形品1の解析モデルを生成するために、樹脂成形品1の3次元の形状データを複数の微小要素(メッシュ)に分割する。また、微小要素分割部11は、樹脂成形品1の各部分の厚み方向において、微小要素の個数が1つとなるように樹脂成形品1を分割する。また、1つの微小要素は、直方体形状を有する。 As shown in FIG. 2, in this embodiment, the microelement division unit 11 of the processor 10 divides the resin molded product 1 into a plurality of microelements. Specifically, the analysis device 100 acquires shape data of the resin molded product 1 by acquiring three-dimensional CAD information of the resin molded product 1 based on an input operation by an operator who analyzes the resin molded product 1. The microelement division unit 11 then divides the three-dimensional shape data of the resin molded product 1 into a plurality of microelements (meshes) to generate an analysis model of the resin molded product 1. The microelement division unit 11 also divides the resin molded product 1 in the thickness direction of each part of the resin molded product 1 so that there is one microelement. Each microelement has a rectangular parallelepiped shape.

たとえば、樹脂成形品1の各部分のうちXY平面に沿って延びるような形状の部分(図2の部分a)では、厚み方向はZ方向となる。したがって、微小要素分割部11は、Z方向における個数が少なくとも1つ(1層)となるように樹脂成形品1を複数の微小要素に分割する。なお、本実施形態の樹脂成形品1は、1つのスプルー部Sから樹脂材料が流動されることにより成形される。したがって、部分aでは、樹脂材料は、X1方向に沿うように流動しながら、樹脂成形品1を成形する金型内に充填される。 For example, in the portion of the resin molded product 1 that is shaped to extend along the XY plane (part a in Figure 2), the thickness direction is the Z direction. Therefore, the microelement division unit 11 divides the resin molded product 1 into multiple microelements so that the number in the Z direction is at least one (one layer). Note that the resin molded product 1 of this embodiment is formed by flowing the resin material from one sprue part S. Therefore, in part a, the resin material flows along the X1 direction as it fills the mold that molds the resin molded product 1.

なお、樹脂成形品1の各部において、樹脂の流動方向をMD(Machine Direction)、樹脂成形品1の表面に沿う面内(厚み方向に対して直交する面内)においてMDに直交する直交方向をTD(Transverse Direction)とする。たとえば、図2の部分aでは、MDは図中のX方向であり、TDは図中のY方向となる。また、プロセッサ10によって、取得された樹脂成形品1の形状データに基づいて、複数の微小要素の各々においてMDおよびTDと、厚み方向とが設定される。 In each part of the resin molded product 1, the flow direction of the resin is defined as MD (Machine Direction), and the orthogonal direction perpendicular to the MD in a plane along the surface of the resin molded product 1 (in a plane perpendicular to the thickness direction) is defined as TD (Transverse Direction). For example, in part a of Figure 2, MD is the X direction in the figure, and TD is the Y direction in the figure. In addition, the processor 10 sets the MD and TD and the thickness direction in each of the multiple microelements based on the acquired shape data of the resin molded product 1.

そして、図3に示すように、本実施形態では、プロセッサ10の仮想節点設定部12は、微小要素分割部11によって分割された複数の微小要素の各々において、樹脂成形品1の厚み方向に沿って複数の仮想節点Pを設定する。具体的には、仮想節点設定部12は、樹脂成形品1の複数の微小要素の各々において、樹脂成形品1の厚み方向に沿って偶数個(例えば20個)の仮想節点Pを設定する。たとえば、図2の部分aに含まれる複数の微小要素のうちの1つの微小要素1aでは、仮想節点設定部12によって、Z方向に沿うように20個の仮想節点P~P20が設定される。ここで、仮想節点Pは、樹脂成形品1においてZ2方向側の表面に位置するように設定される。そして、仮想節点P20は、樹脂成形品1においてZ1方向側の表面に位置するように設定される。 3, in this embodiment, the virtual node setting unit 12 of the processor 10 sets a plurality of virtual nodes P along the thickness direction of the resin molded product 1 in each of the plurality of microelements divided by the microelement dividing unit 11. Specifically, the virtual node setting unit 12 sets an even number (e.g., 20) of virtual nodes P along the thickness direction of the resin molded product 1 in each of the plurality of microelements of the resin molded product 1. For example, in one microelement 1a of the plurality of microelements included in the portion a of FIG. 2, 20 virtual nodes P 1 to P 20 are set by the virtual node setting unit 12 along the Z direction. Here, the virtual node P 1 is set so as to be located on the surface of the resin molded product 1 on the Z2 direction side. And the virtual node P 20 is set so as to be located on the surface of the resin molded product 1 on the Z1 direction side.

〈繊維配向(解析配向)の取得〉
本実施形態では、プロセッサ10の解析配向取得部13は、仮想節点設定部12によって設定された各微小要素の複数の仮想節点Pの各々において、樹脂材料に含まれる繊維の配向を解析する。具体的には、解析配向取得部13は、プログラム21を実行することによって、記憶部20に記憶された解析用データ22に基づいて、複数の仮想節点Pの各々における樹脂材料の流動状態をシミュレートする。そして、解析配向取得部13は、複数の仮想節点Pの各々において、3次元繊維配向テンソルを算出する。MD、TD、および、厚み方向における繊維の配向は、3次元繊維配向テンソルより算出する。詳細には、解析配向取得部13は、プログラム21を実行することによって、射出成形中の樹脂温度、樹脂圧力、樹脂材料の流動速度などの流動データを算出することによって、複数の仮想節点Pの各々において繊維の配向を算出(解析)する。そして、解析配向取得部13は、解析された繊維配向を解析配向Q(図5参照)として取得する。解析された繊維配向である解析配向Qは、実際の樹脂成形品1における繊維配向がシミュレートされたものである。
Obtaining fiber orientation (analysis orientation)
In this embodiment, the analytical orientation acquisition unit 13 of the processor 10 analyzes the orientation of fibers contained in the resin material at each of the multiple virtual nodes P of each microelement set by the virtual node setting unit 12. Specifically, the analytical orientation acquisition unit 13 executes the program 21 to simulate the flow state of the resin material at each of the multiple virtual nodes P based on the analysis data 22 stored in the storage unit 20. Then, the analytical orientation acquisition unit 13 calculates a three-dimensional fiber orientation tensor at each of the multiple virtual nodes P. The fiber orientations in the MD, TD, and thickness directions are calculated from the three-dimensional fiber orientation tensor. In detail, the analytical orientation acquisition unit 13 executes the program 21 to calculate flow data such as the resin temperature, resin pressure, and flow speed of the resin material during injection molding, thereby calculating (analyzing) the fiber orientation at each of the multiple virtual nodes P. Then, the analytical orientation acquisition unit 13 acquires the analyzed fiber orientation as an analytical orientation Q (see FIG. 5). The analyzed fiber orientation, ie, analyzed orientation Q, is a simulation of the fiber orientation in the actual resin molded product 1.

なお、図4に示すように、繊維を含む樹脂成形品1は、樹脂材料の流動方向(MD)、直交方向(TD)、および、厚み方向において、繊維の配向に偏りが生じる。樹脂成形品1に含まれる繊維(強化繊維)は、細長形状を有しており、金型内を流動する際に樹脂材料の剪断力によりMDに配向する比率が高くなる。また、樹脂成形品1において、厚み方向の表層側(表面側)では、MDの配向が大きく(強く)なるとともに、厚み方向における中心側(内側)では、繊維の向き(繊維配向)がランダムな配向となる。たとえば、図3の微小要素1aにおける繊維の配向(繊維配向)では、表層側(表面付近側)においてMDであるX方向に沿うように配置されている繊維の割合が比較的多くなる。 As shown in FIG. 4, the fiber-containing resin molded product 1 has biased fiber orientation in the flow direction (MD), transverse direction (TD), and thickness direction of the resin material. The fibers (reinforcing fibers) contained in the resin molded product 1 have an elongated shape, and the proportion of fibers oriented in the MD is high due to the shear force of the resin material when it flows through the mold. In addition, in the resin molded product 1, the MD orientation is large (strong) on the surface side (surface side) in the thickness direction, and the fiber orientation (fiber orientation) is random on the center side (inside) in the thickness direction. For example, in the fiber orientation (fiber orientation) of the microelement 1a in FIG. 3, the proportion of fibers arranged along the X direction, which is the MD, is relatively high on the surface side (near the surface).

図5に示すように、本実施形態では、解析配向取得部13により取得される解析配向Qは、3次元繊維配向テンソルで表される。MD、TD、および、厚み方向の繊維配向ベクトルは、3次元繊維配向テンソルから計算される。図5の縦軸の値は、樹脂材料における繊維配向ベクトル(MD)の大きさの値を示している。繊維配向テンソルでは、数値が1.0に近づくほどMD(X方向)に沿うように配置されている繊維の割合(確率)が大きいことを示している。そして、数値が0.0に近づくほど、TD(Y方向)、もしくは厚み方向に沿うように配置されている繊維の割合(確率)が大きいことを示している。 As shown in FIG. 5, in this embodiment, the analytical orientation Q acquired by the analytical orientation acquisition unit 13 is expressed as a three-dimensional fiber orientation tensor. The fiber orientation vectors in the MD, TD, and thickness directions are calculated from the three-dimensional fiber orientation tensor. The value on the vertical axis in FIG. 5 indicates the magnitude of the fiber orientation vector (MD) in the resin material. In the fiber orientation tensor, the closer the value is to 1.0, the greater the proportion (probability) of fibers that are arranged along the MD (X direction). And the closer the value is to 0.0, the greater the proportion (probability) of fibers that are arranged along the TD (Y direction) or thickness direction.

〈解析配向の補正〉
ここで、樹脂成形品1の成形において、樹脂成形品1の樹脂材料は、金型内において厚み方向の表層側から順に徐々に冷却されて固化していく。したがって、樹脂成形品1の成形において、表層側の部分が先に冷却されて固化することに起因して、厚み方向における中心側(内側)の部分は、表層側の部分からの影響を受けながら収縮する。そのため、繊維配向(解析配向Q)に基づいて樹脂成形品1の収縮量を解析する場合には、樹脂成形品1の各部における繊維配向を解析することのみならず、表層側から順に樹脂材料が固化することによる影響を考慮して、収縮量を解析する必要がある。
<Analysis Orientation Correction>
Here, in the molding of the resin molded product 1, the resin material of the resin molded product 1 is gradually cooled and solidified in the mold starting from the surface layer side in the thickness direction. Therefore, in the molding of the resin molded product 1, the surface layer side is cooled and solidified first, and the center side (inner side) part in the thickness direction shrinks while being influenced by the surface layer side part. Therefore, when analyzing the shrinkage amount of the resin molded product 1 based on the fiber orientation (analysis orientation Q), it is necessary to not only analyze the fiber orientation in each part of the resin molded product 1, but also to analyze the shrinkage amount taking into account the influence of the resin material solidifying in order from the surface layer side.

そこで、図5および図6に示すように、プロセッサ10の補正配向算出部14は、樹脂成形品1の収縮量を解析するために、解析配向取得部13によって取得された解析配向Qに対して補正を行うことによって、補正された繊維配向である補正配向Rを複数の仮想節点Pの各々において算出するように構成されている。具体的には、本実施形態では、補正配向算出部14は、各微小要素の複数の仮想節点Pの各々において、樹脂成形品1の表層側から順に解析配向Qを補正する計算を行うことによって、収縮量を解析するための補正配向Rを算出する。また、本実施形態では、補正配向算出部14は、各微小要素の複数の仮想節点Pの各々において、自身よりも表層側の仮想節点Pにおける補正結果に基づいて解析配向Qを補正する計算を行う。 As shown in FIG. 5 and FIG. 6, the corrected orientation calculation unit 14 of the processor 10 is configured to calculate the corrected fiber orientation, that is, the corrected orientation R, at each of the multiple virtual nodes P by correcting the analytical orientation Q acquired by the analytical orientation acquisition unit 13 in order to analyze the amount of shrinkage of the resin molded product 1. Specifically, in this embodiment, the corrected orientation calculation unit 14 calculates the corrected orientation R for analyzing the amount of shrinkage at each of the multiple virtual nodes P of each microelement by performing a calculation to correct the analytical orientation Q in sequence from the surface side of the resin molded product 1. Also, in this embodiment, the corrected orientation calculation unit 14 performs a calculation to correct the analytical orientation Q at each of the multiple virtual nodes P of each microelement based on the correction result at the virtual node P on the surface side of the microelement.

図6に示すように、一例として、図3の微小要素1aにおいて、表層からn番目の仮想節点Pにおいて取得された解析配向Qを補正することによって補正配向Rを算出する場合について説明する。なお、n番目の仮想節点Pは、Z2方向から数えて1番目から10番目までの仮想節点P~P10のうちの1つである。すなわち、n-1番目の仮想節点Pn-1は、n番目の仮想節点Pの表層側に隣接する。また、n-1番目の仮想節点Pn-1における解析配向Qを解析配向Qn-1、補正配向Rを補正配向Rn-1とする。なお、仮想節点Pおよび仮想節点Pn-1は、それぞれ、特許請求の範囲における「対象仮想節点」および「隣接仮想節点」の一例である。そして、解析配向Qは、特許請求の範囲における「対象解析配向」の一例である。また、補正配向Rおよび補正配向Rn-1は、それぞれ、特許請求の範囲における「対象補正配向」および「隣接補正配向」の一例である。 As shown in FIG. 6, as an example, in the infinitesimal element 1a of FIG. 3, a case will be described in which the analytical orientation Q n acquired at the n-th virtual node P n from the surface layer is corrected to calculate the corrected orientation R n . The n-th virtual node P n is one of the first to tenth virtual nodes P 1 to P 10 counting from the Z2 direction. That is, the n-1th virtual node P n-1 is adjacent to the surface layer side of the n-th virtual node P n . The analytical orientation Q at the n-1th virtual node P n-1 is the analytical orientation Q n-1 , and the corrected orientation R is the corrected orientation R n-1 . The virtual node P n and the virtual node P n-1 are examples of the "target virtual node" and the "adjacent virtual node" in the claims, respectively. The analytical orientation Q n is an example of the "target analytical orientation" in the claims. Moreover, the corrective orientation R n and the corrective orientation R n-1 are examples of the "target corrective orientation" and the "adjacent corrective orientation" in the claims, respectively.

補正配向算出部14は、仮想節点Pにおける補正配向Rを算出する場合には、仮想節点P自身よりも表層側に隣接する仮想節点Pn-1における補正結果である補正配向Rn-1に基づいて、仮想節点Pにおける解析配向Qを補正する計算を行う。仮想節点Pにおける補正配向Rの算出は、補正係数kと、表層側に隣接する仮想節点Pn-1における補正配向Rn-1とによって、式(1)のように表される。

Figure 0007645207000001

すなわち、本実施形態では、補正配向算出部14は、解析配向Qと、補正係数kによって重み付けされた補正配向Rn-1との平均を計算することによって、解析配向Qを補正する計算を行うことによって、補正配向Rを算出する。 When calculating the correction orientation R n at the virtual node P n , the correction orientation calculation unit 14 performs a calculation to correct the analytical orientation Q n at the virtual node P n based on the correction orientation R n -1, which is the correction result at the virtual node P n-1 adjacent to the virtual node P n on the surface side of the virtual node P n itself. The calculation of the correction orientation R n at the virtual node P n is expressed as shown in formula (1) by the correction coefficient k and the correction orientation R n-1 at the virtual node P n-1 adjacent to the surface side.
Figure 0007645207000001

That is, in this embodiment, the corrected orientation calculation unit 14 calculates the average of the analytical orientation Q n and the corrected orientation R n-1 weighted by the correction coefficient k, thereby performing a calculation to correct the analytical orientation Q n , thereby calculating the corrected orientation R n .

なお、本実施形態では、補正配向算出部14は、解析配向取得部13による流動データ(樹脂材料の解析結果)に基づいて、仮想節点Pn-1と仮想節点Pとの温度差と、仮想節点Pn-1における樹脂材料の弾性率とに基づいて補正係数kを算出するように構成されている。具体的には、本実施形態では、補正係数kは、式(2)の計算によって算出される。

Figure 0007645207000002

固化温度弾性率は、樹脂材料の固化温度における弾性率を示す。また、隣接層弾性率は、樹脂材料の固化温度と隣接層温度差との差の温度における樹脂材料の弾性率を示す。隣接層温度差は、複数の仮想節点Pのうちの隣り合う2つの仮想節点P同士の間の温度差を示す。すなわち、隣接層温度差は、仮想節点Pと仮想節点Pn-1との間の温度差を示す。したがって、補正係数kは、仮想節点Pにおいて樹脂材料が固化する時点での弾性率と、その時点においてすでに固化している仮想節点Pn-1における樹脂材料の弾性率との比率を模擬的に示した値である。 In this embodiment, the correction orientation calculation unit 14 is configured to calculate the correction coefficient k based on the temperature difference between the virtual node P n-1 and the virtual node P n and the elastic modulus of the resin material at the virtual node P n-1 , based on the flow data (resin material analysis results) by the analysis orientation acquisition unit 13. Specifically, in this embodiment, the correction coefficient k is calculated by the calculation of formula (2).
Figure 0007645207000002

The solidification temperature elastic modulus indicates the elastic modulus at the solidification temperature of the resin material. The adjacent layer elastic modulus indicates the elastic modulus of the resin material at the temperature difference between the solidification temperature of the resin material and the adjacent layer temperature difference. The adjacent layer temperature difference indicates the temperature difference between two adjacent virtual nodes P among the multiple virtual nodes P. In other words, the adjacent layer temperature difference indicates the temperature difference between virtual node P n and virtual node P n-1 . Therefore, the correction coefficient k is a value that simulates the ratio between the elastic modulus at the time when the resin material at virtual node P n solidifies and the elastic modulus of the resin material at virtual node P n-1 that has already solidified at that time.

詳細には、解析配向取得部13は、プログラム21を実行することによって、複数の仮想節点Pの各々において、樹脂成形品1の射出成形をシミュレートすることにより成形プロセス中において、保圧が完了した時点、充填工程から保圧工程へと切り替えるVP(velocity pressure)が切り替えられた時点、または、離型した時点、のいずれかにおける樹脂温度を解析(算出)する。そして、複数の仮想節点Pの各々において算出された樹脂温度に基づいて、仮想節点Pと仮想節点Pn-1との間の温度差が隣接層温度差として取得される。また、固化温度は、樹脂材料の特有の固有値であって、記憶部20に記憶されている解析用データ22に含まれる材料物性データに基づいて取得される。また、樹脂材料の弾性率は、同様に、記憶部20に記憶されている解析用データ22に含まれる材料物性データに基づいて算出される。記憶部20に記憶されている解析用データ22には、樹脂成形品1の樹脂材料の弾性率を算出するための条件式が含まれている。なお、樹脂成形品1を形成する樹脂材料の弾性率は、温度の低下に比例して略線形に増加するように値が変化する。したがって、補正係数kは、1よりも大きい値となり、仮想節点Pと仮想節点Pn-1との温度差が大きいほど(弾性率の差が大きいほど)大きい値となる。このように、補正係数kは、仮想節点Pと仮想節点Pn-1との温度差と、隣接層(仮想節点Pn-1)の弾性率を考慮した大きさとなるように算出される。 In detail, the analysis orientation acquisition unit 13 executes the program 21 to simulate the injection molding of the resin molded product 1 at each of the multiple virtual nodes P, thereby analyzing (calculating) the resin temperature at any of the following points during the molding process: when the pressure holding is completed, when the velocity pressure (VP) is switched from the filling process to the pressure holding process, or when the mold is released. Then, based on the resin temperature calculated at each of the multiple virtual nodes P, the temperature difference between the virtual node P n and the virtual node P n-1 is acquired as the adjacent layer temperature difference. In addition, the solidification temperature is a unique value specific to the resin material, and is acquired based on the material property data included in the analysis data 22 stored in the storage unit 20. In addition, the elastic modulus of the resin material is similarly calculated based on the material property data included in the analysis data 22 stored in the storage unit 20. The analysis data 22 stored in the storage unit 20 includes a conditional formula for calculating the elastic modulus of the resin material of the resin molded product 1. The elastic modulus of the resin material forming the resin molded product 1 changes in value so as to increase approximately linearly in proportion to a decrease in temperature. Therefore, the correction coefficient k is a value greater than 1, and increases as the temperature difference between virtual node P n and virtual node P n-1 increases (as the difference in elastic modulus increases). In this way, the correction coefficient k is calculated to have a magnitude that takes into account the temperature difference between virtual node P n and virtual node P n-1 and the elastic modulus of the adjacent layer (virtual node P n-1 ).

なお、仮想節点Pは、仮想節点P~P10のうちの1つとしたが、仮想節点P11~P20においても、仮想節点Pと同様の計算が行われる。その場合には、表層側に隣接する仮想節点は仮想節点Pn+1となる。そして、本実施形態では、補正配向算出部14は、仮想節点Pにおける補正配向Rの算出と同様の計算を、各微小要素ごとの複数の仮想節点Pの各々において、樹脂成形品1の表層側から順に繰り返すことによって、複数の仮想節点Pの各々において補正配向Rを算出するように構成されている。たとえば、補正配向算出部14は、微小要素1aに含まれる20個の仮想節点P~P20に対して、仮想節点P、P、P、・・・、P10の順に解析配向Q~Q10を補正する計算を行うことによって補正配向R~R10を算出するとともに、仮想節点P20、P19、P18、・・・、P11の順に解析配向Q20~Q11を補正する計算を行うことによって、補正配向R20~R11を算出する。なお、仮想節点Pおよび仮想節点Pは、一番表層側であるため、解析配向QおよびQ20が、それぞれ、そのまま補正配向RおよびR20として算出される。そして、表層側から2番目以降の仮想節点Pにおいて、式(1)に示した計算によって、補正配向Rが算出される。 Although the virtual node Pn is one of the virtual nodes P1 to P10 , the same calculation as for the virtual node Pn is performed for the virtual nodes P11 to P20 . In that case, the virtual node adjacent to the surface side is the virtual node Pn +1 . In this embodiment, the corrected orientation calculation unit 14 is configured to calculate the corrected orientation R at each of the multiple virtual nodes P by repeating a calculation similar to the calculation of the corrected orientation Rn at the virtual node Pn for each of the multiple virtual nodes P for each infinitesimal element, in order from the surface side of the resin molded product 1. For example, the correction orientation calculation unit 14 calculates the correction orientations R 1 to R 10 by performing calculations to correct the analytical orientations Q 1 to Q 10 in the order of virtual nodes P 1 , P 2 , P 3 , ..., P 10 for the 20 virtual nodes P 1 to P 20 included in the infinitesimal element 1a, and calculates the correction orientations R 20 to R 11 by performing calculations to correct the analytical orientations Q 20 to Q 11 in the order of virtual nodes P 20 , P 19 , P 18 , ..., P 11. Since the virtual nodes P 1 and P 2 are on the outermost layer side, the analytical orientations Q 1 and Q 20 are calculated as the correction orientations R 1 and R 20 as they are, respectively. Then, the correction orientation R is calculated by the calculation shown in formula (1) for the virtual nodes P that are second or later from the outermost layer side.

〈収縮量の解析〉
そして、本実施形態では、プロセッサ10の収縮量解析部15は、各微小要素の複数の仮想節点Pの各々において算出された補正配向Rの平均に基づいて、複数の微小要素の各々の収縮量を解析する。また、収縮量解析部15は、複数の微小要素の各々において、樹脂材料の流動方向(MD)、直交方向(TD)、および、厚み方向のうちの少なくとも2方向における収縮量をそれぞれ解析するように構成されている。
Analysis of shrinkage amount
In this embodiment, the shrinkage amount analysis unit 15 of the processor 10 analyzes the shrinkage amount of each of the multiple microelements based on the average of the corrected orientations R calculated at each of the multiple virtual nodes P of each microelement. The shrinkage amount analysis unit 15 is also configured to analyze the shrinkage amount in at least two directions of the flow direction (MD), the orthogonal direction (TD), and the thickness direction of the resin material for each of the multiple microelements.

たとえば、収縮量解析部15は、微小要素1aにおける20個の仮想節点P~P20の各々における補正配向R~R20の平均値(平均の繊維配向テンソル)を算出する。そして、収縮量解析部15は、算出された平均値を流動方向(MD)と、直交方向(TD)との収縮量(収縮率)の比率として微小要素1aにおける収縮量を算出する。たとえば、算出された補正配向R~R20の平均値が0.8である場合には、MDにおける収縮量と、TDにおける収縮量との比率が、0.8:0.2となるように微小要素1aにおける収縮量が算出される。 For example, the shrinkage amount analysis unit 15 calculates the average value (average fiber orientation tensor) of the corrected orientations R 1 to R 20 at each of the 20 virtual nodes P 1 to P 20 in the infinitesimal element 1a. The shrinkage amount analysis unit 15 then calculates the shrinkage amount in the infinitesimal element 1a by taking the calculated average value as the ratio of the shrinkage amount (shrinkage rate) in the flow direction (MD) to the orthogonal direction (TD). For example, when the average value of the calculated corrected orientations R 1 to R 20 is 0.8, the shrinkage amount in the infinitesimal element 1a is calculated so that the ratio of the shrinkage amount in the MD to the shrinkage amount in the TD is 0.8:0.2.

そして、収縮量解析部15(プロセッサ10)は、解析された複数の微小要素の各々における流動方向(MD)と直交方向(TD)との収縮量に基づいて、樹脂成形品1の流動方向(MD)と直交方向(TD)と収縮量を解析する。上記のようにして、解析装置100のプロセッサ10(収縮量解析部15)は、樹脂材料の流動による繊維配向と、表層側から順に樹脂材料が固化することによる影響(弾性率の差異)とを考慮した樹脂成形品1の収縮量を解析することによって、樹脂成形品1の形状(そり変形)の解析を行う。 Then, the shrinkage amount analysis unit 15 (processor 10) analyzes the flow direction (MD) and the orthogonal direction (TD) and the shrinkage amount of the resin molded product 1 based on the shrinkage amount in the flow direction (MD) and the orthogonal direction (TD) of each of the analyzed microelements. In this manner, the processor 10 (shrinkage amount analysis unit 15) of the analysis device 100 analyzes the shape (warpage deformation) of the resin molded product 1 by analyzing the shrinkage amount of the resin molded product 1 taking into account the fiber orientation due to the flow of the resin material and the effect (difference in elastic modulus) due to the solidification of the resin material from the surface side.

また、プロセッサ10は、樹脂成形品1の解析結果を表示部30に表示させるように構成されている。たとえば、プロセッサ10は、解析結果として、解析後の形状(変形後の形状)の樹脂成形品1を表示部30に表示させる。また、プロセッサ10は、補正配向算出部14によって算出された補正配向Rを視認可能なように表示部30に表示させるように構成されている。 The processor 10 is also configured to display the analysis results of the resin molded product 1 on the display unit 30. For example, the processor 10 displays the post-analysis shape (post-deformation shape) of the resin molded product 1 on the display unit 30 as the analysis result. The processor 10 is also configured to display the corrected orientation R calculated by the corrected orientation calculation unit 14 on the display unit 30 so that it is visible.

(樹脂成形品解析方法)
次に、図7を参照して、本実施形態の解析装置100による樹脂成形品解析方法について説明する。本実施形態の樹脂成形品解析方法は、繊維を含む樹脂材料によって成形される樹脂成形品1の収縮量を解析する方法である。また、本実施形態による樹脂成形品解析方法は、解析装置100のプロセッサ10にプログラム21を実行させることにより実施される。
(Resin molded product analysis method)
Next, a resin molded product analysis method using the analysis device 100 of this embodiment will be described with reference to Fig. 7. The resin molded product analysis method of this embodiment is a method for analyzing the amount of shrinkage of a resin molded product 1 that is molded from a resin material that contains fibers. The resin molded product analysis method of this embodiment is implemented by causing the processor 10 of the analysis device 100 to execute a program 21.

まず、ステップS1において、樹脂成形品1の形状データが取得されるとともに、収縮量の解析を行うために用いられる各種の解析用データ22が取得(記憶)される。 First, in step S1, shape data of the resin molded product 1 is acquired, and various analysis data 22 used to analyze the amount of shrinkage are acquired (stored).

次に、ステップS2において、樹脂成形品1が複数の微小要素(メッシュ)に分割される。具体的には、樹脂成形品1の形状データが複数の微小要素に分割されることによって、樹脂成形品1の解析モデルが生成される。 Next, in step S2, the resin molded product 1 is divided into a plurality of minute elements (meshes). Specifically, the shape data of the resin molded product 1 is divided into a plurality of minute elements, and an analytical model of the resin molded product 1 is generated.

次に、ステップS3において、複数の微小要素の各々において、樹脂成形品1の厚み方向に沿って複数の仮想節点Pが設定される。たとえば、微小要素1aにおいて、偶数個(20個)の仮想節点P~P20が設定される。 Next, in step S3, in each of the plurality of infinitesimal elements, a plurality of virtual nodes P are set along the thickness direction of the resin molded product 1. For example, in the infinitesimal element 1a, an even number (20) of virtual nodes P 1 to P 20 are set.

次に、ステップS4において、各微小要素の複数の仮想節点Pの各々において、繊維の配向を解析することによって、解析配向Qが取得される。たとえば、微小要素1aにおいて、20個の仮想節点P~P20のうちの1の仮想節点Pにおける解析配向Qが取得される。 Next, in step S4, an analytical orientation Q is obtained by analyzing the orientation of the fiber at each of the multiple virtual nodes P of each microelement. For example, in the microelement 1a, an analytical orientation Qn is obtained at one virtual node Pn among the 20 virtual nodes P1 to P20 .

次に、ステップS5において、各微小要素の複数の仮想節点Pの各々において、樹脂成形品1の表層側から順に解析配向Qを補正する計算を行うことによって、収縮量を解析するための補正配向Rが算出される。具体的には、各微小要素の複数の仮想節点Pの各々において、自身よりも表層側の仮想節点Pにおける補正結果に基づいて解析配向Qを補正する計算を行うことによって、補正配向Rが算出される。なお、解析配向Qの補正方法(補正配向Rの算出方法)の詳細は後述する。 Next, in step S5, a calculation is performed to correct the analytical orientation Q at each of the multiple virtual nodes P of each microelement, starting from the surface side of the resin molded product 1, thereby calculating a corrected orientation R for analyzing the amount of shrinkage. Specifically, a calculation is performed to correct the analytical orientation Q at each of the multiple virtual nodes P of each microelement, based on the correction result at the virtual node P on the surface side of the microelement, thereby calculating the corrected orientation R. The method of correcting the analytical orientation Q (method of calculating the corrected orientation R) will be described in detail later.

次に、ステップS6において、各微小要素の複数の仮想節点Pの各々において算出された補正配向Rの平均に基づいて、複数の微小要素の各々の流動方向(MD)と直交方向(TD)とにおける収縮量が解析される。そして、解析された複数の微小要素の各々におけるMDとTDとの収縮量に基づいて、樹脂成形品1のMDとTDとの収縮量が解析される。 Next, in step S6, the amount of shrinkage in the flow direction (MD) and orthogonal direction (TD) of each of the microelements is analyzed based on the average of the corrected orientations R calculated at each of the multiple virtual nodes P of each microelement. Then, the amount of shrinkage in the MD and TD of the resin molded product 1 is analyzed based on the amount of shrinkage in the MD and TD of each of the analyzed multiple microelements.

次に、ステップS7において、樹脂成形品1の収縮量の解析結果が表示部30に表示される。具体的には、解析後の形状(変形後の形状)の樹脂成形品1が表示部30に表示される。また、算出された補正配向Rが視認可能なように表示部30に表示される。 Next, in step S7, the analysis result of the amount of shrinkage of the resin molded product 1 is displayed on the display unit 30. Specifically, the shape of the resin molded product 1 after analysis (shape after deformation) is displayed on the display unit 30. In addition, the calculated corrected orientation R is displayed on the display unit 30 so that it can be seen.

〈解析配向の補正方法〉
次に、図8を参照して、解析配向Qの補正方法の詳細について説明する。一例として、複数の微小要素のうちの一の微小要素1a(図3参照)における仮想節点Pの解析配向Qの補正方法について説明する。なお、n=1~20である。
<Method of correcting analytical orientation>
Next, the details of the correction method of the analytical orientation Q will be described with reference to Fig. 8. As an example, the correction method of the analytical orientation Qn of the virtual node Pn in one infinitesimal element 1a (see Fig. 3) among the infinitesimal elements will be described. Here, n = 1 to 20.

まず、ステップS51において、微小要素1aの20個の仮想節点P~P20のうちの一番表層側(最も外側)の仮想節点PおよびP20において(n=1および20の場合)、取得された解析配向QおよびQ20が、そのまま、それぞれ補正配向RおよびR20として取得される。 First, in step S51, the analytical orientations Q1 and Q20 obtained at the virtual nodes P1 and P20 on the outermost side (outermost) of the 20 virtual nodes P1 to P20 of the infinitesimal element 1a (when n = 1 and 20 ) are directly obtained as corrective orientations R1 and R20 , respectively.

次に、ステップS52において、n=2として、20個の仮想節点P~P20のうちの表層側から2番目の仮想節点Pにおいて、仮想節点Pの表層側に隣接する仮想節点Pn-1における樹脂材料の弾性率と、仮想節点Pと仮想節点Pn-1との間の温度差とに基づいて、補正係数kが算出される。同様にn=19として、仮想節点Pにおける補正係数kが算出される。なお、n=11~19の場合には、隣接する仮想節点は仮想節点Pn+1となる。 Next, in step S52, for n=2, a correction coefficient k is calculated for the second virtual node Pn from the surface side of the 20 virtual nodes P1 to P20 , based on the elastic modulus of the resin material at the virtual node Pn -1 adjacent to the surface side of the virtual node Pn and the temperature difference between the virtual node Pn and the virtual node Pn-1 . Similarly, for n=19, a correction coefficient k for the virtual node Pn is calculated. Note that for n=11 to 19, the adjacent virtual node is the virtual node Pn +1 .

次に、ステップS53において、仮想節点Pn-1において算出された補正配向Rn-1と、補正係数kとに基づいて、解析配向Qを補正する計算が行われる。具体的には、解析配向Qと、算出された補正係数kによって重み付けされた補正配向Rn-1との平均を計算することによって、解析配向Qを補正する計算が行われる。これにより、補正配向Rn-1と、仮想節点Pn-1における樹脂材料の弾性率とに基づいて、解析配向Qを補正する計算が行われることによって補正配向Rが算出される。 Next, in step S53, a calculation is performed to correct the analytical orientation Q n based on the correction orientation R n-1 calculated at the virtual node P n-1 and the correction coefficient k. Specifically, a calculation is performed to correct the analytical orientation Q n by calculating the average of the analytical orientation Q n and the correction orientation R n-1 weighted by the calculated correction coefficient k. As a result, a calculation is performed to correct the analytical orientation Q n based on the correction orientation R n-1 and the elastic modulus of the resin material at the virtual node P n -1 , thereby calculating the correction orientation R n .

次に、ステップS54において、複数の仮想節点Pの各々において、樹脂成形品1の表層側から順に(nの値を変更しながら)、ステップS52およびステップS53が繰り返される。具体的には、仮想節点P以降、仮想節点P、P、・・・P10の順にステップS52およびステップS53が繰り返されることによって、補正配向R~R10が算出される。また、反対側の表層からも同様に、仮想節点P19以降、仮想節点P18、P17、・・・、P11の順にステップS52およびステップS53が繰り返されることによって、補正配向R11~R19が算出される。 Next, in step S54, steps S52 and S53 are repeated for each of the multiple virtual nodes P in order from the surface layer side of the resin molded product 1 (while changing the value of n). Specifically, steps S52 and S53 are repeated for virtual nodes P2 , P3 , P4 , ..., P10 in this order to calculate corrected orientations R2 to R10 . Similarly, steps S52 and S53 are repeated for virtual nodes P19 , P18 , P17 , ..., P11 in the order from the surface layer on the opposite side to calculate corrected orientations R11 to R19 .

上記のように、ステップS51~ステップS54が実行されることによって、微小要素1aにおける20個の仮想節点P(仮想節点P~P20)の各々における解析配向Q(解析配向Q~Q20)が補正されて、補正配向R(補正配向R~R20)が算出される。樹脂成形品1の複数の微小要素の各々において、同様の処理が実行される。 As described above, by executing steps S51 to S54, the analytical orientation Q n (analytical orientation Q 1 to Q 20 ) at each of the 20 virtual nodes P n (virtual nodes P 1 to P 20 ) in the infinitesimal element 1a is corrected, and the corrected orientation R n (corrected orientation R 1 to R 20 ) is calculated. A similar process is executed in each of the multiple infinitesimal elements of the resin molded product 1.

[本実施形態の効果]
本実施形態では、以下のような効果を得ることができる。
[Effects of this embodiment]
In this embodiment, the following effects can be obtained.

本実施形態の樹脂成形品解析方法では、上記のように、各微小要素の複数の仮想節点Pの各々において、樹脂成形品1の表層側から順に解析配向Qを補正する計算を行うことによって、収縮量を解析するための補正配向Rを算出する。そして、各微小要素の複数の仮想節点Pの各々において、自身よりも表層側の仮想節点Pにおける補正結果(補正配向R)に基づいて解析配向Qを補正する計算を行うことによって、補正配向Rを算出する。これにより、各微小要素の複数の仮想節点Pの各々において解析された繊維配向である解析配向Qに対して、樹脂成形品1の表層側から順に、自身よりも表層側の仮想節点Pにおける補正結果に基づいて補正する計算を行うため、表層側から順に樹脂材料が固化することによる影響を考慮して、解析された繊維配向(解析配向Q)を補正することができる。そのため、表層側から順に樹脂材料が固化することによる影響を考慮して補正された繊維配向の分布を合成することによって、樹脂成形品1の収縮量を解析することができる。その結果、樹脂成形品1の厚み方向における繊維配向の分布を合成することによって、合成された繊維配向に基づいて樹脂成形品1の収縮量を解析する場合に、樹脂成形品1の収縮量を正確に解析することができる。 In the resin molded product analysis method of this embodiment, as described above, a calculation is performed to correct the analytical orientation Q at each of the multiple virtual nodes P of each microelement in order from the surface side of the resin molded product 1, thereby calculating the corrected orientation R for analyzing the amount of shrinkage. Then, at each of the multiple virtual nodes P of each microelement, a calculation is performed to correct the analytical orientation Q based on the correction result (corrected orientation R) at the virtual node P on the surface side of the microelement, thereby calculating the corrected orientation R. As a result, a calculation is performed to correct the analytical orientation Q, which is the fiber orientation analyzed at each of the multiple virtual nodes P of each microelement, in order from the surface side of the resin molded product 1, based on the correction result at the virtual node P on the surface side of the microelement, so that the analyzed fiber orientation (analysis orientation Q) can be corrected taking into account the effect of the resin material solidifying in order from the surface side. Therefore, the amount of shrinkage of the resin molded product 1 can be analyzed by synthesizing the distribution of the fiber orientation corrected in consideration of the effect of the resin material solidifying in order from the surface side. As a result, by synthesizing the distribution of fiber orientation in the thickness direction of the resin molded product 1, the amount of shrinkage of the resin molded product 1 can be accurately analyzed when analyzing the amount of shrinkage of the resin molded product 1 based on the synthesized fiber orientation.

また、本実施形態では、上記のように、複数の仮想節点Pのうちの一の仮想節点P(対象仮想節点)において、仮想節点Pにおける解析配向Q(対象解析配向)を取得する。そして、複数の仮想節点Pのうちの仮想節点Pの表層側に隣接する仮想節点Pn-1(隣接仮想節点)において算出された補正配向Rn-1(隣接補正配向)に基づいて、解析配向Qを補正する計算を行うことによって、仮想節点Pにおける補正配向R(対象補正配向)を算出する。そして、補正配向Rの算出を、複数の仮想節点Pの各々において、樹脂成形品1の表層側から順に繰り返す。これにより、生成される樹脂成形品1の厚み方向において設定された複数の仮想節点Pの各々において、表層側に隣り合う仮想節点Pにおける補正結果(補正配向R)を用いて、1層ずつ順に表層側から内側に向かって補正配向Rを算出することができる。そのため、複数の仮想節点Pの各々において、自身よりも表層側の全ての仮想節点Pにおける補正結果を用いる場合に比べて、補正配向Rの算出に用いる仮想節点Pの数を少なくすることができるので、収縮量を解析するための解析配向Qに対する補正を容易に行うことができる。その結果、表層側に隣り合う仮想節点Pにおける補正結果を用いることによって、合成された補正後の繊維配向に基づいて樹脂成形品1の収縮量を容易にかつ正確に解析することができる。 In addition, in this embodiment, as described above, an analytical orientation Q n (target analytical orientation) at one virtual node P n (target virtual node) among the multiple virtual nodes P is obtained. Then, a calculation is performed to correct the analytical orientation Q n based on a corrected orientation R n-1 (adjacent corrected orientation) calculated at a virtual node P n-1 (adjacent virtual node) adjacent to the surface side of the virtual node P n among the multiple virtual nodes P , thereby calculating a corrected orientation R n ( target corrected orientation) at the virtual node P n. Then, the calculation of the corrected orientation R n is repeated in order from the surface side of the resin molded product 1 at each of the multiple virtual nodes P. As a result, in each of the multiple virtual nodes P set in the thickness direction of the generated resin molded product 1, the corrected orientation R can be calculated from the surface side to the inside, layer by layer, using the correction result (corrected orientation R) at the virtual node P adjacent to the surface side. Therefore, for each of the multiple virtual nodes P, the number of virtual nodes P used to calculate the corrected orientation R can be reduced compared to the case where the correction results for all virtual nodes P on the surface layer side of the virtual node P are used, and therefore correction of the analytical orientation Q for analyzing the amount of shrinkage can be easily performed. As a result, by using the correction results for the adjacent virtual nodes P on the surface layer side, the amount of shrinkage of the resin molded product 1 can be easily and accurately analyzed based on the synthesized corrected fiber orientation.

また、本実施形態では、上記のように、補正配向Rn-1(隣接補正配向)と仮想節点Pn-1(隣接仮想節点)における樹脂材料の弾性率とに基づいて、解析配向Q(対象解析配向)を補正する計算を行うことによって、補正配向R(対象補正配向)を算出する。ここで、表層側に隣接する層が先に冷却されて固化した場合には、すでに固化した隣接する層との弾性率の差異に起因して、実際の収縮量と解析された繊維配向に基づいて算出された収縮量とに差が生じると考えられる。これを考慮して、本実施形態では、補正配向Rn-1と仮想節点Pn-1における樹脂材料の弾性率とに基づいて、解析配向Qを補正する計算を行うことによって、補正配向Rを算出する。これにより、表層側に隣り合う仮想節点Pにおける補正配向Rに加えて、表層側に隣り合う仮想節点Pにおける弾性率をも考慮にいれて、解析配向Qを補正することができる。そのため、収縮量を解析するために解析された繊維配向である解析配向Qをより精度よく補正することができる。その結果、補正された繊維配向である補正配向Rに基づいて収縮量を解析する場合に、収縮量をより精度よく正確に解析することができる。 In addition, in this embodiment, as described above, the analytical orientation Q n (target analytical orientation) is corrected based on the correction orientation R n-1 (adjacent correction orientation) and the elastic modulus of the resin material at the virtual node P n -1 (adjacent virtual node), and the correction orientation R n (target correction orientation) is calculated. Here, when the layer adjacent to the surface layer side is cooled and solidified first, it is considered that a difference occurs between the actual shrinkage amount and the shrinkage amount calculated based on the analyzed fiber orientation due to the difference in elastic modulus with the adjacent layer that has already solidified. In consideration of this, in this embodiment, the correction orientation R n is calculated by performing a calculation to correct the analytical orientation Q n based on the correction orientation R n -1 and the elastic modulus of the resin material at the virtual node P n-1 . As a result, in addition to the correction orientation R at the virtual node P adjacent to the surface layer side, the elastic modulus at the virtual node P adjacent to the surface layer side can also be taken into consideration to correct the analytical orientation Q. Therefore, the analytical orientation Q, which is the fiber orientation analyzed to analyze the shrinkage amount, can be corrected more accurately. As a result, when the amount of shrinkage is analyzed based on the corrected fiber orientation R, it is possible to analyze the amount of shrinkage more accurately and precisely.

また、本実施形態では、上記のように、仮想節点Pn-1(隣接仮想節点)における樹脂材料の弾性率と、仮想節点P(対象仮想節点)と仮想節点Pn-1との間の温度差とに基づいて、補正係数kを算出する。そして、補正配向Rn-1(隣接補正配向)と補正係数kとに基づいて、解析配向Q(対象解析配向)を補正する計算を行う。これにより、樹脂材料の弾性率は温度に依存して変化するため、表層側に隣接する仮想節点Pとの温度差と隣接する仮想節点Pにおける弾性率とに基づいて算出された補正係数kを用いることによって、すでに固化した隣接する層との弾性率の差異に起因する影響を考慮しながら、収縮量を解析するために解析配向Qをより一層精度よく補正することができる。その結果、補正された繊維配向である補正配向Rを合成することによって収縮量を解析する場合に、収縮量をより一層精度よく正確に解析することができる。 In addition, in this embodiment, as described above, the correction coefficient k is calculated based on the elastic modulus of the resin material at the virtual node P n -1 (adjacent virtual node) and the temperature difference between the virtual node P n ( target virtual node) and the virtual node P n-1. Then, a calculation is performed to correct the analysis orientation Q n (target analysis orientation) based on the correction orientation R n-1 (adjacent correction orientation) and the correction coefficient k. As a result, since the elastic modulus of the resin material changes depending on the temperature, by using the correction coefficient k calculated based on the temperature difference with the adjacent virtual node P on the surface layer side and the elastic modulus at the adjacent virtual node P, it is possible to more accurately correct the analysis orientation Q to analyze the amount of shrinkage while considering the influence caused by the difference in the elastic modulus with the adjacent layer that has already solidified. As a result, when the amount of shrinkage is analyzed by combining the correction orientation R, which is the corrected fiber orientation, it is possible to more accurately analyze the amount of shrinkage.

また、本実施形態では、上記のように、解析配向Q(対象解析配向)と、補正係数kによって重み付けされた補正配向Rn-1(隣接補正配向)との平均を計算することによって、解析配向Qを補正する計算を行う。これにより、補正係数kによって重み付けされた補正配向Rn-1との平均を計算することによって、解析配向Qを補正する計算をより容易に行うことができる。その結果、補正された解析配向Q(解析配向Q)である補正配向R(補正配向R)に基づいて、収縮量をより容易に解析することができる。 In addition, in this embodiment, as described above, a calculation for correcting the analytical orientation Q n is performed by calculating the average of the analytical orientation Q n (target analytical orientation) and the corrected orientation R n-1 (adjacent corrected orientation) weighted by the correction coefficient k . This makes it easier to perform a calculation for correcting the analytical orientation Q n by calculating the average with the corrected orientation R n-1 weighted by the correction coefficient k. As a result, the amount of shrinkage can be analyzed more easily based on the corrected orientation R (corrected orientation R n ), which is the corrected analytical orientation Q (analysis orientation Q n ).

また、本実施形態では、上記のように、各微小要素の複数の仮想節点Pの各々において算出された補正配向Rの平均に基づいて、複数の微小要素の各々の収縮量を解析するとともに、解析された複数の微小要素の各々における収縮量に基づいて、樹脂成形品1の収縮量を解析する。これにより、仮想節点Pの各々において算出された補正配向Rの平均を取得することによって、厚み方向において合成された繊維配向を容易に取得することができる。そのため、補正配向Rに基づいて、各微小要素ごとに合成された補正配向Rを容易に取得することができるので、各微小要素ごとの収縮量を容易に解析することができる。その結果、表層側から順に樹脂材料が固化することによる影響を考慮して、樹脂成形品1の収縮量を容易に正確に解析することができる。 In addition, in this embodiment, as described above, the amount of shrinkage of each of the multiple microelements is analyzed based on the average of the corrected orientation R calculated at each of the multiple virtual nodes P of each microelement, and the amount of shrinkage of the resin molded product 1 is analyzed based on the amount of shrinkage of each of the analyzed multiple microelements. In this way, by obtaining the average of the corrected orientation R calculated at each of the virtual nodes P, the fiber orientation synthesized in the thickness direction can be easily obtained. Therefore, the corrected orientation R synthesized for each microelement can be easily obtained based on the corrected orientation R, so that the amount of shrinkage for each microelement can be easily analyzed. As a result, the amount of shrinkage of the resin molded product 1 can be easily and accurately analyzed, taking into account the effect of the resin material solidifying in sequence from the surface layer side.

また、本実施形態では、上記のように、複数の微小要素の各々において、樹脂成形品1の厚み方向に沿って偶数個(例えば20個)の複数の仮想節点Pを設定する。これにより、偶数個の仮想節点Pが設定されるため、厚み方向における一方側からの計算と他方側からの計算とを同じ回数ずつ行うことによって、複数の仮想節点Pの全体において補正配向Rの算出を表層側から順に行うことができる。その結果、奇数個の仮想節点Pが設定されている場合とは異なり、余りがでることなく全ての仮想節点Pにおいて一方側と他方側との両方から順に同様の補正の計算を行うことができるので、全ての仮想節点Pにおいて補正の計算をより容易に行うことができる。 In addition, in this embodiment, as described above, an even number (e.g., 20) of virtual nodes P are set along the thickness direction of the resin molded product 1 in each of the multiple infinitesimal elements. As a result, an even number of virtual nodes P are set, and by performing calculations from one side and the other side in the thickness direction an equal number of times, the calculation of the correction orientation R can be performed for all of the multiple virtual nodes P in sequence starting from the surface layer side. As a result, unlike when an odd number of virtual nodes P are set, similar correction calculations can be performed in sequence from both one side and the other side for all of the virtual nodes P without any remainders, making it easier to perform correction calculations for all of the virtual nodes P.

また、本実施形態では、上記のように、各微小要素の複数の仮想節点Pの各々において算出された補正配向Rに基づいて、複数の微小要素の各々において、樹脂材料の流動方向(MD)と、樹脂成形品1の表面に沿う面内において流動方向と直交する直交方向(TD)との収縮量を解析するとともに、解析された複数の微小要素の各々における流動方向と直交方向との収縮量に基づいて、樹脂成形品1の流動方向と直交方向との収縮量を解析する。これにより、樹脂材料の流動方向と直交方向との2つの方向における繊維配向(解析配向Q)を補正する計算を行うため、流動方向および直交方向に厚み方向をも加えた3つの方向における繊維配向の補正を行う場合に比べて、補正の計算に要する処理負担を軽減することができる。その結果、算出された補正配向Rに基づいて流動方向と直交方向との2つの方向における収縮量の解析を行うことによって、樹脂成形品1における収縮量を解析する場合における処理負担を軽減することができる。 In addition, in this embodiment, as described above, based on the correction orientation R calculated at each of the multiple virtual nodes P of each microelement, the shrinkage amount in the flow direction (MD) of the resin material and the orthogonal direction (TD) perpendicular to the flow direction in the plane along the surface of the resin molded product 1 is analyzed for each of the multiple microelements, and the shrinkage amount in the flow direction and orthogonal direction of the resin molded product 1 is analyzed based on the shrinkage amount in the flow direction and orthogonal direction in each of the analyzed multiple microelements. As a result, since a calculation is performed to correct the fiber orientation (analysis orientation Q) in two directions, the flow direction and the orthogonal direction of the resin material, the processing load required for the correction calculation can be reduced compared to the case where the fiber orientation is corrected in three directions, including the flow direction, the orthogonal direction, and the thickness direction. As a result, by analyzing the shrinkage amount in two directions, the flow direction and the orthogonal direction, based on the calculated correction orientation R, the processing load when analyzing the shrinkage amount in the resin molded product 1 can be reduced.

また、本実施形態では、プログラム21、記録媒体101および、解析装置100は、上記の樹脂成形品解析方法と同様に、樹脂成形品1の厚み方向における繊維配向の分布を合成することによって、合成された繊維配向に基づいて樹脂成形品1の収縮量を解析する場合に、樹脂成形品1の収縮量を正確に解析することができる。 In addition, in this embodiment, the program 21, recording medium 101, and analysis device 100 can accurately analyze the amount of shrinkage of the resin molded product 1 when analyzing the amount of shrinkage of the resin molded product 1 based on the synthesized fiber orientation by synthesizing the distribution of fiber orientation in the thickness direction of the resin molded product 1, similar to the above-mentioned resin molded product analysis method.

[変形例]
なお、今回開示された実施形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施形態の説明ではなく特許請求の範囲によって示され、さらに特許請求の範囲と均等の意味および範囲内でのすべての変更(変形例)が含まれる。
[Modification]
It should be noted that the embodiments disclosed herein are illustrative and not restrictive in all respects. The scope of the present invention is indicated by the claims, not by the description of the embodiments above, and further includes all modifications (variations) within the meaning and scope of the claims.

たとえば、上記実施形態では、複数の微小要素の各々において、樹脂成形品1の厚み方向に沿って複数の仮想節点Pを設定する例を示したが、本発明はこれに限られない。本発明では、図9に示す変形例による仮想節点Paのように、複数の微小要素の頂点を構成する複数の節点201の各々において、樹脂成形品1厚み方向(Z方向)に沿って複数の仮想節点Paを設定するようにしてもよい。 For example, in the above embodiment, an example was shown in which multiple virtual nodes P are set along the thickness direction of the resin molded product 1 in each of multiple microelements, but the present invention is not limited to this. In the present invention, multiple virtual nodes Pa may be set along the thickness direction (Z direction) of the resin molded product 1 in each of multiple nodes 201 that constitute the vertices of multiple microelements, such as the virtual nodes Pa in the modified example shown in Figure 9.

また、上記実施形態では、仮想節点P(対象仮想節点)における解析配向Q(対象解析配向)を補正する場合に、表層側に隣接する仮想節点Pn-1(隣接仮想節点)において算出された補正配向Rn-1(隣接補正配向)に基づいて補正の計算が行われる例を示したが、本発明はこれに限られない。本発明では、表層側に隣接する仮想節点Pのみならず、自身より表層側の複数の仮想節点Pにおける補正結果(補正配向R)に基づいて、解析配向Qを補正する計算を実行するようにしてもよい。 In the above embodiment, when correcting the analytical orientation Q n (target analytical orientation) at the virtual node P n (target virtual node), the correction calculation is performed based on the correction orientation R n -1 (adjacent correction orientation) calculated at the adjacent virtual node P n-1 (adjacent virtual node) on the surface side, but the present invention is not limited to this. In the present invention, the calculation for correcting the analytical orientation Q may be performed based on the correction results (correction orientations R) not only at the virtual node P adjacent to the surface side but also at multiple virtual nodes P on the surface side.

また、上記実施形態では、仮想節点P(対象仮想節点)と仮想節点Pn-1(隣接仮想節点)との温度差に基づいて算出された仮想節点Pn-1における樹脂材料の弾性率に基づいて補正係数kを算出する例を示したが、本発明はこれに限られない。たとえば、補正係数kを、弾性率を算出せず、仮想節点Pと仮想節点Pn-1との温度差に基づいて算出するようにしてもよい。また、補正係数kを、複数の仮想節点Pの各々ごとに算出せずに、1つの微小要素において共通の値を用いるようにしてもよい。また、補正係数kを、樹脂成形品1の全体において、樹脂材料の特性ごとに予め設定された一定の固有値を用いるようにしてもよい。 In the above embodiment, the correction coefficient k is calculated based on the elastic modulus of the resin material at the virtual node P n -1 calculated based on the temperature difference between the virtual node P n (target virtual node) and the virtual node P n-1 (adjacent virtual node), but the present invention is not limited to this. For example, the correction coefficient k may be calculated based on the temperature difference between the virtual node P n and the virtual node P n-1 without calculating the elastic modulus. Also, the correction coefficient k may not be calculated for each of the multiple virtual nodes P, but may be a common value for one infinitesimal element. Also, the correction coefficient k may be a constant characteristic value that is preset for each characteristic of the resin material in the entire resin molded product 1.

また、上記実施形態では、解析配向Q(対象解析配向)と、補正係数kによって重み付けされた補正配向Rn-1(隣接補正配向)との平均を計算することによって、解析配向Qを補正する計算を行う例を示したが、本発明はこれに限られない。たとえば、補正係数kを用いずに解析配向Qと補正配向Rn-1との平均を計算することによって、解析配向Qを補正するようにしてもよい。すなわち、仮想節点P(対象仮想節点)の解析配向Qを補正する場合には、仮想節点P自身よりも表層側の仮想節点Pにおける弾性率、または、温度などの物性値に基づいて表層側から順に補正の計算を行うように構成されていれば、具体的な計算式はどのような形態であってもよい。 In the above embodiment, the analytical orientation Q n is corrected by calculating the average of the analytical orientation Q n (target analytical orientation) and the correction orientation R n-1 (adjacent correction orientation) weighted by the correction coefficient k, but the present invention is not limited to this. For example, the analytical orientation Q n may be corrected by calculating the average of the analytical orientation Q n and the correction orientation R n -1 without using the correction coefficient k. That is, when correcting the analytical orientation Q n of the virtual node P n (target virtual node), the specific calculation formula may be in any form as long as it is configured to perform correction calculation from the surface side based on the elastic modulus or physical property value such as temperature at the virtual node P on the surface side of the virtual node P n itself.

また、上記実施形態では、1つの微小要素において、20個の仮想節点Pを設定するとともに、20個の仮想節点Pの各々において、算出された補正配向Rの平均に基づいて微小要素の収縮量を解析する例を示したが、本発明はこれに限られない。たとえば、20個より少ない10個の仮想節点Pの各々において算出された補正配向Rの平均に基づいて微小要素の収縮量を解析してもよい。また、仮想節点Pを20個より多く設定してもよい。また、複数の微小要素の各々における厚み方向の大きさに応じて、複数の微小要素の各々における仮想節点Pの個数を異ならせるようにしてもよい。 In addition, in the above embodiment, an example was shown in which 20 virtual nodes P were set in one microelement, and the shrinkage amount of the microelement was analyzed based on the average of the calculated correction orientation R for each of the 20 virtual nodes P, but the present invention is not limited to this. For example, the shrinkage amount of the microelement may be analyzed based on the average of the calculated correction orientation R for each of 10 virtual nodes P, which is less than 20. Also, more than 20 virtual nodes P may be set. Also, the number of virtual nodes P in each of the multiple microelements may be made different depending on the size of each of the multiple microelements in the thickness direction.

また、上記実施形態では、複数の微小要素の各々において偶数個の仮想節点Pを設定する例を示したが、本発明はこれに限られない。たとえば、奇数個の仮想節点Pを設定するようにしてもよい。その場合に、表層側の両側から順に解析配向Qを補正する計算を行った場合には、中心の1つの仮想節点Pが余ることとなる。余った中心の1つの仮想節点Pにおいて解析配向Qを補正する計算を行う場合には、いずれか一方側に隣接する仮想節点Pにおける補正配向Rに基づいて解析配向Qを補正する計算を行うようにしてもよいし、両側に隣接する2つの仮想節点Pにおける補正配向Rに基づいて解析配向Qを補正する計算を行うようにしてもよい。 In the above embodiment, an example is shown in which an even number of virtual nodes P are set in each of a plurality of infinitesimal elements, but the present invention is not limited to this. For example, an odd number of virtual nodes P may be set. In that case, when a calculation is performed to correct the analytical orientation Q from both sides of the surface layer side, one virtual node P at the center will be left over. When a calculation is performed to correct the analytical orientation Q at the remaining one virtual node P at the center, a calculation may be performed to correct the analytical orientation Q based on the correction orientation R at the virtual node P adjacent to one side, or a calculation may be performed to correct the analytical orientation Q based on the correction orientation R at the two virtual nodes P adjacent to both sides.

また、上記実施形態では、複数の微小要素の各々において、算出された補正配向Rに基づいて、樹脂材料の流動方向(MD)と直交方向(TD)との2つの方向における収縮量を解析する例を示したが、本発明はこれに限られない。たとえば、厚み方向を加えた3つの方向における繊維配向を解析するとともに、3つの方向における配向を補正することによって、3つの方向における収縮量を解析するように構成してもよい。 In the above embodiment, an example was shown in which the amount of shrinkage in two directions, the flow direction (MD) and the perpendicular direction (TD) of the resin material, was analyzed for each of a plurality of microelements based on the calculated corrected orientation R, but the present invention is not limited to this. For example, the fiber orientation in three directions, including the thickness direction, may be analyzed, and the amount of shrinkage in the three directions may be analyzed by correcting the orientation in the three directions.

また、上記実施形態では、樹脂成形品1の解析モデルとして、樹脂成形品1を六面体(直方体)である複数の微小要素に分割する例を示したが、本発明はこれに限られない。たとえば、微小要素は、四面体であってもよい。また、微小要素は、三角形または四角形などの二次元の要素であってもよい。 In the above embodiment, an example was shown in which the resin molded product 1 was divided into multiple infinitesimal elements that are hexahedrons (rectangular parallelepipeds) as an analysis model of the resin molded product 1, but the present invention is not limited to this. For example, the infinitesimal elements may be tetrahedrons. Furthermore, the infinitesimal elements may be two-dimensional elements such as triangles or rectangles.

また、上記実施形態では、樹脂成形品1の各部分の厚み方向において分割される微小要素の個数が1つ(1層)となるように、樹脂成形品1を複数の微小要素に分割する例を示したが、本発明はこれに限られない。たとえば、樹脂成形品1の厚み方向において、分割される微小要素の個数が複数個となるようにしてもよい。 In addition, in the above embodiment, an example was shown in which the resin molded product 1 was divided into multiple microelements so that the number of microelements divided in the thickness direction of each part of the resin molded product 1 was one (one layer), but the present invention is not limited to this. For example, the number of microelements divided in the thickness direction of the resin molded product 1 may be multiple.

また、上記実施形態では、解析配向Q(対象解析配向)の補正の計算を行うために、対象の仮想節点P(対象仮想節点)と、隣接する仮想節点Pn-1(隣接仮想節点)との温度差を算出する場合に、金型の保圧が完了した時点における樹脂温度をシミュレートすることによって補正係数kを算出する例を示したが、本発明はこれに限られない。たとえば、金型から樹脂成形品1を取り出す時点における樹脂温度をシミュレートすることによって、補正係数kを算出するようにしてもよい。 In the above embodiment, when calculating the temperature difference between the target virtual node P n (target virtual node) and the adjacent virtual node P n-1 (adjacent virtual node) in order to perform a calculation for correcting the analytical orientation Q n (target analytical orientation), an example was shown in which the correction coefficient k is calculated by simulating the resin temperature at the time when the pressure retention of the mold is completed, but the present invention is not limited to this. For example, the correction coefficient k may be calculated by simulating the resin temperature at the time when the resin molded product 1 is removed from the mold.

1 樹脂成形品
11 微小要素分割部
12 仮想節点設定部
13 解析配向取得部
14 補正配向算出部
15 収縮量解析部
21 プログラム
100 解析装置(樹脂成形品解析装置、コンピュータ)
101 記録媒体
REFERENCE SIGNS LIST 1 resin molded product 11 minute element division unit 12 virtual node setting unit 13 analysis orientation acquisition unit 14 correction orientation calculation unit 15 shrinkage amount analysis unit 21 program 100 analysis device (resin molded product analysis device, computer)
101 Recording medium

Claims (11)

繊維を含む樹脂材料によって成形される樹脂成形品の収縮量を解析する樹脂成形品解析方法であって、
前記樹脂成形品を複数の微小要素に分割するステップと、
前記複数の微小要素の各々、もしくは、各前記微小要素の頂点を構成する複数の節点の各々において、前記樹脂成形品の厚み方向に沿って複数の仮想節点を設定するステップと、
各前記微小要素もしくは各前記節点の前記複数の仮想節点の各々において、前記繊維の配向を解析することによって解析配向を取得するステップと、
各前記微小要素もしくは各前記節点の前記複数の仮想節点の各々において、前記樹脂成形品の表層側から順に前記解析配向を補正する計算を行うことによって、前記収縮量を解析するための補正配向を算出するステップと、
各前記微小要素もしくは各前記節点の前記複数の仮想節点の各々において算出された前記補正配向に基づいて、前記複数の微小要素もしくは前記複数の節点の各々の前記収縮量を解析するとともに、解析された前記複数の微小要素もしくは前記複数の節点の各々における前記収縮量に基づいて、前記樹脂成形品の前記収縮量を解析するステップと、を備え、
前記補正配向を算出するステップは、各前記微小要素もしくは各前記節点の前記複数の仮想節点の各々において、自身よりも表層側の前記仮想節点における補正結果に基づいて前記解析配向を補正する計算を行うことによって、前記補正配向を算出するステップを含む、樹脂成形品解析方法。
A method for analyzing a resin molded product, which analyzes a shrinkage amount of a resin molded product molded from a resin material containing fiber, comprising the steps of:
Dividing the resin molded product into a plurality of minute elements;
setting a plurality of virtual nodes along a thickness direction of the resin molded product for each of the plurality of microelements or for each of a plurality of nodes constituting vertices of each of the plurality of microelements;
obtaining an analytical orientation by analyzing the orientation of the fiber at each of the plurality of virtual nodes of each of the microelements or each of the nodes;
A step of calculating a correction orientation for analyzing the amount of shrinkage by performing a calculation to correct the analytical orientation in each of the infinitesimal elements or each of the multiple virtual nodes of each of the nodes, starting from the surface layer side of the resin molded product;
analyzing the amount of shrinkage of each of the microelements or each of the nodes based on the corrected orientation calculated for each of the multiple virtual nodes of each of the microelements or each of the nodes, and analyzing the amount of shrinkage of the resin molded product based on the amount of shrinkage of each of the analyzed multiple microelements or each of the multiple nodes,
The step of calculating the corrected orientation includes a step of calculating the corrected orientation by performing a calculation to correct the analysis orientation at each of the multiple virtual nodes of each of the infinitesimal elements or each of the nodes based on a correction result at the virtual node on the surface side of the infinitesimal element or each of the nodes.
前記解析配向を取得するステップは、前記複数の仮想節点のうちの一の前記仮想節点である対象仮想節点において、前記対象仮想節点における前記解析配向である対象解析配向を取得するステップを含み、
前記補正配向を算出するステップは、
前記複数の仮想節点のうちの前記対象仮想節点の表層側に隣接する隣接仮想節点において算出された前記補正配向である隣接補正配向に基づいて、前記対象解析配向を補正する計算を行うことによって、前記対象仮想節点における前記補正配向である対象補正配向を算出するステップと、
前記対象補正配向を算出するステップを、前記複数の仮想節点の各々において、前記樹脂成形品の表層側から順に繰り返すステップとを含む、請求項1に記載の樹脂成形品解析方法。
The step of acquiring the analytical orientation includes a step of acquiring, at a target virtual node which is one of the virtual nodes among the plurality of virtual nodes, a target analytical orientation which is the analytical orientation at the target virtual node;
The step of calculating the corrected orientation includes:
A step of calculating a target correction orientation, which is the correction orientation at the target virtual node, by performing a calculation to correct the target analysis orientation based on an adjacent correction orientation, which is the correction orientation calculated at an adjacent virtual node adjacent to the surface side of the target virtual node among the plurality of virtual nodes;
2. The method for analyzing a resin molded product according to claim 1, further comprising: repeating the step of calculating the target corrected orientation for each of the plurality of virtual nodes in order from a surface layer side of the resin molded product.
前記対象補正配向を算出するステップは、前記隣接補正配向と前記隣接仮想節点における前記樹脂材料の弾性率とに基づいて、前記対象解析配向を補正する計算を行うことによって、前記対象補正配向を算出するステップを含む、請求項2に記載の樹脂成形品解析方法。 The resin molded product analysis method according to claim 2, wherein the step of calculating the target correction orientation includes a step of calculating the target correction orientation by performing a calculation to correct the target analysis orientation based on the adjacent correction orientation and the elastic modulus of the resin material at the adjacent virtual node. 前記対象補正配向を算出するステップは、
前記隣接仮想節点における前記樹脂材料の弾性率と、前記対象仮想節点と前記隣接仮想節点との間の温度差とに基づいて、補正係数を算出するステップと、
前記隣接補正配向と前記補正係数とに基づいて、前記対象解析配向を補正する計算を行うステップとを含む、請求項3に記載の樹脂成形品解析方法。
The step of calculating the object corrected orientation comprises:
calculating a correction coefficient based on an elastic modulus of the resin material at the adjacent virtual node and a temperature difference between the target virtual node and the adjacent virtual node;
The resin molded product analysis method according to claim 3 , further comprising: performing a calculation to correct the target analysis orientation based on the adjacent correction orientation and the correction coefficient.
前記対象解析配向を補正する計算を行うステップは、前記対象解析配向と、前記補正係数によって重み付けされた前記隣接補正配向との平均を計算することによって、前記対象解析配向を補正する計算を行うステップを含む、請求項4に記載の樹脂成形品解析方法。 The resin molded product analysis method according to claim 4, wherein the step of performing a calculation to correct the target analysis orientation includes a step of performing a calculation to correct the target analysis orientation by calculating an average of the target analysis orientation and the adjacent correction orientation weighted by the correction coefficient. 前記収縮量を解析するステップは、各前記微小要素もしくは各前記節点の前記複数の仮想節点の各々において算出された前記補正配向の平均に基づいて、前記複数の微小要素もしくは前記複数の節点の各々の前記収縮量を解析するとともに、解析された前記複数の微小要素もしくは前記複数の節点の各々における前記収縮量に基づいて、前記樹脂成形品の前記収縮量を解析するステップを含む、請求項1~5のいずれか1項に記載の樹脂成形品解析方法。 The method for analyzing a resin molded product according to any one of claims 1 to 5, wherein the step of analyzing the amount of shrinkage includes a step of analyzing the amount of shrinkage of each of the microelements or each of the nodes based on the average of the corrected orientation calculated for each of the multiple virtual nodes of each of the microelements or each of the nodes, and a step of analyzing the amount of shrinkage of the resin molded product based on the amount of shrinkage of each of the analyzed multiple microelements or each of the multiple nodes. 前記複数の仮想節点を設定するステップは、前記複数の微小要素もしくは前記複数の節点の各々において、前記樹脂成形品の厚み方向に沿って偶数個の前記複数の仮想節点を設定するステップを含む、請求項1~6のいずれか1項に記載の樹脂成形品解析方法。 The method for analyzing a resin molded product according to any one of claims 1 to 6, wherein the step of setting the plurality of virtual nodes includes a step of setting an even number of the plurality of virtual nodes along the thickness direction of the resin molded product in each of the plurality of infinitesimal elements or the plurality of nodes. 前記収縮量を解析するステップは、各前記微小要素もしくは各前記節点の前記複数の仮想節点の各々において算出された前記補正配向に基づいて、前記複数の微小要素もしくは前記複数の節点の各々において、前記樹脂材料の流動方向と、前記樹脂成形品の表面に沿う面内において前記流動方向と直交する直交方向との前記収縮量を解析するとともに、解析された前記複数の微小要素もしくは前記複数の節点の各々における前記流動方向と前記直交方向との前記収縮量に基づいて、前記樹脂成形品の前記流動方向と前記直交方向との前記収縮量を解析するステップを含む、請求項1~7のいずれか1項に記載の樹脂成形品解析方法。 The resin molded product analysis method according to any one of claims 1 to 7, wherein the step of analyzing the shrinkage amount includes a step of analyzing the shrinkage amount between the flow direction of the resin material and an orthogonal direction orthogonal to the flow direction in a plane along the surface of the resin molded product at each of the multiple microelements or multiple nodes based on the corrected orientation calculated at each of the multiple virtual nodes of each of the multiple microelements or each of the nodes, and analyzing the shrinkage amount between the flow direction and the orthogonal direction of the resin molded product based on the shrinkage amount between the flow direction and the orthogonal direction at each of the analyzed multiple microelements or multiple nodes. 請求項1~8のいずれか1項に記載された樹脂成形品解析方法をコンピュータに実行させる、プログラム。 A program that causes a computer to execute the resin molded product analysis method described in any one of claims 1 to 8. 請求項9に記載のプログラムが記録され、前記コンピュータにより読み取り可能な、記録媒体。 A recording medium on which the program according to claim 9 is recorded and which is readable by the computer. 繊維を含む樹脂材料によって成形される樹脂成形品の収縮量を解析する樹脂成形品解析装置であって、
前記樹脂成形品を複数の微小要素に分割する微小要素分割部と、
前記複数の微小要素の各々、もしくは、各前記微小要素の頂点を構成する複数の節点の各々において、前記樹脂成形品の厚み方向に沿って複数の仮想節点を設定する仮想節点設定部と、
各前記微小要素もしくは各前記節点の前記複数の仮想節点の各々において、前記繊維の配向を解析することによって解析配向を取得する解析配向取得部と、
各前記微小要素もしくは各前記節点の前記複数の仮想節点の各々において、前記樹脂成形品の表層側から順に前記解析配向を補正する計算を行うことによって、前記収縮量を解析するための補正配向を算出する補正配向算出部と、
各前記微小要素もしくは各前記節点の前記複数の仮想節点の各々において算出された前記補正配向に基づいて、前記複数の微小要素もしくは前記複数の節点の各々の前記収縮量を解析するとともに、解析された前記複数の微小要素もしくは前記複数の節点の各々における前記収縮量に基づいて、前記樹脂成形品の前記収縮量を解析する収縮量解析部と、を備え、
前記補正配向算出部は、各前記微小要素もしくは各前記節点の前記複数の仮想節点の各々において、自身よりも表層側の前記仮想節点における補正結果に基づいて前記解析配向を補正する計算を行うことによって、前記補正配向を算出するように構成されている、樹脂成形品解析装置。
A resin molded product analysis device that analyzes a shrinkage amount of a resin molded product molded from a resin material containing fiber, comprising:
a microelement dividing unit that divides the resin molded product into a plurality of microelements;
a virtual node setting unit that sets a plurality of virtual nodes along a thickness direction of the resin molded product in each of the plurality of infinitesimal elements or in each of a plurality of nodes that constitute the vertices of each of the infinitesimal elements;
an analytical orientation acquisition unit that acquires an analytical orientation by analyzing the orientation of the fiber in each of the infinitesimal elements or each of the plurality of virtual nodes of each of the nodes;
a correction orientation calculation unit that calculates a correction orientation for analyzing the shrinkage amount by performing a calculation to correct the analytical orientation in each of the infinitesimal elements or each of the multiple virtual nodes of each of the nodes, starting from the surface layer side of the resin molded product;
a shrinkage amount analysis unit that analyzes the amount of shrinkage of each of the plurality of infinitesimal elements or each of the plurality of nodes based on the corrected orientation calculated for each of the plurality of virtual nodes of each of the infinitesimal elements or each of the nodes, and analyzes the amount of shrinkage of the resin molded product based on the amount of shrinkage of each of the analyzed plurality of infinitesimal elements or each of the plurality of nodes,
The corrected orientation calculation unit is configured to calculate the corrected orientation at each of the multiple virtual nodes of each of the infinitesimal elements or each of the nodes by performing calculations to correct the analysis orientation based on the correction result at the virtual node that is on the surface side of the infinitesimal element or each of the nodes, in the plastic molded product analysis device.
JP2022006020A 2022-01-18 2022-01-18 Method, program, recording medium, and apparatus for analyzing resin molded products Active JP7645207B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2022006020A JP7645207B2 (en) 2022-01-18 2022-01-18 Method, program, recording medium, and apparatus for analyzing resin molded products

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2022006020A JP7645207B2 (en) 2022-01-18 2022-01-18 Method, program, recording medium, and apparatus for analyzing resin molded products

Publications (2)

Publication Number Publication Date
JP2023104808A JP2023104808A (en) 2023-07-28
JP7645207B2 true JP7645207B2 (en) 2025-03-13

Family

ID=87379106

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022006020A Active JP7645207B2 (en) 2022-01-18 2022-01-18 Method, program, recording medium, and apparatus for analyzing resin molded products

Country Status (1)

Country Link
JP (1) JP7645207B2 (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002321265A (en) 2001-04-25 2002-11-05 Toray Ind Inc Method for forecasting weldline and system therefor
JP2008120089A (en) 2007-12-11 2008-05-29 Fujitsu Ltd Data converter, flow analyzer, structure analyzer, data conversion program, flow analysis program, structure analysis program, and data conversion method
JP2011201108A (en) 2010-03-25 2011-10-13 Mazda Motor Corp Flow analyzing method, flow analyzer and flow analysis program of resin molding
JP2013082096A (en) 2011-10-06 2013-05-09 Ube Industries Ltd Method for estimating natural frequency of fiber reinforced resin injection molded article
JP6203787B2 (en) 2015-07-03 2017-09-27 本田技研工業株式会社 Computer-aided resin behavior analyzer
WO2018116898A1 (en) 2016-12-19 2018-06-28 東レエンジニアリング株式会社 Resin flow analysis method, resin flow analysis device, program, and computer-readable recording medium
JP2019181801A (en) 2018-04-10 2019-10-24 株式会社豊田中央研究所 Shrinkage rate prediction device, shrinkage rate prediction model learning device, and program
JP2021008081A (en) 2019-07-02 2021-01-28 本田技研工業株式会社 Resin behavior analyzer, resin behavior analyzing method, and resin behavior analyzing program

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002321265A (en) 2001-04-25 2002-11-05 Toray Ind Inc Method for forecasting weldline and system therefor
JP2008120089A (en) 2007-12-11 2008-05-29 Fujitsu Ltd Data converter, flow analyzer, structure analyzer, data conversion program, flow analysis program, structure analysis program, and data conversion method
JP2011201108A (en) 2010-03-25 2011-10-13 Mazda Motor Corp Flow analyzing method, flow analyzer and flow analysis program of resin molding
JP2013082096A (en) 2011-10-06 2013-05-09 Ube Industries Ltd Method for estimating natural frequency of fiber reinforced resin injection molded article
JP6203787B2 (en) 2015-07-03 2017-09-27 本田技研工業株式会社 Computer-aided resin behavior analyzer
WO2018116898A1 (en) 2016-12-19 2018-06-28 東レエンジニアリング株式会社 Resin flow analysis method, resin flow analysis device, program, and computer-readable recording medium
JP2019181801A (en) 2018-04-10 2019-10-24 株式会社豊田中央研究所 Shrinkage rate prediction device, shrinkage rate prediction model learning device, and program
JP2021008081A (en) 2019-07-02 2021-01-28 本田技研工業株式会社 Resin behavior analyzer, resin behavior analyzing method, and resin behavior analyzing program

Also Published As

Publication number Publication date
JP2023104808A (en) 2023-07-28

Similar Documents

Publication Publication Date Title
JP7118941B2 (en) RESIN MOLDING ANALYSIS METHOD, PROGRAM AND RECORDING MEDIUM
US5811133A (en) Injection molding apparatus for predicting deformation amount in injection-molded article
JP3848602B2 (en) Resin molded product design support apparatus and method
CN101685475A (en) Analytical model preparation method, and simulation system method for predicting molding failure
JP2014226871A (en) Method for simulating behavior of filler, and method for analyzing physical property of composite material
JP2006205740A (en) Resin molded product design support apparatus and method
JP7645207B2 (en) Method, program, recording medium, and apparatus for analyzing resin molded products
Tie et al. Three-dimensional finite element method for the filling simulation of injection molding
JP6684166B2 (en) Resin flow analysis method, program, and computer-readable recording medium
US11651124B2 (en) Anti-warping design method for resin molded article, recording medium, and anti-warping design device for resin molded article
JP2001191336A (en) Mold design apparatus and method for designing mold shape
JP6618069B2 (en) Composite molded product design support device, composite molded product manufacturing method, computer software, storage medium
CN115938507A (en) Injection Molding Warpage Prediction
JP2003271678A (en) Numerical analysis method and device
US11602876B2 (en) Machine learning device and design support device
JP6639899B2 (en) Molded article design support method, molded article design support apparatus, computer software, storage medium
JPH08230008A (en) Method and device for estimating warping deformation of injection molded article
JP2002219739A (en) Method for estimating amount of change in shape of injection-molded article
JP4330404B2 (en) Molded product design support apparatus, design support method, and software
JP2007083602A (en) Method for forecasting molding shrinkage rate of injection-molded article
JP6900184B2 (en) Resin flow analysis method, resin flow analysis device, program and computer readable recording medium
JP3582930B2 (en) Manufacturing method for injection molded products
US20180267011A1 (en) Method of predicting deformation of resin molded article
JP2002082996A (en) Method and apparatus for analyzing structure
JP6689015B2 (en) Fiber concentration setting method

Legal Events

Date Code Title Description
RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7426

Effective date: 20220203

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20220124

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20220207

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240502

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20250128

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20250303

R150 Certificate of patent or registration of utility model

Ref document number: 7645207

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150