JP2005283539A - Analysis method for filler resin composition, analysis program and recording medium recording analysis program - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To perform three-dimensional strength analysis of a filler resin composition with high accuracy, and to reduce trouble of calculation necessary for the strength analysis. <P>SOLUTION: By performing two-dimensional CAE flow analysis of the filler resin composition, a 2D flow model divided into a plurality of mesh-like shell elements is configured, and orientation direction data are obtained in each of the plurality of shell elements (a flow analysis processing procedure 11). By assuming that a filler is not filled with respect to the filler resin composition and performing three-dimensional CAE structural analysis, a non-filler 3D structure model divided into a plurality of mesh-like solid elements is configured (a structural analysis processing procedure 12). Each the orientation direction data or the like are converted into an isotropic physical property value, are made to be projected to the non-filler 3D structure model to perform a mapping process (a mapping processing procedure 3), and the strength analysis is performed on the basis of an obtained 3D mapping model. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

  The present invention relates to a filler resin composition analysis method, an analysis program, and a recording medium recorded with the analysis program. For example, in order to increase the shape accuracy of a resin composition formed by injection molding or the like, It is possible to support optimal settings such as design and molding conditions.

  Compositions obtained by injection molding of resin materials (hereinafter referred to as resin compositions) are applicable as the resin materials, molding technology, and fillers (glass fibers, etc.) for reinforcing the molded body are improved. For example, the resin composition can be applied to an engine head cover, an intake pipe, and the like.

With the expansion of the application range of the resin composition, it becomes necessary to analyze the strength of the resin composition. For example, when a resin composition such as an engine head cover is an object to be analyzed, is the durability during engine operation sufficient? It is required to analyze whether or not in advance. As a system used for this analysis, for example, three-dimensional CAE (computer aided engineering) strength analysis using flow analysis, structural analysis, etc. by a finite element method is known (for example, Patent Documents 1 and 2 and Non-Patent Documents). Literature 1, 2, 3).
Japanese Unexamined Patent Publication No. 2002-14021 (for example, paragraphs [0002] to [0007]). Japanese Patent Laid-Open No. 4-102179 (for example, page 1 right column, line 18 to page 3, left column 6th line, FIGS. 1 and 3). Ryo Nakano, Katsuya Sakaba, Satoshi Sawada, Takashi Yuki, Yasuo Suga, “Development of 3D injection molding CAE system” (Japan), Molding, Vol. 15, No. 8, 2003, Japan Society for Plastic Molding, 550-555. Takeshi Oba, Kunio Machida, “Strength Prediction of GF Reinforced PA66 Thick Molded Product” (Japan), Molding '02, published on May 28, 2002, Japan Society for Plastic Molding Processing, pp. 345-346. Kazumi Tada, Masahiro Mori, “Current Sled Analysis” (Japan), Molding '99, published on May 28, 1999, Japan Society for Plastic Molding, 351-352.

  When the object to be analyzed is an anisotropic resin composition, that is, a resin composition using a resin material filled with a filler such as glass fiber (hereinafter referred to as a filler resin composition), the filler resin is simply used. Performing a CAE strength analysis in three dimensions (obtaining a model divided into a plurality of elements (hereinafter referred to as cell elements) in a mesh shape) based on data such as molding conditions and shapes of the composition Calculations required for analysis (particularly, three-dimensional calculations related to the orientation direction of the filler) are considered to be difficult due to enormous reasons.

  As a method for omitting the enormous calculation as described above, for example, the orientation direction (vector) of the filler in the filler resin composition is specified by speculation or specified experimentally (for example, the actual resin composition is cut and the filler A method of identifying (by observing the orientation direction) was known. In this method, the filler resin composition is subjected to a three-dimensional CAE structural analysis without considering the filler (that is, a three-dimensional CAE structural analysis is performed on the resin composition that is assumed not to be filled with filler). To obtain a three-dimensional structural analysis model (hereinafter referred to as a non-filler 3D structural model) divided into a plurality of elements (hereinafter referred to as solid elements) in the form of a mesh. This is a method of performing CAE strength analysis of the filler resin composition by setting the orientation direction data of the identified filler as an anisotropic physical property value (hereinafter referred to as an estimation / experimental method). However, the above-described estimation / experimental method is primitive because the orientation direction of the filler is manually specified, and it takes a lot of time and effort.

  In recent years, for example, by performing a two-dimensional CAE flow analysis on a filler resin composition, a two-dimensional analysis model (hereinafter referred to as a 2D flow model) composed of a plurality of elements (hereinafter referred to as shell elements) in a mesh shape. In general, a method of analyzing the strength of the filler resin composition by simply converting the orientation direction data of the filler relating to the 2D flow model into an anisotropic physical property value is employed.

  As described above, the method of analyzing the strength by converting the orientation direction data of the filler related to the 2D flow model into the anisotropic physical property value corresponds to the method of simply performing the two-dimensional CAE strength analysis of the filler resin composition. Since it is difficult to perform a three-dimensional CAE analysis with respect to the orientation direction, it is possible to reduce the labor required for the calculation as compared with a method in which the filler resin composition is simply subjected to a three-dimensional CAE strength analysis.

  However, for example, in the case of a filler resin composition having a large undulation (for example, uneven thickness) or a relatively large capacity filler resin composition, the orientation direction of the filler in the actual filler resin composition changes variously. It is difficult to accurately specify the orientation direction in three dimensions, and sufficient modeling cannot be achieved. Therefore, the accuracy of strength analysis of the filler resin composition is reduced, that is, a result far from the actual measurement value is obtained. There is a fear.

  The present invention has been made based on the above problems, and is a filler resin capable of analyzing the strength of a filler resin composition three-dimensionally with good accuracy and reducing the labor of calculation required for the strength analysis. The object is to provide a composition analysis method and an analysis program.

  The present invention is intended to solve the above problems, and the inventions according to claims 1 to 8 are based on a filler resin composition containing a filler having anisotropy based on at least molding conditions of the filler resin composition. A three-dimensional strength analysis is performed by creating a 3D mapping model, and a two-dimensional CAE flow analysis is performed on the filler resin composition. By creating a 2D flow model in which the orientation direction data of the filler is specified, and assuming that the filler resin composition is not filled with the filler resin composition, a three-dimensional CAE structure analysis is performed. A non-filler 3D structural model is created, and each orientation data is converted into an anisotropic physical property value for each solid element. It is by mapping process by the projection, characterized by creating a 3D mapping model comprising a plurality of cells elements.

  In addition, the invention according to claims 9 to 13 includes a plurality of filler resin compositions containing an anisotropic filler by performing a two-dimensional CAE flow analysis based on at least the molding conditions of the filler resin composition. A flow analysis processing procedure for creating a 2D flow model in which the orientation data of the filler is specified for each shell element and storing the 2D flow model in the storage means; and the filler resin composition A non-filler 3D structural model composed of a plurality of solid elements is created by performing a three-dimensional CAE structural analysis based on molding conditions assuming that at least the filler is not filled, and the filler 3D structural model Are stored in the storage means, and the orientation direction data are transmitted via the calculation means. A 3D mapping model for creating a 3D mapping model composed of a plurality of cell elements by converting the data into anisotropic physical property values and projecting each of the solid elements to perform a mapping process. A program for performing a three-dimensional intensity analysis based on a mapping model, which is executed via a computer.

  The present invention obtains a 2D flow model that specifies the orientation direction of the filler by performing a two-dimensional CAE flow analysis on the filler resin composition, and maps the orientation direction data of the 2D flow model to a non-filler 3D structure model. Thus, the orientation direction is specified three-dimensionally, and an enormous calculation such as a three-dimensional CAE intensity analysis is not required.

  As described above, in the present invention, the orientation direction is specified three-dimensionally, and an enormous calculation such as a three-dimensional CAE intensity analysis is not required. Therefore, the filler resin composition can be obtained with good accuracy. In addition to being able to perform three-dimensional strength analysis (ie, analysis in the same manner as three-dimensional CAE strength analysis), it is possible to reduce the labor of calculation required for the strength analysis. For example, a resin composition formed by injection molding or the like In order to increase the shape accuracy, it is possible to optimize the design and molding conditions of the mold.

  Hereinafter, a filler resin composition analysis method and analysis program according to an embodiment of the present invention will be described in detail with reference to the drawings.

  The present embodiment does not simply perform a two-dimensional or three-dimensional CAE strength analysis as in the prior art, but first, by performing a two-dimensional CAE flow analysis on the filler resin composition to be analyzed, By obtaining a 2D flow model in which the orientation direction is specified and mapping the orientation direction data to a non-filler 3D structure model, the filler resin composition is three-dimensionally analyzed for strength.

  That is, in the present embodiment, for example, a filler resin is used by using a computer provided with a calculation means, a storage means, and an input means (for example, a user interface for inputting data relating to the filler resin composition, molding condition data, shape data, etc.). A two-dimensional CAE flow analysis is performed on the composition based on at least the molding conditions of the filler resin composition by a finite element method or the like to form a 2D flow model divided into a plurality of shell elements in a mesh shape, and a plurality of Orientation direction data (two-dimensional data) is obtained for each shell element.

  Further, assuming that the filler resin composition is not filled with filler, a three-dimensional CAE structure based on at least molding conditions (molding conditions when filler is not filled) by a finite element method or the like. By performing the analysis, a non-filler 3D structural model is formed which is divided into a plurality of solid elements in a mesh shape.

  Then, the orientation direction data and the like are converted into anisotropic physical property values, projected onto the non-filler 3D structure model, mapped, and the mapped three-dimensional model (that is, obtained by the mapping process). Strength analysis is performed based on a model including a plurality of cell elements (hereinafter referred to as a 3D mapping model).

  In the mapping process, first, the shell element is divided into a plurality of layers (hereinafter referred to as shell division layers) in the thickness direction (for example, the thickness direction of the target filler resin composition), for example, a surface layer, It is divided into a core layer (a layer opposite to the surface layer) and an intermediate layer (a layer located approximately in the middle between the surface layer and the core layer).

  Then, select one of the above-mentioned shell divided layers, for example, in the case of mainly analyzing the behavior of deformation under external stress, select the surface layer, and by thermal expansion according to the ambient temperature used When analyzing the deformation behavior mainly, select an intermediate layer (an intermediate layer relatively thicker than the surface layer), and analyze the behavior of complex deformation due to external stress and thermal expansion. In this case, an average layer taking an average from the surface layer to the core layer is selected, and orientation direction data and the like in the selected shell division layer are converted into anisotropic physical property values and projected.

  Since each shell element and each solid element have different dimensions, shapes, coordinate positions, etc., each shell element and each solid element are appropriately matched and the selected orientation direction data is mapped. In some cases, it may be difficult (that is, orientation direction data or the like is not mapped to a desired coordinate position, and a result far from the actual measurement value may be obtained).

  Moreover, the orientation direction of the filler in an actual filler resin composition is easily affected by molding conditions (injection molding conditions), and the filler resin composition has a shape (for example, a shape with a large undulation) or the filler resin composition. Varies three-dimensionally according to the distribution position (for example, surface layer, intermediate layer) and the like. For shell elements where such changes are significant, there is a large difference in the orientation direction of each shell split layer, and accurate data (that is, the orientation direction that matches or is approximately the same as the measured value) among the orientation direction data. If the mapping process is not performed by selecting (data), there is a case where an error from the actual measurement value becomes large and the accuracy of the intensity analysis is lowered.

  As described above, when there is a possibility that the orientation direction data or the like may not be mapped to the desired coordinate position, first, the center of gravity position of the solid element that is coincident with or closest to the center of gravity of each shell element is calculated, respectively. It is preferable that the mapping process is performed by converting the orientation direction data of the shell element into an anisotropic physical property value and projecting it to the solid element in the combination of the shell element and the solid element.

  In addition to calculating the combination of the shell element and the solid element as described above, the orientation direction data of each shell division layer in the shell element is averaged, and the averaged orientation direction (hereinafter, average orientation) is calculated. By converting the data (referred to as the direction) data, that is, the average physical layer into anisotropic physical property values and projecting them, it is possible to maintain good accuracy in the intensity analysis.

  FIG. 1 is a schematic explanatory diagram showing an example of an analysis method and an analysis program in the present embodiment. In the contents shown in FIG. 1, using a computer equipped with arithmetic means, storage means, input means, etc., software for general two-dimensional CAE flow analysis, software for three-dimensional CAE structure analysis, mapping It shall be executed using processing software.

  In FIG. 1, the code | symbol 11 shows a flow analysis process sequence, The molding conditions (for example, injection speed, resin temperature, holding pressure value, holding pressure time of the filler resin composition which is analysis object via an input means. ), Data of materials used (resin, filler, etc.), shape, etc.) are input, and based on these input values, two-dimensional CAE flow analysis is performed by a calculation means. As a result, a 2D flow model divided into a plurality of shell elements in a mesh shape is obtained, and node information, element information, filler orientation direction (orientation direction for each shell division layer), etc. are two-dimensionally for each shell element. Can be analyzed. Each data obtained by the above analysis is stored in a storage means, for example.

  Reference numeral 12 denotes a structural analysis processing procedure. Assuming that the filler resin composition is not filled with filler, the molding conditions (filler filling is performed) through the input means in the same manner as the flow analysis processing procedure described above. Molding conditions assumed not to be entered; for example, injection speed, resin temperature, holding pressure value, holding pressure time, used material (resin etc.), shape, etc.) data are input and based on these input values A three-dimensional CAE structure analysis is performed by an arithmetic means. Thereby, a non-filler 3D structure model divided into a plurality of solid elements in a mesh shape is obtained, and node information, shape, physical properties, etc. can be analyzed three-dimensionally for each solid element. Each data obtained by the above analysis is stored in a storage means, for example.

  In the flow analysis process procedure 11 and the structure analysis process procedure 12 described above, for example, a constraint condition (for example, a condition to be constrained to a reference position), a boundary condition, or the like in consideration of a subsequent process (for example, a mapping process described later). May be added. The node information includes, for example, the coordinate position of a node with the reference position as the origin (each shell element, the coordinate position of a node of each solid (hereinafter referred to as the node coordinate position), the number of nodes, and the like. Furthermore, the element information includes, for example, the distribution position, shape, physical properties, number of elements, etc. in the filler resin composition of each shell element and each solid element.

  Thereafter, in the two-dimensional data extraction processing procedure indicated by reference numeral 21, the node information data, element information data, orientation direction data, etc. of each shell element stored in the storage means are extracted (for example, written in a text format). In the three-dimensional data extraction processing procedure indicated by 22, node information data, element information data, etc. of each solid element stored in the storage means are extracted (for example, written in a text format). The data of each shell element and each solid element is preferably distinguished by adding a file name for each shell element and each solid element, for example.

  Reference numeral 3 indicates a mapping processing procedure (details will be described later). For example, the data for each shell element is converted into an anisotropic physical property value through software including a general mapping program and arithmetic means. Then, these anisotropic physical property values are projected onto the data for each solid element and sequentially mapped to create a 3D mapping model composed of a plurality of cell elements in a mesh form.

  Then, according to the mapping determination procedure of reference numeral 4, it is confirmed whether or not the mapping process for each shell element and each solid element has been completed (that is, the creation of all cell elements has been completed) (if necessary) It is possible to perform strength analysis (that is, three-dimensional strength analysis) based on the obtained 3D mapping model.

  FIG. 2 is a schematic explanatory diagram illustrating an example of the mapping process procedure 3 of FIG. 1 (for example, an example of using I-DEAS (registered trademark)). The data of each shell element and each solid element, and the orientation direction and average orientation direction data of each shell division layer are the file names suitable for the mapping processing procedure 3 (hereinafter referred to as shell F, solid F, and orientation direction F, respectively). , Referred to as average F), and for example, an extension is added based on the flowchart shown in FIG. 3 and then stored in the storage means or the like. It shall be executed via

  In FIG. 2, reference numeral 30 denotes at least one or more arbitrary shell elements (shell elements including node information data, element information data, orientation tendency data, etc.) among the shell elements described above. A processing procedure for inputting the shell F of the shell element via an input means or the like (hereinafter referred to as a shell specifying processing procedure) is shown.

  Reference numeral 31 denotes at least one arbitrary solid element (node information data, element information data) that approximates the selected shell element (for example, the x and y coordinate positions coincide) among the solid elements. The processing procedure (hereinafter referred to as a solid identification processing procedure) for inputting the solid F of the solid element via the input means or the like is shown.

  Reference numeral 32 selects (or averages) the orientation direction of the shell divided layer according to the molding condition and shape of the filler resin composition, the distribution position of the shell element in the filler resin composition, etc. This shows a processing procedure (hereinafter referred to as an orientation direction selection processing procedure) in which the selected orientation direction F (or average F) is input via an input means or the like.

  Reference numeral 33 denotes a processing procedure (hereinafter referred to as an allowable value input processing procedure) for inputting an allowable value related to the target 3D mapping model by the input means, and each data of the shell F and the solid F is stored in the processing procedure. It is for determining whether the 3D mapping model is effective or not through an arithmetic means or the like, and deleting unnecessary shell F and solid F.

  Reference numeral 34 denotes an output file name (hereinafter referred to as output) for reading the data of the shell F, solid F, and orientation direction F (or average F) in each subsequent processing procedure (reference numerals 35 to 39 described later). A process procedure (hereinafter referred to as an output F definition process procedure) is defined.

  Reference numeral 35 denotes a processing procedure related to the solid element of the input solid F (hereinafter referred to as a solid element processing procedure). The node information data, element information data, etc. in the solid element are read and calculated. The coordinate position of the center of gravity of the solid element (hereinafter referred to as the solid center of gravity coordinates) is calculated through means or the like.

  Reference numeral 36 denotes a processing procedure related to the shell element of the input shell F (hereinafter referred to as a shell element processing procedure). Node information data, element information data, orientation direction data, etc. in the shell element And the coordinate position of the center of gravity of the shell element (hereinafter referred to as the shell center-of-gravity coordinates) is calculated via an arithmetic means or the like.

  Reference numeral 37 denotes a shell element that calculates a distance between each of the calculated solid center-of-gravity coordinate positions and each solid center-of-gravity coordinate position (hereinafter referred to as a distance between the center-of-gravity coordinates). And a solid element (for example, a combination in which a solid barycentric coordinate position and a solid barycentric coordinate position match) are obtained through calculation means (required by repeated calculation if necessary) (hereinafter, referred to as a procedure) , Referred to as a combination definition processing procedure).

  Reference numeral 38 calculates a vector of a line segment (ie, equivalent to a filler) indicated by the orientation direction data (or average orientation direction data) of the inputted orientation direction F (or average F) (ie, anisotropic). A processing procedure (hereinafter referred to as 3D) for specifying the three-dimensional coordinate position in the solid element (solid element related to the combination obtained in the combination definition processing procedure 37) with respect to the orientation data by the calculation means or the like. This is referred to as a coordinate calculation processing procedure). In this 3D coordinate calculation processing procedure 38, when the line segment of the target orientation direction data is AP, the vector of the line segment AP is calculated based on the following calculation formulas (1) to (17). Can do.

First, in the solid element corresponding to the orientation direction to be calculated (that is, the solid element onto which the anisotropic physical property value is projected), as shown in FIG. Are (a ₁ , b ₁ , c ₁ ), (a ₂ , b ₂ , c ₂ ), (a ₃ , b ₃ , c ₃ ), and the end point (ie, P point) of the line segment AP is the line segment. Assuming that it is located on BC, a straight line passing through the point BC can be expressed by the following equation (1).

  When the three-dimensional coordinate position at the end point of the vector AP is (x, y, z), the following equations (2) to (4) are established.

The above equations (2) to (4) are expressed as follows: a ₃ −a ₂ = a ₃₂ , b ₃ −b ₂ = b ₃₂ , c ₃ −c ₂ = c ₃₂ , a ₂ −a ₁ = a ₂₁ , b ₂ When −b ₁ = b ₂₁ and c ₂ −c ₁ = c ₂₁ , they can be expressed by the following equations (5) to (7), respectively.

  Further, the following expression (9) is established from the following expression (8).

Assuming that a ₂₁ ² + b ₂₁ ² + c ₂₁ ² = L ₂₁ , the above formula (9) can be expressed by the following formula (10), and formulas (11) and (12) hold.

  From the above equations (8) to (12), the following equation (13) is established.

The above equation (13) is obtained by the following equation (14), assuming that (a ₃₂ ² + b ₃₂ ² + c ₃₂ ² ) = L ₃₂ , a ₂₁ · a ₃₂ + b ₂₁ · b ₃₂ + c ₂₁ · c ₃₂ = M (15) is established.

The equation (15) _{is, (L 21 · L 32 ·} cos 2 θ-M 2) = X, 2 · L 21 · M (cos 2 θ-1) = Y, L 21 2 · (cos 2 θ- If 1) = Z, it can be expressed by the following equation (16), and equation (17) is established.

  Reference numeral 39 reads, for example, the output F and extracts node information data, element information data, orientation direction data, vector data, etc. input or calculated in the shell specifying processing procedure 30 to 3D coordinate calculation processing procedure 38. , Shows a processing procedure (hereinafter referred to as a mapping cell element processing procedure) for defining a cell element of a target 3D mapping model based on the respective data by an arithmetic means or the like, and various strength analyzes ( For example, analysis based on Young's modulus, Poisson's ratio, etc. can be performed. If necessary, node information data, mapping information data, and strength information (for example, Young's modulus, Poisson's ratio, vector information, etc.) data of the cell element can be extracted and edited via a desired editor. good.

Next, when the average orientation direction is selected in the alignment direction selection processing procedure 32, a method for calculating the average orientation direction will be described with reference to FIGS. 5A and 5B. First, the number of shell division layers of the shell element is N (N is a natural number (1, 2,..., N-1, N)), the orientation direction of each shell division layer is AP ₁ to AP _N , and the target average orientation. The direction is AP _AVE. An average value of angles with respect to the reference position in each of the orientation directions AP _{1 to} AP _N (hereinafter referred to as a representative angle θq) is obtained.

Thereafter, the vector of the representative angle θq is defined as AQ, and the angle difference between the AQ and each of the orientation directions AP _{1 to} AP _N (hereinafter referred to as θ _{1 to} θ _N ) is obtained, respectively. _{In 1 to} AP _N , components in the direction parallel to AQ (that is, AP ₁ cos θ _{1 to} AP _N cos θ _N ) and components in the vertical direction (that is, AP ₁ sin θ _{1 to} AP _N cos θ _N ) are calculated.

Then, by dividing the sum of AP ₁ cos θ _{1 to} AP _N cos θ _N and AP ₁ sin θ _{1 to} AP _N cos θ _N by N, respectively, the component in the direction parallel to the AQ and the component in the vertical direction are Each average value is calculated, and AP _AVE. Can be specified.

[Example]
Next, in the analysis system and analysis method of the present embodiment, verification was performed using the following first and second examples.

(First embodiment)
In the first embodiment, as shown in FIG. 6, a sample of a filler resin composition made of a polyamide resin filled with glass fibers in a shape in which a member protruding in the vertical direction from one end surface side of a substantially flat plate-like member is formed. S1 was created, a 3D mapping model was created in the area surrounded by the dotted line in the figure, and the accuracy of the mapping process was verified.

  First, a two-dimensional CAE flow analysis is performed on the dotted line region of the sample S1 based on the flow analysis processing procedure 11 shown in FIG. In the example, a 2D flow model in which orientation direction data and the like were specified for every 260 shell elements) was obtained. Further, assuming that the filler is not filled with respect to the dotted line region, a three-dimensional CAE structural analysis is performed based on the structural analysis processing procedure 12 shown in FIG. In the example, a non-filler 3D structural model consisting of 537 solid elements) was obtained.

  Then, based on the mapping process procedure 3 shown in FIG. 1 (shell identification process procedure 30 to mapping cell element process procedure 39 shown in FIG. 2) and the mapping determination procedure 4, the orientation direction data of each shell element is obtained. A 3D mapping model related to the sample S1 as shown in FIG. 8 was created by converting it into an anisotropic physical property value, projecting it onto each solid element, and performing mapping processing.

  Here, the angle of the average orientation direction (angle with respect to the reference position) of each shell element (shell element numbers 1 to 260 in Table 1 below) and each cell element in the 3D mapping model corresponding to each average orientation direction The difference (angle difference in Table 1 below) from the angle in the orientation direction (angle relative to the reference position) was determined, and the results are shown in Table 1 below.

  Further, the distance between the center of gravity coordinates between the shell element and the solid element (solid element numbers 1 to 537 in Table 2 below) of each combination identified by the combination definition processing procedure 37 is calculated via a general theoretical formula (for example, The distance between the center of gravity coordinates (theoretical value) calculated using software “EXCEL” manufactured by Microsoft Corporation was compared, and the results are shown in Table 2 below.

  In Table 1, it was determined that the accuracy of the mapping process was good (the symbol “◯” in Table 1 below) when the angle difference was less than 0.01 degrees. The distance between the barycentric coordinates calculated through the above theoretical formula is obtained by paying attention to the fact that the nodes constituting each element have coordinates (X, Y, Z). It shows the distance between the theoretical center-of-gravity position coordinates in each element calculated using (relatively simple calculation formula).

  As shown in Table 1 and Table 2, the difference between the average orientation direction angle in the two-dimensional analysis model and the orientation direction angle in the 3D mapping model is less than 0.01 degrees, respectively. It can be read that each identified combination and the distance between the barycentric coordinates coincide with the theoretical value. From this, it has been found that according to the analysis system and analysis method of the present embodiment, even when the mapping process is performed by applying the average orientation direction, the mapping process can be performed with good accuracy with respect to the filler resin composition.

(Second embodiment)
In the second embodiment, a resin material filled with glass fibers is used, as shown in the schematic explanatory views of FIGS. 9A (schematic plan view) and B (A-A cross-sectional view). A sample S2 of a head cover (in this example, a head cover for a water-cooled in-line three-cylinder SOHC engine) was prepared (in this example, formed using an injection molding machine (made by Nippon Steel Works; 650ton)), and the sample S2 A 3D mapping model is created, and strength analysis (deformation) assuming the reaction force of the rubber gasket (rubber gasket interposed between the cylinder head of the engine and the cylinder head cover) for each part (reference numerals B1 to B36 in FIG. 9) Analysis). The main molding conditions for the sample S2 are shown in Table 3 below.

  First, as in the first embodiment, a two-dimensional CAE flow analysis is performed on the sample S2 based on the flow analysis processing procedure 11 shown in FIG. A 2D flow model in which orientation direction data was specified was obtained. Further, assuming that the filler is not filled based on the structural analysis processing procedure 12 shown in FIG. 1, a non-filler 3D structure model composed of a plurality of cell elements in a three-dimensional mesh shape was obtained.

  Further, based on the mapping processing procedure 3 shown in FIG. 1 (shell identification processing procedure 30 to mapping cell element processing procedure 39 shown in FIG. 3) and the mapping determination procedure 4, the orientation direction data and the like of each shell element is obtained. A 3D mapping model related to the sample S2 was created by converting it into an anisotropic physical property value, projecting it onto each solid element, and performing mapping processing.

  And about each part B1-B37 of the said 3D mapping model, this deformation | transformation amount is calculated | required by performing the deformation | transformation analysis supposing the reaction force of a rubber gasket, The deformation | transformation amount (mm) characteristic with respect to the measurement part of FIG. Shown in the figure. In addition, as a comparative example with the result of the deformation analysis, the strength analysis of the filler resin composition is performed on the sample S2 by converting the orientation direction data of the filler relating to the 2D flow model into simply anisotropic physical property values. The deformation analysis result and the actual measurement value obtained by the method of FIG. 10 are obtained, and these results are also shown in the characteristic diagram of FIG.

As shown in FIG. 10, it can be read that the deformation analysis result obtained by this embodiment approximates the actual measurement value as compared with the conventional method. The Young's modulus in the parallel direction E1 and the vertical direction E2 with respect to the orientation direction was obtained based on the 3D mapping model of the sample S2 obtained in the present embodiment, and about 9555 N (975.0 kgf / mm ² ) and about 4868.6N (507.0 kgf / mm ² ), and the Poisson's ratio ν was 0.36.

  Although the present invention has been described in detail only for the specific examples described above, it is obvious to those skilled in the art that various changes and modifications are possible within the scope of the technical idea of the present invention. Such variations and modifications are naturally within the scope of the claims.

  For example, a computer system including the above-described arithmetic means, storage means, input means, etc. is configured by a well-known CPU, ROM, RAM, computer main body having an I / O interface, keyboard, display, external memory, printer, etc. It is obvious that the above-described processing procedures and determination procedures are specifically realized by the CPU executing various program modules stored in the external memory, and a desired 3D mapping model can be obtained. is there.

  Further, the filler of the filler resin composition is not limited to glass fiber, and it is obvious that an anisotropic filler (for example, carbon, ceramic, etc.) can be applied.

  Further, the filler resin composition is not limited to products such as an engine head cover and an intake pipe, and it is apparent that it can be applied to various resin products.

Schematic explanatory drawing which shows an example of the analysis method and analysis program in this Embodiment. FIG. 6 is a schematic explanatory diagram for illustrating an example of a mapping process procedure 3; FIG. 10 is a schematic explanatory diagram for illustrating an example of a method for adding an extension or the like to a file suitable for the mapping processing procedure 3; Schematic explanatory drawing for showing an example of the 3D coordinate calculation processing procedure 38. Schematic explanatory drawing for showing an example of the calculation method of an average orientation direction. The schematic explanatory drawing of sample S1 of the filler resin composition in 1st Example. The schematic explanatory drawing for showing the 2D flow model in sample S1. Schematic explanatory drawing for showing the 3D mapping model in sample S1. The schematic explanatory drawing of the sample S2 of the engine head cover in 2nd Example. The deformation | transformation amount characteristic view with respect to the measurement location in sample S2.

Explanation of symbols

DESCRIPTION OF SYMBOLS 11 ... Flow analysis processing procedure 12 ... Structural analysis processing procedure 21 ... Two-dimensional data extraction processing procedure 22 ... Three-dimensional data extraction processing procedure 3 ... Mapping processing procedure 32 ... Orientation direction selection processing procedure 33 ... Allowable value input processing procedure 37 ... Combination Definition processing procedure 38 ... 3D coordinate definition processing procedure 4 ... Mapping determination procedure

Claims (13)

A method of conducting a three-dimensional strength analysis by creating a 3D mapping model of a filler resin composition containing a filler having anisotropy based on at least molding conditions of the filler resin composition,
By performing a two-dimensional CAE flow analysis on the filler resin composition, a 2D flow model consisting of a plurality of shell elements and specifying filler orientation direction data for each shell element is created,
Assuming that the filler resin composition is not filled with the filler resin composition, a three-dimensional CAE structure analysis is performed to create a non-filler 3D structural model composed of a plurality of solid elements,
A 3D mapping model composed of a plurality of cell elements is created by converting each orientation direction data into an anisotropic physical property value and projecting it onto each solid element to perform mapping processing. Analysis method of filler resin composition.   The center-of-gravity position of each shell element and the center-of-gravity position of the solid element closest to the center-of-gravity position are calculated, and the orientation data of the shell element is anisotropic from the calculated shell element and solid element. 2. The method for analyzing a filler resin composition according to claim 1, wherein each mapping process is performed by converting into physical property values and projecting them onto solid elements. Dividing the shell element in the thickness direction to obtain a plurality of shell division layers,
One of the shell division layers is selected, filler orientation direction data in the selected shell division layer is converted into an anisotropic physical property value, and projected onto the solid element for mapping processing. The analysis method of the filler resin composition of Claim 1 or 2.   4. The selected shell dividing layer is selected according to at least molding conditions and shape of the filler resin composition and a distribution position of shell elements in the filler resin composition. The analysis method of the filler resin composition of description.   The average orientation direction data of each shell division layer is averaged, the averaged average orientation direction data is converted into anisotropic physical property values, projected onto the solid elements, and mapped. The method for analyzing a filler resin composition according to claim 3. The average orientation direction is AP _AVE. , The number of shell division layers of the shell element is N, and the orientation directions of the shell division layers are AP ₁ to AP _N.
Obtains a representative angle θq as an average value of the angle with respect to the reference position in each orientation direction AP ₁ ~AP _N above,
The angle differences θ _{1 to} θ _N between the vector AQ of the representative angle θq and the orientation directions AP _{1 to} AP _N are respectively determined.
In each orientation direction AP ₁ ~AP _N above, it calculates the direction of the component AP ₁ cosθ ₁ ~AP _N cosθ _N parallel to AQ, perpendicular direction component _{AP 1 sinθ 1 ~AP N cosθ N} , respectively ,
By dividing the sum of the AP ₁ cos θ _{1 to} AP _N cos θ _{N and} the sum of AP ₁ sin θ _{1 to} AP _N cos θ _N by N, the average for the component in the direction parallel to the AQ and the component in the vertical direction Calculate each value,
6. The method for analyzing a filler resin composition according to claim 5, wherein the average orientation direction AP _AVE. Is specified by the calculated average value and the representative angle θq. The shell element includes at least node information data, element information data, orientation tendency data,
The filler resin composition analysis method according to claim 1, wherein the solid element includes at least node information data and element information data.   8. The method for analyzing a filler resin composition according to claim 1, wherein the filler comprises at least one of anisotropic glass fiber, carbon, and ceramic. By performing a two-dimensional CAE flow analysis on the filler resin composition containing the filler having anisotropy based on at least the molding conditions of the filler resin composition, the shell resin is composed of a plurality of shell elements and at least a filler for each shell element. A flow analysis processing procedure for creating a 2D flow model in which the orientation direction data is specified and storing the 2D flow model in the storage means;
A non-filler 3D structural model composed of a plurality of solid elements is created by performing a three-dimensional CAE structural analysis based on molding conditions when it is assumed that the filler resin composition is not filled with at least the filler resin composition. , A structural analysis processing procedure for storing the filler 3D structural model in the storage means,
A 3D mapping model composed of a plurality of cell elements is obtained by converting each orientation direction data into anisotropic physical property values and projecting them onto the respective solid elements and performing mapping processing via the arithmetic means. Mapping process procedure to create
A filler resin composition analysis program, comprising: a computer program for performing a three-dimensional strength analysis based on the 3D mapping model. The mapping process procedure is as follows:
Shell element processing procedure for extracting node information data, element information data, orientation direction data in at least one shell element among each of the shell elements, and calculating shell centroid coordinates via a calculation means,
Among the solid elements, the node information data and element information data are extracted for at least one solid element whose x and y coordinate positions match or approximate to the shell element, and the solid barycentric coordinates are calculated via the calculation means. Solid element processing procedure to
A combination obtained by calculating the distance between the center of gravity coordinates between each calculated solid center of gravity coordinate position and each solid center of gravity coordinate position, and the combination calculation means of the shell element and the solid element having the minimum distance between the center of gravity coordinates Definition processing procedure,
A vector of a line segment indicated by any one orientation direction data is calculated from the shell division layers obtained by dividing the shell element according to the combination into a plurality of three-dimensional coordinate positions in the solid element with respect to the orientation direction data. 3D coordinate calculation processing procedure for specifying
Based on the node information data, element information data, orientation direction data, and vector data extracted or calculated by the shell element processing procedure, solid element processing procedure, combination definition processing procedure, 3D coordinate calculation processing procedure, the target 3D mapping model A mapping cell element processing procedure for defining a cell element of
The filler resin composition analysis program according to claim 9, comprising: The orientation direction data specified in the 3D coordinate calculation processing procedure is as follows.
Among the shell divided layers, at least the molding condition and shape of the filler resin composition, and the orientation direction data of the shell divided layer selected according to the distribution position of the shell elements in the filler resin composition The analysis program for a filler resin composition according to claim 10, wherein   11. The filler according to claim 10, wherein the orientation direction data specified by the 3D coordinate calculation processing procedure is average orientation direction data obtained by averaging the orientation direction data of the shell divided layers. Analysis program for resin composition. Each procedure in the analysis program of the filler resin composition according to claim 9 to 12 is recorded on a recording medium readable by the computer as a program to be executed by the computer.
The recording medium which recorded the analysis program of the filler resin composition characterized by the above-mentioned.

Images