CN109635439A - Network chain effect research method is produced under bumpy weather based on transmission dynamics - Google Patents
Network chain effect research method is produced under bumpy weather based on transmission dynamics Download PDFInfo
- Publication number
- CN109635439A CN109635439A CN201811521510.8A CN201811521510A CN109635439A CN 109635439 A CN109635439 A CN 109635439A CN 201811521510 A CN201811521510 A CN 201811521510A CN 109635439 A CN109635439 A CN 109635439A
- Authority
- CN
- China
- Prior art keywords
- node
- network
- manufacture
- disturbance
- production
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- General Factory Administration (AREA)
Abstract
The present invention relates to workshop chain effect analysis technical fields, it is that network chain effect research method is produced under a kind of bumpy weather based on transmission dynamics, using the manufacture node in production system as node, process relation between production activity is to connect side, the parameter in propagation model is determined according to the production network of foundation, the coupling parameter of task inflow and outflow, using the linear theory of 5W2H analytic method combined process, assessment judgement is carried out to strength of turbulence according to delay time at stop and economic loss, disturbing influence degree is divided into slight perturbations, moderate disturbance, severe disturbs and four kinds of serious disturbance, node different in network is selected as source node is propagated to be activated, correlation between production network topology structure and propagation of disturbance kinetic characteristics is furtherd investigate, it is unique chain to disclose job shop under bumpy weather Propagation law.The present invention provides effectively reference and foundation for the even running control of workshop.
Description
Technical field
The present invention relates to workshop chain effect analysis technical fields, are a kind of bumpy weathers based on transmission dynamics
Lower production network chain effect research method.
Background technique
Personalized customization reduces cost as manufacturing enterprise, the important means of raising competitiveness is fallen over each other to adopt by vast enterprise
With.But the change of this production method also will directly induce rush order insertion, material covers shortage together, the customer delivery phase is changed etc.
Disturbance event occurrence frequency sharply increases.Even if working condition is by small disturbance during workshop operation, with life
The propulsion of production process, these disturbances also can constantly be passed, add up, amplifying, and eventually lead to the very big delay of manufacturing schedule, production
The rapid decrease of stability also will cause the confusion and interruption of production process when serious.Therefore, research manufacturing shop is in disturbance ring
Transmission dynamics behavior and its evolving trend under border for accurately implementing prevention, reclamation activities, and then guarantee workshop
Even running is particularly important, is increasingly becoming the focal issue of recent domestic focus of attention.
Due to the complexity of production system, it is uncertain the features such as, under bumpy weather, the exception of any manufacture link is made
At disturbance can all depend on certain carrier and propagated and spread, cause production system from the mistake operated normally to collapse interrupt
Journey can undergo intermediate state several different, this brings to production process propagation of disturbance mechanism study model foundation and its solution
Great challenge.Most of researchs, which primarily focus on, at present inquires under bumpy weather, and certain production attributes are by external disturbance
Or when Self-variation, fluctuating characteristic and changing rule of the manufacturing parameter under the effect of technique incidence relation, for disturbing to system
Whether the influence and production system chain effect for making location mode, which can allow production system to be continually maintained in stable state, still lacks
In-depth study.
Therefore complex network transmission dynamics theory is mainly applied to production system by the present invention, to production network topology knot
Correlation between structure and propagation of disturbance kinetic characteristics is furtherd investigate, unique to disclose job shop under bumpy weather
Chain propagation law.Effectively reference and foundation are provided for the even running control of workshop.
Summary of the invention
In order to solve the problems existing in the prior art, the invention proposes raw under a kind of bumpy weather based on transmission dynamics
Produce network chain effect research method.
The present invention is realized by following technical measures:
Network chain effect research method is produced under a kind of bumpy weather based on transmission dynamics, chain propagation scale is calculated
Steps are as follows for method:
It is as follows to construct propagation model:
Step 1: using the manufacture node in production system as node, the process relation between production activity is to connect side, processing
Time is side right wij, determine adjacency matrix A, constitute the production system network under complex network visual angle, generate network；
Step 2: the parameter in propagation model is determined according to the production network of foundation, the coupling parameter of task inflow and outflow,
ε1、ε2, manufacture Node B threshold Yi, calculate node parameter
Step 3: the linear theory of 5W2H analytic method combined process is used, it is strong to disturbing according to delay time at stop and economic loss
Degree carry out assessment judgement, by disturbing influence degree be divided into slight perturbations, moderate disturbance, severe disturbance and four kinds of serious disturbance, it is right
I, II, III, IV 4 grade is answered, node different in network is selected as source node is propagated and is activated；
Step 4: all manufacture node states being updated by the parameter setting in propagation model formula (2) and step 2, if xi
(t+1)≥Yi, then manufacture node and be activated and switch to propagate state, while renewal time step-length；
Step 4.1: initialization node state value sums to the side right of node i and precedence activities, obtains the original negative of node
It carries L (i), the working ability of node is that the build-in attribute of node is Capality (i), and definition node original state is
Step 4.2: in addition t0Moment adds the load Δ L (i) of node, recalculates the state of node
The downstream adjacent node for finding node i, obtains downstream node set U (i), obtain node set interior joint with thereon
Swim the connection weight w of nodeji+With enter intensityAnd the connection weight w with node downstreami+jIntensity outAnd handle will
Process is saved in interim set Temp downstream, uses formula
It calculates
In t0The state value of+1 moment downstream node.When the node traverses of set U (i) are completed, using Temp as the influence of next stage
Node set, until whole network state reaches stable state；
Step 4.3: node being repaired, according to the working ability Capality (i) of node itself and all nodes
The maximum CapalityMax (i) of middle working ability ratio as the probability repaired, the working ability of node is improved 0
Random number between~20%；
Step 4.4: the quantity of failure node in each moment network is counted,
Step 4.5: choosing the angle value of node, intensity, betweenness and gather four network attributes of coefficient, in same disturbance situation
Under, it is arranged by the descending of each network attribute, identical disturbance successively is added to node, analysis heterogeneous networks characteristic is for the scale that fails
Influence, under different disturbances, analyzing the failure scale of network is influenced by disturbing；
Step 5: entire step 4 is recycled, until the state value of all nodes is respectively less than Node B threshold, i.e. x in networki(t+
1)<Yi, i=1,2 ..., n, the chain propagation effect of production system caused by disturbance factor terminates；
Wherein: xi(t+1) production status of the manufacture node i in moment t+1, x are indicatedi(t) indicate manufacture node i in moment t
Production status, ε1∈ (0,1) is the stiffness of coupling that task flows into, ε2∈ (0,1), for the stiffness of coupling of task outflow, ε1+ε2<
1, function f (x) indicate the individual node own dynamics behavior in network, and physical significance may be characterized as resource node production energy
The Evolution of force constraint selects Logistic mapping: f (x)=4x (1-x) herein, as 0≤x≤1,0≤f (x)≤1,
N1For the sum of manufacture node in-degrees all in network, N2For the sum of the out-degree of manufacture nodes all in network, A=a (i, j)N×NFor
Adjacency matrix describes the coupled relation that node is manufactured in network, if having Bian Xianglian, a between manufacture node i and manufacture node jij
=1, otherwise aij=aji=0, do not allow to manufacture node and is connected with itself, wijNeighbours, which are directed toward, for manufacture node i manufactures node j's
Side right,To enter intensity, indicate manufacture node i enters the sum of side right, To go out intensity, manufacture node i is indicated
The sum of side right out,Absolute value guarantees to manufacture in formula node production status not less than 0.
In the present invention, manufacture node i is disturbed in t moment, and i is activated at the t+1 moment, manufactures the association of node i
Manufacture node, which is manufactured node i state by the t+1 moment, to be influenced, and the state of these manufacture nodes involved is based on formula (2)
It calculates, state value is if it is greater than or equal to threshold value Yi, will lead to more manufacture node failures, such cycle calculations workshop is in bumpy weather
Under sprawling and chain effect, the polymorphism of production system is described, four kinds of states of production system under bumpy weather are portrayed
Evolution process and disturbing influence under manufacture the state of node and constitute, disclose disturbance factor propagation row in process of production
For.Embodiment:
(1) data preparation
1) by taking certain agricultural machinery product enterprise C workshop data as an example, the chain effect that may cause under bumpy weather is ground
Study carefully.The workshop shares 36 manufacturing cells, mainly include fiber cutter, 3cm (6cm) cutter, 6cm bending machine, press machine,
Numerical control press, common punch press etc., are expressed as Ri, i={ 1,2 ..., 36 }.When table 1 is the processing and manufacturing unit of each workpiece, processing
Between and time.
Table 1 work time piece period t/min and time r/min
2) actual job-shop scheduling problem is solved, is obtained using makespan as the Optimized Operation side of target
Case: include machine sum m, workpiece sum n, the process letter processed on s-th of order of k-th of machine in initial schedule scheme
Breath: operation number, workpiece number, process time started, process end time；The process of each workpiece is made of m procedure.
(2) production system network model is established
According to Optimized Operation scheme, process is obtainedManufacturing recourses m and machine dimension successor activities machine dimensionWith
The successor activities of workpiece dimensionProcess is mapped to manufacturing recourses, occupies manufacture node j to manufacture the task of node i output
Time as the weight between node, to manufacture threshold value of the working ability as node of node, being formed with manufacturing cell is section
Point, process constraint are the production system network on even side.
(3) network parameter of network model is obtained
According to established network model, the network attribute of statistics manufacture node, node degree indicates the number of node i associated nodes
Amount,N is number of nodes, aijIndicate that node i, with the presence or absence of being associated with, exists for a with jij=1, otherwise aij=0；Section
Point intensitywijFor node i and j side right.Node betweenness indicates any two node shortest path warp in network
The number of the node i is crossed,δtShow all shortest path numbers from vertex s to vertex t, and δst(r) it indicates
From vertex s to vertex t by all shortest path numbers of vertex r.C (i)=2Ei/(ki(ki-1)),kiIt indicates to close with node i
The number of nodes of connection, EiIndicate that there are associated quantity for node i associated nodes.
Table 2 produces network part node topological property statistics
(4) the chain communication process analysis of network is produced
The CML established the production chain propagation model of network is subjected to programming using python, from propagation of disturbance source
Node type, initial disturbance intensity analyze propagation of disturbance effect.
1) influence of the strength of turbulence size to the chain propagation of network
Consider network topology structure, side right and industry characteristics tripartite's region feature, chooses angle value maximum node, enters maximum intensity
Node goes out maximum intensity node and random node as propagation of disturbance source, while level of disruption is arranged and disturbs from slight perturbations, moderate
Dynamic, severe is disturbed to serious disturbance and is successively increased, and analyzes its chain propagation law.
When disturb initial value it is smaller when, propagate source node itself repair power can control disturbance propagate, fail scale domination
Below 10%.As disturbed value increases, propagates the chain propagation scale that source node causes for four kinds and show a increasing trend.
When maximum intensity node is as source node is propagated out, causes the failure of network size to be significantly larger than other nodes and make
To propagate destruction caused by source node.Speed is spread in the chain reaction that node (fiber cutter) failure of maximum intensity causes out
Degree is most fast, it is largest to involve, and followed by enters the node (machining center) of maximum intensity.Different type is propagated caused by source node
Failure scale and manufacture node itself capacity constraints and with the close phase of characteristic attributes such as the differentiation of neighbor node business association
It closes.Therefore in production activity, the big node of intensity is especially gone out to key node should put into more guarantees, reclamation activities, with
Containment disturbance is a wide range of to propagate, and causes heavy losses.
2) influence of the characteristic attribute of source node to the chain propagation of network is propagated
In order to which the discussion being more clear manufactures node unique characteristics attribute to the scale and process of the production chain failure of network
Influence, be provided with the experiment of three groups of simulation calculations, respectively analysis node angle value, enter intensity, go out intensity and the production chain mistake of network
The correlation of effect scale.
Influence of the angle value to the production chain failure process of network:
When identical disturbed value acts on different angle value nodes, the identical node of most of angle value causes the deviation of failure scale
± 10% or so.And as angle value increases, the chain propagation scale of initiation shows a increasing trend, and the big propagation source node of angle value is more
It is also easy to produce the propagation of the large area of disturbance, production system is caused to be paralysed, meets universal law.But its interior joint 2 and node 29 are
Apparent abnormal point.Analyze from network perspective: node 2 is in entirely production network, although its in-degree, the local features category such as out-degree
Property does not protrude, but has significant global property index, i.e., highest to gather coefficient.This shows node 2 in overall network
In criticality it is very high, and node is more downstream, and business association is strong, which changes the shadow to its associated nodes
It is big to ring power, therefore node 2 is big as caused network chain propagation scale when propagating source node, the very Gao Shui in whole network
It is flat.And for node 29, although its angle value is in by-level in overall network, its global convergence factor, betweenness are located
In reduced levels, and in entirely production network, subsequent association node only has node 50, and business association is weak, the influence of disturbance
Range is small, therefore node 29 is used as caused chain propagation small scale when propagating source node.
From production angle analysis: node 2 is machining center, and node 29 is boring machine.In actual production, since machining center has
Have preferable manufacture flexible, therefore the equipment be in a plurality of process route and flows through intertwined point in production, can and other equipment production
Raw very strong business association, and its occupancy is up to 93%, show the equipment resist disturbance buffer capacity it is weaker.Therefore,
Once such resource is disturbed, it can greatly may cause the large area sprawling of disturbance and propagate, reduce the production of entire producing line
Efficiency.And boring machine is limited by working ability and business association limits, under same disturbance, influence power is smaller, therefore caused
Chain propagation scale is limited.
Enter influence of the intensity to the production chain failure process of network:
Enter the network failure scale that intensity node causes to identical disturbance, difference to emulate, identical disturbed value acts on
When difference enters intensity node, with the increase for entering intensity value, the chain propagation scale of initiation shows a increasing trend, that is, undertakes task amount
After higher node is disturbed, the chain propagation scale of initiation also can be higher, this meets common cognition rule, but its interior joint
1 and node 20 be obvious abnormal point, it is only 80 that node 1, which enters intensity, but causes 85% failure scale, and node 20 enters intensity and is
178, cause the failure scale of network but very small, only 23.5%.
It is analyzed from network perspective: as can be seen from Table 2, intensity is 230 out, node betweenness value this is because 1 out-degree of node is 12
Highest level in the whole network.This shows that the depended on amount vector of the node propagation of disturbance is more, with downstream node business association
By force, processing tasks flow through a possibility that node, and it is very high with the totality of non-neighbor node connect level.So under identical disturbance
Its influence power is larger, can by the propagation of disturbance into network most nodes.Although and node 20 enter intensity compared to node 1 compared with
It greatly, is 178, but its out-degree is only 1, downstream node only has node 8, and gathering coefficient is only 0.002, therefore, it is possible to judge that should
Node belongs to fringe node in a network, and the influence power of generation is smaller, therefore causes network failure small scale.
From production angle analysis: node 1 is fiber cutter, and there are more associations to manufacture node propagation disturbance, therefore
Influence is wider, and chain propagation is larger.Slotting machine is limited because association manufacture node is few by disturbance propagation path, chain propagation rule
Mould is smaller；Production analysis result is consistent with simulation analysis result.
Influence of the intensity to the production chain failure process of network out:
Identical disturbance, the network failure scales that intensity node causes different out are emulated, out the higher section of intensity value
Point more easily takes the lead in that network large area is caused to fail；Node under identical intensity out, there is also not on causing network failure scale
With the difference of degree, but as ascendant trend is integrally presented in the increase of intensity out, network failure scale, this meets universal law,
But there is obvious exception in its interior joint 29.As shown in Table 2, it is 66 that node 29, which goes out intensity, gathers that coefficient is low, pass in a network
Connection degree is low, and downstream node only has node 50, is limited by propagation path, and communication process is relatively slow.Therefore, even if to such
Node applies very big disturbance, and the entire network that produces also can finally tend towards stability in communication process, caused chain propagation rule
Mould is limited.
Can be seen that propagating source node diagnostic attribute difference from above-mentioned simulation result can generate difference to propagation of disturbance effect
Influence, and show different propagation laws.Wherein, when potential key node is as source node is propagated, failure procedure is more acute
It is strong, the chain effect for involving the whole network may be caused.Case verification also further illustrates proposed method and causes network for identification
The validity of the potential key node of chain effect.
The present invention is established based on the chain propagation algorithm for improving coupled map lattice systems, explores production system in bumpy weather
The chain propagation characteristic of lower formation carries out validation verification to above-mentioned theory and method in conjunction with agricultural machinery production practices, obtains as follows
Conclusion:
1) generation disturbed can break original equilibrium state of manufacturing process and lose a certain number of manufacture nodes
Effect, strength of turbulence and chain propagation scale are positively correlated trend.
2) angle value and intensity for propagating source node are in non-linear effects to the chain communication process of production system and scale.
3) propagating source node diagnostic attribute difference can generate different influences to propagation of disturbance effect, and show different
Propagation law, wherein when potential key node is as source node is propagated, failure procedure is more violent, may cause and involve the whole network
Chain effect, therefore, how to effectively control the propagation effect of disturbance bring risk in a network, this is also follow-up study work
Emphasis.
Detailed description of the invention
Fig. 1: for the polymorphic differentiation relational graph of production system of the invention.
Fig. 2: for chain propagation algorithm flow chart of the invention.
Fig. 3: production network topology structure figure of the invention.
Fig. 4: node failure I-R figure.
Fig. 5: degree-scale figure.
Fig. 6: in_strength-scale figure.
Fig. 7: out_strength-scale figure.
Specific embodiment
As shown in Figure 1, 2, 3, network chain effect research method is produced under a kind of bumpy weather based on transmission dynamics,
Chain propagation scale algorithm steps are as follows:
It is as follows to construct propagation model:
Step 1: using the manufacture node in production system as node, the process relation between production activity is to connect side, processing
Time is side right wij, determine adjacency matrix A, constitute the production system network under complex network visual angle, generate network；
Step 2: the parameter in propagation model is determined according to the production network of foundation, the coupling parameter of task inflow and outflow,
ε1、ε2, manufacture Node B threshold Yi, calculate node parameter
Step 3: the linear theory of 5W2H analytic method combined process is used, it is strong to disturbing according to delay time at stop and economic loss
Degree carry out assessment judgement, by disturbing influence degree be divided into slight perturbations, moderate disturbance, severe disturbance and four kinds of serious disturbance, it is right
I, II, III, IV 4 grade is answered, node different in network is selected as source node is propagated and is activated；
Step 4: all manufacture node states being updated by the parameter setting in propagation model formula (2) and step 2, if xi
(t+1)≥Yi, then manufacture node and be activated and switch to propagate state, while renewal time step-length；
Step 4.1: initialization node state value sums to the side right of node i and precedence activities, obtains the original negative of node
It carries L (i), the working ability of node is that the build-in attribute of node is Capality (i), and definition node original state is
Step 4.2: in addition t0Moment adds the load Δ L (i) of node, recalculates the state of node
The downstream adjacent node for finding node i, obtains downstream node set U (i), obtain node set interior joint with thereon
Swim the connection weight w of nodeji+With enter intensityAnd the connection weight w with node downstreami+jIntensity outAnd handle will
Process is saved in interim set Temp downstream, uses formula
It calculates
In t0The state value of+1 moment downstream node.When the node traverses of set U (i) are completed, using Temp as the influence of next stage
Node set, until whole network state reaches stable state；
Step 4.3: node being repaired, according to the working ability Capality (i) of node itself and all nodes
The maximum CapalityMax (i) of middle working ability ratio as the probability repaired, the working ability of node is improved 0
Random number between~20%；
Step 4.4: the quantity of failure node in each moment network is counted,
Step 4.5: choosing the angle value of node, intensity, betweenness and gather four network attributes of coefficient, in same disturbance situation
Under, it is arranged by the descending of each network attribute, identical disturbance successively is added to node, analysis heterogeneous networks characteristic is for the scale that fails
Influence, under different disturbances, analyzing the failure scale of network is influenced by disturbing；
Step 5: entire step 4 is recycled, until the state value of all nodes is respectively less than Node B threshold, i.e. x in networki(t+
1)<Yi, i=1,2 ..., n, the chain propagation effect of production system caused by disturbance factor terminates；
Wherein: xi(t+1) production status of the manufacture node i in moment t+1, x are indicatedi(t) indicate manufacture node i in moment t
Production status, ε1∈ (0,1) is the stiffness of coupling that task flows into, ε2∈ (0,1), for the stiffness of coupling of task outflow, ε1+ε2<
1, function f (x) indicate the individual node own dynamics behavior in network, and physical significance may be characterized as resource node production energy
The Evolution of force constraint selects Logistic mapping: f (x)=4x (1-x) herein, as 0≤x≤1,0≤f (x)≤1,
N1For the sum of manufacture node in-degrees all in network, N2For the sum of the out-degree of manufacture nodes all in network, A=a (i, j)N×NFor
Adjacency matrix describes the coupled relation that node is manufactured in network, if having Bian Xianglian, a between manufacture node i and manufacture node jij
=1, otherwise aij=aji=0, do not allow to manufacture node and is connected with itself, wijNeighbours, which are directed toward, for manufacture node i manufactures node j's
Side right,To enter intensity, indicate manufacture node i enters the sum of side right, To go out intensity, manufacture node i is indicated
The sum of side right out,Absolute value guarantees to manufacture in formula node production status not less than 0.
In the present invention, manufacture node i is disturbed in t moment, and i is activated at the t+1 moment, manufactures the association of node i
Manufacture node, which is manufactured node i state by the t+1 moment, to be influenced, and the state of these manufacture nodes involved is based on formula (2)
It calculates, state value is if it is greater than or equal to threshold value Yi, will lead to more manufacture node failures, such cycle calculations workshop is in bumpy weather
Under sprawling and chain effect, the polymorphism of production system is described, four kinds of states of production system under bumpy weather are portrayed
Evolution process and disturbing influence under manufacture the state of node and constitute, disclose disturbance factor propagation row in process of production
For.
Embodiment:
(1) data preparation
1) by taking certain agricultural machinery product enterprise C workshop data as an example, the chain effect that may cause under bumpy weather is ground
Study carefully.The workshop shares 36 manufacturing cells, mainly include fiber cutter, 3cm (6cm) cutter, 6cm bending machine, press machine,
Numerical control press, common punch press etc., are expressed as Ri, i={ 1,2 ..., 36 }.When table 1 is the processing and manufacturing unit of each workpiece, processing
Between and time.
Table 1 work time piece period t/min and time r/min
3) actual job-shop scheduling problem is solved, is obtained using makespan as the Optimized Operation side of target
Case: include machine sum m, workpiece sum n, the process letter processed on s-th of order of k-th of machine in initial schedule scheme
Breath: operation number, workpiece number, process time started, process end time；The process of each workpiece is made of m procedure.
(2) production system network model is established
According to Optimized Operation scheme, process is obtainedManufacturing recourses m and machine dimension successor activities machine dimensionWith
The successor activities of workpiece dimensionProcess is mapped to manufacturing recourses, occupies manufacture node j to manufacture the task of node i output
Time as the weight between node, to manufacture threshold value of the working ability as node of node, being formed with manufacturing cell is section
Point, process constraint are the production system network on even side.
(3) network parameter of network model is obtained
According to established network model, the network attribute of statistics manufacture node, node degree indicates node i associated nodes
Quantity,N is number of nodes, aijIndicate that node i, with the presence or absence of being associated with, exists for a with jij=1, otherwise aij
=0；Node strengthwijFor node i and j side right.Node betweenness indicates any two node in network
Number of the shortest path Jing Guo the node i,δtShow all shortest path numbers from vertex s to vertex t, and
δst(r) it indicates from vertex s to vertex t by all shortest path numbers of vertex r.C (i)=2Ei/(ki(ki-1)),kiIt indicates
With the associated number of nodes of node i, EiIndicate that there are associated quantity for node i associated nodes.
Table 2 produces network part node topological property statistics
(4) the chain communication process analysis of network is produced
The CML established the production chain propagation model of network is subjected to programming using python, from propagation of disturbance source
Node type, initial disturbance intensity analyze propagation of disturbance effect.
1) influence of the strength of turbulence size to the chain propagation of network
Consider network topology structure, side right and industry characteristics tripartite's region feature, chooses angle value maximum node, enters maximum intensity
Node goes out maximum intensity node and random node as propagation of disturbance source, while level of disruption is arranged and disturbs from slight perturbations, moderate
Dynamic, severe is disturbed to serious disturbance and is successively increased, and analyzes its chain propagation law.
When disturb initial value it is smaller when, propagate source node itself repair power can control disturbance propagate, fail scale domination
Below 10%.As disturbed value increases, propagates the chain propagation scale that source node causes for four kinds and show a increasing trend.
When maximum intensity node is as source node is propagated out, causes the failure of network size to be significantly larger than other nodes and make
To propagate destruction caused by source node.Speed is spread in the chain reaction that node (fiber cutter) failure of maximum intensity causes out
Degree is most fast, it is largest to involve, and followed by enters the node (machining center) of maximum intensity.Different type is propagated caused by source node
Failure scale and manufacture node itself capacity constraints and with the close phase of characteristic attributes such as the differentiation of neighbor node business association
It closes.Therefore in production activity, the big node of intensity is especially gone out to key node should put into more guarantees, reclamation activities, with
Containment disturbance is a wide range of to propagate, and causes heavy losses.
2) influence of the characteristic attribute of source node to the chain propagation of network is propagated
In order to which the discussion being more clear manufactures node unique characteristics attribute to the scale and process of the production chain failure of network
Influence, be provided with the experiment of three groups of simulation calculations, respectively analysis node angle value, enter intensity, go out intensity and the production chain mistake of network
The correlation of effect scale.
Influence of the angle value to the production chain failure process of network:
When identical disturbed value acts on different angle value nodes, the identical node of most of angle value causes the deviation of failure scale
± 10% or so, and as angle value increases, the chain propagation scale of initiation shows a increasing trend, and the big propagation source node of angle value is more
It is also easy to produce the propagation of the large area of disturbance, production system is caused to be paralysed, meets universal law.But its interior joint 2 and node 29 are
Apparent abnormal point.Analyze from network perspective: node 2 is in entirely production network, although its in-degree, the local features category such as out-degree
Property does not protrude, but has significant global property index, i.e., highest to gather coefficient.This shows node 2 at whole network
In criticality it is very high, and node is more downstream, and business association is strong, which changes the shadow to its associated nodes
It is big to ring power, therefore node 2 is big as caused network chain propagation scale when propagating source node, the very Gao Shui in whole network
It is flat.And for node 29, although its angle value is in by-level in overall network, its global convergence factor, betweenness are located
In reduced levels, and in entirely production network, subsequent association node only has node 50, and business association is weak, the influence of disturbance
Range is small, therefore node 29 is used as caused chain propagation small scale when propagating source node.
From production angle analysis: node 2 is machining center, and node 29 is boring machine.In actual production, since machining center has
Have preferable manufacture flexible, therefore the equipment be in a plurality of process route and flows through intertwined point in production, can and other equipment production
Raw very strong business association, and its occupancy is up to 93%, show the equipment resist disturbance buffer capacity it is weaker.Therefore,
Once such resource is disturbed, it can greatly may cause the large area sprawling of disturbance and propagate, reduce the production of entire producing line
Efficiency.And boring machine is limited by working ability and business association limits, under same disturbance, influence power is smaller, therefore caused
Chain propagation scale is limited.
Enter influence of the intensity to the production chain failure process of network:
Enter the network failure scale that intensity node causes to identical disturbance, difference to emulate, identical disturbed value acts on
When difference enters intensity node, with the increase for entering intensity value, the chain propagation scale of initiation shows a increasing trend, that is, undertakes task amount
After higher node is disturbed, the chain propagation scale of initiation also can be higher, this meets common cognition rule, but its interior joint
1 and node 20 be obvious abnormal point, it is only 80 that node 1, which enters intensity, but causes 85% failure scale, and node 20 enters intensity and is
178, cause the failure scale of network but very small, only 23.5%.
It is analyzed from network perspective: as can be seen from Table 2, intensity is 230 out, node betweenness value this is because 1 out-degree of node is 12
Highest level in the whole network.This shows that the depended on amount vector of the node propagation of disturbance is more, with downstream node business association
By force, processing tasks flow through a possibility that node, and it is very high with the totality of non-neighbor node connect level.So under identical disturbance
Its influence power is larger, can by the propagation of disturbance into network most nodes.Although and node 20 enter intensity compared to node 1 compared with
It greatly, is 178, but its out-degree is only 1, downstream node only has node 8, and gathering coefficient is only 0.002, therefore, it is possible to judge that should
Node belongs to fringe node in a network, and the influence power of generation is smaller, therefore causes network failure small scale.
From production angle analysis: node 1 is fiber cutter, and there are more associations to manufacture node propagation disturbance, therefore
Influence is wider, and chain propagation is larger.Slotting machine is limited because association manufacture node is few by disturbance propagation path, chain propagation rule
Mould is smaller；Production analysis result is consistent with simulation analysis result.
Influence of the intensity to the production chain failure process of network out:
Identical disturbance, the network failure scales that intensity node causes different out are emulated, out the higher section of intensity value
Point more easily takes the lead in that network large area is caused to fail；Node under identical intensity out, there is also not on causing network failure scale
With the difference of degree, but as ascendant trend is integrally presented in the increase of intensity out, network failure scale, this meets universal law,
But there is obvious exception in its interior joint 29.As shown in Table 2, it is 66 that node 29, which goes out intensity, gathers that coefficient is low, pass in a network
Connection degree is low, and downstream node only has node 50, is limited by propagation path, and communication process is relatively slow.Therefore, even if to such
Node applies very big disturbance, and the entire network that produces also can finally tend towards stability in communication process, caused chain propagation rule
Mould is limited.
Can be seen that propagating source node diagnostic attribute difference from above-mentioned simulation result can generate difference to propagation of disturbance effect
Influence, and show different propagation laws.Wherein, when potential key node is as source node is propagated, failure procedure is more acute
It is strong, the chain effect for involving the whole network may be caused.Case verification also further illustrates proposed method and causes network for identification
The validity of the potential key node of chain effect.
The present invention is established based on the chain propagation algorithm for improving coupled map lattice systems, explores production system in bumpy weather
The chain propagation characteristic of lower formation carries out validation verification to above-mentioned theory and method in conjunction with agricultural machinery production practices, obtains as follows
Conclusion:
1) generation disturbed can break original equilibrium state of manufacturing process and lose a certain number of manufacture nodes
Effect, strength of turbulence and chain propagation scale are positively correlated trend.
2) angle value and intensity for propagating source node are in non-linear effects to the chain communication process of production system and scale.
3) propagating source node diagnostic attribute difference can generate different influences to propagation of disturbance effect, and show different
Propagation law, wherein when potential key node is as source node is propagated, failure procedure is more violent, may cause and involve the whole network
Chain effect, therefore, how to effectively control the propagation effect of disturbance bring risk in a network, this is also follow-up study work
Emphasis.
The above technical characteristic constitutes the optimal embodiment of the present invention, and there is stronger adaptability and optimal implementation to imitate
Fruit can increase and decrease non-essential technical characteristic, according to actual needs to meet the needs of different situations.
Claims (1)
1. producing network chain effect research method under a kind of bumpy weather based on transmission dynamics, which is characterized in that chain
Propagation scale algorithm steps are as follows:
It is as follows to construct propagation model:
Step 1: using the manufacture node in production system as node, the process relation between production activity is to connect side, process time
For side right wij, determine adjacency matrix A, constitute the production system network under complex network visual angle, generate network；
Step 2: the parameter in propagation model, the coupling parameter of task inflow and outflow, ε are determined according to the production network of foundation1、ε2,
Manufacture Node B threshold Yi, calculate node parameter
Step 3: use the linear theory of 5W2H analytic method combined process, according to delay time at stop and economic loss to strength of turbulence into
Row assessment determines, disturbing influence degree is divided into slight perturbations, moderate disturbance, severe disturbance and four kinds of serious disturbance, corresponding I,
Tetra- grades of II, III, IV select node different in network as source node is propagated and are activated；
Step 4: all manufacture node states being updated by the parameter setting in propagation model formula (2) and step 2, if xi(t+1)
≥Yi, then manufacture node and be activated and switch to propagate state, while renewal time step-length；
Step 4.1: initialization node state value sums to the side right of node i and precedence activities, obtains the initial load L of node
(i), it is Capality (i) that the working ability of node, which is the build-in attribute of node, and definition node original state is
Step 4.2: in addition t0Moment adds the load Δ L (i) of node, recalculates the state of nodeThe downstream adjacent node for finding node i obtains downstream node set U (i), obtains and saves in node set
The connection weight w of point and its upstream nodeji+With enter intensityAnd the connection weight w with node downstreami+jIntensity outAnd that process will be saved in interim set Temp downstream, formula is used It calculates in t0When+1
Carve the state value of downstream node；When the node traverses of set U (i) are completed, using Temp as the influence node set of next stage,
Until whole network state reaches stable state；
Step 4.3: node being repaired, according in the working ability Capality (i) of node itself and all nodes
The ratio of the maximum CapalityMax (i) of working ability as the probability repaired, by the working ability of node improves 0~
Random number between 20%；
Step 4.4: the quantity of failure node in each moment network is counted,
Step 4.5: it chooses the angle value of node, intensity, betweenness and gathers four network attributes of coefficient, in same disturbance,
It is arranged by the descending of each network attribute, identical disturbance successively is added to node, analysis heterogeneous networks characteristic is for failure scale
It influences, under different disturbances, analyzing the failure scale of network is influenced by disturbing；
Step 5: entire step 4 is recycled, until the state value of all nodes is respectively less than Node B threshold, i.e. x in networki(t+1) <
Yi, i=1,2, L, n, the chain propagation effect of production system caused by disturbance factor terminates；
Wherein: xi(t+1) production status of the manufacture node i in moment t+1, x are indicatedi(t) indicate manufacture node i in the life of moment t
Occurrence state, ε1∈ (0,1) is the stiffness of coupling that task flows into, ε2∈ (0,1), for the stiffness of coupling of task outflow, ε1+ε2< 1,
Function f (x) indicates the individual node own dynamics behavior in network, and physical significance may be characterized as resource node production capacity
The Evolution of constraint selects Logistic mapping: f (x)=4x (1-x), as 0≤x≤1,0≤f (x)≤1, N herein1
For the sum of manufacture node in-degrees all in network, N2For the sum of the out-degree of manufacture nodes all in network, A=a (i, j)N×NFor neighbour
Matrix is connect, the coupled relation for manufacturing node in network is described, if having Bian Xianglian, a between manufacture node i and manufacture node jij=
1, otherwise aij=aji=0, do not allow to manufacture node and is connected with itself, wijThe side that neighbours manufacture node j is directed toward for manufacture node i
Power,To enter intensity, indicate manufacture node i enters the sum of side right, To go out intensity, indicate that manufacture node i goes out
The sum of side right,Absolute value guarantees to manufacture in formula node production status not less than 0.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811521510.8A CN109635439B (en) | 2018-12-12 | 2018-12-12 | Production network linkage effect research method based on propagation dynamics in disturbance environment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811521510.8A CN109635439B (en) | 2018-12-12 | 2018-12-12 | Production network linkage effect research method based on propagation dynamics in disturbance environment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109635439A true CN109635439A (en) | 2019-04-16 |
| CN109635439B CN109635439B (en) | 2023-05-05 |
Family
ID=66073369
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811521510.8A Active CN109635439B (en) | 2018-12-12 | 2018-12-12 | Production network linkage effect research method based on propagation dynamics in disturbance environment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109635439B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115099581A (en) * | 2022-05-25 | 2022-09-23 | 北京航天飞行控制中心 | Dynamic task planning method and device for satellite, electronic equipment and medium |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6466935B1 (en) * | 1999-06-04 | 2002-10-15 | International Business Machines Corporation | Applying relational database technology to process control in manufacturing processes |
| US20060235670A1 (en) * | 2003-08-01 | 2006-10-19 | Helios Biosciences | Method and system for selecting therapeutic targets using molecular interaction dynamic networks |
| CN101364107A (en) * | 2008-09-24 | 2009-02-11 | 永凯软件技术（上海）有限公司 | Dynamic dispatching method with prediction in complex large system environment |
| CN101587486A (en) * | 2009-06-23 | 2009-11-25 | 四川大学 | Method for organizing and applying design knowledge supporting product innovation |
| DE202009005097U1 (en) * | 2009-07-31 | 2009-12-24 | Soontai Tech Co., Ltd., Yongkang City | Multimedia network distributor |
| CN101916404A (en) * | 2010-08-06 | 2010-12-15 | 沈阳工业大学 | Multi-factory cooperative scheduling optimization method during equipment manufacturing |
| CN102129522A (en) * | 2011-03-17 | 2011-07-20 | 北京交通大学 | Method for quickly identifying and eliminating transportation capacity bottleneck of high-speed railway |
| CN103870642A (en) * | 2014-03-05 | 2014-06-18 | 湖州师范学院 | Complex supply network robust performance analyzing method based on topological structure |
| US20150112498A1 (en) * | 2013-09-23 | 2015-04-23 | Sureshchandra B Patel | Methods of Patel Loadflow Computation for Electrical Power System |
| CN108809946A (en) * | 2018-05-18 | 2018-11-13 | 哈尔滨工业大学深圳研究生院 | Vehicle-mounted net Malware spread restraint based on node influence power |
-
2018
- 2018-12-12 CN CN201811521510.8A patent/CN109635439B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6466935B1 (en) * | 1999-06-04 | 2002-10-15 | International Business Machines Corporation | Applying relational database technology to process control in manufacturing processes |
| US20060235670A1 (en) * | 2003-08-01 | 2006-10-19 | Helios Biosciences | Method and system for selecting therapeutic targets using molecular interaction dynamic networks |
| CN101364107A (en) * | 2008-09-24 | 2009-02-11 | 永凯软件技术（上海）有限公司 | Dynamic dispatching method with prediction in complex large system environment |
| CN101587486A (en) * | 2009-06-23 | 2009-11-25 | 四川大学 | Method for organizing and applying design knowledge supporting product innovation |
| DE202009005097U1 (en) * | 2009-07-31 | 2009-12-24 | Soontai Tech Co., Ltd., Yongkang City | Multimedia network distributor |
| CN101916404A (en) * | 2010-08-06 | 2010-12-15 | 沈阳工业大学 | Multi-factory cooperative scheduling optimization method during equipment manufacturing |
| CN102129522A (en) * | 2011-03-17 | 2011-07-20 | 北京交通大学 | Method for quickly identifying and eliminating transportation capacity bottleneck of high-speed railway |
| US20150112498A1 (en) * | 2013-09-23 | 2015-04-23 | Sureshchandra B Patel | Methods of Patel Loadflow Computation for Electrical Power System |
| CN103870642A (en) * | 2014-03-05 | 2014-06-18 | 湖州师范学院 | Complex supply network robust performance analyzing method based on topological structure |
| CN108809946A (en) * | 2018-05-18 | 2018-11-13 | 哈尔滨工业大学深圳研究生院 | Vehicle-mounted net Malware spread restraint based on node influence power |
Non-Patent Citations (1)
| Title |
|---|
| 丁剑 等: "基于贝叶斯网络的扰动后预想事故分析方法", 《电力系统自动化》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115099581A (en) * | 2022-05-25 | 2022-09-23 | 北京航天飞行控制中心 | Dynamic task planning method and device for satellite, electronic equipment and medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109635439B (en) | 2023-05-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103716202B (en) | A kind of intelligent maintenance strategy management method for power communication | |
| CN105975797B (en) | A kind of product initial failure root primordium recognition methods based on Fuzzy data processing | |
| CN106168799A (en) | A kind of method carrying out batteries of electric automobile predictive maintenance based on big data machine learning | |
| Sahu et al. | Benchmarking CNC machine tool using hybrid-fuzzy methodology: a multi-indices decision making (MCDM) approach | |
| CN105069521A (en) | Photovoltaic power plant output power prediction method based on weighted FCM clustering algorithm | |
| CN106355308B (en) | A method of wind power integration system core equipment is recognized based on decision tree | |
| CN107153845A (en) | A kind of isolated island detection method of the adaptive grid-connected photovoltaic system based on machine learning | |
| Kobbacy et al. | A survey of AI in operations management from 2005 to 2009 | |
| CN103729553A (en) | Classification control method for urban safety complex events on basis of Bayesian network learning | |
| CN107544450A (en) | Process industry network model construction method and system based on data | |
| Petković | Adaptive neuro-fuzzy optimization of the net present value and internal rate of return of a wind farm project under wake effect | |
| CN104506137A (en) | Equipment fault diagnosis method and apparatus | |
| CN114066122A (en) | Scheduling method based on multi-strategy water wave optimization algorithm | |
| CN109635439A (en) | Network chain effect research method is produced under bumpy weather based on transmission dynamics | |
| CN106408016A (en) | Distribution network power outage time automatic identification model construction method | |
| Rosli et al. | Predictive maintenance of air booster compressor (ABC) motor failure using artificial neural network trained by particle swarm optimization | |
| CN111160386B (en) | Data-driven power island detection method based on active reactive power disturbance | |
| Shahrabi et al. | Application of FMEA and AHP in Lean Maintenance | |
| CN117273428A (en) | Industrial chain toughness assessment method and system based on multi-chain coupling | |
| CN114118592B (en) | Smart power grids power consumption end short-term energy consumption prediction system | |
| WO2009082042A1 (en) | Nonlinear time series prediction method for weighted average defuzzification based on newfm | |
| CN114567562B (en) | Method for identifying key nodes of coupling network of power grid and communication network | |
| Chang et al. | A study of application of assessment modes in development of smart automatic production systems in Industry 4.0 | |
| CN114202205A (en) | HPLC product performance multidimensional evaluation method based on fuzzy analytic hierarchy process | |
| CN106033217B (en) | Dispatching system and work dispatching method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |