CN108563890A - A kind of quartz sand gas liquid inclusion three dimensional PIC code system - Google Patents
A kind of quartz sand gas liquid inclusion three dimensional PIC code system Download PDFInfo
- Publication number
- CN108563890A CN108563890A CN201810368507.0A CN201810368507A CN108563890A CN 108563890 A CN108563890 A CN 108563890A CN 201810368507 A CN201810368507 A CN 201810368507A CN 108563890 A CN108563890 A CN 108563890A
- Authority
- CN
- China
- Prior art keywords
- module
- model
- weights
- point
- triangle
- 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.)
- Pending
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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/02—CAD in a network environment, e.g. collaborative CAD or distributed simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention belongs to field of material technology, disclose a kind of quartz sand gas liquid inclusion three dimensional PIC code system, and the quartz sand gas liquid inclusion three dimensional PIC code system includes：Pressure detecting module, temperature detecting module, cubing module, main control module, coordinate curve drafting module, three-dimensional simulation module, data memory module, display module.The present invention accurately measures quartz sand gas liquid inclusion volume by cubing module, and measurement accuracy is significantly improved compared with previous technology；Observation quartz sand gas liquid inclusion and its stressing conditions more three-dimensional, comprehensive by three-dimensional simulation module simultaneously, is conducive to accurately analyze.
Description
Technical field
The invention belongs to field of material technology more particularly to a kind of quartz sand gas liquid inclusion three dimensional PIC code systems.
Background technology
Currently, the prior art commonly used in the trade is such：
Quartz sand is quartz through the broken quartz particles being process.Quartz is a kind of nonmetallic mineral, is one
Hard, wear-resisting, stable chemical performance the silicate mineral of kind.The color of quartz sand is milky or semi-transparent clear shape, Mohs
Hardness 7.Quartz sand is important industrial mineral raw material, and non-chemical dangerous material are widely used in glass, casting, ceramics and fire proofing wood
Material, smelting duriron, flux for metallurgy, metallurgy, building, chemical industry, plastics, rubber, abrasive material, the industry such as filtrate.However, existing to quartz
Sand gas liquid inclusion volume, performance measurement precision be not high, impact analysis data result；It is analyzed simultaneously by flat image, it cannot
Comprehensive, three-dimensional displaying analysis.
With the fast development of 3-D scanning hardware and the raising of three-dimensional modeling level, the fineness of threedimensional model is increasingly
Height, at the same time, the data volume of model also become very huge.It is on the one hand difficult to meet model real-time Transmission and processing in this way
It needs, large effect also on the other hand is caused to scene rendering efficiency, therefore must be carried out in such applications to model
Abbreviation.
It is also more and more for the demand for carrying out texture mapping on Irregular Boundary Surface in the every field that we live.Than
Such as, free form surface is commonly used in industrial modeling to carry out the configuration design of aircraft, steamer, automobile, is needed in laser marking machine
To carrying out laser marking on irregular mold cavity surface, usually reality simulation is carried out to virtual scene in 3d gaming,
Wallpaper or furniture effect are simulated in advance in intelligently decorateing, preview etc..Therefore mapping technology industrial design,
The fields such as the design of 3D scene of game, video display cartoon making and virtual emulation have a very wide range of applications.Irregular Boundary Surface
The research of texture mapping has certain actual application value.
The texture mapping of triangle mesh curved surface refers to that 2 d texture image is mapped to the curved surface of three-dimensional triangulation grid expression
On, that is, find a relationship mapped one by one of 2 d texture plane and triangle mesh curved surface.By this relationship, by two-dimentional line
The face of patting fits on three-dimensional Irregular Boundary Surface, it is made to have more the sense of reality.
Nowadays, as geometry obtains equipment (such as：Depth camera, spatial digitizer etc.) development, triangle gridding has become
To indicate a kind of main stream approach of body surface.For set curved surface, triangle gridding quantity is more, the curved surface more light of expression
It is sliding.And when directly doing texture mapping to the curved surface of quantity numerous triangle griddings expression, demand goes out each triangle gridding vertex
Texture coordinate can reduce the time performance of method in this way.
In conclusion problem of the existing technology is：
It is existing not high to quartz sand gas liquid inclusion volume, performance measurement precision, impact analysis data result；Pass through simultaneously
Flat image is analyzed, and displaying that cannot be comprehensive, three-dimensional is analyzed.
It is existing directly to do texture mapping method with the curved surface that triangle gridding indicates, it needs to find out out each triangle gridding vertex
Texture coordinate so that time performance is low, affects the practicability of analysis.
The prior art predicts that accuracy is poor, can not adapt to various change caused by different conditions to the prediction of volume.
Invention content
In view of the problems of the existing technology, the present invention provides a kind of quartz sand gas liquid inclusion three dimensional PIC code systems
System.
The invention is realized in this way a kind of quartz sand gas liquid inclusion three dimensional PIC code system, including：
Cubing module, connect with main control module, the volume data for measuring gas liquid inclusion；Cubing module
Measurement method include：
Relevant database data, sensing data and controller data that cubing module embeds are passed through into Sqoop
It is uploaded to distributed file system HDFS, and is stored into NoSQL databases；Using MapReduce Computational frames to relationship type
Database data, sensing data and controller data carry out mining analysis, and NoSQL databases are written in the data analyzed,
And it is shown by Web；
Apriori association rule minings are used under MapReduce frames, are specifically included：
The set L of frequent 1 item collection is obtained using MapReduce computation module1, generate the set C of candidate's k item collectionsk(k≥
2)；
In Map function processing stages, Transaction Information of each Map task computations handled by it is concentrated in each transaction journal
Included in CkIn Item Sets occurrence number, for each Map tasks, if some item collection of candidate's k item collections (includes
K project) it appears in a transaction journal, then Map functions are generated and are exported<Some item collection, 1>Key-value pair is to Combiner
Function gives Reduce functions after being handled by Combiner functions；
In Reduce function processing stages, Reduce functions add up CkIn Item Sets occurrence number, obtain all items
The support frequency of mesh collection, all minimum Item Sets composition frequent item set L for supporting frequency for supporting frequency >=settingkSet,
If the maximum iterations of k ＜ and not be sky, execute k++, be transferred to step Map function processing stages；Otherwise, terminate fortune
Row；
Main control module, with pressure detecting module, temperature detecting module, cubing module, coordinate curve drafting module, three
Analog module, data memory module, display module connection are tieed up, for dispatching modules normal work；
Coordinate curve drafting module, connect with main control module, for being depicted as sitting by mathematical software by the data of detection
Mark song line；
Three-dimensional simulation module, connect with main control module, for gas liquid inclusion to be fabricated to threedimensional model into Mobile state point
Analysis；It specifically includes：
Step 1: simplifying to triangle mesh model, to be simplified triangle mesh model；
Step 2: using the obtained simplified triangle mesh model of step 1, spring-mass model is established, is simplified
Final position of the triangle gridding by spring-mass model parameterization after rectangular domain；
Step 3: using the process for being simplified triangle mesh model, calculates and be deleted point in parameterizing rectangular domain
Initial position thus obtains position of all vertex in parameterizing rectangular domain in initial model；
Step 4: establishing spring-mass model to parametrization initial model, it is motionless to fix its border vertices, cause in point
Iteration obtains the parametrization coordinate of spring-mass Modifying model until interior position reaches balance；
Step 5:The parametrization coordinate of the spring-mass Modifying model obtained to step 4 is normalized, as
The texture coordinate of grid vertex carries out texture mapping.
Further, the measurement method of cubing module further includes：
Generate neural network model BP initial weights：By random number of the weight initialization between [a, b], wherein a, b are
Meet the integer of following equation：
Wherein H is network node in hidden layer；
The weights and threshold value of neural network model BP are dynamically refined, dynamic neural network model DBP is obtained, are generated
The weights and threshold value of dynamic neural network model DBP；
With adaptive immune genetic AIGA algorithm optimizations dynamic neural network model DBP, prediction model AIGA- is obtained
DBP calculates volume predictions value according to prediction model AIGA-DBP；
The weights and threshold value of neural network model BP are dynamically refined, dynamic neural network model DBP is obtained, are generated
The weights and threshold value of dynamic neural network model DBP specifically include following steps：
Adjust the weight w between neural network model BP hidden layers and output layerkj；
Adjust wkjPurpose be desirable to the new output o of output node j* pjO is exported than currentlypjCloser to desired value tpj, fixed
Justice：
Wherein α represents the degree of approach, is remained unchanged in each cycle of training, and becomes smaller with the adjustment of node in hidden layer H,
Do not consider threshold value, then has：
Wherein wk jAnd w* k jRespectively update front and back weights, ypkIt is exported for hidden layer, △ wkjFor wkjKnots modification；
Obtain △ wkjSolution equation：
Wherein,
Equation is solved according to least square and error principleObtain △ wkjApproximate solution：
It is connected to the hidden layer node k of output node j to each, calculates the weights variation △ w between k and jkj, update
Weights simultaneously calculate error of sum square E, then in one optimal k of k ∈ [1, H] interval selection so that E is minimum；
Adjust the weights v between neural network model BP input layers and hidden layerik；
Adjust vikPurpose be once neural network algorithm is absorbed in local minimum point, it is minimum that modification weights can jump out this
Point judges that condition that neural network algorithm is absorbed in local minimum point is the change rate △ E=0 of error E, and E>0；
Do not consider that threshold value, the change of the weights of hidden layer node k are solved by following equation：
Wherein δpj=f-1(ypk+Δypk)-f-1(ypk), M is natural number, then hidden layer exports ypkSolution formula is：
Wherein △ ypkFor ypkKnots modification, then have：
According to least square and error principle solution formulaThe matrix equation of structure, can
To calculate：
Aggregative formulaWithIt calculates implicit
The consecutive mean knots modification of weights between layer and output layer
Calculate the consecutive mean knots modification of weights between input layer and hidden layer
M takes the natural number between 10~20, the consecutive mean weights of neural network model BP is obtained, according to neural network mould
The consecutive mean weights of type BP obtain dynamic neural network model DBP.
It is planted with the power of adaptive immune genetic AIGA algorithm optimization dynamic neural network models DBP and threshold value, acquisition is pre-
It surveys model AIGA-DBP and specifically includes following steps：
Antigen recognizing, the i.e. identification of mean square error are carried out first, generate initial antibody (dynamic neural network model DBP
Weights and threshold value) after, carry out antibody coding, then calculating antibody fitness and antibody concentration are carried out according to the result of calculating
Adaptive antibody selection operation based on concentration obtains suitable weights and threshold value if this operation meets condition, and by result
Output；If this operation is unsatisfactory for condition, need to carry out adaptive crossover operation and TSP question operation, and carry out weights and threshold
The update of value, then antibody fitness and antibody concentration are recalculated, weights and threshold value until obtaining the condition that meets, last basis
Meet the weights and threshold value of condition, obtains prediction model AIGA-DBP.
Further, described Step 1: simplify to triangle mesh model, to be simplified triangle mesh model, packet
Include following steps：
(1) triangle grid model initialization is carried out, every directed edge in triangle gridding is traversed, if this side is only by one three
Angular to include, then two vertex for illustrating to connect this side are boundary point；
(2) input needs the number of vertices x retained and angle threshold value θ；
(3) judge whether existing number of vertices is more than vertex threshold value x；If so, carrying out (4)；If it is not, then carrying out (9)；
(4) length of side on each side, and the normal vector of each triangle are found out respectively, are calculated using following formula oriented
Collapse cost cost (the e on sideij)；And find out the directed edge e of collapse cost minimumuv, then carried out (5)；
cost(eij)=‖ eij‖·w(eij) (2)
Wherein w (eij) indicate directed edge eijLocal curvature, TiIt indicates to include vertex viTriangle sets, TeijIt indicates
Including side eijTriangle sets, f.normal refers to the normal vector in the faces f, and n.normal refers to the normal vector on the faces n；It is indicated
Include directed edge eijTri patch and include vertex viTri patch normal direction maximum value；
(5) judge whether u, v are boundary point.If so, carrying out (6)；If it is not, then carrying out (7)；
(6) judge whether the rim angle at boundary point u is more than threshold θ；If so, carrying out (8)；If it is not, then carrying out (3)；
(7) it is boundary point to judge whether u, v meet u, and v is not boundary point；If satisfied, then carrying out (3)；If not satisfied, then
It carries out (8)；
(8) vertex u is deleted, the frontier juncture being connect with vertex u is connected into vertex v.And it is impacted to mark 2 points of u, v
Point carries out (3)；
(9) terminate；
It is described Step 2: using the obtained simplified triangle mesh model of step 1, establish spring-mass model, obtain
Simplify final position of the triangle gridding by spring-mass model parameterization after rectangular domain, includes the following steps：
(a) triangle grid model is read in, point mass m, device of spring stiffness coefficient k, time step Δ t, iteration domain are initialized
Value e；
(b) triangle grid data pretreatment is carried out, triangle grid model boundary point sequence is obtained, to there is the adjoining on each vertex
Point sequence；
(c) perspective plane is selected according to formula (3), (4), triangle grid model is projected, established two dimension on perspective plane and sit
Mark system calculates the two-dimensional coordinate (x after the projection of triangle gridding vertexi,yi)；
a(x-x0)+b(y-y0)+c(z-z0)=0 (3)
Wherein, (a, b, c) is the normal vector on perspective plane, x0,y0,z0It is a fixed point, (a in spacei,bi,ci) it is i-th three
Angular normal vector, n are the numbers of Irregular Boundary Surface on piece triangle；
(d) four are chosen according to the structure feature of the triangle gridding after projection in four vertex of Selecting All Parameters rectangular area
A boundary point is fixed on the vertex of given rectangular domain；
(e) due to the movement of point, the stressing conditions for causing on-fixed point change；According to formula (5), on-fixed point is calculated
Suffered resultant force f (Pi)；If for arbitrary point Pi, suffered resultant force is respectively less than given threshold value, i.e. f (Pi)≤e, then terminate current iteration,
Turn (i) when terminating iteration for the first time, turns (j) when terminating iteration for the second time；Otherwise, continue；
Wherein, M is the mass matrix of spring-mass model system, and q is particle coordinate, and K is stiffness matrix；
(f) according to formula (6), by resultant force f (Pi) calculate the acceleration each put
Wherein, m is point mass,For in t moment particle PiAcceleration；
(g) by accelerationAnd given time step Δ t, calculate the displacement each put according to formula (7), (8)；
Wherein,WithThe respectively particle P in t and t+ time ΔtsiSpeed,WithRespectively in t and
T+ time Δt particles PiPosition；
(h) by the displacement in (g), each point is calculated after the Δ t times, the current location of pointTurn (e)；
(i) by loose boundary point according in space lattice between boundary point and boundary point broken line ratio, be fixed to
Four sides of given rectangular domain；It is further continued for carrying out (i)；
(j) parameterized results are obtained.
Further, the quartz sand gas liquid inclusion three dimensional PIC code system further includes：
Pressure detecting module, connect with main control module, the number pressure for detecting gas liquid inclusion by pressure sensor
According to；
Temperature detecting module is connect with main control module, the temperature number for detecting gas liquid inclusion by temperature sensor
According to；
Data memory module is connect with main control module, is stored for that will analyze data；
Display module is connect with main control module, for showing detection data and threedimensional model.
Further, the cubing module measurement method is as follows：
First, it is scanned using laser confocal microscope, obtains quartz sand gas liquid inclusion three-dimensional image；
Secondly, the boundary of three-dimensional coordinate calibration three-dimensional image liquid phase is utilized under transmitted light and fluorescence；
Then, according to the liquid phase boundary of calibration, the volume of quartz sand gas liquid inclusion liquid phase is calculated；
Then, the gas phase diameter d in quartz sand gas liquid inclusion is accurately obtained using laser co-focusing micro-fluorescence spectrum
And gas phase radius R；
Finally, the volume V of gas phase in inclusion enclave is calculated according to formula V=(4/3) π R3.
Further, the three-dimensional simulation module includes modelling module, parameter setting module, deduction module；
Modelling module, for gas liquid inclusion to be fabricated to threedimensional model by 3d prototype softwares；
Parameter setting module, for gas liquid inclusion to be arranged by force parameter；
Deduction module, for the dynamic state of parameters deduction gas liquid inclusion stressing conditions by being arranged.
Advantages of the present invention and good effect are：
The present invention accurately measures quartz sand gas liquid inclusion volume, measurement accuracy and previous skill by cubing module
Art is compared and is significantly improved；Observation quartz sand gas liquid inclusion more three-dimensional, comprehensive by three-dimensional simulation module simultaneously
And its stressing conditions, be conducive to accurately analyze.
The present invention need not find out the texture coordinate on each triangular network lattice vertex, can directly do texture mapping, can be with
The largely time complexity of reduction method, improves time performance, makes parametrization by the application of spring-mass model
Triangle gridding afterwards maintains the geometry of original triangle gridding, and the application of spring-mass model can keep triangle gridding
Original geometry is characteristically uniformly distributed, and will not cause texture convergence and texture moderate finite deformation phenomenon, therefore, texture mapping result
In deformation texture it is smaller.
Then the present invention excavates total complete time influence factor with association rule algorithm, and build neural network model
BP is dynamically refined the weights and threshold value of neural network model BP, to obtain dynamic neural network model DBP, then transports
Prediction model AIGA-DBP is obtained with adaptive immune genetic AIGA algorithm optimization dynamic neural network models DBP, is finally used
Prediction model AIGA-DBP calculates volume predictions value, and the dynamic neural network model DBP in the present invention can adapt to because of difference
Various change caused by state.
Big data analysis technology is used in the present invention so that the excavation of cubing influence factor is highly efficient and accurate
Really, influence factor considers more comprehensive, effectively improves the accuracy of prediction.
Description of the drawings
Fig. 1 is quartz sand gas liquid inclusion three dimensional PIC code system structure diagram provided in an embodiment of the present invention.
In figure：1, pressure detecting module；2, temperature detecting module；3, cubing module；4, main control module；5, coordinate is bent
Line drafting module；6, three-dimensional simulation module；7, data memory module；8, display module.
Fig. 2 is the rough schematic view of boundary edge provided in an embodiment of the present invention.
Fig. 3 is spring-mass model schematic provided in an embodiment of the present invention.
Fig. 4 is simplification provided in an embodiment of the present invention determines schematic diagram by the parametrization initial position for deleting a little.
Specific implementation mode
In order to further understand the content, features and effects of the present invention, the following examples are hereby given, and coordinate attached drawing
Detailed description are as follows.
As shown in Figure 1, quartz sand gas liquid inclusion three dimensional PIC code system provided in an embodiment of the present invention includes：Pressure
Detection module 1, temperature detecting module 2, cubing module 3, main control module 4, coordinate curve drafting module 5, three-dimensional simulation mould
Block 6, data memory module 7, display module 8.
Pressure detecting module 1 is connect with main control module 4, the pressure for detecting gas liquid inclusion by pressure sensor
Data；
Temperature detecting module 2 is connect with main control module 4, the temperature for detecting gas liquid inclusion by temperature sensor
Data；
Cubing module 3 is connect with main control module 4, the volume data for measuring gas liquid inclusion；
Main control module 4 draws mould with pressure detecting module 1, temperature detecting module 2, cubing module 3, coordinate curve
Block 5, three-dimensional simulation module 6, data memory module 7, display module 8 connect, for dispatching modules normal work；
Coordinate curve drafting module 5 is connect with main control module 4, for being depicted as the data of detection by mathematical software
Coordinate curve；
Three-dimensional simulation module 6 is connect with main control module 4, for gas liquid inclusion to be fabricated to threedimensional model into Mobile state
Analysis；
Data memory module 7 is connect with main control module 4, is stored for that will analyze data；
Display module 8 is connect with main control module 4, for showing detection data and threedimensional model.
3 measurement method of cubing module provided by the invention is as follows：
First, it is scanned using laser confocal microscope, obtains quartz sand gas liquid inclusion three-dimensional image；
Secondly, the boundary of three-dimensional coordinate calibration three-dimensional image liquid phase is utilized under transmitted light and fluorescence；
Then, according to the liquid phase boundary of calibration, the volume of quartz sand gas liquid inclusion liquid phase is calculated；
Then, the gas phase diameter d in quartz sand gas liquid inclusion is accurately obtained using laser co-focusing micro-fluorescence spectrum
And gas phase radius R；
Finally, according to formula V=(4/3) π R3Calculate the volume V of gas phase in inclusion enclave.
Three-dimensional simulation module 6 provided by the invention includes modelling module, parameter setting module, deduction module；
Modelling module, for gas liquid inclusion to be fabricated to threedimensional model by 3d prototype softwares；
Parameter setting module, for gas liquid inclusion to be arranged by force parameter；
Deduction module, for the dynamic state of parameters deduction gas liquid inclusion stressing conditions by being arranged.
When the present invention works,
The pressure data of gas liquid inclusion is detected by pressure detecting module 1；Gas-liquid packet is detected by temperature detecting module 2
Wrap up in the temperature data of body；The volume data of gas liquid inclusion is measured by cubing module 3；It is bent that main control module 4 dispatches coordinate
Line drafting module 5 draws mathematical coordinates curve；Gas liquid inclusion is fabricated to threedimensional model into action by three-dimensional simulation module 6
State is analyzed；It is stored by data memory module 7 by data are analyzed；Detection data and three-dimensional mould are shown by display module 8
Type.
The measurement method of cubing module of the present invention includes：
Relevant database data, sensing data and controller data that cubing module embeds are passed through into Sqoop
It is uploaded to distributed file system HDFS, and is stored into NoSQL databases；Using MapReduce Computational frames to relationship type
Database data, sensing data and controller data carry out mining analysis, and NoSQL databases are written in the data analyzed,
And it is shown by Web；
Apriori association rule minings are used under MapReduce frames, are specifically included：
The set L of frequent 1 item collection is obtained using MapReduce computation module1, generate the set C of candidate's k item collectionsk(k≥
2)；
In Map function processing stages, Transaction Information of each Map task computations handled by it is concentrated in each transaction journal
Included in CkIn Item Sets occurrence number, for each Map tasks, if some item collection of candidate's k item collections (includes
K project) it appears in a transaction journal, then Map functions are generated and are exported<Some item collection, 1>Key-value pair is to Combiner
Function gives Reduce functions after being handled by Combiner functions；
In Reduce function processing stages, Reduce functions add up CkIn Item Sets occurrence number, obtain all items
The support frequency of mesh collection, all minimum Item Sets composition frequent item set L for supporting frequency for supporting frequency >=settingkSet,
If the maximum iterations of k ＜ and not be sky, execute k++, be transferred to step Map function processing stages；Otherwise, terminate fortune
Row；
Main control module, with pressure detecting module, temperature detecting module, cubing module, coordinate curve drafting module, three
Analog module, data memory module, display module connection are tieed up, for dispatching modules normal work；
Coordinate curve drafting module, connect with main control module, for being depicted as sitting by mathematical software by the data of detection
Mark song line；
Three-dimensional simulation module, connect with main control module, for gas liquid inclusion to be fabricated to threedimensional model into Mobile state point
Analysis；It specifically includes：
Step 1: simplifying to triangle mesh model, to be simplified triangle mesh model；
Step 2: using the obtained simplified triangle mesh model of step 1, spring-mass model is established, is simplified
Final position of the triangle gridding by spring-mass model parameterization after rectangular domain；
Step 3: using the process for being simplified triangle mesh model, calculates and be deleted point in parameterizing rectangular domain
Initial position thus obtains position of all vertex in parameterizing rectangular domain in initial model；
Step 4: establishing spring-mass model to parametrization initial model, it is motionless to fix its border vertices, cause in point
Iteration obtains the parametrization coordinate of spring-mass Modifying model until interior position reaches balance；
Step 5:The parametrization coordinate of the spring-mass Modifying model obtained to step 4 is normalized, as
The texture coordinate of grid vertex carries out texture mapping.
Further, the measurement method of cubing module further includes：
Generate neural network model BP initial weights：By random number of the weight initialization between [a, b], wherein a, b are
Meet the integer of following equation：
Wherein H is network node in hidden layer；
The weights and threshold value of neural network model BP are dynamically refined, dynamic neural network model DBP is obtained, are generated
The weights and threshold value of dynamic neural network model DBP；
With adaptive immune genetic AIGA algorithm optimizations dynamic neural network model DBP, prediction model AIGA- is obtained
DBP calculates volume predictions value according to prediction model AIGA-DBP；
The weights and threshold value of neural network model BP are dynamically refined, dynamic neural network model DBP is obtained, are generated
The weights and threshold value of dynamic neural network model DBP specifically include following steps：
Adjust the weight w between neural network model BP hidden layers and output layerkj；
Adjust wkjPurpose be desirable to the new output o of output node j* pjO is exported than currentlypjCloser to desired value tpj, fixed
Justice：
Wherein α represents the degree of approach, is remained unchanged in each cycle of training, and becomes smaller with the adjustment of node in hidden layer H,
Do not consider threshold value, then has：
Wherein wkjAnd w* kjRespectively update front and back weights, ypkIt is exported for hidden layer, △ wkjFor wkjKnots modification；
Obtain △ wkjSolution equation：
Wherein,
Equation is solved according to least square and error principleObtain △ wkjApproximate solution：
It is connected to the hidden layer node k of output node j to each, calculates the weights variation △ w between k and jkj, update
Weights simultaneously calculate error of sum square E, then in one optimal k of k ∈ [1, H] interval selection so that E is minimum；
Adjust the weights v between neural network model BP input layers and hidden layerik；
Adjust vikPurpose be once neural network algorithm is absorbed in local minimum point, it is minimum that modification weights can jump out this
Point judges that condition that neural network algorithm is absorbed in local minimum point is the change rate △ E=0 of error E, and E>0；
Do not consider that threshold value, the change of the weights of hidden layer node k are solved by following equation：
Wherein δpj=f-1(ypk+Δypk)-f-1(ypk), M is natural number, then hidden layer exports ypkSolution formula is：
Wherein △ ypkFor ypkKnots modification, then have：
According to least square and error principle solution formulaThe matrix equation of structure, can
To calculate：
Aggregative formulaWithCalculate hidden layer
The consecutive mean knots modification of weights between output layer
Calculate the consecutive mean knots modification of weights between input layer and hidden layer
M takes the natural number between 10~20, the consecutive mean weights of neural network model BP is obtained, according to neural network mould
The consecutive mean weights of type BP obtain dynamic neural network model DBP.
It is planted with the power of adaptive immune genetic AIGA algorithm optimization dynamic neural network models DBP and threshold value, acquisition is pre-
It surveys model AIGA-DBP and specifically includes following steps：
Antigen recognizing, the i.e. identification of mean square error are carried out first, generate initial antibody (dynamic neural network model DBP
Weights and threshold value) after, carry out antibody coding, then calculating antibody fitness and antibody concentration are carried out according to the result of calculating
Adaptive antibody selection operation based on concentration obtains suitable weights and threshold value if this operation meets condition, and by result
Output；If this operation is unsatisfactory for condition, need to carry out adaptive crossover operation and TSP question operation, and carry out weights and threshold
The update of value, then antibody fitness and antibody concentration are recalculated, weights and threshold value until obtaining the condition that meets, last basis
Meet the weights and threshold value of condition, obtains prediction model AIGA-DBP.
Further, described Step 1: simplify to triangle mesh model, to be simplified triangle mesh model, packet
Include following steps：
(1) triangle grid model initialization is carried out, every directed edge in triangle gridding is traversed, if this side is only by one three
Angular to include, then two vertex for illustrating to connect this side are boundary point；
(2) input needs the number of vertices x retained and angle threshold value θ；
(3) judge whether existing number of vertices is more than vertex threshold value x；If so, carrying out (4)；If it is not, then carrying out (9)；
(4) length of side on each side, and the normal vector of each triangle are found out respectively, are calculated using following formula oriented
Collapse cost cost (the e on sideij)；And find out the directed edge e of collapse cost minimumuv, then carried out (5)；
cost(eij)=‖ eij‖·w(eij) (2)
Wherein w (eij) indicate directed edge eijLocal curvature, TiIt indicates to include vertex viTriangle sets,Indicate packet
E containing sideijTriangle sets, f.normal refers to the normal vector in the faces f, and n.normal refers to the normal vector on the faces n；It is illustrated
Including directed edge eijTri patch and include vertex viTri patch normal direction maximum value；
(5) judge whether u, v are boundary point.If so, carrying out (6)；If it is not, then carrying out (7)；
(6) judge whether the rim angle at boundary point u is more than threshold θ；If so, carrying out (8)；If it is not, then carrying out (3)；
(7) it is boundary point to judge whether u, v meet u, and v is not boundary point；If satisfied, then carrying out (3)；If not satisfied, then
It carries out (8)；
(8) vertex u is deleted, the frontier juncture being connect with vertex u is connected into vertex v.And it is impacted to mark 2 points of u, v
Point carries out (3)；
(9) terminate；
It is described Step 2: using the obtained simplified triangle mesh model of step 1, establish spring-mass model, obtain
Simplify final position of the triangle gridding by spring-mass model parameterization after rectangular domain, includes the following steps：
(a) triangle grid model is read in, point mass m, device of spring stiffness coefficient k, time step Δ t, iteration domain are initialized
Value e；
(b) triangle grid data pretreatment is carried out, triangle grid model boundary point sequence is obtained, to there is the adjoining on each vertex
Point sequence；
(c) perspective plane is selected according to formula (3), (4), triangle grid model is projected, established two dimension on perspective plane and sit
Mark system calculates the two-dimensional coordinate (x after the projection of triangle gridding vertexi,yi)；
a(x-x0)+b(y-y0)+c(z-z0)=0 (3)
Wherein, (a, b, c) is the normal vector on perspective plane, x0,y0,z0It is a fixed point, (a in spacei,bi,ci) it is i-th three
Angular normal vector, n are the numbers of Irregular Boundary Surface on piece triangle；
(d) four are chosen according to the structure feature of the triangle gridding after projection in four vertex of Selecting All Parameters rectangular area
A boundary point is fixed on the vertex of given rectangular domain；
(e) due to the movement of point, the stressing conditions for causing on-fixed point change；According to formula (5), on-fixed point is calculated
Suffered resultant force f (Pi)；If for arbitrary point Pi, suffered resultant force is respectively less than given threshold value, i.e. f (Pi)≤e, then terminate current iteration,
Turn (i) when terminating iteration for the first time, turns (j) when terminating iteration for the second time；Otherwise, continue；
Wherein, M is the mass matrix of spring-mass model system, and q is particle coordinate, and K is stiffness matrix；
(f) according to formula (6), by resultant force f (Pi) calculate the acceleration each put
Wherein, m is point mass,For in t moment particle PiAcceleration；
(g) by accelerationAnd given time step Δ t, calculate the displacement each put according to formula (7), (8)；
Wherein,WithThe respectively particle P in t and t+ time ΔtsiSpeed,WithRespectively in t and t+ Δs
T moment particle PiPosition；
(h) by the displacement in (g), each point is calculated after the Δ t times, the current location of pointTurn (e)；
(i) by loose boundary point according in space lattice between boundary point and boundary point broken line ratio, be fixed to
Four sides of given rectangular domain；It is further continued for carrying out (i)；
(j) parameterized results are obtained.
Fig. 2 is the rough schematic view of boundary edge provided in an embodiment of the present invention.
Fig. 3 is spring-mass model schematic provided in an embodiment of the present invention.
Fig. 4 is simplification provided in an embodiment of the present invention determines schematic diagram by the parametrization initial position for deleting a little.
With reference to concrete analysis, the invention will be further described.
Spring-mass model foundation
Spring-mass model is often used as simulating and controlling the deformation of triangle gridding in computer graphics.In reality
Spring-Mass Model is multi-purpose in the application of border makees simulation fabric movement [38] and complex plane expansion etc..Its basic principle is：By three
Point in the grid of angle regards particle as, while regarding spring as；While the original that the length in original triangle gridding corresponds to spring is long, side
Long variation causes particle movement corresponding to the stretching or compression of spring to which spring generates pulling force or pressure to particle.Bullet
Spring-Mass Model schematic diagram is as shown in figure 3, wherein Fig. 3 (a) indicates particle P in triangle griddingiAnd the initial bit on one field side
It sets, Fig. 3 (b) indicates particle PiStress balance situation after position changes.Define particle PiInstitute's stress is：
Wherein, C is the stiffness factor of spring, | PiPik| andRespectively PiPoint original state and after changing with its
K-th of abutment points PikThe distance between,For by particle Pi *It is directed towardDirection, j indicates particle PiThe number of abutment points.
It is constrained when to certain particles in spring-mass model, then passes through the fortune of spring-mass model system particle
It is dynamic, when reaching balance with joint efforts suffered by each particle, it will be able to obtain the equilbrium position of each particle in spring-mass model system.
And the movement of particle is described using lagrange equation of motion in entire spring-mass model, i.e.,：
Wherein, M, D and K are respectively the mass matrix, damping matrix and stiffness matrix of spring-mass model system, q
For particle coordinate.gqFor force vector in the conjunction of whole system, fqFor the bonding force vector of system.In spring-mass modeling three
When the grid deformation of angle, gqAnd fqIt is zero.In normal circumstances, we do not consider damping term generally, in this way, above-mentioned (3.4)
Formula can abbreviation be：
Wherein, Kq is elastic force.Above-mentioned (5) formula is usually solved with Euler's method.Entire spring-mass model system
Balance is made of the balance of each particle.When Δ t is smaller, particle PiThe motion process of t to t+ time Δts can be by upper
The simplification Lagrange's equation of (5) after simplification is stated to indicate：
Wherein, m is point mass,For in t moment particle PiAcceleration,WithRespectively in t and t+ Δ t
Particle P when quarteriSpeed,WithRespectively in t and t+ time Δt particles PiPosition.
Because by formula (6-1) it was found that the speed in moment t can generate certain influence to the speed of moment t+ Δ t,
It is this to influence the inertia that namely we often say.From initial motion to terminating to move, will once be used to every time Δt respectively
Property effect, this creates the terminal accumulative effects.When particle is when equilbrium position reaches maximum speed, if acceleration and speed
Direction on the contrary, particle still can continue to move according to the direction of original speed this moment, rather than at once towards the side of acceleration
To movement.When this happens, entire spring-mass model system just will appear larger oscillation.And this oscillation
In the presence of situation, spring-mass model system just needs more times to reach the balance of each particle, greatly reduces algorithm
Speed.In addition to this, due to not accounting for damping factor in the algorithm, this accumulative effect is it is also possible to cause to send out
It dissipates, more cannot get the balance of particle.
It is solved the above problems using the method for ignoring initial velocity of the propositions such as hair.Give tacit consent to each moment t, particle just
Speed is zero, can prevent the influence of the effect of inertia of particle in this way, then above-mentioned formula (6-1), (6-2) are adjusted to：
In addition, in this algorithm, point mass takes constant.
After the spring-mass system of triangle grid model is established as stated above, triangle gridding mould is artificially chosen
Four boundary points of type, and it is first fixed on to four vertex of some rectangle.After fixing four vertex, remaining particle
Suffered resultant force is not zero, i.e. spring-mass system disequilibrium.Then at this point, providing a threshold θ again, with time interval Δ t.
We haggle over resultant force suffered by each primary particle every time Δt.The resultant force suffered by each particle in this system is
When less than this threshold θ, this spring-mass system is just claimed to reach balance.At this point, further according in triangle gridding archetype
Remaining boundary is pressed length of side ratio and is fixed on given square boundary by boundary edge company commander's ratio.Above-mentioned change is carried out once again
Generation.When all particles reach balance, i.e., when entire spring-mass model system reaches balance, terminate iteration.At this moment particle
Position is exactly to simplify final position of the triangle gridding by spring-mass model parameterization after rectangular domain.
Parameterized results amendment
Vertex parameter initial position determines
Due in the parameterized procedure, after spring-mass model is established in the projection of original pessimistic concurrency control, four sides are specified
Four vertex that boundary vertex is stretched to specified rectangular domain are fixed, remaining point force unbalance is caused, thus given time step-length,
Remaining point is iterated, and specified rectangular domain is significantly larger than view field, the parameter that can be just more uniformly distributed in this way during this
Change result.Therefore vertex is larger to equilbrium position displacement by initial position, and because time step cannot be too big, too big system hair
It dissipates, so the time that iteration is spent is longer.It is spent in conclusion to reduce the algorithm time it is necessary to reduce particle to the greatest extent by initial
Displacement of the position to equilbrium position.
After model simplification, retain the parametric method that preceding Zhang Suoshu is passed through on vertex, obtains the two dimension in specified rectangular domain
Coordinate.Using this as the initial position for retaining vertex parameter in this algorithm.The initial position on vertex is deleted in archetype
It is the parametrization initial position center of one field adjacent vertex.Then an areas is vi1vi2…vim(m is adjacent vertex
Number) vertex viParametrization initial coordinate position (xi,yi) calculation formula such as formula (9) shown in：
Wherein, (xik,yik) (k=1...m) vertex v for beingiA field adjacent vertex vikParametrization initial position.
The determination for being deleted vertex parameter initial position is illustrated with, as shown in figure 4, wherein Fig. 4 (a)-Fig. 4
(d) it is that Triangular Mesh Simplification deletes point process figure, Fig. 4 (e)-Fig. 4 (f) is to be deleted a parametrization location determination process.
In Fig. 4, (a)-(d) is that Triangular Mesh Simplification deletes point process, vertex vi、vj、vkSuccessively by quilt during simplification
It deletes.(e)-(f) is to be deleted a parametrization location determination process.Rectangular domain intrinsic parameter is carried out to grid after simplification, if top
Point v1v2v3v4v5v6In the region surrounded, interior point obtains the 4th width figure in (a)-(d), the threedimensional model after simplifying by simplifying
By parameterized procedure, the first width figure in (e)-(f) is obtained, at this time vertex v1、v2、v3、v4、v5、v6Parametrization initial position
It has been obtained that, respectively (x1,y1)、(x2,y2)、(x3,y3)、(x4,y4)、(x5,y5)、(x6,y6).Record current state is finally deleted
Vertex vkAnd delete vkIt is the second width figure v in (e)-(f) beforekAn areas v1v3v4v5.So vkParameter it is initial
Position (xk,yk) just by its adjacent vertex v1、v3、v4、v5Parametrization initial position (x1,y1)、(x3,y3)、(x4,y4)、(x5,
y5) determined by formula：
vkParameterizing initial coordinate is deleted in back when the average value of previous field apex coordinate by it.vkGinseng
After numberization initial coordinate determines, v is foundkThe last one deleted vertex v before being deletedj, then as stated above, find out vjQuilt
An areas v before deletion1v2v3vk, at this point, one field vertex (including vk) parametrization initial coordinate be it is known that by
The above method finds out vjParametrization initial coordinate (xj,yj)；vjParametrization initial coordinate determine after, as stated above, can transport
With the v found outk、vjPoint parametrization initial coordinate finds out vjThe last one deleted vertex v before being deletediParametrization at the beginning of
Beginning position.
Spring-mass model foundation
After obtaining the parametrization initial position of all the points, vertex is connected by the topological structure of original triangle grid model
It connects, obtains initial parameter model.Regard vertex as particle on initial parameter model, Grid Edge is regarded as to the bullet of massless
Spring is grown with a length of spring original of original on side in initial three-dimensional triangle gridding, using the current length of side as the length at spring current time, is established
Spring-mass model.The movement of particle is still described using the simplification lagrange equation of motion for ignoring initial velocity.
After establishing spring-mass model, the boundary point of fixed model is motionless, given time step-length, allows after simplifying and retains
Vertex and deleted vertex are iterated, until all vertex stress balances, terminate iteration.The vertex warp retained after simplification
The process that iteration reaches balance is crossed, the amendment for parameterizing initial results to it is equivalent to.At this point, obtaining institute in original triangle gridding
There is the final rectangular domain intrinsic parameter coordinate on vertex.The parameterized results that this method can be worked well, but if
Need the time performance of algorithm more excellent, i.e., interior point iterations are less, then just needing interior number of reduction participation iteration.
Interior point is divided into affected point and unaffected point, affected point refer to during Triangular Mesh Simplification, this
The connected side of point, which has, to be folded, then this point is exactly affected point.Unaffected point refers to that this puts connected side
It is not folded.During simplification, deleted vertex must be affected point, retain vertex and have plenty of affected point, have
Be unaffected point.After establishing spring-mass model, fixed boundary point and unaffected point are motionless, allow affected point
It is iterated, until all vertex stress balances, terminate iteration.It is motionless by being put in fixed part, you can to improve efficiency of algorithm.
The above is only the preferred embodiments of the present invention, and is not intended to limit the present invention in any form,
Every any simple modification made to the above embodiment according to the technical essence of the invention, equivalent variations and modification, belong to
In the range of technical solution of the present invention.
Claims (6)
1. a kind of quartz sand gas liquid inclusion three dimensional PIC code system, which is characterized in that the quartz sand gas liquid inclusion three
Tieing up mechanical analysis system includes：
Cubing module, connect with main control module, the volume data for measuring gas liquid inclusion；The survey of cubing module
Amount method includes：
Relevant database data, sensing data and controller data that cubing module embeds are uploaded by Sqoop
To distributed file system HDFS, and store into NoSQL databases；Using MapReduce Computational frames to relational data
Library data, sensing data and controller data carry out mining analysis, NoSQL databases are written in the data analyzed, and lead to
Cross Web displayings；
Apriori association rule minings are used under MapReduce frames, are specifically included：
The set L of frequent 1 item collection is obtained using MapReduce computation module1, generate the set C of candidate's k item collectionsk(k≥2)；
In Map function processing stages, Transaction Information of each Map task computations handled by it, which is concentrated in each transaction journal, includes
In CkIn Item Sets occurrence number, for each Map tasks, if some item collection of candidate's k item collections (includes k
Project) it appears in a transaction journal, then Map functions are generated and are exported<Some item collection, 1>Key-value pair gives Combiner letters
Number gives Reduce functions after being handled by Combiner functions；
In Reduce function processing stages, Reduce functions add up CkIn Item Sets occurrence number, obtain all items collection
Support frequency, all minimum Item Sets composition frequent item set L for supporting frequency for supporting frequency >=settingkSet, if k ＜
Maximum iterations and not be sky, then execute k++, be transferred to step Map function processing stages；Otherwise, terminate operation；
Main control module, with pressure detecting module, temperature detecting module, cubing module, coordinate curve drafting module, three-dimensional mould
Quasi- module, data memory module, display module connection, for dispatching modules normal work；
Coordinate curve drafting module, connect with main control module, for the data of detection to be depicted as coordinate song by mathematical software
Line；
Three-dimensional simulation module, connect with main control module, and dynamic analysis are carried out for gas liquid inclusion to be fabricated to threedimensional model；Tool
Body includes：
Step 1: simplifying to triangle mesh model, to be simplified triangle mesh model；
Step 2: using the obtained simplified triangle mesh model of step 1, spring-mass model is established, triangle is simplified
Final position of the grid by spring-mass model parameterization after rectangular domain；
Step 3: using the process for being simplified triangle mesh model, it is initial in parameterizing rectangular domain to calculate deleted point
Position thus obtains position of all vertex in parameterizing rectangular domain in initial model；
Step 4: establishing spring-mass model to parametrization initial model, it is motionless to fix its border vertices, cause in put iteration,
Until interior position reaches balance, the parametrization coordinate of spring-mass Modifying model is obtained；
Step 5:The parametrization coordinate of the spring-mass Modifying model obtained to step 4 is normalized, as grid
The texture coordinate on vertex carries out texture mapping.
2. quartz sand gas liquid inclusion three dimensional PIC code system as described in claim 1, which is characterized in that
The measurement method of cubing module further includes：
Generate neural network model BP initial weights：By random number of the weight initialization between [a, b], wherein a, b are to meet
The integer of following equation：
Wherein H is network node in hidden layer；
The weights and threshold value of neural network model BP are dynamically refined, dynamic neural network model DBP is obtained, generate dynamic
The weights and threshold value of neural network model DBP；
With adaptive immune genetic AIGA algorithm optimizations dynamic neural network model DBP, prediction model AIGA-DBP, root are obtained
It is predicted that model AIGA-DBP calculates volume predictions value；
The weights and threshold value of neural network model BP are dynamically refined, dynamic neural network model DBP is obtained, generate dynamic
The weights and threshold value of neural network model DBP specifically include following steps：
Adjust the weight w between neural network model BP hidden layers and output layerkj；
Adjust wkjPurpose be desirable to the new output of output node jO is exported than currentlypjCloser to desired value tpj, definition：
Wherein α represents the degree of approach, is remained unchanged in each cycle of training, and becomes smaller with the adjustment of node in hidden layer H, does not examine
Consider threshold value, then has：
Wherein wkjWithRespectively update front and back weights, ypkIt is exported for hidden layer, △ wkjFor wkjKnots modification；
Obtain △ wkjSolution equation：
Wherein,
Equation is solved according to least square and error principleObtain △ wkjApproximate solution：
It is connected to the hidden layer node k of output node j to each, calculates the weights variation △ w between k and jkj, update weights
And error of sum square E is calculated, then in one optimal k of k ∈ [1, H] interval selection so that E is minimum；
Adjust the weights v between neural network model BP input layers and hidden layerik；
Adjust vikPurpose be once neural network algorithm is absorbed in local minimum point, modification weights can jump out the minimal point, sentence
The condition that disconnected neural network algorithm is absorbed in local minimum point is the change rate △ E=0 of error E, and E>0；
Do not consider that threshold value, the change of the weights of hidden layer node k are solved by following equation：
Wherein δpj=f-1(ypk+Δypk)-f-1(ypk), M is natural number, then hidden layer exports ypkSolution formula is：
Wherein △ ypkFor ypkKnots modification, then have：
According to least square and error principle solution formulaThe matrix equation of structure can be calculated
Go out：
Aggregative formulaWithCalculate hidden layer with it is defeated
Go out the consecutive mean knots modification of weights between layer
Calculate the consecutive mean knots modification of weights between input layer and hidden layer
M takes the natural number between 10~20, the consecutive mean weights of neural network model BP is obtained, according to neural network model BP
Consecutive mean weights obtain dynamic neural network model DBP.
It is planted with the power of adaptive immune genetic AIGA algorithm optimization dynamic neural network models DBP and mould is predicted in threshold value, acquisition
Type AIGA-DBP specifically includes following steps：
Antigen recognizing, the i.e. identification of mean square error are carried out first, generate the initial antibody (power of dynamic neural network model DBP
Value and threshold value) after, antibody coding is carried out, then calculating antibody fitness and antibody concentration, are based on according to the result of calculating
The adaptive antibody selection operation of concentration obtains suitable weights and threshold value, and result is defeated if this operation meets condition
Go out；If this operation is unsatisfactory for condition, need to carry out adaptive crossover operation and TSP question operation, and carry out weights and threshold value
Update, then recalculate antibody fitness and antibody concentration, weights and threshold value until obtaining the condition that meets are last according to full
The weights and threshold value of sufficient condition obtain prediction model AIGA-DBP.
3. quartz sand gas liquid inclusion three dimensional PIC code system as described in claim 1, which is characterized in that it is described Step 1:
Triangle mesh model is simplified, to be simplified triangle mesh model, is included the following steps：
(1) triangle grid model initialization is carried out, every directed edge in triangle gridding is traversed, if this side is only by a triangle
Including two vertex for then illustrating to connect this side are boundary point；
(2) input needs the number of vertices x retained and angle threshold value θ；
(3) judge whether existing number of vertices is more than vertex threshold value x；If so, carrying out (4)；If it is not, then carrying out (9)；
(4) length of side on each side, and the normal vector of each triangle are found out respectively, and directed edge is calculated using following formula
Collapse cost cost (eij)；And find out the directed edge e of collapse cost minimumuv, then carried out (5)；
cost(eij)=| | eij||·w(eij) (2)
Wherein w (eij) indicate directed edge eijLocal curvature, TiIt indicates to include vertex viTriangle sets,It indicates to include side
eijTriangle sets, f.normal refers to the normal vector in the faces f, and n.normal refers to the normal vector on the faces n；It, which is illustrated, includes
Directed edge eijTri patch and include vertex viTri patch normal direction maximum value；
(5) judge whether u, v are boundary point.If so, carrying out (6)；If it is not, then carrying out (7)；
(6) judge whether the rim angle at boundary point u is more than threshold θ；If so, carrying out (8)；If it is not, then carrying out (3)；
(7) it is boundary point to judge whether u, v meet u, and v is not boundary point；If satisfied, then carrying out (3)；If not satisfied, then carrying out
(8)；
(8) vertex u is deleted, the frontier juncture being connect with vertex u is connected into vertex v.And it is affected point to mark 2 points of u, v,
It carries out (3)；
(9) terminate；
It is described Step 2: using the obtained simplified triangle mesh model of step 1, establish spring-mass model, be simplified
Final position of the triangle gridding by spring-mass model parameterization after rectangular domain, includes the following steps：
(a) triangle grid model is read in, point mass m, device of spring stiffness coefficient k, time step Δ t, iteration thresholding e are initialized；
(b) triangle grid data pretreatment is carried out, triangle grid model boundary point sequence is obtained, to there is the abutment points sequence on each vertex
Row；
(c) perspective plane is selected according to formula (3), (4), triangle grid model is projected, two-dimensional coordinate system is established on perspective plane,
Calculate the two-dimensional coordinate (x after the projection of triangle gridding vertexi,yi)；
a(x-x0)+b(y-y0)+c(z-z0)=0 (3)
Wherein, (a, b, c) is the normal vector on perspective plane, x0,y0,z0It is a fixed point, (a in spacei,bi,ci) it is i-th of triangle
Normal vector, n is the number of Irregular Boundary Surface on piece triangle；
(d) four vertex of Selecting All Parameters rectangular area four sides are chosen according to the structure feature of the triangle gridding after projection
Boundary's point, is fixed on the vertex of given rectangular domain；
(e) due to the movement of point, the stressing conditions for causing on-fixed point change；According to formula (5), calculate suffered by on-fixed point
Resultant force f (Pi)；If for arbitrary point Pi, suffered resultant force is respectively less than given threshold value, i.e. f (Pi)≤e, then terminate current iteration, and first
Turn (i) when secondary termination iteration, turns (j) when terminating iteration for the second time；Otherwise, continue；
Wherein, M is the mass matrix of spring-mass model system, and q is particle coordinate, and K is stiffness matrix；
(f) according to formula (6), by resultant force f (Pi) calculate the acceleration each put
Wherein, m is point mass,For in t moment particle PiAcceleration；
(g) by accelerationAnd given time step Δ t, calculate the displacement each put according to formula (7), (8)；
Wherein,WithThe respectively particle P in t and t+ time ΔtsiSpeed,WithRespectively in t and t+ Δ t
Carve particle PiPosition；
(h) by the displacement in (g), each point is calculated after the Δ t times, the current location of pointTurn (e)；
(i) by loose boundary point according in space lattice between boundary point and boundary point broken line ratio, be fixed to given
Four sides of rectangular domain；It is further continued for carrying out (i)；
(j) parameterized results are obtained.
4. quartz sand gas liquid inclusion three dimensional PIC code system as described in claim 1, which is characterized in that the quartz sand gas
Liquid inclusion enclave three dimensional PIC code system further includes：
Pressure detecting module, connect with main control module, the pressure data for detecting gas liquid inclusion by pressure sensor；
Temperature detecting module is connect with main control module, the temperature data for detecting gas liquid inclusion by temperature sensor；
Data memory module is connect with main control module, is stored for that will analyze data；
Display module is connect with main control module, for showing detection data and threedimensional model.
5. quartz sand gas liquid inclusion three dimensional PIC code system as described in claim 1, which is characterized in that the cubing
Module measurement method is as follows：
First, it is scanned using laser confocal microscope, obtains quartz sand gas liquid inclusion three-dimensional image；
Secondly, the boundary of three-dimensional coordinate calibration three-dimensional image liquid phase is utilized under transmitted light and fluorescence；
Then, according to the liquid phase boundary of calibration, the volume of quartz sand gas liquid inclusion liquid phase is calculated；
Then, gas phase diameter d and the gas in quartz sand gas liquid inclusion are accurately obtained using laser co-focusing micro-fluorescence spectrum
Phase radius R；
Finally, the volume V of gas phase in inclusion enclave is calculated according to formula V=(4/3) π R3.
6. quartz sand gas liquid inclusion three dimensional PIC code system as described in claim 1, which is characterized in that the three-dimensional simulation
Module includes modelling module, parameter setting module, deduction module；
Modelling module, for gas liquid inclusion to be fabricated to threedimensional model by 3d prototype softwares；
Parameter setting module, for gas liquid inclusion to be arranged by force parameter；
Deduction module, for the dynamic state of parameters deduction gas liquid inclusion stressing conditions by being arranged.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810368507.0A CN108563890A (en) | 2018-04-23 | 2018-04-23 | A kind of quartz sand gas liquid inclusion three dimensional PIC code system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810368507.0A CN108563890A (en) | 2018-04-23 | 2018-04-23 | A kind of quartz sand gas liquid inclusion three dimensional PIC code system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108563890A true CN108563890A (en) | 2018-09-21 |
Family
ID=63536160
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810368507.0A Pending CN108563890A (en) | 2018-04-23 | 2018-04-23 | A kind of quartz sand gas liquid inclusion three dimensional PIC code system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108563890A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109632016A (en) * | 2019-02-20 | 2019-04-16 | 湖北理工学院 | Rock And Soil adit digging and surrouding rock stress, strain monitoring experimental rig and its method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104636479A (en) * | 2015-02-15 | 2015-05-20 | 西安电子科技大学 | Industrial big data driven total completion time prediction method |
| CN104680573A (en) * | 2015-01-20 | 2015-06-03 | 西安科技大学 | Texture mapping method based on triangular mesh simplification |
| CN104849256A (en) * | 2015-04-15 | 2015-08-19 | 中国地质大学（武汉） | Method for obtaining trapping pressure of pure methane inclusion |
| CN204917995U (en) * | 2015-08-26 | 2015-12-30 | 湖南工业大学 | Device of gas -liquid inclusion enclave in rejection glass sand |
| WO2017052911A1 (en) * | 2015-09-21 | 2017-03-30 | Siemens Product Lifecycle Management Software Inc. | Method, system, and non-transitory computer readable medium for distributing multiple layers of a composite within a structural volume containing an inclusion |
| CN106705835A (en) * | 2015-07-14 | 2017-05-24 | 中国石油化工股份有限公司 | Method of precisely measuring volume of hydrocarbon inclusion |
-
2018
- 2018-04-23 CN CN201810368507.0A patent/CN108563890A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104680573A (en) * | 2015-01-20 | 2015-06-03 | 西安科技大学 | Texture mapping method based on triangular mesh simplification |
| CN104636479A (en) * | 2015-02-15 | 2015-05-20 | 西安电子科技大学 | Industrial big data driven total completion time prediction method |
| CN104849256A (en) * | 2015-04-15 | 2015-08-19 | 中国地质大学（武汉） | Method for obtaining trapping pressure of pure methane inclusion |
| CN106705835A (en) * | 2015-07-14 | 2017-05-24 | 中国石油化工股份有限公司 | Method of precisely measuring volume of hydrocarbon inclusion |
| CN204917995U (en) * | 2015-08-26 | 2015-12-30 | 湖南工业大学 | Device of gas -liquid inclusion enclave in rejection glass sand |
| WO2017052911A1 (en) * | 2015-09-21 | 2017-03-30 | Siemens Product Lifecycle Management Software Inc. | Method, system, and non-transitory computer readable medium for distributing multiple layers of a composite within a structural volume containing an inclusion |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109632016A (en) * | 2019-02-20 | 2019-04-16 | 湖北理工学院 | Rock And Soil adit digging and surrouding rock stress, strain monitoring experimental rig and its method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | A deformable surface model for real-time water drop animation | |
| US9953460B2 (en) | Garment simulation using thread and data level parallelism | |
| CN104680573B (en) | A kind of texture mapping method based on Triangular Mesh Simplification | |
| Ng et al. | Computer graphics techniques for modeling cloth | |
| CN106407605B (en) | A kind of particlized computer dynamic stimulating method of three-dimensional garment | |
| Iben et al. | Generating surface crack patterns | |
| CN105793866B (en) | Method and system for Intelligent crane lifting | |
| CN101944144B (en) | Meshless cloth-based simulation method | |
| CN107590853A (en) | A kind of high validity methods of exhibiting of architecture ensemble earthquake | |
| CN106646645B (en) | A kind of gravity forward modeling accelerated method | |
| CN103080941A (en) | Apparatus for generating computational data, method for generating computational data, and program for generating computational data | |
| CN108491619A (en) | Based on physics and the non-physical complex scene fluid structurecoupling efficient analogy method mixed | |
| CN108053483A (en) | A kind of Wei Nuotu three-dimensional grid reconstructing methods accelerated based on GPU | |
| CN106934192B (en) | A kind of shallow water equations model water body modeling method of parameter optimization | |
| CN108776745A (en) | A kind of complex topography wind field analogy method and device based on WRF and Fluent couplings | |
| Ruspini et al. | Haptic display for human interaction with virtual dynamic environments | |
| Feng et al. | Learning-based bending stiffness parameter estimation by a drape tester | |
| Larionov et al. | Estimating cloth elasticity parameters using position-based simulation of compliant constrained dynamics | |
| CN108563890A (en) | A kind of quartz sand gas liquid inclusion three dimensional PIC code system | |
| Chen et al. | Real-time continuum grass | |
| CN108171745A (en) | A kind of three-dimensional plant plant projected area computational methods | |
| CN104299265A (en) | Group behavior control technology under fluid environment influence | |
| Rudomin et al. | Gpu generation of large varied animated crowds | |
| O'Brien | Graphical modeling and animation of fracture | |
| CN114912303A (en) | Three-dimensional fluid display method and device, electronic equipment and storage medium |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180921 |