JP6420881B1 - Compression molding analysis system, compression molding analysis method, and compression molding analysis program - Google Patents

Abstract

To provide a compression molding analysis system, a compression molding analysis method, and a compression molding analysis program capable of performing appropriate compression molding analysis even when fiber reinforced resin is compression molded. A compression molding analysis system for a fiber reinforced resin in which fibers are mixed with a resin, the fiber model unit 11 modeling a fiber with a beam element, and a node not shared with the node of the beam element. Resin model unit 12 that models resin with solid elements, coupled analysis unit 13 that performs coupled analysis of beam element movement and solid element movement, and large deformation that performs large deformation analysis on beam and solid elements And an analysis unit 14. The coupling analysis unit 13 performs coupled analysis based on the movement of the beam element and the solid element during the large deformation of the beam element and the solid element by the large deformation analysis unit 14. . [Selection] Figure 2

Description

  The present invention relates to a compression molding analysis system, a compression molding analysis method, and a compression molding analysis program.

Countermeasures for global warming are an urgent issue, and automobile manufacturers and the like are required to significantly improve fuel efficiency by reducing the weight of automobiles in order to reduce CO ₂ emissions. Steel materials are used in conventional automobile parts and the like, but in recent years, fiber reinforced resin has been studied as an alternative material. Among them, the carbon fiber reinforced resin has a mechanical property equal to or higher than that of a steel material, but its weight can be suppressed to ¼ that of a steel material.

  The fiber reinforced resin is classified into a discontinuous short fiber reinforced resin having a fiber length of 2 mm or less, a discontinuous long fiber reinforced resin having a fiber length of 2 mm to 50 mm, and a continuous fiber reinforced resin having a fiber length of 50 mm or more. The fiber reinforced resin is known to have better mechanical properties as the fiber length is longer, and exhibits high strength when the fiber length is 10 mm or more and high toughness when the fiber length is 40 mm or more.

  Although the discontinuous short fiber reinforced resin can be molded into a part having a complicated shape such as a rib by injection molding, it is difficult to have strength characteristics equivalent to those of a steel material because the fiber length is short. Continuous fibers can be formed into high-strength and high-toughness parts by press molding, but complex shapes such as ribs cannot be formed. Among them, the discontinuous long fiber reinforced resin can be molded into a complicated shape by compression molding, and can be applied to various parts. In addition, when carbon fiber is used, it can have high strength and high toughness equivalent to or higher than that of steel, so it can be applied to car body structural parts in general, and it Weight reduction is achieved

  When the discontinuous long fiber reinforced resin is compression-molded, the fibers and the resin flow to be molded into a desired part shape. However, since fiber orientation (orientation) and fiber density occur during compression molding, it is known that regions having different rigidity and strength exist in the molded part. If there is a portion having low rigidity or strength, the performance of the parts is greatly affected, and therefore, a numerical simulation technique for predicting fiber orientation (direction) and fiber density is required.

  Conventionally, the fiber orientation was predicted by a fiber orientation algorithm such as the ARD (Anisotropic Rotary Diffusion) method based on the Folgar-Tucker Model used in the flow simulation of injection molding of discontinuous short fiber reinforced resin ( For example, Patent Document 1).

US Patent No. 8,571,828

  However, in the ARD method based on the Folgar-Tucker Model such as Patent Document 1, the fiber is treated as a rigid body that does not deform, and only the translational movement and rotation of the center of gravity of the fiber are calculated, and the deformation and stress of the fiber are calculated. Cannot be calculated. On the other hand, in actual compression molding, the fiber is deformed to generate stress, and when the fiber length exceeds 10 mm, deformation such as large bending or undulation occurs in the fiber and affects the molding behavior.

  The present invention has been made to solve the above-described problem, and is a compression molding analysis system capable of performing an appropriate compression molding analysis even when a fiber reinforced resin such as a discontinuous long fiber reinforced resin is compression molded. An object of the present invention is to provide a compression molding analysis method and a compression molding analysis program.

  In order to solve the above-described problems, one aspect of the compression molding analysis system of the present invention is a fiber reinforced resin compression molding analysis system in which fibers are mixed with a resin, and the fibers are modeled as fiber models using beam elements. A fiber model part to be converted, a resin model part that models the resin as a resin model with a solid element so as to have a node that is not shared with the node of the beam element, and the movement of the beam element and the movement of the solid element A coupled analysis unit that performs coupled analysis of the beam element and a large deformation analysis unit that performs large deformation analysis on the beam element and the solid element, and the coupled analysis unit includes the beam element by the large deformation analysis unit. The coupled analysis is performed based on the movement of the beam element during the large deformation of the solid element and the movement of the solid element.

  In this specification, “fiber reinforced resin” is a concept including a discontinuous long fiber reinforced resin, a discontinuous short fiber reinforced resin, a continuous fiber reinforced resin, and the like. In this specification, “resin” is a concept including a thermoplastic resin, a thermosetting resin, and the like.

  According to the present invention, in the fiber reinforced resin compression molding analysis system, the fiber model unit models the fiber with a beam element. The resin model unit models the resin with a solid element so as to have a node that is not shared with the node of the beam element. The coupled analysis unit performs coupled analysis between the movement of the beam element and the movement of the solid element. The large deformation analysis unit performs large deformation analysis on the beam element and the solid element. The coupled analysis by the coupled analysis unit is performed based on the movement of the beam element and the solid element during the large deformation of the beam element and the solid element by the large deformation analysis unit. Therefore, the behavior of the fiber and the resin, the orientation (orientation) of the fiber, and the density of the fiber are predicted with high accuracy.

  In another aspect of the compression molding analysis system of the present invention, the large deformation analysis unit may perform large deformation analysis on the solid element by an adaptive remeshing method.

  According to this aspect, since the large deformation analysis unit performs large deformation analysis on the solid element by the adaptive remeshing method, for example, even in a large deformation portion where a rib is formed, the accuracy is stable and accurate. Large deformation analysis is performed.

  According to another aspect of the compression molding analysis system of the present invention, the coupled analysis unit is configured to detect the movement of the beam element and the movement of the solid element during the large deformation of the beam element and the solid element by the large deformation analysis unit. For the beam element and the solid element based on, the relative displacement in the axial direction of the beam element is acquired as a slip amount, a resistance force is calculated from the acquired slip amount, and the large deformation analysis unit includes the beam element The resistance may be applied as a load to the nodes of the solid element.

  According to this aspect, the coupled analysis unit performs the beam element axial direction on the beam element and the solid element based on the movement of the beam element and the solid element during the large deformation of the beam element and the solid element by the large deformation analysis unit. The relative displacement of is acquired as the slip amount. The coupled analysis unit calculates a resistance force from the acquired slip amount, and applies the resistance force as a load to the nodes of the beam element and the solid element. Therefore, the phenomenon that the fiber moves in the resin when the fiber reinforced resin is deformed is accurately reproduced.

  According to another aspect of the compression molding analysis system of the present invention, the coupled analysis unit is configured to detect the movement of the beam element and the movement of the solid element during the large deformation of the beam element and the solid element by the large deformation analysis unit. For the beam element and the solid element based on the slip amount, which is a relative displacement in the axial direction of the beam element, and the temperature calculated at the node of the solid element, interpolation of the node of the beam element is calculated. The temperature at the position is acquired, and the resistance force is calculated from the acquired slip amount and the temperature. The large deformation analysis unit may apply the resistance force to the nodes of the beam element and the solid element as a load. Good.

  According to this aspect, the coupled analysis unit performs the beam element axial direction on the beam element and the solid element based on the movement of the beam element and the solid element during the large deformation of the beam element and the solid element by the large deformation analysis unit. Is obtained by interpolating from the temperature calculated at the node of the solid element and the temperature at the node position of the beam element calculated by interpolation. The coupled analysis unit calculates a resistance force from the acquired slip amount and temperature, and applies the resistance force as a load to the nodes of the beam element and the solid element. Therefore, even when the fiber reinforced resin is a thermoplastic resin, the fiber moves in the thermoplastic resin when the thermoplastic resin heated to a high temperature is deformed by compression molding while changing the temperature. Is accurately reproduced.

  In another aspect of the compression molding analysis system of the present invention, the coupled analysis unit is configured such that, in the direction perpendicular to the axial direction of the beam element, the movement of the nodal point of the beam element and the solid element are determined by a constraint method. In the axial direction of the beam element, the resistance force is calculated and the resistance force is applied to the node of the beam element and the node of the solid element as a load. Synthetic analysis may be performed.

  According to this aspect, the coupled analysis unit performs coupled analysis between the movement of the node of the beam element and the movement of the node of the solid element by the constraint method in the direction perpendicular to the axial direction of the beam element. In addition, the coupled analysis unit calculates a resistance force in the axial direction of the beam element, and applies the resistance force as a load to the node of the beam element and the node of the solid element. Therefore, the phenomenon that the fiber moves in the resin when the fiber reinforced resin is deformed is accurately reproduced.

  In another aspect of the compression molding analysis system of the present invention, the fiber model unit sets a material property of a fiber measured from an actual object in the beam element, and a fiber volume fraction included in the resin model of the fiber model. The cross-sectional area of the beam element is set so that becomes the same as the actual fiber volume fraction, and the resin model portion adds (1-fiber volume fraction) to the material properties measured from the actual solid element. You may set the multiplied value.

  In this embodiment, “material properties” is a concept including mass density, Young's modulus, yield stress, viscosity, and the like. “Fiber volume fraction” represents the proportion of the volume of fibers contained in the fiber reinforced resin.

  According to this aspect, in the fiber model portion, the fiber volume fraction included in the resin model of the fiber model by the beam element in which the material properties of the fiber measured from the real thing are set to be the same as the real fiber volume fraction. The cross-sectional area of the beam element is set so that Therefore, modeling with a fiber model using beam elements is performed with a smaller number than the actual fibers. The resin model unit sets a value obtained by multiplying the solid element by the material property measured from the actual product by (1-fiber volume fraction). Therefore, the resin model is provided with material characteristics excluding the volume of the fiber model.

  In another aspect of the compression molding analysis system of the present invention, the fiber model unit models the fiber with a two-node beam element having cross-sectional shape information, and the resin model unit converts the resin into four nodes. You may model with a tetrahedral solid element.

According to this aspect, the fiber model part models the fiber with a two-node beam element having cross-sectional shape information, and the resin model part models the resin with a solid element of a four-node tetrahedron. Appropriate compression molding analysis of fiber reinforced resin is performed.

  In another aspect of the compression molding analysis system of the present invention, the fiber may be a fiber having any length.

  According to this aspect, since the fiber has a fiber length of any length, compression molding analysis is appropriately performed for all types of fiber reinforced resin.

  In another aspect of the compression molding analysis system of the present invention, the resin may be a thermoplastic resin or a thermosetting resin.

  According to this aspect, since the resin is a thermoplastic resin or a thermosetting resin, compression molding analysis is appropriately performed on various types of fiber reinforced resins.

  In order to solve the above-described problems, one aspect of the compression molding analysis method of the present invention is a compression molding analysis method for fiber reinforced resin in which fibers are mixed with resin, and the fibers are modeled as fiber models using beam elements. The resin is modeled as a resin model with a solid element so as to have a node that is not shared with the node of the beam element, a large deformation analysis is performed on the beam element and the solid element, and the beam element and the solid element A coupled analysis between the movement of the beam element and the movement of the solid element is performed based on the movement of the beam element and the movement of the solid element during large deformation.

  According to one aspect of the compression molding analysis method of the present invention, the fibers are modeled with beam elements. The resin is modeled with solid elements to have nodes that are not shared with the beam element nodes. Large deformation analysis is performed on beam elements and solid elements. Based on the movement of the beam element and the movement of the solid element during the large deformation of the beam element and the solid element, a coupled analysis of the movement of the beam element and the movement of the solid element is performed. Therefore, the behavior of the fiber and the resin, the orientation (orientation) of the fiber, and the density of the fiber are predicted with high accuracy.

  In order to solve the above-described problem, one aspect of the compression molding analysis program of the present invention is a compression molding analysis program for causing a computer to perform a compression molding analysis of a fiber reinforced resin obtained by mixing a fiber with a resin. Modeling a fiber as a fiber model with a beam element; modeling the resin as a resin model with a solid element so as to have nodes that are not shared with the nodes of the beam element; and the beam element and the solid element Performing a large deformation analysis for the beam element and the solid element based on the beam element and the solid element movement during the large deformation of the beam element and the solid element. And a step of performing a coupled analysis.

  According to one aspect of the compression molding analysis program of the present invention, the fibers are modeled with beam elements. The resin is modeled with solid elements to have nodes that are not shared with the beam element nodes. Large deformation analysis is performed on beam elements and solid elements. Based on the movement of the beam element and the movement of the solid element during the large deformation of the beam element and the solid element, a coupled analysis of the movement of the beam element and the movement of the solid element is performed. Therefore, the behavior of the fiber and the resin, the orientation (orientation) of the fiber, and the density of the fiber are predicted with high accuracy.

  ADVANTAGE OF THE INVENTION According to this invention, the deformation | transformation behavior of the fiber and resin, the orientation (direction) of a fiber, and the density of a fiber in the compression molding of a fiber reinforced resin can be predicted with high accuracy.

It is a figure showing a schematic structure of a compression molding analysis system of one embodiment concerning the present invention. It is a figure which shows the functional block of a compression molding analysis system. It is a figure which shows the image which image | photographed the real thing of discontinuous long fiber reinforced resin. It is the figure which expanded the upper surface of discontinuous long fiber reinforced resin. It is a figure for demonstrating the lamination | stacking state in the side view of discontinuous long fiber reinforced resin. (A) is a figure which shows the tetra-solid element which comprises a resin model, (B) is a figure which shows the beam element which comprises a fiber model, and the solid element which comprises a resin model. It is a figure which shows a fiber model. It is a figure which shows a resin model. It is a figure which shows the metal mold | die used for compression molding. (A) to (E) is a figure explaining the change of the shape of the discontinuous long fiber reinforced resin by compression molding. It is a figure for demonstrating the fiber reinforced resin analysis model of compression molding analysis simulation. It is a flowchart which shows the operation | movement of a compression molding analysis simulation. It is a figure which shows the mass of the node of a beam element, and each node of a solid element. It is a figure which shows the velocity vector of the node of a beam element, and each node of a solid element. It is a figure which shows the acceleration vector of the node of a beam element, and each node of a solid element. It is a figure which shows the local coordinate system comprised by the axial direction and the orthogonal | vertical direction of a beam element. It is a figure for demonstrating a beam solid constraint coupling function, (A) is a figure which shows the velocity vector of the node of a beam element, and each node of a solid element, (B) is an axial direction of a beam element And (C) is a diagram in which the velocity vector of each node of the beam element and each node of the solid element is decomposed into an axial component and a vertical component according to a local coordinate system composed of FIG. 4D is a diagram showing only the axial component of the velocity vector of each node of the solid element, and FIG. 4D is a diagram showing only the vertical component of the velocity vector of the beam element and each node of the solid element. It is a figure for demonstrating a beam solid constraint coupling function, (A) is a figure which shows the mass of the node of a beam element, and each node of a solid element, (B) is the mass of the node of a beam element It is a figure which shows the mass of each node of a solid element calculated by distributing to each node of a solid element. It is a figure for demonstrating a beam solid constraint coupling function, (A) is a figure which shows the momentum vector in the perpendicular direction of the node of a beam element and each node of a solid element, (B) is a beam element It is a figure which shows the momentum vector in the vertical direction of each node of a solid element calculated by distributing the momentum vector in the vertical direction of each node to each node of the solid element. It is a figure for demonstrating a beam solid constraint coupling function, (A) is a figure which shows the momentum vector in the perpendicular direction of each node of a solid element, (B) is the perpendicular | vertical of each node of a solid element It is a figure which shows that the velocity vector in a direction is updated. It is a figure for demonstrating a beam solid constraint coupling function, (A) is the velocity vector in the vertical direction of a solid element in the position of the node of a beam element in the vertical direction of each node of a solid element. FIG. 7B is a diagram for explaining that the velocity vector is calculated by interpolation using the shape function from the velocity vector. FIG. 5B shows that the velocity vector calculated by the interpolation is updated as the velocity vector in the vertical direction of the node of the beam element. FIG. It is a figure for demonstrating a beam solid constraint coupling function. It is a figure which shows that the acceleration vector of the node of a beam element and each node of a solid element is updated by the same method as the update of the velocity vector in a perpendicular direction. It is a figure for demonstrating calculation of the resistance force of the beam element and solid element in the axial direction of a beam element, (A) is a figure which shows the displacement increment vector in the axial direction of the node of a beam element and each node of a solid element. , (B) is for explaining that the displacement increment vector in the axial direction of the solid element at the position of the node of the beam element is calculated from the displacement increment vector in the axial direction of each node of the solid element by interpolation using a shape function. FIG. It is a figure for demonstrating calculation of the resistance force of the beam element and solid element in the axial direction of a beam element, (A) is a figure which shows the slip amount increment vector in the position of the node of a beam element, (B) is a beam It is a figure which shows the slip amount in the position of the node of an element. It is a figure for demonstrating calculation of the resistance force of the beam element and solid element in the axial direction of a beam element. It is a figure for demonstrating the resistance force model by an elastic perfect plastic body. It is a flowchart which shows the calculation method of resistance force. It is a figure for demonstrating the resistance force model by an elastic perfect plastic body. It is a figure which shows the relationship between a slip amount and resistance. It is a figure for demonstrating the resistance force model by an elastic perfect plastic body, (A) is a figure which shows the state before a beam element moves, (B) is the state which moved the beam element to the positive direction of the X-axis. (C) is a diagram showing a state in which the beam element is further moved in the positive direction of the X axis, (D) is a diagram showing a state in which the beam element is moved in the negative direction of the X axis, (E) Is a diagram showing a state in which the beam element is further moved in the negative direction of the X axis, (F) is a diagram showing a state in which the beam element is further moved in the negative direction of the X axis, and (G) is a resistance force against the slip amount. It is a figure which shows the change of. It is a figure explaining a compression test about the disk test piece of the discontinuous long fiber reinforced resin of different temperature, (A) is so that the disk test piece of the discontinuous long fiber reinforced resin of diameter d1 may become height h4. The figure which shows the state piled up, (B) is a figure which shows the state which compressed the flat plate at the same speed as compression molding, heating a disk test piece. It is a figure which shows the result of the compression test done about the disk test piece of the discontinuous long fiber reinforced resin of different temperature. It is a figure for demonstrating the resistance force model by a temperature dependence elastic perfect plastic body. It is a flowchart which shows the calculation method of temperature-dependent resistance. It is a figure for demonstrating the resistance force model by a temperature dependence elastic perfect plastic body. It is a figure which shows calculating the temperature of the solid element in the position of the node of a beam element from the temperature of each node of a solid element by interpolation. It is a figure for demonstrating the resistance force model by a temperature dependence elastic perfect plastic body. It is a figure which shows the relationship between a slip amount and temperature-dependent resistance. It is a figure which shows the temperature distribution of the resin model using a heat conduction / structure coupling analysis function step by step according to a time series, and (C) to (C) show the heat conduction / structure coupling analysis function in this embodiment. It is a figure which shows the temperature distribution of the used resin model in steps according to a time series. It is a figure which shows the result of the compression molding test of discontinuous long fiber reinforced resin, and the result of a compression molding analysis simulation, (A)-(C) performed the compression molding test with respect to discontinuous long fiber reinforced resin. The figure which shows the image which image | photographed stepwise the shape of the discontinuous long fiber reinforced resin in the case according to a time series, (D) to (F) by the compression molding analysis simulation of this embodiment corresponding to each step of the compression molding test It is a figure which shows the shape of a resin model. It is a figure which shows the result of the compression molding test of discontinuous long fiber reinforced resin, and the result of a compression molding analysis simulation, (A) and (B) performed the compression molding test with respect to discontinuous long fiber reinforced resin. The figure which shows the image which image | photographed stepwise the shape of the discontinuous long fiber reinforced resin in the case according to a time series, (C) And (D) is by the compression molding analysis simulation of this embodiment corresponding to each step of a compression molding test It is a figure which shows the shape of a resin model. It is a figure which shows the result of the compression molding analysis simulation of a discontinuous long fiber reinforced resin, (A) to (C) is a step-by-step process of forming ribs in the fiber model by the compression molding analysis simulation of this embodiment in time series. FIGS. 4D to 4F are diagrams showing the rib formation process in the resin model according to the compression molding analysis simulation of the present embodiment in a time-series manner in correspondence with (A) to (C). It is a figure which shows the result of the compression molding test of a discontinuous long fiber reinforced resin, and the result of a compression molding analysis simulation, (A) is the discontinuity at the time of performing the compression molding test with respect to a discontinuous long fiber reinforced resin. The figure which shows the image which image | photographed the shape of the outer edge part of long fiber reinforced resin, (B) is a figure which shows the shape of the outer edge part of the fiber model by the compression molding analysis simulation of this embodiment, (C) is discontinuous long fiber reinforcement The figure which shows the image which image | photographed the stage in the middle of the rib of the discontinuous long fiber reinforced resin at the time of performing the compression molding test with respect to resin, (D) is in the middle by the compression molding analysis simulation of this embodiment It is a figure which shows the rib on the fiber model in a stage. It is a figure which shows the result of the compression molding test of a discontinuous long fiber reinforced resin, and the result of a compression molding analysis simulation, (A) is the discontinuity at the time of performing the compression molding test with respect to a discontinuous long fiber reinforced resin. The figure which shows the image which image | photographed the cross section of the rib of long fiber reinforced resin, (B) is a figure which shows the cross section of the rib on the fiber model by the compression molding analysis simulation of this embodiment.

  Hereinafter, a compression molding analysis system according to an embodiment of the present invention will be described in detail with reference to the drawings. FIG. 1 is a diagram showing a schematic configuration of a compression molding analysis system 100 according to the present embodiment. As shown in FIG. 1, the compression molding analysis system 100 of this embodiment includes a central processing unit 1, a display device 2, a storage device 3, an input device 4, and an output device 5.

  The compression molding analysis system 100 of this embodiment is a system that performs an analysis when discontinuous long fiber reinforced resin is used as a fiber reinforced resin and compression molding is used as a molding method. In this embodiment, an example in which a flow simulation is performed when a discontinuous long fiber reinforced resin is compression molded using the compression molding analysis system 100 will be described.

  The central processing unit 1 is a device such as a personal computer that can execute a program, and includes a CPU, a memory, and the like. The display device 2 is a device such as a liquid crystal display capable of displaying characters and images. The storage device 3 is a device capable of storing a program and data such as an HDD (Hard Disk Drive), and may use an external database server or the like. The program of the present invention is stored in the storage device 3. The input device 4 is a device such as a keyboard that allows a user to input data or instructions. The output device 5 is a device such as a printer that can output characters and images. In the compression molding analysis system 100 of the present embodiment, the output device 5 may be omitted.

  FIG. 2 is a diagram showing functional blocks that function when the central processing unit 1 executes the program of the present invention. As shown in FIG. 2, the central processing unit 1 functions as a control unit 10. Moreover, the control part 10 functions as the fiber model part 11, the resin model part 12, the coupled analysis part 13, and the large deformation | transformation analysis part 14 according to the program of this invention.

  The fiber model unit 11 models the discontinuous long fiber with a beam element of a two-node finite element method having cross-sectional shape information, and sets the material properties of the fiber measured from the actual object. The fiber model unit 11 is modeled by beam elements with a smaller number than the actual discontinuous long fibers. Therefore, the fiber model unit 11 sets the cross-sectional area of the beam element so that the fiber model of the beam element has the same fiber volume fraction contained in the resin model as the actual fiber volume fraction. Details of the fiber model by the beam element will be described later. By modeling with fewer beam elements than actual discontinuous long fibers, the calculation time in large deformation analysis and coupled analysis can be shortened.

  The resin model unit 12 models the resin with a tetra-solid element of a four-node tetrahedron finite element method. The resin model unit 12 is set so that the nodes of the solid element are not shared with the nodes of the beam element. Further, the resin model unit 12 has a value obtained by multiplying the material characteristic measured from the actual resin by (1-fiber volume fraction) in order to give the resin model the material characteristic excluding the volume of the fiber model. Set. Details of the resin model using solid elements will be described later.

  The coupled analysis unit 13 performs coupled analysis between the movement of the beam element and the movement of the solid element. Coupled analysis means that analysis is performed while two different phenomena influence each other through interaction. The coupled analysis unit 13 performs coupled analysis on the vertical movement of the beam element by the beam solid constraint coupling function, and detects the beam element in the axial direction of the beam element. The sliding amount of the solid element is acquired, the resistance is calculated from the acquired sliding amount, and the coupled analysis is performed. This reproduces the phenomenon that the discontinuous long fibers move in the resin when the discontinuous long fiber reinforced resin is deformed. Details of the beam solid constraint coupling function and the calculation of the slip amount and the resistance force will be described later.

  The large deformation analysis unit 14 performs large deformation analysis on the beam elements constituting the fiber model and the solid elements constituting the resin model. The large deformation analysis unit 14 performs large deformation analysis on the solid element while recreating the solid element so that the shape of the solid element does not become an irregular shape by using the adaptive remeshing function. Further, the large deformation analysis unit 14 performs large deformation analysis on the beam element without using the adaptive remeshing function because the beam element does not have an irregular shape.

(Discontinuous fiber reinforced resin)
Next, the discontinuous fiber reinforced resin modeled by the fiber model part 11 and the resin model part 12 will be described.
Fiber reinforced resin (FRP (Fiber Reinforced Plastics)) is a composite material in which fibers are mixed with a resin, and glass fibers (Glass) or carbon fibers (Carbon) are often used as the fibers. Although glass fiber is lighter than steel, it has the characteristics that rigidity and intensity | strength are lower than steel. Carbon fiber is a material that can replace steel because it is lighter than steel and has rigidity and strength equal to or greater than that of steel. However, carbon fiber is more expensive than glass fiber or steel.

  The resin includes a thermosetting resin and a thermoplastic resin. Examples of thermosetting resins include epoxy resins and phenol resins. Thermosetting resins are generally soft at room temperature and harden by a chemical reaction at high temperatures, which solidifies, but the chemical reaction takes time. The thermosetting resin remains solid even after cooling to room temperature after curing, and has irreversibility. Because of this property, it is generally difficult to recycle.

  Examples of the thermoplastic resin include polypropylene and polyamide. Thermoplastic resins are generally hard at room temperature, soft when heated to a temperature higher than the melting temperature of the resin, and hard when returned to room temperature, and therefore have reversibility. Because of this property, it is generally easy to recycle. The thermoplastic resin has an advantage of high productivity because it does not require a chemical reaction for curing the resin unlike the thermosetting resin.

  The fibers used for the fiber reinforced resin include continuous fibers and discontinuous fibers. The continuous fiber is a fiber having a fiber length of 50 mm or more and is very long. A discontinuous fiber is called a discontinuous long fiber when the fiber length is 2 mm to 50 mm, and is called a discontinuous short fiber when the fiber length is 2 mm or less. The longer the fiber, the better the strength characteristics. When the fiber length is 10 mm or more, high strength is exhibited, and when the fiber length is 40 mm or more, high toughness is exhibited.

  In the present embodiment, as an example, a flow simulation is performed when a product name is Tepex (registered trademark) flowcore and a discontinuous long fiber reinforced resin manufactured by Bond-Laminates GmbH is compression-molded. FIG. 3 is a view showing an image obtained by photographing the actual discontinuous long fiber reinforced resin. FIG. 4 is an enlarged view of the upper surface of the discontinuous long fiber reinforced resin. In the discontinuous long fiber reinforced resin 20 shown in FIG. 3, the fiber is glass, and the fiber length is 30 to 50 mm. As shown in FIG. 4, the fiber orientation of the fibers 21 is random in two dimensions. As the resin, thermoplastic polyamide nylon 6 is used.

  FIG. 5 is a view for explaining a laminated state of the discontinuous long fiber reinforced resin in a side view. The plate thickness d of the discontinuous long fiber reinforced resin 20 shown in FIG. 5 is 2 mm, and about 40 fiber layers 22 are laminated. The fiber volume content is 0.47 (47%). The fiber volume fraction is the ratio of the volume of fibers contained in the fiber reinforced resin.

  The discontinuous long fiber reinforced resin 20 as described above can be formed into a complicated shape by compression molding. Compression molding is a technique in which a material is pressed and molded using a mold composed of a male mold and a female mold. A feature of compression molding is that the fibers contained in the discontinuous long fiber reinforced resin 20 can be molded into a desired shape while maintaining the long fibers without breaking, and thus can have high strength characteristics even after molding. Further, in the case of the thermoplastic discontinuous long fiber reinforced resin 20, the discontinuous long fiber reinforced resin heated to a high temperature is molded while being cooled in a low temperature mold and cured in about several minutes. The advantage is that the time to complete the parts is short.

(Fiber model and resin model)
Next, the fiber model and resin model by the fiber model part 11 and the resin model part 12 are demonstrated. FIG. 6A is a diagram showing tetrasolid elements constituting a resin model. FIG. 6B is a diagram showing beam elements constituting the fiber model and tetrasolid elements constituting the resin model. In FIG. 6B, tetrasolid elements are schematically represented in two dimensions. The same applies to the tetrasolid elements shown in FIG. 13 and subsequent drawings. FIG. 7 is a diagram showing a fiber model. FIG. 8 is a diagram showing a resin model.

  The length of the beam element 30 shown in FIG. 6 is 1 mm as an example. Therefore, when the fiber length of the discontinuous long fiber reinforced resin 20 is 40 mm, 40 beam elements 30 are connected. Each beam element 30 is connected by a node 31. The node 31 of the beam element 30 is set so as not to be shared with the node 41 of the solid element 40.

  The material model of the beam element 30 is an elastic body or an elastic-plastic body. The material properties (mass density, Young's modulus, Poisson's ratio, yield stress, work hardening coefficient) are set to the same values as those of the actual fiber 21. In addition, since the bending rigidity of the real fiber 21 is very small, the bending rigidity of the beam element 30 is also set to be very small or set to zero. The number of fibers 21 contained in the actual discontinuous long fiber reinforced resin 20 is enormous, and it is not practical to create a fiber model with the same number of beam elements because it takes too much calculation time. Therefore, as shown in FIG. 7, the fiber model 50 is set to be smaller than the number of real fibers 21. However, the cross-sectional area of the cross-sectional shape of the beam element 30 is set so that the fiber volume fraction is the same as that of the actual discontinuous long fiber reinforced resin 20.

  A comparison between the actual discontinuous long fiber reinforced resin 20 and the fiber model 50 is as follows. For example, in the case of the discontinuous long fiber reinforced resin 20 having a length of 100 mm, a width of 100 mm, a plate thickness of 2 mm, a fiber volume fraction of 47%, and a fiber length of 40 mm, the actual discontinuous long fiber reinforced resin 20 The fiber model 50 is compared as follows.

<Real discontinuous long fiber reinforced resin 20>
Fiber diameter: 5-10 μm
Number of fibers: 50,000,000-120,000,000
Number of layers: 32 to 48 layers

<Fiber model 50>
Fiber diameter: 100-500 μm
Number of fibers: 50,000 to 1,200,000
Number of layers: 2 to 12 layers

  A resin model 60 shown in FIG. 8 models the resin with tetrasolid elements (four-node tetrahedral elements) 40. As shown in FIG. 6, each surface of the solid element 40 is connected by a node 41. The material model of the resin model 60 is a temperature-dependent viscoelastic body. The material characteristics (mass density, Young's modulus at each temperature, Poisson's ratio at each temperature, viscosity at each temperature) are set to values obtained by multiplying the characteristics of the actual resin by (1-fiber volume fraction). By setting in this way, the resin model can have material characteristics excluding the volume of the fiber model.

(Compression molding)
Next, compression molding for performing compression molding analysis simulation of the present embodiment will be described. FIG. 9 is a diagram showing a mold used for compression molding. FIG. 10 is a diagram for explaining a change in the shape of the discontinuous long fiber reinforced resin due to compression molding.

  As shown in FIG. 9, the mold used for compression molding is divided into a male mold 70 and a female mold 80. FIG. 9 shows a state in which the male mold 70, the discontinuous long fiber reinforced resin 20, and the female mold 80 are cut on the ZX plane. As shown in FIG. 9, the male mold 70 has a flat plate portion 71 having a width W2 of 130 mm in the X direction and the Y direction, and a width W1 in the X direction and the Y direction is formed at the central portion of the flat plate portion 71. A convex portion 72 of 68 mm is formed. Inside the convex part 72, a rib forming part 73 having a height h1 of 15 mm is formed. The male mold 70 is controlled to be kept at a constant temperature after being heated to 160 ° C.

  Similar to the male mold 70, the female mold 80 has a flat plate portion 81 with a width W2 in the X direction and the Y direction of 130 mm, and the central portion of the flat plate portion 81 has a width W1 in the X direction and the Y direction. A 68 mm concave portion 82 is formed. The depth h2 of the concave portion 82 is 25.4 mm. Similar to the male mold 70, the female mold 80 is controlled to be maintained at a constant temperature after being heated to 160 ° C.

  As the discontinuous long fiber reinforced resin 20, a discontinuous long fiber reinforced resin having a thickness of 7 mm is prepared, heated to 290 ° C. as an initial temperature, placed on the female mold 80, and the male mold 70 is replaced with the female mold 80. By moving in the direction, compression molding is performed. The height h3 from the surface of the discontinuous long fiber reinforced resin 20 to the flat plate portion 71 of the male mold 70 is 25.4 mm. After the compression molding, the discontinuous long fiber reinforced resin 20 has a thickness of 2 to 3 mm.

  As shown in FIGS. 10A to 10E, the discontinuous long fiber reinforced resin 20 is deformed by the movement of the male mold 70 and the female mold 80, and the discontinuous long fiber reinforced resin is formed in the rib forming portion 73. As the fluid 20 flows, the ribs 24 are formed.

(Compression molding analysis simulation)
FIG. 11 is a diagram for explaining a fiber-reinforced resin analysis model of a compression molding analysis simulation in the present embodiment. In this embodiment, the flow simulation at the time of compression molding of the discontinuous long fiber reinforced resin 20 as described above is performed by the large deformation analysis of the fiber model 50 and the large deformation analysis of the resin model 60 shown in FIG. At this time, since the fiber model 50 and the resin model 60 affect the fluidity of each other, the movement of the beam element 30 constituting the fiber model 50 and the movement of the solid element 40 constituting the resin model 60 are coupled. By performing the analysis, a flow simulation of the analysis model 90 is performed.

  In this embodiment, LS, which is commercial coded by Livermore Software Technology Corporation as software for performing compression molding analysis simulation as described above, is widely used in academia and industry as an advanced nonlinear finite element solver. -DYNA was used.

  In LS-DYNA, the phenomenon time is divided into small intervals at a certain time interval, and the deformation state at each time is subjected to sequential deformation analysis by a direct integration method called explicit method or implicit method. Each time is called a step.

  LS-DYNA has a function of performing large deformation analysis of beam elements.

  LS-DYNA has an adaptive remeshing function for tetra-solid elements developed for performing metal forging. In a large deformation analysis such as metal forging, solid elements become irregular shapes, causing problems such as deterioration in analysis accuracy, increase in calculation time, and analysis calculation being stopped halfway. To avoid this problem, the solid element is recreated before the solid element becomes distorted, and the stress and strain that the solid element had before it was recreated are automatically recreated. An adaptive remeshing function that performs analysis calculation while taking over is used. Furthermore, the adaptive remeshing function can be coupled with heat conduction analysis, and can perform hot forging analysis of metals. In the hot forging analysis of a metal, a metal model that is heated to a high temperature and has a low Young's modulus and yield stress and improved formability is formed while being cooled to a low-temperature mold. In the present invention, a 3D adaptive EFG (Element Free Galerkin) function having excellent calculation accuracy is used among the adaptive remeshing functions. The 3D adaptive EFG function is an adaptive remeshing function for tetra-solid elements based on the mesh-free Galerkin method. An adaptive remeshing function other than the 3D adaptive EFG function can also be used.

  LS-DYNA has a beam, solid, constraint, and coupling function developed to perform the structural analysis function of reinforced concrete. In the reinforced concrete analysis function, the reinforcing bar is modeled as a beam element and the concrete is modeled as a solid element, and the movement of the beam element and the movement of the solid element are calculated by the coupled analysis using the beam solid constraint coupling function.

  In this embodiment, the large deformation behavior of the beam element 30 of the fiber model is calculated by the large deformation analysis function of the beam element. In the present embodiment, as described above, the adaptive remeshing function generally used for metal forging analysis is applied to the compression molding analysis of the resin model, and the large deformation behavior of the solid element 40 of the resin model is Calculate using the adaptive remeshing function. Furthermore, in the present embodiment, as described above, the beam solid constraint coupling function generally used for structural analysis of reinforced concrete is analyzed by compression molding analysis of a fiber reinforced resin composed of a resin model and a fiber model. And the coupled analysis of the movement of the fiber element 30 of the fiber model and the movement of the solid element 40 of the resin model is calculated by the beam solid constraint coupling function.

  In the structural analysis of reinforced concrete, the reinforcing bars are firmly constrained in the concrete and hardly slip, but in this embodiment, the fiber model beam element is a soft (low viscosity) resin model solid element. By setting it so that it can move while receiving resistance, it succeeded in reproducing the deformation behavior of fiber reinforced resin during compression molding.

  FIG. 12 is a flowchart showing the operation of the compression molding analysis simulation in the present embodiment. First, the fiber model unit 11 models the discontinuous long fiber reinforced resin 20 with the beam element 30 (S1). Further, the resin model unit 12 models the resin with the solid element 40 so as to have the node 41 that is not shared with the node 31 of the beam element 30 (S1).

Next, the control unit 10 performs a flow simulation (S2). At this time, the large deformation analysis unit 14 performs large deformation analysis on the beam element 30 without using the adaptive remeshing function (S3). The large deformation analysis unit 14 performs large deformation analysis on the solid element 40 using the adaptive remeshing function (S5). Then, the coupled analysis unit 13 performs coupled analysis between the movement of the beam element 30 and the movement of the solid element 40 by using the beam / solid / constraint / coupling function and the slip resistance (S4).
When the flow simulation ends, the control unit 10 performs output processing using the output device 5 (S5).

(Beam, solid, constraint, coupling function)
Next, the beam solid constraint coupling function will be described. 13 to 22 are views for explaining the beam solid constraint coupling function.

  As shown in FIGS. 13 to 15, the node 31 of the beam element 30 and the node 41 of the solid element 40 have a mass, a velocity vector, and an acceleration vector.

In FIG. 13, m ₁ represents the mass of the node 31 of the beam element 30, and M _1, M _2, M _{3, and} M ₄ represent the mass of the node 41 of the solid element 40. In FIG. 14, v ₁ represents the velocity vector of the node 31 of the beam element 30, and V _1, V _2, V _{3, and} V ₄ represent the velocity vector of the node 41 of the solid element 40. In FIG. 15, a ₁ represents the acceleration vector of the node 31 of the beam element 30, and A _1, A _2, A _{3, and} A ₄ represent the acceleration vector of the node 41 of the solid element 40.

  Note that the subscript numbers of m, M, v, V, a, and A indicate the numbers given to the nodes 31 and 41 for explanation. The same applies to the subscript numbers below.

Here, as shown in FIG. 16, a local coordinate system configured by the axial direction LCx of the beam element 30 and the vertical direction LCy perpendicular to the axial direction LCx is set. The velocity vector v ₁ of the node 31 of the beam element 30 shown in FIG. 17A is perpendicular to the axial component v _{T, 1 of} the beam element 30 according to the local coordinate system, as shown in FIG. It is decomposed into a directional component v _{N, 1} . Incidentally, the subscript T represents Tangential which means the tangential direction as the direction along the axis, and N represents Normal which means the vertical direction.

As shown in FIG. 17B, the velocity vectors V _1, V _2, V _{3, and} V _{4 of the} nodes 41 of the solid element 40 are also expressed in the axial direction components V _{T, 1} , V _{T, 2} , V _{T, 3} , V _{T, 4} and the vertical components V _{N, 1} , V _{N, 2} , V _{N, 3} , V _{N, 4} are decomposed.

FIG. 17C shows the axial component v of the velocity vector v ₁ of the node 31 of the beam element 30 and the velocity vectors V _1, V _2, V _{3, and} V _{4 of the} nodes 41 of the solid element 40. _{T, 1} , V _{T, 1} , V _{T, 2} , V _{T, 3} , V _{T, 4} are shown.

FIG. 17D shows the vertical component v of the velocity vector v ₁ of the node 31 of the beam element 30 and the velocity vectors V _1, V _2, V _{3, and} V _{4 of the} nodes 41 of the solid element 40. _{N, 1} , _{VN, 1} , _{VN, 2} , _{VN, 3} , _{VN, 4} are shown.

  In the present embodiment, in the vertical direction of the beam element 30, the velocity coupling calculation and the acceleration coupling calculation are performed using the constraint coupling method, and the nodal point 31 of the beam element 30 and the solid element 40. The movement of the node 41 is calculated. Further, with respect to the axial direction of the beam element 30, the resistance force is calculated according to the relative displacement (slip amount) between the beam element 30 and the solid element 40, and the node 31 of the beam element 30 and the node of the solid element 40 are calculated. 41 movements are calculated. Hereinafter, calculation of the movement of the node 31 in the beam element 30 and the movement of the node 41 in the solid element 40 will be described separately for the vertical direction of the beam element 30 and the axial direction of the beam element 30.

<Coupling calculation of beam element vertical velocity and acceleration>
First, the vertical velocity coupling calculation and acceleration coupling calculation of the beam element 30 will be described.
In the velocity coupling calculation, the mass m ₁ of the node 31 of the beam element 30 shown in FIG. 18 (A) is determined by the shape functions N ₁ , N ₂ , N ₃ , and N ₄ of the solid element 40 as shown in FIG. As shown in FIG. The masses M ′ _1, M ′ _2, M ′ _{3 and} M ′ ₄ of the nodes 41 of the solid element 40 after distribution are expressed as follows.

M ′ ₁ = M ₁ + N ₁ · m ₁
M ′ ₂ = M ₂ + N ₂ · m ₁
M ′ ₃ = M ₃ + N ₃ · m ₁
M ′ ₄ = M ₄ + N ₄ · m ₁

Similarly, the momentum vector m ₁ · v _{N, 1} in the vertical direction of the node 31 of the beam element 30 shown in FIG. 19A is obtained by the shape function N ₁ , N ₂ , N ₃ , N ₄ of the solid element 40. 19B, it is distributed to the nodes 41 of the solid element 40. Momentum vectors (M ₁ · V _{N, 1} ) ', (M ₂ · V _{N, 2} )', (M ₃ · V _{N, 3} ) ', () in the vertical direction of the nodes 41 of the solid element 40 after distribution M ₄ · V _{N, 4} ) ′ is expressed as follows.

(M ₁ · V _{N, 1} ) '= M ₁ · V _{N, 1} + N ₁ (m ₁ · V _{N, 1} )
(M ₂ · V _{N, 2} ) '= M ₂ · V _{N, 2} + N ₂ (m ₁ · v _{N, 1} )
(M ₃ · V _{N, 3} ) '= M ₃ · V _{N, 3} + N ₃ (m ₁ · V _{N, 1} )
(M ₄ · V _{N, 4} ) '= M ₄ · V _{N, 4} + N ₄ (m ₁ · v _{N, i} )

Next, momentum vectors (M ₁ · V _{N, 1} ) ′, (M ₂ · V _{N, 2} ) ′, (M ₃ · V) in the vertical direction of the nodes 41 of the solid element 40 shown in FIG. _{N, 3} ) ′, (M ₄ · V _{N, 4} ) ′, the mass M of the node 41 of the solid element 40 to which the mass m ₁ of the node 31 of the beam element 30 is distributed as shown in FIG. _'1, M' _2, M divided by _'3, M' _4, the velocity vector in the vertical direction of the nodal 41 of solid element _{40 V N, 1, V N} , 2, V N, 3, V N, Update ₄ . The updated velocity vectors V ′ _{N, 1,} V ′ _{N, 2,} V ′ _{N, 3,} V ′ _{N, 4} are expressed as follows.

V ′ _{N, 1} = [M ₁ · V _{N, 1} + N ₁ (m ₁ · v _{N, 1} )] / [M ₁ + N ₁ · m ₁ ]
V ′ _{N, 2} = [M ₂ · V _{N, 2} + N ₂ (m ₁ · v _{N, 1} )] / [M ₂ + N ₂ · m ₁ ]
V ′ _{N, 3} = [M ₃ · V _{N, 3} + N ₃ (m ₁ · v _{N, 1} )] / [M ₃ + N ₃ · m ₁ ]
V ′ _{N, 4} = [M ₄ · V _{N, 4} + N ₄ (m ₁ · v _{N, 1} )] / [M ₄ + N ₄ · m ₁ ]

Next, as shown in FIG. 21 (A), the velocity vector in the vertical direction of the solid element 40 at the position of the node 31 of the beam element 30 is used as the velocity in the vertical direction of the node 41 of the updated solid element 40. Calculate from the vectors V ′ _{N, 1,} V ′ _{N, 2,} V ′ _{N, 3,} V ′ _{N, 4} by interpolation using the shape function N ₁ , N ₂ , N ₃ , N ₄ of the solid element 40 . As shown in FIG. 21B, the velocity vector v _{N, 1} in the vertical direction of the node 31 of the beam element 30 is updated according to the value. The updated velocity vector v ′ _{N, 1} is expressed as follows.

v ′ _{N, 1} = ΣN _k · V ′ _{N, k} (k indicates the number of the node 41 of the solid element 40)

As in the case of the velocity vector, the acceleration vectors A _{N, 1,} A _{N, 2,} A _{N, 3,} A _{N, 4} in the vertical direction of the node 41 in the solid element 40 are updated. Next, the acceleration vector in the vertical direction of the solid element 40 at the position of the node 31 of the beam element 30 is used as the updated acceleration vector A ′ _{N, 1,} A ′ _N in the vertical direction of the node 41 of the solid element 40. _{, 2,} A ′ _{N, 3 and} A ′ _{N, 4} are calculated by interpolation using the shape functions N ₁ , N ₂ , N ₃ and N ₄ of the solid element 40. Next, as shown in FIG. 22, the acceleration vector a _{N, 1} in the vertical direction of the node 31 of the beam element 30 is updated with the value. The updated acceleration vector a ′ _{N, 1} is expressed as follows.

a ′ _{N, 1} = ΣN _k · A ′ _{N, k} (k indicates the number of the node 41 of the solid element 40)

<Calculation of resistance in the beam element axis direction>
Next, calculation of the resistance force between the beam element 30 and the solid element 40 in the axial direction of the beam element 30 will be described. 23 to 25 are diagrams for explaining the calculation of the resistance force between the beam element 30 and the solid element 40 in the axial direction of the beam element 30. FIG. By multiplying the velocity vector in the axial direction of the node 31 of the beam element 30 and the node 41 of the solid element 40 shown in FIG. 17C by a time increment ΔT described later, the node of the beam element 30 shown in FIG. 31 of the displacement increment vector Δd _{T, 1 in the} axial direction and the displacement increment vectors ΔD _{T, 1} , ΔD _{T, 2} , ΔD _{T, 3} , ΔD _{T, 4} of the node 41 of the solid element 40 in the axial direction. Calculated.

Next, as shown in FIG. 23B, the displacement increment vector in the axial direction of the solid element 40 at the position of the node 31 of the beam element 30 is replaced with the displacement increment vector in the axial direction of the node 41 of the solid element 40. Calculation is performed from ΔD _{T, 1} , ΔD _{T, 2} , ΔD _{T, 3} , ΔD _{T, 4} by interpolation using the shape functions N ₁ , N ₂ , N ₃ , N ₄ of the solid element 40. The axial displacement increment vector ΔD ′ _T of the solid element 40 at the calculated position of the node 31 of the beam element 30 is expressed as follows.

ΔD ′ _T = ΣN _k · ΔD _{T, k} (k indicates the number of the node 41 of the solid element 40)

FIG. 23B shows the displacement increment vector Δd _{T, 1} in the axial direction of the node 31 of the beam element 30 and the displacement increment vector ΔD ′ in the axial direction of the solid element 40 at the position of the node 31 of the beam element 30. _T is shown.

Relative increment vector Δd _{T, 1} in the axial direction of the node 31 of the beam element 30 minus the displacement increment vector ΔD ′ _T in the axial direction of the solid element 40 at the position of the node 31 of the beam element 30 is relative. Represents a typical displacement increment, resulting in a slip increment vector ΔS ₁ . The slip amount increment vector ΔS ₁ is expressed as follows.

ΔS ₁ = Δd _{T, 1} −ΔD ′ _T

FIG. 24A shows a slip increment vector ΔS ₁ at the position of the node 31 of the beam element 30.

The slip amount S ₁ is calculated by integrating the slip amount increment vector ΔS ₁ for a predetermined period.

FIG. 24B shows the slip amount S ₁ at the position of the node 31 of the beam element 30.

The resistance force in the axial direction of the beam element 30 is calculated from the slip amount increment vector ΔS ₁ by a “resistance force model by an elastic perfect plastic body” and a “resistance force model by a temperature-dependent elastic perfect plastic body” described later. . As shown in FIG. 25, the resistance force is given to the node 31 of the beam element 30 as a resistance vector Fb ₁ , and the shape function N ₁ , N ₂ , N ₃ , N of the solid element 40 is applied to the node 41 of the solid element 40. _The resistance force distributed by ₄ is given in the opposite direction as resistance vector Fs1 _, Fs2 _, Fs3 _, Fs4. Resistance force vector _{Fs 1, Fs 2, Fs 3} , Fs 4 is expressed as follows.

Fs ₁ = N ₁ · Fb ₁
Fs ₂ = N ₂ · Fb ₁
Fs ₃ = N ₃ · Fb ₁
Fs ₄ = N ₄ · Fb ₁

  That is, in the present embodiment, the slip amount increment vector between the beam element 30 and the solid element 40 at the position of the node 31 of the beam element 30 is calculated in the axial direction of the beam element 30, and the node 31 of the beam element 30 is If it is recognized that slip has occurred in the axial direction, a resistance force is applied in the opposite direction. And the process which gives resistance force as reaction force to the node 41 of the solid element 40 is performed so that the resistance force may be offset.

  The axial movement of the beam element 30 at the node 31 of the beam element 30 and the node 41 of the solid element 40 is obtained by calculating the equation of motion in which the above-described velocity vector, acceleration vector, and resistance force vector are applied. .

  By performing such treatment, in compression molding where the resin of the discontinuous fiber reinforced resin is soft (low viscosity), the phenomenon that the discontinuous fiber moves while receiving resistance in the resin Can be reproduced well.

(Resistance model with elastic perfect plastic)
Next, a resistance force model using an elastic perfect plastic body will be described. FIG. 26 to FIG. 28 are diagrams for explaining a resistance force model by an elastic perfect plastic body. Note that, as described above, in the compression analysis simulation of this embodiment, the phenomenon time is finely divided by a certain time increment ΔT, and the deformation state at each time is sequentially subjected to deformation analysis by a direct integration method called explicit method or implicit method. Each time is called a step. The calculation of the resistance force described here is also performed for each step.

To calculate the resistance force, first, as shown in FIG. 26, the coupled analysis unit 13 calculates the component ΔS of the slip increment vector generated between the previous step and the current step, and the resistance force of the previous step. Fb _T-ΔT is read from the storage device 3 (S10). Note that the resistance force Fb _T-ΔT of the previous step is stored in the storage device 3 in advance.

Next, the coupled analysis unit 13 calculates the _trial resistance force Fb _trial (S11). The trial resistance force Fb _trial is calculated by the following equation.

Fb _trial = Fb _T-ΔT + K · ΔS

Here, K is the slope of the elastic region in the elastic perfect plastic slip resistance model. FIG. 27 shows an elastic perfect plastic slip resistance model. As shown in FIG. 27, in the elastic perfect plastic slip resistance model, when the slip amount is within a certain range, the resistance force Fb changes in the elastic region. However, when the slip amount S becomes S1 and the resistance force Fb reaches the maximum resistance value Fb _max , it becomes a complete plastic region, and thereafter the resistance force Fb does not change even if the slip amount increases. In such an elastic perfect plastic slip resistance model, the slope K of the elastic region is expressed by the following equation.

K = Fb _max / S1

  The inventors of the present invention have found that the behavior of the node 31 of the beam element 30 in the axial direction of the beam element 30 is considered that the resistance to the slip amount has the characteristics of an elastic perfect plastic body. The resistance force was calculated by a slip resistance model and a resistance force model based on a temperature-dependent elastic perfect plastic body, which will be described later.

  In the discontinuous long fiber reinforced resin 20 of the present embodiment, since the resin is a thermoplastic resin, the resistance force is calculated by a resistance model using a temperature-dependent elastic perfect plastic body. Before explaining the model, here, calculation of the resistance force using a resistance force model based on an elastic perfect plastic body independent of temperature will be described.

Next, the coupled analysis unit 13 calculates the yield strength in order to determine whether the trial resistance force Fb _trial is the elastic force of the elastic perfect plastic slip resistance force model or the resistance force of the complete plastic region. (S12). Yield calculation is performed by the following formula.

| Fb _trial | <Fb _max

If the absolute value of the _trial resistance force Fb _trial is smaller than the maximum resistance value Fb _{max as} a result of the yield calculation (S12: YES), the coupled analysis unit 13 determines that the trial resistance force Fb _trial is in the elastic range. Then, the trial resistance force Fb _trial is obtained as the resistance force Fb _T of this step (S13).

However, if the absolute value of the _trial resistance force Fb _trial is greater than or equal to the maximum resistance value Fb _{max as} a result of the yield calculation (S12: NO), the coupled analysis unit 13 determines that the trial resistance force Fb _trial is in the complete plastic range. It is determined that there is a resistance force Fb _T by the following equation (S14).

Fb _T = Fb _max · (Fb _trial / | Fb _trial |)

The coupled analysis unit 13 applies the resistance force Fb _T of this step obtained as described above as a load to the node 31 of the beam element 30 and the node 41 of the solid element 40 (S15).

  FIG. 28 explains how the resistance force changes with respect to the slip amount of the node 31 of the beam element 30 in the axial direction of the beam element 30 when the elastic perfect plastic slip resistance model is used. FIG. In this example, a case where the beam element 30 is forcibly moved in the fixed solid element 40 will be described in order to simplify the description.

  FIG. 28A shows a state before the beam element 30 moves. In this case, the slip amount is zero, which corresponds to the point (a) shown in FIG. As is clear from FIG. 28G, the resistance force in this case is zero.

Next, as shown in FIG. 28 (B), the beam element 30 is moved in the direction of the arrow indicated by the dotted line, that is, in the positive direction of the X axis (rightward in FIG. 28), and the slip amount is as shown in FIG. Consider a case where the point changes from point (a) to point (b). In this case, as shown in FIG. 28 (G), the resistance force becomes Fb ₁ , and the resistance force Fb ₁ is given in the direction opposite to the moving direction of the beam element 30 (left direction in FIG. 28 (B)). .

Further, as shown in FIG. 28 (C), when the beam element 30 is moved in the positive direction of the X axis, the slip amount changes from the point (b) shown in FIG. 28 (G) to the point (c). The resistance force reaches the maximum resistance force Fb _max . Thereafter, when the beam element 30 is moved in the positive direction of the X axis, the slip amount changes from the point (c) to the point (d) shown in FIG. 28 (G), but the resistance force is the maximum resistance force Fb _max. It is given to the node 31 of the beam element 30 as it is. FIG. 28C shows a state where the maximum resistance force Fb _max is applied to the node 31 of the beam element 30.

  Next, as shown in FIG. 28 (D), when the beam element 30 is moved in the negative direction of the X axis (leftward in FIG. 28), the slip amount starts from the point (d) shown in FIG. 28 (G). It changes toward the point (e), and the resistance decreases. When the slip amount reaches the point (e) shown in FIG. 28G, the resistance force becomes zero. FIG. 28D shows a state in which the resistance force at the node 31 of the beam element 30 has become zero.

Further, as shown in FIG. 28 (E), when the beam element 30 is moved in the negative direction of the X axis, a resistance force is generated in a direction opposite to the movement direction of the beam element 30 (right direction in FIG. 28). When the slip amount reaches the point (f) shown in FIG. 28G, the resistance force becomes the maximum resistance force Fb _max . FIG. 28E shows a state in which the maximum resistance force Fb _max is applied to the node 31 of the beam element 30.

Thereafter, as shown in FIG. 28 (F), when the beam element 30 is moved in the negative direction of the X axis, the amount of slip changes from the point (f) shown in FIG. 28 (G) to the point (g). However, the resistance force is applied to the node 31 of the beam element 30 with the maximum resistance force Fb _max as shown in FIG.

  Thereafter, if the node 31 of the beam element 30 is moved in the positive direction of the X-axis, the slip amount changes as shown by the dotted arrow in FIG. 28 (G), and the resistance according to the elastic perfect plastic slip resistance model is obtained. The power will change.

  When the resin in the discontinuous long fiber reinforced resin is a thermosetting resin, since there is no temperature change in compression molding, the node 31 of the beam element 30 is obtained by using the above-described elastic complete plastic slip resistance model. It is possible to calculate a resistance force when the beam element 30 is moved in the axial direction of the beam element 30.

  The node 31 of the beam element 30 is a coupling point for performing a coupling calculation with the node 41 of the solid element 40. However, any number of coupling points other than the node 31 may be set on the beam element 30. Coupled analysis can also be performed.

(Thermoplastic resin)
In the above-described elastic perfect plastic slip resistance model, the temperature parameter is not taken into consideration. However, when the resin is a thermoplastic resin, such as the discontinuous long fiber reinforced resin 20 of the present embodiment, a temperature change occurs in compression molding, and thus the hardness (viscosity) changes.

  FIG. 29 is a diagram for explaining a compression test performed on disc test pieces of discontinuous long fiber reinforced resin at different temperatures. FIG. 30 is a diagram showing the results of a compression test. In the compression test, as shown in FIG. 29A, the disc test pieces of the discontinuous long fiber reinforced resin 20 having a diameter d1 of 75 mm are stacked so that the height h4 is 6 mm, and these disc tests are performed. The pieces were heated to temperatures different from 160 ° C., 250 ° C., and 280 ° C., respectively. Then, as shown in FIG. 29B, the flat plate was compressed in the thickness direction at the same speed as the compression molding. In the following description, the direction in which the disk specimen is displaced is the thickness direction. The melting temperature of the thermoplastic polyamide nylon 6 which is the resin of the discontinuous long fiber reinforced resin 20 is 225 ° C.

  As shown in FIG. 30, when the temperature is 160 ° C., which is lower than the melting temperature, a load of about 2000 N is required to displace the disk specimen by 0.2 mm, and after that, the disk specimen is 2 mm. It can be seen that a load of about 5000 N is required for displacement.

  However, when the temperature is 250 ° C., the load necessary for displacing the disk specimen by about 0.2 mm is about 500 N, and after that, the load necessary for displacing the disk specimen by about 2 mm is about It turns out that it is about 1000N.

  Furthermore, when the temperature is 280 ° C., the load necessary to displace the disk specimen by about 0.2 mm is about 100 N, and the load necessary to displace the disk specimen by 2 mm is about It turns out that it is about 300N.

  As described above, when the resin is thermoplastic, the hardness (viscosity) varies depending on the temperature. Therefore, in the present embodiment, the following temperature-dependent elastic perfect plastic slip resistance model is used.

(Temperature dependent elastic perfect plastic slip resistance model)
FIG. 31 to FIG. 33 are diagrams for explaining the temperature-dependent elastic perfect plastic slip resistance model.

In order to calculate the resistance force, first, as shown in FIG. 31, the coupled analysis unit 13 calculates the component ΔS of the slip increment vector generated between the previous step and the present step, and the nodal point of the beam element 30. The temperature Ts of the solid element 40 at the position 31 is calculated (S20). Further, the coupled analysis unit 13 reads the resistance force Fb _T-ΔT of the previous step from the storage device 3 (S20). Note that the resistance force Fb _T-ΔT of the previous step is stored in the storage device 3 in advance.

As shown in FIG. 32, the temperature Ts of the solid element at the position of the node 31 of the beam element 30 is the temperature Ts, _K (K is the number of the node 41 of the solid element 40). ) Is calculated by interpolation with a shape function. As shown in FIG. 32, when the temperature of the node 41 of the solid element 40 is 270 ° C., 280 ° C., 250 ° C., and 240 ° C., the temperature of the solid element at the node 31 of the beam element 30 is 260 ° C. Become. In this case, the temperature Ts of the solid element 40 at the position of the node 31 of the beam element 30 is expressed as follows.

Ts = ΣN _K · Ts, _K (K represents the number of the node 41 of the solid element 40)

Next, the coupled analysis unit 13 calculates the _trial resistance force Fb _trial (S21). The trial resistance force Fb _trial is calculated by the following equation.

Fb _trial = Fb _T-ΔT + K (Ts) · ΔS

Here, K (Ts) is the slope of the elastic region in the temperature-dependent elastic perfect plastic slip resistance model. K (Ts) indicates that K is a function of temperature Ts. FIG. 33 shows a temperature-dependent elastic perfect plastic slip resistance model. As shown in FIG. 33, in the temperature-dependent elastic perfect plastic slip resistance force model, when the slip amount is within a certain range, the resistance force Fb changes in the elastic region. However, when the slip amount S becomes S1 and the resistance force Fb reaches the maximum resistance value Fb _max (Ts), the resistance force Fb is obtained even if the slip amount increases at the same temperature after that, when the slip temperature becomes the same. Does not change. The maximum resistance value Fb _max (Ts) indicates that Fb _max is a function of the temperature Ts. In such a temperature-dependent elastic perfect plastic slip resistance model, the slope K (Ts) of the elastic region is expressed by the following equation.

K (Ts) = Fb _max (Ts) / S1

Next, the coupled analysis unit 13 determines whether the trial resistance force Fb _trial is a resistance force in an elastic region or a resistance force in a complete plastic region in a temperature-dependent elastic perfect plastic slip resistance model. Yield calculation is performed (S22). Yield calculation is performed by the following formula.

| Fb _trial | <Fb _max (Ts)

If the absolute value of the _trial resistance force Fb _trial is smaller than the maximum resistance value Fb _max (Ts) as a result of the yield calculation, the coupled analysis unit 13 determines that the trial resistance force Fb _trial is in the elastic range. It is determined that there is a trial resistance force Fb _trial as the resistance force Fb _T of this step (S23).

However, if the absolute value of the _trial resistance force Fb _trial is greater than or equal to the maximum resistance value Fb _max (Ts) as a result of the yield calculation (S22: NO), the coupled analysis unit 13 completes the _trial resistance force Fb _trial. It determines that the plastic region, determine the resistance force Fb _T by the following equation (S24).

Fb _T = Fb _max (Ts) ・ (Fb _trial / | Fb _trial |)

The coupled analysis unit 13 applies the resistance force Fb _T of this step obtained as described above as a load to the node 31 of the beam element 30 and the node 41 of the solid element 40 (S25).

(Heat conduction / structure coupled analysis)
The calculation was performed using the heat conduction / structure coupled analysis function implemented in the compression molding analysis simulation software of this embodiment. FIGS. 34A, 34B, and 34C are diagrams showing the temperature distribution of the resin model 60 using the heat conduction / structure coupling analysis function in this embodiment step by step in time series.

  In FIGS. 34 (A), (B), and (C), a portion where the temperature of the resin model 60 is high is represented by a light color, and a portion where the temperature is low is represented by a dark color. As shown in FIG. 34A, the temperature of the resin model 60 is high before the start of compression molding. However, as shown in FIGS. 34 (B) and 34 (C), when compression molding is started, the temperature of the resin model 60 decreases from the portion in contact with the male mold 70 and the female mold 80. I understand.

  The behavior of the resin model 60 and the fiber model 50 is calculated by the above-described temperature-dependent elastic perfect plastic slip resistance force model, where the portion where the temperature is low is hard (high viscosity) and the high region is soft (low viscosity). .

  As described above, according to the present embodiment, by making the resistance force dependent on the temperature in the heat conduction / structure coupled analysis, the discontinuous fibers are discontinuous among the discontinuous long fiber reinforced resins when the resin is thermoplastic. The phenomenon of moving while receiving resistance can be accurately reproduced.

(Comparison of compression molding test of discontinuous long fiber reinforced resin and compression molding analysis simulation)
Next, a comparison between a compression molding test of discontinuous long fiber reinforced resin and a compression molding analysis simulation will be described.
35 to 39 are diagrams showing the results of a compression molding test of discontinuous long fiber reinforced resin and the results of a compression molding analysis simulation.
35 (A), (B), (C), and FIGS. 36 (A), (B) are discontinuous long fiber reinforced resins when a compression molding test is performed on the discontinuous long fiber reinforced resin 20. FIG. It is a figure which shows the image which image | photographed 20 shapes in steps according to the time series.
35 (D), (E), (F), and FIGS. 36 (C), (D) show the shape of the resin model 60 by the compression molding analysis simulation of this embodiment corresponding to each stage of the compression molding test. FIG.

  As is clear from FIGS. 35 and 36, the compression molding analysis simulation of the present embodiment accurately represents the shape change of the discontinuous long fiber reinforced resin 20.

  In particular, the rib 24 formed by the discontinuous long fiber reinforced resin 20 flowing in the rib forming portion 73 of the male mold 70, and the rib 61 on the resin model 60 by the compression molding analysis simulation of the present embodiment. It can be seen that the shapes of these are consistent and the molding process is represented with high accuracy.

FIGS. 37A, 37B, and 37C are diagrams showing the formation process of the rib 61 in the fiber model 50 by the compression molding analysis simulation of the present embodiment step by step in time series.
37 (D), (E), and (F) correspond to FIGS. 37 (A), (B), and (C) showing the formation process of the rib 61 in the resin model 60 by the compression molding analysis simulation of the present embodiment. It is a figure shown in steps according to time series.

  As can be understood by comparing FIGS. 37 (A), (B), and (C) with FIGS. 37 (D), (E), and (F), in this embodiment, the above-described constraint coupling calculation is performed. Large deformation analysis of beam element and solid element is performed while performing coupled analysis with sliding resistance force.

  Further, in compression molding, a recess called a weld line is formed, but a weld line 62 is represented on the resin model 60 as shown in FIG. When the place where the weld line 62 is formed is viewed in the fiber model 50, the orientation (orientation) of the fiber can be confirmed as shown by an arrow in FIG.

FIG. 38A is a diagram illustrating an image obtained by photographing the shape of the outer edge portion of the discontinuous long fiber reinforced resin 20 when the compression molding test is performed on the discontinuous long fiber reinforced resin 20. FIG. 38B is a diagram showing the shape of the outer edge portion of the fiber model 50 by the compression molding analysis simulation of the present embodiment.
FIG. 38C is a diagram illustrating an image obtained by photographing a stage in the middle of forming the ribs 24 of the discontinuous long fiber reinforced resin 20 when the compression molding test is performed on the discontinuous long fiber reinforced resin 20. is there. FIG. 38D is a diagram showing the rib 61 on the fiber model 50 in the intermediate stage according to the compression molding analysis simulation of the present embodiment.

  As can be seen by comparing FIG. 38 (A) and FIG. 38 (B), the shape of the wrinkles generated at the outer edge of the actual discontinuous long fiber reinforced resin 20 and the shape of the wrinkles in the fiber model 50 And the real molding state can be accurately reproduced by the compression molding analysis simulation of this embodiment.

  Further, as shown in FIG. 38C, the discontinuous long fiber reinforced resin 20 flows in the direction of the white arrow to form the rib 24, but the rib 61 on the fiber model 50 in FIG. Also, the shape in the middle of formation matches, and it can be seen that the actual molding process can be accurately reproduced by the compression molding analysis simulation of the present embodiment.

  FIG. 39A is a diagram illustrating an image obtained by photographing a cross section of the rib 24 of the discontinuous long fiber reinforced resin 20 when a compression molding test is performed on the discontinuous long fiber reinforced resin 20. FIG. 39B is a view showing a cross section of the rib 61 on the fiber model 50 by the compression molding analysis simulation of the present embodiment.

  As can be clearly understood by comparing FIG. 39A and FIG. 39B, the fiber layer constituting the fiber model 50 is thicker than the fibers of the actual discontinuous long fiber reinforced resin 20, and the number thereof is also Although there are few, it turns out that the orientation (direction) of a fiber and the density of a fiber are expressed with sufficient precision. The fiber model 50 has a cylindrical shape with a set cross-sectional area, which is only drawn to make the analysis result easier to see. In other words, in the analysis, the calculation is performed with the beam element of two nodes.

  As described above, according to the present embodiment, in the compression molding analysis simulation of the discontinuous long fiber reinforced resin, the large deformation analysis of the beam element modeling the discontinuous long fiber and the solid element modeling the resin is performed. In this case, coupled analysis is performed based on the movement of the beam element and the movement of the solid element during the large deformation. Therefore, the behavior of the discontinuous long fibers and the resin, the orientation (direction) of the discontinuous long fibers, and the density of the discontinuous long fibers can be accurately predicted.

  In particular, according to the present embodiment, the movement of the beam element during the large deformation of the beam element and the solid element and the slip amount of the beam element and the solid element based on the movement of the solid element are obtained, and the resistance force is obtained from the obtained slip amount. It calculates and gives resistance as a load to the node of a beam element and a solid element. Therefore, the phenomenon that the discontinuous long fibers move in the resin when the discontinuous long fiber reinforced resin is deformed can be accurately reproduced.

  Further, according to the present embodiment, the temperature of the beam element calculated by interpolating from the temperature calculated at the node of the solid element is obtained, the resistance force is calculated from the slip amount and the temperature, the beam element and Resistive force is applied as a load to the nodes of solid elements. Therefore, even if the resin of the discontinuous long fiber reinforced resin is a thermoplastic resin, the phenomenon that the discontinuous long fibers move in the thermoplastic resin when the thermoplastic resin deforms can be accurately reproduced. Can do.

  The above embodiments are merely examples, and various modifications can be made without departing from the scope of the present invention.

  In the above-described embodiment, the discontinuous long fiber reinforced resin is used as an example of the fiber reinforced resin. The same applies to the case where a resin is used. Moreover, even if it is continuous fiber reinforced resin, it is applicable similarly.

  The fiber is not limited to the discontinuous fiber having the fiber length shown in the embodiment described above, and the present invention can be applied to fibers having any fiber length. For example, it may be a continuous fiber.

  In the above-described embodiment, an analysis simulation when performing compression molding has been described as an example of compression molding. However, the present invention is not limited to such an embodiment, and is also used for analysis simulation when performing press molding. Applicable.

  In the above-described embodiment, the mode in which the resin is modeled by the tetrasolid element having a four-node tetrahedron has been described. However, the present invention is not limited to such a mode, and the solid element other than the tetrasolid element is also used. Applicable.

  The compression molding analysis system program according to the above aspect can be provided in a form stored in a computer-readable recording medium and installed in the computer. The recording medium is, for example, a non-transitory recording medium, and an optical recording medium such as a CD-ROM is a good example, but a known arbitrary format such as a semiconductor recording medium or a magnetic recording medium is used. A recording medium may be included. It is also possible to provide the aforementioned program in the form of distribution via a communication network and install it on a computer.

  As described above, the compression molding analysis system, the compression molding analysis method, and the compression molding analysis program according to the embodiment of the present invention have been described, but the present invention is not limited to this, and does not depart from the gist of the present invention. Various modifications are possible.

DESCRIPTION OF SYMBOLS 11 Fiber model part 12 Resin model part 13 Coupled analysis part 14 Large deformation analysis part 100 Compression molding analysis system

Claims (12)

A compression molding analysis system for fiber reinforced resin in which fibers are mixed with resin,
A fiber model unit that models the fiber as a fiber model with a beam element;
A resin model portion that models the resin as a resin model with a solid element so as to have a node that is not shared with the node of the beam element;
A coupled analysis unit that performs coupled analysis of the movement of the beam element and the movement of the solid element;
A large deformation analysis unit for performing large deformation analysis on the beam element and the solid element,
The coupled analysis unit, on the basis of the motion of said beam elements in the large deformation of the beam element and the solid element by the large deformation analysis unit and the motion of the solid element, have rows the coupled analysis,
The coupled analysis unit includes the beam element for the beam element and the solid element based on the movement of the beam element and the movement of the solid element during the large deformation of the beam element and the solid element by the large deformation analysis unit. The relative displacement in the axial direction is acquired as a slip amount, and the resistance force is calculated from the acquired slip amount,
The large deformation analysis unit applies the resistance force as a load to the nodes of the beam element and the solid element.
A compression molding analysis system characterized by this. A compression molding analysis system for fiber reinforced resin in which fibers are mixed with resin,
A fiber model unit that models the fiber as a fiber model with a beam element;
A resin model portion that models the resin as a resin model with a solid element so as to have a node that is not shared with the node of the beam element;
A coupled analysis unit that performs coupled analysis of the movement of the beam element and the movement of the solid element;
A large deformation analysis unit for performing large deformation analysis on the beam element and the solid element,
The coupled analysis unit performs the coupled analysis based on the movement of the beam element and the movement of the solid element during the large deformation of the beam element and the solid element by the large deformation analysis unit,
The coupled analysis unit includes the beam element for the beam element and the solid element based on the movement of the beam element and the movement of the solid element during the large deformation of the beam element and the solid element by the large deformation analysis unit. The slip amount, which is a relative displacement in the axial direction, and the temperature at the node position of the beam element calculated by interpolation from the temperature calculated at the node of the solid element are acquired, and the acquired slip Calculate the resistance from the amount and the temperature,
The large deformation analysis unit applies the resistance force as a load to the nodes of the beam element and the solid element .
Compression forming analysis system that is characterized in that. The large deformation analysis unit performs large deformation analysis on the solid element by an adaptive remeshing method.
The compression molding analysis system according to claim 1 or 2, characterized by the above. The coupled analysis unit performs coupled analysis between the movement of the node of the beam element and the movement of the node of the solid element by a constraint method in the direction perpendicular to the axial direction of the beam element, As for the axial direction, a coupled analysis is performed by calculating the resistance force and applying the resistance force as a load to the node of the beam element and the node of the solid element.
Compression molding analysis system according to any one of claims 1 to 3, characterized in that. The fiber model unit sets the material characteristics of the fiber measured from the real to the beam element, and the fiber volume fraction contained in the resin model of the fiber model is the same as the fiber volume fraction of the real thing. Setting the cross-sectional area of the beam element;
The resin model part sets a value obtained by multiplying the solid material by measuring (1-fiber volume fraction) the resin material property measured from the real thing.
Compression molding analysis system according to any one of claims 1 to 4, characterized in that. The fiber model unit models the fiber with a beam element of two nodes having cross-sectional shape information,
The resin model unit models the resin with a solid element of a four-node tetrahedron.
The compression molding analysis system according to any one of claims 1 to 5, wherein The fiber can be set to any length fiber length,
The compression molding analysis system according to any one of claims 1 to 6, wherein the compression molding analysis system is characterized in that: The resin is a thermoplastic resin or a thermosetting resin.
The compression molding analysis system according to any one of claims 1 to 7, wherein the compression molding analysis system is characterized in that: A compression molding analysis method for a fiber reinforced resin in which fibers are mixed with a resin,
The fiber model unit models the fiber as a fiber model with a beam element,
With the resin model part, the resin is modeled as a resin model with a solid element so as to have a node that is not shared with the node of the beam element,
The large deformation analysis unit performs large deformation analysis on the beam element and the solid element,
Based on the movement of the beam element and the movement of the solid element based on the movement of the beam element and the movement of the solid element during the large deformation of the beam element and the solid element, the coupled analysis unit And
The coupled analysis unit causes the beam element and the solid element based on the movement of the beam element and the solid element during the large deformation of the beam element and the solid element to be relative to each other in the axial direction of the beam element. A large amount of displacement as a slip amount, calculate a resistance force from the acquired slip amount,
The large deformation analysis unit applies the resistance force as a load to the nodes of the beam element and the solid element.
A compression molding analysis method characterized by the above . A compression molding analysis method for a fiber reinforced resin in which fibers are mixed with a resin,
The fiber model unit models the fiber as a fiber model with a beam element,
With the resin model part, the resin is modeled as a resin model with a solid element so as to have a node that is not shared with the node of the beam element,
The large deformation analysis unit performs large deformation analysis on the beam element and the solid element,
Based on the movement of the beam element and the movement of the solid element based on the movement of the beam element and the movement of the solid element during the large deformation of the beam element and the solid element, the coupled analysis unit And
The coupled analysis unit causes the beam element and the solid element based on the movement of the beam element and the solid element during the large deformation of the beam element and the solid element to be relative to each other in the axial direction of the beam element. A slip amount that is a large displacement and a temperature at the node position of the beam element calculated by interpolation from the temperature calculated at the node of the solid element, and resistance from the acquired slip amount and the temperature Calculate the force,
The large deformation analysis unit applies the resistance force as a load to the nodes of the beam element and the solid element .
A compression molding analysis method characterized by the above . A compression molding analysis program for causing a computer to perform compression molding analysis of a fiber reinforced resin obtained by mixing fibers with resin.
Modeling the fiber as a fiber model with a beam element;
Modeling the resin as a resin model with solid elements so as to have nodes that are not shared with the nodes of the beam elements;
Performing a large deformation analysis on the beam element and the solid element;
Performing a coupled analysis of the beam element movement and the solid element movement based on the beam element movement and the solid element movement during the large deformation of the beam element and the solid element;
The relative displacement in the axial direction of the beam element is obtained as a slip amount for the beam element and the solid element based on the movement of the beam element and the movement of the solid element during the large deformation of the beam element and the solid element. And calculating a resistance force from the acquired slip amount;
Applying the resistive force as a load to the nodes of the beam element and the solid element,
A compression molding analysis program characterized by that. A compression molding analysis program for causing a computer to perform compression molding analysis of a fiber reinforced resin obtained by mixing fibers with resin.
Modeling the fiber as a fiber model with a beam element;
Modeling the resin as a resin model with solid elements so as to have nodes that are not shared with the nodes of the beam elements;
Performing a large deformation analysis on the beam element and the solid element;
Performing a coupled analysis of the beam element movement and the solid element movement based on the beam element movement and the solid element movement during the large deformation of the beam element and the solid element;
A sliding amount that is a relative displacement in the axial direction of the beam element with respect to the beam element and the solid element based on the movement of the beam element during the large deformation of the beam element and the solid element and the movement of the solid element; Obtaining a temperature at the position of the node of the beam element calculated by interpolation from the temperature calculated at the node of the solid element, and calculating a resistance force from the acquired amount of the slip and the temperature;
Applying the resistive force as a load to the nodes of the beam element and the solid element,
A compression molding analysis program characterized by that.