Polymer model generation method, system, and program

Abstract

To provide a method, a system and a program for polymer model generation, being usable for highly accurate evaluation even in a case of a ratio of monomers having no experimental data.SOLUTION: A system 1 for polymer model generation, includes an angle potential calculation section 10 of determining angle potentials of first to third type monomer models 5-7 so as to match experimental sites. A monomer allocation section 11 allocates any of the first to third type monomer models 5-7 to a plurality of coarse-grained particles of the polymer model. A potential setting section 12 sets a potential of the allocated model to each coarse-grained particle. A cross-linking processing execution section 13 so sets that each coarse-grained particle can bind or cannot bind to a cross-linking particle model according to the allocated monomers and executes cross-linking processing of binding the bindably-set first type monomer model 5 to the cross-linking particle model.SELECTED DRAWING: Figure 1

Description

[Technical field]

[0001] The present disclosure relates to a method, system, and program for generating a polymer model. [Background technology]

[0002] Hydrogenated SBR (styrene butadiene rubber) is known to have better abrasion resistance than non-hydrogenated SBR. In order to understand the reasons for this superiority, it is necessary to create a molecular dynamics simulation model that reproduces hydrogenated SBR. [Prior art documents] [Non-patent literature]

[0003] [Non-Patent Document 1] Simulation of SS curves using entanglement density control model, Hiroki Yamamura, Tetsuo Tominaga, Takumi Adachi, Hirofumi Senga, JSR TECHNICAL REVIEW No.126 / 2019 Summary of the Invention [Problem to be solved by the invention]

[0004] Non-Patent Document 1 discloses a method of using the Kremer-Grest model to model styrene and butadiene together without distinguishing between them. Specifically, an angle potential is added to the Kremer-Grest model, which is a styrene-butadiene model, and the angle potential is determined so as to match the experimental viscoelasticity (equilibrium elastic modulus, Young's modulus). After that, the number of crosslinks in the styrene-butadiene model is determined so that the hydrogenation rate matches the experimental value of the SS curve (stress-strain curve) of the known hydrogenated SBR.

[0005] However, this method uses experimental data with known hydrogenation rates and cannot generate models with unknown hydrogenation rates. In addition, a method is also required to generate polymer models that are not limited to hydrogenated SBR, but in which the number of crosslinks with crosslinked particles can be changed depending on the proportion of monomers.

[0006] The present disclosure provides a method, system, and program for generating a polymer model that can be used for accurate evaluation even for monomer ratios for which no experimental data is available. [Means for solving the problem]

[0007] The method for generating a polymer model according to the present disclosure is a method for generating a polymer model executed by one or more processors, the polymer model including a plurality of polymer models having a plurality of coarse-grained particles connected in a linear or branched manner, and a bridge particle model, the plurality of coarse-grained particles representing any of a plurality of types of monomer models including a first type monomer model and a second type monomer model, the first type monomer model and the second type monomer model each having a plurality of particles connected in a linear or branched manner and each having an angular potential representing an interaction between three particles, the first type monomer model can be set to be bondable with the bridge particle model, and the second type monomer model can be set to be bondable with the bridge particle model, and the bridge particle model can be set to be bonded with the bridge particle model. the first type monomer model and the second type monomer model are set to be unable to be bonded to a crosslinked particle model, and the method includes: calculating the angular potential of each of the first type monomer model and the second type monomer model based on an experimental value; assigning one of the plurality of types of monomer models to each of the plurality of coarse-grained particles constituting the polymer model such that a ratio of the second type monomer model to the first type monomer model is a specified ratio; setting a potential corresponding to the assigned monomer model for each of the plurality of coarse-grained particles; and performing a crosslinking process to bond the first type monomer model and the crosslinked particle model in the polymer model. [Brief description of the drawings]

[0008] [Figure 1] FIG. 1 is a block diagram showing a system according to an embodiment of the present invention. [Diagram 2] 4 is a flowchart showing a process executed by the system. [Diagram 3] FIG. 2 is a schematic diagram showing a polymer model (SBR model) generated in this embodiment. [Figure 4] 4 is an explanatory diagram of an angular potential (angle parameter) and a characteristic ratio. [Diagram 5] 11A and 11B are diagrams showing the relationship between the characteristic ratio and the angle parameter, and between the angular potential and the angle parameter. [Figure 6] FIG. 11 is an explanatory diagram relating to a process of assigning a monomer model to each particle of a polymer model. [Figure 7] FIG. 2 is an explanatory diagram of four polymer models created in this embodiment. [Figure 8] FIG. 2 is a diagram showing SS curves of four polymer models created in this embodiment. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0009] Hereinafter, an embodiment of the present disclosure will be described with reference to the drawings.

[0010] [Model generation system] The system 1 (apparatus) of this embodiment generates a polymer model. In this embodiment, the polymer model is hydrogenated SBR. SBR is a copolymer of butadiene and styrene, and the butadiene portion is crosslinked with a crosslinking agent such as sulfur. When SBR is hydrogenated, hydrogen is added to the double bond of the butadiene portion, resulting in an ethylene structure. This structural change is modeled.

[0011] As shown in FIG. 1, the system 1 includes an angular potential calculation unit 10, a monomer allocation unit 11, a potential setting unit 12, and a crosslinking process execution unit 13. The system 1 may include an extension process execution unit 14 that stretches the generated polymer model. Each of these units 10 to 14 is realized by software and hardware working together as the processor 1a executes a processing routine shown in FIG. 2 that is stored in advance in a computer equipped with a processor 1a, a memory 1b, various interfaces, etc. In this embodiment, the processor 1a in one device realizes each unit, but this is not limited to this. For example, the system may be configured to be distributed using a network, with multiple processors executing the processing of each unit. That is, one or multiple processors execute the processing. The memory 1b stores the generated polymer model 2, multiple types of monomer models (5, 6, 7), and analysis conditions (predetermined temperature, predetermined pressure, etc.) for executing molecular dynamics calculations.

[0012] As shown in FIG. 3, the polymer model 2 has a plurality of polymer models 3 and a crosslinked particle model 4. The polymer model 3 has a plurality of coarse-grained particles 30 that are linearly or branched. The crosslinked particle model 4 bonds (crosslinks) with the particles 30 of the polymer model 3 that are set to be crosslinkable. The coarse-grained particles 30 represent any of a plurality of types of monomer models including a first type monomer model 5 and a second type monomer model 6. In FIG. 3, each of the coarse-grained particles that constitute the chain of the polymer 3 is depicted as if it further contains a monomer having a plurality of particles, but this is depicted for ease of understanding. The coarse-grained particles 30 do not have particles that represent monomers therein, but various potentials (described later) for representing the characteristics of the monomers are set, and the number of particles is reduced to reduce calculation costs. In this embodiment, the multiple types of monomer models further include a third type monomer model 7, the first type monomer model 5 represents butadiene, the second type monomer model 6 represents ethylene, and the third type monomer model 7 represents styrene. The first type monomer model 5 represents butadiene, and can be set to be bondable with the crosslinked particle model 4. This does not mean that all of the first type monomer models 5 need to be set to be bondable with the crosslinked particle model 4, but a part of the first type monomer model 5 may be set to be bondable with the crosslinked particle model 4, and the remaining part may be set to be unable to be bonded with the crosslinked particle model 4. The second type monomer model 6 represents ethylene, and is set to be unable to be bonded with the crosslinked particle model 4. In addition, the third type monomer model 7 represents styrene, and is set to be unable to be bonded with the crosslinked particle model 4.

[0013] In each of the multiple types of monomer models, multiple particles are connected in a linear or branched manner. In each of the multiple types of monomer models, a bonding potential and a non-bonding potential are set for each particle. In the first type monomer model 5 in this embodiment, multiple particles 50 are connected in a linear or branched manner by a bonding potential, and a non-bonding potential is set for each particle 50. In the second type monomer model 6, multiple particles 60 are connected in a linear or branched manner by a bonding potential, and a non-bonding potential is set for each particle 60. In the third type monomer model 7, multiple particles 70 are connected in a linear or branched manner by a bonding potential, and a non-bonding potential is set for each particle 70. The bonding potential is a potential acting between a bonded partner. The non-bonding potential is a potential acting between other particles.

[0014] The multiple types of monomer models 5 to 7 of this embodiment use the Kremer-Grest model, and the non-bonded potential (U non-bond ) and the binding potential of the Kremer-Grest model (U bond) and the Kremer-Grest model has a unique characteristic ratio that represents the spread and entanglement of the polymer chain, so the potential of multiple types of monomer models is further increased by an angular potential (U angle ) has been added. The characteristic ratio can be changed by the angular potential, which changes depending on the angle with the bonded molecule. The non-bonded potential can be expressed by equation (1), and the bonded potential can be expressed by equation (2). The angular potential can be expressed by equation (3). This equation is just one example.

number

[0015] The characteristic ratio is the average of all particle chains (square of the distance between the ends) / (number of particles × square of the distance between particles) = <(r1-r N ) 2 〉 / (NL 2 ) where N represents the number of particles constituting one particle chain. As shown in FIG. 4 by way of example in the first type monomer model 5, the particle chain is connected (bonded) in a linear or branched manner. The particle chain is divided into a first end r1 and a second end r2, and the particle chain is connected (bonded) in a linear or branched manner. N The distance between is |r1-r N If the particle's identification number is i, then the position is expressed as r i The position of the particle with identification number 1 is r1, the position of the particle with identification number 2 is r2, and the position of the particle with identification number N is r N The distance between two adjacent particles 50 in a bonded state is represented by L. The characteristic ratio is a value obtained from light scattering in an experiment. When the particle chain is branched, the one with the longest distance between the above ends is used.

[0016] The angular potential calculation unit 10 calculates the angular potential U angle (Angle parameter (K angle)) are calculated based on experimental values. In this embodiment, the various parameters for determining the potential are unified to K=30.0, R0=1.5, ε=1.0, and σ=1.0 for all models. In this embodiment, the first type monomer model 5 is butadiene, and the characteristic ratio described in the literature (Polymer, Vol. 48, September issue, p. 737 (1999), Society of Polymer Science) is 5.4. The second type monomer model 6 is ethylene, and the characteristic ratio described in the above literature is 6.8. The third type monomer model 7 is styrene, and the characteristic ratio described in the above literature is 10.8. Each monomer model in this embodiment is a model in which 200 particles are connected in a linear chain using the Kremer-Grest model, and the inter-particle distance L is 0.97. An angle parameter (K angle ), the relationship between the angular potential U(θ) and the angle θ can be calculated as shown in the lower diagram of Figure 5. angle = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10) to the first type monomer model 5. If the molecular dynamics calculation is continued until the equilibrium is reached, the end-to-end distance |r1-r N | can be calculated, and the angle parameter (K angle ) and the characteristic ratio. angle ) = 2 to 7 data is approximated by the formula [characteristic ratio = α × K angle +β] was fitted using the least squares method to calculate the coefficients α and β of the approximation equation. α = 1.7288, β = -0.0513, and the sum of squares of the residuals (coefficient of determination) R 2 =0.9987. The angular potential calculation unit 10 calculates the angular parameter (K angle ) and determine the angular potential. Using the above approximation, the angular parameter (K angle ) is 3.1, and the angular parameter of the second kind monomer model (K angle ) is 4.0, and the angular parameter of the third kind monomer model (K angle ) calculates to 6.3. From the above, the angular potential (parameters that determine the angular potential) of each monomer model 5 can be calculated so as to match the characteristic ratio of each monomer.

[0017] As shown in FIG. 6, the monomer allocation unit 11 allocates one of the multiple types of monomer models 5 to 7 to each of the multiple coarse-grained particles 30 constituting the polymer model 3 (hydrogenated SBR). In FIG. 6, the particle 30 to which the first type monomer model 5 (butadiene) is allocated is represented as B, the particle 30 to which the second type monomer model 6 (ethylene) is allocated is represented as E, and the particle 30 to which the third type monomer model (styrene) is allocated is represented as S. At this time, the ratio of the second type monomer model 6 (ethylene) to the first type monomer model 5 (butadiene) is set to a specified ratio. The specified ratio is a ratio designated by the user and can be freely set. For example, SBR is 23.5% styrene and the remaining 76.5% butadiene. When hydrogen is added, hydrogen is added to the double bond of butadiene to form an ethylene structure, so the ratio of butadiene decreases and the ratio of ethylene increases instead. If the hydrogenation rate is Y [%], the rate [%] of butadiene is 76.5 × (100-Y), and the rate [%] of ethylene is 76.5 × Y. The monomer allocation unit 11 determines the number of coarse-grained particles 30 to be allocated to the first type monomer model 5 and the number of coarse-grained particles 30 to be allocated to the second type monomer model 6 so as to achieve the specified rate. In this embodiment, one polymer model 3 has 200 coarse-grained particles, and three types of monomer models 5 to 7 are randomly allocated so as to match the allocation number of each monomer model. The number of particles allocated to the third type monomer model 7 (styrene) is 47. If the hydrogenation rate is 0%, the number of particles allocated to the first type monomer model 5 (butadiene) is 153. If the hydrogenation rate is 10%, the number of particles allocated to the first type monomer model 5 (butadiene) is about 137, and the number of particles allocated to the second type monomer model 6 (ethylene) is about 16.

[0018] The potential setting unit 12 sets a potential corresponding to the assigned monomer models 5 to 7 for each of the multiple coarse-grained particles 30. In this embodiment, the potentials of the formulas (1) and (2) are the same regardless of the type of monomer model. The potential setting unit 12 sets an angle parameter (K angle ) was set to 3.1, and the angular parameter (K angle ) was set to 4.0, and the angular parameter (K angle ) is set to 6.3 and the angular potential is set.

[0019] The crosslinking process execution unit 13 places a designated number of polymer models 3 and crosslinked particle models 4 in space, performs an equilibration process, and then executes a crosslinking process to bond the first type monomer model 5 and the crosslinked particle model 4 in the polymer model 3 in an equilibrium state. The equilibration process executes molecular dynamics calculations at a predetermined pressure and a predetermined temperature until an equilibrium state is reached. When the first type monomer model 5 and the crosslinked particle model 4 approach each other within a predetermined distance, the crosslinking process execution unit 13 bonds them with a predetermined probability. The crosslinking process execution unit 13 repeatedly executes the process until all the crosslinked particle models 4 have been bonded. This completes the generation of the hydrogenated SBR model (polymer model 2). An example of an equilibrium state is a state in which the fluctuation in the potential of the entire model over a certain period of time falls within a certain range.

[0020] Although all the coarse-grained particles 30 of the first type monomer model 5 (butadiene) may be set to be bondable to the crosslinked particle model 4, since the number of the first type monomer models 5 (butadiene) in the polymer model 3 is relatively large compared to other monomer models, it is preferable to set only a part of the coarse-grained particles 30 of the first type monomer model 5 (butadiene) to be bondable to the crosslinked particle model 4. Since it is preferable that the crosslinked particle model 4 is bonded to the first type monomer model 5 of a certain polymer model 3 and to the first type monomer model 5 of another polymer model 3, it is preferable that a predetermined number of coarse-grained particles centered on the coarse-grained particle 30 set to be bondable to the crosslinked particle model 4 are set to be unable to be bonded to the crosslinked particle model 4. This makes it possible to prevent one crosslinked particle model 4 from being bonded to multiple first type monomer models 5 (butadiene) of the same polymer model 3, thereby improving the uniformity of the crosslinking. In addition, since the number of first type monomer models 5 (butadiene) decreases as the hydrogenation rate increases, it is preferable to reduce the number of coarse-grained particles 30 (first type monomer models 5) that can be bonded to the crosslinked particle models 4 in accordance with the reduction in the number of first type monomer models 5 (butadiene). This makes it possible to improve the uniformity of crosslinking.

[0021] [How to generate polymer models] A method for generating a polymer model executed by the above system 1 will be described with reference to Fig. 2. Here, an example is given of generating a hydrogenated SBR model as a polymer model. A predetermined pressure and a predetermined temperature used in the molecular dynamics calculation are set and stored in advance in the memory 1b. Also, the characteristic ratios of styrene, butadiene, and ethylene are set and stored in advance in the memory 1b.

[0022] First, in step ST1, the angular potential calculation unit 10 calculates the angular potentials (angular parameters) of the first to third monomer models (butadiene, ethylene, and styrene) based on experimental values. In this embodiment, the angular parameter (K angle ) was calculated to be 3.1, and the angle parameter (K angle) was calculated to be 4.0, and the angle parameter (K angle ) is calculated to be 6.3.

[0023] 6, in the next step ST2, the monomer allocation unit 11 allocates any one of the first to third types of monomer models to each of the multiple coarse-grained particles 30 constituting the polymer model 3. The ratio of each is set to a designated ratio.

[0024] In step ST2, when the assignment of the first to third types of monomer models to all the coarse-grained particles 30 is completed, in the next step ST3, the potential setting unit 12 sets a potential corresponding to the assigned monomer model for each of the multiple coarse-grained particles 30. In this embodiment, since the bond potentials other than the angular potential are common, only the angular potential is set to the potential corresponding to the assigned monomer model.

[0025] In the next step ST4, the system 1 sets the coarse-grained particles 30 assigned to the second type monomer model (ethylene) and the third type monomer model 7 (styrene) as unable to bond with the crosslinked particle model 4. In addition, some of the coarse-grained particles 30 of the first type monomer model (butadiene) are set as able to bond with the crosslinked particle model 4. The remaining coarse-grained particles 30 of the first type monomer model (butadiene) are set as unable to bond with the crosslinked particle model 4. In this embodiment, about 13% of the number of the coarse-grained particles 30 assigned to the first type monomer model (butadiene) are set as able to bond with the crosslinked particle model 4, but the proportion is not limited to this.

[0026] In the next step ST5, the crosslinking process execution unit 13 executes a crosslinking process to bond the first type monomer model (butadiene) set to be bondable in the polymer model 3 to the crosslinked particle model 4. This completes the generation of the model. Before the crosslinking process, it is recommended to execute an equilibration process to stabilize the uncrosslinked polymer model. Once the crosslinking process is completed (when all the crosslinked particle models have been bonded), it is recommended to execute an equilibration process.

[0027] [Tensile test simulation] A tensile test simulation was performed to confirm whether the created polymer model exhibits qualitative behavior. A known method described in JP 2017-96871 A is used, so detailed description is omitted. The extension processing execution unit 14 shown in FIG. 1 places a polymer model in a virtual space, gradually extends the virtual space in the z-axis direction under constant volume and temperature conditions, and calculates the extension ratio and stress. Since it is not possible to extend the polymer model to the desired extension ratio all at once, a process of slightly extending the polymer model and a relaxation process (equilibrium process) are repeatedly performed. For example, the polymer model is extended to about 1.0139 times, then a relaxation process is performed, and a process of extending the polymer model to 1.0139 times again is repeatedly performed. If the process of extending the polymer model to about 1.0139 times in the z-axis direction is repeated 100 times, the extension ratio can be obtained from 1 to 7, and the stress at each time point can be obtained, and the strain can also be obtained.

[0028] Working Example As shown in FIG. 7, four polymer models (hydrogenated SBR) with hydrogenation rates of 0%, 30%, 50%, and 70% were generated. In each model, the chain length of one polymer model 3 is 200 (200 coarse-grained particles 30), the number of polymer models 3 is 200, and 200 crosslinked particle models 4 are arranged as single particles. As shown in FIG. 7, with an increase in the hydrogenation rate, the number of first-type monomer models 5 (butadiene) that are set to be bondable with the crosslinked particle model 4 (number of bondable points) is reduced. This is to improve the uniformity of the crosslinking. FIG. 7 shows the four polymer models and their total particle number, butadiene particle number, the ratio set to the bondable points, the number of bondable points, and the ratio of the number of bondable points to the total particles. Figure 8 shows the SS curves for the four hydrogenated SBR models (polymer models) shown in Figure 7. These were obtained by the tensile simulation described above. In all of the models, the stress increases as the hydrogenation rate increases, and it can be seen that the hydrogenated SBR can be qualitatively reproduced. The stress in Figure 8 increases upward. Each grid on the vertical axis shows the same amount.

[0029] As described above, the method for generating a polymer model according to the present embodiment is a method for generating a polymer model executed by one or more processors, in which the polymer model 2 includes a plurality of polymer models 3 having a plurality of coarse-grained particles 30 connected in a linear or branched manner, and a crosslinked particle model 4, the plurality of coarse-grained particles 30 represent any one of a plurality of types of monomer models 5 to 7 including a first type monomer model 5 and a second type monomer model 6, the first type monomer model 5 and the second type monomer model 6 each have a plurality of particles (50, 60) connected in a linear or branched manner and have an angular potential representing an interaction between the three particles, the first type monomer model 5 can be set to be bondable with the crosslinked particle model 4, and the second type monomer model 6 can be set to be bondable with the crosslinked particle model 4. The polymer model 6 is set so as not to be able to bond with the crosslinked particle model 4, and the method may include: calculating the angular potentials of the first type monomer model 5 and the second type monomer model 6 based on experimental values; assigning one of multiple types of monomer models 5 to 7 to each of the multiple coarse-grained particles 30 constituting the polymer model 3 so that the ratio of the second type monomer model 6 to the first type monomer model 5 is a specified ratio; setting a potential corresponding to the assigned monomer model for each of the multiple coarse-grained particles 30; and performing a crosslinking process to bond the first type monomer model 5 and the crosslinked particle model 4 in the polymer model 3.

[0030] In this way, instead of determining the angle parameters of a model in which the first type monomer (butadiene) and the second type monomer (ethylene) are mixed without distinction, the first type monomer model 5 and the second type monomer model 6 are distinguished and the angular potentials of each monomer model 5 to 7 are determined. Each monomer model 5 to 7 is assigned to a coarse-grained particle that constitutes the polymer model. Therefore, the ratio of the first type monomer and the second type monomer can be freely controlled, and the degree of entanglement can be controlled by controlling the ratio. Since the characteristic ratio can be reproduced by the degree of entanglement, a model that can be used for accurate evaluation can be generated even if the ratio is unknown.

[0031] As in this embodiment, the multiple types of monomer models 5 to 7 may include a third type monomer model 7 that cannot be bonded to the crosslinked particle model 4. This makes it possible to increase the types of polymers that can be reproduced.

[0032] As in this embodiment, the first monomer model 5 may be a butadiene model, the second monomer model 6 may be an ethylene model, the third monomer model 7 may be a styrene model, and the polymer model 2 may be a hydrogenated SBR model. This makes it possible to generate a model of hydrogenated SBR, where the ratio is the hydrogenation rate.

[0033] As in this embodiment, the method includes generating a first polymer model (hydrogenation rate 0%) and a second polymer model (hydrogenation rate 30%) as the polymer model 2, and the ratio of the second type monomer model 6 to the first type monomer model 5 in the second polymer model (hydrogenation rate 30%) is higher than the ratio of the second type monomer model 6 to the first type monomer model 5 in the first polymer model (hydrogenation rate 0%), and the number of first type monomer models set to be crosslinkable in the second polymer model (hydrogenation rate 30%) is set to be less than the number of first type monomer models set to be crosslinkable in the first polymer model. According to this configuration, when the ratio of the second type monomer model to the first type monomer model becomes high, the number of first type monomer models that can be set as crosslinkable decreases, and by reducing the number of first type monomer models that are set as crosslinkable accordingly, unintended crosslinking can be prevented and crosslinking uniformity can be improved.

[0034] The program according to this embodiment is a program for causing one or more computers to execute the above method. By executing these programs, it is possible to obtain the effects of the above-mentioned methods.

[0035] Although the embodiments of the present disclosure have been described above based on the drawings, the specific configurations should not be considered to be limited to these embodiments. The scope of the present disclosure is indicated not only by the description of the above-mentioned embodiments but also by the claims, and further includes all modifications within the meaning and scope equivalent to the claims.

[0036] (1) Although the polymer model 2 in the above embodiment has the third monomer model 7 as styrene, the third monomer model 7 can be omitted.

[0037] (2) In the above embodiment, the polymer model is SBR, but is not limited thereto. As long as it has two types of monomers whose ratios are in a competitive relationship (negative correlation) with each other, it can be used to reproduce various polymers.

[0038] The structures employed in the above-described embodiments may be employed in any other embodiment. The specific configurations of the components are not limited to the above-described embodiments, and may be modified in various ways without departing from the spirit of the present disclosure.

[0039] For example, the order of execution of each process, such as operations, procedures, steps, and stages, in the devices, systems, programs, and methods shown in the claims, specifications, and drawings, can be arbitrary, as long as the output of a previous process is not used in a subsequent process. Even if the flow in the claims, specifications, and drawings is described using "first," "next," and the like for convenience, this does not mean that it is essential to execute the processes in this order.

[0040] Each unit shown in Fig. 1 is realized by executing a predetermined program on one or more processors, but each unit may be configured with a dedicated memory or a dedicated circuit. In the system 1 of the above embodiment, each unit is implemented in the processor 1a of one computer, but each unit may be distributed and implemented on multiple computers or in the cloud. In other words, the above method may be executed on one or more processors.

[0041] The system 1 includes a processor 1a. For example, the processor 1a can be a central processing unit (CPU), a microprocessor, or other processing unit capable of executing computer-executable instructions. The system 1 also includes a memory 1b for storing data for the system 1. In one example, the memory 1b includes a computer storage medium, including RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, DVD or other optical disk storage, magnetic cassette, magnetic tape, magnetic disk storage or other magnetic storage device, or any other medium that can be used to store desired data and that can be accessed by the system 1. [Explanation of symbols]

[0042] 3...polymer model, 4...crosslinked particle model, 5...first type monomer model, 6...second type monomer model, 7...third type monomer model, 30...coarse-grained particle, 50, 60, 70...particles.

Claims

1. 1. A method for generating a polymer model, executed by one or more processors, comprising: The polymer model includes a plurality of polymer models having a plurality of coarse-grained particles connected in a linear or branched manner, and a bridged particle model; the plurality of coarse-grained particles represent any one of a plurality of types of monomer models including a first type monomer model and a second type monomer model; The first type monomer model and the second type monomer model each have a plurality of particles connected in a linear or branched chain form and an angular potential representing an interaction between three particles, The first type monomer model can be set to be capable of being bound to the bridge particle model, The second type monomer model is set to be incapable of bonding with the bridge particle model, The method comprises: Calculating the angular potentials of the first type monomer model and the second type monomer model based on experimental values; Allocating any one of the plurality of types of monomer models to each of the plurality of coarse-grained particles constituting the polymer model such that a ratio of the second type monomer model to the first type monomer model is a designated ratio; Setting a potential corresponding to the assigned monomer model for each of the plurality of coarse-grained particles; and performing a cross-linking process to bond the first type monomer model and the cross-linked particle model in the polymer model.

2. The method of claim 1 , wherein the plurality of types of monomer models includes a third type of monomer model that cannot be bonded to the bridge particle model.

3. 3. The method of claim 2, wherein the first monomer model is a butadiene model, the second monomer model is an ethylene model, the third monomer model is a styrene model, and the polymer model is a hydrogenated SBR model.

4. The method includes generating a first polymer model and a second polymer model as the polymer model; a ratio of the second type monomer model to the first type monomer model in the second polymer model is higher than a ratio of the second type monomer model to the first type monomer model in the first polymer model; The method according to any one of claims 1 to 3, wherein the number of the first type monomer models set to be crosslinkable in the second polymer model is set to be less than the number of the first type monomer models set to be crosslinkable in the first polymer model.

5. A system comprising one or more processors for carrying out the method according to any one of claims 1 to 4.

6. A program for causing one or more processors to execute the method according to any one of claims 1 to 4.

Images