JP2016151480A - Composite material analysis method, computer program for analyzing composite material, evaluation method of analysis result for composite material, and computer program for evaluating analysis result for composite material - Google Patents

Abstract

PROBLEM TO BE SOLVED: To provide a composite material analysis method capable of analyzing local change in a molecular state of polymer molecules, a computer program for analyzing the composite material, an evaluation method of an analysis result for the composite material, and a computer program for evaluating the analysis result for the composite material.SOLUTION: A composite material analysis method includes: a first step of preparing an analysis model of composite material containing a polymer model and a filler model; a second step of cross-linking the polymer model by performing cross-linking analysis; a third step of forming segments from at least one kind selected from the group consisting of polymer atoms, polymer particles and bond chains, contained in the polymer model after cross-linking, and creating a segment list in which a plurality of segments as analysis targets are specified; a fourth step of setting interaction in the analysis model of the composite material and analyzing a molecular state of the polymer model for each specified segment by numerical analysis.SELECTED DRAWING: Figure 6

Description

  The present invention relates to a composite material analysis method, a composite material analysis computer program, a composite material analysis result evaluation method, and a composite material analysis result evaluation computer program. The present invention relates to a composite material analysis method capable of analysis, a composite material analysis computer program, a composite material analysis result evaluation method, and a composite material analysis result evaluation computer program.

  Conventionally, a method for creating a model of a composite material including a modified polymer and a filler used in an automobile tire or the like has been proposed (see, for example, Patent Document 1). In this model creation method, the interaction between the modified polymer and the filler is made larger than the interaction between other particles, and the modified polymer and the filler are dispersed in the composite material. In addition, the interaction between the modified polymer and the filler is made smaller than the interaction between other particles, and the end of the modified polymer reacts with the filler to enable analysis of the material properties of the composite material. Create a model.

JP 2012-177609 A

  By the way, in a composite material containing two or more kinds of substances such as rubber containing a polymer and a filler, it is predicted that the molecular motion of the polymer existing around the filler has a great influence on the material properties of the compound. In order to elucidate the mechanism of expression of the material properties of such a composite material, it is desirable to accurately analyze changes in molecular state such as molecular orientation and molecular spread accompanying compound deformation.

  However, in the analysis using the conventional composite material analysis model, since the movement of the polymer is evaluated for each large polymer molecule after crosslinking, only the change in the molecular state of the entire polymer molecule after crosslinking can be evaluated. For this reason, in the conventional analysis, for example, the local state change of the polymer molecule is accurately analyzed, such as the state change of the region near the filler in the polymer molecule and the state change of the polymer molecular state in the region away from the filler. There are circumstances that cannot be evaluated.

  The present invention has been made in view of such circumstances, and a composite material analysis method capable of analyzing and evaluating a local molecular state change of a polymer molecule, a composite material analysis computer program, and a composite material An object is to provide an analysis result evaluation method and a computer program for evaluation of an analysis result of a composite material.

  The composite material analysis method of the present invention is a composite material analysis method using a composite material analysis model created by a molecular dynamics method using a computer, and includes analysis of a composite material including a polymer model and a filler model. At least one selected from the group consisting of a polymer atom, a polymer particle, and a bond chain included in the polymer model after crosslinking, a first step of creating a model for use, a second step of crosslinking the polymer model by crosslinking analysis, and A third step of creating a segment list in which a seed is a segment and a plurality of segments to be analyzed is specified, and the polymer model for each segment specified by numerical analysis by setting an interaction in the analysis model of the composite material And a fourth step of analyzing the molecular state of

  According to the composite material analysis method of the present invention, it is possible to perform numerical analysis according to the interaction set for each segment of the polymer model after cross-linking. Therefore, the molecular state of the polymer model is analyzed for each segment of the polymer model. It becomes possible to do. As a result, the composite material analysis method can analyze the change in the molecular state of the polymer model in the region near the filler model and the molecular state of the polymer model in the region away from the filler model. Analysis of local molecular state changes of molecules becomes possible.

  The computer program for analyzing a composite material according to the present invention causes a computer to execute the above-described composite material analysis method.

  According to the computer program for analyzing a composite material of the present invention, it is possible to perform a numerical analysis according to the interaction set for each segment of the polymer model after crosslinking, so that the molecular state of the polymer model is determined for each segment of the polymer model. Can be analyzed. As a result, the computer program for analyzing the composite material can analyze the change in the molecular state of the polymer model in the region near the filler model and the molecular state of the polymer model in the region away from the filler model, respectively. This makes it possible to analyze local changes in the molecular state of polymer molecules.

  The analysis method of the analysis result of the composite material of the present invention is an evaluation method of the analysis result of a composite material analyzed using a model for analysis of a composite material created by a molecular dynamics method using a computer, the polymer model and A first step of creating a model for analyzing a composite material including a filler model, a second step of cross-linking the polymer model by cross-linking analysis, and polymer atoms, polymer particles, and bonding chains included in the polymer model after cross-linking A third step of creating a segment list in which at least one selected from the group consisting of a segment and a plurality of segments to be evaluated is specified, and an analysis model of a composite material including a polymer model and a filler model are interacted with Set and solve the molecular state of the polymer model for each segment specified by numerical analysis. Characterized in that it comprises a fourth step for to evaluate the.

  According to the evaluation method of the analysis result of the composite material of the present invention, numerical analysis corresponding to the interaction set for each segment of the polymer model after crosslinking is possible. It is possible to analyze and evaluate the state. Thereby, the evaluation method of the analysis result of the composite material can analyze the change in the molecular state of the polymer model in the region near the filler model and the molecular state of the polymer model in the region away from the filler model, respectively. Therefore, it is possible to evaluate the local molecular state change of the polymer molecule.

  In the evaluation method of the analysis result of the composite material of the present invention, it is preferable that in the fourth step, a numerical value representing a state change of a specific segment included in the segment list is calculated. By this method, it is possible to quantitatively evaluate the state change of a specific segment included in the segment list.

  In the composite material analysis result evaluation method of the present invention, in the third step, a distance between the polymer model and the filler model is added to the segment list, and in the fourth step, the polymer model is added. It is preferable to evaluate the relationship between the distance between the filler model and the numerical value representing the state change of a specific segment included in the segment list. This method makes it possible to evaluate the relationship between the distance from the filler model and the state change of a specific segment of the polymer model.

  In the composite material analysis result evaluation method of the present invention, in the fourth step, a numerical value representing a state change of a specific segment included in the segment list is calculated in a plurality of analysis times, and the calculated numerical value and the analysis time are calculated. It is preferable to evaluate the relationship. By this method, the relationship between a specific segment that changes at each analysis time and a physical quantity can be evaluated, so that a more detailed analysis of a local molecular state change of a polymer molecule becomes possible.

  In the evaluation method of the analysis result of the composite material of the present invention, in the fourth step, it is preferable to evaluate the relationship between the numerical value representing the state change of the segment and the strain or analysis time. By this method, it is possible to evaluate the relationship between bond length of polymer molecules and polymer particle velocity and strain or analysis time in a specific segment that changes every analysis time, so that more detailed analysis of local molecular state changes of polymer molecules Is possible.

  In the analysis method of the analysis result of the composite material of the present invention, in the fourth step, in the fourth step, the transformation analysis is performed to evaluate the relationship between the numerical value representing the state change of the segment and the pressure or the analysis time. It is preferable. By this method, it is possible to evaluate the relationship between the bond length of polymer molecules and the velocity of polymer particles in a specific segment that changes at each analysis time, and the pressure or analysis time. Is possible.

  In the evaluation method of the analysis result of the composite material of the present invention, it is preferable that in the fourth step, the temperature change analysis is executed to evaluate the relationship between the numerical value indicating the state change of the segment and the temperature or the analysis time. This method can evaluate the relationship between polymer molecule bond length and polymer particle velocity and temperature or analysis time in a specific segment that changes with analysis time, allowing more detailed analysis of local molecular state changes of polymer molecules. Is possible.

  In the evaluation method of the analysis result of the composite material of the present invention, it is preferable to visualize a specific segment included in the segment list based on the numerical analysis result in the fourth step. By this method, since a specific segment is visualized, the state change of the segment of the characteristic of the polymer model can be drawn and evaluated. For example, by continuously displaying the visualized result for each analysis time, it is possible to evaluate the state change of a specific segment as a moving image.

  The computer program for evaluating the analysis result of the composite material according to the present invention causes the computer to execute the evaluation method of the analysis result of the composite material.

  According to the computer program for evaluating the analysis result of the composite material of the present invention, it is possible to perform numerical analysis according to the interaction set for each segment of the polymer model after cross-linking. It is possible to analyze and evaluate the molecular state. As a result, the computer program for evaluating the analysis result of the composite material can analyze the change in the molecular state of the polymer model in the region near the filler model and the molecular state of the polymer model in the region away from the filler model, respectively. Therefore, it is possible to realize a computer program for evaluating the analysis result of the composite material capable of evaluating the local molecular state change of the polymer molecule.

  According to the present invention, a composite material analysis method capable of analyzing and evaluating local molecular state changes of polymer molecules, a composite material analysis computer program, a composite material analysis result evaluation method, and a composite material analysis result The evaluation computer program can be realized.

FIG. 1 is a conceptual diagram showing an example of an analysis model used in the composite material analysis method according to the present embodiment. FIG. 2A is a diagram illustrating an example in which polymer particles in a region near a filler model in a polymer model molecule after crosslinking are visualized as segments. FIG. 2B is a diagram showing an example in which polymer particles in a region near a filler model in a polymer model molecule after crosslinking are visualized as segments. FIG. 2C is a diagram showing an example in which polymer particles in a region near a filler model in a polymer model molecule after crosslinking are visualized as segments. FIG. 3 is a diagram showing an example in which polymer particles in a region in the vicinity of the filler model in the polymer model molecule after crosslinking are visualized as a segment after extension analysis. FIG. 4A is a diagram showing an example in which polymer particles in a region away from a filler model in a polymer model molecule after crosslinking are visualized as segments. FIG. 4B is a diagram showing an example in which polymer particles in a region away from the filler model in the polymer model molecule after crosslinking are visualized as segments. FIG. 4C is a diagram showing an example in which polymer particles in a region away from the filler model in the polymer model molecule after crosslinking are visualized as segments. FIG. 5 is a diagram showing an example in which polymer particles in a region away from a filler model in a polymer model molecule after crosslinking are visualized as a segment after extension analysis. FIG. 6 is a flowchart of the evaluation method of the analysis result of the composite material according to the present embodiment. FIG. 7A is a conceptual diagram of an uncrosslinked polymer model. FIG. 7B is a conceptual diagram of a polymer model after crosslinking. FIG. 8A is a diagram illustrating an example of a segment. FIG. 8B is a diagram illustrating an example of a segment. FIG. 8C is a diagram illustrating an example of a segment. FIG. 8D is a diagram illustrating an example of a segment. FIG. 9 is a functional block diagram of an analysis apparatus that executes a composite material analysis result evaluation method and a composite material simulation method according to an embodiment of the present invention. FIG. 10 is a diagram illustrating the relationship between the distance from the filler and the inertial radius.

  Hereinafter, an embodiment of the present invention will be described in detail with reference to the accompanying drawings. In addition, this invention is not limited to the following embodiment, It can implement by changing suitably.

  The composite material analysis method according to the present embodiment is a composite material analysis method using a composite material analysis model created by a molecular dynamics method using a computer, and includes a polymer model and a filler model. A first step of creating a model for analyzing the material, a second step of cross-linking the polymer model by cross-linking analysis, and at least selected from the group consisting of polymer atoms, polymer particles and bonding chains included in the polymer model after cross-linking The third step of creating a segment list that identifies a plurality of segments to be analyzed, with one type as a segment, and the polymer model for each segment identified by numerical analysis by setting an interaction in the composite material analysis model A fourth step of analyzing the molecular state. The evaluation method of the analysis result of the composite material according to the present embodiment is an evaluation method of the local molecular state change of the polymer molecule using the analysis method of the composite material according to the present embodiment. First, an outline of a composite material analysis model used in the composite material analysis result evaluation method according to the present embodiment will be described.

  FIG. 1 is a conceptual diagram showing an example of an analysis model 1 used in the evaluation method of the analysis result of the composite material according to the present embodiment. As shown in FIG. 1, the composite material analysis model 1 according to the present embodiment includes two filler models 11 in a model creation region A, which is a substantially cubic virtual space whose one side is a distance L1. A composite material to be analyzed in which a plurality of polymer models 12 are arranged is modeled. Note that in this embodiment, an example in which a composite material to be analyzed contains a filler and a polymer will be described. However, the present invention can be applied to various composite materials in addition to a composite material containing two kinds of substances. is there.

  In the composite material analysis model 1, examples of the filler include carbon black, silica, and alumina. The filler model 11 is modeled as a substantially spherical body in which a plurality of filler atoms and filler particles 11a as an aggregate of a plurality of filler atoms are aggregated. The relative position of the filler particles 11a is specified by a bond chain (not shown) between the plurality of filler particles 11a. This binding chain (not shown) has a function as a spring in which an equilibrium length, which is a binding distance between filler particles 11a, and a spring constant are defined, and binds between the filler particles 11a. The bond chain is a bond in which the relative position of the filler particle 11a and the potential at which force is generated by twisting, bending, or the like are defined. The filler model 11 is numerical data (including the mass, volume, diameter, initial coordinates, and the like of the filler particles 11a) for handling the filler by molecular dynamics. Numerical data of the filler model 11 is input to the computer.

  Examples of the polymer include rubber, resin, and elastomer. A modifier is added to the polymer as needed to enhance the affinity with the filler. Examples of the modifier include a hydroxyl group, a carbonyl group, and a functional group of an atomic group. The polymer model 12 is modeled by filling a plurality of polymer atoms and polymer particles 12a, which are aggregates of a plurality of polymer atoms, into the model creation region A at a predetermined density. The relative position of the polymer particle 12a is specified by the bonding chain 12b between the plurality of polymer particles 12a. The binding chain 12b functions as a spring in which an equilibrium length, which is a binding distance between the polymer particles 12a, and a spring constant are defined, and restrains the polymer particles 12a. The bond chain 12b is a bond in which the relative position of the polymer particle 12a and the potential at which force is generated by twisting, bending, or the like are defined. The polymer model 12 is numerical data (including the mass, volume, diameter, initial coordinates, and the like of the polymer particle 12a) for handling the polymer by molecular dynamics. Numerical data of the polymer model 12 is input to a computer.

  Next, with reference to FIG. 2A to FIG. 5, the evaluation method of the analysis result of the composite material according to the present embodiment will be described in detail. In the evaluation method of the analysis result of the composite material according to the present embodiment, a part of the polymer model 12 molecule after crosslinking is specified as a segment of the polymer model 12 and analyzed and evaluated. FIG. 2A to FIG. 2C are diagrams showing examples in which polymer particles 12a in the vicinity of the filler model in the polymer model 12 molecule after crosslinking are visualized as segments 121, and FIG. 3 is a diagram in the polymer model 12 molecule after crosslinking. It is a figure which shows the example visualized after the expansion | extension analysis as the polymer particle 12a of the area | region of the filler model 11 vicinity as the segment 121. FIG. FIGS. 4A to 4C are diagrams showing an example in which polymer particles 12a in a region away from the filler model 11 in the polymer model 12 molecule after crosslinking are visualized as segments 122, and FIG. 5 is a diagram illustrating the polymer model 12 after crosslinking. It is a figure which shows the example visualized after the extension analysis as the polymer particle 12a of the area | region away from the filler model 11 in a molecule | numerator as the segment 122. FIG. 2A and 4A show the state of the filler particles 11a and the polymer particles 12a viewed from the X direction, and in FIGS. 2B and 4B, the state of the filler particles 11a and the polymer particles 12a viewed from the Y direction. 2C, FIG. 3, FIG. 4C, and FIG. 5 show the state of the filler particles 11a and the polymer particles 12a as viewed from the Z direction.

  In the example shown in FIGS. 2A to 2C, the polymer model 12 is a region in the vicinity of the two filler models 11, which is a part of the entire polymer model 12 molecule in which a plurality of polymer models 12 are cross-linked by a bonding chain 12 b. Only the segment 121 of the polymer particle 12a is visualized. As shown in FIG. 3, the distance between the two filler models 11 is expanded after the extension analysis, and the segment 121 of the polymer particle 12 a in the region in the vicinity of the two filler models 11 interacts with the filler model 11. As a result, the state is gathered in a region between the two filler models 11 in the vicinity of the filler model 11.

  In the example shown in FIGS. 4A to 4C, the polymer model 12 is a region away from the two filler models 11 that are a part of the entire polymer model 12 molecule in which a plurality of polymer models 12 are cross-linked by a bonding chain 12 b. Only the segment 121 of the polymer particle 12a is visualized. As shown in FIG. 5, when the extension analysis is performed in the same manner as in the example shown in FIG. 3, the segment 122 of the polymer particle 12a in the region away from the two filler models 11 has an interaction with the filler shown in FIG. It is relatively small with respect to the example shown in (2) and is in a state of being distributed between the two filler models 11.

  Thus, in the present embodiment, a part of the polymer particles 12a in the entire large polymer model 12 molecule in which a plurality of polymer models 12 are cross-linked by the bonding chains 12b are analyzed as the segments 121 and 122, so that the filler model It is possible to evaluate changes in the orientation and spread of some molecular states in the entire molecule of the polymer model 12 due to the interaction depending on the distance between the polymer model 12 and the polymer model 12. As a result, the evaluation method of the analysis result of the composite material according to the present embodiment accurately analyzes and evaluates the change in the molecular state of the polymer accompanying the deformation of the compound, which contributes to the influence on the material properties of the compound. It becomes possible to do. In the above-described embodiment, the example in which the polymer particle 12a is specified as a segment has been described. However, the bonding chain 12b may be specified as a segment.

  Next, a composite material analysis method and a composite material analysis result evaluation method according to the present embodiment will be described in more detail with reference to FIG. FIG. 6 is a flowchart of the evaluation method of the analysis result of the composite material according to the present embodiment.

  As shown in FIG. 6, in the first step, an uncrosslinked polymer model 12 is created and a filler model 11 is created (step ST11 and step ST12). FIG. 7A is a conceptual diagram of the uncrosslinked polymer model 12. As shown in FIG. 7A, the uncrosslinked polymer models 12-1 and 12-2 are formed by connecting a plurality of polymer particles 12a with a bonding chain 12b. Next, as shown in FIG. 1 in these filler models 11, polymer models 12-1 and 12-2 are arranged to create a composite material analysis model 1 (step ST13).

  Next, in the second step, cross-linking points are created by cross-linking analysis of the polymer models 12-1 and 12-2 in the composite material analysis model 1 (step ST14). FIG. 7B is a conceptual diagram of polymer models 12-1 and 12-2 after crosslinking. As shown in FIG. 7B, in the crosslinking analysis, predetermined polymer particles 12a-1 and 12a-2 in the uncrosslinked polymer models 12-1 and 12-2 are specified to create a crosslinking point 12b-1. As a result, in the composite material analysis model 1, a large polymer model 12 after crosslinking is created by the plurality of polymer models 12. Here, the term “crosslinking” refers to forming the polymer model 12 including the polymer particles 12a formed by connecting three or more bonding chains 12b.

  Next, in the third step, at least one selected from the group consisting of the polymer atom of the polymer model 12 after cross-linking, the polymer particle 12a and the bonding chain 12b is defined as a segment, and is specified for each of a plurality of segments to be evaluated. Create a segment list. Next, a polymer molecule list for each evaluation segment to be evaluated is created (step ST15). Next, a segment to be evaluated is selected from the prepared polymer molecule list (step ST16). 8A to 8C are diagrams illustrating examples of segments. As shown in FIG. 8A, the segments may specify the entire polymer models 12-1 and 12-2 before cross-linking as one segment S1 and S2, respectively. Further, as shown in FIG. 8B, the segments are identified as four segments S3 to S6 between the crosslinking points 12b-1 and 12b-2 of the polymer models 12-1 and 12-2 after crosslinking. May be. Moreover, as shown to FIG. 8C, a segment may specify the range containing the crosslinking points 12b-1 and 12b-2 of the polymer models 12-1 and 12-2 after crosslinking as six segments S7 to S12. . Further, as shown in FIG. 8D, the segment may be specified as a segment S13 including a crosslinking point 12a-5 and a free end 12c of the polymer model 12-3 after crosslinking, and only the free end is specified as a segment. Also good. The free end means that only one end of the segment is a crosslinking point. Further, in the segment, the crosslinking point 12a-6 and the bonding chain 12d bonded to the crosslinking point 12a-6 in a loop shape may be specified as the segment S14, or only the bonding chain 12d may be specified as the segment. The segment list includes at least one of a segment number, a polymer particle number constituting the segment, and a bond chain number. Only the particle number or the bond number may be used.

  Next, in the fourth step, an interaction is set in the composite material analysis model including the polymer model 12 and the filler model 11 (step ST17). As the interaction here, for example, the interaction between the filler particles 11a, the polymer particles 12a, the interaction between the filler particles 11a and the polymer particles 12a, and the state in which the filler particles 11a and the polymer particles 12a are bonded with a bond chain. Interactions can be mentioned, but not all of these need to be set. Further, when the polymer model 12 includes a plurality of types of polymer models, the polymer model 12 may be set between the polymer particles 12a, the filler particles 11a, and the polymer particles 12a constituting each polymer model 12. Furthermore, in this case, depending on the type of the polymer model 12 to be created, for example, the first interaction between the polymer particles 12aA and the filler particles 11a constituting the first polymer model 12 and the second polymer The second interaction between the polymer particles 12aB constituting the model 12 and the filler particles 11a may be set as a different interaction.

  Next, numerical analysis is performed to analyze and evaluate the molecular state of the polymer model 12 for each identified segment (step ST18). As the numerical analysis here, for example, a physical quantity for each analysis time that changes for each analysis time is evaluated. Thereby, since the physical quantity of the whole analysis model which changes for every analysis time can be evaluated, a more detailed analysis of a local molecular state change of a polymer molecule becomes possible. Examples of numerical analysis include deformation analysis, transformation analysis, and temperature change analysis. This numerical analysis represents changes in the state of segments such as the bond length and polymer particle speed of the entire analysis model, which varies with the analysis time, the speed or bond length between the crosslinking points and the free ends, and physical quantities such as orientation. The relationship between the numerical value and strain, the bond length of the polymer molecule that changes with each analysis time and the state of the segment such as the polymer particle velocity and the pressure or analysis time, and the polymer molecule that changes with the analysis time Since the relationship between the numerical value representing the state change of the segment such as the bond length and the polymer particle velocity and the temperature or the analysis time can be evaluated, a more detailed analysis of the local molecular state change of the polymer molecule becomes possible. Here, evaluation may be made by taking an average of a plurality of polymer molecule segments. In this case, for example, the average of the free end segments of the polymer model 12 may be evaluated, and a part or all of the segments between the cross-linking points of the entire polymer model 12 may be evaluated.

  Finally, the segment of the selected polymer model 12 is evaluated, and the evaluation of the analysis result is finished (step ST19). Here, for example, among the segments of the polymer model 12 in the polymer molecule list for each evaluation segment, for example, the region near the filler model 11 of the polymer model 12 after crosslinking is separated from the segment 121 and the filler model 11 of the polymer particle 12a. By evaluating the segments 122 of the polymer particles 12a selected from the regions, the regions of the polymer particles 12a of the polymer model 12 in the region near the filler model 11 and the region away from the filler model 11 shown in FIGS. Each molecular motion can be evaluated.

  In the fourth step, it is possible to calculate a numerical value representing a state change of a specific segment included in the created segment list. Examples of the numerical value calculated here include the orientation, inertia radius, mean square displacement, end-to-end distance, average bond length, acceleration, and velocity of the polymer particles 12a included in the segment list. By calculating the numerical value in this way, it becomes possible to quantitatively evaluate the state change of the polymer particles 12a belonging to a specific segment included in the segment list.

  In the fourth step, the distance between the polymer model 12 and the filler model 11 calculated in advance in the third step is added to the segment list, and the distance between the polymer model 12 and the filler model 11 and the segment list are added. You may evaluate the relationship with the numerical value showing the state change of the specific segment contained in. In this case, for example, the relationship between the distance between the polymer model 12 and the filler model 11 and the numerical value indicating the state change of a specific segment included in the segment list may be evaluated as a table. You may evaluate as. FIG. 10 shows the relationship between the distance from the filler and the radius of inertia as an example of a numerical value representing the state change of a specific segment included in the segment list. In the example shown in FIG. 10, it can be seen that the radius of inertia increases in a range where the distance from the filler is short, but rapidly decreases when a certain distance from the filler leaves. Thus, by evaluating, it becomes possible to evaluate the relationship between the distance from a filler, and the state change of the specific segment of a polymer.

  Here, the radius of inertia when the analysis model 1 of the composite material including the polymer model 12 and the filler model 11 after crosslinking is expanded by 200% by extension analysis will be described with reference to Table 1 below. In Example 1 in Table 1 below, the radius of inertia was evaluated using one molecule of the polymer model 12 in a state in the vicinity of the filler model 11 included in the large one molecule polymer model 12 after crosslinking as a segment. An example is shown. Example 2 in Table 1 below is an example in which the radius of inertia was evaluated using one molecule of the polymer model 12 in a state before the crosslinking in a region away from the filler model 11 included in the large one molecule polymer model 12 after crosslinking as a segment. Is shown. The comparative examples in Table 1 below show examples in which the radius of inertia of the entire polymer model 12 having a large molecule after crosslinking was evaluated.

  As can be seen from Table 1, when the inertial radius was evaluated using one molecule of the polymer model 12 in the state before crosslinking as a segment, the value of the inertial radius before and after deformation by extension analysis was a large molecule of polymer after crosslinking. It can be seen that there is a large difference from the case where the inertia radius of the entire model 12 is evaluated (Examples 1 and 2 and Comparative Example). And the value of the inertial radius before and after deformation by extension analysis is the case where one molecule of the polymer model 12 in a state before cross-linking in a region away from the filler model 11 is a segment (Example 1), and in the vicinity of the filler model 11 It can be seen that there is a great difference between the case where one molecule of the polymer model 12 in the state before crosslinking of the region is segmented (Example 2). From these results, when the radius of inertia of the entire polymer model 12 after crosslinking is evaluated, only the radius of inertia of the entire molecule of the large polymer model 12 after crosslinking accompanying the deformation of the entire model 1 for analysis of a composite material can be obtained. It turns out that it has not been evaluated. Then, by evaluating the radius of inertia using one molecule of the polymer model 12 before cross-linking as a segment, the local spread (deformation) of the polymer molecule corresponding to the distance between the filler model 11 and the polymer model 12 is evaluated. I understand that I can do it.

  In the composite material analysis result evaluation method according to the present embodiment, it is preferable to visualize a specific segment included in the segment list based on the numerical analysis result in the fourth step. When a specific segment is visualized, only the specific segment may be visualized and displayed, or the specific segment may be displayed in a different color from the other segments. Moreover, when visualizing a specific segment, you may visualize and display the said specific segment with different materials, such as a filler and another polymer. When a specific segment is visualized, contour display may be performed according to a numerical value representing a state change of the specific segment included in the calculated segment list. In this way, by visualizing a specific segment, it is possible to describe and evaluate the state change of the polymer particle 12a. Further, by visualizing and continuously displaying the state change of the polymer particle 12a for each analysis time, it is possible to make the state change of the polymer particle 12a into a moving image.

  Next, a method for evaluating the analysis result of the composite material and a computer program for evaluating the analysis result of the composite material according to the present embodiment will be described in detail. FIG. 9 is a functional block diagram of an analysis apparatus that executes a composite material analysis result evaluation method and a composite material simulation method according to an embodiment of the present invention.

  As shown in FIG. 9, the analysis method 50 and the composite material simulation method according to the present embodiment are realized by an analysis device 50 that is a computer including a processing unit 52 and a storage unit 54. This analysis device 50 is electrically connected to an input / output device 51 having an input means 53. The input means 53 processes various physical property values of the polymer and filler for which a composite material analysis model is to be created, the actual measurement result of the extension test result using the composite material containing the polymer and filler, and the boundary conditions in the analysis. Input to the unit 52 or the storage unit 54. As the input means 53, for example, an input device such as a keyboard and a mouse is used.

  The processing unit 52 includes, for example, a central processing unit (CPU: CentralL1 Processing Unit) and a memory. The processing unit 52 reads a computer program from the storage unit 54 and develops it in a memory when executing various processes. The computer program expanded in the memory executes various processes. For example, the processing unit 52 expands data relating to various processes stored in advance from the storage unit 54 to an area allocated to itself on the memory as necessary, and analyzes the composite material based on the expanded data. Various processes related to the simulation of the composite material using the model creation and the composite material analysis model are executed.

  The processing unit 52 includes a model creation unit 52a, a condition setting unit 52b, and an analysis unit 52c. The model creation unit 52a is based on data stored in the storage unit 54 in advance, and the number of particles, the number of molecules, and the molecular weight of the composite material such as filler and polymer when creating the composite material analysis model 1 by the molecular dynamics method. , Molecular chain length, number of molecular chains, branching, shape, size, reaction time, reaction conditions, arrangement of components such as the number of molecules included in the model for analysis to be created, setting and number of calculation steps, etc. The initial conditions of various calculation parameters such as the coarse-grained model and the interaction between molecular chains are set.

As calculation parameters for adjusting the interaction between the filler particles 11a and the interaction between the polymer particles 12a, σ and ε of Leonard-Jones potential expressed by the following formula (1) are used, and these are adjusted. By increasing the upper limit distance (cutoff distance) for calculating the potential, it is possible to adjust the attractive force and repulsive force that worked to a long distance. It is preferable that the interaction parameter between the filler particles 11a and the interaction between the polymer particles 12a are sequentially reduced until the interaction between the filler particles 11a and the interaction between the polymer particles 12a reach a constant value. By gradually bringing the σ and ε of the Leonard-Jones potential closer to the original values from large values, it is possible to approach the particles at a gentle speed that does not lead the molecule to an unnatural state. Further, by gradually reducing the cut-off distance, the attractive force and the repulsive force can be adjusted within an appropriate range.

  The model creation unit 52a performs balancing calculation after setting initial conditions. In the equilibration calculation, molecular dynamics calculation is performed at a predetermined temperature, density, and pressure for a predetermined time for various components after the initial setting to reach an equilibrium state. The model creation unit 52a includes the polymer particles 12a and the bonding chains 12b in the model creation region A for creating a composite material analysis model set in the calculation region after the initial condition setting and equilibration calculation processing. A filler model 11 including a polymer model 12 and filler particles 11a is created. Moreover, the model creation part 52a may mix | blend modifiers, such as a hydroxyl group, a carbonyl group, and a functional group of an atomic group which raise affinity with a filler to a polymer as needed.

  The model creation unit 52a creates a crosslinking point by specifying a predetermined polymer particle 12a in the uncrosslinked polymer model 12 in the composite material analysis model by crosslinking analysis, and connects three or more bonding chains 12b. A polymer model 12 including the polymer particles 12a is created.

  The condition setting unit 52b sets various conditions for executing a motion simulation (analysis) by a molecular dynamics method using the composite material analysis model 1 created by the model creation unit 52a. The condition setting unit 52 b sets various conditions based on the input from the input unit 53 and the information stored in the storage unit 54. As various conditions, the position and number of the filler model 11 for performing the analysis, the position and number of the filler atom, the filler atomic group, the filler particle 11a and the filler particle group, the filler particle 11a number, the position and number of the molecular chain of the polymer, Position and number of polymer atom, polymer atomic group, polymer particle 12a and polymer particle group, polymer particle number, position and number of bond chain 12b and bond chain 12b, bond chain 12b number, stress strain curve which is a preset physical quantity history And fixed values that do not change the conditions.

  The analysis unit 52c uses a molecular dynamics method for relaxation calculation, elongation analysis, deformation analysis such as shear analysis, etc. using the model for analysis 1 of the composite material including the filler model 11 and the polymer model 12 created by the model creation unit 52a. Various physical quantities are acquired by executing the motion simulation. Examples of the physical quantity here include a motion displacement obtained as a result of simulation and a nominal strain obtained by calculating a nominal stress or motion displacement.

  In addition, the analysis unit 52c executes calculation using the molecular dynamics method using the polymer model 12 after cross-linking that is cross-linked by the cross-linking analysis. The analysis unit 52c executes a calculation using a molecular dynamics method using the interaction between the analyzed filler and polymer.

  The analysis unit 52c creates at least one type selected from the group consisting of polymer atoms, polymer particles, and bonding chains of the polymer model 12 after crosslinking as a segment, and creates a segment list specified for each of a plurality of segments to be evaluated. . Here, the analysis unit 52c may specify the entirety of the polymer model 12 before cross-linking as one segment, or may specify between the cross-linking points 12b of the polymer model 12 after cross-linking as segments. A range including the crosslinking point 12b of the later polymer model 12 may be specified as a segment. Here, the analysis unit 52c creates a segment list that is specified for each of a plurality of segments for which at least one of a segment number, a polymer particle number and a bond chain number constituting the segment is an evaluation target. The analysis unit 52c creates a polymer molecule list for each evaluation segment to be evaluated.

  The analysis unit 52 c sets an interaction in the composite material analysis model 1 including the polymer model 12 and the filler model 11. Here, the analysis unit 52c sets interactions between the filler particles 11a, between the polymer particles 12a, between the filler particles 11a and the polymer particles 12a, and the like. Moreover, the analysis part 52c sets between the polymer particle 12a which comprises each polymer model 12, the filler particle 11a, and the polymer particle 12a, when the polymer model 12 contains multiple types of polymer model.

  Furthermore, the analysis unit 52c performs numerical analysis, analyzes and evaluates the molecular state of the polymer model 12 for each identified segment. The analysis unit 52c evaluates the selected segment of the polymer model 12. Here, the analysis unit 52c is separated from the filler model 11 and the segment of the polymer particle 12a in which the region in the vicinity of the filler model 11 of the polymer model 12 after crosslinking is selected as the segment in the polymer molecule list for each evaluation segment, for example. By evaluating the segments of the polymer particles 12a selected from the regions, it is possible to evaluate the molecular motion of the polymer particles 12a of the polymer model 12 in the region near the filler model 11 and the region away from the filler model 11, respectively. Become.

  In addition, the analysis unit 52c calculates a numerical value representing a state change of a specific segment included in the created segment list. Here, the analysis unit 52c calculates numerical values such as the orientation, inertia radius, mean square displacement, end-to-end distance, average bond length, acceleration, and velocity of the segment of the polymer particle 12a included in the segment list.

  Furthermore, the analysis unit 52c adds the distance between the polymer model 12 and the filler model 11 calculated in advance to the segment list, and the specific distance included in the segment list and the distance between the polymer model 12 and the filler model 11 Evaluate the relationship with the numerical value representing the state change of the segment. Here, the analysis unit 52c evaluates the relationship between the distance between the polymer model 12 and the filler model 11 and the numerical value indicating the state change of a specific segment included in the segment list as a table, a graph, or the like.

  The analysis unit 52c visualizes a specific segment included in the segment list based on the numerical analysis result. Here, when visualizing a specific segment, the analysis unit 52c displays only the specific segment or displays the specific segment in a different color from the other segments. Moreover, the analysis part 52c visualizes and displays the specific segment with different materials, such as a filler and another polymer, when visualizing a specific segment. In addition, when the specific segment is visualized, the analysis unit 52c displays a contour according to a numerical value representing a state change of the specific segment included in the calculated segment list. The analysis unit 52 c stores the analysis result of the analyzed composite material in the storage unit 54.

  The storage unit 54 is a non-volatile memory that is a readable recording medium such as a hard disk device, a magneto-optical disk device, a flash memory, and a CD-ROM, and a read / write operation such as a RAM (Random Access Memory). A volatile memory which is a possible recording medium is appropriately combined.

  In the storage unit 54, data on fillers such as rubber carbon black, silica, and alumina, which are data for creating a model for analysis of a composite material to be analyzed via the input means 53, rubber, resin, and elastomer Of data such as polymer data, stress-strain curve, which is a preset physical quantity history, and a method for creating a composite material analysis model, a composite material simulation method, and a composite material analysis result evaluation method A computer program or the like is stored. This computer program may be capable of realizing the composite material simulation method according to the present embodiment in combination with a computer program already recorded in a computer or computer system. The “computer system” here includes hardware such as an OS (Operating System) and peripheral devices.

  The display means 55 is a display device such as a liquid crystal display device. The storage unit 54 may be in another device such as a database server. For example, the analysis device 50 may access the processing unit 52 and the storage unit 54 by communication from a terminal device including the input / output device 51.

  As described above, according to the evaluation method of the analysis result of the composite material according to the above embodiment, it is possible to perform numerical analysis according to the interaction set for each segment of the polymer model 12 after crosslinking. It becomes possible to analyze and evaluate the molecular state of the polymer model 12 for each segment of the polymer model 12. Thereby, the evaluation method of the analysis result of the composite material is to analyze the state change of the molecular state of the polymer model 12 in the region near the filler model 11 and the molecular state of the polymer model 12 in the region far from the filler model 11. Therefore, it is possible to realize a method for evaluating the analysis result of a composite material capable of evaluating a local molecular state change of a polymer molecule.

DESCRIPTION OF SYMBOLS 1 Model for analysis 11 Filler model 11a Filler particle 12,12-1,12-2 Polymer model 12a, 12a-1,12a-2 Polymer particle 12b, 12b-1,12b-2 Crosslinking point 121,122 Segment 50 Analysis apparatus 51 I / O Device 52 Processing Unit 52a Model Creation Unit 52b Condition Setting Unit 52c Analysis Unit 53 Input Unit 54 Storage Unit 55 Display Unit A Area L1 Distance S1 to S12 Segment

Claims (11)

A composite material analysis method using a composite material analysis model created by a molecular dynamics method using a computer,
A first step of creating a model for analysis of a composite material including a polymer model and a filler model;
A second step of crosslinking the polymer model by crosslinking analysis;
A third step of creating a segment list in which at least one selected from the group consisting of polymer atoms, polymer particles, and bonding chains included in the polymer model after crosslinking is a segment, and a plurality of segments to be analyzed are specified; ,
A fourth step of analyzing the molecular state of the polymer model for each segment specified by numerical analysis by setting an interaction in the analysis model of the composite material;
A method for analyzing a composite material, comprising:   A computer program for analyzing a composite material, which causes a computer to execute the composite material analysis method according to claim 1. A method for evaluating the analysis result of a composite material analyzed using a model for analysis of a composite material created by a molecular dynamics method using a computer,
A first step of creating a model for analysis of a composite material including a polymer model and a filler model;
A second step of crosslinking the polymer model by crosslinking analysis;
A third step of creating a segment list in which at least one selected from the group consisting of polymer atoms, polymer particles, and bonding chains included in the polymer model after crosslinking is a segment, and a plurality of segments to be evaluated are specified; ,
A fourth step of setting an interaction in an analysis model of a composite material including a polymer model and a filler model, and analyzing and evaluating a molecular state of the polymer model for each segment specified by numerical analysis;
A method for evaluating an analysis result of a composite material, comprising:   The evaluation method of the analysis result of the composite material according to claim 3, wherein in the fourth step, a numerical value representing a state change of a specific segment included in the segment list is calculated.   In the third step, a distance between the polymer model and the filler model is added to the segment list, and in the fourth step, a distance between the polymer model and the filler model is added to the segment list. The evaluation method of the analysis result of the composite material of Claim 4 which evaluates the relationship with the numerical value showing the state change of the specific segment contained.   5. The composite according to claim 4, wherein in the fourth step, a numerical value representing a state change of a specific segment included in the segment list is calculated in a plurality of analysis times, and a relationship between the calculated numerical value and the analysis time is evaluated. Evaluation method for analysis results of materials.   The method for evaluating an analysis result of a composite material according to claim 6, wherein in the fourth step, deformation analysis is executed to evaluate a relationship between a numerical value representing a state change of the segment and strain or analysis time.   The method for evaluating an analysis result of a composite material according to claim 6, wherein in the fourth step, a transformation analysis is executed to evaluate a relationship between a numerical value representing a state change of the segment and a pressure or an analysis time.   The evaluation method of the analysis result of the composite material according to claim 6, wherein in the fourth step, a temperature change analysis is executed to evaluate a relationship between a numerical value representing a state change of the segment and a temperature or an analysis time.   The evaluation method of the analysis result of the composite material according to any one of claims 3 to 9, wherein in the fourth step, a specific segment included in the segment list is visualized based on a numerical analysis result.   A computer program for evaluating the analysis result of a composite material, which causes a computer to execute the method for evaluating the analysis result of a composite material according to any one of claims 3 to 10.

Images