Patent Application
20150302180
The present invention relates to a simulation method for a high-polymer material which is useful for estimating and improving the dispersion of filler.
In a high-polymer material (e.g. rubber) containing filler (e.g. carbon black, silica and the like), the strength is largely affected by the dispersion of the filler.
In recent years, in order to improve and develop compositions of rubber compounds, there have been proposed various computer-implemented numerical simulation methods for estimating the dispersion of filler contained in a high-polymer material.
In this kind of numerical simulation method, a numerical model of the filler is defined by a plurality of particles (hereinafter “F-particle”), and a numerical model of the high-polymer material is defined by a plurality of particles (hereinafter “P-particle”). Then, the P-particles are bound with the F-particles, and a molecular dynamics calculation is performed by a computer.
In the above-described simulation method, P-particles to be bound with each F-particle are randomly selected by a computer. Consequently, the number of the P-particles bound with one F-particle varies largely. If the number becomes large, as the total force of the P-particles exerted on the F-particle becomes large, the molecular dynamics calculation has a high probability of causing a calculation failure.
It is therefore, a primary object of the present invention to provide a simulation method for a high-polymer material in which, a calculation failure during a molecular dynamics calculation is efficiency prevented by limiting the number of P-particles to be bound with each F-particle not to exceed a predetermined upper limit value.
According to one aspect of the present invention, a simulation method for a high-polymer material is a computer-implemented method for estimating the dispersion of filler contained in the high-polymer material, and comprises:
a process for defining filler models of the filler, wherein each filler model is a numerical model made up of a plurality of F-particles,
a process for defining polymer models of the high-polymer material, wherein each polymer models is a numerical model made up of a plurality of P-particles, and
a simulation process for executing a molecular dynamics calculation by the use of the filler models and the polymer models disposed in a predetermined virtual space,
the simulation process including a bonding process for binding bonding-target P-particles selected from the P-particles with bonding-target F-particles selected from the F-particles, wherein the number of the bonding-target P-particle(s) to be bound with each bonding-target F-particle satisfies a predetermined upper limit value Su.
Here, the term “F-particle” means a particle of a filler model. The term “P-particle” means a particle of a polymer model.
The upper limit value Su may satisfy the following expression:
Na/Nb−1=<Su=<Na/Nb+1
wherein
Na is a total number of the bonding-target P-particles existing in the virtual space, and
Nb is a total number of the bonding-target F-particles existing in the virtual space.
The bonding process may includes:
a process for defining the upper limit value Su,
a selecting process for selecting the bonding-target P-particle(s) to be bound with each bonding-target F-particle, and
a process for binding the selected bonding-target P-particle(s) with the concerned bonding-target F-particle.
The selecting process may includes:
a 1st process for associating each bonding-target P-particle with a bonding-target F-particle to which interparticle distance is smallest,
a 2nd process for selecting the bonding-target P-particle(s) to be bound with each bonding-target F-particle from a group of the bonding-target P-particle(s) associated with the concerned bonding-target F-particle in the ascending order of the interparticle distance to satisfy the upper limit value Su,
a third process for associating each bonding-target P-particle not selected in the 2nd process with a bonding-target F-particle to which interparticle distance is second smallest next to said bonding-target F-particle to which interparticle distance is smallest, and
a fourth process for selecting the bonding-target P-particle(s) to be bound with each bonding-target F-particle from a group of the bonding-target P-particle(s) associated in the third process, in the ascending order of the interparticle distance to satisfy the upper limit value su, and
the third process and the fourth process are repeated until all of the bonding-target P-particles are selected for the bonding-target F-particles.
Preferably, the F-particles of each filler model include those disposed at apexes of a polyhedron, and the bonding-target F-particles are the F-particles disposed at the apexes.
According to another aspect of the present invention, a simulation method for a high-polymer material is a computer-implemented method for estimating the dispersion of filler contained in the high-polymer material together with a coupling agent for coupling the filler with polymer by the use of a computer, comprises:
a process for defining filler models of the filler, wherein each filler model is a numerical model made up of a plurality of F-particles,
a process for defining polymer models of the high-polymer material, wherein each polymer models is a numerical model made up of a plurality of P-particles, and
a process for defining coupling agent models of the coupling agent, wherein each coupling agent model is a numerical model made up of at least one C-particle,
a simulation process for executing a molecular dynamics calculation by the use of the filler models, the coupling agent models and the polymer models disposed in a predetermined virtual space,
the simulation process including a bonding process in which bonding-target P-particles selected from the P-particles and bonding-target C-particles selected from the C-particles are bound with bonding-target F-particles selected from the F-particles, wherein
the number of the bonding-target P-particle(s) and/or the bonding-target C-particle(s) to be bound with each bonding-target F-particle satisfies a predetermined upper limit value Su. Here, The term “F-particle” means a particle of a filler model. The term “P-particle” means a particle of a polymer model.
The term “C-particle” means a particle of a coupling agent model.
The upper limit value Su may satisfy the following expression:
Ns=[(Na+Nc)/Nb+1]
wherein
Na is a total number of the bonding-target P-particles existing in the virtual space,
Nb is a total number of the bonding-target F-particles existing in the virtual space, and
Nc is a total number of the bonding-target C-particles existing in the virtual space.
The bonding process may includes:
a process for defining the upper limit value Ns,
a selecting process for selecting the bonding-target P-particle(s) and/or the bonding-target C-particle(s) to be bound with each bonding-target F-particle, and
a process for binding the selected bonding-target P-particle(s) and/or bonding-target C-particle(s) with the concerned bonding-target F-particle.
The selecting process may includes:
a 1st process for associating each bonding-target P-particle with a bonding-target F-particle to which interparticle distance is smallest, and associating each bonding-target C-particle with a bonding-target F-particle to which interparticle distance is smallest,
a 2nd process for selecting the bonding-target P-particle(s) and/or the bonding-target C-particle(s) to be bound with each bonding-target F-particle from the associated bonding-target P-particle(s) and/or bonding-target C-particle(s) in the ascending order of the interparticle distance to satisfy the upper limit value Ns,
a third process for associating each bonding-target P-particle not selected in the 2nd process with a bonding-target F-particle to which interparticle distance is second smallest next to said bonding-target F-particle to which interparticle distance is smallest, and
associating each bonding-target C-particle not selected in the 2nd process with a bonding-target F-particle to which interparticle distance is second smallest next to said bonding-target F-particle to which interparticle distance is smallest, and
a fourth process for selecting the bonding-target P-particle(s) and/or the bonding-target C-particle(s) to be bound with each bonding-target F-particle from the bonding-target P-particle(s) and/or the bonding-target C-particle(s) associated in the third process in the ascending order of the interparticle distance to satisfy the upper limit value Ns, and
the third process and the fourth process are repeated until all of the bonding-target P-particles and the bonding-target C-particles are selected for the bonding-target F-particles.
Preferably, the F-particles of each filler model include those disposed at apexes of a polyhedron, and the bonding-target F-particles are the F-particles disposed at the apexes.
According to the present invention, the number of bonding-target P-particle(s) and/or C-particle(s) to be bound with each bonding-target F-particle does not exceed the predetermined upper limit value to prevent an excessive increase in the summation of forces of the bonding-target P-particle(s) and/or C-particle(s) exerted on each bonding-target F-particle, therefore, it is possible to perform the molecular dynamics calculation stably without causing calculation failure.
a) is a diagram showing a state where the C-particles are not yet bound with the P-particles.
b) is a diagram showing a state where a P-particle is bound with a C-particle.
Embodiments of the present invention will now be described in detail in conjunction with accompanying drawings.
The simulation method for a high-polymer material according the present invention (shortly, the simulation method) is a computer-implemented method for estimating the dispersion of filler contained in the high-polymer material.
For example, the filler may be carbon black, silica, alumina or the like. For example, the high-polymer material may be rubber, resin, elastomer or the like.
As shown in
Process S1
In the simulation method, as shown in
Each filler model 3 is a particle model made up of F-particles 4 as shown in
The term “F-particle” means a particle of a filler model 3. The F-particles 4 are arranged to accord with a predetermined polyhedron 13.
In this example, the polyhedron 13 is a regular hexahedron.
More specifically, the F-particles 4 of each filler model 3 are
a single center F-particle 4c disposed at the center of gravity of the polyhedron 13, and
a plurality of surface F-particles 4s disposed at the surface of the polyhedron 13.
The surface F-particles 4s in this example are
apex F-particles 7 disposed at the apexes 13a of the polyhedron 13, and
middle F-particles 8 disposed between the apex F-particles 7.
The center F-particle 4c and the surface F-particles 4s are each defined as a sphere having a certain diameter.
Between the center F-particle 4c and each apex F-particle 7 (surface F-particle 4s), and
between each apex F-particle 7 (surface F-particle 4s) and each of the adjacent middle F-particles 8,
links 4j having respective equilibrium lengths are defined.
Incidentally, the filler models 3 are numerical data processable with the computer 1 and stored in the computer 1. Such numerical data may include the mass, volume and diameter of each F-particle 4, and initial coordinate values of the F-particles 4.
The equilibrium lengths are equal to the bond distances between the center F-particle 4c and each surface F-particles 4s and between the surface F-particles 4s and 4s, at which the relative positions of the surface F-particles 4s on the polyhedron 13 are stably maintained.
Thereby, the arrangement of the center F-particle 4c and surface F-particles 4s of each filler model 3 can be stably maintained in the form of the polyhedron 13.
Process S2
In the simulation method, as shown in
As shown in
The term “P-particle” means a particle of a polymer model 5. The P-particles 6 are each defined as a sphere having a certain diameter.
In this example, the P-particles 6 include a native P-particle 6a and a denatured P-particle 6b for which different potentials are defined.
Between the native P-particles 6a and 6a, and between the native P-particle 6a and the denatured P-particle 6b, links 6j having respective equilibrium lengths are defined so that the polymer model 5 has a three-dimensional straight-chain structure.
Incidentally, the polymer model 5 is numerical data processable with the computer 1 and stored in the computer 1.
The numerical data may include the mass, volume, and diameter of each P-particle 6, and initial coordinate values of the P-particles 6.
Condition Setting Process S3
In the simulation method, simulation conditions are defined. (Condition setting process S3)
Process S31
In the condition setting process S3, potentials U are defined between the F-particles (4c and 4s) and the P-particles (6a and 6b) as shown in
The potential U is a function of the interparticle distance between the concerned two particles and used to calculate a force exerted between the concerned two particles.
If the interparticle distance rij is less than a predetermined cutoff distance rc, the potential U is defined by the following expression (1):
wherein
rij is the interparticle distance between the centers of the two particles concerned,
rc is the cutoff distance between the centers of the two particles concerned, and
aij is a coefficient for the intensity of the potential defined between the two particles concerned.
Thus, the potential U exerts an interaction force (in this embodiment, a repulsive force).
If the interparticle distance rij is equal to or more than the predetermined cutoff distance rc, the potential U is defined as zero, and does not exert a force (repulsive force) between the two particles.
In the first embodiment, the following potentials U1-U10 are defined between different models (3 or 5) and (3 or 5).
U1: between center F-particle 4c and native P-particle 6a
U2: between center F-particle 4c and denatured P-particle 6b
U3: between center F-particle 4c and surface F-particles 4s
U4: between surface F-particles 4s and native P-particle 6a
U5: between surface F-particles 4s and denatured P-particle 6b
U6: between native P-particle 6a and denatured P-particle 6b
U7: between center F-particle 4c and center F-particle 4c
U8: between surface F-particles 4s and surface F-particles 4s
U9: between native P-particle 6a and native P-particle 6a
U10: between denatured P-particle 6b and denatured P-particle 6b
In the first embodiment, the values of the coefficient aij (for example, 25) for the potential U2 between the denatured P-particle 6b and the center F-particle 4c, and the potential U5 between the denatured P-particle 6b and the surface F-particles 4s are set to be less than the values of the coefficient aij (for example 72) for the potential U1 between the native P-particle 6a and the center F-particle 4c, and the potential U4 between the native P-particle 6a and the surface F-particles 4s.
More specifically, the values of the coefficient aij for the potentials U1-U10 are set as follows.
potential U1: aij=72
potential U2: aij=25
potential U3: aij=50
potential U4: aij=72
potential U5: aij=25
potential U6: aij=72
potential U7: aij=50
potential U8: aij=50
potential U9: aij=50
potential U10: aij=50
These values were determined by reference to treatise—J. Chem Phys. 107(11) 4423-4435 (1997) and treatise—Macromolcule vol. 39 6744(2006).
Thereby, in comparison with the native P-particle 6a, the repulsive force exerted on the denatured P-particle 6b becomes small. Thus, the denatured P-particle 6b is defined as having a higher affinity for the F-particle (4c, 4s) than the native P-particle 6a. Therefore, the denatured P-particle 6b can model a denatured agent added in the actual high-polymer material. Accordingly, similarly to denaturalizing of the polymer in the actual high-polymer material, by incorporating the denatured P-particle 6b in the polymer model 5, it is possible to change the dispersion of the filler models 3 in the polymer models 5. Thus, the flexibility of the simulation method is increased.
The cutoff distances rc for the potentials U1-U10 are set as follows.
potential U1: rc=3
potential U2: rc=3
potential U3: rc=3
potential U4: rc=1
potential U5: rc=1
potential U6: rc=1
potential U7: rc=5
potential U8: rc=1
potential U9: rc=1
potential U10: rc=1
Each of the cutoff distances rc for the potentials (for example, potential U1) relating to the center F-particle 4c of a filler model 3 is set to be larger than the cutoff distances rc for the potentials (for example, potential U4) relating to the surface F-particles 4s of the filler model 3.
In the first embodiment, further, each of the cutoff distances rc for the potentials (for example, U1) relating to the center F-particle 4c is set to be larger than the sum (rc+Lc) of each of the cutoff distances rc for the potentials (for example, U4) relating to each surface F-particle 4s and the interparticle distance Lc between the above-mentioned surface F-particle 4s and the center F-particle 4c.
Therefore, by defining the cutoff distance rc as above in a filler model 3, the potentials (for example, U1-U3 and U7) relating to the center F-particle 4c can exert interaction forces prior to the potentials (for example, U4, U5 and U8) relating to the surface F-particles 4s.
In addition, as the center F-particle 4c is defined as a sphere having a certain diameter, the potentials (U1-U3 and U7) exert radially or in all directions from the center F-particle 4c. Therefore, in a molecular dynamics calculation performed by the computer 1, a filler model 3 is treated as being a sphere approximate to an actual shape of a filler particle, and as a result, the simulation accuracy is improved.
The computer 1 can perform the molecular dynamics calculation for the filler models 3 by the use of substantively only the potentials (for example, U1-U3 and U7) relating to the center F-particle 4c of each filler model 3, excluding particles entered into the cutoff distances rc of the surface F-particles 4s thereof.
Therefore, in the simulation method according to the present invention, the computational efficiency can be improved. Incidentally, numerical data of such potentials U are stored in the computer 1.
Process S32
Next, with the computer 1, the filler models 3 and the polymer models 5 are randomly disposed (or so defined) within a virtual space 9 having a predetermined volume. (Process S32)
As shown in
In this example, the virtual space 9 is defined as being a regular hexahedron whose each side has a length L1 of 30[σ] for example.
For example, 100 filler models 3 and 2000 polymer models 5 are randomly disposed in the virtual space 9.
Process S33
Next, in order to numerically model a compact cluster of filler particles which is liable to occur in the rubber polymer as the analysis object, a compact cluster of the filler models 3 is defined. (Process S33) In this process S33, the compact cluster of at least two filler models 3 is formed by placing the filler models 3 so that the above-mentioned interaction force between the particles occurs between adjacent filler models 3, in other words, the interparticle distance becomes less than the cutoff distance rc.
Simulation Process S4
Next, as shown in
Filler Restraining Process S41
In the simulation process S4, firstly, as shown in
Process S42
Then, targeting only the polymer models 5, the molecular dynamics calculation is performed by the computer 1.
In the molecular dynamics calculation, the Newton's equation of motion is applied to all of the polymer models 5 excluding the immovable filler models 3 on the assumption that they accord with the classical dynamics for a specified period of time in the virtual space 9.
And, the P-particles 6a and 6b are traced to obtain their coordinates at each clock time.
During the molecular dynamics calculation, conditions such as the number of each of the F-particles 4c and 4s and P-particles 6a and 6b and the volume and temperature of the virtual space 9, are kept unchanged.
Thus, in the process S42 in the first embodiment, only the polymer models 5 are distributed in the virtual space 9, while preserving the compact cluster of the filler models 3, and a stable arrangement of the polymer models 5 can be obtained.
Process S43
Next, the computer 1 checks whether or not the polymer models 5 have been sufficiently distributed. (Process S43)
If sufficiently distributed (“Y” in the process S43), then the next process S44 is performed.
If not yet sufficiently distributed (“N” in the process S43), then a unit time step is incremented (Process S45) and the process S42 is repeated.
In order to obtain the stable arrangement of the polymer models 5 through the process S42, the number of steps of the molecular dynamics calculation is preferably not less than 100, more preferably not less than 1000 but preferably not more than 1000000 if the time interval of one step is 0.05[τ] ([τ] is a unit of time). If less than 100 steps, it becomes difficult to sufficiently distribute the polymer models 5.
Therefore, the polymer models 5 are sufficiently or uniformly distributed through the simulation process S4.
Process S44
Next, all of the restrained F-particles 4 of the filler models 3 are released in order to allow the filler models 3 to move within the virtual space 9 in the after-mentioned process S47. (Process S44)
In this process S44, the computer 1 unfixes the coordinates of the F-particles 4.
Bonding Process S46
Next, the F-particles 4 and the P-particles 6 are bound.
In the bonding process S46, bonding-target F-particles 15 are selected from the F-particles 4, and bonding-target P-particles 16 are selected from the P-particles 6.
The bonding-target F-particles 15 are bound with the bonding-target P-particles 16 via links to simulate chemical bond.
In this example, the bonding-target F-particles 15 of each filler model 3 are eight F-particles 4 which are the apex F-particles 7 positioned at the apexes 13a of the polyhedron 13 as shown in
The bonding-target P-particle 16 of each polymer model 5 is, as shown in
Process S11
In the bonding process S46, an upper limit value Su for the number Sa of the bonding-target P-particles 16 to be bound with each bonding-target F-particle 15 is defined.
The upper limit value Su may be arbitrarily determined, for example, based on the total number Nb of the bonding-target F-particles 15 existing in the virtual space 9 and the total number Na of the bonding-target P-particles 16 existing in the virtual space 9.
In this example, the upper limit value Su is defined to satisfy the following conditional expression (2):
Na/Nb−1=<Su=<Na/Nb+1 (2).
In this example, each polymer model 5 has one bonding-target P-particle 16, therefore, the total number Na is equal to the total number (for example, 2000) of the polymer models 5 existing in the virtual space 9.
Each filler model 3 has eight bonding-target F-particles 15, therefore, the total number Nb is equal to the total number (for example 100) of the filler models 3 existing in the virtual space 9 multiplied by eight.
The Na/Nb means the number Sa when the bonding-target P-particles 16 are evenly bound with the bonding-target F-particles 15.
If Na=2000 and Nb=800, Na/Nb=2.5 (2.5−1=<Su=<2.5+1).
Since the Su is a positive integer, the upper limit value Su for each bonding-target F-particle 15 is defined as 2 or 3.
Selecting Process S12
In the bonding process S46, next, the bonding-target P-particles 16 to be bound with the bonding-target F-particles 15 are selected by the computer 1. (Selecting process S12)
1st Process S21
In the selecting process S12, firstly, each bonding-target P-particle 16 is associated with a bonding-target F-particle 15 to which interparticle distance Lp therefrom is smallest. (1st process S21)
In the 1st process S21, for each bonding-target P-particle 16 in the virtual space 9, a bonding-target F-particle 15 to which interparticle distance Lp from the concerned bonding-target P-particle 16 is smallest, is searched. Then, the searched-out bonding-target F-particle 15 is associated with the concerned bonding-target P-particle 16. Here, the interparticle distance Lp is the distance between centers of the concerned particles 15 and 16.
Data about such associations, for example, identification numbers of each bonding-target P-particle 16 and the bonding-target F-particle 15 to be associated therewith, and their coordinate values in relation to the virtual space 9, are stored in the computer 1.
In an example shown in
With respect to the first bonding-target F-particle 15a, three bonding-target P-particles 16 are associated therewith, which are a first bonding-target P-particle 16a, a second bonding-target P-particle 16b and a third bonding-target P-particle 16c in ascending order of the interparticle distance Lp from the concerned bonding-target F-particle 15.
With respect to the second bonding-target F-particle 15b, a single bonding-target P-particle 16 is associated therewith, which is a fourth bonding-target P-particle 16d.
In
2nd Process S22
In the selecting process S12, next, for each bonding-target F-particle 15, one or more bonding-target P-particles 16 to be bound therewith are selected from a group of the bonding-target P-particle(s) 16 associated therewith (hereinafter, the “bonding-target P-particle group”). (2nd process S22)
The selection of the bonding-target P-particle(s) 16 is made in the ascending order of the interparticle distance Lp from the concerned bonding-target F-particle 15 so that the number of the selected bonding-target P-particle(s) 16 does not exceed the above-mentioned upper limit value Su for the number Sa.
As explained above in connection with
Therefore, from the group of the bonding-target P-particles (16a, 16b and 16c) associated with the first bonding-target F-particle 15a, a number (which is, in this example, equal to the upper limit value Su) of the bonding-target P-particles (16a and 16b) are selected in their ascending order of the interparticle distance Lp.
Then, the selected bonding-target P-particles (16a and 16b) are defined as particles to be bound with the first bonding-target F-particle 15a as shown in
On the other hand, each bonding-target P-particle not selected from the group (for example, the third bonding-target P-particle 16c to which interparticle distance Lp is largest) is defined as a particle not to be bound with the first bonding-target F-particle 15a.
In the case of the second bonding-target F-particle 15b, only one bonding-target P-particle 16 (fourth bonding-target P-particle 16d) is associated therewith, and
the number of which is less than the upper limit value Su (=2). Therefore, the only one P-particle 16d is defined as a particle to be bound with the second bonding-target F-particle 15b.
In the 2nd process S22, as described above, the bonding-target P-particles 16 are selected so that the number Sa of the bonding-target P-particle(s) 16 to be bound with each bonding-target F-particle 15 does not exceed the upper limit value Su.
Data about such selections, for example, identification numbers of each bonding-target F-particle 15 and the bonding-target P-particle(s) 16 selected therefor and their coordinate values in relation to the virtual space 9, are stored in the computer 1.
Third Process S23
In the selecting process S12, next, each bonding-target P-particle 16 not selected in the 2nd process S22 is associated with another bonding-target F-particle 15. (Third process S23)
In the third process S23, firstly, for each bonding-target P-particle 16 not selected in the 2nd process S22, there is searched a bonding-target F-particle 15 to which interparticle distance Lp from the concerned not-selected bonding-target P-particle 16 is second smallest next to the smallest interparticle distance Lp of the bonding-target F-particle 15 associated, in the 1st process S21, with the concerned not-selected bonding-target P-particle 16, and the searched-out bonding-target F-particle 15 is associated with the concerned bonding-target P-particle 16.
Data about such association are stored in the computer 1.
In an example shown in
In the third process S23, there is searched a bonding-target F-particle 15 to which interparticle distance Lp from the third bonding-target P-particle 16c is second smallest next to the smallest interparticle distance Lp from the third bonding-target P-particle 16c to the first bonding-target F-particle 15a.
In this example, the second bonding-target F-particle 15b is searched out.
Then, the second bonding-target F-particle 15b is associated with the third bonding-target P-particle 16c.
Fourth Process S24
In the selecting process S12, next, for each bonding-target F-particle 15, bonding-target P-particles 16 to be bound therewith are selected from a group of the bonding-target P-particle(s) associated with the concerned bonding-target F-particle 15 in the third process S23. (Fourth process S24)
In the same way as in the 2nd process S22, the selection of the bonding-target P-particle(s) 16 is made in the ascending order of the interparticle distance Lp so that the number of the selected bonding-target P-particle(s) 16 satisfies the above-mentioned upper limit value Su for the number Sa.
In this example, for the second bonding-target F-particle 15b, the single bonding-target P-particle 16 (fourth bonding-target P-particle 16d) has been already selected as explained above.
Therefore, in order that the number Sa not exceed the upper limit value Su (=2), only the third bonding-target P-particle 16c is associated with the second bonding-target F-particle 15b. And the third bonding-target P-particle 16c is defined as a particle to be bound with the second bonding-target F-particle 15b.
If more than two bonding-target P-particles 16 are associated with the second bonding-target F-particle 15b, the selection from the group of the more than two bonding-target P-particles 16 is made in the ascending order of the interparticle distance Lp so that the resultant number Sa does not exceed the upper limit value Su.
Thus, each bonding-target P-particle 16 (not shown) to which interparticle distance Lp is more than the second smallest, is defined as a particle not bound with the second bonding-target F-particle 15b.
In the fourth process S24, as described above, the bonding-target P-particle(s) 16 to be bound with each bonding-target F-particle 15 is(are) selected from a group of the bonding-target P-particle(s) 16 associated therewith through the third process S23 so that the resultant number Sa does not exceed the upper limit value Su.
Data about such selections are stored in the computer 1.
Process S25
In the selecting process S12, next, it is checked whether or not all of the bonding-target P-particles 16 have been selected for the bonding-target F-particles 15.
If all of the bonding-target P-particles 16 have been selected, then the next process S13 is performed.
If not yet selected, then with respect to each bonding-target P-particle 16 not selected, the third process S23 and the fourth process S24 are repeated.
Therefore, through the selecting process S12, all of the bonding-target P-particles 16 are selected for the bonding-target F-particles 15.
But, in this embodiment, there is a possibility of existence of a bonding-target F-particle 15 for which no bonding-target P-particle 16 is selected finally.
Process S13
In the bonding process S46, next, the bonding-target P-particles 16 selected through the selecting process S12 are bound with the bonding-target F-particles 15 via links 10 in order that the filler models 3 and the polymer models 5 can simulate chemical bond between the filler and the polymer.
Process S47
In the simulation process S4, next, with respect to the filler models 3 and the polymer models 5, a molecular dynamics calculation is performed. (Process S47)
In the molecular dynamics calculation performed in the process S47, the Newton's equation of motion is applied similarly to the process S42.
And, the F-particles 4c and 4s of each filler model 4 (
Through the process S47, as shown in
In the bonding process S46, as explained, the number Sa of the bonding-target P-particle(s) 16 to be bound with each bonding-target F-particle 15 is limited to a value not more than the predetermined upper limit value Su, so that the summation of the forces of the bonding-target P-particles 16 exerted on each bonding-target F-particle 15 does not increase excessively during the molecular dynamics calculation, therefore, a stable molecular dynamics calculation without failure is possible.
In this embodiment, the upper limit value Su is set so as to satisfy the above-mentioned conditional expression (2), based on the number Sa of the bonding-target P-particles 16 per a bonding-target F-particle 15 when the bonding-target P-particles 16 are evenly bound with the bonding-target F-particles 15. Consequently, the bonding-target P-particles 16 can be evenly bound with the bonding-target F-particles 15, and it is possible to prevent an excessive increase in the summation of the forces of the bonding-target P-particles 16 exerted on each bonding-target F-particle 15.
In this embodiment, since the bonding-target F-particles 15 are only the apex F-particles 7 disposed at the apexes 13a of the polyhedron 13, it is possible to provide a space between the bonding-target P-particle(s) 16 bound with a bonding-target F-particle 15 and the bonding-target P-particle(s) 16 bound with an adjacent bonding-target F-particle 15.
Thereby, it is possible to prevent the interaction between the bonding-target P-particles 16 from becoming excessively increased, and it is possible to stably perform the molecular dynamics calculation.
Process S48
Next, the computer 1 judges whether or not the filler models 3 and the polymer models 5 have been sufficiently distributed. (Process S48)
If sufficiently distributed (“Y” in the process S48), then the next estimating process S5 (
If not yet sufficiently distributed (“N” in the process S48), then a unit time step is incremented (Process S49), and the process S47 is repeated.
In the simulation process S4, therefore, a structurally-relaxed high-polymer material model 11 in which the filler models 3 and the polymer models 5 are fully distributed, can be defined.
In order to effectually distribute the filler models 3 and the polymer models 5, the number of steps of the molecular dynamics calculation performed in the process S47, is preferably not less than 1000, more preferably not less than 5000 but preferably not more than 1000000 if the time interval of one step is 0.05[τ].
If the number of steps is less than 1000, there is a possibility that the filler models 3 and the polymer models 5 are not sufficiently distributed.
Estimating Process S5
As shown in
(Estimating process S5)
In this process S5, similarly to the patent document 1, with respect to only the center F-particle 4c of each filler model 3, the radial distribution function is computed.
The radial distribution function is a function which expresses a probability density that another center F-particle 4c can exist at a distance r from a certain center F-particle 4c (
In the first embodiment, the dispersion state of the filler models 3 is obtained by computing the radial distribution functions with respect to only the center F-particles 4c. Therefore, it is possible to control an increase in the computational cost.
More specifically, in this estimating process S5, with the computer 1, it is checked whether the results of the radial distribution functions are within the predetermined permissible limits or not.
If within the permissible limits, namely, the dispersion state is good, then according to the high-polymer material model 11, an actual high-polymer material is manufactured. (Process S6)
If outside the permissible limits, namely, the dispersion state is not good, then the conditions previously set to the filler models 3 and the polymer models 5 are changed, for example, taking the results of the radial distribution functions into account (process S7), and then the processes S1 to S5 are repeated.
Thereby, the structure of the high-polymer material model 11 (shown in
In the simulation method according to the first embodiment, since the structures of the filler models 3 and the polymer models 5 and the potentials can be arbitrarily defined, it is possible to estimate the properties and performances of a nonexistent unknown high-polymer material.
In the simulation method according to the first embodiment, the dispersion of filler contained in a high-polymer material is estimated, using the filler models 3 and the polymer models 5 disposed in the virtual space 9.
But, the high-polymer material is not limited to this type. The high-polymer material can be a high-polymer material further containing a coupling agent for coupling the filler with the polymer as follows.
In the following description, with respect to the processes, models and the like of the second embodiment which are the same as those in the first embodiment, the same reference numbers and characters are assigned and the descriptions are omitted in certain instances.
In the simulation method according to the second embodiment, as described above in connection with the first embodiment, the process S1 for defining filler models 3 and the process S2 for defining the polymer models 5 are performed.
Process S8
In the second embodiment, coupling agent models 17 of the coupling agent are defined. (Process S8)
The coupling agent model 17 is a particle model made up of at least one, in this example, a plurality of C-particles 18. The term “C-particle” means a particle of a coupling agent model 17.
The C-particles 18 are each defied as a sphere having a certain diameter.
Between the C-particles 18, links 21 having respective equilibrium lengths are defined.
In this example, each coupling agent model 17 has a three-dimensional straight-chain structure.
The number of the C-particles 18 of each coupling agent model 17 is smaller than the number of the P-particle 6 of a polymer model 5 in order to simulate the coupling agent having a molecular weight smaller than that of a molecular chain of the high-polymer material.
Incidentally, if a coupling agent model 17 is made up of a single C-particle 18, there is no link 21.
In this example, all of the coupling agent models 17 are made up of the identical number of the C-particles 18.
But, it is also possible to employ plural kinds of coupling agent models 17 made up of different numbers of the C-particles 18, for example, in order to estimate the effect of variations of the lengths of the links of the coupling agent on the dispersion of the filler.
Incidentally, the coupling agent models 17 are numerical data processable with the computer 1 including the mass, volume, diameter, initial coordinate values and the like of each C-particle 18 and stored in the computer 1.
Condition Setting Process S3
In the simulation method according to the second embodiment, next, simulation conditions are defined.
In this process S3, in the same way as in the first embodiment, the above-mentioned process S31 for defining potentials between the F-particles 4c and 4s and the P-particles 6a and 6b is performed.
The potentials U1-U10 between the F-particles 4c and 4s and the P-particles 6a and 6b shown in
Process S34
In the condition setting process S3 in the second embodiment, next, potentials U are defined
between the C-particles 18 and the F-particles 4c and 4s, and
between the C-particles 18 and the P-particles 6a and 6b.
The potential U is defined by the above-mentioned expression (1).
As shown, potentials U11-U14 are defined as follows.
U11: between C-particle 18 and center F-particle 4c
U12: between C-particle 18 and surface F-particles 4s
U13: between C-particle 18 and native P-particle 6a
U14: between C-particle 18 and denatured P-particle 6b
In this example, the coefficient aij for each potential U11-U14 is defined as follows.
potential U11: aij=50
potential U12: aij=50
potential U13: aij=50
potential U14: aij=50
The reason for setting the same value to the potentials U11-U14 is to even the repulsive force of the coupling agent model 17 exerted on the filler model 3 and the repulsive force of the coupling agent model 17 exerted on the polymer model 5, namely, to simulate a coupling agent in which the affinity for the filler is the same as the affinity for the polymer.
In the second embodiment, it is also possible that the coefficients aij for the potentials U11-U14 are set to different values in order to simulate a coupling agent in which the affinity for the filler is different from the affinity for the polymer, for example as follows.
potential U11: aij=10
potential U12: aij=10
potential U13: aij=100
potential U14: aij=100
In this example, the repulsive force between the filler model 3 and the coupling agent model 17 can be smaller than the repulsive force between the polymer model 5 and the coupling agent model 17, and the affinity of the coupling agent model 17 becomes higher for the filler model 3 than the polymer model.
Further, the cutoff distance rc in the expression (1) is set as follows.
potential U11: rc=3
potential U12: rc=1
potential U13: rc=1
potential U14: rc=1
The cutoff distance rc of the potential U11 relating to the center F-particle 4c of the filler model 3 is set to be larger than the cutoff distance rc of the potential U12 relating to the surface F-particles 4s of the filler model 3, and the cutoff distance rc of the potential U11 is set to be larger than the summation (rc+Lc) of the cutoff distance rc of the potential U12 and the interparticle distance Lc (see
Further, as the center F-particle 4c is expressed by the sphere having the certain diameter, it is possible to exert the potential U11 radially. Therefore, in the simulation process S4, the filler model 3 can be treated as a sphere approximate to the actual filler, and the simulation accuracy can be improved.
Process S32
In the condition setting process S3 in the second embodiment, a plurality of the filler models 3, a plurality of polymer models 5 and a plurality of coupling agent models 17 are disposed in the predetermined virtual space 9. (process S32)
As to the virtual space 9, the virtual space show in
For example, in the virtual space 9, 100 filler models 3, 1000 polymer models 5, and 1000 coupling agent models 17 are randomly disposed.
Process S33
The process S33 for forming a compact cluster of the filler models 3 is performed in the same way as in the first embodiment. (Process S33)
As to the technique for forming a compact cluster of the filler models 3, that employed in the first embodiment is employed.
Simulation Process S4
In the simulation method according to the second embodiment, next, with the computer 1, a molecular dynamics calculation is performed using the filler models 3, the polymer models 5 and the coupling agent models 17 disposed in the virtual space 9. (Simulation process S4)
Filler Restraining Process S41
In the simulation process S4 in the second embodiment, first, the filler restraining process S41 for restraining the motion of each of the F-particles 4c and 4s (shown in
Process S51
Then, with respect to the polymer models 5 and the coupling agent models 17, a molecular dynamics calculation is performed. (Process S51)
In the molecular dynamics calculation in this process S51, the Newton's equation of motion is applied to all of the polymer models 5 and the coupling agent models 17 excluding the immovable filler models 3 on the assumption that they accord with the classical dynamics for a specified period of time in the virtual space 9.
And the P-particles 6a and 6b and the C-particles 18 are traced to obtain their coordinates at each clock time.
During the molecular dynamics calculation, the conditions, such as the numbers of the F-particles 4c and 4s, P-particles 6a and 6b and C-particles 18 and the volume and temperature of the virtual space 9 are kept unchanged.
Thus, in the process S51, only the polymer models 5 and the coupling agent models 17 are distributed within the virtual space 9, while keeping the compact cluster of the filler models 3 unchanged, and a stable arrangement of the polymer models 5 and the coupling agent models 17 can be obtained.
Process S52
In the simulation process S4 in the second embodiment, next, with the computer 1, it is judged whether the polymer models 5 and the coupling agent models 17 are sufficiently distributed or not. (Process S52)
If sufficiently distributed (“Y” in the process S52), then the next process S44 is performed.
If not yet sufficiently distributed (“N” in the process S52), then a unit time step is incremented (Process S45), and the process S51 is repeated.
Therefore, in the simulation process S4 in the second embodiment, the polymer models 5 and the coupling agent models 17 can be effectually distributed.
In order to obtain a stable arrangement of the polymer models 5 and the coupling agent models 17, similarly to the first embodiment, the number of steps of the molecular dynamics calculation is preferably set to not less than 1000 if the time interval of one step is 0.05[τ].
Process S44
In the process S44, all of the restrained F-particles 4 of the filler models 3 are released in the same way as in the first embodiment.
a) and
Process S53
In the simulation process S4, next, the P-particles 6 are bound with the C-particles. (Process S53)
In this process S53, a P-particle 6 and a C-particle 18 are bound with each other via a link 20 when the interparticle distance therebetween becomes a predetermined value L3 or less to simulate a chemical bond between the coupling agent and the polymer.
Preferably, the predetermined value L3 is 0% to 200% of the cutoff distance rc (
If, with respect to one coupling agent model 17, a plurality of P-particles 6 whose interparticle distances are not more than L3 exist, then the coupling agent model 17 is bound with only one P-particle 6 of which interparticle distance is not more than the other(s), which usually has the smallest interparticle distance Lp in order that the forces of a plurality of P-particles 6 are not accumulated on one C-particle 18, and thereby to avoid calculation failure.
Bonding Process S54
In the simulation process S4 in the second embodiment, next, the F-particles 4 are bound with the P-particles 6 and the C-particles to simulate chemical bond between the filler, polymer and coupling agent. (Bonding process S54)
In this process S54, each bonding-target F-particle 15 selected from a plurality of the F-particles 4 is bound with a bonding-target P-particle 16 selected from a plurality of the P-particles 6 via a link and/or bound with a bonding-target C-particle 19 selected from a plurality of the C-particle 18 via a link.
The bonding-target F-particle(s) 15 may be arbitrarily selected from the F-particles 4 of each filler model 3.
As shown in
The bonding-target P-particle(s) 16 may be arbitrarily selected from the P-particles 6 of each polymer model 5.
As shown in
In the second embodiment, as shown in
Process S61
In the bonding process S54, there is defined an upper limit value Ns for the sum (Sa+sb) of the number Sa of the bonding-target P-particle(s) 16 and the number sb of the bonding-target C-particle(s) 19 to be bound with each bonding-target F-particle 15. (Process S61)
The upper limit value Ns may be arbitrarily determined, for example, depending on
the total number Na of the bonding-target P-particle 16 existing in the virtual space 9,
the total number Nb of the bonding-target F-particles 15 existing in the virtual space 9, and
the total number Nc of the bonding-target C-particle 19 existing in the virtual space 9.
In this example, the upper limit value Ns is defined by the following expression (3).
Ns=[(Na+Nc)/Nb+1] (3)
The total number Na is equal to the total number (in this example, 1000) of the polymer models 5 disposed in the virtual space 9.
The total number Nb is equal to the total number of the filler models 3 (=100) multiplied by the number of the bonding-target F-particles 15 per a filler model 3 (=8), namely 800.
In the second embodiment, the total number Nc is equal to the total number of the bonding-target C-particles 19 disposed in the virtual space 9, namely, 1000, since each coupling agent model 17 has a single bonding-target C-particle 19 as shown in
Thus, the (Na+Nc)/Nb means the number of the particles (16, 19) per one bonding-target F-particle 15 when the bonding-target P-particles 16 and the bonding-target C-particles 19 are evenly bound with the bonding-target F-particles 15.
In the expression (3), the [(Na+Nc)/Nb+1] means an integer obtained by rounding up the value (Na+Nc)/Nb.
For example, if Na=1000, Nb=800 and Nc=1000, then (Na+Nc)/Nb=2.5 and Ns=[(Na+Nc)/Nb+1]=3.
Selecting Process S62
In the bonding process S54 in the second embodiment, the bonding-target P-particle(a) 16 and/or bonding-target C-particle(s) 19 to be bound with each bonding-target F-particle 15 are selected. (selecting process S62)
1st Process S71
In the selecting process S62, first, each of the bonding-target P-particles 16 and bonding-target C-particles 19 is associated with one bonding-target F-particle 15 to which interparticle distance Lp is smallest. (1st process S71)
In
In the 1st process S71, for each bonding-target P-particle 16 in the virtual space 9 (shown in
a bonding-target F-particle 15 to which interparticle distance Lp from the concerned bonding-target P-particle 16 is smallest is searched, and
the searched-out bonding-target F-particle 15 is associated with the concerned bonding-target P-particle 16.
In the 1st process S71, further, for each bonding-target C-particle 19 in the virtual space 9 (shown in
a bonding-target F-particle 15 to which interparticle distance Lp from the concerned bonding-target C-particle 19 is smallest is searched, and
the searched-out bonding-target F-particle 15 is associated with
the concerned bonding-target C-particle 19.
The interparticle distance Lp is the distance between the centers of the particles concerned.
Data about such associations, for example, identification numbers of the bonding-target P-particles 16, the bonding-target C-particles 19 and the bonding-target F-particles 15, and their coordinate values in relation to the virtual space 9, are stored in the computer 1.
In the example shown in
As to the first bonding-target F-particle 15a, three bonding-target P-particles 16 and two bonding-target C-particles 19 are associated therewith.
The three bonding-target P-particles 16 are a first bonding-target P-particle 16a, a second bonding-target P-particle 16b, and a third bonding-target P-particle 16c in the ascending order of the interparticle distance Lp from the bonding-target F-particle 15a.
The two bonding-target C-particles 19 are a first bonding-target C-particle 19a, and a second bonding-target C-particle 19b in the ascending order of the interparticle distance Lp from the bonding-target F-particle 15a.
As to the second bonding-target F-particle 15b, a single bonding-target P-particle 16 (referred to as the fourth bonding-target P-particle 16d) is associated therewith.
In
2nd Process S72
In the selecting process S62 in the second embodiment, next, the bonding-target P-particle(s) 16 and/or bonding-target C-particle(s) 19 to be bound with each bonding-target F-particle 15 are selected from the bonding-target P-particle(s) 16 and/or the bonding-target C-particle(s) 19 associated with the concerned bonding-target F-particle 15. (2nd process S72)
The selection is made in the ascending order of the interparticle distance Lp so that the sum (Sa+Sb) of the number Sa of the selected bonding-target P-particle(s) 16 and the number sb of the selected bonding-target C-particle(s) 19 does not exceed the upper limit value Ns.
In this example, the upper limit value Ns is 3.
In the example shown in
Therefore, from the bonding-target P-particles 16 and bonding-target C-particles 19 associated with the first bonding-target F-particle 15a, regardless of whether P-particle 16 or C-particle 19, particles (16, 19) are selected in the ascending order of the interparticle distance Lp so that the number of the selected particles (16, 19) does not exceed the upper limit value Ns (in this example, 3).
In this example, the number of the selected particles (16, 19) is equal to Ns, and the first bonding-target P-particle 16a, the second bonding-target P-particle 16b, and the first bonding-target C-particle 19a are selected.
The selected particles (16a, 16b and 19a) are defined as particles to be bound with the first bonding-target F-particle 15a.
In
The not-selected particles (16c and 19b) are defined as particles not to be bound with the first bonding-target F-particle 15a.
As to the second bonding-target F-particle 15b, on the other hand, the number (in this example, =1) of the fourth bonding-target P-particle 16d associated with the second bonding-target F-particle 15b is under the upper limit value Ns (=3).
Therefore, as shown in
In the 2nd process S72, therefore, the bonding-target P-particle(s) 16 and/or bonding-target C-particle(s) 19 is/are selected so that the number of the bonding-target P-particle(s) 16 and/or bonding-target C-particle(s) 19 to be bound with each bonding-target F-particle 15 does not exceed the upper limit value Ns.
Data about such selections, for example, identification numbers of the bonding-target P-particles 16, the bonding-target C-particles 19 and the bonding-target F-particles 15 and their coordinate values in relation to the virtual space 9, are stored in the computer 1.
Third Process S73
In the selecting process S62 in the second embodiment, next, each of the bonding-target P-particles 16 and bonding-target C-particles 19 not selected in the 2nd process S72 is associated with another bonding-target F-particle 15.
In the process S73, for each of the bonding-target P-particles 16 and bonding-target C-particles 19 not selected in the 2nd process S72,
a bonding-target F-particle 15 to which interparticle distance Lp from the concerned not-selected particle (16, 19) is second smallest next to the interparticle distance Lp of the bonding-target F-particle 15 associated with the concerned not-selected particle (16, 19) through the 1st process S71,
is searched, and
the searched-out bonding-target F-particle 15 is associated with the concerned not-selected particle (16, 19).
Data about such association are stored in the computer 1.
In the example shown in
Therefore, as shown in
Fourth Process S74
In the selecting process S62 in the second embodiment, next, the bonding-target P-particle(s) 16 and/or bonding-target C-particle(s) 19 to be bound with each bonding-target F-particle 15 is(are) selected
from the bonding-target P-particle(s) 16 and/or bonding-target C-particle(s) 19 associated, through the third process S73, with the concerned bonding-target F-particle 15.
In the fourth process S74, similarly to the 2nd process S72, the selection is made in the ascending order of the interparticle distance Lp so that the sum (Sa+Sb) of the number Sa of the selected bonding-target P-particle(s) 16 and the number sb of the selected bonding-target C-particle(s) 19 satisfies the upper limit value Ns, namely,
does not render the total number of the resultant bound particles over the upper limit value Ns.
In the case of the second bonding-target F-particle 15b, therefore, the sum (Sa+Sb) satisfies the upper limit value Ns (=3 in this example).
Consequently, as shown in
If three or more bonding-target particles 16 and/or 19 are associated with the second bonding-target F-particle 15b, they are selected in the ascending order of the interparticle distance Lp so that the sum (Sa+Sb) satisfies the upper limit value Ns (=3 in this example).
Therefore, each of the P-particle(s) 16 and/or the C-particle(s) 19 to which interparticle distance Lp is the third smallest or more, is not defined as a particle to be bound with the second bonding-target F-particle 15b.
Thus, in the fourth process S74, the bonding-target P-particle(s) 16 and bonding-target C-particle(s) 19 to be bound with each bonding-target F-particle 15 are selected from the bonding-target P-particle(s) 16 and/or the bonding-target C-particle(s) 19 associated through the third process S73 so that the sum (Sa+sb) does not render the total number of the resultant bound particles over the upper limit value Ns.
Data about the selections of the bonding-target P-particles 16 and the bonding-target C-particles 19 are stored in the computer 1.
Process S75
In the selecting process S62 in the second embodiment, next, it is checked whether or not all of the bonding-target P-particles 16 and the bonding-target C-particles 19 have been selected for the bonding-target F-particles 15. (Process S75)
If all of them have been selected, the next process S63 is performed.
If not selected, with respect to the not selected particle(s) 16 and/or 19, the third process S73 and the fourth process S74 are repeated.
Thus, in the selecting process S62, all of the bonding-target P-particles 16 and the bonding-target C-particles 19 are selected for the bonding-target F-particles 15.
But, in this embodiment, there is a possibility of existence of a bonding-target F-particle 15 for which none of the bonding-target P-particles 16 and the bonding-target C-particles 19 is selected finally.
Process S63
Next, the bonding-target P-particle(s) 16 and/or bonding-target C-particle(s) 19 selected in the bonding process S46 are bound with the bonding-target F-particle 15.
In the process S63, each selected bonding-target P-particle 16 is bound with the bonding-target F-particle 15 via a link 22 so that the filler models 3 and the polymer models 5 can simulate chemical bond between the filler and the polymer. Further, each selected bonding-target C-particle 19 is bound with the bonding-target F-particle 15 via a link 22 so that the filler models 3 and the coupling agent models 17 can simulate chemical bond between the filler and the coupling agent.
Process S55
In the simulation process S4, next, a molecular dynamics calculation is performed. (Process S55)
In the process S55, the Newton's equation of motion is applied to all of the filler models 3, the polymer models 5 and the coupling agent models 17, on the assumption that they accord with the classical dynamics for a specified period of time in the virtual space 9.
And the F-particles 4c and 4s (shown in
Through the process S55, the filler models 3, the polymer models 5 and the coupling agent models 17 are distributed within the virtual space 9 as shown in
As explained above, in the bonding process S46 in the simulation process S4 in the second embodiment, with respect to each bonding-target F-particle 15, the sum (Sa+Sb) of the number Sa of the bonding-target P-particle(s) 16 to be bound therewith and the number sb of the bonding-target C-particle(s) 19 to be bound therewith is set to a value not more than the predetermined upper limit value Ns. Therefore, it is possible to prevent the occurrence of large variation in the number of the P-particle(s) 6 and C-particle(s) 18 bound with each F-particle 4. Thereby, in the simulation process S4 in the second embodiment, it is possible to prevent an excessive increase in the summation of the forces of the bonding-target P-particle(s) 16 and bonding-target C-particle(s) 19 exerted on each bonding-target F-particle 15 during the molecular dynamics calculation. Accordingly, in the simulation method according to the second embodiment, it is possible to avoid the occurrence of calculation failure to continue the molecular dynamics calculation stably.
In the second embodiment, as the upper limit value Ns is defined so as to satisfy the expression (3), based on the number of the particles (16, 19) per one bonding-target F-particle 15 when the bonding-target P-particles 16 and the bonding-target C-particles 19 are evenly bound with the bonding-target F-particles 15, the bonding-target P-particles 16 and the bonding-target C-particles 19 can be evenly bound with the bonding-target F-particles 15, and it is possible to prevent an excessive increase in the summation of the forces exerted on each bonding-target F-particle 15 during the molecular dynamics calculation.
In the second embodiment, since the bonding-target F-particles 15 are only the apex F-particles 7 disposed at the apexes 13a of the polyhedron 13 (shown in
Process S56
In the simulation process S4, next, it is judged with the computer 1 whether the filler models 3, the polymer models 5 and the coupling agent models 17 have been sufficiently distributed or not. (Process S56)
If sufficiently distributed (“Y” in the process S56), the next estimating process S5 (
If not yet sufficiently distributed (“N” in the process S56), then a unit time step is incremented (process S49), and the process S55 is repeated.
Therefore, in the simulation process S4 in the second embodiment, a structurally-relaxed high-polymer material model 11 in which the filler models 3 and the polymer models 5 are effectually distributed, can be defined.
In order to effectually disperse the filler models 3, the polymer models 5 and the coupling agent models 17, the number of steps of the molecular dynamics calculation performed in the process S55, is preferably not less than 1000, more preferably not less than 5000 but preferably not more than 1000000 if the time interval of one step is 0.05[τ].
Estimating Process S5
In the simulation method according to the second embodiment, next, it is checked whether or not the dispersion state of the filler models 3 is goof or not as shown in
In the estimating process S5, similarly to the estimating process S5 in the first embodiment, the radial distribution function is computed with respect to the center F-particle 4c of each filler model 3, and
it is checked with the computer 1 whether the results of the radial distribution functions are within the predetermined permissible limits or not.
If within the permissible limits, namely, the dispersion state is good, then according to the high-polymer material model 11, an actual high-polymer material is manufactured. (Process S6)
If outside the permissible limits, namely, the dispersion state is not good, then the conditions previously set to the filler models 3 and the polymer models 5 are changed, for example, taking the results of the radial distribution functions into account (process S7), and then the processes S1 to S5 are repeated.
Thereby, in the simulation method according to the second embodiment, the structure of the high-polymer material model 11 (shown in
In the simulation method according to the second embodiment, the structures of the filler models 3 and the polymer models 5 and the potentials can be arbitrarily defined, therefore, it is possible to estimate the properties and performances of a nonexistent unknown high-polymer material.
While description has been made of particularly preferable embodiments of the present invention, the illustrated embodiments should not be construed as to limit the scope of the present invention; various modifications are possible without departing from the scope of the present invention.
According to the procedures shown in
For comparison, bonding-target P-particles were bound with bonding-target F-particles randomly without the limitation to the upper limit value Su, and the molecular dynamics calculation was performed.
Each of a simulation method as Embodiment 1 and a simulation method as comparative example 1 was executed thirty times in order to check whether or not calculation failure occurred during the molecular dynamics calculation.
Each parameter was as is written in the description.
Common specifications are as follows:
filler models:
upper limit value Su: 2
In the simulation method as Comparative example 1, calculation failure was occurred during 15-time executions per 30-time executions.
In the simulation method as Embodiment 1, calculation failure was not occurred during 30-time executions and it was confirmed that the molecular dynamics calculation can be stably performed.
According to the procedures shown in
In the Embodiments 2-4, the numbers of the C-particles and the coefficients aij for the potentials U11 to U14 were as follows.
number of C-particles: 5
coefficients aij for potentials U11 to U14: 50
number of C-particles: 1
coefficients aij for potentials U11 to U14: 50
number of C-particles: 5
coefficient aij for potential U11: 10
coefficient aij for potential U12: 10
coefficient aij for potential U13: 100
coefficient aij for potential U14: 100
Thus, in Embodiment 4, the affinity of the C-particle for the F-particle was defined as being higher than the affinity of C-particle for the P-particle.
For comparison, bonding-target P-particles and bonding-target C-particles were bound with bonding-target F-particles randomly without the limitation to the upper limit value Su and the molecular dynamics calculation was performed. (comparative example 2)
Each of simulation methods as Embodiments 2 to 4 and a simulation method as comparative example 2 was executed thirty times in order to check whether or not calculation failure occurred during the molecular dynamics calculation.
Except for the above specifications for the C-particle, each parameter was as is written in the description.
Common specifications are as follows:
total number Na of bonding-target P-particles: 1000
total number Nb of bonding-target F-particles: 800
total number Nc of bonding-target C-particles: 2500
upper limit value Ns: 6
In the embodiments 2 to 4, calculation failure was not occurred and the molecular dynamics calculation was made stably. In the comparative example 2, the molecular dynamics calculation could not be completed due to calculation failure.
Radial distribution functions of Embodiments 2 and 3 are shown in
Radial distribution functions of Embodiments 2 and 4 are shown in
In
It has been known in the art that, in general, if the chain length of a coupling agent is longer, the degree of dispersion of filler becomes higher.
Therefore, the simulation method according to the present invention can accurately simulate and estimate the influence of the chain length of the coupling agent on the dispersion of the filler.
In
It has been known in the art that, in general, if a coupling agent whose affinity for the filler is the same as that for the polymer, the degree of dispersion of filler becomes high. Therefore, the simulation method according to the present invention can accurately simulate and estimate the influence of the interaction of the coupling agent on the dispersion of the filler.
| Number | Date | Country | Kind |
|---|---|---|---|
| 2014-086665 | Apr 2014 | JP | national |
| 2014-231208 | Nov 2014 | JP | national |
