JP2017013437A - Computer-aided resin behavior analysis device - Google Patents

Computer-aided resin behavior analysis device Download PDF

Info

Publication number
JP2017013437A
JP2017013437A JP2015134616A JP2015134616A JP2017013437A JP 2017013437 A JP2017013437 A JP 2017013437A JP 2015134616 A JP2015134616 A JP 2015134616A JP 2015134616 A JP2015134616 A JP 2015134616A JP 2017013437 A JP2017013437 A JP 2017013437A
Authority
JP
Japan
Prior art keywords
fibers
calculated
bending
fiber
computer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2015134616A
Other languages
Japanese (ja)
Other versions
JP6203787B2 (en
Inventor
功司 団
Koji Dan
功司 団
正俊 小林
Masatoshi Kobayashi
正俊 小林
剛志 馬場
Tsuyoshi Baba
剛志 馬場
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.)
Honda Motor Co Ltd
Original Assignee
Honda Motor 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 Honda Motor Co Ltd filed Critical Honda Motor Co Ltd
Priority to JP2015134616A priority Critical patent/JP6203787B2/en
Priority to US15/096,168 priority patent/US10223481B2/en
Publication of JP2017013437A publication Critical patent/JP2017013437A/en
Application granted granted Critical
Publication of JP6203787B2 publication Critical patent/JP6203787B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Casting Or Compression Moulding Of Plastics Or The Like (AREA)
  • Moulding By Coating Moulds (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a computer-aided resin behavior analysis device enabling prediction of deterioration in a structural strength of a molded product due to fiber bending when analyzing a behavior when a fiber-reinforced resin is molded in a mold on a predetermined molding condition through a simulation program stored in a computer.SOLUTION: An analysis device (S10-S42) which analyzes a behavior when a resin to which a plurality of continuous fibers and discontinuous long fibers are mixed on a predetermined molding condition through a simulation program stored in a computer includes making the simulation program select at least a part of fiber group out of the plurality of fibers when an analysis condition including a joint F of the fiber is input, calculating a bending ratio Af of the respective fibers constituting the fiber group on the basis of the analysis condition, evaluating the bending of the respective fibers on the basis of the bending ratio Af, and evaluating the strength of the fiber group selected on the basis of the evaluation result (S44).SELECTED DRAWING: Figure 2

Description

この発明はコンピュータ支援の樹脂挙動解析装置に関し、より詳しくは連続繊維および不連続な長繊維が混入された樹脂を金型内で成形するときの挙動をコンピュータに格納されたシミュレーションプログラムを介して解析するコンピュータ支援(Computer Aided Engineering)の樹脂挙動解析装置に関する。   The present invention relates to a computer-aided resin behavior analysis apparatus, and more specifically, analyzes a behavior when molding a resin mixed with continuous fibers and discontinuous long fibers in a mold through a simulation program stored in a computer. The present invention relates to a resin behavior analysis apparatus of Computer Aided Engineering.

従来から、繊維強化樹脂のプレス成形時に発生する繊維曲がりに起因する成形品の構造強度低下を低減する技術が知られている。例えば、特許文献1記載の技術は、キャビティ内に樹脂を充填して繊維強化樹脂材を製造するとき、下型の内部に摺動自在に設けられた中子の上に連続繊維補強材を載置し、それを押上げて樹脂内に埋め込むようにしている。   2. Description of the Related Art Conventionally, a technique for reducing a decrease in structural strength of a molded product due to fiber bending that occurs during press molding of a fiber reinforced resin is known. For example, in the technique described in Patent Document 1, when a fiber reinforced resin material is manufactured by filling a resin in a cavity, a continuous fiber reinforcing material is placed on a core slidably provided in a lower mold. And push it up to embed it in the resin.

また、特許文献2記載の技術は、樹脂基材に繊維強化樹脂材を補強材として配置して樹脂成形体を製造するとき、下型において繊維強化樹脂に可動ピンを貫通させ、樹脂基材を貫通孔に導入しながら成形するようにしている。   In addition, the technique described in Patent Document 2 is such that when a fiber reinforced resin material is disposed on a resin base material as a reinforcing material to manufacture a resin molded body, a movable pin is passed through the fiber reinforced resin in the lower mold, Molding is performed while being introduced into the through hole.

特許第5418684号公報Japanese Patent No. 5418684 特開2013−022852号公報JP2013-022852A

ところで、特許文献1,2記載の技術は研究開発にあたって試作品(成形品)の構造強度を確認するが、実際の成形品が必要となるため大きな手間がかかる。これを解決するには樹脂挙動解析を用いて成形品の構造強度または弾性率(剛性)を予測することが考えられるが、繊維曲がりに起因する成形品の構造強度低下または弾性率(剛性)低下を対象とした解析技術は存在していない。   By the way, although the technique of patent documents 1 and 2 confirms the structural strength of a prototype (molded product) in research and development, since an actual molded product is required, it takes a lot of work. To solve this, it is conceivable to predict the structural strength or elastic modulus (rigidity) of the molded product using resin behavior analysis. However, the structural strength or elastic modulus (rigidity) of the molded product is reduced due to fiber bending. There is no analysis technology for.

従って、この発明の目的は上記した課題を解決し、連続繊維および不連続な長繊維が複数本混入された繊維強化樹脂を金型内で所定の成形条件で成形するときの挙動をコンピュータに格納されたシミュレーションプログラムを介して解析するとき、繊維曲がりに起因する成形品の構造強度低下または弾性率(剛性)低下の予測を可能としたコンピュータ支援の樹脂挙動解析装置を提供することにある。   Therefore, the object of the present invention is to solve the above-mentioned problems and store the behavior when a fiber reinforced resin mixed with a plurality of continuous fibers and discontinuous long fibers is molded in a mold under predetermined molding conditions in a computer. It is an object of the present invention to provide a computer-aided resin behavior analysis apparatus capable of predicting a decrease in structural strength or a decrease in elastic modulus (rigidity) of a molded product due to fiber bending when an analysis is performed via a simulation program.

上記した課題を解決するために、請求項1にあっては、連続繊維および不連続な長繊維が複数本混入された樹脂を金型内で所定の成形条件で成形するときの前記繊維の曲がり挙動をコンピュータに格納されたシミュレーションプログラムを介して解析する樹脂挙動解析装置において、前記シミュレーションプログラムが、前記複数本の繊維の複数個の節点Fを少なくとも含む解析条件が入力されるとき、前記複数本の繊維のうちの少なくとも一部の繊維群を選択し、前記入力された解析条件に基づいて前記選択された繊維群を構成するそれぞれの繊維の曲がり率Afを算出し、前記算出されたそれぞれの繊維の曲がり率Afに基づいて前記それぞれの繊維の曲がりを評価し、前記得られたそれぞれの繊維の曲がり評価結果に基づいて前記選択された繊維群の強度または弾性率(剛性)を評価するステップからなる如く構成した。   In order to solve the above-described problem, in Claim 1, the bending of the fiber when a resin in which a plurality of continuous fibers and discontinuous long fibers are mixed is molded in a mold under predetermined molding conditions. In the resin behavior analysis apparatus for analyzing behavior via a simulation program stored in a computer, when the simulation program receives an analysis condition including at least a plurality of nodes F of the plurality of fibers, the plurality of And selecting a fiber group of at least a part of the fibers, calculating a bending rate Af of each of the fibers constituting the selected fiber group based on the input analysis condition, The bending of each of the fibers is evaluated based on the bending rate Af of the fiber, and the selection is performed based on the bending evaluation result of each of the obtained fibers. It was constructed as consisting of the step of evaluating the strength or elastic modulus of the fiber group a (rigid).

請求項2に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記シミュレーションプログラムは、前記得られた繊維群の強度または弾性率(剛性)評価結果に基づいて前記複数本の繊維の全体の強度または弾性率(剛性)を評価するステップを含む如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 2, the simulation program calculates the overall strength of the plurality of fibers based on the strength or elastic modulus (rigidity) evaluation result of the obtained fiber group. Or it comprised so that the step which evaluates an elastic modulus (rigidity) might be included.

請求項3に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がりを評価するステップは、前記算出されたそれぞれの繊維の曲がり率Afを所定のしきい値と比較して所定の範囲ごとに分類して評価するステップからなり、前記選択された繊維群の強度または弾性率(剛性)を評価するステップは、前記所定の範囲に属する繊維の割合を算出し、前記算出された繊維の割合に基づいて前記選択された繊維群の強度または弾性率(剛性)評価値を算出するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 3, the step of evaluating the bending of each of the fibers compares the calculated bending rate Af of each of the fibers with a predetermined threshold value. The method comprises the steps of classifying and evaluating for each predetermined range, and the step of evaluating the strength or elastic modulus (rigidity) of the selected fiber group calculates the ratio of fibers belonging to the predetermined range and calculates And a step of calculating a strength or elastic modulus (stiffness) evaluation value of the selected fiber group based on the ratio of the fibers.

請求項4に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記選択された繊維群の強度または弾性率(剛性)を評価するステップは、前記所定の範囲に属する繊維の割合を前記所定の範囲ごとに表示して評価するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 4, the step of evaluating the strength or elastic modulus (rigidity) of the selected fiber group may include the ratio of the fibers belonging to the predetermined range. It comprised so that it might consist of the step displayed and evaluated for every range.

請求項5に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記選択された繊維群の強度または弾性率(剛性)を評価するステップは、前記算出された繊維群の強度または弾性率(剛性)評価値から予め設定された特性に基づいて前記選択された繊維群の強度または弾性率(剛性)を算出して評価するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 5, the step of evaluating the strength or elastic modulus (rigidity) of the selected fiber group includes the calculated strength or elastic modulus (rigidity) of the fiber group. ) It is configured to include a step of calculating and evaluating the strength or elastic modulus (rigidity) of the selected fiber group based on a preset characteristic from the evaluation value.

請求項6に係るコンピュータ支援の樹脂挙動解析装置にあっては、連続繊維および不連続な長繊維が複数本混入された樹脂を金型内で所定の成形条件で成形するときの前記繊維の曲がり挙動をコンピュータに格納されたシミュレーションプログラムを介して解析する樹脂挙動解析装置において、前記シミュレーションプログラムが、前記複数本の繊維の複数個の節点Fを少なくとも含む解析条件が入力されるとき、前記複数本の繊維のうちの少なくとも一部の繊維群を選択し、前記入力された解析条件に基づいて前記選択された繊維群を構成するそれぞれの繊維の曲がり率Afを算出し、前記算出されたそれぞれの繊維の曲がり率Afを表示色に関連付けて前記それぞれの繊維の曲がりを評価し、前記得られたそれぞれの繊維の曲がり評価結果に基づいて前記選択された繊維群の強度または弾性率(剛性)を評価するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 6, the fiber bends when a resin in which a plurality of continuous fibers and discontinuous long fibers are mixed is molded in a mold under predetermined molding conditions. In the resin behavior analysis apparatus for analyzing behavior via a simulation program stored in a computer, when the simulation program receives an analysis condition including at least a plurality of nodes F of the plurality of fibers, the plurality of And selecting a fiber group of at least a part of the fibers, calculating a bending rate Af of each of the fibers constituting the selected fiber group based on the input analysis condition, The bending rate Af of the fiber is related to the display color to evaluate the bending of each of the fibers, and the bending evaluation result of each of the obtained fibers is evaluated. It was constructed as consisting of the step of evaluating the strength or elastic modulus (stiffness) the selected fiber group based on.

請求項7に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記シミュレーションプログラムは、前記得られた繊維群の強度または弾性率(剛性)評価結果に基づいて前記複数本の繊維の全体の強度または弾性率(剛性)を評価するステップを含む如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 7, the simulation program calculates the overall strength of the plurality of fibers based on the strength or elastic modulus (rigidity) evaluation result of the obtained fiber group. Or it comprised so that the step which evaluates an elastic modulus (rigidity) might be included.

請求項8に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がり率Afを算出するステップは、前記複数個の節点Fの中の少なくとも一つの節点Fnから得られる評価長さに対する前記成形条件から予測される前記それぞれの繊維の曲がり率Afを算出するステップからなり、前記それぞれの繊維の曲がりを評価するステップは、前記算出された評価長さに対する曲がり率Afに基づいて前記それぞれの繊維の曲がりを評価するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 8, the step of calculating the bending rate Af of each fiber is an evaluation length obtained from at least one node Fn among the plurality of nodes F. The step of calculating the bending rate Af of each of the fibers predicted from the molding conditions with respect to the thickness, and the step of evaluating the bending of each of the fibers is based on the bending rate Af with respect to the calculated evaluation length. It comprised so that it might consist of the step which evaluates the bending of each said fiber.

請求項9に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がり率Afを算出するステップは、前記複数個の節点Fの少なくとも一つの節点Fnを曲線または直線で接続して生成される三次元図形に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 9, the step of calculating the bending rate Af of each of the fibers includes connecting at least one node Fn of the plurality of nodes F with a curve or a straight line. The method comprises a step of calculating the bending rate Af of each of the fibers based on the three-dimensional figure generated in this manner.

請求項10に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がり率Afを算出するステップは、前記複数個の節点Fの少なくとも一つの節点Fnを曲線または直線で接続して生成される三次元図形上の点を曲線または直線で接続して生成される第2の三次元図形に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 10, the step of calculating the bending rate Af of each of the fibers includes connecting at least one node Fn of the plurality of nodes F with a curve or a straight line. The method comprises a step of calculating the bending rate Af of each of the fibers based on a second three-dimensional graphic generated by connecting points on the three-dimensional graphic generated by a curve or a straight line.

請求項11に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記第2の三次元図形は三角形からなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 11, the second three-dimensional figure is configured to be a triangle.

請求項12に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がり率Afを算出するステップは、前記節点Fnとそれに隣接する節点Fn-1,Fn+1とをそれぞれ結ぶ線FnFn-1FnFn+1を斜辺とすると共に、隣接する節点Fn-1,Fn+1を結ぶ線Fn-1Fn+1を底辺とする前記三角形を設定し、前記節点Fnあるいはその付近の任意の点から前記底辺に向けられた垂線の長さh1を算出し、前記算出された垂線の長さh1を前記斜辺FnFn-1FnFn+1の長さの合計値で除算してアスペクト比As1を算出し、前記算出されたアスペクト比As1に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 12, the step of calculating the bending rate Af of each of the fibers connects the node Fn and the nodes Fn-1 and Fn + 1 adjacent thereto. The triangle having the line FnFn-1FnFn + 1 as the hypotenuse and the line Fn-1Fn + 1 connecting the adjacent nodes Fn-1 and Fn + 1 as the base is set, and the node Fn or any point in the vicinity thereof is set. To calculate the aspect ratio As1 by dividing the length h1 of the perpendicular directed to the base from the total length of the hypotenuses FnFn-1FnFn + 1, Based on the calculated aspect ratio As1, the bending rate Af of each of the fibers is calculated.

請求項13に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がり率Afを算出するステップは、前記節点Fn+1とそれに隣接する節点Fn,Fn+2とをそれぞれ結ぶ線Fn+1FnFn+1Fn+2を斜辺とすると共に、隣接する節点Fn,Fn+2を結ぶ線FnFn+2を底辺とする前記三角形を設定し、前記節点Fn+1あるいはその付近の任意の点から前記底辺に向けられた垂線の長さh2を算出し、前記算出された垂線の長さh2を前記斜辺Fn+1FnFn+1Fn+2の長さの合計値で除算してアスペクト比As2を算出し、前記算出されたアスペクト比As2と前記節点Fnについて算出されたアスペクト比As1とに基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 13, the step of calculating the bending rate Af of each fiber connects the node Fn + 1 and the nodes Fn and Fn + 2 adjacent thereto. The triangle having the line Fn + 1FnFn + 1Fn + 2 as the hypotenuse and the line FnFn + 2 connecting the adjacent nodes Fn and Fn + 2 as the base is set, and the node Fn + 1 or any point in the vicinity thereof is set. To calculate the aspect ratio As2 by dividing the length h2 of the perpendicular directed to the bottom side by the total length of the hypotenuses Fn + 1FnFn + 1Fn + 2 And it comprised so that it might consist of the step which calculates the bending rate Af of each said fiber based on the calculated aspect-ratio As2 and the aspect-ratio As1 calculated about the said node Fn.

請求項14に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がり率Afを算出するステップは、前記節点Fn+2とそれに隣接する/しない節点Fn+2,Fn+mとをそれぞれ結ぶ線Fn+2Fn+1Fn+2Fn+mを斜辺とすると共に、隣接する節点Fn+1,Fn+mを結ぶ線Fn+1Fn+mを底辺とする前記三角形を設定し、前記節点Fn+2あるいはその付近の任意の点から前記底辺に向けられた垂線の長さhmを算出し、前記算出された垂線の長さhmを前記斜辺Fn+2Fn+1Fn+2Fn+mの長さの合計値で除算してアスペクト比Asmを算出し、前記算出されたアスペクト比Asmと前記節点Fnと前記節点Fn+1について算出されたアスペクト比As1,As2とに基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 14, the step of calculating the bending rate Af of each of the fibers includes the node Fn + 2 and the nodes Fn + 2 and Fn + m adjacent to / not adjacent to the node Fn + 2. And the line Fn + 2Fn + 1Fn + 2Fn + m that connects to each other is set as the hypotenuse and the triangle with the line Fn + 1Fn + m connecting adjacent nodes Fn + 1 and Fn + m as the base is set. The length hm of the perpendicular line directed to the base from Fn + 2 or an arbitrary point in the vicinity thereof is calculated, and the calculated length hm of the perpendicular line is the length of the hypotenuse Fn + 2Fn + 1Fn + 2Fn + m The aspect ratio Asm is calculated by dividing by the total value of the above, and the bending of each of the fibers based on the calculated aspect ratio Asm, the aspect ratios As1 and As2 calculated for the node Fn and the node Fn + 1 It consists of steps to calculate rate Af .

請求項15に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がり率Afを算出するステップは、前記節点Fnとそれに隣接する節点Fn-1,Fn+1とをそれぞれ結ぶ線FnFn-1FnFn+1上のいずれかの任意の点Pと隣接する節点Fn-1,Fn+2を結ぶ線PFn-1PFn+2を斜辺とすると共に、節点Fn-1,Fn+2を結ぶ線Fn-1Fn+2を底辺とする前記三角形を設定し、前記任意の点Pから前記底辺に向けられた垂線の長さh1を算出し、前記算出された垂線の長さh1を前記斜辺PFn-1PFn+2の長さの合計値で除算してアスペクト比As1を算出し、次いで他の節点について前記アスペクト比Asnを算出し、前記算出されたアスペクト比As1,Asnに基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 15, the step of calculating the bending rate Af of each of the fibers connects the node Fn and the nodes Fn-1 and Fn + 1 adjacent thereto. The line PFn-1 PFn + 2 connecting any arbitrary point P on the line FnFn-1FnFn + 1 and the adjacent nodes Fn-1, Fn + 2 is set as the hypotenuse and the nodes Fn-1, Fn + 2 are connected. The triangle whose base is the line Fn-1Fn + 2 is set, the length h1 of the perpendicular directed from the arbitrary point P to the base is calculated, and the length h1 of the calculated perpendicular is set to the hypotenuse PFn The aspect ratio As1 is calculated by dividing by the total length of -1 PFn + 2, then the aspect ratio Asn is calculated for other nodes, and the respective fibers are calculated based on the calculated aspect ratios As1 and Asn. From the step of calculating the bending rate Af of It was constructed as.

請求項16に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がり率Afを算出するステップは、前記算出された値の平均値に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 16, the step of calculating the bending rate Af of each of the fibers includes the bending rate Af of each of the fibers based on an average value of the calculated values. It comprised so that it might consist of the step which calculates.

請求項17に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がり率Afを算出するステップは、前記算出された値の最大値に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 17, the step of calculating the bending rate Af of each of the fibers includes the bending rate Af of each of the fibers based on a maximum value of the calculated values. It comprised so that it might consist of the step which calculates.

請求項18に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記それぞれの繊維の曲がり率Afを算出するステップは、前記金型のキャビティにおける前記繊維の三次元空間での座標位置に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 18, the step of calculating the bending rate Af of each of the fibers is based on a coordinate position of the fibers in a three-dimensional space in the cavity of the mold. It comprised so that it might consist of the step which calculates the bending rate Af of each said fiber.

請求項19に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記評価長さは前記三次元図形の始点から終点までの長さからなる如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 19, the evaluation length is configured to be a length from a start point to an end point of the three-dimensional figure.

請求項20に係るコンピュータ支援の樹脂挙動解析装置にあっては、前記シミュレーションプログラムは、分析器を介して得られる前記複数本の繊維の複数個の節点Fを少なくとも含むデータが入力されるとき、前記複数本の繊維のうちの少なくとも一部の繊維群を選択し、前記入力されたデータに基づいて前記選択された繊維群を構成するそれぞれの繊維の曲がり率Afを算出し、前記算出されたそれぞれの繊維の曲がり率Afに基づいて前記それぞれの繊維の曲がりを評価し、前記得られたそれぞれの繊維の曲がり評価結果に基づいて前記選択された繊維群の強度または弾性率を評価するステップを含む如く構成した。   In the computer-aided resin behavior analysis apparatus according to claim 20, the simulation program receives data including at least a plurality of nodes F of the plurality of fibers obtained through an analyzer. Selecting at least a part of the plurality of fibers, calculating a bending rate Af of each fiber constituting the selected fiber group based on the input data, and calculating the calculated Evaluating the bending of each of the fibers based on the bending rate Af of each of the fibers, and evaluating the strength or elastic modulus of the selected fiber group based on the bending evaluation result of each of the obtained fibers. Constructed to include.

請求項1に係るコンピュータ支援の樹脂挙動解析装置にあっては、金型内で所定の成形条件で成形するときの繊維混入樹脂の挙動を解析するシミュレーションプログラムが、複数本の繊維の複数個の節点Fを少なくとも含む解析条件が入力されるとき、複数本の繊維のうちの少なくとも一部の繊維群を選択し、入力された解析条件に基づいて選択された繊維群を構成するそれぞれの繊維の曲がり率Afを算出し、算出されたそれぞれの繊維の曲がり率Afに基づいてそれぞれの繊維の曲がりを評価し、得られたそれぞれの繊維の曲がり評価結果に基づいて選択された繊維群の強度または弾性率(剛性)を評価するステップからなる如く構成したので、それぞれの繊維の曲がり率Afに基づいて選択された繊維群、即ち、成形品において強度または弾性率(剛性)を評価する範囲の強度または弾性率(剛性)を定量的に評価することができる。   In the computer-aided resin behavior analysis apparatus according to claim 1, a simulation program for analyzing the behavior of the fiber-mixed resin when molding in a mold under predetermined molding conditions includes a plurality of fibers of a plurality of fibers. When an analysis condition including at least the node F is input, at least a part of the plurality of fibers is selected, and each of the fibers constituting the fiber group selected based on the input analysis condition is selected. The bending rate Af is calculated, the bending of each fiber is evaluated based on the calculated bending rate Af of each fiber, and the strength of the fiber group selected based on the obtained bending evaluation result of each fiber or Since the elastic modulus (stiffness) is composed of steps, the fiber group selected based on the bending rate Af of each fiber, that is, the strength of the molded product. It may be assessed range strength or modulus of evaluating modulus (stiffness) and (stiffness) quantitatively.

請求項2に係るコンピュータ支援の樹脂挙動解析装置にあっては、シミュレーションプログラムは、得られた繊維群の強度または弾性率(剛性)評価結果に基づいて複数本の繊維の全体の強度または弾性率(剛性)を評価するステップを含む如く構成したので、上記した効果に加え、複数本の繊維の全体、即ち、複数本の繊維が混入された樹脂を金型内で成形するときの成形品の強度または弾性率(剛性)を定量的に評価することができる。   In the computer-aided resin behavior analysis apparatus according to claim 2, the simulation program calculates the overall strength or elastic modulus of the plurality of fibers based on the obtained strength or elastic modulus (rigidity) evaluation result of the fiber group. Since it is configured to include a step for evaluating (rigidity), in addition to the above-described effects, the whole of a plurality of fibers, that is, a molded product when a resin in which a plurality of fibers are mixed is molded in a mold. The strength or elastic modulus (rigidity) can be quantitatively evaluated.

請求項3に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がりを評価するステップは、算出されたそれぞれの繊維の曲がり率Afを所定のしきい値と比較して所定の範囲ごとに分類して評価するステップからなり、選択された繊維群の強度または弾性率(剛性)を評価するステップは、所定の範囲に属する繊維の割合を算出し、算出された繊維の割合に基づいて選択された繊維群の強度または弾性率(剛性)評価値を算出するステップからなる如く構成したので、上記した効果に加え、成形品において強度または弾性率(剛性)を評価する範囲の強度または弾性率(剛性)を、強度低下または弾性率(剛性)低下に寄与する曲がり率などの特定の曲がり率の繊維割合に基づいて定量的に評価することができる。   In the computer-aided resin behavior analysis apparatus according to claim 3, the step of evaluating the bending of each fiber is performed by comparing the calculated bending rate Af of each fiber with a predetermined threshold value. The method comprises the steps of classifying and evaluating each range, and the step of evaluating the strength or elastic modulus (rigidity) of the selected fiber group calculates the proportion of fibers belonging to a predetermined range, and calculates the calculated proportion of fibers. In addition to the above-described effects, the strength of the range in which the strength or elastic modulus (rigidity) is evaluated in the molded product, because it is configured to include the step of calculating the evaluation value of the strength or elastic modulus (rigidity) of the selected fiber group. Alternatively, the elastic modulus (rigidity) can be quantitatively evaluated based on the fiber ratio of a specific bending rate such as a bending rate that contributes to a decrease in strength or a reduction in elastic modulus (rigidity).

請求項4に係るコンピュータ支援の樹脂挙動解析装置にあっては、選択された繊維群の強度または弾性率(剛性)を評価するステップは、所定の範囲に属する繊維の割合を所定の範囲ごとに表示して評価するステップからなる如く構成したので、上記した効果に加え、成形品において強度または弾性率(剛性)を評価する範囲の強度または弾性率(剛性)を数値的に表示して評価することができ、成形条件を変更した結果を比較して評価することが容易になる。   In the computer-aided resin behavior analysis apparatus according to claim 4, the step of evaluating the strength or elastic modulus (rigidity) of the selected fiber group is configured to calculate a ratio of fibers belonging to a predetermined range for each predetermined range. Since it is configured to include the step of displaying and evaluating, in addition to the above-described effects, the strength or elastic modulus (rigidity) in the range in which the strength or elastic modulus (rigidity) is evaluated in the molded product is numerically displayed and evaluated. It is possible to compare and evaluate the results of changing the molding conditions.

請求項5に係るコンピュータ支援の樹脂挙動解析装置にあっては、選択された繊維群の強度または弾性率(剛性)を評価するステップは、算出された繊維群の強度または弾性率(剛性)評価値から予め設定された特性に基づいて選択された繊維群の強度または弾性率(剛性)を算出して評価するステップからなる如く構成したので、上記した効果に加え、成形品において強度または弾性率(剛性)を評価する範囲の強度または弾性率(剛性)を精度良く算出することができ、要求強度または弾性率(剛性)と比較して強度または弾性率(剛性)が十分か確認することができる。   In the computer-aided resin behavior analysis apparatus according to claim 5, the step of evaluating the strength or elastic modulus (rigidity) of the selected fiber group is performed by evaluating the calculated strength or elastic modulus (rigidity) of the fiber group. In addition to the above-described effects, the strength or elastic modulus of the molded product is obtained by calculating and evaluating the strength or elastic modulus (rigidity) of the fiber group selected based on the preset characteristics from the values. It is possible to calculate the strength or elastic modulus (rigidity) within the range to evaluate (rigidity) with high accuracy and confirm whether the strength or elastic modulus (rigidity) is sufficient compared to the required strength or elastic modulus (rigidity). it can.

請求項6に係るコンピュータ支援の樹脂挙動解析装置にあっては、金型内で所定の成形条件で成形するときの繊維混入樹脂の挙動を解析するシミュレーションプログラムが、複数本の繊維の複数個の節点Fを少なくとも含む解析条件が入力されるとき、複数本の繊維のうちの少なくとも一部の繊維群を選択し、入力された解析条件に基づいて選択された繊維群を構成するそれぞれの繊維の曲がり率Afを算出し、算出されたそれぞれの繊維の曲がり率Afを表示色に関連付けて前記それぞれの繊維の曲がりを評価し、得られたそれぞれの繊維の曲がり評価結果に基づいて選択された繊維群の強度または弾性率(剛性)を評価するステップからなる如く構成したので、それぞれの繊維の曲がり率Afに基づいて成形品において強度または弾性率(剛性)を評価する範囲の強度または弾性率(剛性)を視覚的に表示して評価することができ、成形条件を変更した結果を比較して評価することが容易になる。   In the computer-aided resin behavior analysis apparatus according to claim 6, a simulation program for analyzing the behavior of the fiber-mixed resin when molding in a mold under a predetermined molding condition includes a plurality of fibers of a plurality of fibers. When an analysis condition including at least the node F is input, at least a part of the plurality of fibers is selected, and each of the fibers constituting the fiber group selected based on the input analysis condition is selected. The bending rate Af is calculated, the calculated bending rate Af of each fiber is correlated with the display color, the bending of each of the fibers is evaluated, and the selected fiber based on the obtained bending evaluation result of each fiber Since it comprises the step which evaluates the intensity | strength or elastic modulus (rigidity) of a group, based on the bending rate Af of each fiber, strength or elastic modulus in a molded article Strength or elastic modulus in the range of evaluating stiffness) (stiffness) visually displayed to be able to evaluate, it is easy to evaluate and compare the result of changing the molding conditions.

請求項7に係るコンピュータ支援の樹脂挙動解析装置にあっては、シミュレーションプログラムは、得られた繊維群の強度または弾性率(剛性)評価結果に基づいて複数本の繊維の全体の強度または弾性率(剛性)を評価するステップを含む如く構成したので、上記した効果に加え、成形品の強度または弾性率(剛性)を定量的に評価することができる。   In the computer-aided resin behavior analysis apparatus according to claim 7, the simulation program calculates the overall strength or elastic modulus of the plurality of fibers based on the obtained strength or elastic modulus (rigidity) evaluation result of the fiber group. Since it is configured to include a step of evaluating (rigidity), in addition to the above-described effects, the strength or elastic modulus (rigidity) of the molded product can be quantitatively evaluated.

請求項8に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がり率Afを算出するステップは、前記複数個の節点Fの中の少なくとも一つの節点Fnから得られる評価長さに対する前記成形条件から予測される前記それぞれの繊維の曲がり率Afを算出するステップからなり、前記それぞれの繊維の曲がりを評価するステップは、前記算出された評価長さに対する曲がり率Afに基づいて前記それぞれの繊維の曲がりを評価するステップからなる如く構成したので、上記した効果に加え、評価長さに対する曲がり率Afに基づいてそれぞれの繊維の曲がりを評価することで細かいうねり(小さい曲がり)は小さい評価長さで、大きい曲がりは大きい評価長さで評価することも可能となり、それによって曲がりの形態ごとに定量的に評価することができ、細かいうねりと大きい曲がりで構成される複合的な曲がりも定量的に評価することができる。   In the computer-aided resin behavior analysis apparatus according to claim 8, the step of calculating the bending rate Af of each fiber is an evaluation length obtained from at least one node Fn among the plurality of nodes F. The step of calculating the bending rate Af of each of the fibers predicted from the molding conditions for the step of evaluating the bending of each of the fibers based on the bending rate Af with respect to the calculated evaluation length. Since it comprises so that it might consist of the step which evaluates the bending of each fiber, in addition to the above-mentioned effect, a small wave | undulation (small bending) is small by evaluating the bending of each fiber based on the bending rate Af with respect to evaluation length. With an evaluation length, a large bend can also be evaluated with a large evaluation length, which makes it possible to form a bend. Bets on can be quantitatively evaluated, it can also be quantitatively evaluated complex bends composed of bend larger and fine waviness.

また、それぞれの繊維の曲がりを適正に判定することが可能となることで、それぞれの繊維の曲がりが少なくなる条件を見つけることが容易となり、成形条件を変更して製品の強度または弾性率(剛性)を上げることも可能となる。   In addition, since it is possible to determine the bending of each fiber properly, it is easy to find the conditions that reduce the bending of each fiber, and the molding conditions can be changed to change the strength or elastic modulus (rigidity of the product). ) Can also be raised.

請求項9に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がり率Afを算出するステップは、複数個の節点Fの少なくとも一つの節点Fnを曲線または直線で接続して生成される三次元図形に基づいてそれぞれの繊維の曲がり率Afを算出するステップからなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afを簡易かつ精度よく算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 9, the step of calculating the bending rate Af of each fiber is generated by connecting at least one node Fn of the plurality of nodes F with a curve or a straight line. Since the configuration includes the step of calculating the bending rate Af of each fiber on the basis of the three-dimensional figure to be formed, in addition to the above-described effects, the bending rate Af of the fiber as a numerical value can be calculated easily and accurately. .

請求項10に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がり率Afを算出するステップは、複数個の節点Fの少なくとも一つの節点Fnを曲線または直線で接続して生成される三次元図形上の点を曲線または直線で接続して生成される第2の三次元図形に基づいてそれぞれの繊維の曲がり率Afを算出するステップからなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afを一層簡易に算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 10, the step of calculating the bending rate Af of each fiber is generated by connecting at least one node Fn of the plurality of nodes F with a curve or a straight line. The above-described effect is obtained because it comprises a step of calculating the bending rate Af of each fiber based on the second three-dimensional figure generated by connecting the points on the three-dimensional figure by a curve or a straight line. In addition, the bending rate Af of the fiber as a numerical value can be calculated more easily.

請求項11に係るコンピュータ支援の樹脂挙動解析装置にあっては、第2の三次元図形は三角形からなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afをより一層簡易に算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 11, since the second three-dimensional figure is configured to be a triangle, the fiber bending rate Af as a numerical value can be further simplified in addition to the effects described above. Can be calculated.

請求項12に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がり率Afを算出するステップは、節点Fnとそれに隣接する節点Fn-1,Fn+1とをそれぞれ結ぶ線FnFn-1FnFn+1を斜辺とすると共に、隣接する節点Fn-1,Fn+1を結ぶ線Fn-1Fn+1を底辺とする三角形を設定し、節点Fnあるいはその付近の任意の点から底辺に向けられた垂線の長さh1を算出し、算出された垂線の長さh1を斜辺FnFn-1FnFn+1の長さの合計値で除算してアスペクト比As1を算出し、算出されたアスペクト比As1に基づいてそれぞれの繊維の曲がり率Afを算出するステップからなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afを簡易に算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 12, the step of calculating the bending rate Af of each fiber includes a line FnFn connecting the node Fn and the adjacent nodes Fn-1, Fn + 1. -1FnFn + 1 is set as the hypotenuse, and a triangle having a base Fn-1Fn + 1 connecting adjacent nodes Fn-1 and Fn + 1 is set, and the node Fn or an arbitrary point near it is directed to the base The length h1 of the calculated perpendicular is calculated, the aspect ratio As1 is calculated by dividing the calculated length h1 of the perpendicular by the total length of the hypotenuses FnFn-1FnFn + 1, and the calculated aspect ratio As1. Since the configuration includes the step of calculating the bending rate Af of each fiber based on the above, in addition to the above effects, the bending rate Af of the fiber as a numerical value can be easily calculated.

請求項13に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がり率Afを算出するステップは、節点Fn+1とそれに隣接する節点Fn,Fn+2とをそれぞれ結ぶ線Fn+1FnFn+1Fn+2を斜辺とすると共に、隣接する節点Fn,Fn+2を結ぶ線FnFn+2を底辺とする三角形を設定し、節点Fn+1あるいはその付近の任意の点から底辺に向けられた垂線の長さh2を算出し、算出された垂線の長さh2を斜辺Fn+1FnFn+1Fn+2の長さの合計値で除算してアスペクト比As2を算出し、算出されたアスペクト比As2と節点Fnについて算出されたアスペクト比As1とに基づいてそれぞれの繊維の曲がり率Afを算出するステップからなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afを簡易かつ精度良く算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 13, the step of calculating the bending rate Af of each fiber is a line Fn connecting the node Fn + 1 and the nodes Fn and Fn + 2 adjacent thereto. Set a triangle with + 1FnFn + 1Fn + 2 as the hypotenuse and the base of line FnFn + 2 connecting adjacent nodes Fn and Fn + 2, from the node Fn + 1 or any nearby point toward the base The length h2 of the calculated perpendicular is calculated, the aspect ratio As2 is calculated by dividing the calculated length h2 of the perpendicular by the total length of the hypotenuses Fn + 1FnFn + 1Fn + 2, and the calculated aspect ratio Since it comprises the step which calculates the bending rate Af of each fiber based on As2 and the aspect ratio As1 calculated about the node Fn, in addition to the above-mentioned effect, the bending rate Af of the fiber as a numerical value can be simplified and Spirit It can be well calculated.

請求項14に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がり率Afを算出するステップは、節点Fn+2とそれに隣接する/しない節点Fn+2,Fn+mとをそれぞれ結ぶ線Fn+2Fn+1Fn+2Fn+mを斜辺とすると共に、隣接する節点Fn+1,Fn+mを結ぶ線Fn+1Fn+mを底辺とする三角形を設定し、節点Fn+2あるいはその付近の任意の点から底辺に向けられた垂線の長さhmを算出し、算出された垂線の長さhmを斜辺Fn+2Fn+1Fn+2Fn+mの長さの合計値で除算してアスペクト比Asmを算出し、算出されたアスペクト比Asmと節点Fnと節点Fn+1について算出されたアスペクト比As1,As2とに基づいてそれぞれの繊維の曲がり率Afを算出するステップからなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afを簡易かつ一層精度良く算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 14, the step of calculating the bending rate Af of each fiber includes the node Fn + 2 and the nodes Fn + 2 and Fn + m adjacent to / not from the node Fn + 2. Set a triangle with the line Fn + 2Fn + 1Fn + 2Fn + m connecting to each other as the hypotenuse and the line Fn + 1Fn + m connecting adjacent nodes Fn + 1 and Fn + m as the base, and the node Fn + 2 or Calculate the length hm of the perpendicular from the arbitrary point in the vicinity to the base, and divide the calculated length hm by the total length of the hypotenuse Fn + 2Fn + 1Fn + 2Fn + m The aspect ratio Asm is calculated, and the bending ratio Af of each fiber is calculated based on the calculated aspect ratio Asm and the aspect ratios As1 and As2 calculated for the node Fn and the node Fn + 1. So, in addition to the effects described above, the fiber as a numerical value And simple bending ratio Af can be calculated more accurately.

請求項15に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がり率Afを算出するステップは、節点Fnとそれに隣接する節点Fn-1,Fn+1とをそれぞれ結ぶ線FnFn-1FnFn+1上のいずれかの任意の点Pと隣接する節点Fn-1,Fn+2を結ぶ線PFn-1PFn+2を斜辺とすると共に、節点Fn-1,Fn+2を結ぶ線Fn-1Fn+2を底辺とする三角形を設定し、任意の点Pから底辺に向けられた垂線の長さh1を算出し、算出された垂線の長さh1を斜辺PFn-1PFn+2の長さの合計値で除算してアスペクト比As1を算出し、次いで他の節点についてアスペクト比Asnを算出し、算出されたアスペクト比As1,Asnに基づいてそれぞれの繊維の曲がり率Afを算出するステップからなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afを簡易かつ一層精度良く算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 15, in the step of calculating the bending rate Af of each fiber, the line FnFn connecting the node Fn and the adjacent nodes Fn-1, Fn + 1, respectively. A line PFn-1PFn + 2 connecting any arbitrary point P on -1FnFn + 1 and the adjacent nodes Fn-1, Fn + 2 is a hypotenuse and a line Fn connecting the nodes Fn-1, Fn + 2 Set a triangle with -1Fn + 2 as the base, calculate the length h1 of the perpendicular from any point P to the base, and calculate the length h1 of the perpendicular as the length of the hypotenuse PFn-1PFn + 2 The aspect ratio As1 is calculated by dividing by the sum of the values, then the aspect ratio Asn is calculated for the other nodes, and the bending ratio Af of each fiber is calculated based on the calculated aspect ratios As1 and Asn. In addition to the effects described above, And simple bending ratio Af of the fiber as the value can be calculated more accurately.

請求項16に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がり率Afを算出するステップは、算出された値の平均値に基づいてそれぞれの繊維の曲がり率Afを算出するステップからなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afを簡易かつ一層精度良く算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 16, the step of calculating the bending rate Af of each fiber calculates the bending rate Af of each fiber based on the average value of the calculated values. Since the configuration is made up of steps, in addition to the above-described effects, the fiber bending rate Af as a numerical value can be calculated easily and more accurately.

請求項17に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がり率Afを算出するステップは、算出された値の最大値に基づいてそれぞれの繊維の曲がり率Afを算出するステップからなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afを簡易かつ一層精度良く算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 17, the step of calculating the bending rate Af of each fiber calculates the bending rate Af of each fiber based on the maximum value of the calculated values. Since the configuration is made up of steps, in addition to the above-described effects, the fiber bending rate Af as a numerical value can be calculated easily and more accurately.

請求項18に係るコンピュータ支援の樹脂挙動解析装置にあっては、それぞれの繊維の曲がり率Afを算出するステップは、金型のキャビティにおける繊維の三次元空間での座標位置に基づいてそれぞれの繊維の曲がり率Afを算出するステップからなる如く構成したので、数値としての繊維の曲がり率Afを一層簡易に算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 18, the step of calculating the bending rate Af of each fiber includes the step of calculating each fiber based on the coordinate position in the three-dimensional space of the fiber in the cavity of the mold. Therefore, the fiber bending rate Af as a numerical value can be calculated more easily.

請求項19に係るコンピュータ支援の樹脂挙動解析装置にあっては、評価長さは三次元図形の始点から終点までの長さからなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afを一層精度良く算出することができる。   In the computer-aided resin behavior analysis apparatus according to claim 19, since the evaluation length is configured to be the length from the start point to the end point of the three-dimensional figure, in addition to the above-described effects, the numerical value of the fiber The bending rate Af can be calculated with higher accuracy.

請求項20に係るコンピュータ支援の樹脂挙動解析装置にあっては、シミュレーションプログラムは、分析器を介して得られる複数本の繊維の複数個の節点Fを少なくとも含むデータが入力されるとき、複数本の繊維のうちの少なくとも一部の繊維群を選択し、入力されたデータに基づいて選択された繊維群を構成するそれぞれの繊維の曲がり率Afを算出し、算出されたそれぞれの繊維の曲がり率Afに基づいてそれぞれの繊維の曲がりを評価し、得られたそれぞれの繊維の曲がり評価結果に基づいて選択された繊維群の強度または弾性率を評価するステップからなる如く構成したので、実際の成形品を分析器、例えばX線CTスキャンで分析して得られるデータ(実際の成形品中の繊維の形状、節点)をシミュレーションデータと同様の手法で解析して比較することができる。   In the computer-aided resin behavior analysis apparatus according to claim 20, when the simulation program receives data including at least a plurality of nodes F of a plurality of fibers obtained through the analyzer, a plurality of simulation programs are input. And selecting a fiber group of at least some of the fibers, calculating a bending rate Af of each fiber constituting the selected fiber group based on the input data, and calculating the calculated bending rate of each fiber Since the bending of each fiber is evaluated based on Af and the strength or elastic modulus of the selected fiber group is evaluated based on the obtained bending evaluation result of each fiber, the actual molding is performed. Data obtained by analyzing the product with an analyzer such as an X-ray CT scan (the shape of the fiber in the actual molded product, nodes) is the same as the simulation data. It can be compared and analyzed by law.

この発明の実施形態に係るコンピュータ支援の樹脂挙動解析装置を用いて行われる製品設計から量産までの工程を示す説明図である。It is explanatory drawing which shows the process from product design performed using the computer-aided resin behavior analysis apparatus based on embodiment of this invention to mass production. 図1に示す装置のコンピュータに格納されるシミュレーションプログラムの処理を示すフロー・チャートである。It is a flowchart which shows the process of the simulation program stored in the computer of the apparatus shown in FIG. 図2フロー・チャートのS38の処理のサブ・ルーチン・フロー・チャートである。3 is a sub-routine flow chart of the process of S38 of the flow chart. 図3の曲がりの評価で対象とする1本の繊維を、繊維01から04までの4例について示す説明図である。It is explanatory drawing which shows one fiber made into object by evaluation of the bending | flexion of FIG. 3 about four examples from fiber 01 to 04. FIG. 図3フロー・チャートのS102の処理のサブ・ルーチン・フロー・チャートである。3 is a sub-routine flow chart of the process of S102 of the flow chart. 図5フロー・チャートの処理の説明図である。5 is an explanatory diagram of the processing of the flow chart. 同様に図5フロー・チャートの処理の説明図である。Similarly, it is explanatory drawing of the process of FIG. 5 flowchart. 同様に図5フロー・チャートの処理の説明図である。Similarly, it is explanatory drawing of the process of FIG. 5 flowchart. 同様に図5フロー・チャートの処理の説明図である。Similarly, it is explanatory drawing of the process of FIG. 5 flowchart. 同様に図5フロー・チャートの処理の説明図である。Similarly, it is explanatory drawing of the process of FIG. 5 flowchart. 同様に図5フロー・チャートの処理の説明図である。Similarly, it is explanatory drawing of the process of FIG. 5 flowchart. 同様に図5フロー・チャートの処理の説明図である。Similarly, it is explanatory drawing of the process of FIG. 5 flowchart. この実施形態の変形例を示す、図5と同様のフロー・チャートである。It is a flowchart similar to FIG. 5 which shows the modification of this embodiment. 図13フロー・チャートの処理の説明図である。13 is an explanatory diagram of the processing of the flow chart. 同様に図13フロー・チャートの処理の説明図である。Similarly, it is explanatory drawing of the process of FIG. 13 flowchart. この実施形態の変形例を示す、図5と同様のフロー・チャートである。It is a flowchart similar to FIG. 5 which shows the modification of this embodiment. 図16フロー・チャートの処理の説明図である。FIG. 16 is an explanatory diagram of processing of the flowchart. 同様に図16フロー・チャートの処理の説明図である。Similarly, it is explanatory drawing of the process of FIG. 16 flowchart. この実施形態の変形例を示す、図16に類似するフロー・チャートである。It is a flowchart similar to FIG. 16 which shows the modification of this embodiment. 同様にこの実施形態の変形例を示す、図16に類似するフロー・チャートである。FIG. 17 is a flow chart similar to FIG. 16, similarly showing a modification of this embodiment. 図20フロー・チャートの処理の説明図である。FIG. 20 is an explanatory diagram of the processing of the flow chart. この実施形態の変形例を示す説明図である。It is explanatory drawing which shows the modification of this embodiment. 図16などの処理によって得られる図4に示す4種の繊維についての算出結果を示す説明図である。It is explanatory drawing which shows the calculation result about four types of fibers shown in FIG. 4 obtained by processes, such as FIG. 図3フロー・チャートのS104の処理のサブ・ルーチン・フロー・チャートである。3 is a sub-routine flow chart of the process of S104 of the flow chart of FIG. 図2フロー・チャートのS44の処理のサブ・ルーチン・フロー・チャートである。3 is a sub-routine flow chart of the process of S44 of the flow chart of FIG. 図25フロー・チャートの強度または弾性率(剛性)算出マップの例を示す説明図である。FIG. 25 is an explanatory diagram showing an example of a strength or elastic modulus (rigidity) calculation map of the flow chart. 図25の処理によって得られる繊維割合の算出結果を示す説明図である。It is explanatory drawing which shows the calculation result of the fiber ratio obtained by the process of FIG. この実施形態の変形例を示す、図3と同様のフロー・チャートである。It is the same flow chart as FIG. 3 which shows the modification of this embodiment. この実施形態に係る解析結果の一例を示す写真である。It is a photograph which shows an example of the analysis result which concerns on this embodiment. 同様にこの実施形態に係る解析結果の一例を示す写真である。It is a photograph which similarly shows an example of the analysis result which concerns on this embodiment.

以下、添付図面を参照してこの発明に係るコンピュータ支援の樹脂挙動解析装置を実施するための形態について説明する。   DESCRIPTION OF EMBODIMENTS Hereinafter, an embodiment for implementing a computer-aided resin behavior analysis apparatus according to the present invention will be described with reference to the accompanying drawings.

図1はこの発明の実施形態に係るコンピュータ支援の樹脂挙動解析装置を用いて行われる製品設計から量産までの工程を示す説明図である。   FIG. 1 is an explanatory diagram showing steps from product design to mass production performed using a computer-aided resin behavior analysis apparatus according to an embodiment of the present invention.

図1において符号10は樹脂挙動解析装置を示し、装置10はコンピュータ12を備え、コンピュータによる開発支援装置(Computer Aided Engineering)、即ち、コンピュータ支援の樹脂挙動解析装置として構成される。   In FIG. 1, reference numeral 10 denotes a resin behavior analysis apparatus. The apparatus 10 includes a computer 12 and is configured as a computer-aided development support apparatus (Computer Aided Engineering), that is, a computer-aided resin behavior analysis apparatus.

コンピュータ12はバスを介して互いにデータ通信可能に構成されるCPU(Central Processing Unit)、メモリ、入出力回路などからなると共に、ディスプレイ12aと、キーボード、マウス、タッチパネルなどからなる入力デバイス12bを備える。   The computer 12 includes a central processing unit (CPU), a memory, an input / output circuit, and the like configured to be able to communicate with each other via a bus, and includes a display 12a and an input device 12b including a keyboard, a mouse, a touch panel, and the like.

コンピュータ12のメモリには金型14内の樹脂16の挙動を解析するための対話形式のシミュレーションプログラム20が格納され、CPUはメモリに格納されたシミュレーションプログラム20を実行する。ディスプレイ12aはその結果を表示すると共に、入力デバイス12bは設計者(エンジニア)の操作や指示を受付ける。   An interactive simulation program 20 for analyzing the behavior of the resin 16 in the mold 14 is stored in the memory of the computer 12, and the CPU executes the simulation program 20 stored in the memory. The display 12a displays the result, and the input device 12b accepts a designer (engineer) operation and instruction.

シミュレーションプログラム20の解析対象である金型14について説明すると、図1に示すように、金型14は上型14aと下型14bを備え、その間にキャビティ14cが形成される。キャビティ14cには樹脂16が投入される。樹脂16は炭素繊維などからなる連続繊維および不連続な長繊維(以下「繊維」という)16aが混入されてシート状を呈する。   The mold 14 to be analyzed by the simulation program 20 will be described. As shown in FIG. 1, the mold 14 includes an upper mold 14a and a lower mold 14b, and a cavity 14c is formed therebetween. Resin 16 is put into the cavity 14c. The resin 16 has a sheet shape in which continuous fibers made of carbon fiber or the like and discontinuous long fibers (hereinafter referred to as “fibers”) 16 a are mixed.

コンピュータ12は、設計者による入力デバイス12bの操作によって所定の成形条件が入力されるとき、樹脂16を金型キャビティ14cに投入して所定の成形条件で流動させて成形品(製品あるいは半製品)を成形、より具体的にはプレス成形するときの樹脂16の挙動を格納された対話形式のシミュレーションプログラム20によるシミュレーションを介して解析する。   When a predetermined molding condition is input by an operation of the input device 12b by the designer, the computer 12 puts the resin 16 into the mold cavity 14c and causes it to flow under the predetermined molding condition to form a molded product (product or semi-finished product). , More specifically, the behavior of the resin 16 at the time of press molding is analyzed through simulation by the stored simulation program 20.

樹脂16の挙動の解析は具体的には製品設計から量産までの工程の一環として行われ、設計者(エンジニア)は、入力デバイス12bを介してデータを入力し、シミュレーションプログラム20に格納された指示に従い、対話形式で製品モデル22を設計する。   Specifically, the analysis of the behavior of the resin 16 is performed as part of the process from product design to mass production, and the designer (engineer) inputs data via the input device 12b and stores the instructions stored in the simulation program 20 The product model 22 is designed in an interactive manner.

装置10においては、金型14を用いて製品製造を行う場合、設計者は先ず製品設計工程で製品モデル22を設計し、作成した製品モデル22を使用して金型設計工程で金型モデル24を設計する。   In the apparatus 10, when a product is manufactured using the mold 14, the designer first designs the product model 22 in the product design process, and uses the created product model 22 to mold the mold model 24 in the mold design process. To design.

次いで、設計者は、作成した金型モデル24を使用して金型加工データを作成し、そのデータを使用してNC加工装置26などによって金型14を製造し、製造した金型14を使用して成形を行うことによって製品(あるいは半製品)30を製造する。   Next, the designer creates die machining data using the created die model 24, uses the data to produce the die 14 using the NC processing device 26, and uses the produced die 14. Then, the product (or semi-finished product) 30 is manufactured by molding.

図2はこの発明の実施形態によるコンピュータ支援の樹脂挙動解析装置10の処理(動作)を示すフロー・チャート、より具体的には装置10のコンピュータ12に格納されたシミュレーションプログラム20の処理(ステップ)を示すフロー・チャートである。   FIG. 2 is a flowchart showing the processing (operation) of the computer-aided resin behavior analysis apparatus 10 according to the embodiment of the present invention, more specifically, the processing (step) of the simulation program 20 stored in the computer 12 of the apparatus 10. It is a flowchart which shows.

以下説明すると、S10において樹脂充填解析を開始する。具体的には、S12において設計者によって所定の成形条件が入力される。所定の成形条件は、樹脂(シート状)16の物性(粘度、熱伝導率など)、初期樹脂充填(チャージ)位置と樹脂16の寸法、金型材特性、プレス機特性(最大型締め力Nmax、型締め速度V1など)と、樹脂16に混入される繊維16aの形状、分割長さなどを含む。   To explain below, resin filling analysis is started in S10. Specifically, predetermined molding conditions are input by the designer in S12. Predetermined molding conditions include physical properties (viscosity, thermal conductivity, etc.) of the resin (sheet-like) 16, initial resin filling (charge) position and dimensions of the resin 16, die material characteristics, press machine characteristics (maximum clamping force Nmax, The mold clamping speed V1 and the like, and the shape and split length of the fiber 16a mixed in the resin 16.

次いでS14に進んで初期値を設定(リセット)する。即ち、時間tを0とすると共に、型締め速度Vを成形条件のV1とする。S12とS14の処理が条件入力に相当する。   Next, the process proceeds to S14 to set (reset) the initial value. That is, the time t is set to 0, and the mold clamping speed V is set to V1 of the molding condition. The processing of S12 and S14 corresponds to condition input.

次いでS16に進んで時間tを単位時間Δtだけインクリメントし、S18に進んで型締め速度V1一定としたときの単位時間Δtの間のプレスによる樹脂16の流動を解析する。即ち、金型キャビティ14cに相当する三次元空間における樹脂16の流動速度分布を算出する。   Next, the process proceeds to S16, where the time t is incremented by the unit time Δt, and the process proceeds to S18, where the flow of the resin 16 due to pressing during the unit time Δt when the mold clamping speed V1 is constant is analyzed. That is, the flow velocity distribution of the resin 16 in the three-dimensional space corresponding to the mold cavity 14c is calculated.

次いで、S20に進み、時間tからt+Δtまでの間の圧縮力Nを算出し、S22に進み、算出された圧縮力Nが最大型締め力Nmax未満か否か判断する。S22で肯定されるときはS16に戻って上記した処理を繰り返す。S16からS22までの処理が充填解析に相当する。   Next, the process proceeds to S20, where the compression force N from time t to t + Δt is calculated, and the process proceeds to S22, where it is determined whether or not the calculated compression force N is less than the maximum clamping force Nmax. When the result in S22 is affirmative, the process returns to S16 and the above processing is repeated. The processing from S16 to S22 corresponds to the filling analysis.

一方、S22で否定されるときは圧縮力Nが最大型締め力Nmaxに達したと判断できるため、S24に進んで圧縮力Nを一定としたときの単位時間Δtの間のプレスによる樹脂16の流動を解析する。即ち、圧縮力Nを最大値としたときの単位時間Δtの間に充填される際の樹脂16の流動速度分布を算出する。   On the other hand, when the result in S22 is negative, it can be determined that the compression force N has reached the maximum clamping force Nmax. Therefore, the process proceeds to S24 and the resin 16 is pressed by the unit during the unit time Δt when the compression force N is constant. Analyze flow. That is, the flow velocity distribution of the resin 16 when it is filled during the unit time Δt when the compression force N is the maximum value is calculated.

次いで、S26に進んで時間tからt+Δtまでの間の型締め速度Vを算出し、S28に進み、算出された型締め速度Vが0より大きいか否か判断する。S28で肯定されるときはS30に進んで時間tを単位時間ΔtだけインクリメントしてS24に戻り、上記した処理を繰り返す。   Next, the process proceeds to S26, where the mold clamping speed V from time t to t + Δt is calculated, and the process proceeds to S28, where it is determined whether or not the calculated mold clamping speed V is greater than zero. When the result in S28 is affirmative, the program proceeds to S30, in which the time t is incremented by the unit time Δt, the process returns to S24, and the above processing is repeated.

他方、S28で否定されるときは型締め速度Vが0になって樹脂16の充填が終了したと判断できるため、S32に進んで樹脂充填解析を終了する。S24からS32までの処理も充填解析に相当する。   On the other hand, when the result in S28 is negative, it can be determined that the mold clamping speed V has become 0 and the filling of the resin 16 has been completed, so the process proceeds to S32 and the resin filling analysis is completed. The processing from S24 to S32 also corresponds to the filling analysis.

次いでS34に進み、繊維挙動解析を開始する。   Next, in S34, the fiber behavior analysis is started.

即ち、S36に進んで時間tを0にリセットし、次いでS38に進み、各繊維16a、即ち、繊維16aの1本ごとの挙動を解析する。具体的には、時間t=0(t:樹脂16の充填解析時間)において、各繊維16aに三次元座標を付与し、付与された座標とS18,S24で算出された樹脂16の流動速度分布とに基づき、樹脂16の流動に伴う各繊維16aの挙動を解析する、即ち、各繊維16aの各節点Fの座標を算出して各繊維16aの曲がりを評価することで挙動を解析する。   That is, the process proceeds to S36, the time t is reset to 0, and then the process proceeds to S38, in which the behavior of each fiber 16a, that is, each fiber 16a, is analyzed. Specifically, at time t = 0 (t: resin 16 filling analysis time), a three-dimensional coordinate is assigned to each fiber 16a, and the flow rate distribution of the resin 16 calculated in S18 and S24 is given. Based on the above, the behavior of each fiber 16a accompanying the flow of the resin 16 is analyzed, that is, the behavior is analyzed by calculating the coordinates of each node F of each fiber 16a and evaluating the bending of each fiber 16a.

次いで、S40に進み、時間を単位時間Δtだけインクリメントし、S42に進み、時間tが樹脂16の充填終了時間に達したか否か判断する。S42で否定されるときはS38に戻る一方、肯定されるときはS44に進んで成形品の構造強度または弾性率(剛性)を評価し、S46に進んで繊維挙動解析を終了する。S34からS46までの処理が繊維挙動解析に相当する。   Next, the process proceeds to S40, the time is incremented by the unit time Δt, the process proceeds to S42, and it is determined whether or not the time t has reached the filling end time of the resin 16. When the result in S42 is negative, the process returns to S38. When the result is positive, the process proceeds to S44 to evaluate the structural strength or elastic modulus (rigidity) of the molded product, and the process proceeds to S46 to complete the fiber behavior analysis. The processing from S34 to S46 corresponds to fiber behavior analysis.

図3は図2フロー・チャートのS38の処理のサブ・ルーチン・フロー・チャートである。S100において、樹脂16に混入される複数本の繊維16aの複数個の節点Fを少なくとも含む解析条件が設計者によって入力デバイス12bから入力されるとき、複数本の繊維16aのうち少なくとも一部の繊維群を選択、即ち、成形品において強度または弾性率(剛性)を評価する範囲に相当する繊維群を選択する。次いで、S102において、選択された繊維群を構成するそれぞれの繊維16aに着目し、それぞれの繊維16aの複数個の節点Fの中の少なくとも一つの節点Fnから得られる評価長さに対する、上記した成形条件から予測される1本の繊維16aの曲がり率Afを算出する。解析条件は繊維16aの長さ、節点Fの間の距離、製品30の形状などからなる。   FIG. 3 is a sub-routine flow chart of the process of S38 in the flow chart of FIG. In S100, when an analysis condition including at least the plurality of nodes F of the plurality of fibers 16a mixed in the resin 16 is input from the input device 12b by the designer, at least some of the plurality of fibers 16a. A group is selected, that is, a fiber group corresponding to a range in which strength or elastic modulus (rigidity) is evaluated in a molded product is selected. Next, in S102, paying attention to each fiber 16a constituting the selected fiber group, the above-described forming for the evaluation length obtained from at least one node Fn among the plurality of nodes F of each fiber 16a. The bending rate Af of one fiber 16a predicted from the conditions is calculated. The analysis conditions include the length of the fiber 16a, the distance between the nodes F, the shape of the product 30, and the like.

図4は曲がりの評価で対象とする1本の繊維16aを、繊維01から04までの4例について示す説明図である。この実施形態では4例のいずれの場合でも、繊維16aの長さ(繊維長)は20mm、分割数10個、分割長さ一定とする。即ち、4例の繊維16aはいずれもその長さを10個に均等に分割して得られる11個の節点Fを備える。実施形態では4例の全てについて繊維16aの1本ごとに曲がりを評価するようにした。   FIG. 4 is an explanatory diagram showing one example of the fibers 16a to be subjected to the bending evaluation for the fibers 01 to 04. In this embodiment, in any of the four examples, the length of the fiber 16a (fiber length) is 20 mm, the number of divisions is 10, and the division length is constant. That is, each of the four examples of the fibers 16a includes eleven nodes F obtained by equally dividing the length into ten. In the embodiment, the bending is evaluated for each of the fibers 16a in all four examples.

繊維16aの曲がりを評価するには繊維16aの形状(位置)を把握する必要があるが、図1の構成にあってはシミュレーションにおいて樹脂16が金型14のキャビティ14cに投入されるとき、そこに混入された多数の繊維16aのそれぞれは金型キャビティ14cにおける三次元空間での座標位置で指定されることから、繊維16aの位置はその座標位置を流用して特定(検出)する。   In order to evaluate the bending of the fiber 16a, it is necessary to grasp the shape (position) of the fiber 16a. In the configuration of FIG. 1, when the resin 16 is introduced into the cavity 14c of the mold 14 in the simulation, there is Since each of the many fibers 16a mixed in is designated by the coordinate position in the three-dimensional space in the mold cavity 14c, the position of the fiber 16a is specified (detected) using the coordinate position.

図5は図3フロー・チャートのS102の処理のサブ・ルーチン・フロー・チャート、図6から図12は図5フロー・チャートの処理の説明図である。   FIG. 5 is a sub-routine flow chart of the process of S102 of FIG. 3 and FIGS. 6 to 12 are explanatory diagrams of the process of the flowchart of FIG.

まずS200において複数個の節点Fの少なくとも一つの節点Fnを曲線または直線で接続して三次元図形を生成する。各節点Fは、図2フロー・チャートで規定される時間tにおける、各繊維16aの各節点Fの座標としてコンピュータ12のメモリに格納されたデータから読み出される。   First, in S200, at least one node Fn of the plurality of nodes F is connected by a curve or a straight line to generate a three-dimensional figure. Each node F is read from the data stored in the memory of the computer 12 as the coordinates of each node F of each fiber 16a at the time t specified in the flow chart of FIG.

次いでS202において生成された三次元図形に基づいて繊維16aの曲がり率Afを算出する。   Next, the bending rate Af of the fiber 16a is calculated based on the three-dimensional figure generated in S202.

図6に示すように、曲がり率Afは繊維16aの「評価長さ」と複数個の節点Fから選択された節点Fnを接続して生成される三次元図形の「評価高さ」の比に基づいて算出する。また、評価長さは選択された節点Fnの全部または一部を曲線または直線で接続して生成される三次元図形の始点から終点までの長さとして算出する。   As shown in FIG. 6, the bending rate Af is the ratio of the “evaluation length” of the fiber 16a and the “evaluation height” of the three-dimensional figure generated by connecting the node Fn selected from the plurality of nodes F. Calculate based on The evaluation length is calculated as the length from the start point to the end point of the three-dimensional figure generated by connecting all or part of the selected nodes Fn with curves or straight lines.

尚、図6およびそれ以降の図に記載される図形はxyz三次元空間における三次元図形をxy平面上に投影して得た二次元図形であるが、図示の便宜のため、この二次元図形を請求項でいう「三次元図形」を意味するものとして使用する。   6 and the subsequent figures are two-dimensional figures obtained by projecting a three-dimensional figure in the xyz three-dimensional space onto the xy plane. For convenience of illustration, this two-dimensional figure is used. Is used to mean a “three-dimensional figure” in the claims.

図7に示すように、繊維16aの曲がりは細かいうねり(小さい曲がり)、大きい曲がり、小さい曲がりと大きい曲がりの複合曲がりに大きく分類することができる。   As shown in FIG. 7, the bending of the fiber 16a can be roughly classified into fine undulations (small bending), large bendings, small bendings, and compound bendings of large bendings.

このように形態の異なる曲がりを定量的に評価するには、形態ごとに切り分けて評価する必要がある。そのため、曲がりの大きさ(曲がり率Af)を特定の長さ(評価長さ)ごとに評価することとした。即ち、小さい曲がりは小さい評価長さで、大きい曲がりは大きい評価長さで曲がり率Afを算出することとした。   In order to quantitatively evaluate the bends having different forms as described above, it is necessary to evaluate them separately for each form. For this reason, the magnitude of bending (bending rate Af) is evaluated for each specific length (evaluation length). That is, the bending rate Af is calculated with a small evaluation length for small bends and a large evaluation length for large bends.

上記したように、この実施形態およびその変形例において、評価長さは、三次元図形の始点と終点との直線距離ではなく、始点から終点までの繊維16aの長さとした。これは、曲がりの強いものと弱いものの横軸(評価長さ)を揃えるためである。   As described above, in this embodiment and its modifications, the evaluation length is not the linear distance between the start point and the end point of the three-dimensional figure, but the length of the fiber 16a from the start point to the end point. This is to align the horizontal axis (evaluation length) of the strong and weak bends.

続いて、図8を参照して図6の評価高さの算出方法について説明する。繊維16aの曲がりの形状を正弦曲線と仮定するとき、評価高さhは正弦曲線の振幅として算出することができる。   Next, the evaluation height calculation method of FIG. 6 will be described with reference to FIG. When it is assumed that the bending shape of the fiber 16a is a sine curve, the evaluation height h can be calculated as the amplitude of the sine curve.

即ち、三次元図形の始点を(0,0,0)、終点を(π,0,0)とする新たな三次元座標空間に点H(π/2,h,0)を取り、各節点Fnと曲線y=h・sin(x)(z=0)との距離の総和が最小になるよう、最小二乗法によって評価高さhを求める。   That is, a point H (π / 2, h, 0) is taken in a new three-dimensional coordinate space where the start point of the three-dimensional figure is (0, 0, 0) and the end point is (π, 0, 0), and each node The evaluation height h is obtained by the least square method so that the sum of the distances between Fn and the curve y = h · sin (x) (z = 0) is minimized.

しかしながら、上記のように正弦曲線を仮定して評価高さhを算出する方法は計算負荷が高く現実的ではない。そのため、評価高さhは、すでに座標が確定している節点Fnを連結して生成される三次元図形の形状に基づいて算出することとした。   However, the method of calculating the evaluation height h assuming a sine curve as described above has a high calculation load and is not practical. Therefore, the evaluation height h is calculated based on the shape of the three-dimensional figure generated by connecting the nodes Fn whose coordinates are already determined.

具体的には、図9に示すような形状に基づいて算出することとした。即ち、節点Fnを直線で接続して三次元図形を生成し、生成された三次元図形上を点Pが始点から終点までの全部または一部を動くとき、点Pから始点と終点とを結んだ直線に下した垂線が通過する部分の面積Sを、垂線の足が通過した部分の長さLhで除算して、評価高さhを算出するように構成した。   Specifically, the calculation is based on the shape as shown in FIG. That is, the node Fn is connected by a straight line to generate a three-dimensional figure, and when the point P moves all or part of the generated three-dimensional figure from the start point to the end point, the start point and the end point are connected from the point P. The evaluation height h is calculated by dividing the area S of the portion through which the perpendicular line passes through the length Lh of the portion through which the vertical leg passes.

より具体的には、隣接する二つの節点間を一つの区間とし、各区間について点Pを隣接する二つの節点Fnを通る直線上の点として数式で表す。次いで、始点と終点を通る直線と点Pの距離を数式で表し、それを積分することでその区間の面積を求め、すべての区間の面積の合計値を求める。   More specifically, between two adjacent nodes is defined as one section, and for each section, the point P is expressed as a point on a straight line passing through the two adjacent nodes Fn. Next, the distance between the straight line passing through the start point and the end point and the point P is expressed by a mathematical expression, and by integrating the distance, the area of the section is obtained, and the total value of the areas of all the sections is obtained.

尚、数式で表すことができれば、節点Fnは曲線で接続しても良いことになる。節点Fの番号はコンピュータ12のメモリに格納されたデータから繊維16aの節点(あるいは一端(始端)と他端(終端))の座標(x,y,zの三次元座標)を読み出し、ナンバリングすることで行う。   If it can be expressed by a mathematical expression, the node Fn may be connected by a curve. As for the number of the node F, the coordinates (three-dimensional coordinates of x, y, z) of the node (or one end (start end) and the other end (end)) of the fiber 16a are read from the data stored in the memory of the computer 12 and numbered. Do that.

図10にこの実施形態の変形例を示す。この変形例では、節点Fnを直線で接続して三次元図形を生成し、生成された三次元図形上を点Pが始点から終点までの全部または一部を動くとき、点Pから始点と終点とを結んだ直線に下した垂線が通過する部分の重心Gから始点と終点とを結んだ直線に下した垂線の長さとして、評価高さhを算出するように構成した。節点Fnは曲線で接続しても良い。   FIG. 10 shows a modification of this embodiment. In this modified example, the node Fn is connected by a straight line to generate a three-dimensional figure, and when the point P moves all or part of the point from the start point to the end point on the generated three-dimensional figure, the start point and the end point from the point P The evaluation height h is calculated as the length of the perpendicular drawn from the center of gravity G of the part through which the perpendicular drawn from the straight line connecting the starting point and the ending point passes. The node Fn may be connected by a curve.

上記した垂線の中点と重心Gとの各軸方向の変位と垂線の長さとの積は、三次元図形の始点と終点とを結んだ直線の方向に積分するとき、すべて0となる。従って、始点を(0,0,0)、終点をx軸上に置いた新たな三次元座標空間において、各軸方向の変位と垂線の長さとの積をxの関数として表し、それぞれx軸方向に積分したとき0となるように重心Gの座標を求める。節点Fnは曲線で接続しても良い。   The product of the displacement of each axial direction between the midpoint of the perpendicular line and the center of gravity G and the length of the perpendicular line becomes 0 when integrating in the direction of a straight line connecting the start point and the end point of the three-dimensional figure. Therefore, in a new three-dimensional coordinate space where the start point is (0, 0, 0) and the end point is on the x-axis, the product of the displacement in the direction of each axis and the length of the perpendicular is expressed as a function of x, respectively. The coordinates of the center of gravity G are obtained so that it becomes 0 when integrated in the direction. The node Fn may be connected by a curve.

図11にこの実施形態の変形例を示す。この変形例では、節点Fnを直線で接続して三次元図形を生成し、生成された三次元図形上に複数個の点Pを生成し、各点Pから始点と終点とを結んだ直線に下した垂線の長さの最大値として、評価高さhを算出するように構成した。節点Fnは曲線で接続しても良い。   FIG. 11 shows a modification of this embodiment. In this modified example, the node Fn is connected by a straight line to generate a three-dimensional figure, a plurality of points P are generated on the generated three-dimensional figure, and a straight line connecting the start point and the end point from each point P. The evaluation height h was calculated as the maximum value of the length of the dropped perpendicular. The node Fn may be connected by a curve.

図12にこの実施形態の変形例を示す。この変形例では、節点Fnを直線で接続して三次元図形を生成し、生成された三次元図形上に複数個の点Pを生成し、各点Pから始点と終点とを結んだ直線に下した垂線の長さの平均値として、評価高さhを算出するように構成した。節点Fnは曲線で接続しても良い。   FIG. 12 shows a modification of this embodiment. In this modified example, the node Fn is connected by a straight line to generate a three-dimensional figure, a plurality of points P are generated on the generated three-dimensional figure, and a straight line connecting the start point and the end point from each point P. The evaluation height h was calculated as an average value of the lengths of the dropped vertical lines. The node Fn may be connected by a curve.

図13はこの実施形態の変形例を示す、図5と同様のフロー・チャート、図14、図15は図13フロー・チャートの処理の説明図である。   FIG. 13 is a flow chart similar to FIG. 5 showing a modification of this embodiment, and FIGS. 14 and 15 are explanatory diagrams of the processing of the flow chart of FIG.

図5と同様にS300において複数個の節点Fの少なくとも一つの節点Fnを曲線または直線で接続して三次元図形を生成し、S302において生成された三次元図形上の点を曲線または直線で接続して第2の三次元図形を生成し、S304において生成された第2の三次元図形に基づいて繊維16aの曲がり率Afを算出するように構成した。   As in FIG. 5, in S300, at least one node Fn of the plurality of nodes F is connected by a curve or a straight line to generate a three-dimensional figure, and the points on the three-dimensional figure generated in S302 are connected by a curve or a straight line. Thus, the second three-dimensional figure is generated, and the bending rate Af of the fiber 16a is calculated based on the second three-dimensional figure generated in S304.

曲がり率Afは、繊維16aの評価長さと複数個の節点Fから選択された節点Fnを接続して生成される三次元図形上の点を接続して生成される第2の三次元図形の評価高さの比に基づいて算出する。評価長さは三次元図形の始点から終点までの長さとして算出する。   The curvature Af is an evaluation of the second 3D graphic generated by connecting the points on the 3D graphic generated by connecting the evaluation length of the fiber 16a and the node Fn selected from the plurality of nodes F. Calculated based on the height ratio. The evaluation length is calculated as the length from the start point to the end point of the three-dimensional figure.

図14に示すように、節点Fnの中から複数個の点Pを選択し、複数個の点Pを直線で接続して生成される第2の三次元図形を生成し、生成された第2の三次元図形の重心Gから、三次元図形の始点と終点とを結んだ直線に下した垂線の長さとして、評価高さhを算出するように構成した。重心Gの座標は複数個の点Pの座標の平均値として算出する。   As shown in FIG. 14, a plurality of points P are selected from the nodes Fn, a second three-dimensional figure generated by connecting the plurality of points P with straight lines is generated, and the generated second From the center of gravity G of the three-dimensional figure, the evaluation height h is calculated as the length of the perpendicular line drawn from the straight line connecting the start point and the end point of the three-dimensional figure. The coordinates of the center of gravity G are calculated as an average value of the coordinates of a plurality of points P.

図15にこの実施形態の変形例を示す。この変形例では、節点Fnを直線で接続して三次元図形を生成し、生成された三次元図形上に複数個の点Pを生成し、複数個の点Pを直線で接続して第2の三次元図形を生成し、生成された第2の三次元図形の重心Gから、三次元図形の始点と終点とを結んだ直線に下した垂線の長さとして、評価高さhを算出するように構成した。重心Gの座標は複数個の点Pの座標の平均値として算出する。節点Fnは曲線で接続しても良い。   FIG. 15 shows a modification of this embodiment. In this modification, the node Fn is connected by a straight line to generate a three-dimensional figure, a plurality of points P are generated on the generated three-dimensional figure, and the plurality of points P are connected by a straight line to form a second one. And the evaluation height h is calculated as the length of the perpendicular line drawn from the center of gravity G of the generated second 3D figure to the straight line connecting the start point and the end point of the 3D figure. It was configured as follows. The coordinates of the center of gravity G are calculated as an average value of the coordinates of a plurality of points P. The node Fn may be connected by a curve.

図16はこの実施形態の変形例を示す、図5と同様のフロー・チャート、図17、図18は図16フロー・チャートの処理の説明図である。この変形例ではさらに第2の三次元図形を点Pおよび三次元図形の始点と終点を頂点とする三角形として生成し、生成された三角形の高さを評価高さhとして算出するように構成した。   FIG. 16 is a flow chart similar to FIG. 5 showing a modification of this embodiment, and FIGS. 17 and 18 are explanatory diagrams of processing of the flow chart of FIG. In this modification, the second three-dimensional figure is further generated as a triangle having the point P and the start and end points of the three-dimensional figure as vertices, and the height of the generated triangle is calculated as the evaluation height h. .

まずS400において節点Fnとそれに隣接する節点Fn-1,Fn+1とをそれぞれ結ぶ線FnFn-1FnFn+1を斜辺(評価長さに相当)とすると共に、隣接する節点Fn-1,Fn+1を結ぶ線Fn-1Fn+1を底辺とする三角形を設定する。尚、図18においては図示の便宜のため、斜辺と底辺をa,b,cで示した。   First, in S400, lines FnFn-1FnFn + 1 connecting the nodes Fn and the adjacent nodes Fn-1, Fn + 1 are set as hypotenuses (corresponding to evaluation lengths), and adjacent nodes Fn-1, Fn + 1 are used. A triangle with the base of line Fn-1Fn + 1 connecting is established. In FIG. 18, the hypotenuse and the base are indicated by a, b, and c for convenience of illustration.

次いでS402において節点Fnあるいはその付近の任意の点から底辺に向けられた垂線の長さh1を算出する。図示例では垂線は節点Fnから底辺に向けられるが、節点Fn付近の任意の点から底辺に向けられても良い。   In step S402, the length h1 of the perpendicular line directed from the node Fn or an arbitrary point in the vicinity thereof to the base is calculated. In the illustrated example, the vertical line is directed from the node Fn to the base, but may be directed from any point near the node Fn to the base.

次いでS404において算出された垂線の長さh1を斜辺FnFn-1FnFn+1の長さの合計値(評価長さ)で除算してアスペクト比As1を算出する。このようにアスペクト比は3個の節点を用いて算出する。   Next, the aspect ratio As1 is calculated by dividing the length h1 of the perpendicular calculated in S404 by the total value (evaluation length) of the lengths of the hypotenuses FnFn-1FnFn + 1. As described above, the aspect ratio is calculated using three nodes.

次いでS406に進み、隣接する節点に対して同様の処理を行う。   In step S406, the same processing is performed on adjacent nodes.

即ち、節点Fn+1とそれに隣接する節点Fn,Fn+2とをそれぞれ結ぶ線Fn+1FnFn+1Fn+2を斜辺とすると共に、隣接する節点Fn,Fn+2を結ぶ線FnFn+2を底辺とする三角形を設定し、節点Fn+1あるいはその付近の任意の点から斜辺に向けられた垂線の長さh2を算出し、算出された垂線の長さh2を斜辺Fn+1FnFn+1Fn+2の長さの合計値で除算してアスペクト比As2を算出し、以降、同様の処理を繊維16aの終端まで繰り返す。   That is, the line Fn + 1FnFn + 1Fn + 2 connecting the node Fn + 1 and the adjacent nodes Fn and Fn + 2 is the hypotenuse, and the line FnFn + 2 connecting the adjacent nodes Fn and Fn + 2 is the base. Is set, and the length h2 of the perpendicular line directed to the hypotenuse from the node Fn + 1 or an arbitrary point near the node Fn + 1 is calculated, and the calculated length h2 of the perpendicular is the hypotenuse Fn + 1FnFn + 1Fn + 2 Then, the aspect ratio As2 is calculated by dividing by the total value of the lengths, and thereafter the same processing is repeated until the end of the fiber 16a.

具体的には、それ以降の節点(あるいは他端)のいずれかをmと総称するとき、節点Fn+2とそれに隣接する/しない節点Fn+1,Fn+mとをそれぞれ結ぶ線Fn+2Fn+1Fn+2Fn+mを斜辺とすると共に、隣接する節点Fn+1,Fn+mを結ぶ線Fn+1Fn+mを底辺とする三角形を設定し、節点Fn+2あるいはその付近の任意の点から斜辺に向けられた垂線の長さhmを算出し、算出された垂線の長さhmを斜辺Fn+2Fn+1Fn+2Fn+mの長さの合計値で除算してアスペクト比Asmを算出する。   Specifically, when any one of the subsequent nodes (or the other end) is collectively referred to as m, a line Fn + 2Fn connecting the node Fn + 2 and the adjacent nodes Fn + 1 and Fn + m which are not adjacent to each other. Set a triangle with + 1Fn + 2Fn + m as the hypotenuse and a line Fn + 1Fn + m connecting adjacent nodes Fn + 1 and Fn + m as the base, and any point in the vicinity of node Fn + 2 To calculate the aspect ratio Asm by dividing the length hm of the perpendicular directed to the hypotenuse from the total length of the hypotenuse Fn + 2Fn + 1Fn + 2Fn + m .

次いでS408に進み、算出されたアスペクト比As2と、節点Fnについて算出されたアスペクト比As1と、それ以降の終端までについて算出されたアスペクト比Asmとに基づいて繊維の曲がり率Afを繊維16aの長さに対して算出する。   Next, the process proceeds to S408, where the fiber bending rate Af is calculated based on the calculated aspect ratio As2, the aspect ratio As1 calculated for the node Fn, and the aspect ratio Asm calculated up to the end of the node Fn. Calculate for the length.

具体的にはアスペクト比As1,As2,..Asmまでの合計値を算出すると共に、その平均値を算出し、その平均値に基づいて繊維16aの曲がり率Afを算出、より具体的にはその平均値を繊維16aの曲がり率Afとする。   Specifically, the aspect ratios As1, As2,. . The total value up to Asm is calculated, the average value is calculated, the bending rate Af of the fiber 16a is calculated based on the average value, and more specifically, the average value is set as the bending rate Af of the fiber 16a.

図17と図18に図16フロー・チャートに示す処理を示す。   17 and 18 show the processing shown in the flowchart of FIG.

図19はこの実施形態の変形例を示す、図16に類似するフロー・チャートである。   FIG. 19 is a flow chart similar to FIG. 16, showing a modification of this embodiment.

この変形例の場合、S500からS506までに図16と同様の処理を行ってS508に進み、アスペクト比As1,As2,..Asmを算出すると共に、その中の最大値を選択し、その最大値に基づいて繊維16aの曲がり率Afを算出、より具体的にはその最大値を繊維16aの曲がり率とするようにした。   In the case of this modification, the same processing as in FIG. 16 is performed from S500 to S506, and the processing proceeds to S508, where the aspect ratios As1, As2,. . Asm was calculated, the maximum value among them was selected, the bending rate Af of the fiber 16a was calculated based on the maximum value, and more specifically, the maximum value was used as the bending rate of the fiber 16a.

図20はこの実施形態の変形例を示す、図16に類似するフロー・チャート、図21は図20フロー・チャートの処理の説明図である。   FIG. 20 is a flowchart similar to FIG. 16 showing a modification of this embodiment, and FIG. 21 is an explanatory diagram of the processing of the flowchart of FIG.

以下説明すると、S600において節点Fnとそれに隣接する節点Fn-1,Fn+1とをそれぞれ結ぶ線FnFn-1FnFn+1上のいずれかの任意の点Pと隣接する節点Fn-1,Fn+2を結ぶ線PFn-1PFn+2を斜辺とすると共に、節点Fn-1,Fn+2を結ぶ線Fn-1Fn+2を底辺とする三角形を設定する。図示例では任意の点Pは中点とするが、それに限られるものではない。   In the following description, in S600, the node Fn-1, Fn + 2 adjacent to any point P on the line FnFn-1 FnFn + 1 connecting the node Fn and the adjacent nodes Fn-1, Fn + 1 respectively. A triangle is set with the line PFn-1 PFn + 2 connecting the lines Fn-1 and Fn + 2 connecting the nodes Fn-1 and Fn + 2 as the base. In the illustrated example, the arbitrary point P is a middle point, but is not limited thereto.

次いでS602において任意の点Pから底辺に向けられた垂線の長さh1を図18と同様な手法で算出する。   Next, in S602, the length h1 of the perpendicular line directed from the arbitrary point P to the bottom is calculated by the same method as in FIG.

次いでS604において算出された垂線の長さh1を斜辺PFn-1PFn+2の長さの合計値で除算してアスペクト比As1を算出する。即ち、図示例では4個の節点を用いてアスペクト比を算出する。   Next, the aspect ratio As1 is calculated by dividing the length h1 of the perpendicular calculated in S604 by the total value of the lengths of the hypotenuses PFn-1PFn + 2. That is, in the illustrated example, the aspect ratio is calculated using four nodes.

次いでS606において他の節点についてS600からS604に述べた手法でアスペクト比Asnを順次算出する。アスペクト比の算出自体は図18と同様なので、図21においてAsnの図示は省略する。   Next, in step S606, the aspect ratio Asn is sequentially calculated for the other nodes by the method described in steps S600 to S604. Since the aspect ratio calculation itself is the same as in FIG. 18, the illustration of Asn is omitted in FIG.

次いでS608において算出されたアスペクト比As1,Asnに基づいて繊維16aの曲がり率Afを算出する。具体的には算出されたアスペクト比As1,Asnの平均値を算出して曲がり率Afとする、あるいは算出されたアスペクト比As1,Asnの中の最大値を算出して曲がり率とする。   Next, the bending rate Af of the fiber 16a is calculated based on the aspect ratios As1 and Asn calculated in S608. Specifically, the average value of the calculated aspect ratios As1 and Asn is calculated to be the bending rate Af, or the maximum value of the calculated aspect ratios As1 and Asn is calculated to be the bending rate.

尚、この実施形態の変形例として、図22に示す如く、中点を節点にしない場合、同図(a)(b)に示すように周囲の点から内部補間(平均)した値を節点Fn+1としても良い。   As a modification of this embodiment, as shown in FIG. 22, when the midpoint is not a node, a value obtained by internal interpolation (average) from surrounding points as shown in FIGS. It is good also as +1.

図3フロー・チャートのS102の処理においてそれぞれの繊維16aの曲がり率Afを評価長さに対して算出するが、図23に図4に示す4種の繊維についての算出結果を示す。数値で算出することにより、4種の例でいえば、繊維01(図3)は4mmのところにのみ曲がり(うねり)があり、全体として細かい曲がり(うねり)があると評価でき、繊維02は全体的に曲がり(うねり)なしと評価することができる。   In the process of S102 in the flow chart of FIG. 3, the bending rate Af of each fiber 16a is calculated with respect to the evaluation length. FIG. 23 shows the calculation results for the four types of fibers shown in FIG. By calculating numerically, it can be evaluated that the fiber 01 (FIG. 3) has a bend (swell) only at 4 mm and a fine bend (swell) as a whole. Overall, it can be evaluated that there is no bend.

また、繊維03は長さが短い領域も長い領域も曲がり率(アスペクト比)が大きく、細かい曲がり(うねり)だけではなく、大きな曲がりも存在する繊維と評価でき、繊維04は長さが長い領域では曲がり率が大きいことから、細かい曲がり(うねり)はなく、大きな曲がりが存在する繊維と評価することができる。また、繊維16aの曲がりを適正に判定することが可能となることで、成形条件を変更して製品30の強度または弾性率(剛性)を上げることも可能となる。   In addition, the fiber 03 has a large curvature ratio (aspect ratio) in both a short region and a long region, and can be evaluated as a fiber having not only a fine bend (swell) but also a large bend, and the fiber 04 is a long region. Then, since the bending rate is large, there is no fine bending (undulation), and it can be evaluated as a fiber having a large bending. Moreover, since it becomes possible to determine appropriately the bending of the fiber 16a, it is also possible to change the molding conditions and increase the strength or elastic modulus (rigidity) of the product 30.

図3フロー・チャートの説明に戻ると、次いでS104に進み、算出された評価長さに対する曲がり率Afに基づいてそれぞれの繊維16aの曲がりを評価する。   Returning to the description of the flow chart of FIG. 3, the process proceeds to S104, and the bending of each fiber 16a is evaluated based on the bending rate Af with respect to the calculated evaluation length.

図24は図3フロー・チャートのS104の処理のサブ・ルーチン・フロー・チャートである。   FIG. 24 is a sub-routine flowchart of the process of S104 in the flowchart of FIG.

以下説明すると、S700において算出された繊維16aの評価長さに対する曲がり率Af(アスペクト比As)を所定のしきい値TH1と比較し、曲がり率Afがしきい値未満か否か判断し、肯定されるときはS702に進んでその評価長さに対する繊維16aの曲がりなしと判断する。   Explaining below, the bending rate Af (aspect ratio As) with respect to the evaluation length of the fiber 16a calculated in S700 is compared with a predetermined threshold value TH1, and it is determined whether or not the bending rate Af is less than the threshold value. If so, the process proceeds to S702 and it is determined that the fiber 16a is not bent with respect to the evaluation length.

一方、S700で否定されるときはS704に進んでその評価長さに対する繊維16aの曲がりありと判断し、S706に進んで曲がり率Afを適宜設定する第2、第3、・・・の所定のしきい値TH2,TH3,・・・と比較して繊維16を所定の範囲ごとに分類する。   On the other hand, when the result in S700 is negative, the process proceeds to S704, where it is determined that the fiber 16a is bent with respect to the evaluation length, and the process proceeds to S706 to appropriately set the bending rate Af. Compared with the threshold values TH2, TH3,..., The fibers 16 are classified into predetermined ranges.

図2フロー・チャートのS44の処理において成形品の構造強度または弾性率(剛性)を評価するが、具体的には、得られたそれぞれの繊維16aの曲がり評価結果に基づいて選択された繊維群の強度または弾性率(剛性)を評価し、その強度または弾性率(剛性)評価結果に基づいて複数本の繊維16aの全体の強度または弾性率(剛性)を評価する。   The structural strength or elastic modulus (rigidity) of the molded product is evaluated in the processing of S44 in the flow chart of FIG. 2, and specifically, the fiber group selected based on the bending evaluation result of each obtained fiber 16a. The strength or elastic modulus (rigidity) of the plurality of fibers 16a is evaluated based on the strength or elastic modulus (rigidity) evaluation result.

図25は図2フロー・チャートのS44の処理のサブ・ルーチン・フロー・チャートである。   FIG. 25 is a sub-routine flowchart of the process in S44 of the flowchart of FIG.

以下説明すると、S800において所定の範囲(曲がりなし、または曲がりありの所定の範囲)に属する繊維16aの割合(繊維割合)Rを算出する。即ち、選択された繊維群を構成する繊維16aの総本数に対する所定の範囲に属する繊維16aの本数の割合を算出する。   In the following, in S800, the ratio (fiber ratio) R of the fibers 16a belonging to a predetermined range (no bending or a predetermined range with bending) is calculated. That is, the ratio of the number of fibers 16a belonging to a predetermined range to the total number of fibers 16a constituting the selected fiber group is calculated.

次いでS802に進み、算出された繊維割合Rに基づいて選択された繊維群の強度または弾性率(剛性)評価値を算出する。具体的には、S800において1つの繊維群に対して所定の範囲の数だけ算出された繊維割合Rのうち、特に成形品の構造強度または弾性率(剛性)低下に寄与する範囲の繊維割合Rを強度または弾性率(剛性)評価値として算出する。   Next, in S802, the strength or elastic modulus (rigidity) evaluation value of the fiber group selected based on the calculated fiber ratio R is calculated. Specifically, among the fiber ratios R calculated for a single fiber group in a predetermined range in S800, the fiber ratio R in a range particularly contributing to a decrease in the structural strength or elastic modulus (rigidity) of the molded product. Is calculated as an evaluation value of strength or elastic modulus (rigidity).

より具体的には、図24フロー・チャートのS702において曲がりありと判断される所定のしきい値TH1以上のすべての範囲(TH1〜TH2,TH2〜TH3,・・・)の繊維割合Rの合計値を強度または弾性率(剛性)評価値として算出する。   More specifically, the sum of the fiber ratios R in all ranges (TH1 to TH2, TH2 to TH3,...) Greater than or equal to the predetermined threshold value TH1 that is determined to be bent in S702 of the flow chart of FIG. The value is calculated as an evaluation value of strength or elastic modulus (rigidity).

次いでS804に進み、算出された繊維群の強度または弾性率(剛性)評価値から予め設定された特性(強度または弾性率(剛性)算出マップ)に基づいて選択された繊維群の強度または弾性率(剛性)を算出する。即ち、図26に示すような強度または弾性率(剛性)算出マップを参照して強度または弾性率(剛性)評価値に対応する強度または弾性率(剛性)を算出し、S806に進んで要求強度または弾性率(剛性)を満たしているか否か判断する。強度または弾性率(剛性)評価マップは試験などにより事前に作成する。   Next, the process proceeds to S804, and the strength or elastic modulus of the fiber group selected based on the characteristic (strength or elastic modulus (rigidity) calculation map) set in advance from the calculated strength or elastic modulus (rigidity) evaluation value of the fiber group. (Rigidity) is calculated. That is, the strength or elastic modulus (rigidity) corresponding to the strength or elastic modulus (rigidity) evaluation value is calculated with reference to the strength or elastic modulus (rigidity) calculation map as shown in FIG. Alternatively, it is determined whether or not the elastic modulus (rigidity) is satisfied. A strength or elastic modulus (stiffness) evaluation map is created in advance by a test or the like.

S806で肯定されるときはS808に進んでまたは弾性率(剛性)十分と判断する一方、否定されるときはS810に進んで強度または弾性率(剛性)不足と判断する。強度または弾性率(剛性)不足の場合は製品30の強度または弾性率(剛性)が向上するように成形条件を変更する。   When the result in S806 is affirmative, the process proceeds to S808 or it is determined that the elastic modulus (rigidity) is sufficient, while when the result is negative, the process proceeds to S810 and it is determined that the strength or elastic modulus (rigidity) is insufficient. When the strength or elastic modulus (rigidity) is insufficient, the molding conditions are changed so that the strength or elastic modulus (rigidity) of the product 30 is improved.

図27は、成形条件変更の前後において強度または弾性率(剛性)評価値の算出結果を比較した例を示す、繊維16aの曲がり率Afの所定の範囲ごとに繊維割合Rを表示したヒストグラムである。   FIG. 27 is a histogram displaying the fiber ratio R for each predetermined range of the bending rate Af of the fiber 16a, showing an example of comparison of calculation results of strength or elastic modulus (rigidity) evaluation values before and after the molding condition change. .

繊維16aの曲がり率Afが0〜TH1の範囲(曲がりなし)の繊維割合Rは約85%から約97%に増加した一方、TH1以上の範囲(曲がりあり)の繊維割合Rはいずれも減少、即ち、TH1以上の範囲の繊維割合Rの合計値は約15%から約3%に減少し、選択された繊維群(成形品において強度または弾性率(剛性)を評価する範囲に相当する繊維群)全体として曲がりが抑制され、成形品の強度または弾性率(剛性)が向上していることを数値的に確認することができる。   While the fiber ratio R of the fiber 16a in the bending range Af of 0 to TH1 (no bending) increased from about 85% to about 97%, the fiber ratio R in the range of TH1 or more (with bending) decreased. That is, the total value of the fiber ratio R in the range of TH1 or more is reduced from about 15% to about 3%, and the selected fiber group (fiber group corresponding to the range for evaluating the strength or elastic modulus (rigidity) in the molded product). ) It is possible to numerically confirm that bending is suppressed as a whole and the strength or elastic modulus (rigidity) of the molded product is improved.

図28はこの実施形態の変形例を示す、図3と同様のフロー・チャートである。   FIG. 28 is a flow chart similar to FIG. 3, showing a modification of this embodiment.

この変形例の場合、S900からS902までに図3と同様の処理を行ってS904に進み、算出された繊維16aの評価長さに対する曲がり率Afをそれに応じた表示色に関連付けて繊維16aの曲がりを評価する。   In the case of this modification, the same processing as in FIG. 3 is performed from S900 to S902, and the process proceeds to S904, where the bending rate Af with respect to the calculated evaluation length of the fiber 16a is associated with the display color corresponding thereto and the bending of the fiber 16a is performed. To evaluate.

図2フロー・チャートのS44の処理において成形品の構造強度または弾性率(剛性)を評価するが、具体的には、得られたそれぞれの繊維16aの曲がり評価結果に基づいて選択された繊維群の強度または弾性率(剛性)を評価し、その強度または弾性率(剛性)評価結果に基づいて複数本の繊維16aの全体の強度または弾性率(剛性)を評価する。   The structural strength or elastic modulus (rigidity) of the molded product is evaluated in the processing of S44 in the flow chart of FIG. 2, and specifically, the fiber group selected based on the bending evaluation result of each obtained fiber 16a. The strength or elastic modulus (rigidity) of the plurality of fibers 16a is evaluated based on the strength or elastic modulus (rigidity) evaluation result.

図29、図30はそれぞれ成形条件の変更前、変更後における、この実施形態に係る解析結果の一例を示す写真である。図29と図30は複数本の繊維16aの全体、即ち、複数本の繊維16aが混入された樹脂16を金型14内で成形するときの成形品の全体像である。また、図中に繊維群、即ち、成形品において強度または弾性率(剛性)を評価する範囲に相当する繊維群の例を符号A,B,Cで示す。   29 and 30 are photographs showing examples of analysis results according to this embodiment before and after changing the molding conditions, respectively. FIG. 29 and FIG. 30 are the entire images of the molded product when the whole of the plurality of fibers 16 a, that is, the resin 16 mixed with the plurality of fibers 16 a is molded in the mold 14. Further, in the figure, reference numerals A, B, and C denote examples of fiber groups, that is, fiber groups corresponding to a range for evaluating strength or elastic modulus (rigidity) in a molded product.

繊維16aの曲がりの有無や大小を、表示色を用いてコンター表示しているため、成形品における繊維曲がりの深刻な箇所を容易に把握することができ、成形条件変更の前後で選択された繊維群全体として曲がりが抑制され、成形品の強度または弾性率(剛性)が向上していることを視覚的に確認することができる。   Since the presence / absence or size of the fiber 16a is contoured using the display color, it is possible to easily grasp the severe part of the fiber bending in the molded product, and the fiber selected before and after the molding condition change It is possible to visually confirm that the bending of the entire group is suppressed and the strength or elastic modulus (rigidity) of the molded product is improved.

上記した如く、この発明の実施形態およびその変形例にあっては、連続繊維および不連続な長繊維16aが複数本混入された樹脂16を金型14内で所定の成形条件で成形するときの前記繊維16aの曲がり挙動をコンピュータ12に格納されたシミュレーションプログラム20を介して解析する樹脂挙動解析装置10において、前記シミュレーションプログラム20が、前記複数本の繊維16aの複数個の節点Fを少なくとも含む解析条件が入力されるとき、前記複数本の繊維16aのうちの少なくとも一部の繊維群を選択し(S100)、前記入力された解析条件に基づいて前記選択された繊維群を構成するそれぞれの繊維16aの曲がり率Afを算出し(S102)、前記算出されたそれぞれの繊維16aの曲がり率Afに基づいて前記それぞれの繊維16aの曲がりを評価し(S104)、前記得られたそれぞれの繊維16aの曲がり評価結果に基づいて前記選択された繊維群の強度または弾性率(剛性)を評価するステップ(S44)からなる如く構成したので、それぞれの繊維16aの曲がり率Afに基づいて選択された繊維群、即ち、成形品において強度または弾性率(剛性)を評価する範囲の強度または弾性率(剛性)を定量的に評価することができる。   As described above, in the embodiment of the present invention and its modification, the resin 16 in which a plurality of continuous fibers and discontinuous long fibers 16a are mixed is molded in the mold 14 under predetermined molding conditions. In the resin behavior analysis apparatus 10 that analyzes the bending behavior of the fiber 16a via a simulation program 20 stored in the computer 12, the simulation program 20 includes an analysis including at least a plurality of nodes F of the plurality of fibers 16a. When a condition is input, at least a part of the plurality of fibers 16a is selected (S100), and each fiber constituting the selected fiber group based on the input analysis condition A bending rate Af of 16a is calculated (S102), and based on the calculated bending rate Af of each fiber 16a From the step (S44) of evaluating the bending of each fiber 16a (S104) and evaluating the strength or elastic modulus (rigidity) of the selected fiber group based on the obtained bending evaluation result of each fiber 16a. Since it is configured as described above, the strength or elastic modulus (rigidity) in a range in which the strength or elastic modulus (rigidity) is evaluated in the fiber group selected based on the bending rate Af of each fiber 16a, that is, the molded product, is quantitative. Can be evaluated.

また、前記シミュレーションプログラム20は、前記得られた繊維群の強度または弾性率(剛性)評価結果に基づいて前記複数本の繊維16aの全体の強度または弾性率(剛性)を評価するステップ(S44)を含む如く構成したので、上記した効果に加え、複数本の繊維16aの全体、即ち、複数本の繊維16aが混入された樹脂16を金型14内で成形するときの成形品の強度または弾性率(剛性)を定量的に評価することができる。   Further, the simulation program 20 evaluates the overall strength or elastic modulus (rigidity) of the plurality of fibers 16a based on the obtained strength or elastic modulus (rigidity) evaluation result of the fiber group (S44). In addition to the effects described above, the strength or elasticity of the molded product when the whole of the plurality of fibers 16a, that is, the resin 16 mixed with the plurality of fibers 16a is molded in the mold 14 is added. The rate (rigidity) can be quantitatively evaluated.

また、前記それぞれの繊維16aの曲がりを評価するステップは、前記算出されたそれぞれの繊維16aの曲がり率Afを所定のしきい値(TH1,TH2,・・・)と比較して所定の範囲(TH1〜TH2,TH2〜TH3,・・・)ごとに分類して評価するステップ(S700からS706)からなり、前記選択された繊維群の強度または弾性率(剛性)を評価するステップは、前記所定の範囲に属する繊維16aの割合(繊維割合)Rを算出し(S800)、前記算出された繊維16aの割合(繊維割合)Rに基づいて前記選択された繊維群の強度または弾性率(剛性)評価値(TH1以上の範囲の繊維割合Rの合計値)を算出するステップ(S802)からなる如く構成したので、上記した効果に加え、成形品において強度または弾性率(剛性)を評価する範囲の強度または弾性率(剛性)を、強度または弾性率(剛性)低下に寄与する曲がり率Af(TH1以上)などの特定の曲がり率Afの繊維割合Rに基づいて定量的に評価することができる。   Further, the step of evaluating the bending of each of the fibers 16a compares the calculated bending rate Af of each of the fibers 16a with a predetermined threshold (TH1, TH2,...) In a predetermined range ( TH1 to TH2, TH2 to TH3,...) And evaluation (S700 to S706), and the step of evaluating the strength or elastic modulus (rigidity) of the selected fiber group includes the predetermined step The ratio (fiber ratio) R of the fibers 16a belonging to the range is calculated (S800), and the strength or elastic modulus (rigidity) of the selected fiber group based on the calculated ratio (fiber ratio) R of the fibers 16a. Since the evaluation value (the total value of the fiber ratio R in the range of TH1 or more) is calculated (S802), in addition to the above-described effects, the molded product has strength or elasticity. Based on the fiber ratio R of a specific bending rate Af such as a bending rate Af (TH1 or more) that contributes to a decrease in strength or elastic modulus (rigidity) It can be evaluated quantitatively.

また、前記選択された繊維群の強度または弾性率(剛性)を評価するステップは、前記所定の範囲(TH1〜TH2,TH2〜TH3,・・・)に属する繊維16aの割合(繊維割合)Rを前記所定の範囲ごとに表示して評価するステップ(S800、図27)からなる如く構成したので、上記した効果に加え、成形品において強度を評価する範囲の強度を数値的に表示して評価することができ、成形条件を変更した結果を比較して評価することが容易になる。   In the step of evaluating the strength or elastic modulus (rigidity) of the selected fiber group, the ratio (fiber ratio) R of the fibers 16a belonging to the predetermined range (TH1 to TH2, TH2 to TH3,...) In addition to the above-described effects, the strength of the range in which the strength is evaluated in the molded product is numerically displayed and evaluated. This makes it easy to compare and evaluate the results of changing the molding conditions.

また、前記選択された繊維群の強度または弾性率(剛性)を評価するステップは、前記算出された繊維群の強度または弾性率(剛性)評価値(TH1以上の範囲の繊維割合Rの合計値)から予め設定された特性(強度または弾性率(剛性)算出マップ)に基づいて前記選択された繊維群の強度または弾性率(剛性)を算出して評価するステップ(S804からS810)からなる如く構成したので、上記した効果に加え、成形品において強度または弾性率(剛性)を評価する範囲の強度または弾性率(剛性)を精度良く算出することができ、要求強度または弾性率(剛性)と比較して強度または弾性率(剛性)が十分か確認することができる。   Further, the step of evaluating the strength or elastic modulus (rigidity) of the selected fiber group includes the calculated strength or elastic modulus (rigidity) evaluation value of the fiber group (the total value of the fiber ratio R in the range of TH1 or more). ) To calculate and evaluate the strength or elastic modulus (rigidity) of the selected fiber group based on a preset characteristic (strength or elastic modulus (rigidity) calculation map) (S804 to S810). In addition to the effects described above, the strength or elastic modulus (rigidity) within the range for evaluating the strength or elastic modulus (rigidity) of the molded product can be accurately calculated, and the required strength or elastic modulus (rigidity) In comparison, it can be confirmed whether the strength or elastic modulus (rigidity) is sufficient.

また、連続繊維および不連続な長繊維16aが複数本混入された樹脂16を金型14内で所定の成形条件で成形するときの前記繊維16aの曲がり挙動をコンピュータ12に格納されたシミュレーションプログラム20を介して解析する樹脂挙動解析装置10において、前記シミュレーションプログラム20が、前記複数本の繊維16aの複数個の節点Fを少なくとも含む解析条件が入力されるとき、前記複数本の繊維16aのうちの少なくとも一部の繊維群を選択し(S900)、前記入力された解析条件に基づいて前記選択された繊維群を構成するそれぞれの繊維16aの曲がり率Afを算出し(S902)、前記算出されたそれぞれの繊維16aの曲がり率Afを表示色に関連付けて前記それぞれの繊維16aの曲がりを評価し(S904、図29、図30)、前記得られたそれぞれの繊維16aの曲がり評価結果に基づいて前記選択された繊維群の強度または弾性率(剛性)を評価するステップ(S44)からなる如く構成したので、それぞれの繊維16aの曲がり率Afに基づいて成形品において強度または弾性率(剛性)を評価する範囲の強度または弾性率(剛性)を視覚的に表示して評価することができ、成形条件を変更した結果を比較して評価することが容易になる。   A simulation program 20 stored in the computer 12 shows the bending behavior of the fiber 16a when the resin 16 in which a plurality of continuous fibers and discontinuous long fibers 16a are mixed is molded in the mold 14 under predetermined molding conditions. In the resin behavior analysis apparatus 10 that performs analysis through the simulation program 20, when an analysis condition including at least a plurality of nodes F of the plurality of fibers 16a is input, of the plurality of fibers 16a. At least some of the fiber groups are selected (S900), and the bending rate Af of each of the fibers 16a constituting the selected fiber group is calculated based on the input analysis conditions (S902), and the calculated The bending rate Af of each fiber 16a is associated with the display color to evaluate the bending of each fiber 16a (S9). 4, FIG. 29, FIG. 30), and a step (S44) of evaluating the strength or elastic modulus (rigidity) of the selected fiber group based on the obtained bending evaluation result of each fiber 16a. Therefore, it is possible to visually display and evaluate the strength or elastic modulus (rigidity) in the range in which the strength or elastic modulus (rigidity) is evaluated in the molded product based on the bending rate Af of each fiber 16a. It becomes easy to compare and evaluate the results of changing the.

また、前記シミュレーションプログラム20は、前記得られた繊維群の強度または弾性率(剛性)評価結果に基づいて前記複数本の繊維16aの全体の強度または弾性率(剛性)を評価するステップ(S44)を含む如く構成したので、上記した効果に加え、成形品の強度または弾性率(剛性)を定量的に評価することができる。   Further, the simulation program 20 evaluates the overall strength or elastic modulus (rigidity) of the plurality of fibers 16a based on the obtained strength or elastic modulus (rigidity) evaluation result of the fiber group (S44). In addition to the above-described effects, the strength or elastic modulus (rigidity) of the molded product can be quantitatively evaluated.

また、前記それぞれの繊維16aの曲がり率Afを算出するステップは、前記複数個の節点Fの中の少なくとも一つの節点Fnから得られる評価長さに対する前記成形条件から予測される前記それぞれの繊維16aの曲がり率Afを算出するステップ(S102,S902)からなり、前記それぞれの繊維16aの曲がりを評価するステップは、前記算出された評価長さに対する曲がり率Afに基づいて前記それぞれの繊維16aの曲がりを評価するステップ(S104,S904)からなる如く構成したので、上記した効果に加え、評価長さに対する曲がり率Afに基づいてそれぞれの繊維16aの曲がりを評価することで細かいうねり(小さい曲がり)は小さい評価長さで、大きい曲がりは大きい評価長さで評価することも可能となり、それによって曲がりの形態ごとに定量的に評価することができ、細かいうねりと大きい曲がりで構成される複合的な曲がりも定量的に評価することができる。また、それぞれの繊維16aの曲がりを適正に判定することが可能となることで、それぞれの繊維16aの曲がりが少なくなる条件を見つけることが容易となり、成形条件を変更して製品30の強度または弾性率(剛性)を上げることも可能となる。   The step of calculating the bending rate Af of each of the fibers 16a includes the step of calculating each of the fibers 16a predicted from the molding condition with respect to the evaluation length obtained from at least one node Fn among the plurality of nodes F. The step (S102, S902) of calculating the bending rate Af of each of the fibers 16a is evaluated, and the step of evaluating the bending of each of the fibers 16a is based on the bending rate Af with respect to the calculated evaluation length. In addition to the effects described above, the fine undulations (small bends) can be obtained by evaluating the bends of the fibers 16a based on the bend rate Af with respect to the evaluation length in addition to the above-described effects (S104, S904). With a small evaluation length, it becomes possible to evaluate a large bend with a large evaluation length. Can be quantitatively assessed for each form of bending by Les, complex bends composed bend as large as fine undulations can also be quantitatively evaluated. Moreover, since it becomes possible to determine appropriately the bending of each fiber 16a, it becomes easy to find the conditions under which the bending of each fiber 16a is reduced, and the strength or elasticity of the product 30 is changed by changing the molding conditions. It is also possible to increase the rate (rigidity).

また、前記それぞれの繊維16aの曲がり率Afを算出するステップは、前記複数個の節点Fの少なくとも一つの節点Fnを曲線または直線で接続して生成される三次元図形に基づいて前記それぞれの繊維16aの曲がり率Afを算出するステップ(S200,S202)からなる如く構成したので、上記した効果に加え、数値としての繊維16aの曲がり率Afを簡易かつ精度よく算出することができる。   The step of calculating the bending rate Af of each of the fibers 16a includes the step of calculating each of the fibers based on a three-dimensional figure generated by connecting at least one node Fn of the plurality of nodes F with a curve or a straight line. Since the configuration includes the steps (S200, S202) of calculating the bending rate Af of 16a, the bending rate Af of the fiber 16a as a numerical value can be calculated easily and accurately in addition to the above-described effects.

また、前記それぞれの繊維16aの曲がり率Afを算出するステップは、前記複数個の節点Fの少なくとも一つの節点Fnを曲線または直線で接続して生成される三次元図形上の点を曲線または直線で接続して生成される第2の三次元図形に基づいて前記それぞれの繊維16aの曲がり率Afを算出するステップ(S300からS304)からなる如く構成したので、上記した効果に加え、数値としての繊維16aの曲がり率Afを一層簡易に算出することができる。   The step of calculating the bending rate Af of each of the fibers 16a may be obtained by connecting a point on a three-dimensional figure generated by connecting at least one node Fn of the plurality of nodes F with a curve or a line. In addition to the above-described effects, the numerical value as the numerical value is calculated by the step (S300 to S304) of calculating the bending rate Af of each of the fibers 16a based on the second three-dimensional figure generated by connecting with The bending rate Af of the fiber 16a can be calculated more easily.

また、前記第2の三次元図形は三角形からなる如く構成したので、上記した効果に加え、数値としての繊維16aの曲がり率Afをより一層簡易に算出することができる。   Further, since the second three-dimensional figure is constituted by a triangle, in addition to the above-described effects, the bending rate Af of the fiber 16a as a numerical value can be calculated more easily.

また、前記それぞれの繊維16aの曲がり率Afを算出するステップは、前記節点Fnとそれに隣接する節点Fn-1,Fn+1とをそれぞれ結ぶ線FnFn-1FnFn+1を斜辺とすると共に、隣接する節点Fn-1,Fn+1を結ぶ線Fn-1Fn+1を底辺とする前記三角形を設定し(S400,S500)、前記節点Fnあるいはその付近の任意の点から前記底辺に向けられた垂線の長さh1を算出し(S402,S502)、前記算出された垂線の長さh1を前記斜辺FnFn-1FnFn+1の長さの合計値で除算してアスペクト比As1を算出し(S404,S504)、前記算出されたアスペクト比As1に基づいて前記それぞれの繊維の曲がり率Afを算出するステップ(S408,S508)からなる如く構成したので、上記した効果に加え、数値としての繊維16aの曲がり率Afを簡易に算出することができる。   Further, the step of calculating the bending rate Af of each of the fibers 16a uses the line FnFn-1FnFn + 1 connecting the node Fn and the adjacent nodes Fn-1, Fn + 1 as hypotenuses and is adjacent thereto. The triangle having the base Fn-1Fn + 1 connecting the nodes Fn-1 and Fn + 1 is set (S400, S500), and a perpendicular line directed to the base from the node Fn or any nearby point is set. A length h1 is calculated (S402, S502), and an aspect ratio As1 is calculated by dividing the calculated perpendicular length h1 by the total length of the hypotenuses FnFn-1FnFn + 1 (S404, S504). In addition to the above-described effects, the fiber as a numerical value is constituted by the steps (S408, S508) of calculating the bending rate Af of each of the fibers based on the calculated aspect ratio As1. Bending ratio Af of 6a can be calculated easily.

また、前記それぞれの繊維16aの曲がり率Afを算出するステップは、前記節点Fn+1とそれに隣接する節点Fn,Fn+2とをそれぞれ結ぶ線Fn+1FnFn+1Fn+2を斜辺とすると共に、隣接する節点Fn,Fn+2を結ぶ線FnFn+2を底辺とする前記三角形を設定し(S406,S506)、前記節点Fn+1あるいはその付近の任意の点から前記底辺に向けられた垂線の長さh2を算出し、前記算出された垂線の長さh2を前記斜辺Fn+1FnFn+1Fn+2の長さの合計値で除算してアスペクト比As2を算出し(S406,S506)、前記算出されたアスペクト比As2と前記節点Fnについて算出されたアスペクト比As1とに基づいて前記それぞれの繊維の曲がり率Afを算出するステップ(S408,S508)からなる如く構成したので、上記した効果に加え、数値としての繊維16aの曲がり率Afを簡易かつ精度良く算出することができる。   The step of calculating the bending rate Af of each of the fibers 16a includes a line Fn + 1FnFn + 1Fn + 2 connecting the node Fn + 1 and the adjacent nodes Fn, Fn + 2 as hypotenuses, respectively. The triangle having a base FnFn + 2 connecting adjacent nodes Fn and Fn + 2 is set (S406, S506), and a perpendicular line directed to the base from the node Fn + 1 or an arbitrary point in the vicinity thereof is set. The length h2 is calculated, the aspect ratio As2 is calculated by dividing the calculated length h2 of the perpendicular by the total length of the hypotenuse Fn + 1FnFn + 1Fn + 2 (S406, S506), and the calculation The above-described aspect ratio As2 and the aspect ratio As1 calculated for the node Fn are used to calculate the bending rate Af of each of the fibers (S408, S508). In addition to the effects, the bending ratio Af of the fibers 16a of the numeric can be calculated easily and accurately.

また、前記それぞれの繊維16aの曲がり率Afを算出するステップは、前記節点Fn+2とそれに隣接する/しない節点Fn+2,Fn+mとをそれぞれ結ぶ線Fn+2Fn+1Fn+2Fn+mを斜辺とすると共に、隣接する節点Fn+1,Fn+mを結ぶ線Fn+1Fn+mを底辺とする前記三角形を設定し、前記節点Fn+2あるいはその付近の任意の点から前記底辺に向けられた垂線の長さhmを算出し、前記算出された垂線の長さhmを前記斜辺Fn+2Fn+1Fn+2Fn+mの長さの合計値で除算してアスペクト比Asmを算出し、前記算出されたアスペクト比Asmと前記節点Fnと前記節点Fn+1について算出されたアスペクト比As1,As2とに基づいて前記それぞれの繊維の曲がり率Afを算出するステップ(S400からS408,S500からS508)からなる如く構成したので、上記した効果に加え、数値としての繊維16aの曲がり率Afを簡易かつ一層精度良く算出することができる。   Further, the step of calculating the bending rate Af of each of the fibers 16a includes a line Fn + 2Fn + 1Fn + 2Fn + m connecting the node Fn + 2 and the adjacent nodes Fn + 2 and Fn + m. Is set to the hypotenuse, and the triangle having the base Fn + 1Fn + m connecting adjacent nodes Fn + 1 and Fn + m is set, and the node Fn + 2 or an arbitrary point in the vicinity thereof is set to the base. Calculate the length hm of the directed perpendicular, and divide the calculated length hm of the perpendicular by the total length of the hypotenuse Fn + 2Fn + 1Fn + 2Fn + m to calculate the aspect ratio Asm; Steps of calculating the bending ratios Af of the respective fibers based on the calculated aspect ratio Asm, the aspect ratios As1 and As2 calculated for the node Fn and the node Fn + 1 (S400 to S408, S500 to S508) ) , In addition to the effects mentioned above, the bending ratio Af of the fibers 16a of the numeric can be easily and calculated more accurately.

また、前記それぞれの繊維16aの曲がり率Afを算出するステップは、前記節点Fnとそれに隣接する節点Fn-1,Fn+1とをそれぞれ結ぶ線FnFn-1FnFn+1上のいずれかの任意の点Pと隣接する節点Fn-1,Fn+2を結ぶ線PFn-1PFn+2を斜辺とすると共に、節点Fn-1,Fn+2を結ぶ線Fn-1Fn+2を底辺とする前記三角形を設定し(S600)、前記任意の点Pから前記底辺に向けられた垂線の長さh1を算出し(S602)、前記算出された垂線の長さh1を前記斜辺PFn-1PFn+2の長さの合計値で除算してアスペクト比As1を算出し(S604)、次いで他の節点について前記アスペクト比Asnを算出し(S606)、前記算出されたアスペクト比As1,Asnに基づいて前記それぞれの繊維の曲がり率Afを算出するステップ(S608)からなる如く構成したので、上記した効果に加え、数値としての繊維16aの曲がり率Afを簡易かつ一層精度良く算出することができる。   The step of calculating the bending rate Af of each of the fibers 16a includes any point on the line FnFn-1FnFn + 1 connecting the node Fn and the adjacent nodes Fn-1, Fn + 1. Set the triangle with the line PFn-1 PFn + 2 connecting the nodes Fn-1 and Fn + 2 adjacent to P as the hypotenuse and the line Fn-1Fn + 2 connecting the nodes Fn-1 and Fn + 2 as the base (S600), the length h1 of the perpendicular line directed from the arbitrary point P to the base side is calculated (S602), and the calculated length h1 of the perpendicular line is set to the length of the oblique side PFn-1PFn + 2. The aspect ratio As1 is calculated by dividing by the total value (S604), then the aspect ratio Asn is calculated for other nodes (S606), and the bending of the respective fibers is calculated based on the calculated aspect ratios As1 and Asn. Calculating the rate Af (S60) Since it is configured as consisting), in addition to the effects mentioned above, and simplify the bending ratio Af of the fibers 16a of the numeric can be calculated more accurately.

また、前記それぞれの繊維16aの曲がり率Afを算出するステップは、前記算出された値(アスペクト比As1,As2;As1,As2,Asm;As1,Asm)の平均値に基づいて前記それぞれの繊維の曲がり率Afを算出するステップ(S408)からなる如く構成したので、上記した効果に加え、数値としての繊維16aの曲がり率Afを簡易かつ一層精度良く算出することができる。   The step of calculating the bending rate Af of each of the fibers 16a is based on the average value of the calculated values (aspect ratio As1, As2; As1, As2, Asm; As1, Asm). Since the configuration includes the step (S408) of calculating the bending rate Af, in addition to the above-described effects, the bending rate Af of the fiber 16a as a numerical value can be calculated easily and more accurately.

また、前記それぞれの繊維16aの曲がり率Afを算出するステップは、前記算出された値(アスペクト比As1,As2;As1,As2,Asm;As1,Asm)の最大値に基づいて前記それぞれの繊維の曲がり率Afを算出するステップ(S508)からなる如く構成したので、上記した効果に加え、数値としての繊維16aの曲がり率Afを簡易かつ一層精度良く算出することができる。   The step of calculating the bending rate Af of each of the fibers 16a is based on the maximum value of the calculated values (aspect ratio As1, As2; As1, As2, Asm; As1, Asm). Since the configuration includes the step (S508) of calculating the bending rate Af, in addition to the above-described effects, the bending rate Af of the fiber 16a as a numerical value can be calculated easily and more accurately.

また、前記それぞれの繊維16aの曲がり率Afを算出するステップは、前記金型14のキャビティ14cにおける前記繊維16aの三次元空間での座標位置x,y,zに基づいて前記それぞれの繊維16aの曲がり率Afを算出するステップ(S102)からなる如く構成したので、数値としての繊維16aの曲がり率Afを一層簡易に算出することができる。   The step of calculating the bending rate Af of each of the fibers 16a is based on the coordinate positions x, y, and z of the fibers 16a in the three-dimensional space in the cavity 14c of the mold 14. Since the configuration includes the step (S102) of calculating the bending rate Af, the bending rate Af of the fiber 16a as a numerical value can be calculated more easily.

また、前記評価長さは前記三次元図形の始点から終点までの長さからなる如く構成したので、上記した効果に加え、数値としての繊維の曲がり率Afを一層精度良く算出することができる。   Further, since the evaluation length is configured to be the length from the start point to the end point of the three-dimensional figure, in addition to the effects described above, the fiber bending rate Af as a numerical value can be calculated with higher accuracy.

また、前記シミュレーションプログラム20は、分析器を介して得られる前記複数本の繊維16aの複数個の節点Fを少なくとも含むデータが入力されるとき、前記複数本の繊維16aのうちの少なくとも一部の繊維群を選択し(S100)、前記入力されたデータに基づいて前記選択された繊維群を構成するそれぞれの繊維16aの曲がり率Afを算出し(S102)、前記算出されたそれぞれの繊維16aの曲がり率Afに基づいて前記それぞれの繊維16aの曲がりを評価し(S104)、前記得られたそれぞれの繊維16aの曲がり評価結果に基づいて前記選択された繊維群の強度または弾性率(剛性)を評価するステップ(S44)からなる如く構成したので、実際の成形品を分析器、例えばX線CTスキャンで分析して得られるデータ(実際の成形品中の繊維の形状、節点)をシミュレーションデータと同様の手法で解析して比較することができる。   The simulation program 20 receives at least a part of the plurality of fibers 16a when data including at least a plurality of nodes F of the plurality of fibers 16a obtained through an analyzer is input. A fiber group is selected (S100), a bending rate Af of each fiber 16a constituting the selected fiber group is calculated based on the input data (S102), and the calculated fiber 16a of each fiber 16a is calculated. The bending of each of the fibers 16a is evaluated based on the bending rate Af (S104), and the strength or elastic modulus (rigidity) of the selected fiber group is calculated based on the bending evaluation result of each of the obtained fibers 16a. Since it comprises so that the evaluation step (S44) may be comprised, it is obtained by analyzing an actual molded product with an analyzer, for example, an X-ray CT scan. Chromatography data (actual shape of the fibers in the molded article, nodes) may be to analyze in the same manner and simulation data comparing.

尚、この実施形態およびその変形例においては、評価長さは一定の値(4mmなど)としたが、不定の値(4mm,6mm,...,20mmなど)としても良い。   In this embodiment and its modification, the evaluation length is a constant value (4 mm, etc.), but may be an indefinite value (4 mm, 6 mm,..., 20 mm, etc.).

また、繊維16aの分割数は10に限られるものではなく、さらに分割長さも不定でも良いことはいうまでもない。また、繊維16aの一端(始端)から始めて他端(終端)まで節点の全てあるいは一部についてアスペクト比Asnを算出したが、繊維16aの途中から始めても良く、節点の一部についてのみ算出するようにしても良い。また、アスペクト比を算出するのに3個あるいは4個の節点を用いたが、それに限られるものではない。また、繊維16aは1本に限られるものではなく、複数本であっても良い。   Needless to say, the number of divisions of the fibers 16a is not limited to 10, and the division length may be indefinite. In addition, the aspect ratio Asn is calculated for all or part of the nodes starting from one end (starting end) to the other end (ending) of the fiber 16a. However, the aspect ratio Asn may be started in the middle of the fiber 16a, and only for a part of the nodes. Anyway. Further, although three or four nodes are used to calculate the aspect ratio, the present invention is not limited to this. Further, the number of fibers 16a is not limited to one, and a plurality of fibers may be used.

また、上記において装置構成も開示したものに限られるものではないことはいうまでもない。   In addition, it goes without saying that the device configuration is not limited to that disclosed above.

10 コンピュータ支援の樹脂挙動解析装置、12 コンピュータ、12a ディスプレイ、12b 入力デバイス、14 金型、14a 上型、14b 下型、14c キャビティ、16 樹脂、16a 繊維、20 シミュレーションプログラム、22 製品モデル、24 金型モデル、26 NC加工機、30 製品   10 Computer Assisted Resin Behavior Analyzer, 12 Computer, 12a Display, 12b Input Device, 14 Mold, 14a Upper Mold, 14b Lower Mold, 14c Cavity, 16 Resin, 16a Fiber, 20 Simulation Program, 22 Product Model, 24 Gold Model, 26 NC processing machines, 30 products

Claims (20)

連続繊維および不連続な長繊維が複数本混入された樹脂を金型内で所定の成形条件で成形するときの前記繊維の曲がり挙動をコンピュータに格納されたシミュレーションプログラムを介して解析する樹脂挙動解析装置において、前記シミュレーションプログラムが、前記複数本の繊維の複数個の節点Fを少なくとも含む解析条件が入力されるとき、前記複数本の繊維のうちの少なくとも一部の繊維群を選択し、前記入力された解析条件に基づいて前記選択された繊維群を構成するそれぞれの繊維の曲がり率Afを算出し、前記算出されたそれぞれの繊維の曲がり率Afに基づいて前記それぞれの繊維の曲がりを評価し、前記得られたそれぞれの繊維の曲がり評価結果に基づいて前記選択された繊維群の強度または弾性率を評価するステップからなることを特徴とするコンピュータ支援の樹脂挙動解析装置。   Resin behavior analysis that analyzes the bending behavior of the fiber mixed with a plurality of continuous fibers and discontinuous long fibers in a mold under predetermined molding conditions through a simulation program stored in a computer In the apparatus, when an analysis condition including at least a plurality of nodes F of the plurality of fibers is input, the simulation program selects at least a part of the plurality of fibers, and the input The bending rate Af of each of the fibers constituting the selected fiber group is calculated based on the analysis condition thus calculated, and the bending of each of the fibers is evaluated based on the calculated bending rate Af of each of the fibers. Evaluating the strength or elastic modulus of the selected fiber group based on the bending evaluation result of each of the obtained fibers. Resin behavior analysis device of a computer assisted characterized by comprising. 前記シミュレーションプログラムは、前記得られた繊維群の強度または弾性率評価結果に基づいて前記複数本の繊維の全体の強度または弾性率を評価するステップを含むことを特徴とする請求項1記載のコンピュータ支援の樹脂挙動解析装置。   2. The computer according to claim 1, wherein the simulation program includes a step of evaluating an overall strength or elastic modulus of the plurality of fibers based on the obtained strength or elastic modulus evaluation result of the fiber group. Supporting resin behavior analyzer. 前記それぞれの繊維の曲がりを評価するステップは、前記算出されたそれぞれの繊維の曲がり率Afを所定のしきい値と比較して所定の範囲ごとに分類して評価するステップからなり、前記選択された繊維群の強度または弾性率を評価するステップは、前記所定の範囲に属する繊維の割合を算出し、前記算出された繊維の割合に基づいて前記選択された繊維群の強度または弾性率評価値を算出するステップからなることを特徴とする請求項1または2記載のコンピュータ支援の樹脂挙動解析装置。   The step of evaluating the bending of each of the fibers includes the step of comparing the calculated bending rate Af of each of the fibers with a predetermined threshold value and classifying the calculated bending rate for each predetermined range, and selecting the selected bending rate. The step of evaluating the strength or elastic modulus of the selected fiber group calculates the ratio of the fibers belonging to the predetermined range, and the strength or elastic modulus evaluation value of the selected fiber group based on the calculated ratio of the fibers The computer-aided resin behavior analysis apparatus according to claim 1 or 2, characterized by comprising the steps of: 前記選択された繊維群の強度または弾性率を評価するステップは、前記所定の範囲に属する繊維の割合を前記所定の範囲ごとに表示して評価するステップからなることを特徴とする請求項3記載のコンピュータ支援の樹脂挙動解析装置。   4. The step of evaluating the strength or elastic modulus of the selected fiber group includes the step of displaying and evaluating the proportion of fibers belonging to the predetermined range for each predetermined range. Computer-aided resin behavior analyzer. 前記選択された繊維群の強度または弾性率を評価するステップは、前記算出された繊維群の強度または弾性率評価値から予め設定された特性に基づいて前記選択された繊維群の強度または弾性率を算出して評価するステップからなることを特徴とする請求項1から4のいずれかに記載のコンピュータ支援の樹脂挙動解析装置。   The step of evaluating the strength or elastic modulus of the selected fiber group includes the step of evaluating the strength or elastic modulus of the selected fiber group based on a property preset from the calculated strength or elastic modulus evaluation value of the fiber group. 5. The computer-aided resin behavior analysis apparatus according to claim 1, further comprising a step of calculating and evaluating the value. 連続繊維および不連続な長繊維が複数本混入された樹脂を金型内で所定の成形条件で成形するときの前記繊維の曲がり挙動をコンピュータに格納されたシミュレーションプログラムを介して解析する樹脂挙動解析装置において、前記シミュレーションプログラムが、前記複数本の繊維の複数個の節点Fを少なくとも含む解析条件が入力されるとき、前記複数本の繊維のうちの少なくとも一部の繊維群を選択し、前記入力された解析条件に基づいて前記選択された繊維群を構成するそれぞれの繊維の曲がり率Afを算出し、前記算出されたそれぞれの繊維の曲がり率Afを表示色に関連付けて前記それぞれの繊維の曲がりを評価し、前記得られたそれぞれの繊維の曲がり評価結果に基づいて前記選択された繊維群の強度または弾性率を評価するステップからなることを特徴とするコンピュータ支援の樹脂挙動解析装置。   Resin behavior analysis that analyzes the bending behavior of the fiber mixed with a plurality of continuous fibers and discontinuous long fibers in a mold under predetermined molding conditions through a simulation program stored in a computer In the apparatus, when an analysis condition including at least a plurality of nodes F of the plurality of fibers is input, the simulation program selects at least a part of the plurality of fibers, and the input The bending rate Af of each of the fibers constituting the selected fiber group is calculated based on the analyzed condition, and the bending rate of each of the fibers is related to the display color by associating the calculated bending rate Af of each of the fibers. And evaluate the strength or elastic modulus of the selected fiber group based on the obtained bending evaluation result of each fiber. Resin behavior analysis device of a computer assisted, characterized in that it consists of steps. 前記シミュレーションプログラムは、前記得られた繊維群の強度または弾性率評価結果に基づいて前記複数本の繊維の全体の強度または弾性率を評価するステップを含むことを特徴とする請求項6記載のコンピュータ支援の樹脂挙動解析装置。   7. The computer according to claim 6, wherein the simulation program includes a step of evaluating an overall strength or elastic modulus of the plurality of fibers based on the obtained strength or elastic modulus evaluation result of the fiber group. Supporting resin behavior analyzer. 前記それぞれの繊維の曲がり率Afを算出するステップは、前記複数個の節点Fの中の少なくとも一つの節点Fnから得られる評価長さに対する前記成形条件から予測される前記それぞれの繊維の曲がり率Afを算出するステップからなり、前記それぞれの繊維の曲がりを評価するステップは、前記算出された評価長さに対する曲がり率Afに基づいて前記それぞれの繊維の曲がりを評価するステップからなることを特徴とする請求項1から7のいずれかに記載のコンピュータ支援の樹脂挙動解析装置。   The step of calculating the bending rate Af of each of the fibers includes the bending rate Af of each of the fibers predicted from the molding condition with respect to an evaluation length obtained from at least one node Fn among the plurality of nodes F. The step of evaluating the bending of each of the fibers comprises the step of evaluating the bending of each of the fibers based on a bending rate Af with respect to the calculated evaluation length. The computer-aided resin behavior analysis apparatus according to claim 1. 前記それぞれの繊維の曲がり率Afを算出するステップは、前記複数個の節点Fの少なくとも一つの節点Fnを曲線または直線で接続して生成される三次元図形に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなることを特徴とする請求項8記載のコンピュータ支援の樹脂挙動解析装置。   The step of calculating the bending rate Af of each of the fibers includes the bending rate of each of the fibers based on a three-dimensional figure generated by connecting at least one node Fn of the plurality of nodes F with a curve or a straight line. 9. The computer-aided resin behavior analysis apparatus according to claim 8, further comprising a step of calculating Af. 前記それぞれの繊維の曲がり率Afを算出するステップは、前記複数個の節点Fの少なくとも一つの節点Fnを曲線または直線で接続して生成される三次元図形上の点を曲線または直線で接続して生成される第2の三次元図形に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなることを特徴とする請求項8記載のコンピュータ支援の樹脂挙動解析装置。   In the step of calculating the bending rate Af of each of the fibers, the points on the three-dimensional figure generated by connecting at least one node Fn of the plurality of nodes F with a curve or a line are connected with the curve or the line. 9. The computer-aided resin behavior analysis apparatus according to claim 8, further comprising a step of calculating a bending rate Af of each of the fibers based on the second three-dimensional figure generated in the above-described manner. 前記第2の三次元図形は三角形からなることを特徴とする請求項10記載のコンピュータ支援の樹脂挙動解析装置。   The computer-aided resin behavior analysis apparatus according to claim 10, wherein the second three-dimensional figure is a triangle. 前記それぞれの繊維の曲がり率Afを算出するステップは、前記節点Fnとそれに隣接する節点Fn-1,Fn+1とをそれぞれ結ぶ線FnFn-1FnFn+1を斜辺とすると共に、隣接する節点Fn-1,Fn+1を結ぶ線Fn-1Fn+1を底辺とする前記三角形を設定し、前記節点Fnあるいはその付近の任意の点から前記底辺に向けられた垂線の長さh1を算出し、前記算出された垂線の長さh1を前記斜辺FnFn-1FnFn+1の長さの合計値で除算してアスペクト比As1を算出し、前記算出されたアスペクト比As1に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなることを特徴とする請求項11記載のコンピュータ支援の樹脂挙動解析装置。   The step of calculating the bending rate Af of each of the fibers includes a line FnFn-1FnFn + 1 connecting the node Fn and the adjacent nodes Fn-1 and Fn + 1, respectively, and an adjacent node Fn- 1. Set the triangle whose base is the line Fn-1Fn + 1 connecting 1, Fn + 1, calculate the length h1 of the perpendicular line directed to the base from the node Fn or an arbitrary point in the vicinity thereof, The aspect ratio As1 is calculated by dividing the calculated length h1 of the perpendicular line by the total length of the hypotenuses FnFn-1FnFn + 1, and the bending rate of each fiber based on the calculated aspect ratio As1. The computer-aided resin behavior analysis apparatus according to claim 11, comprising a step of calculating Af. 前記それぞれの繊維の曲がり率Afを算出するステップは、前記節点Fn+1とそれに隣接する節点Fn,Fn+2とをそれぞれ結ぶ線Fn+1FnFn+1Fn+2を斜辺とすると共に、隣接する節点Fn,Fn+2を結ぶ線FnFn+2を底辺とする前記三角形を設定し、前記節点Fn+1あるいはその付近の任意の点から前記底辺に向けられた垂線の長さh2を算出し、前記算出された垂線の長さh2を前記斜辺Fn+1FnFn+1Fn+2の長さの合計値で除算してアスペクト比As2を算出し、前記算出されたアスペクト比As2と前記節点Fnについて算出されたアスペクト比As1とに基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなることを特徴とする請求項12記載のコンピュータ支援の樹脂挙動解析装置。   The step of calculating the bending rate Af of each of the fibers includes a line Fn + 1FnFn + 1Fn + 2 connecting the node Fn + 1 and the nodes Fn and Fn + 2 adjacent to the node Fn + 1 and an adjacent node. The triangle having the base FnFn + 2 connecting Fn and Fn + 2 is set, and the length h2 of the perpendicular line directed to the base from the node Fn + 1 or any point nearby is calculated. The aspect ratio As2 is calculated by dividing the calculated length h2 of the perpendicular by the total length of the hypotenuse Fn + 1FnFn + 1Fn + 2, and the aspect ratio As2 and the node Fn are calculated. 13. The computer-aided resin behavior analysis apparatus according to claim 12, comprising a step of calculating a bending rate Af of each of the fibers based on the aspect ratio As1. 前記それぞれの繊維の曲がり率Afを算出するステップは、前記節点Fn+2とそれに隣接する/しない節点Fn+2,Fn+mとをそれぞれ結ぶ線Fn+2Fn+1Fn+2Fn+mを斜辺とすると共に、隣接する節点Fn+1,Fn+mを結ぶ線Fn+1Fn+mを底辺とする前記三角形を設定し、前記節点Fn+2あるいはその付近の任意の点から前記底辺に向けられた垂線の長さhmを算出し、前記算出された垂線の長さhmを前記斜辺Fn+2Fn+1Fn+2Fn+mの長さの合計値で除算してアスペクト比Asmを算出し、前記算出されたアスペクト比Asmと前記節点Fnと前記節点Fn+1について算出されたアスペクト比As1,As2とに基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなることを特徴とする請求項13記載のコンピュータ支援の樹脂挙動解析装置。   In the step of calculating the bending rate Af of each of the fibers, a line Fn + 2Fn + 1Fn + 2Fn + m connecting the node Fn + 2 and the nodes Fn + 2 and Fn + m adjacent to / not adjacent to the node Fn + 2 At the same time, the triangle having a base Fn + 1Fn + m connecting adjacent nodes Fn + 1 and Fn + m is set and directed from the node Fn + 2 or an arbitrary point in the vicinity thereof to the base. The vertical length hm is calculated, and the calculated vertical length hm is divided by the total length of the hypotenuses Fn + 2Fn + 1Fn + 2Fn + m to calculate the aspect ratio Asm. 14. The method according to claim 13, further comprising: calculating a bending rate Af of each of the fibers based on the aspect ratio Asm, the aspect ratios As1 and As2 calculated for the node Fn and the node Fn + 1. Computer-aided resin behavior analyzer. 前記それぞれの繊維の曲がり率Afを算出するステップは、前記節点Fnとそれに隣接する節点Fn-1,Fn+1とをそれぞれ結ぶ線FnFn-1FnFn+1上のいずれかの任意の点Pと隣接する節点Fn-1,Fn+2を結ぶ線PFn-1PFn+2を斜辺とすると共に、節点Fn-1,Fn+2を結ぶ線Fn-1Fn+2を底辺とする前記三角形を設定し、前記任意の点Pから前記底辺に向けられた垂線の長さh1を算出し、前記算出された垂線の長さh1を前記斜辺PFn-1PFn+2の長さの合計値で除算してアスペクト比As1を算出し、次いで他の節点について前記アスペクト比Asnを算出し、前記算出されたアスペクト比As1,Asnに基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなることを特徴とする請求項11記載のコンピュータ支援の樹脂挙動解析装置。   The step of calculating the bending rate Af of each of the fibers is adjacent to any point P on the line FnFn-1FnFn + 1 connecting the node Fn and the adjacent nodes Fn-1, Fn + 1. The triangles having the line PFn-1 PFn + 2 connecting the nodes Fn-1 and Fn + 2 to be processed as the hypotenuse and the line Fn-1 Fn + 2 connecting the nodes Fn-1 and Fn + 2 as the base are set. The length h1 of the perpendicular line directed from the arbitrary point P to the base is calculated, and the calculated length h1 of the perpendicular line is divided by the total length of the oblique sides PFn-1PFn + 2 to obtain the aspect ratio As1. Then, the aspect ratio Asn is calculated for other nodes, and the bending ratio Af of each of the fibers is calculated based on the calculated aspect ratios As1 and Asn. 11. Computer-aided resin behavior solution according to 11 Apparatus. 前記それぞれの繊維の曲がり率Afを算出するステップは、前記算出された値の平均値に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなることを特徴とする請求項13から15のいずれかに記載のコンピュータ支援の樹脂挙動解析装置。   16. The step of calculating the bending rate Af of each of the fibers comprises the step of calculating the bending rate Af of each of the fibers based on an average value of the calculated values. The computer-aided resin behavior analysis apparatus according to any one of the above. 前記それぞれの繊維の曲がり率Afを算出するステップは、前記算出された値の最大値に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなることを特徴とする請求項13から15のいずれかに記載のコンピュータ支援の樹脂挙動解析装置。   16. The step of calculating the bending rate Af of each of the fibers comprises the step of calculating the bending rate Af of each of the fibers based on the maximum value of the calculated values. The computer-aided resin behavior analysis apparatus according to any one of the above. 前記それぞれの繊維の曲がり率Afを算出するステップは、前記金型のキャビティにおける前記繊維の三次元空間での座標位置に基づいて前記それぞれの繊維の曲がり率Afを算出するステップからなることを特徴とする請求項1から17のいずれかに記載のコンピュータ支援の樹脂挙動解析装置。   The step of calculating the bending rate Af of each of the fibers includes the step of calculating the bending rate Af of each of the fibers based on a coordinate position of the fiber in a three-dimensional space in the cavity of the mold. A computer-aided resin behavior analysis apparatus according to any one of claims 1 to 17. 前記評価長さは前記三次元図形の始点から終点までの長さからなることを特徴とする請求項8から18のいずれかに記載のコンピュータ支援の樹脂挙動解析装置。   The computer-aided resin behavior analysis apparatus according to claim 8, wherein the evaluation length is a length from a start point to an end point of the three-dimensional figure. 前記シミュレーションプログラムは、分析器を介して得られる前記複数本の繊維の複数個の節点Fを少なくとも含むデータが入力されるとき、前記複数本の繊維のうちの少なくとも一部の繊維群を選択し、前記入力されたデータに基づいて前記選択された繊維群を構成するそれぞれの繊維の曲がり率Afを算出し、前記算出されたそれぞれの繊維の曲がり率Afに基づいて前記それぞれの繊維の曲がりを評価し、前記得られたそれぞれの繊維の曲がり評価結果に基づいて前記選択された繊維群の強度または弾性率を評価するステップを含むことを特徴とする請求項1から19のいずれかに記載のコンピュータ支援の樹脂挙動解析装置。   The simulation program selects at least some of the plurality of fibers when the data including at least the plurality of nodes F of the plurality of fibers obtained through the analyzer is input. The bending rate Af of each of the fibers constituting the selected fiber group is calculated based on the input data, and the bending of each of the fibers is calculated based on the calculated bending rate Af of each of the fibers. 20. The method according to claim 1, further comprising a step of evaluating and evaluating a strength or an elastic modulus of the selected fiber group based on a bending evaluation result of each of the obtained fibers. Computer-aided resin behavior analyzer.
JP2015134616A 2015-04-14 2015-07-03 Computer-aided resin behavior analyzer Expired - Fee Related JP6203787B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2015134616A JP6203787B2 (en) 2015-07-03 2015-07-03 Computer-aided resin behavior analyzer
US15/096,168 US10223481B2 (en) 2015-04-14 2016-04-11 Computer-aided resin behavior analyzer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015134616A JP6203787B2 (en) 2015-07-03 2015-07-03 Computer-aided resin behavior analyzer

Publications (2)

Publication Number Publication Date
JP2017013437A true JP2017013437A (en) 2017-01-19
JP6203787B2 JP6203787B2 (en) 2017-09-27

Family

ID=57827611

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015134616A Expired - Fee Related JP6203787B2 (en) 2015-04-14 2015-07-03 Computer-aided resin behavior analyzer

Country Status (1)

Country Link
JP (1) JP6203787B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6420881B1 (en) * 2017-10-31 2018-11-07 株式会社Jsol Compression molding analysis system, compression molding analysis method, and compression molding analysis program
WO2020111039A1 (en) * 2018-11-29 2020-06-04 東レエンジニアリング株式会社 Movement analysis device for fibers in resin, movement analysis method for fibers in resin, program, and storage medium
WO2021002210A1 (en) * 2019-07-02 2021-01-07 本田技研工業株式会社 Resin behavior analysis device, resin behavior analysis method, and resin behavior analysis program

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05314091A (en) * 1992-05-01 1993-11-26 Toyota Central Res & Dev Lab Inc Method for analyzing motion of grain in fluidized substrate
JP2005250767A (en) * 2004-03-03 2005-09-15 Yokohama Rubber Co Ltd:The Dynamic characteristic simulation method for composite material and dynamic behavior simulation device for composite material
JP2012131220A (en) * 2010-11-29 2012-07-12 Toray Ind Inc Strength analysis device, and strength analysis program
JP2013097521A (en) * 2011-10-31 2013-05-20 Toray Ind Inc Strength analysis method for composite material
JP2013166282A (en) * 2012-02-15 2013-08-29 Japan Polypropylene Corp Method of analyzing flow process of fiber filler-containing thermoplastic resin material
JP2014108612A (en) * 2012-12-04 2014-06-12 Ehime Univ Analysis method in molding carbon fiber reinforced plastic, and analyzer
JP2014226871A (en) * 2013-05-24 2014-12-08 東レエンジニアリング株式会社 Method for simulating behavior of filler, and method for analyzing physical property of composite material
US20160342716A1 (en) * 2015-04-14 2016-11-24 Honda Motor Co., Ltd. Computer-aided resin behavior analyzer
JP2016198997A (en) * 2015-04-14 2016-12-01 本田技研工業株式会社 Computer-aided resin behavior analysis device

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05314091A (en) * 1992-05-01 1993-11-26 Toyota Central Res & Dev Lab Inc Method for analyzing motion of grain in fluidized substrate
JP2005250767A (en) * 2004-03-03 2005-09-15 Yokohama Rubber Co Ltd:The Dynamic characteristic simulation method for composite material and dynamic behavior simulation device for composite material
JP2012131220A (en) * 2010-11-29 2012-07-12 Toray Ind Inc Strength analysis device, and strength analysis program
JP2013097521A (en) * 2011-10-31 2013-05-20 Toray Ind Inc Strength analysis method for composite material
JP2013166282A (en) * 2012-02-15 2013-08-29 Japan Polypropylene Corp Method of analyzing flow process of fiber filler-containing thermoplastic resin material
JP2014108612A (en) * 2012-12-04 2014-06-12 Ehime Univ Analysis method in molding carbon fiber reinforced plastic, and analyzer
JP2014226871A (en) * 2013-05-24 2014-12-08 東レエンジニアリング株式会社 Method for simulating behavior of filler, and method for analyzing physical property of composite material
US20160342716A1 (en) * 2015-04-14 2016-11-24 Honda Motor Co., Ltd. Computer-aided resin behavior analyzer
JP2016198997A (en) * 2015-04-14 2016-12-01 本田技研工業株式会社 Computer-aided resin behavior analysis device

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6420881B1 (en) * 2017-10-31 2018-11-07 株式会社Jsol Compression molding analysis system, compression molding analysis method, and compression molding analysis program
JP2019082929A (en) * 2017-10-31 2019-05-30 株式会社Jsol Compression molding analysis system, compression molding analysis method, and compression molding analysis program
WO2020111039A1 (en) * 2018-11-29 2020-06-04 東レエンジニアリング株式会社 Movement analysis device for fibers in resin, movement analysis method for fibers in resin, program, and storage medium
WO2021002210A1 (en) * 2019-07-02 2021-01-07 本田技研工業株式会社 Resin behavior analysis device, resin behavior analysis method, and resin behavior analysis program
JP2021008081A (en) * 2019-07-02 2021-01-28 本田技研工業株式会社 Resin behavior analyzer, resin behavior analyzing method, and resin behavior analyzing program
JP7134924B2 (en) 2019-07-02 2022-09-12 本田技研工業株式会社 Resin behavior analysis device, resin behavior analysis method, and resin behavior analysis program

Also Published As

Publication number Publication date
JP6203787B2 (en) 2017-09-27

Similar Documents

Publication Publication Date Title
Bahloul et al. A study on optimal design of process parameters in single point incremental forming of sheet metal by combining Box–Behnken design of experiments, response surface methods and genetic algorithms
JP6203787B2 (en) Computer-aided resin behavior analyzer
WO2010073756A1 (en) Springback occurrence cause analyzing method, springback occurrence cause analyzing device, springback occurrence cause analyzing program, and recording medium
JP2016198997A (en) Computer-aided resin behavior analysis device
CN108108582A (en) A kind of method for numerical simulation of curved-surface piece flexible rolling forming process
Minguez et al. Framework for verification of positional tolerances with a 3D non-contact measurement method
JP2007229784A (en) Method for preparing correction shape data of press die
Shim Optimal preform design for the free forging of 3D shapes by the sensitivity method
JP2018183885A (en) Mold CAD model data creation apparatus and mold CAD model data creation method
JP2018069263A (en) Press working condition determination method
US10223481B2 (en) Computer-aided resin behavior analyzer
JP7407642B2 (en) Design support device and design support method
JP5969886B2 (en) Corresponding point calculation system and program, mold shape generation system and program
JP2014026301A (en) Setting method for material anisotropic information and plate thickness information for analytical model of molded article, stiffness analysis method and collision analysis method
KR101359661B1 (en) conversion system and method for 3D object represented by triangle mesh to 3D object represented by DoSurface
Bambach et al. Quantitative validation of FEM simulations for incremental sheet forming using optical deformation measurement
JP2016020000A (en) Creation method of metal mold expected shape and device
JP2000122996A (en) Method for evaluating shape defect of work forming surface
JP2011145876A (en) Die production method and creation method for die machining data
Cekic et al. Production of optimized layered products using intelligent support
JP2015056027A (en) Method for calculating contact surface pressure between two objects, and computer program for calculating contact surface pressure
US20220277117A1 (en) Resin behavior analysis apparatus, resin behavior analysis method and resin behavior analysis program
JP6265811B2 (en) Draw model generation method and draw model generation system
Dan et al. Direct generation of die surfaces from measured data points based on springback compensation
JP6303815B2 (en) Forging crack prediction method, forging crack prediction program, and recording medium

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170411

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170606

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170727

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170830

R150 Certificate of patent or registration of utility model

Ref document number: 6203787

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees