CN103955592A - Method for establishing multi-scale model of medicament particles during powder compaction process of long thin metal pipe - Google Patents

Abstract

The invention discloses a method for establishing a multi-scale model of medicament particles during a powder compaction process of a long thin metal pipe, and belongs to the field of computer simulation analysis. The method specifically comprises the following steps of (1) establishing a long thin metal pipe model; (2) establishing a medicament particle adding model; (3) establishing a medicament particle model; (4) establishing a long thin metal pipe compaction action model; (5) carrying out compaction simulation by utilizing discrete meta models which are established from the step (1) to the step (4). According to the method for establishing the multi-scale model of the medicament particles during the powder compaction process of the long thin metal pipe, disclosed by the invention, through establishing a multi-scale compaction model of the powder particle, the movement of the medicament particles at the inner part of the long thin metal pipe and a density changing process of the medicament particles in a compaction action are simulated by utilizing a computer, and a feasible theory basis is provided for a compaction technology of a compaction device based on impact vibration.

Description

A kind of multiple dimensioned method for establishing model of slender metal pipe powder compacting process pharmacy particle

Technical field

The present invention relates to computer simulation analysis field, particularly a kind of multiple dimensioned method for establishing model of slender metal pipe powder compacting process pharmacy particle.

Background technology

The compacting of slender metal pipe powder granule is the produce difficult point of priming system powder charge.Because priming system is higher to the density of pharmacy particle in slender metal pipe, uniformity requirement, the compacted density of slender metal pipe pharmacy particle becomes the key factor that affects priming system quality.The priming system of China still adopts the manual production pattern of single mode single-shot at present, under this production model, exists following problem: (1) homogeneity of product is poor; (2) labour intensity is high, production efficiency is low; (3) be there is to health and safety hidden danger in operating personnel.For these problems, take vertical impact vibrating compacting principle as basis, developed the tamping unit based on impact shock, this device utilizes the self gravitation of pharmacy particle to load, by the vibratory impulse effect in vertical direction, pharmacy particle is rearranged, to reach the object that improves powder charge packing.This tamping unit has solved the poor problem of homogeneity of product, has reduced labour intensity, has improved production efficiency, has realized man-machine isolation completely simultaneously, has ensured operating personnel's health and safety.

But in tamping unit operational process, cannot observe motion, the pharmacy particle variable density process of the inner pharmacy particle of slender metal pipe in compacting action.And the multi-scale Modeling method of pharmacy particle is the angle from particle, set up the compacting model of pharmacy particle, by foundation and the emulation of model, disclose macroscopic motion process, for addressing the above problem, provide effective solution route.

At present, multi-scale Modeling method has application widely in fields such as ground, mining and metallurgy, agricultural, chemical industry, pharmacy and environment, but also waits research for slender metal pipe powder compacting process pharmacy particle multi-scale Modeling.

Summary of the invention

The object of the invention is to overcome above-mentioned deficiency, the multiple dimensioned method for establishing model that a kind of slender metal pipe powder compacting process pharmacy particle is provided, it can simulate motion and the pharmacy particle variable density process of the inner pharmacy particle of slender metal pipe in compacting action.

A multiple dimensioned method for establishing model for slender metal pipe powder compacting process pharmacy particle, concrete steps are as follows:

1) set up long and thin metal tube model, adopt hollow form right cylinder metope to simulate slender metal pipe, and to set right cylinder metope inwall be effective side; In right cylinder metope bottom, bottom metope is set plug is simulated, and the upper surface of setting bottom metope is effective side;

2) set up the reinforced model of pharmacy particle, adopt funnel-form metope to simulate the reinforced model of pharmacy particle, funnel-form metope is positioned at the top of hollow form right cylinder metope, and setting funnel-form metope inwall is effective side, and sets the boundary condition that it generates for pharmacy particle;

3) set up pharmacy particle model, according to step 2) in the pharmacy particle that obtains generate boundary condition, set up linear contact stiffness model pharmacy particle model simulated;

4) setting up slender metal pipe compacting action model simulates compacting action;

5) utilize step 1) to step 4) the discrete element model set up tamps emulation, particle porosity while observing inner pharmacy particle motion and differing heights in simulation process, and adopting factor of porosity move the parameter index of rear pharmacy particle even density degree as weighing compacting, the inner respectively pharmacy particle density uniformity highly of rear slender metal pipe tamp in checking.

Further, the boundary condition judgment formula that step 2) Chinese medicine particle generates is:

R_xy+r≤R_lim；

-R_upper≤x≤R_upper；

-R_upper≤y≤R_upper；

H_cylinder≤z≤H+H_cylinder；

And $R\_\lim =\frac{(z-H\_\mathrm{cylinder})*(R\_\mathrm{upper}-R\_\mathrm{bottom})}{H}+R\_\mathrm{bottom};$

Wherein, generating particle center-of-mass coordinate is (x, y, z), particle radius is r, the reinforced model height of pharmacy particle is H, upper opening radius is R_upper, and lower openings radius is R_bottom, and long and thin metal tube model height is H_cylinder, particle barycenter and z axle base are R_xy, and the particle of same z axle height is allowed to the limit radius R _ lim generating.

Further, step 3) described in, setting up linear contact stiffness model to the concrete grammar of pharmacy particle simulation is: utilize spring simulation particulate units between elasticity, non-resilient between damper simulation particulate units, uses with the slide block of friction factor and simulates and between particle, have friction

Utilize power-displacement Indentation Law can calculate the contact force of two Interaction between particles:

$F_n^c=K_n{U}_{n}n$

In formula, for normal direction contact force; K _nfor normal stiffness coefficient; U _nfor normal direction contact displacement; N is unit normal vector;

${\mathrm{\ΔF}}_s^c=k_s{\mathrm{\ΔU}}_s$

$F_s^c=K_{\mathrm{os}}^c+\mathrm{\Δ}{F}_s^c$

In formula, K _sfor tangential contact stiffness; Δ U _sfor relative tangential displacement increment; for tangential contact force increment; tangential contact force for current time step; walk tangential contact force when last;

So the suffered F that makes a concerted effort of particle is:

${F=F}_n^c+F_s^c+\mathrm{mg}$

According to the situation of making a concerted effort of particle, can utilize Newton second law to set up law of motion equation:

F＝ma。

Further, step 4), compacting action is simulated and is included two stages:

4-1) the freely falling body stage of pharmacy particle;

4-2) identification bottom pharmacy particle layer, and it is applied to compacting acting force;

And compacting action frequency, freely falling body height and compacting power are all adjustable.

Emulation and the observation of further, step 5), tamping action include two stages:

5-1) initial equilibrium conditions emulation and observation, initial equilibrium conditions is for having fed in raw material rear initial equilibrium state;

5-2) tamp balance emulation and the observation after having moved.

Further, the emulation of described compacting action comprises loop calculation, between following two steps, circulates:

6-1) utilize power-displacement law to calculate, apply to contacting of particle and body of wall, obtain contact force;

6-2) utilize the law of motion to calculate, apply to, on each particle, obtain the displacement of particle and body of wall, form new contact.

Owing to having adopted technique scheme, the present invention has advantages of as follows:

1, the present invention is by setting up the multiple dimensioned compacting model of powder particle, utilize motion and the pharmacy particle variable density process of the inner pharmacy particle of slender metal pipe in computer simulation compacting action, for the tamping unit compacting technique based on impact shock provides performable theory foundation;

2, having realized is that the simulation that in compacting action, the inner pharmacy particle of slender metal pipe moves is visual;

3,, by emulation and observation, can improve the efficiency of actual compacting technique, thereby determine best compacting scheme.

Other advantages of the present invention, target and feature will be set forth to a certain extent in the following description, and to a certain extent, based on will be apparent to those skilled in the art to investigating below, or can be instructed from the practice of the present invention.Target of the present invention and other advantages can be realized and be obtained by instructions and claims below.

Accompanying drawing explanation

In order to make the object, technical solutions and advantages of the present invention clearer, below in conjunction with accompanying drawing, the present invention is described in further detail, wherein:

Fig. 1 is the schematic flow sheet that the multiple dimensioned model of slender metal pipe powder compacting process pharmacy particle is set up;

Fig. 2 is slender metal pipe model demonstration figure;

Fig. 3 is long and thin metal tube model;

Fig. 4 is the reinforced model of pharmacy particle;

Fig. 5 is the reinforced model boundary condition of pharmacy particle;

Fig. 6 is the reinforced model boundary condition sectional drawing of pharmacy particle;

Fig. 7 is pharmacy particle model demonstration figure;

Fig. 8 is pharmacy particle model;

Fig. 9 is computation cycles process;

Figure 10 is for repeatedly feeding in raw material and tamped rear granular model;

Figure 11 is 4 and measures ball factor of porosity change curve.

Embodiment

Below in conjunction with drawings and Examples, the invention will be further described.

Fig. 1 is the schematic flow sheet that the multiple dimensioned model of slender metal pipe powder of the present invention compacting process pharmacy particle is set up, and as shown in the figure, mainly comprises the steps: 1) foundation of long and thin metal tube model; 2) foundation of the reinforced model of pharmacy particle; 3) foundation of pharmacy particle model; 4) foundation of slender metal pipe compacting action model; 5) the multiple dimensioned model that application is set up carries out emulation and observes motion, the pharmacy particle variable density process of the inner pharmacy particle of slender metal pipe in compacting action.

Described compacting action is mainly Vertical Free falling.By vertical cylinder, there is the slender metal pipe of sheath to be promoted to certain altitude in cover, by gas pawl retaining sleeve, have the slender metal pipe of sheath, guarantee the elemental height of freely falling body, unclamp the freely falling body that can realize slender metal pipe after gas pawl.

The present invention can specifically set forth from following five parts:

1), the foundation of long and thin metal tube model

Fig. 2 is slender metal pipe model demonstration figure, Fig. 3 is long and thin metal tube model, as shown in the figure, adopts hollow form right cylinder metope to simulate slender metal pipe, and to set right cylinder metope inwall be effective side, realize and cylindrical plot metope and internal particle between interaction; In right cylinder metope bottom, bottom metope is set plug is simulated, and the upper surface of setting bottom metope is significant surface, realize and bottom metope and internal particle between interaction; Slender metal pipe is yielding, so that its rigidity arranges numbered magnitude is low compared with bottom surface infinitepiston.

2), the foundation of the reinforced model of pharmacy particle

Fig. 4 is the reinforced model of pharmacy particle, as shown in the figure, adopts funnel-form metope to simulate the reinforced model of pharmacy particle, is positioned at long and thin metal tube model top; The effective side of funnel-form metope is funnel-form metope inwall, the interaction between realization and funnel-form metope and internal particle; The reinforced model of pharmacy particle is the boundary condition that pharmacy particle generates, and border determination methods is as follows:

If generating particle center-of-mass coordinate is (x, y, z), particle radius is r, and the reinforced model height of pharmacy particle is H, upper opening radius is R_upper, lower openings radius is R_bottom, and long and thin metal tube model height is H_cylinder, and particle barycenter and z axle base are R_xy, the particle of same z axle height is allowed to the limit radius R _ lim generating, as shown in Figures 4 and 5:

Wherein particle barycenter and z axle base are:

$R\_\mathrm{xy}=\sqrt{x^2+y^2}$

Its midsagittal plane trapezoidal area S is:

$S=\frac{H(R\_\mathrm{upper}+R\_\mathrm{bottom})}{2}$

Its midsagittal plane trapezoidal upper area S1 is:

$S1=\frac{(H+H\_\mathrm{cylinder}-z)(R\_\mathrm{upper}+R\_\lim )}{2}$

Its midsagittal plane trapezoidal upper area S2 is:

$S2=\frac{(z+H\_\mathrm{cylinder})(R\_\lim +R\_\mathrm{bottom})}{2}$

According to S1, equate with S with S2 area sum:

S＝S1+S2

Can derive R_lim:

$R\_\lim =\frac{(z-H\_\mathrm{cylinder})*(R\_\mathrm{upper}-R\_\mathrm{bottom})}{H}+R\_\mathrm{bottom}$

If generating particle coordinate and radius meets:

R_xy+r≤R_lim

-R_upper≤x≤R_upper

-R_upper≤y≤R_upper

H_cylinder≤z≤H+H_cylinder

Accept, judge that this particle allows to generate;

If generate any one that particle coordinate and radius do not meet above-mentioned condition, refusal, judges that this particle does not allow to generate.

3), the foundation of pharmacy particle model

Fig. 7 is pharmacy particle model demonstration figure, with reference to Fig. 7, sets up linear contact stiffness model and realizes the simulation to pharmacy particle model, realizes between particle and particle, the interaction between particle and body of wall.

In Fig. 7, with spring and damper, represent the elasticity and the stiff character that between unit, contact.Elasticity between the particle of spring representative unit, damper represents non-resilient between particulate units, uses slide block with friction factor to represent to exist between particle friction.K in figure _nb1, K _nb2represent particle normal stiffness coefficient, K _nwrepresent body of wall normal stiffness coefficient; K _sb1, K _sb2represent particle shear stiffness coefficient, K _swrepresent body of wall shear stiffness coefficient; F represents the friction factor of particle, body of wall; C represents spring damping, and spring damping can make spring vibration stop, and does not need special setting in model emulation.

In linear contact stiffness model, suppose the effect of connecting with the rigidity of body of wall of two contacts particles or particles, so normal stiffness COEFFICIENT K _nbb, K _nbwwith tangential contact stiffness K _sbb, K _sbwbe calculated as:

$K_{\mathrm{nbb}}=\frac{K_{\mathrm{nb}1}*K_{\mathrm{nb}2}}{K_{\mathrm{nb}1}+K_{\mathrm{nb}2}}$

$K_{\mathrm{sbb}}=\frac{K_{\mathrm{sb}1}*K_{\mathrm{sb}2}}{K_{\mathrm{sb}1}+K_{\mathrm{sb}2}}$

$K_{\mathrm{nbw}}=\frac{K_{\mathrm{nb}1}*K_{\mathrm{nw}}}{K_{\mathrm{nb}1}+K_{\mathrm{nw}}}$

$K_{\mathrm{sbw}}=\frac{K_{\mathrm{sb}1}*K_{\mathrm{sw}}}{K_{\mathrm{sb}1}+K_{\mathrm{sw}}}$

Make firmly-displacement Indentation Law can calculate the contact force of two Interaction between particles:

$F_n^c=K_n{U}_{n}n$

In formula, for normal direction contact force; K _nfor normal stiffness coefficient; U _nfor normal direction contact displacement; N is unit normal vector.

Tangential contact force is calculated with the form of increment.When contact forms, total tangential contact force is initialized as zero, then each relative tangential displacement increment can produce the tangential contact force increment of elasticity, the tangential contact force existing when step starts when new tangential contact force equals current and tangential elastic connecting touch increment sum,

${\mathrm{\ΔF}}_s^c=k_s{\mathrm{\ΔU}}_s$

$F_s^c=K_{\mathrm{os}}^c+\mathrm{\Δ}{F}_s^c$

In formula, K _sfor tangential contact stiffness; Δ U _sfor relative tangential displacement increment; for tangential contact force increment; tangential contact force for current time step; walk tangential contact force when last, so the suffered F that makes a concerted effort of particle is:

${F=F}_n^c+F_s^c+\mathrm{mg}$

According to the situation of making a concerted effort of particle, can use Newton second law to set up law of motion equation:

F＝ma。

Fig. 8 is pharmacy particle model, and according to the production technology of tamping unit, pharmacy particle each time generates just in time for metal tube allows spatial altitude.

4), slender metal pipe compacting action model is set up

Simulation to slender metal pipe compacting action is divided into two stages: the 1. freely falling body of pharmacy particle after promoting; 2. identify bottom pharmacy particle layer, and it is applied to compacting acting force.It tamps action frequency, freely falling body height, compacting power is all adjustable.

5) the multiple dimensioned model emulation of particle that, application is set up is also observed

The multiple dimensioned model of particle set up of application carries out emulation, and simulation process is divided into two stages: 1. initial balance process; 2. tamp the equilibrium process after having moved.Fig. 9 is simulation calculation cyclic process, comprising: 1) utilize power-displacement law to calculate, apply to contacting of particle and body of wall, obtain contact force; 2) utilize the law of motion to calculate, apply to, on each particle, obtain the displacement of particle and body of wall, form new contact.Between above two steps, cycle calculations is the motion process of the inner pharmacy particle of observable.

The size of factor of porosity can embody the compaction rate of pharmacy particle, so select factor of porosity to move the parameter index of rear pharmacy particle even density degree as weighing compacting.In simulation process in order to reflect the homogeneity of inner each height pharmacy particle compacted density of right cylinder internal face, in 4 height, place 4 and measure ball, and choose two of bottoms and measure ball as the basis of variable element compacting effect research, by change, tamp number of times and tamp and highly carry out difference simulation, obtaining tamping the highly impact effect to slender metal pipe inside pharmacy particle density of number of times and compacting.

Factor of porosity refers in pharmacy particle, the ratio of volume of voids and whole prose style free from parallelism volume, as formula:

$n=\frac{V_0}{V_0+V_1}$

In above formula: n is factor of porosity; V ₀for volume of voids, the m of unit ³; V _lfor prose style free from parallelism volume, the m of unit ³.

Figure 10 is for repeatedly feeding in raw material and tamped rear model.Figure 11 is 4 measurement ball factor of porosity change curves, and wherein respectively highly measuring ball porosity curve marks respectively 1,2,3,4 for 0.05m, 0.13m, 0.21m, 0.28m.

Embodiment mono-:

Step 1: the foundation of long and thin metal tube model

Set up hollow form right cylinder metope slender metal pipe is simulated, the effective side of right cylinder metope is set and points to cylindrical metope inwall, this inwall metope is for No. ID 1; In cylindrical metope bottom, set up bottom metope plug is simulated, effective side that bottom metope is set is upper surface, and this upper surface metope is for No. ID 2; Slender metal pipe is yielding, so its rigidity magnitude setting level is low compared with ground infinitepiston.As shown in Figure 3, its design parameter is as follows for long and thin metal tube model:

Cylindrical metope height is H_cylinder=0.5m;

Cylindrical metope bottom surface radius is rad=0.05m;

Step 2: the foundation of the reinforced model of pharmacy particle

Set up funnel-form metope the reinforced model of pharmacy particle is simulated, its effective side is set and points to funnel-form metope inwall, this inwall metope is for No. ID 3; As shown in Figure 4, its design parameter is as follows for the reinforced model of pharmacy particle:

The high H=0.2m of funnel-form metope;

Upper opening radius is R_upper=0.1m;

Lower openings radius is R_bottom=0.04m;

The limit radius R_lim=0.3* (z-0.5)+0.04 of the particle that particle (x, y, z) being allowed on z axle sustained height generates; So the boundary condition that pharmacy particle generates is following formula:

R_xy+r≤0.3*(z-0.5)+0.04

-0.1≤x≤0.1

-0.1≤y≤0.1

0.5≤z≤0.7

Pharmacy particle generates in the reinforced model scope of above-mentioned pharmacy particle, completes altogether repeatedly reinforced action.

Step 3: the foundation of pharmacy particle model

Set up linear contact stiffness model pharmacy particle model is simulated, the particle attribute repeatedly generating is identical, pharmacy particle model initial density and actual loose shape pharmacy particle density is set and approaches, and meet the boundary condition in step 2, and it increases progressively for No. ID successively.As shown in Figure 8, its attribute design parameter is as follows for pharmacy particle model:

Pharmacy particle model radius is r=0.004m;

Pharmacy particle model initial density is dens=500kg/m ³.

Step 4: slender metal pipe compacting action model is set up

The problem of modelling of compacting action is divided into two stages: the 1. freely falling body of pharmacy particle after promoting, by the initial velocity of the certain step number of pharmacy particle model to hollow form right cylinder inside, realize; 2. identify bottom pharmacy particle layer, radius r=0.004m due to generated granular model, so judgement granular model barycenter is at z direction of principal axis particle swarm for contact bottom surface infinitepiston in the scope of 0≤z≤0.008, naming this particle swarm is z_bottom, and it is applied to compacting acting force.But after completing once compacting action, due to the motion of particle, need to delete this particle swarm at every turn, until next time, when moving, compacting redefines particle swarm, if delete not in time, this particle swarm quantity can increase along with the increase of compacting number of times, causes abnormal compacting power transmission.The common simulation realizing slender metal pipe compacting action of these two processes.It tamps action frequency, freely falling body height, compacting power is all adjustable.

Step 5: the multiple dimensioned model of particle that application is set up carries out emulation and observation

The multiple dimensioned model of particle set up of application carries out emulation, and simulation process is divided into two stages: 1. initial balance process; 2. tamp the equilibrium process after having moved.Figure 10 has tamped rear granular model for repeatedly reinforced, with reference to Figure 10, and in simulation process, the motion process of the inner pharmacy particle of observable after cycle calculations.4 differing heights select respectively 1 pharmacy particle as measuring ball, its z direction of principal axis height is respectively 0.05m, 0.13m, 0.21m, 0.28m, the measurement radius of a ball is 0.04m, and that respectively highly measures ball is respectively 1,2,3,4 No. ID, and corresponding porosity curve marks respectively 1,2,3,4.Select bottom 0.05m, two of 0.13m to measure the basis that ball carries out the research of variable element compacting effect, by change, tamp number of times and tamp and highly carry out difference simulation, obtaining tamping the highly impact effect to slender metal pipe inside pharmacy particle density of number of times and compacting.

4 homogeneitys of measuring the inner differing heights pharmacy particle of ball factor of porosity reflection cylinder wall compacted density, as shown in figure 11.

Figure 11 is 4 and measures ball factor of porosity change curve, with reference to Figure 11, after having fed in raw material, due to the effect of institute's adding medicine particle gravity, can make bottom pharmacy particle factor of porosity further reduce at every turn; Once compacting action of the most advanced and sophisticated expression of each curve, after each compacting has been moved, pharmacy particle factor of porosity all reduces, i.e. and pharmacy particle density increase shows that compacting action can make pharmacy particle density raise; Z direction of principal axis height is larger compared with low measurement ball factor of porosity than z direction of principal axis height compared with high measurement ball factor of porosity, shows that pharmacy particle compacted density successively decreases from bottom to top successively.The factor of porosity data that derive after reinforced compacting are as shown in table 1:

The factor of porosity tables of data that table 1 is reinforced after having tamped

Measure ball height	0.05m	0.13m	0.21m	0.28m
					Factor of porosity	0.318	0.336	0.341	0.347

Wherein z direction of principal axis height is that the factor of porosity average of 0.05m, 0.13m, 0.21m, 0.28m is 0.3355, variance is 0.0001175, variance is less, showing to feed in raw material, it is more even to have tamped the factor of porosity that rear z direction of principal axis height is 0.05m, 0.13m, 0.21m, 0.28m, and each z direction of principal axis height pharmacy particle density is comparatively even.

Finally explanation is, above embodiment is only unrestricted in order to technical scheme of the present invention to be described, although the present invention is had been described in detail with reference to preferred embodiment, those of ordinary skill in the art is to be understood that, can modify or be equal to replacement technical scheme of the present invention, and not departing from aim and the scope of the technical program, it all should be encompassed in the middle of claim scope of the present invention.

Claims (6)

1. a multiple dimensioned method for establishing model for slender metal pipe powder compacting process pharmacy particle, is characterized in that, concrete steps are as follows:

1) set up long and thin metal tube model, adopt hollow form right cylinder metope to simulate slender metal pipe, and to set right cylinder metope inwall be effective side; In right cylinder metope bottom, bottom metope is set plug is simulated, and the upper surface of setting bottom metope is effective side;

2) set up the reinforced model of pharmacy particle, adopt funnel-form metope to simulate the reinforced model of pharmacy particle, funnel-form metope is positioned at the top of hollow form right cylinder metope, and setting funnel-form metope inwall is effective side, and sets the boundary condition that it generates for pharmacy particle;

3) set up pharmacy particle model, according to step 2) in the pharmacy particle that obtains generate boundary condition, set up linear contact stiffness model pharmacy particle model simulated;

4) setting up slender metal pipe compacting action model simulates compacting action;

5) utilize step 1) to step 4) the discrete element model set up tamps emulation, particle porosity while observing inner pharmacy particle motion and differing heights in simulation process, and adopting factor of porosity move the parameter index of rear pharmacy particle even density degree as weighing compacting, the inner respectively pharmacy particle density uniformity highly of rear slender metal pipe tamp in checking.

2. the multiple dimensioned method for establishing model of a kind of slender metal pipe powder compacting process pharmacy particle as claimed in claim 1, is characterized in that step 2) the boundary condition judgment formula that generates of Chinese medicine particle is:

R_xy+r≤R_lim；

-R_upper≤x≤R_upper；

-R_upper≤y≤R_upper；

H_cylinder≤z≤H+H_cylinder；

And $R\_\lim =\frac{(z-H\_\mathrm{cylinder})*(R\_\mathrm{upper}-R\_\mathrm{bottom})}{H}+R\_\mathrm{bottom};$

Wherein, generating particle center-of-mass coordinate is (x, y, z), particle radius is r, the reinforced model height of pharmacy particle is H, upper opening radius is R_upper, and lower openings radius is R_bottom, and long and thin metal tube model height is H_cylinder, particle barycenter and z axle base are R_xy, and the particle of same z axle height is allowed to the limit radius R _ lim generating.

3. a kind of slender metal pipe powder as claimed in claim 1 is tamped the multiple dimensioned method for establishing model of process pharmacy particle, it is characterized in that, step 3) described in, setting up linear contact stiffness model to the concrete grammar of pharmacy particle simulation is: utilize spring simulation particulate units between elasticity, non-resilient between damper simulation particulate units, use with the slide block of friction factor and simulate between particle and exist and rub

Utilize power-displacement Indentation Law can calculate the contact force of two Interaction between particles:

$F_n^c=K_n{U}_{n}n$

In formula, for normal direction contact force; K _nfor normal stiffness coefficient; U _nfor normal direction contact displacement; N is unit normal vector;

${\mathrm{\ΔF}}_s^c=k_s{\mathrm{\ΔU}}_s$

$F_s^c=K_{\mathrm{os}}^c+\mathrm{\Δ}{F}_s^c$

In formula, K _sfor tangential contact stiffness; Δ U _sfor relative tangential displacement increment; for tangential contact force increment; tangential contact force for current time step; walk tangential contact force when last;

So the suffered F that makes a concerted effort of particle is:

${F=F}_n^c+F_s^c+\mathrm{mg}$

According to the situation of making a concerted effort of particle, can utilize Newton second law to set up law of motion equation:

F＝ma。

4. the multiple dimensioned method for establishing model of a kind of slender metal pipe powder compacting process pharmacy particle as claimed in claim 1, is characterized in that step 4) described in compacting action is simulated and is included two stages:

4-1) the freely falling body stage of pharmacy particle;

4-2) identification bottom pharmacy particle layer, and it is applied to compacting acting force;

And compacting action frequency, freely falling body height and compacting power are all adjustable.

5. the multiple dimensioned method for establishing model of a kind of slender metal pipe powder compacting process pharmacy particle as claimed in claim 1, is characterized in that step 5) in tamp emulation and the observation of moving and includes two stages:

5-1) initial equilibrium conditions emulation and observation, initial equilibrium conditions is for having fed in raw material rear initial equilibrium state;

5-2) tamp balance emulation and the observation after having moved.

6. the multiple dimensioned method for establishing model of a kind of slender metal pipe powder compacting process pharmacy particle as claimed in claim 1, is characterized in that, the emulation of described compacting action comprises loop calculation, between following two steps, circulates:

6-1) utilize power-displacement law to calculate, apply to contacting of particle and body of wall, obtain contact force;

6-2) utilize the law of motion to calculate, apply to, on each particle, obtain the displacement of particle and body of wall, form new contact.