JP6061982B2 - Method for structural analysis of anisotropic resin molding - Google Patents

Description

  The present invention relates to a method for creating a model for structure analysis of an anisotropic resin molded body, a method for structure analysis, and a program for creating a model for structure analysis.

  Anisotropic resin moldings obtained by injection molding etc. of anisotropic resin compositions containing thermoplastic resins and anisotropic fillers are widely used as industrial products, and shape analysis by numerical analysis is not possible. It is used a lot. By the way, physical strength values such as elastic modulus, Poisson's ratio, and linear expansion coefficient are required for the strength design of the resin molded body by numerical analysis. When the resin molding contains an anisotropic filler, it is not easy to actually measure and predict these physical property values.

  The elastic modulus calculation method described in Patent Literature 1 stores in advance the first elastic modulus in the fiber orientation direction of the object, the second elastic modulus in the direction orthogonal to the fiber orientation, and the load direction when performing the structural analysis. The first elastic modulus is obtained when the fiber orientation direction and the load direction are parallel for any element of the object, and the first elastic modulus is obtained when the fiber orientation direction and the load direction are orthogonal to each other. 2 elastic modulus is set respectively, and when the fiber orientation direction and the load direction are neither parallel nor orthogonal, the angle between the load direction and the fiber orientation direction and either the first or second elastic modulus The elastic modulus is calculated based on the elastic modulus, and the elastic modulus is set.

  In Patent Document 2, there is a technique of obtaining a tensor that expresses the physical properties of the entire composite material by adding and calculating the physical properties of each single fiber and the virtual strain tensor for handling the contained fibers as equivalent inclusions equivalent to a matrix. It has been introduced.

  In Patent Document 3, the flow orientation analysis in the evaluation shape is performed using the mechanical properties of the reinforcing fiber and the base resin alone contained, and the mechanical properties calculated from the calculated orientation parameters of the resin molded product are actually converted into the evaluation shapes. Compare the measured physical properties of the molded product with the measured values, identify the base resin that the calculated value is equal to the actual measured value, and the mechanical properties of the reinforcing fiber alone, and use the identified mechanical properties to deform the molded part A method for evaluating the quantity with high accuracy is described.

  However, in any of the methods, there is a numerical discrepancy between the predicted value of the mechanical property prediction method of the molded product by fiber orientation and the actual physical property value of the molded product, and further improvement has been demanded.

  As a cause of this discrepancy, there is a difference between an assumption and a case where it is obtained by calculation. For example, the anisotropic filler breaks during injection filling, and the value of the aspect ratio of the anisotropic filler used in the calculation is different from the actual value, the fiber content is uneven, Insufficient predictive power of flow orientation when anisotropic filler is contained due to interference between the isotropic fillers. In other words, in addition to the fact that the calculation assumptions can not express the actual phenomenon sufficiently, the physical property values of each component are inaccurate, so the error of the input physical property value is added to the error at the time of building the analysis model and compounded There was an error, and even if the mechanical properties were predicted from the filler orientation calculation results, this prediction was insufficient.

  In order to improve the accuracy of prediction, in Patent Document 4, a plurality of slice images are acquired by X-ray CT measurement or the like for a resin molded product obtained by molding a resin composition containing a filler at a predetermined ratio, and slice images are obtained. A finite element method model that reflects the orientation of the filler in the resin molded product is created based on the information divided into pixels with a small area, and mechanical properties are identified through finite element method analysis. A method for evaluating the strength and the deformation amount with high accuracy is described.

  Moreover, in patent document 5, the three-dimensional orientation data of the filler in the injection molded product containing the filler is acquired from a three-dimensional image by X-ray CT, and homogenized based on the integral value of the shear stress and the data of the molecular orientation state. The linear expansion coefficient is obtained through the method, and the hot warpage of the injection molded product is obtained by structural analysis.

  However, the methods described in Patent Documents 4 and 5 require a great deal of experimentation and measurement in order to obtain the physical property values of the constituent elements, and the calculation time for obtaining them is also great. Thereby, the cost required for calculation when determining the physical property value is enormous.

Japanese Patent Laid-Open No. 9-237267 Japanese Patent Laid-Open No. 7-304056 JP 2004-25796 A JP 2011-758 A JP 2014-1000087 A

  The present invention has been made in order to solve the above-described problems, and in the structural analysis of an anisotropic resin molded body that is considered to be difficult to analyze the structure, the accuracy of the structural analysis is improved. An object of the present invention is to provide a structural analysis model that can significantly reduce the load on the computer in a short time with the analysis result.

  The inventors of the present invention have made extensive studies to solve the above problems. As a result, assuming a plurality of types of virtual molded bodies having different orientation states from the anisotropic resin molded body, for each of the plurality of types of virtual molded bodies, physical property information including at least one of the Poisson's ratio and the linear expansion coefficient is previously stored. By creating and applying this physical property information to each of the plurality of second regions divided into elements for structural analysis, it has been found that the above problems can be solved, and the present invention has been completed. More specifically, the present invention provides the following.

[1] The present invention is a structural analysis model creation method for creating a structural analysis model for performing structural analysis of an anisotropic resin molded body containing an anisotropic filler, the anisotropic resin molding For each of a plurality of types of virtual molded bodies having different orientation states from the body, physical property information including at least one of the Poisson's ratio defined by the following formula (1) and the linear expansion coefficient defined by the following formula (2) is created. A physical property information creating step, a first region orientation state calculating step for calculating a first region orientation state in each of the plurality of first regions divided in the resin flow analysis of the anisotropic resin molding, A search step for searching for the first region closest to the second region for each of the plurality of second regions divided into elements for the structural analysis of the anisotropic resin molded body, and a search step that is closest in the search step 1 area A second region orientation state setting step in which the first region orientation state in the second region orientation state is set to each second region orientation state, and the physical property information is referred to, and each of the plurality of second regions is changed to the second region orientation state. And a second region physical property information setting step in which the corresponding physical property information is used as each second region physical property information.
(1)
(In Formula (1),
ν ₂₃ is a Poisson's ratio to two directions orthogonal to the principal axis direction of the target material,
ν _f is the Poisson's ratio of the anisotropic filler contained in the target material,
ν _m is the Poisson's ratio of the resin composition constituting the target material,
_Vf is the volume content of the anisotropic filler contained in the target material with respect to the target material. )
(In Formula (2),
α _ν is the volume expansion coefficient of the target material,
α ₁₁ is the linear expansion coefficient in the principal axis direction of the target material,
α ₂₂ is a linear expansion coefficient in a direction orthogonal to the material principal axis direction,
α ₃₃ is a linear expansion coefficient in a direction orthogonal to the material principal axis direction,
α ₂₂ ≧ α ₃₃ )

[2] Further, the present invention provides the structural analysis model creation method according to [1], in which the α ₁₁ is obtained by the following equation (3) and the α ₂₂ is obtained by the following equation (4). is there.
(In Formula (3),
α _s is the linear expansion coefficient in the principal axis direction of a virtual anisotropic resin molded body that is completely oriented, that is, a virtual anisotropic resin molded body having an orientation degree of 1,
α _r is the linear expansion coefficient in the principal axis direction of the virtual anisotropic resin molded body in which the orientation is random, and is defined by the above formula (3) ′,
λ ₁₁ is the degree of orientation of the target material in the principal axis direction,
x is a parameter indicating the degree of nonlinearity in each direction with respect to the degree of orientation. )
(In Formula (3) ′, α _ν is the volume expansion coefficient of the target material.)
(In Formula (4),
α _b is a linear expansion coefficient in a direction orthogonal to the principal axis direction of a virtual anisotropic resin molded body that is completely oriented, that is, a virtual anisotropic resin molded body having an orientation degree of 1,
α _r is the linear expansion coefficient in the principal axis direction of the virtual anisotropic resin molded body in which the orientation is random, and is defined by the following formula (4) ′.
λ ₁₁ is the degree of orientation of the target material in the principal axis direction,
x is a parameter indicating the degree of nonlinearity in each direction with respect to the degree of orientation. )
(In Expression (4) ′, α _ν is the volume expansion coefficient of the target material.)

[3] The present invention is the structural analysis model creation method according to [1] or [2], wherein the physical property information further includes an elastic modulus defined by the following formula (5).
(In Formula (5),
E ₁₁ is a principal axis direction of the elastic modulus of the target material,

_Em is the elastic modulus of the molded body made of the resin composition constituting the target material,
ζ ₁₁ is a value represented by equation (5) ′,
η is a value represented by the formula (5) '',
_Vf is the volume content of the anisotropic filler with respect to the target material. )
(In Formula (5) ′,
λ ₁₁ is the degree of orientation of the target material in the principal axis direction,
L / d is the aspect ratio of the anisotropic filler contained in the target material. )
(In Formula (5) ″, E _f is the elastic modulus of the filler contained in the target material, and E _m is the elastic modulus of the molded body made of the resin composition constituting the target material.)

[4] Further, the present invention is the structural analysis model creation method according to any one of [1] to [3], wherein the physical property information further includes a shear elastic modulus defined by Expression (6). .
(In Formula (6),
G ₁₂ is a principal axis direction of the shear modulus of the target material,
G _m is the shear modulus of the molded body made of the resin composition constituting the target material,
ζ ₁₁ is a value represented by equation (6) ′,
η _g is a value represented by equation (6) ″,
_Vf is the volume content of the anisotropic filler with respect to the target material. )
(In Formula (6) ′,
λ ₁₁ is the degree of orientation of the target material in the principal axis direction,
L / d is the aspect ratio of the anisotropic filler contained in the target material. )
(In Formula (6) ″, G _f is the shear elastic modulus of the filler contained in the target material, and G _m is the shear elastic modulus of the molded body made of the resin composition constituting the target material.)

[5] Further, according to the present invention, the structural analysis model according to any one of [1] to [4], wherein α _ν is an actual measurement value of a volume expansion coefficient of a target material obtained by PVT resin characteristic analysis. It is a creation method.

  [6] Further, the present invention is the structural analysis model creation method according to any one of [1] to [5], wherein the anisotropic resin molded body includes a weld portion.

  [7] In the present invention, the search step includes a first centroid position deriving step of deriving a first centroid position for each of the plurality of first areas, and a second centroid of each of the plurality of second areas. A second centroid position deriving step for deriving a position; and for each of the second regions, a first centroid position closest to the second centroid position is searched for, and the first centroid position closest to the second centroid position is provided. The method for creating a model for structural analysis according to any one of [1] to [6], further including a shortest first region setting step in which a first region is a first region closest to the second region.

  [8] Further, in the present invention, the orientation state is an orientation degree in a principal axis direction of a target material, and the plurality of types of virtual molded bodies are in a range of 10 types to 1000 types according to the range of the orientation degree. The structural analysis model creation method according to any one of [1] to [7], set in step (1).

  [9] Further, the present invention is a structural analysis method for structurally analyzing an anisotropic resin molded body using the structural analysis model creation method according to any one of [1] to [8], First to obtain first region orientation information including orientation state calculation basic information for performing resin flow analysis on the anisotropic resin molded body and calculating an orientation state in each of the plurality of first regions divided into elements Separately from the region orientation information acquisition step and the first region orientation information acquisition step, the anisotropic resin molded body is divided into the plurality of second regions for structural analysis of the anisotropic resin molded body. A structural analysis element splitting step; and a structural analysis step of performing structural analysis of the anisotropic resin molding based on structural analysis model information including information on the second region orientation state and the second region physical property information. , The first region The orientation state calculating step is a step of calculating a first region orientation state corresponding to an orientation state in each of the plurality of first regions based on the orientation state calculation basic information, and the searching step is for the structural analysis It is a structure analysis method of an anisotropic resin molding which is a step which searches for the 1st field nearest to the 2nd field about each of a plurality of 2nd fields in a model.

  [10] Further, the present invention is a warpage deformation prediction method for predicting warpage deformation of the anisotropic resin injection molded body by using the structural analysis method according to [9].

[11] The present invention is a structural analysis model creation program for causing a computer to create a structural analysis model for structural analysis of an anisotropic resin molded body containing an anisotropic filler. For each of a plurality of types of virtual molded bodies having different orientation states from the anisotropic resin molded body, the Poisson's ratio defined by the following formula (1) and the linear expansion coefficient defined by the following formula (2) A physical property information creating step for creating physical property information including at least one, and a first region for calculating a first region orientation state in each of the plurality of first regions divided into elements in the resin flow analysis of the anisotropic resin molded body An orientation state calculation step, and a search step for searching for the first region closest to the second region for each of the plurality of second regions divided for the structural analysis of the anisotropic resin molding, A second region orientation state setting step in which the first region orientation state in the first region determined to be the closest in the search step is set to each second region orientation state, and the physical property information is referred to, the plurality of second regions A program for causing a computer to execute a second region physical property information setting step in which each physical property information corresponding to the second region orientation state is used as each second region physical property information for each region.
(In Formula (1),
ν ₂₃ is a Poisson's ratio to two directions orthogonal to the principal axis direction of the target material,
ν _f is the Poisson's ratio of the anisotropic filler contained in the target material,
ν _m is the Poisson's ratio of the resin composition constituting the target material. )
(In Formula (2),
α _ν is the volume expansion coefficient of the target material,
α ₁₁ is the linear expansion coefficient in the principal axis direction of the target material,
α ₂₂ is a linear expansion coefficient in a direction orthogonal to the material principal axis direction,
α ₃₃ is a linear expansion coefficient in a direction orthogonal to the material principal axis direction,
α ₂₂ ≧ α ₃₃ )

  According to the present invention, in analyzing the structure of an anisotropic resin molded body, highly accurate results can be obtained in a short time, and the load on the computer can be greatly reduced.

It is a flowchart which shows an example of the structural analysis method of the anisotropic resin molding which concerns on this invention. It is a flowchart which shows an example of the model creation model for structural analysis in FIG. It is a flowchart which shows an example of the search step in FIG. An example of hardware resources for realizing a series of analysis programs is shown. The shape of the test piece A after element division in the resin flow analysis is shown. The orientation direction distribution of the anisotropic filler in the test piece A obtained by resin flow analysis is shown. The shape of the test piece A after dividing into elements for structural analysis is shown. An example is shown in which the orientation direction distribution of the anisotropic filler in the test piece A obtained by the flow analysis is applied to a finite element division model used for structural analysis. The principal stress distribution of the test piece A obtained when the test piece A is structurally analyzed by the method of Example 1 is shown. The principal stress distribution of the test piece A obtained when the structure of the test piece B is analyzed by the method of Comparative Example 1 is shown. The shape of the test piece B after the element division in the resin flow analysis is shown. The orientation direction distribution of the anisotropic filler in the test piece B obtained by resin flow analysis is shown. The shape of the test piece B after dividing into elements for structural analysis is shown. The relationship between the temperature and specific volume of reinforced polybutylene terephthalate containing 30% by weight of glass fiber is shown. An example when the orientation direction distribution of the anisotropic filler in the specimen B obtained by the flow analysis is applied to a finite element division model used for structural analysis is shown. The position which evaluated the linear expansion coefficient in the test piece B is shown. The principal stress distribution of the test piece A obtained when the structure of the test piece B is analyzed by the method of Example 2 is shown. The principal stress distribution of the test piece A obtained when the structure of the test piece B is analyzed by the method of Comparative Example 2 is shown.

  Hereinafter, specific embodiments of the present invention will be described in detail. However, the present invention is not limited to the following embodiments, and may be implemented with appropriate modifications within the scope of the object of the present invention. can do.

<Structural analysis method of anisotropic resin molding>
FIG. 1 is a flowchart showing an example of a method for analyzing the structure of an anisotropic resin molded body according to the present invention. The structural analysis method according to the present invention is an orientation state for performing resin flow analysis on an anisotropic resin molded body containing an anisotropic filler and calculating an orientation state in each of a plurality of first regions divided into elements. Separately from the first region orientation information acquisition step (S1) for acquiring the first region orientation information including the calculation basic information and the first region orientation information acquisition step (S1), for structural analysis of the anisotropic resin molded body A structural analysis element dividing step (S2) for dividing the anisotropic resin molded body into a plurality of second regions, and a structural analysis for creating a structural analysis model for performing the structural analysis of the anisotropic resin molded body Model creation step (S3) and a structural analysis step (S4) for performing a structural analysis of the anisotropic resin molding based on the structural analysis model.

[First Region Orientation Information Acquisition Step (S1)]
First, the first region orientation information acquisition step (S1) will be described. In obtaining the orientation state of the anisotropic resin molded body, it may be obtained by observing an actual anisotropic resin molded body obtained by processing such as injection molding using X-ray CT or the like. However, in the design stage, which is the pre-trial stage, it is difficult to actually measure the orientation state because an actual anisotropic resin molded body has not been prototyped. Therefore, it is preferable to simulate the orientation state in each of the plurality of first regions divided into elements by using resin flow analysis.

  The method of resin flow analysis is not particularly limited, and examples thereof include resin flow analysis using a finite element method. When obtaining the orientation state by resin flow analysis using the finite element method, it is preferable to divide the element into five or more parts in the thickness direction of the anisotropic resin molded body. When the degree of element division is 4 or less, the fiber orientation in the thickness direction is averaged, and the accuracy of the structural analysis can be lowered, which is not preferable.

  The kind of orientation state is not specifically limited, For example, an orientation degree, an orientation direction, etc. are mentioned.

  Orientation degree information relating to the degree of orientation is usually expressed in a 3 × 3 matrix. From this matrix, eigenvectors and eigenvalues in three directions are obtained. The orientation direction is obtained from each eigenvector. Further, the degree of orientation can be obtained from each eigenvalue.

  Although it is not essential, if a two-point gate is used during injection molding, welds occur and stress concentrates on the welds. Therefore, it is preferable to model the anisotropic resin molded product as including the welds. .

  In addition, it is preferable to model the gate, runner, insert portion and the like in order to increase the accuracy of the flow analysis and, consequently, the final accuracy of the structural analysis.

[Structural analysis element division step (S2)]
Next, the structural analysis element dividing step (S2) will be described. The structural analysis element dividing step (S2) separates the anisotropic resin molded body into a plurality of second regions for structural analysis of the anisotropic resin molded body separately from the first region orientation information acquisition step (S1). This is the step of element division.

  The element division step for structural analysis (S2) is not particularly limited as long as it performs element division separately from the element division in the first region orientation information acquisition step (S1), and the first region orientation information is obtained. It may be performed after step (S1), may be performed before the first region orientation information acquisition step (S1), or may be performed in parallel with the first region orientation information acquisition step (S1).

  The method for dividing the element is not particularly limited. As an example, first, import the shape of the anisotropic resin molding into a personal computer using a CAD interface or the like, or create the shape of the anisotropic resin molding using a CAD system and set the modeling range. To do. Next, the anisotropic resin molded body is divided into a plurality of regions by performing element division such as a finite element method with an element division preprocessor or the like.

  The shape of the element is not particularly limited, and a tetrahedral primary element, a tetrahedral secondary element, a hexahedral primary element, a hexahedral secondary element, and the like can be selected. What is necessary is just to select suitably according to the specification of a system, calculation cost, etc.

  The number of elements is not particularly limited, and may be appropriately selected in consideration of calculation accuracy, calculation time, and the like.

[Structural analysis element division step (S3)]
The structural analysis model creation step (S3) is a step of creating a structural analysis model as a preliminary stage for performing the structural analysis of the anisotropic resin molded body. Hereinafter, the structural analysis model creation step (S3) will be described in detail with reference to FIG.

  FIG. 2 is a flowchart showing an example of the structural analysis model creation step (S3) in FIG. 1, and corresponds to the structural analysis model creation method according to the present invention. The structural analysis model creation step (S3) includes physical property information including at least one of a predetermined Poisson's ratio and a predetermined linear expansion coefficient for each of a plurality of types of virtual molded bodies having different orientation states from the anisotropic resin molded body. A physical property information creation step (S31) to be created, and a first region orientation state calculation step for calculating a first region orientation state in each of the plurality of first regions divided in the resin flow analysis of the anisotropic resin molded body ( S32), a search step (S33) for searching for the first region closest to the second region for each of the plurality of second regions divided for the structural analysis of the anisotropic resin molded body, and a search step A second region alignment state setting step (S34) in which the first region alignment state in the first region determined to be the closest in (S33) is set to each second region alignment state, and see physical property information , For each of a plurality of second regions, and a second region property information setting step (S35) and include that the property information corresponding to the second region alignment state and the second region property information for each.

[Physical property information creation step (S31)]
The physical property information creation step (S31) is a step of creating physical property information (mechanical property values) corresponding to the orientation state exemplified by the degree of orientation in advance, specifically, the anisotropic resin molded body and the orientation state. Although different from each other, each of a plurality of types of virtual molded bodies having the same other conditions such as the composition of the resin constituting the anisotropic resin molded body and the shape of the resin molded body is defined by the following formula (1). This is a step of creating physical property information including at least one of the Poisson's ratio and the linear expansion coefficient defined by the following formula (2). By creating the physical property information (mechanical property value) in advance, there is an advantage that confirmation and correction of the input physical property information can be facilitated.

  The creation information of the physical property information is not particularly limited, and when creating the physical property information constituting the rigidity matrix, it is possible to use a model that is partially improved based on the Halpin-Tsai equation. The created physical property information is stored in a hardware recording area exemplified by a hard disk, a memory, and the like.

In Expression (1), ν ₂₃ is a Poisson's ratio with respect to two directions orthogonal to the principal axis direction of the target material. Further, ν _f is the Poisson's ratio of the anisotropic filler contained in the target material. Further, ν _m is the Poisson's ratio of the resin composition that constitutes the target material. V _f is the volume content of the anisotropic filler contained in the target material with respect to the target material.

  In the present invention, the method for obtaining the Poisson's ratio is improved. In general, as a method for experimentally determining the Poisson's ratio in each direction, the method described in Japanese Patent Application Laid-Open No. 2009-128033 is known, but it is extremely difficult to determine the orientation degree dependency. Actually, it is obtained by model calculation. At that time, depending on the calculation method, it tends to be an actually impossible value. In this case, the positive definite value becomes negative, and calculation itself may be impossible depending on the structure analysis software. By adopting the formula (1), a positive definite value in the material physical property is likely to be positive, and thus there is an advantage that the general-purpose structural calculation software can be easily handled.

In the formula (2), α _ν is a volume expansion coefficient of the target material. Α ₁₁ is the linear expansion coefficient of the target material in the principal axis direction. Further, α _{22 and} α ₃₃ are linear expansion coefficients in a direction orthogonal to the material principal axis direction, and have a relationship of α22 ≧ α33.

  In formula (2), the relationship between the linear expansion coefficient and the volume expansion coefficient in each direction was used for the calculation. If this relationship is not satisfied, the volume expansion / contraction amount contradicts the expansion / contraction amount in each direction. . However, the general Schapery model, the Hashin-Shtrikman model, etc. do not always satisfy this relationship in the method for obtaining the linear expansion coefficient, which has caused the accuracy of the calculation result to be poor.

α ₁₁ is preferably obtained by the following equation (3), and α ₂₂ is preferably obtained by the following equation (4).

In the formula (3), α _s is a linear expansion coefficient in the main axis direction of a virtual anisotropic resin molded body having a complete orientation, that is, a virtual anisotropic resin molded body having an orientation degree of 1. Further, α _r is a linear expansion coefficient in the principal axis direction of the virtual anisotropic resin molded body in which the orientation is random, and is defined by the expression (3) ′. λ ₁₁ is the degree of orientation of the target material in the principal axis direction. x is a parameter indicating the degree of nonlinearity in each direction with respect to the degree of orientation.

Moreover, in Formula (3) ', (alpha) _v is the volume expansion coefficient of object material.

In the formula (4), α _b is a linear expansion coefficient in a direction perpendicular to the principal axis direction of a virtual anisotropic resin molded body having a complete orientation, that is, a virtual anisotropic resin molded body having an orientation degree of 1. α _r is a linear expansion coefficient in the principal axis direction of the virtual anisotropic resin molded body in which the orientation is in a random state, and is defined by Expression (4) ′. λ ₁₁ is the degree of orientation of the target material in the principal axis direction. x is a parameter indicating the degree of nonlinearity in each direction with respect to the degree of orientation.

Further, in the formula (4) ′, α _ν is a volume expansion coefficient of the target material.

In the embodiment of the present invention, material property values are not continuously given according to the degree of orientation of the virtual anisotropic resin molded body in the principal axis direction, but are handled piecewise. Specifically, since the orientation degree λ _{11 in} the principal axis direction of the virtual anisotropic resin molded body can take a value of 1/3 or more and 1 or less, for example, λ ₁₁ is in the range of 0.9 or more and 1 or less. The case shall have one kind of material property. As for this section, the smaller the section, the more continuous and the higher the analysis accuracy. However, since the number of sections increases, confirmation / change of material properties becomes difficult. In consideration of both analysis accuracy in structural analysis and ease of confirmation / change of material properties, a plurality of types of virtual molded bodies are set in a range of 10 to 1000 types according to the range of the degree of orientation. It is preferable.

  As physical property information, in addition to the Poisson's ratio represented by the above formula (1) and the linear expansion coefficient represented by the formula (2), the elastic modulus defined by the following formula (5), defined by the following formula (6) Shear modulus of elasticity and the like.

In Equation (5), E ₁₁ is the elastic modulus in the principal axis direction of the target material, E _m is the elastic modulus of the molded body made of the resin composition constituting the target material, and ζ ₁₁ is ) ′, Η is a value represented by the formula (5) ″, and V _f is a volume content of the anisotropic filler with respect to the target material.

In Formula (5) ′, λ ₁₁ is the degree of orientation in the principal axis direction of the target material, and L / d is the aspect ratio of the anisotropic filler contained in the target material.

In Formula (5) ″, E _f is the elastic modulus of the filler contained in the target material, and E _m is the elastic modulus of the molded body made of the resin composition constituting the target material.

In the formula (6), G ₁₂ is the main axis of the shear modulus of the subject material, G _m is the shear modulus of the molded article comprising the resin composition constituting the target material, zeta ₁₁ has the formula (6) ′, η _g is a value represented by the formula (6) ″, and V _f is a volume content of the anisotropic filler with respect to the target material.

In Equation (6) ′, λ ₁₁ is the degree of orientation of the target material in the principal axis direction, and L / d is the aspect ratio of the anisotropic filler contained in the target material.

In Formula (6) ″, G _f is the shear modulus of the filler contained in the target material, and G _m is the shear modulus of the molded body made of the resin composition that constitutes the target material.

  With respect to the formula (5) ′ and the formula (6) ′, a technique for obtaining the elastic modulus and the like based on the orientation state exemplified by the orientation degree and the like has been known as a homogenization method. However, since the homogenization method has a complicated calculation method for obtaining the elastic modulus and the like, a calculation error, a program coding error, and the like may occur. In order to facilitate the calculation procedure and the creation of the calculation program, it is preferable to employ the above formulas (5) ′ and (6) ′ as a model that can take into account the degree of orientation in a microscopic region.

  In the formula (5) ′ and the formula (6) ′, it is preferable to adopt an actual measurement value for the aspect ratio L / d. However, the aspect ratio of the filler contained in the anisotropic resin molded body depends on the damage of the filler and the like. It has a distribution and changes from the state before injection molding. Therefore, the aspect ratio may be treated as a parameter for the operation of the expressions (5) ′ and (6) ′.

  In addition, as physical property information, the following formulas (7) to (9) can be exemplified.

In Expression (7), ν ₁₂ and ν ₁₃ are Poisson's ratios to the direction orthogonal to the principal axis direction of the target material. V _f is the volume content of the anisotropic filler with respect to the target material, and ν _f is the Poisson's ratio of the anisotropic filler included in the target material. V _m is the volume content of the resin composition with respect to the target material, and ν _m is the Poisson's ratio of the resin composition constituting the target material.

In Formula (8-1) and Formula (8-2), E ₂₂ and E ₃₃ are elastic moduli in a direction orthogonal to the principal axis direction of the target material, and E _m is from the resin composition constituting the target material. It is an elastic modulus of a formed product. Also, ζ ₂₂ is 2, η is a value represented by the above formula (5) ″, and V _f is a volume content of the anisotropic filler with respect to the target material.

In the formula (9-1) and formula (9-2), G ₂₃ is the direction of the shear modulus perpendicular to the main axis direction of the subject material, G _m is made of a resin composition constituting the target material forming This is the shear modulus of the body. Further, ζ ₂₂ is 2, η _g is a value represented by the formula (6) ″, and V _f is a volume content of the anisotropic filler with respect to the target material.

[First Region Orientation State Calculation Step (S32)]
A 1st area | region orientation state calculation step (S32) is a step which calculates the 1st area | region orientation state in each of the some 1st area | region divided into the element in the resin flow analysis of an anisotropic resin molding. In the present embodiment, the alignment state is described as the alignment direction and the alignment degree, but is not limited thereto.

  Orientation degree information relating to the degree of orientation is usually expressed in a 3 × 3 matrix. From this matrix, eigenvectors and eigenvalues are obtained. The orientation direction is obtained from the eigenvector. The degree of orientation can be obtained from the eigenvalue. The three directions of the eigenvector are converted into a coordinate system and given to each element in a later step. Further, the physical property information obtained in the physical property information creation step (S31) is set in each element in a later step according to the degree of orientation. Since these settings can be easily confirmed before and after the structural analysis, there is an effect that it is easy to confirm a defect in the physical property setting.

[Search Step (S33)]
The searching step (S33) is a step of searching for the first region closest to the second region for each of the plurality of second regions divided into elements for the structural analysis of the anisotropic resin molded body. The mode of the search step (S33) is not particularly limited. For example, as shown in FIG. 3, a first centroid position deriving step (S331) for deriving a first centroid position for each of the plurality of first regions; A second barycentric position deriving step (S332) for deriving a second barycentric position for each of the plurality of second areas, and searching for the first barycentric position closest to the second barycentric position for each of the second areas, A mode having a shortest first region setting step (S333) in which the first region having the first centroid position closest to the two centroid positions is set as the first region closest to the second region can be exemplified.

(First gravity center derivation step (S331))
A specific example of the search step (S33) will be described with reference to FIG. In the first center-of-gravity position deriving step (S331), the first center-of-gravity position is obtained from the position of the node constituting the element for each element related to the first region used in the calculation in the first region orientation information acquisition step (S1). . The first centroid position is compared with the second centroid position of the structural analysis model to be calculated next.

(Second center of gravity position deriving step (S332))
In the second centroid position deriving step (S332), as in the first centroid position deriving step (S331), the nodes constituting the elements for each element of the second region divided in the structural analysis element dividing step (S2) The position of the center of gravity is obtained from the position of.

(Shortest first area setting step (S333))
In the shortest first region setting step (S333), for each element of the second region divided in the structural analysis element dividing step (S2), the second center of gravity position obtained in the second center of gravity position deriving step (S332), The distance from the first centroid position obtained in the first centroid position derivation step (S331) is calculated, and the closest position information is applied. When the shortest first region setting step (S333) is completed, the process proceeds to the second region orientation state setting step (S34) in FIG.

[Second Region Orientation State Setting Step (S34)]
Returning to FIG. The second region orientation state setting step (S34) is a step in which the first region orientation state in the first region determined to be the closest in the search step (S33) is set as each second region orientation state. In the case of this method, the element division in the resin flow analysis used in the first region orientation information acquisition step (S1) is rough, and the second of the elements in the second region divided in the structural analysis element division step (S2). Even when there is no element related to the first region used in the first region orientation information acquisition step (S1) at the center of gravity position, the information on the first region orientation state can be applied to the second region.

  According to this method, it is possible to apply the information on the first region orientation state to the second region regardless of the element division state. This fitting can be obtained by interpolating the orientation state information using the shape function of the element, but if the element division is different, there may be a situation where the interpolation cannot be performed, so the orientation state information is obtained. In some cases, restrictions are imposed on the information density and the element division of the structural analysis, and the restrictions may affect the accuracy of the structural analysis itself.

[Second Area Physical Property Information Setting Step (S35)]
The second area physical property information setting step (S35) refers to the physical property information created in advance in the physical property information creating step (S31) and stored in the recording area of the hardware. In this step, the physical property information corresponding to the two-region orientation state is used as the second region physical property information. In this step, the coordinate system based on the information on the second region orientation state obtained in the second region orientation state setting step (S34) and the material properties corresponding to the degree of orientation are set as structural analysis elements. When this second region property information setting step (S35) is completed, the structure analysis can be executed, and the process proceeds to the structure analysis step (S4) in FIG.

[Structural analysis step (S4)]
The structural analysis step (S4) is a step of performing structural analysis of the anisotropic resin molded body based on the structural analysis model. In this step, a structure including information on the second region orientation state set in the second region orientation state setting step (S34) and second region property information set in the second region property information setting step (S35) Using the analysis model information, the structural analysis software program is executed to execute the structural analysis. Various conditions such as constraint conditions, load conditions, and calculation conditions in the structural analysis may be in accordance with the specifications of the structural analysis software, and are not particularly limited.

  Structural analysis can simulate deformation, stress, and the like of the anisotropic resin molded body. From this simulation result, it is possible to determine the presence or absence of breakage in the anisotropic resin molding.

  As a specific aspect, warp deformation of an anisotropic resin injection molded body can be analyzed by setting a load condition in a structural analysis to a temperature load and a residual stress load. Therefore, the warp deformation of the anisotropic resin injection molded body can be predicted by using the structural analysis method according to the present invention.

<Analysis program>
A series of analysis programs including the structural analysis model creation program according to the present invention is realized by cooperation of software and hardware resources.

  FIG. 4 shows an example of a hardware resource H for realizing a series of analysis programs. The hardware resource H includes an information processing apparatus 1, an input apparatus 2 that receives various requests from a designer, and an output apparatus 3 that outputs an analysis result performed by the information processing apparatus 1. The information processing apparatus 1 is connected to the CAD apparatus 4 via a network NW such as a LAN (Local Area Network).

  The information processing device 1 includes a CPU (Central Processing Unit) 10, a RAM (Random Access Memory), etc., a main storage device 20, an I / O interface 30, an auxiliary storage device 40 constituted by a hard disk, etc. And a network interface 50 that controls data exchange with a device connected to the network NW.

  The auxiliary storage device 40 stores an analysis program 41 for causing the hardware device 1 to execute the series of steps described above. The analysis program 41 includes a first region orientation information acquisition program 41A for causing the hardware device 1 to execute a first region orientation information acquisition step (S1) by flow analysis, and a structure analysis element division step (S2). Structural analysis element division program 41B for causing the apparatus 1 to execute, structural analysis model creation program 41C for causing the hardware apparatus 1 to execute the structural analysis model creation step (S3), and structural analysis step (S4) And a structural analysis program 41D for causing the hardware device 1 to execute the above. The structural analysis method according to the present invention is realized by the CPU 10 loading the analysis program 41 stored in the auxiliary storage device 40 into the main storage device 20 and executing it. In particular, in the anisotropic resin molded body structural analysis model creation method according to the present invention, the CPU 10 loads the structural analysis model creation program 41C stored in the auxiliary storage device 40 to the main storage device 20 and executes it. Is realized.

  In the above embodiment, each step from the first region orientation information acquisition step (S1) to the structure analysis step (S4) is performed by executing a combination of a plurality of programs. The program may be a program built as a whole from the beginning, and the form, scale, installation location, etc. of the computer to be executed are not limited.

  EXAMPLES Hereinafter, although an Example and a comparative example are shown and this invention is demonstrated concretely, this invention is not limited to these Examples.

<Test Example 1>
In Test Example 1, the load resistance at the time of the three-point bending test was examined using a test piece A obtained by injection molding having a shape in which a rib is provided on the long piece side of the flat plate portion.

[Example 1]
First, the test piece A was subjected to a bending test based on JIS178. The interval between the fixed jigs was set to 62 mm, an indenter was placed in the center of the test piece, the indenter was moved at a strain rate of 1% / min, and the load was measured with a load sensor mounted on the test equipment. As a result, the test piece was damaged at a load of 272N.

[First Region Orientation Information Acquisition Step (S1)]
Then, the shape of the test piece A was created as CAD data, and a runner, a gate, etc. were added and the finite element division | segmentation was performed. For the resin flow analysis, AUTODESK SIMULATION MODFLOW INSIGHT was used. The conditions for resin flow analysis are shown below.
Resin: Reinforced polybutylene terephthalate containing 30% by weight of glass fiber Resin temperature: 260 ° C
Mold temperature: 80 ℃
Injection flow rate: 57.7 cm ³ / s
Holding pressure: 70 MPa
Holding time: 15 seconds Cooling time: 10 seconds Element: Tetrahedral primary element (number of divisions 450538)

  FIG. 5 shows the shape of the test piece A after element division in the resin flow analysis, and FIG. 6 shows the orientation direction distribution of the anisotropic filler in the test piece A obtained by the resin flow analysis.

[Structural analysis element division step (S2)]
Separately from the finite element model for flow analysis, the specimen A was divided into a plurality of second regions for structural analysis of the anisotropic resin molding. The shape of the test piece A after element division is shown in FIG. In FIG. 7, the runner, the gate, and the like are omitted, and a test piece support and an indenter are added instead. For the element, a tetrahedral secondary element was used, and the number of divisions was 66565. The test piece support was completely restrained, and a load 272N, which is the load resistance of the test piece, was given to the indenter, and the maximum principal stress and the generation position were examined.

[Structural analysis model creation step (S3)]
First, in the physical property information creation step (S31), physical property information is created according to the above formulas (1) to (9) for each of a plurality of types of virtual molded bodies having different orientation states from the test piece A, and a predetermined computer Stored in storage area. The physical property information obtained according to the formulas (1) to (9) was as follows.
Elastic modulus E _{f of} filler contained in target material: 72000 MPa
Elastic modulus E _{m of a} molded body made of a resin composition constituting the target material: 2500 MPa
Volume content V _{f with} respect to target material of anisotropic filler V _f : 18 vol%
Volume content V _{m with} respect to target material of resin composition: 82 vol%
Shear elastic modulus G _{f of the} filler contained in the target material: 29500 MPa
Shear modulus of the molded article comprising the resin composition constituting the target material _G m: 926MPa
Poisson's ratio ν _{f of} anisotropic filler contained in target material: 0.22
Poisson's ratio ν _m of the resin composition constituting the target material: 0.35
Aspect ratio L / d of anisotropic filler contained in target material: 14

  Then, a first region orientation state calculating step (S32) for calculating a first region orientation state in each of the plurality of first regions divided into elements in the resin flow analysis of the test piece A, and for the structural analysis of the test piece A A search step (S33) for searching for the first region closest to the second region for each of the plurality of second regions divided into elements in the first region in the first region determined to be closest in the search step (S33) The second region orientation state setting step (S34) in which the region orientation state is set to each second region orientation state, and the physical property information for each of the plurality of second regions with reference to the physical property information. The structure analysis model is completed through a second region property information setting step (S35) in which each is a second region property information. The completed structural analysis model is shown in FIG. Compared to the original distribution state of FIG. 6, the positions indicating the directions are different, but the directions of the two coincide.

[Structural analysis step (S4)]
Calculation is performed using the obtained model for structural analysis. For the calculation, structural analysis software NX I-DEAS (manufactured by Siemens PLM Software) was used.

  FIG. 9 shows a principal stress distribution obtained by structural analysis of the test piece A by the method of the first embodiment. The position showing the maximum stress is slightly shifted from the center part of the test piece. Since the two-point gate was used at the time of molding the test piece, it was considered that weld was generated and stress was concentrated on the weld. Moreover, the maximum principal stress in Example 1 was 139 MPa (Table 1).

  Moreover, the measured value of the tensile fracture strength of the test piece was 140 MPa, and the weld portion was the starting point of fracture (Table 1).

  Thus, the result of Example 1 is very close to the actual measurement, and the accuracy of the present invention is very good in Example 1.

[Comparative Example 1]
For the same test piece as in Example 1, 6320 MPa, which is a value obtained by averaging the elastic modulus in the flow direction and the elastic modulus in the flow vertical direction, is used as the elastic modulus. The structural analysis of the specimen A was performed in the same manner as in Example 1 except that.

  FIG. 10 shows a principal stress distribution obtained by structural analysis of the test piece A by the method of Comparative Example 1. The position showing the maximum stress was the center of the test piece, which was a position different from the weld generation position. The maximum principal stress in Comparative Example 1 was 106 MPa (Table 1).

  As described above, the result of Comparative Example 1 is different from the actual measurement, and it can be said that the accuracy of Comparative Example 1 is inferior to that of Example 1.

<Test Example 2>
In Test Example 2, a thermal expansion behavior was examined using a test piece B obtained by injection molding, which is a shape in which a flat plate portion is provided with a box shape, a rib, a cylinder, and the like.

[Example 2]
First, the test piece B obtained by injection molding was cut out, annealed in a thermostatic chamber at 120 ° C. for about 2 hours, and then allowed to cool. Then, it installed in the linear expansion coefficient measuring apparatus, temperature was raised from -30 degreeC to room temperature, and the linear expansion coefficient was calculated | required from the dimensional change amount in four places. The results are summarized in Table 2 below.

[Fiber orientation information acquisition step (S1)]
Subsequently, the finite element division of the test piece B was performed in the same manner as in Example 1. The conditions for resin flow analysis are shown below.
Resin: Reinforced polybutylene terephthalate containing 30% by weight of glass fiber Resin temperature: 250 ° C
Mold temperature: 60 ℃
Injection flow rate: 32 cm ³ / s
Holding pressure: 78.4 MPa
Holding pressure time: 10 seconds Cooling time: 15 seconds Element: Tetrahedral primary element (number of divisions 486749)

  The shape of the test piece B after element division in the resin flow analysis is shown in FIG. 11, and the orientation direction distribution of the anisotropic filler in the test piece B obtained by the resin flow analysis is shown in FIG.

[Structural analysis element division step (S2)]
Separately from the finite element model for flow analysis, the test piece B was divided into a plurality of second regions for structural analysis of the anisotropic resin molding. The shape of the test piece B after element division is shown in FIG. In FIG. 13, the cut shape is used. For the element, a tetrahedral secondary element was used, and the number of divisions was 66565. XYZ direction, YZ direction, and Z direction were constrained at three corners of the model, and a temperature load of 100 ° C. was applied.

[Material property setting step (S3)]
First, in the physical property information creation step (S31), physical property information is created according to the above formulas (1) to (9) for each of a plurality of types of virtual molded bodies having different orientation states from the test piece A, and a predetermined computer Stored in storage area. The physical property information obtained according to the formulas (1) to (9) was as follows.
The linear expansion coefficient α _{s in} the principal axis direction of a completely anisotropic virtual anisotropic resin molded body, that is, a virtual anisotropic resin molded body having an orientation degree of 1: 2.4 × 10 ⁻⁵ / ° C.
The linear expansion coefficient α _b in the direction orthogonal to the principal axis direction of the completely anisotropic virtual anisotropic resin molded body, ie, the virtual anisotropic resin molded body having an orientation degree of 1, is 6.2 × 10 ⁻⁵ / ° C.
Volume expansion coefficient α _v of target material: 1.53 × 10 ⁻⁴ / ° C.
Non-linearity x in each direction with respect to orientation degree: 15.5 (dimensionless)

The volume expansion coefficient α _ν was determined by determining the volume change of 30% by weight of glass fiber-reinforced reinforced polybutylene terephthalate (reinforced PBT) when the pressure was 0 MPa using a PVT measuring device. Is a value obtained by calculating.

  And like Example 1, 1st area | region orientation state calculation step (S32), search step (S33), 2nd area | region orientation state setting step (S34), 2nd area | region physical property information setting step (S35) After that, the structural analysis model is completed. The completed structural analysis model is shown in FIG. Compared to the original distribution situation of FIG. 12, the positions indicating the directions are different, but the directions of both are the same.

[Structural analysis step (S4)]
Subsequently, the structure of the test piece B was analyzed by the same method as in Example 1.

  From the dimensional changes at positions A, B, C, and D shown in FIG. 16, the linear expansion coefficient at each position was obtained. The results are shown in Table 2. FIG. 17 shows a deformation distribution obtained by structural analysis of the test piece B by the method of the second embodiment. From FIG. 17, it is suggested that the test piece B is deformed so that the flat portion is distorted.

  Further, in the test piece B, the linear expansion coefficients at the positions A, B, C, and D were measured. The results are shown in Table 2. During the actual measurement, it was visually confirmed that the flat portion of the test piece B was deformed so as to be distorted.

  Thus, it can be said that the change of the linear expansion coefficient difference between the positions A, B, C, and D in Example 2 reflects the actual measurement value, and it can be said that it reflects the distortion of the actual test piece B. . Therefore, the result of Example 2 is very close to the actual measurement, and it can be said that the accuracy of the present invention is very good in Example 2.

[Comparative Example 2]
For the same specimen B as in Example 2, 4.3 × 10 ⁻⁵ / ° C., which is a value obtained by averaging the linear expansion coefficient in the flow direction and the linear expansion coefficient in the flow vertical direction, was used as the linear expansion coefficient. Except that, the structural analysis of the test piece B was performed in the same manner as in Example 2.

  From the dimensional changes at positions A, B, C, and D shown in FIG. 16, the linear expansion coefficient at each position was obtained. The results are shown in Table 2. In the case of Comparative Example 2, there is no significant difference in linear expansion coefficient between the positions A, B, C, and D. FIG. 18 shows a deformation distribution obtained by structural analysis of the test piece B by the method of Comparative Example 2. From FIG. 18, it is clear that the flat portion of the test piece B remains flat and does not reflect the actual distortion of the test piece B.

  Thus, it can be said that the results of Comparative Example 2 are inferior to those of Example 2 in the structural analysis accuracy of Comparative Example 2 unlike actual measurement.

H Hardware resource 1 Information processing device 2 Input device 3 Output device 4 CAD device 10 CPU
20 Main storage device 30 I / O interface 40 Auxiliary storage device 41 Analysis program 41A First region orientation information acquisition program 41B Structural analysis element division program 41C Structural analysis model creation program 41D Structural analysis program 50 Network interface NW network S1 First Region orientation information acquisition step S2 Structural analysis element division step S3 Structural analysis model creation step S4 Structural analysis step S31 Physical property information creation step S32 First region orientation state calculation step S33 Search step S34 Second region orientation state setting step S35 Second Area physical property information setting step S331 First centroid position derivation step S332 Second centroid position derivation step S333 Shortest first area setting step

Claims (11)

A structural analysis model creation method for creating a structural analysis model for performing a structural analysis of an anisotropic resin molded article containing an anisotropic filler,
At least one of a Poisson's ratio defined by the following formula (1) and a linear expansion coefficient defined by the following formula (2) for each of a plurality of types of virtual molded bodies having different orientation states from the anisotropic resin molded body A physical property information creation step for creating physical property information including
A first region orientation state calculating step for calculating a first region orientation state in each of the plurality of first regions divided into elements in the resin flow analysis of the anisotropic resin molded body;
A search step of searching for the first region closest to the second region for each of the plurality of second regions divided into elements for the structural analysis of the anisotropic resin molded body,
A second region orientation state setting step in which the first region orientation state in the first region determined to be the closest in the search step is the second region orientation state;
Referring to the physical property information, and for each of the plurality of second regions, a second region physical property information setting step in which the physical property information according to the second region orientation state is used as each second region physical property information, Model creation method for structural analysis of anisotropic resin moldings.
(In Formula (1),
ν ₂₃ is a Poisson's ratio to two directions orthogonal to the principal axis direction of the target material,
ν _f is the Poisson's ratio of the anisotropic filler contained in the target material,
ν _m is the Poisson's ratio of the resin composition constituting the target material,
_Vf is the volume content of the anisotropic filler contained in the target material with respect to the target material. )
(In Formula (2),
α _ν is the volume expansion coefficient of the target material,
α ₁₁ is the linear expansion coefficient in the principal axis direction of the target material,
α ₂₂ is a linear expansion coefficient in a direction orthogonal to the material principal axis direction,
α ₃₃ is a linear expansion coefficient in a direction orthogonal to the material principal axis direction,
α ₂₂ ≧ α ₃₃ ) The structural analysis model creation method according to claim 1, wherein the α ₁₁ is obtained by the following equation (3), and the α ₂₂ is obtained by the following equation (4).
(In Formula (3),
α _s is the linear expansion coefficient in the principal axis direction of a virtual anisotropic resin molded body that is completely oriented, that is, a virtual anisotropic resin molded body having an orientation degree of 1,
α _r is the linear expansion coefficient in the principal axis direction of the virtual anisotropic resin molded body in which the orientation is random, and is defined by the above formula (3) ′,
λ ₁₁ is the degree of orientation of the target material in the principal axis direction,
x is a parameter indicating the degree of nonlinearity in each direction with respect to the degree of orientation. )
(In Formula (3) ′, α _ν is the volume expansion coefficient of the target material.)
(In Formula (4),
α _b is a linear expansion coefficient in a direction orthogonal to the principal axis direction of a virtual anisotropic resin molded body that is completely oriented, that is, a virtual anisotropic resin molded body having an orientation degree of 1,
α _r is the linear expansion coefficient in the principal axis direction of the virtual anisotropic resin molded body in which the orientation is random, and is defined by the following formula (4) ′.
λ ₁₁ is the degree of orientation of the target material in the principal axis direction,
x is a parameter indicating the degree of nonlinearity in each direction with respect to the degree of orientation. )
(In Expression (4) ′, α _ν is the volume expansion coefficient of the target material.) The model for structural analysis according to claim 1 or 2, wherein the physical property information further includes an elastic modulus defined by the following formula (5).
(In Formula (5),
E ₁₁ is a principal axis direction of the elastic modulus of the target material,
_Em is the elastic modulus of the molded body made of the resin composition constituting the target material,
ζ ₁₁ is a value represented by equation (5) ′,
η is a value represented by the formula (5) '',
_Vf is the volume content of the anisotropic filler with respect to the target material. )
(In Formula (5) ′,
λ ₁₁ is the degree of orientation of the target material in the principal axis direction,
L / d is the aspect ratio of the anisotropic filler contained in the target material. )
(In Formula (5) ″, E _f is the elastic modulus of the filler contained in the target material, and E _m is the elastic modulus of the molded body made of the resin composition constituting the target material.) The structural property model creation method according to any one of claims 1 to 3, wherein the physical property information further includes a shear elastic modulus defined by Expression (6).
(In Formula (6),
G ₁₂ is a principal axis direction of the shear modulus of the target material,
G _m is the shear modulus of the molded body made of the resin composition constituting the target material,
ζ ₁₁ is a value represented by equation (6) ′,
η _g is a value represented by equation (6) ″,
_Vf is the volume content of the anisotropic filler with respect to the target material. )
(In Formula (6) ′,
λ ₁₁ is the degree of orientation of the target material in the principal axis direction,
L / d is the aspect ratio of the anisotropic filler contained in the target material. )
(In Formula (6) ″, G _f is the shear elastic modulus of the filler contained in the target material, and G _m is the shear elastic modulus of the molded body made of the resin composition constituting the target material.) 5. The structural analysis model creation method according to claim 1, wherein α _ν is an actual measurement value of a volume expansion coefficient of a target material obtained by PVT resin characteristic analysis.   The method for creating a structural analysis model according to claim 1, wherein the anisotropic resin molded body includes a weld portion. The searching step includes
A first centroid position deriving step of deriving a first centroid position for each of the plurality of first regions;
A second centroid position deriving step of deriving a second centroid position for each of the plurality of second regions;
For each of the second regions, the first centroid position closest to the second centroid position is searched, and the first region having the first centroid position closest to the second centroid position is the closest to the second region. The structural analysis model creation method according to claim 1, further comprising a shortest first region setting step for one region. The orientation state is the degree of orientation in the principal axis direction of the target material,
The structural analysis model creation method according to any one of claims 1 to 7, wherein the plurality of types of virtual molded bodies are set in a range of 10 types to 1000 types according to the range of the degree of orientation. A structural analysis method for structural analysis of an anisotropic resin molded body using the structural analysis model creation method according to any one of claims 1 to 8,
Resin flow analysis is performed on the anisotropic resin molded body, and first region orientation information including orientation state calculation basic information for calculating an orientation state in each of the plurality of first regions divided into elements is obtained. One region orientation information acquisition step;
Apart from the first region orientation information acquisition step, a structural analysis element dividing step for dividing the anisotropic resin molded body into the plurality of second regions for structural analysis of the anisotropic resin molded body; ,
A structural analysis step for performing a structural analysis of the anisotropic resin molding based on structural analysis model information including the information on the second region orientation state and the second region physical property information,
The first region orientation state calculation step is a step of calculating a first region orientation state corresponding to an orientation state in each of the plurality of first regions based on the orientation state calculation basic information.
The method for analyzing the structure of an anisotropic resin molded body, wherein the searching step is a step of searching for the first region closest to the second region for each of the plurality of second regions in the structural analysis model.   A warpage deformation prediction method for predicting warpage deformation of the anisotropic resin injection molded body by using the structure analysis method according to claim 9. A structural analysis model creation program for causing a computer to create a structural analysis model for structural analysis of an anisotropic resin molded body containing an anisotropic filler,
At least one of a Poisson's ratio defined by the following formula (1) and a linear expansion coefficient defined by the following formula (2) for each of a plurality of types of virtual molded bodies having different orientation states from the anisotropic resin molded body A physical property information creation step for creating physical property information including
A first region orientation state calculating step for calculating a first region orientation state in each of the plurality of first regions divided into elements in the resin flow analysis of the anisotropic resin molded body;
A search step of searching for the first region closest to the second region for each of the plurality of second regions divided into elements for the structural analysis of the anisotropic resin molded body,
A second region orientation state setting step in which the first region orientation state in the first region determined to be the closest in the search step is the second region orientation state;
A second region property information setting step of referring to the physical property information and using the physical property information corresponding to the second region orientation state as the second region physical property information for each of the plurality of second regions; A program to be executed.
(In Formula (1),
ν ₂₃ is a Poisson's ratio to two directions orthogonal to the principal axis direction of the target material,
ν _f is the Poisson's ratio of the anisotropic filler contained in the target material,
ν _m is the Poisson's ratio of the resin composition constituting the target material. )
(In Formula (2),
α _ν is the volume expansion coefficient of the target material,
α ₁₁ is the linear expansion coefficient in the principal axis direction of the target material,
α ₂₂ is a linear expansion coefficient in a direction orthogonal to the material principal axis direction,
α ₃₃ is a linear expansion coefficient in a direction orthogonal to the material principal axis direction,
α ₂₂ ≧ α ₃₃ )