CN103377306A - Systems and methods for providing a bonded-particle model in computer aided engineering analysis - Google Patents

Abstract

Systems and methods of providing bonded-particle model amongst a plurality of discrete particles representing a physical domain in a time-marching simulation to obtain numerically simulated continuum physical phenomena are disclosed. A physical domain is represented by a plurality of discrete particles. A domain of influence is assigned to each discrete particle and a bonded-particle model is created for the discrete particles. Respective bonds are established to connect each discrete particle to all other discrete particles within its domain of influence. The bonded-particle model further defines a rule for breakage of a bond. Continuum physical phenomena of the physical domain are numerically represented through a set of formula such that a time-marching simulation of the physical domain can be conducted. Physical properties include material properties and fracture energy release rate. Finally, the bonded-particle model allows size and orientation changes of each discrete particle.

Description

The system and method for particle link model is provided in computer-aided engineering analysis

Technical field

The present invention relates generally to a kind of computer-aided engineering analysis, more specifically, relate to a kind of particle link model that is provided in a plurality of discrete particles, the representative of this particle link model advances the physical definition territory that is made of hard brittle material in the simulation in the time, is used for obtaining the physical phenomenon of simulating on the numerical value.

Background technology

Many modern project analyses are carried out under the assistance of computer system.A kind of discrete element method (DEM) that is known as during this computer assisted engineering (CAE) is analyzed, it is generally used for from the motion of a large amount of discrete particle of numerical value simulation.Along with the computing power that is used for nearest neighbo(u)r classification and the progress of numerical algorithm, become possibility from the ten hundreds of particle of numerical value simulation.Nowadays, DEM is just accepted widely as a kind of effective ways that solve the engineering problem in granular and discontinuous material, particularly in grain flow, bulk material mechanics and rock mechanics.

Classical mechanics solves partial differential equation (PDE) for basic take the field of definition of continuous distribution at the hypothesis material, comprises finite element method, boundary integral method, without gridding method etc.In Other subjects, molecular dynamics (MD) has been used to determine atom and the molecule of acting force and energy, be used for the simulation from the Nano grade to the microcosmic point, but this simulation for macroscopic aspect is unaccommodated.

On the contrary, DEM provides the distinct methods that does not need for the PDE formula of continuum mechanics.Yet, have in the conventional method defective or shortcoming.Particularly, do not exist for the integrated technology that links continuum mechanics and crushed particles after being compromised at continuous medium.Having proposed many special methods, is gratifying but there is not one in these prior aries.For example, the acting force that one of them prior art supposition acts on the particle only is axial, therefore can not correctly simulate to have any relative detrusion or the field of definition of transversely deforming.

Therefore, need to provide a kind of model of the improvement in a plurality of discrete particles, the representative of this model advances the physical definition territory that is made of hard brittle material in the simulation in the time, is used for obtaining the continuous physical phenomenon of simulation on the numerical value.

Summary of the invention

The invention discloses a kind of particle link model in a plurality of discrete particles, the representative of this particle link model advances the physical definition territory that is made of hard brittle material in the simulation in the time, is used for obtaining the continuous physical phenomenon of simulation on the numerical value.

According to an aspect of the present invention, the physical definition territory is represented by a plurality of discrete particles.The domain of influence is specified by each discrete particle, and is formed for the particle link model of this discrete particle.Set up each node so that each discrete particle is being connected on the every other discrete particle within its domain of influence.This particle link model has also defined the rule of the fracture of node.

By one group of formula representative (that is, with this particle link model domination discrete particle), simulate and to be performed so that the time in physical definition territory advances on the continuous physical phenomenon in physical definition territory (for example, the mechanical property) numerical value.Physical characteristics comprises material behavior and energy to failure rate of release.This applies a formula for the acting force that calculates between every pair of discrete particle, and the potential energy in physical definition territory.Momentum conservation and the energy equilibrium in physical definition territory are saved.Finally, the particle link model allows size and the directed conversion of each discrete particle.

According to the particle link model of an embodiment of the invention take by the physical characteristics of the material of a plurality of discrete particle representatives with following characteristics as the basis:

1) characteristic of the node between every pair of discrete particle is determined by material constant, comprising: bulk modulus, modulus of shearing and density.The intensity of node is take the fracture toughness of the material of the form of energy to failure rate of release as the basis;

2) all discrete particles can be freely mobile separately in this physical definition territory;

3) in the particle link model, there is not difference operation (that is, not having PDE); And

4) for counting yield, in the particle link model, there is not integral operation.

Target of the present invention, feature and advantage will become apparent after the detailed description to its embodiment checking by reference to the accompanying drawings hereinafter.

Description of drawings

By explanation, accessory claim and accompanying drawing hereinafter, these and other features of the present invention, aspect and advantage will be better understood, in the accompanying drawing:

Fig. 1 is the process flow diagram that illustrates according to the illustrative process of the particle link model among a plurality of discrete particles an embodiment of the invention, that be provided for numerical value patrix quasi-continuous physical phenomenon.

Fig. 2 is the two-dimensional representation that the variation of size according to the embodiment of the present invention, the schematic discrete particle in time propelling simulation and orientation is shown.

Fig. 3 is the two-dimensional representation that two schematic positions according to the embodiment of the present invention, the schematic discrete particle in time propelling simulation are shown;

Fig. 4 is the two-dimensional representation that illustrates according to schematic domain of influence an embodiment of the invention, the first schematic discrete particle;

Fig. 5 is the two-dimensional representation that illustrates according to another schematic domain of influence an embodiment of the invention, the second schematic discrete particle;

Fig. 6 A be illustrate according to the embodiment of the present invention, based on the two-dimensional representation of the first schematic configuration of the discrete particle of particle link model;

Fig. 6 B is the two-dimensional representation that the second schematic configuration of the discrete particle among Fig. 6 A is shown;

Fig. 7 illustrates according to diagram an embodiment of the invention, schematic three-dimensional discrete particle (that is, spheroid);

Fig. 8 is the two-dimensional representation that illustrates according to a plurality of discrete particles in the schematic physical definition of the representative of an embodiment of the invention territory;

Fig. 9 is the calcspar that illustrates according to the projecting part of the calculation element of the embodiment that is implemented of the present invention.

Embodiment

At first referring to Fig. 1, illustrate according to the embodiment of the present invention, be provided at the process flow diagram that particle link model among a plurality of discrete particles is used for obtaining on the numerical value illustrative process 100 of continuous physical phenomenon, these a plurality of discrete particles representatives are by the physical definition territory that advances the hard brittle material in the simulation to consist of in the time.Process 100 realizes in software and preferably is understood in conjunction with other accompanying drawings.

Process 100 begins with the definition of a plurality of discrete particles of reception in the step 102, the physical definition territory that these a plurality of discrete particle representatives are made by hard brittle material.The physical definition territory can have any size or shape.Fig. 8 illustrates represents physical field of definition 800(random geometry) a plurality of discrete particle 802(circle or dish).For convenience of explanation, as used herein all examples be two dimension and particle for circular or spherical.Yet the present invention can be applied in two dimension or three-dimensional physical definition territory.And discrete particle can have except circular or spherical other different geometric configuratioies (for example, rectangle, cube etc.).Although the discrete particle 802 of Fig. 8 arranges equably, this is not necessarily in the present invention.Any other of discrete particle placed equally can the represents physical field of definition.The definition of discrete particle comprises initial position, orientation and the size of each discrete particle.

Fig. 2 illustrates to have the first directed 202(and is represented by solid arrow) and the second directed 212(represented by dotted arrow) schematic discrete particle (that is, the solid line circle 200).The first and second directed two the different orientations of discrete particle in time propelling simulation that represent, it begins zero at first constantly, and stops at future time.For example, the first orientation can be at zero " t constantly ₀" initial orientation, and the second orientation can for another constantly " t " and orientation.Perhaps, they can be at two different constantly " t ₁" and " t ₂" two different orientations.Rotation angle ω 222 is between the first and second orientations relatively.

The first and second sizes of schematic discrete particle also are shown in Fig. 2.First size illustrates with solid line circle 200 and the second size illustrates with broken circle 210.In this example, the second size 210 is than first size 200 large (that is, expanding).The second size can less than first size (that is, shrinking) (not shown).The difference that expansion/contraction angle θ 200 is illustrated between the first and second sizes.

Shown in Figure 3 is primary importance 300 and the second place 310 of the schematic discrete particle in global coordinate system.Can use this position of known design liaison.For example, vector r301 is the primary importance 300 of measuring in cartesian coordinate system 330, and motion vector u311 representative is with respect to the second place 310 of this primary importance 300.

Referring again to Fig. 1, in step 104, the domain of influence refers to each discrete particle.This can be by finishing for each discrete particle specific characteristic size.For circular or spherical field of definition, characteristic dimension can be radius.

Fig. 4 illustrates, the first schematic discrete particle 400(solid black circle) the appointment domain of influence (broken circle zone 404).The discrete particle that is arranged in the domain of influence 404 is solid line circle 402, is broken circle 412 and be arranged on the domain of influence 404 outer discrete particles.In this example, characteristic dimension is the radius 414 of the domain of influence 404.Should be pointed out that the domain of influence is volume rather than the zone in three-dimensional.The domain of influence 404 is used for the particle link model with the quantity of the contiguous discrete particle that limits any given discrete particle.Characteristic dimension can be for constant or be different for each independent discrete particle for all discrete particles.Schematic 3 d-dem particle shown in Figure 7 (that is, spherical 700).

The second schematic discrete particle 500 specifies in the domain of influence 504 among Fig. 5.Similarly, being arranged on field of definition 504 outer discrete particles is broken circle 512, and the discrete particle that is arranged in the field of definition 504 is solid line circle 502.Should be pointed out that one group of discrete particle in Fig. 4 and Fig. 5 is identical.Each discrete particle (for example, particle 400 and 500) has it self the domain of influence.

After the domain of influence was designated, the particle link model that is used for discrete particle was related with the physical characteristics in physical definition territory in step 106.The particle link model is set up each node so that each discrete particle is connected with every other discrete particle in its domain of influence.Determine node by material behavior, for example, bulk modulus, modulus of shearing, density of material and fracture toughness.In addition, the particle link model has defined the rule of node fracture.The physical characteristics in physical definition territory comprises material behavior and the energy to failure rate of release of the material in this physical definition territory.

Fig. 6 A illustrate according to the embodiment of the present invention, based on the first schematic configuration of the discrete particle interested " Pi " 620 of particle link model.For ease of explanation, the first schematic configuration has been described to be connected to three of discrete particle " Pi " 620 other discrete particles " Pj " 628 via each node 621.In the domain of influence of discrete particle " Pi " 620, can there be other discrete particles " Pj " 628 that are greater than or less than three.The position of discrete particle " Pi " 620 is by vector " r _i" 622 specify, and the position of each other discrete particle is by vector " r _j" 624 appointments.The first structure can advance for the time in the physical definition territory fragment of the special time in the simulation.

In the second schematic configuration shown in Fig. 6 B.This second schematic configuration is in the time place after the first schematic configuration.Corresponding discrete particle interested " Pi " 630 is connected with each node 631 with the discrete particle of being connected " Pj " 638.Two position vector " r _i" 632 and " r _j" 634 be used for new position.The orientation that should be pointed out that each discrete particle is different with size (being specified by dotted line in Fig. 6 B) in the first and second structures.Node 621,631 correspondingly is conditioned and obeys the fracture rule that defines in the mathematical formulae group of setting forth in below the paragraph.

Each discrete particle has its initial position as original state, orientation and size (for example, the volume in three-dimensional or the zone in two dimension).For every pair of particle in the domain of influence 404 forms node.For example use circle, the radius R 414 of the domain of influence is characteristic dimension or affects distance.

After the distortion in time propelling simulation, each discrete particle may move to the second place (motion u311), the second orientation (rotation angle ω 222) and/or the second size (expansion/contraction θ 220).

Defined hereinafter the particle link model for physical definition territory (that is, special material):

1, material behavior:

2, the discrete particle of each in the domain of influence " i " (" Pi " 620) and his contiguous discrete particle " j " (" Pj " 628)

Affect distance R

(characteristic dimension 414)

The initial volume that affects $\∀j:\left|\right|r_j-r_i\left|\right|\≤R,V_i^0=\underset{j\∈{\mathrm{\Ω}}_i^0}{\mathrm{\Σ}}{V}_j$

The current volume that affects $\∀j:\left|\right|r_j-r_i\left|\right|\≤R,V_i^t=\underset{j\∈{\mathrm{\Ω}}_i^t}{\mathrm{\Σ}}{V}_j$

Particle constant (3-D) c ^k=8 (3K-2G)/R

c ^s=8G/R

Particle constant (2-D) c ^k=6E/ (1-v)/R

c ^s=12G/R

Transition energy density (3-D) w ^Cs=5G _c/ R

w ^ck=w ^cs/4

Transition energy density (2-D) w ^Cs=3/2 π G _c/ R

w ^ck=w ^cs/5

Wherein,

Be the domain of influence of initial discrete particle " i ", and

Be the domain of influence in the time afterwards.V _jFor being arranged on the volume of the various discrete particle in the domain of influence.

1, the interaction between contiguous discrete particle comprises:

1) expands

${\mathrm{\θ}}_i=\frac{1}{V_i^t}\underset{j\∈{\mathrm{\Ω}}_i^t}{\mathrm{\Σ}}{V}_j\frac{(r_j-r_i)\·(u_j-u_i)}{(r_j-r_i)\·(r_j-r_i)}$

2) rotation

${\mathrm{\ω}}_i=\frac{1}{V_i^t}\underset{j\∈{\mathrm{\Ω}}_i^t}{\mathrm{\Σ}}{V}_j\frac{(r_j-r_i)\×(u_j-u_i)}{(r_j-r_i)\×(r_j-r_i)}$

3) distortion is decomposed

$u_{\mathrm{ij}}^k=({\mathrm{\θ}}_j+{\mathrm{\θ}}_i)(r_j-r_i)/2$

$u_{\mathrm{ij}}^r=({\mathrm{\ω}}_j+{\mathrm{\ω}}_i)\×(r_j-r_i)/2$

$u_{\mathrm{ij}}^p=u_j-u_i-u_{\mathrm{ij}}^r$

$u_{\mathrm{ij}}^n=\frac{u_{\mathrm{ij}}^p\·(r_j-r_i)}{(r_j-r_i)\·(r_j-r_i)}(r_j-r_i)$

$u_{\mathrm{ij}}^s=u_{\mathrm{ij}}^p-u_{\mathrm{ij}}^n$

4) acting force on each discrete particle

$f_i=\underset{j\∈{\mathrm{\Ω}}_i}{\mathrm{\Σ}}{f}_{\mathrm{ij}}^k+f_{\mathrm{ij}}^n+f_{\mathrm{ij}}^s$

$f_{\mathrm{ij}}^k=\frac{\sqrt{c_i^k{c}_j^k}}{\sqrt{V_i^0{V}_j^0}}\frac{V_i{V}_j{u}_{\mathrm{ij}}^k}{\left|\right|(r_j-r_i)\left|\right|}$

$f_{\mathrm{ij}}^n=\frac{\sqrt{c_i^s{c}_j^s}}{\sqrt{V_i^0{V}_j^0}}\frac{V_i{V}_j{u}_{\mathrm{ij}}^n}{\left|\right|(r_j-r_i)\left|\right|}$

$f_{\mathrm{ij}}^s=\frac{\sqrt{c_i^s{c}_j^s}}{\sqrt{V_i^0{V}_j^0}}\frac{V_i{V}_j{u}_{\mathrm{ij}}^s}{\left|\right|(r_j-r_i)\left|\right|}$

5) deformation energy density

$w_{\mathrm{ij}}^k=\frac{1}{2}(\sqrt{c_i^k{c}_j^k}{u}_{\mathrm{ij}}^k\·u_{\mathrm{ij}}^k)/[(r_j-r_i)\·(r_j-r_i)]$

$w_{\mathrm{ij}}^n=\frac{1}{2}(\sqrt{c_i^s{c}_j^s}{u}_{\mathrm{ij}}^n\·u_{\mathrm{ij}}^n)/[(r_j-r_i)\·(r_j-r_i)]$

$w_{\mathrm{ij}}^s=\frac{1}{2}(\sqrt{c_i^s{c}_j^s}{u}_{\mathrm{ij}}^s\·u_{\mathrm{ij}}^s)/[(r_j-r_i)\·(r_j-r_i)]$

6) if the rule of the fracture of node---when arriving following critical value, the node fracture:

$w_{\mathrm{ij}}^k/\sqrt{w_i^{\mathrm{ck}}{w}_j^{\mathrm{ck}}}+w_{\mathrm{ij}}^n/\sqrt{w_i^{\mathrm{cs}}{w}_j^{\mathrm{cs}}}+w_{\mathrm{ij}}^s/\sqrt{w_i^{\mathrm{cs}}{w}_j^{\mathrm{cs}}}\≥1$

When node is stretched $(u_{\mathrm{ij}}^k+u_{\mathrm{ij}}^n)\·(r_j-r_i)\≥0$

$w_{\mathrm{ij}}^s/\sqrt{w_i^{\mathrm{cs}}{w}_j^{\mathrm{cs}}}\≥1$

When node is compressed $(u_{\mathrm{ij}}^k+u_{\mathrm{ij}}^n)\&CenterDot;(r_j-r_i)<0$

Energy to failure release rate Gc is the energy that dissipates during the fracture of each unit in the new break surface zone that forms and can obtains via known technology, for example, obtains via the material behavior test.

Referring again to Fig. 1, in step 108, the continuous physical phenomenon that the numerical value that process 100 obtains the physical definition territory under loading environment is simulated (for example, breaking of hard brittle material), this loading environment is by using a plurality of discrete particles and relevant particle link model execution time to advance simulation to obtain.This advances simulation to begin (that is, starting condition) constantly and advance until satisfy predetermined condition along with time increment in each of many solutions cycle zero time.The corresponding specific time in each solution cycle.The behavior of discrete particle is by the domination of particle link model.

According to an aspect, the present invention is intended to one or more ofly can realize functional computer system described here.The example of computer system 900 shown in Figure 9.Computer system 900 comprises one or more processor, for example, and processor 904.Processor 904 is connected to inside computer system communication bus 902.For this exemplary computer system the various software embodiment has been described.After reading this explanation, it will be apparent to those skilled in the art that computer system and/or computer organization how to use other realize the present invention.

Computer system 900 also comprises primary memory 908, preferably, is random access memory (RAM), and can comprises supplementary storage 910.Supplementary storage 910 can comprise, for example, and one or more hard disk drive 912 and/or one or more mobile storage driver 914, for example floppy disk driver, tape drive and CD drive.Mobile storage driver 914 is in known manner from mobile storage unit 918 reading out datas or to mobile storage unit 918 data writings.Mobile storage unit 918, take flexible plastic disc, tape, CD etc. as representative, it reads or data writing by mobile storage driver 914.As will be understood, mobile storage unit 918 comprises computer-usable storage medium, and it has the computer software and/or data that is stored in wherein.

In substituting embodiment, supplementary storage 910 can comprise for other the similar device that allows computer program or other instruction load to computer system 900.This device can comprise, for example, and mobile storage unit 922 and interface 920.The example (for example can comprise programming box and cartridge interface (for example device seen in the video game device), mobile memory sheet, the read-only storage of erasable programmable (EPROM), USB (universal serial bus) (USB) flash memory, or PROM) with relevant socket, and other permission software and data are passed to mobile storage unit 922 and the interface 920 of computer system 900 from mobile storage unit 922.Usually, computer system 900 is by operating system (OS) software control and coordination, and operating system software carries out for example process scheduling, memory management and network and I/O service.

Also can there be the communication interface 924 that is connected to bus 902.Communication interface 924 allows software and data to transmit between computer system 900 and external unit.The example of communication interface 924 can comprise modulator-demodular unit, network interface (for example, Ethernet card), communication port, PC memory card international association (PCMCIA) slot and card etc.The software that transmits via communication interface 924 and data are with electronics, electromagnetism, light or other the form of the signal that can be received by communication interface 924.Computing machine 900 is communicated by letter in the data network based on special rules (that is, agreement) with other computing equipments.A kind of in the agreement is TCP/IP(transmission control protocol/Internet Protocol of generally using in the Internet).Usually, communication interface 924 management are maybe ressembled the data file assembling for the less packet that transmits in data network and are raw data file with the packet that receives.In addition, communication interface 924 is processed the address portion of each packet, so that its arrives the packet that computing machine 900 is gone in correct destination or interception.In this document, term " Computer Processing medium " and " computer usable medium " are generally used for referring to for example be installed in mobile storage driver 914 in the hard disk drive 912 and/or the media of hard disk.These computer programs are for being used for providing to computer system 900 device of software.The present invention is intended to this computer program.

Computer system 900 also can comprise I/O (I/O) interface 930, and it provides access display, keyboard, mouse, printer, scanner, plotting apparatus etc. to computer system 900.

Computer program (being also referred to as computer control logic) is stored in primary memory 908 and/or the supplementary storage 910 as application module 906.Computer program also can receive via communication interface 924.This computer program, when being performed, can be so that computer system 900 be carried out feature of the present invention discussed herein.Particularly, computer program, when being performed, can be so that processor 904 be carried out feature of the present invention.Correspondingly, this computer program represents the controller of computer system 900.

Realizing in the embodiments of the present invention that by use software software can be stored in the computer program and by use mobile storage driver 914, hard disk drive 912 or communication interface 924 and be loaded in the computer system 900.Application module 906 is when being performed by processor 904, so that processor 904 is carried out function described here of the present invention.

Primary memory 908 for example can load one or more application module 906(, the discrete element method), application module 906 can be carried out by one or more processor 904 in the situation that has or do not have the user's input by the I/O interface, to realize the task of expectation.In operation, in the time of one at least one processor 904 execution application module 906, the result is calculated and is stored in (that is, hard disk drive 912) in the supplementary storage 910.The state of finite element analysis reports to user with text or diagram from the detecting device that is attached to computing machine via I/O interface 930.

Although with reference to the specific embodiment of the present invention it is described, these embodiments only are schematically, not as limitation of the present invention.Those skilled in the art will contemplate multiple modification or the variation for concrete disclosed exemplary embodiment.Although for convenience of explanation, briefly shown the discrete particle of two-dimensional representation, the present invention can be applied on the three dimensional particles, for example, and spheroid.In a word, scope of the present invention is not subject to concrete exemplary embodiment disclosed herein, and the modification that can easily expect of persons skilled in the art all should be included in the scope of the spirit and scope of application and accessory claim.

Claims (14)

1. method that is provided at the particle link model among a plurality of discrete particles that the time advances represents physical field of definition in the simulation, in order to the continuous physical phenomenon that obtains to simulate on the numerical value, described method comprises:

Be received in the definition of a plurality of discrete particles of the represents physical field of definition of being made by hard brittle material in the computer system, described computer system has application module mounted thereto;

Specify the domain of influence to each discrete particle;

Will be related with the physical characteristics in described physical definition territory for the particle link model of a plurality of discrete particles, described particle link model defines each node, in order to another or more discrete particle in the domain of influence of described each discrete particle and described each discrete particle that is connected are connected, and described particle link model definition is used for determining the rule of node fracture, wherein, described physical characteristics comprises deformation and fracture energy release rate; And

By using described discrete particle to carry out the continuous physical phenomenon of simulating on the time propelling simulation acquisition numerical value in the described physical definition territory of obeying loading environment, the behavior of described discrete particle is arranged by described particle link model.

2. method according to claim 1, wherein, the definition of described a plurality of discrete particles comprises the initial position of described each discrete particle, orientation and size.

3. method according to claim 1, wherein, described particle link model also is included in the described time and advances during the simulation, and described each discrete particle is moved to the second place from primary importance; Described each particle is expanded or is contracted to the second size from first size; With the center of described each particle from the first directed spinning to the second orientation.

4. method according to claim 1, wherein, described energy to failure release rate obtains from the material behavior test of the sample of described hard brittle material.

5. method according to claim 1, wherein, described loading environment comprises tensile force, in order to pull apart the crackle in described physical definition territory.

6. method according to claim 5, wherein, described continuous physical phenomenon is included in the Crack Extension in the described physical definition territory.

7. method according to claim 1,, wherein, described rule can density based on one group of critical strain of deriving from the relative deformation between described every pair of discrete particle.

8. method according to claim 7, wherein, described one group of critical strain can density also relies on characteristic dimension and the energy to failure release rate of described each discrete particle, and wherein, described characteristic dimension is used for defining the domain of influence.

9. system that is used for the particle link model among a plurality of discrete particles that the time that is provided at advances simulation represents physical field of definition, in order to the continuous physical phenomenon that obtains to simulate on the numerical value, described system comprises:

Primary memory is used for the computer-readable code that storage is used for application module;

Be attached at least one processor of described primary memory, described at least one processor is carried out described computer-readable code in described primary memory, with so that described application module according to the following methods executable operations:

The definition of a plurality of discrete particles of the represents physical field of definition that reception is made by hard brittle material;

Specify the domain of influence to each discrete particle;

Will be related with the physical characteristics in described physical definition territory for the particle link model of a plurality of discrete particles, described particle link model defines each node, in order to another or more discrete particle in the domain of influence of described each discrete particle and described each discrete particle that is connected are connected, and described particle link model definition is used for determining the rule of node fracture, wherein, described physical characteristics comprises deformation and fracture energy release rate; And

By using described discrete particle to carry out the continuous physical phenomenon of simulating on the time propelling simulation acquisition numerical value in the described physical definition territory of obeying loading environment, the behavior of described discrete particle is arranged by described particle link model.

10. system according to claim 9, wherein, the definition of described a plurality of discrete particles comprises the initial position of described each discrete particle, orientation and size.

11. system according to claim 9, wherein, described particle link model also is included in during the described time propelling simulation, and described each discrete particle is moved to the second place from primary importance; Described each particle is expanded or is contracted to the second size from first size; With the center of described each particle from the first directed spinning to the second orientation.

12. system according to claim 9, wherein, described energy to failure release rate obtains from the material behavior test of the sample of described hard brittle material.

13. system according to claim 9, wherein, described rule can density based on one group of critical strain of deriving from the relative deformation between described every pair of discrete particle.

14. system according to claim 13, wherein, described one group of critical strain energy density also relies on characteristic dimension and the energy to failure release rate of described each discrete particle, and wherein, described characteristic dimension is used for defining the domain of influence.

Images