CN110334414A - The calculation method and device of polymer degradation intensity based on intensity phasor - Google Patents
The calculation method and device of polymer degradation intensity based on intensity phasor Download PDFInfo
- Publication number
- CN110334414A CN110334414A CN201910525798.4A CN201910525798A CN110334414A CN 110334414 A CN110334414 A CN 110334414A CN 201910525798 A CN201910525798 A CN 201910525798A CN 110334414 A CN110334414 A CN 110334414A
- Authority
- CN
- China
- Prior art keywords
- intensity
- slight crack
- state
- cellular
- region
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Image Analysis (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
The embodiment of the present invention discloses the calculation method and device of a kind of polymer degradation intensity based on intensity phasor, which comprises step 1, reads the intensity phasor of polymer to be calculated；Step 2, it is begun stepping through from the coordinate origin of the intensity phasor, judges the state of strength of the cellular on each pixel；Step 3, according to neighbour-border extension algorithm, region composed by the pixel to the same intensity state in the phasor carries out hole, slight crack identification, obtains the region inside aperture of varying strength state, the region area and concentration of slight crack；Step 4, according to the region inside aperture of the varying strength state, the region area and concentration of slight crack, and the strength model built up in advance, the intensity of the cellular of varying strength state is calculated；Step 5, according to the intensity of the cellular of the varying strength state, the overall strength of the polymer is calculated.
Description
Technical field
The present invention relates to degradable high polymer material field more particularly to a kind of polymer degradation based on intensity phasor are strong
The calculation method and device of degree.
Background technique
Closely now, tend to be mature along with medicine continuous development, degradable high polymer material also grows up therewith.Pass through
The new medical product that PLA and its copolymer are made are gushed, and excellent performance is especially shown on organization bracket,
Approved by society research personnel and doctor as organization bracket, support strength when degradation will be with the growth rate of tissue
Match, therefore the research about material degradation Strength Changes has become the most active topic of current research.For material
The research of degradation process be it is gradually deepening, not accomplish in one move.Zygourakis is proposed to utilize the randomness of Monte Carlo, be built
The drug release that vertical mathematical model be attached in degradation material is simulated, and is for the first time used for monte carlo method high poly-
In the random degenerate of object.Thombre and Himinelstein is the influence for being put forward for the first time reaction environment to degradation, if acidity
When, ester bond meets sour water solution, and hydrionic concentration increases so that degradation system polarity increases, the self-catalysis ability of degradation material compared with
By force, accelerated degradation carries out.Et al. mainly the influence to the diffusion of short molecular chain to degradation carried out point
Analysis.Professor Pan proposes the model of elastic modulus change for the hydrolysis of hypocrystalline degradation material and the diffusion of short chain,
Think that the intensity of crystallization is bigger than hemicrystalline intensity, and empty intensity is zero.But the equation does not take into account crystallization.
Calculation formula of professor Shirazi based on the proposed elasticity modulus of Flory, by molecular weight model and elasticity modulus model phase coupling
It closes, is to study one of model the most accurate in high polymer Strength Changes up to now.The size shadow of the discovery field trash such as Zhang
Intensity needed for ringing interface peel, the size of intensity are related with the size of interface impurity size.Tanaka etc. is adulterated in the material
Spherical impurity not of uniform size, discovery contaminant size is bigger, and stress needed for destroying interface is smaller.Numerous scientific research personnel not
It is disconnected to explore under research, it is still not yet found so far about crystallization, the influence in gap, molecular weight and constituent size on strength
Effect, therefore, analysis are faced with huge difficulty with intensity is calculated, and have extremely challenging.
Summary of the invention
In view of this, the embodiment of the present invention provide it is a kind of based on intensity phasor polymer degradation intensity calculation method,
Device can calculate polymer degradation intensity according to approximating method.
A kind of calculation method of the polymer degradation intensity based on intensity phasor, comprising:
Step 1, the intensity phasor of polymer to be calculated is read；
Step 2, it is begun stepping through from the coordinate origin of the intensity phasor, judges the intensity shape of the cellular on each pixel
State；
Step 3, according to neighbour-border extension algorithm, the pixel of the same intensity state in the phasor is formed
Region carry out hole, slight crack identification, obtain the region inside aperture of varying strength state, the region area of slight crack and intensive journey
Degree；
Step 4, according to the region inside aperture of the varying strength state, the region area and concentration of slight crack, and
The strength model built up in advance calculates the intensity of the cellular of varying strength state；
Step 5, according to the intensity of the cellular of the varying strength state, the overall strength of the polymer is calculated.
A kind of computing device of the polymer degradation intensity based on intensity phasor, comprising:
Reading unit reads the intensity phasor of polymer to be calculated；
Judging unit begins stepping through from the coordinate origin of the intensity phasor, judges the strong of the cellular on each pixel
Degree state；
Acquiring unit, according to neighbour-border extension algorithm, to the pixel institute group of the same intensity state in the phasor
At region carry out hole, slight crack identification, obtain the region inside aperture of varying strength state, the region area of slight crack and intensive journey
Degree；
First computing unit, according to the region inside aperture of the varying strength state, the region area of slight crack and intensive journey
Degree, and the strength model built up in advance calculate the intensity of the cellular of varying strength state；
Second computing unit calculates the total of the polymer according to the intensity of the cellular of the varying strength state
Intensity.
In the present invention, polymer degradation intensity can be calculated according to approximating method, comparison for calculation methods is simple.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
Other attached drawings are obtained according to these attached drawings.
Fig. 1 is the schematic diagram of the calculation method of polymer degradation intensity of the embodiment of the present invention based on intensity phasor；
Fig. 2 is the signal of the calculation method of polymer degradation intensity of the another embodiment of the present invention based on intensity phasor
Figure；
Fig. 3 is another embodiment of the present invention neighbour-border extension algorithm schematic diagram；
Fig. 4 is the definition of slight crack in the present invention；
Fig. 5 is the definition of neighbour-border extension algorithm hole；
Fig. 6 is the definition of neighbour-border extension algorithm complicated shape；
The mixing macromolecule degradable material degradation Strength Changes figure that Fig. 7 is 60%PLLA and 40%PVA, yellow intensity state
For amorphous phase, red color intensity state is crystalline phase, and ash, black intensity state are Kongxiang；
Fig. 8 is the mould of the molecular weight of 60%PLLA and 40%PVA material provided in an embodiment of the present invention, crystallinity, intensity
Analog values are compared with experiment value, and green is crystallinity, and blue is normalized intensity, and red is normalization molecular weight；
Fig. 9 is PLA degradable material degradation Strength Changes figure, and yellow intensity state is amorphous phase, and red strong 0 degree of state is
Crystalline phase, ash, black intensity state are Kongxiang
Molecular weight, crystallinity, the analogue value of intensity and the experiment value ratio of Figure 10 PLA material provided in an embodiment of the present invention
Right, green is crystallinity, and blue is normalized intensity, and red is normalization molecular weight；
Figure 11 is the connection figure of the computing device of polymer degradation intensity of the embodiment of the present invention based on intensity phasor.
Specific embodiment
The embodiment of the present invention is described in detail with reference to the accompanying drawing.
It will be appreciated that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Base
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts it is all its
Its embodiment, shall fall within the protection scope of the present invention.
For convenience of description, description apparatus above is to be divided into various units/modules with function to describe respectively.Certainly, exist
Implement to realize each unit/module function in the same or multiple software and or hardware when the present invention.
As shown in Figure 1, for a kind of calculation method of the polymer degradation intensity based on intensity phasor of the present invention, packet
It includes:
Step 1, the intensity phasor of polymer to be calculated is read；The step 1 specifically:
The intensity phasor of polymer to be calculated is read, the pixel of the intensity phasor is L*L；L is greater than predetermined value.L
It may be set according to actual conditions.
Step 2, it is begun stepping through from the coordinate origin of the intensity phasor, judges the intensity shape of the cellular on each pixel
State；
Step 3, according to neighbour-border extension algorithm, the pixel of the same intensity state in the phasor is formed
Region carry out hole, slight crack identification, obtain the region inside aperture of varying strength state, the region area of slight crack and intensive journey
Degree；
The step 3 includes:
S301 is successively read the phasor according to time sequencing, carries out hole, slight crack identification；
S302 is begun stepping through from the coordinate origin of cellular, and the region composed by the object pixel is horizontal and vertical both greater than
When equal to N, it is judged as hole；The region composed by the object pixel only meet longitudinal length be greater than/or be equal to N, or it is only full
When sufficient lateral length is greater than or equal to N, it is judged as slight crack；
S303 first counts the number of hole and area；If target pixel points belong to the hole, accessed time
Going through traffic sign placement is 1；
S304, then slight crack is traversed, specifically: longitudinal traversal is carried out to slight crack first, if pixel access time
It goes through and is identified as 0, and belong to slight crack, its access flag is arranged to 1, slight crack is added in the set of storage slight crack；Longitudinal traversal
After, then lateral traversal is carried out to slight crack, it is identical as longitudinal method, and then respectively obtain the set of slight crack and hole；
S305, according to phasor picture or state found out come slight crack, hole set, loop through；Again by current
The coordinate of slight crack or hole finds out minimum and maximum coordinate value, determines the range of current region；
S306 judges whether the coordinate in current slight crack is labeled；If marked, perform the next step rapid；No
Label, then carry out S304；
S307 carries out outside rectangle expansion to the range of current slight crack, increases N number of pixel for the first time, second of increase
2*N and so on, current extent is traversed, checks in slight crack or hole set whether to have to be present in currently being traversed
In the range of；If it does, adjacent slight crack or hole set is added, and current adjacent slight crack number adds 1；If except when
Preceding slight crack or hole are adjacent without other regions, or when constituting the gross area in region and accounting for less than expand area 1/3, then stop
Expand, above-mentioned adjacent slight crack is stored in set, circulation terminates；
S308 successively traverses each of slight crack set slight crack or hole, until terminating, obtains varying strength state
Region inside aperture, the region area of slight crack and concentration.
Step 4, according to the region inside aperture of the varying strength state, the region area and concentration of slight crack, and
The strength model built up in advance calculates the intensity of the cellular of varying strength state；
The step 4 includes:
When the state of strength of the cellular is crystalline phase, the intensity of the cellular is calculated specifically:
In formula, σZFor crystalline phase intensity value, αzhFor the number of slight crack in region, αzcFor the number of region inside aperture, SzhIt is
The area summation of slight crack, S in regionzcIt is the area summation of region inside aperture, φ,It is predetermined ginseng related with material properties
Number, l × l are cellular number；
When the state of strength of the cellular is amorphous phase, the intensity of the cellular is calculated specifically:
In formula, σaFor crystalline phase intensity value, αahFor the number of slight crack in region, αacFor the number of region inside aperture, SahIt is
The area summation of slight crack, S in regionacIt is the area summation of the region inside aperture, γ, η are predetermined ginsengs related with material properties
Number, l × l are cellular number；
When the state of strength of the cellular is Kongxiang, the intensity of the cellular is calculated specifically:
In formula, σeFor Kongxiang intensity value, αehFor the number of slight crack in region, αecFor the number of region inside aperture, SehIt is area
The area summation of slight crack, S in domainecIt is the area summation of region inside aperture, μ, ν are predefined parameters related with material properties, l ×
L is cellular number.
The method for pre-establishing model includes:
Step 21, the intensity phasor of the polymer of modeling is read；
Step 22, it is begun stepping through from the coordinate origin of the intensity phasor of modeling, judges the member on each pixel
The state of strength of born of the same parents；
Step 23, according to neighbour-border extension algorithm, to the pixel of same intensity state of the modeling in the phasor
Composed region carries out hole, slight crack identification, obtains the region inside aperture of varying strength state, the region area of slight crack and close
Collection degree；
Step 24, the described of varying strength state is calculated according to the initial value of model according to the state of strength of the cellular
The intensity of cellular；
Step 25, according to the intensity of the cellular of the varying strength state, the modeling overall strength of polymer is calculated；
Step 26, modeling is compared with overall strength with the experiment value of the overall strength of the polymer of modeling, is constantly adjusted
The parameter of model, obtains optimal models, as the model built up.
Step 5, according to the intensity of the cellular of the varying strength state, the overall strength of the polymer is calculated.
The step 5 includes:
σ (t)=XA(t)σa+XZ(t)×σZ+XE(t)×σe
In formula,For overall strength, XA(t)、XZ(t)、XE(t) under respectively amorphous, crystallization, three kinds of Kongxiang state
Cellular accounts for the ratio of cellular sum；σZFor crystalline phase intensity value, σaFor crystalline phase intensity value, σeFor Kongxiang intensity value.
The present invention provides a kind of polymer degradation intensity modeling and calculation method based on phasor image procossing, bio-medical
Degradable high polymer material is because its good mechanical property is widely used with degradation property in medicine, intensity in degradation process
The matching of variation and degradation rate is the critical consideration in the design of high polymer device.Establish polymer degradation processes
Strength model can be developed with the intensity in quantitative research degradation process.Corroding structure change in degradation process with polymer body is to cut
Access point reflects the differentiation of the phase of intensity with the variation of structure, the concept of intensity phasor is proposed, in the method for image procossing to strong
It spends phasor and carries out feature identification, propose that neighbour-border extension algorithm (NNBE) identification has the feature of hole, slight crack shape, with
This feature establishes out of phase strength model.Model is calculated and is compared with experimental data, pattern die analog values and experiment value are quasi-
It is preferable to close effect, illustrates that strength model and feature of phase diagram extraction algorithm are feasible.
A kind of polymer degradation intensity modeling and calculation method based on phasor image procossing,
The cellular of high polymer discretization indicates；
When initial, each member is intracellular, and there are a macromolecular chains；
Material degradation, macromolecular chain are broken；
The macromolecular chain being broken, it may occur however that recrystallization, and diffusion can occur for short chain；
Average molecular weight, macromolecule chain quantity, crystallinity of cellular etc. can all make state of strength change, thus
Obtained phasor is different；
State of strength under different moments, discribed phasor are different；
Phasor is read, different pixels point represents different phasor states；
Under same intensity state, to region composed by pixel, hole, slight crack identification are carried out；
Calculate the region area of hole, slight crack；
The concentration of hole, slight crack is differentiated；
According to hole, the region area of slight crack and concentration, intensity is calculated；
The polymer degradation intensive analysis based on phasor image procossing and computation model are for predicting cellular in t moment
Under intensity value, and be considered as under t moment, the intensity value of material entirety.
In the present invention, different phasors is obtained from the definition of intensity angle, using this phasor as hole, slight crack recognizer
Input, begun stepping through from cellular coordinate origin.When the definition of the region hole composed by the object pixel, mark it for hole
Hole；When meeting the definition of slight crack, it is labeled as slight crack；It is respectively stored into the set of hole and slight crack.
In the present invention, pixel is traversed, the area (use number of pixels) of record hole, slight crack, area is bigger, hole slight crack
More intensive, the influence to intensity is bigger, further determines that the concentration of hole, slight crack.
The gross area occupied area domain gross area of the area of region inseam, seam is bigger, bigger to intensity effect；Because of region inner hole
The gross area occupied area domain gross area of area, hole is bigger, bigger to intensity effect.
In the present invention, according under moment t, cellular number-average molecular weight, macromolecule chain number, the state of strength of cellular are discribed
Phasor.Read phasor, varying strength state use different strength calculation formulas, finally by it is stateful under intensity value phase
Add, i.e. overall strength of the material under t moment.
In the present invention, the diffusion formula of short chain (unit number is less than 8 on macromolecular chain) is calculated are as follows:
Wherein, CeFor ester bond concentration；ColFor low short chain concentration；r1Under the conditions of normal；Hydrolytic rate constant；r2For acidity
Under the conditions of, hydrolytic rate constant；D is diffusion coefficient；Grad is gradient；Div is divergence.
In the present invention, each state is calculated using the strength formula under the state according to cellular state of strength under t moment
Intensity include:
First, crystalline phase strength model
Based on phasor, proposed according to crystalline phase shape using the hole of proposition, slight crack recognizer with the face of crystalline phase
Influence of the long-pending and concentration to intensity, calculates the intensity value of crystalline phase, such as following formula:
σz=σzh+σzc
In formula, σzFor crystalline phase intensity value, σzhInfluence for slight crack to intensity, σzcInfluence for hole to intensity.
Second, amorphous phase strength model
Strength Changes calculation formula under amorphous is as follows:
σa=σaA+σac
In formula, σaFor crystalline phase intensity value, σahInfluence for slight crack to intensity, σacInfluence for hole to intensity.
Third, Kongxiang strength model
Strength Changes calculation formula under Kongxiang is as follows:
σe=σeh+σec
In formula, σeFor Kongxiang intensity value, σehInfluence for slight crack to intensity, σecInfluence for hole to intensity.
In the present invention, the integral strength changed by each cellular intensity statistics polymer with degradation time includes:
Intensity value under three kinds of states is added up, overall strength is calculated:
In formula,For required overall strength, XA(t)、XZ(t)、XEIt (t) is amorphous, crystallization, three kinds of Kongxiang state
Under cellular account for the ratio of cellular sum.
The division number accuracy of calculated result is related with the size of number L is divided in the present invention, when L value is larger,
Time required for computer simulation polymer degradation process is long, occupies resource；When L value is smaller, calculated result accuracy
It is low, the accuracy of model can not be embodied.
Now still not yet find the influence about crystallization, gap and molecular weight to intensity so far in the art, because
This, analysis is faced with huge difficulty with intensity is calculated.The present invention can exist for simulation and forecast high polymer is carried out based on phasor
The variation of mechanical property in degradation process provides a kind of polymer degradation intensive analysis based on phasor image procossing and calculates
Method, the method by reading phasor, calculates the intensity value under high polymer varying strength phase state.
As shown in Fig. 2, the present invention provides a kind of polymer degradation intensive analysis based on phasor image procossing and calculating
Method, comprising:
1) under different moments, intensity phasor is different, by reading phasor, it can be deduced that state of strength；
2) state of strength is divided into three classes: crystalline phase, amorphous phase, Kongxiang；
3) under different conditions, the calculation formula of intensity is different；
4) in the case where determining to state of strength, hole, slight crack identification are carried out；
5) area and number of hole, slight crack are calculated；
6) it when hole, slight crack are intensive, is affected to intensity, therefore, concentration is analyzed；
7) by concentration and area, strength calculation formula is substituted into；
8) intensity value under three kinds of state of strength is added, i.e. overall strength.
Above-mentioned 8 points are specifically addressed below.
Firstly, the variation for the phase that polymer degradation processes occur is inquired into meso-scale, polymer is divided into Jie
The cellular of sight, it is assumed that each cellular contains a macromolecular chain when degrading initial, and to simulate chain break, chain break causes point
The reduction of son amount；The oligomer of fracture diffuses out matrix, forms hole；Chain is broken the formation for also resulting in recrystallization；Go out as a result,
Structural form different in degradation process is showed.From the angle of structural form can separate unformed phase, crystalline phase and Kongxiang it
The phase of class.It has been investigated that influence of the sufficiently small unformed phase cellular of number molecular weight to matrix strength and little, thus from
The angle of intensity can treat as the small unformed phase of these molecular weight as Kongxiang, for this purpose, we are again fixed from the angle of intensity
Justice phase is respectively as follows: degradation phase, crystalline phase and intensity Kongxiang.
1) degradation phase
Refer to that chain is broken, molecular weight becomes smaller, but molecular weight does not have small arrive to be lower than some critical value, this critical value is known as
Strength Failure critical molecular weight；
2) crystalline phase
The cellular of all crystallizations (including the recrystallization generated in crystallization existing before degrading and degradation process) shows knot
Crystal phase.
3) intensity Kongxiang
Oligomer diffuse out the hole cellular presented outside matrix and due to caused by degradation molecular weight lower than Strength Failure
The cellular of critical molecular weight, define these cellulars composition is mutually intensity Kongxiang.
Further, the process degraded with polymer, these three mutually also dynamically carry out evolution, to intensity
Influence the also dynamic change with degradation process.Therefore, by capturing the phase of dynamic change, and quantitative research is unfolded from the angle of phase
To disclose micro-, mesoscopic structure variation to the influence of intensity, the quantitative research to phasor is unfolded using image processing techniques.
It is further carried out using the method for image procossing, phasor is read, according to pixel institute's generation of different colours in phasor
The difference of table meaning differentiates the state of strength of cellular.It is the cellular of crystalline phase or Kongxiang to state of strength, according to compositing area shape
The difference of shape, can be divided into hole, two kinds of shape areas of slight crack are calculated, and variation of the concentration in region to intensity value
It has a huge impact
Further, hole, slight crack are defined, as follows:
1) definition of slight crack
The cellular combination of the adjacent formation strip of one direction (laterally or longitudinally) is defined as slight crack, such as Fig. 1 (a), (b) institute
Show.
2) definition of hole
The adjacent cellular combination of both direction (horizontal and vertical) is defined as hole, as shown in Figure 5.
Wherein for irregular shape, as shown in Fig. 5 (b), (c), it is equal to Fig. 5 (a) processing.
3) complicated shape
For complicated shape, as shown in Fig. 6 (a), the assembly of hole and slight crack can be considered, such figure is split as hole
Hole (Fig. 6 (b)) and slight crack (Fig. 6 (c)), to handled respectively according to hole and slight crack the Influencing Mechanism of intensity.
Further, different phasors is obtained from the definition of intensity angle, using this phasor as hole, slight crack recognizer
Input, begun stepping through from cellular coordinate origin.When the definition of the region hole composed by the object pixel, mark it for hole
Hole；When meeting the definition of slight crack, it is labeled as slight crack；It is respectively stored into the set of hole and slight crack.With traversal pixel, record
The area (using number of pixels here) of hole, slight crack, area is bigger, hole slight crack is more intensive, and the influence to intensity is bigger.
The mechanical strength for further calculating each cellular, in synchronization, it is total that material has crystalline state, amorphous state, hole state
The case where depositing, the intensity value of three kinds of varying strength states, which calculates, uses different calculation formula.
1) crystalline phase strength model
Based on phasor, proposed according to crystalline phase shape using the hole of proposition, slight crack recognizer with the face of crystalline phase
Influence of the long-pending and concentration to intensity, calculates the intensity value of crystalline phase, such as following formula:
In formula, σzFor crystalline phase intensity value, αzhFor the number of slight crack, αzcFor the number of hole, SzhIt is slight crack in the region
Area summation, SzcIt is the area summation of the region inside aperture, φ,It is parameter related with material properties, l × l is cellular
Number.
2) amorphous phase strength model
Strength Changes calculation formula under amorphous is as follows:
In formula, σaFor crystalline phase intensity value, αahFor the number of slight crack, αacFor the number of hole, SahIt is slight crack in the region
Area summation, SacIt is the area summation of the region inside aperture, γ, η are parameters related with material properties, and l × l is cellular
Number.
3) Kongxiang strength model
The cellular of Kongxiang includes that two kinds of forms exist, and one is the cellular that state is hole state, another kind is that state is drop
It solves under state or non-degradable state, but relative molecular mass is less than defined molecular weight threshold, it is believed that it loses supporting role
Cellular, its state of strength are referred to as Kongxiang.According to the form of Kongxiang, phasor Hole, the region area of slight crack and concentration,
Intensity connection is established, proposes Kongxiang strength model are as follows:
In formula, σeFor Kongxiang intensity value, αehFor the number of slight crack, αecFor the number of hole, SehIt is slight crack in the region
Area summation, SecIt is the area summation of hole in the region, μ, ν are parameters related with material properties, and l × l is cellular
Number.
Further, by amorphous phase, crystalline phase, three kinds of Kongxiang strength type, using different strength calculation formulas into
Row calculates, and can be expressed as follows:
The intensity value under three kinds of states is added up in region, calculates overall strength:
In formula,For required overall strength, XA(t)、XZ(t)、XEIt (t) is amorphous, crystallization, three kinds of Kongxiang state
Under cellular account for the ratio of cellular sum.
Further, the phasor being successively read under different moments calculates intensity value.
The advantageous effects of the above technical solutions of the present invention are as follows:
In above scheme, in raw polymer degradation processes, high-molecular copolymer is separated into cellular, and the intensity shape of cellular
State is segmented into three kinds: amorphous phase, crystalline phase, Kongxiang.With degradation process variation constantly occurs for phasor, needs to variation
Phasor carry out Dynamical capture and identification, this be carry out quantitative research premise.By taking intensity Kongxiang as an example, for Kongxiang phasor spy
Kongxiang phasor, can be divided into the features such as hole, slight crack by point and the influence to intensity.The cellular point seen that is situated between can be considered image
Pixel is used for reference the technology of image procossing, is analyzed and processed to phasor with this.Distinguished using different phase strength models
The intensity for calculating the cellular of different conditions counts the intensity of all cellulars up to the elastic strength of material entirety, intends obtaining material
The variation tendency of intensity in degradation process is the precision medicinal equipment of high molecular degradable polymer in mechanical property
The optimization design of aspect provides numerical basis.
Application scenarios of the invention are described below.
The present invention is directed to existing high molecular polymer degradation model, is still not yet found so far about crystallization, gap, is divided
The influence of son amount and constituent size on strength provides a kind of polymer degradation intensity based on phasor image procossing point
Analysis and calculation method.
As shown in figure 3, it is provided in an embodiment of the present invention by the polymer degradation intensive analysis of phasor image procossing with based on
Calculation method proposes the feature that neighbour-border extension algorithm (NNBE) identification has hole, slight crack shape, comprising:
S101 is successively read intensity phasor according to time sequencing, carries out hole, slight crack identification；
S102 is begun stepping through from cellular coordinate origin, and the region composed by the object pixel is horizontal and vertical both greater than etc.
When N, name it for hole；When composition region, which meets longitudinal length, is more than or equal to N or lateral length more than or equal to N, life
Entitled slight crack；
S103, computer first count the number of hole and area, if target pixel points belong to the hole, by it
Access traversal traffic sign placement is 1；
S104, then slight crack is traversed, longitudinal traversal is carried out to slight crack first, if the pixel belongs to slight crack, by it
Access flag is arranged to 1, slight crack is added in the set of storage slight crack, after longitudinal traversal, then carries out lateral time to slight crack
It goes through, it is identical as longitudinal method, and then the set of storage slight crack and hole can be respectively obtained.
S105, according to picture or state found out come slight crack, hole set, loop through.Again by current slight crack
The coordinate of (hole) finds out minimum and maximum coordinate value, determines the range of current region；
S106 judges whether the coordinate in current slight crack is labeled, if marked, without step 4, does not mark
Note, then carrying out step 4；
S107 carries out outside rectangle expansion to the range of current slight crack, increases N number of pixel for the first time, second of increase
Whether 2*N, traverses current extent, check slight crack, have to be present in hole set and currently traversed
In the range of.If it does, adjacent slight crack (hole) set is added, and current adjacent slight crack number adds 1, if except when
Preceding slight crack (hole) without other regions it is adjacent, or constitute region the gross area account for less than expand area 1/3 when, then stop
Expand, above-mentioned adjacent slight crack is stored in set, circulation terminates；
S108 successively traverses each of slight crack set slight crack (hole), until terminating；
S109 carries out Strength co-mputation, obtains intensity value according to concentration and area.
In the present embodiment, the calculating of each cellular intensity is needed using different phase polymer degradation strength models.
Cellular (i, j) may be expressed as: in the intensity of t moment according to its state difference
Intensity value under three kinds of states is added up, overall strength is calculated:
In formula,For required overall strength, XA(t)、XZ(t)、XEIt (t) is amorphous, crystallization, three kinds of Kongxiang state
Under cellular account for the ratio of cellular sum.
To sum up, by high polymer material, it is divided into L*L small grid using cellular automata, each small grid is called
One cellular, neighborhood use four neighborhood of von Neumann.The numerical value of L is bigger, and it is more to represent grid division quantity, to calculate also more
Accurately.But it is there is also a disadvantage, and when the number of grids of division is more, calculation amount is bigger, and time-consuming, may to test into
Degree has an impact.Therefore, a suitable L is selected, is vital.Classified according to state of strength, be segmented into cellular (i,
J) in t moment, four kinds of states: amorphous state, empty state, crystalline state are segmented into.Currently, the forecasting research to intensity mostly stops
In the level of Amorphous Polymers, or only consider its constitutive relation, is not related to the polymer of multiphase coexistence.It is before right
Intensity phase state is the cellular of empty state, it is believed that influence of its intensity to integral strength is defaulted as 0.But in fact, such cellular exists
It is presented in the form of slight crack or hole in phasor, its presence can weaken whole intensity, and viewpoint needs to carry out more before
Just.Therefore, the intensity value under four or three different conditions is solved, add up summation, and the overall strength of material can be obtained.For high score
Optimization design of the precision medicinal equipment of sub- degradable polymer in terms of mechanical property provides numerical basis.The present invention
Different application field suitable for degradable macromolecule equipment.
In order to better understand the present embodiment propose shown polymer degradation intensity modeling and simulation method, be with material
For the mixing macromolecule degradable material and original state of the PVA of 60% PLLA and 40% is the PLA material of nodeless mesh, to this
The high molecular polymer degradation intensity modeling and simulation method that embodiment proposes is described in detail:
Fig. 4 is the definition of slight crack in the present invention；Fig. 5 is the definition of neighbour-border extension algorithm hole；Fig. 6 is neighbour-side
The definition of boundary's expansion algorithm complicated shape；The mixing macromolecule degradable material degradation intensity that Fig. 7 is 60%PLLA and 40%PVA becomes
Change figure, yellow intensity state is amorphous phase, and red color intensity state is crystalline phase, and ash, black intensity state are Kongxiang；Fig. 8 is
Molecular weight, crystallinity, the analogue value of intensity and the experiment value ratio of 60%PLLA and 40%PVA material provided in an embodiment of the present invention
Right, green is crystallinity, and blue is normalized intensity, and red is normalization molecular weight；Fig. 9 is PLA degradable material degradation intensity
Variation diagram, yellow intensity state are amorphous phase, and red strong 0 degree of state is crystalline phase, and ash, black intensity state are Kongxiang；Figure
Molecular weight, crystallinity, the analogue value of intensity and the experiment value of 10 PLA materials provided in an embodiment of the present invention compare, and green is knot
Brilliant degree, blue are normalized intensity, and red is normalization molecular weight.It is described below in conjunction with each figure.
It is the experiment that Tsuji and Muramatsu are done, the material used is 60% PLLA and 40% in the present embodiment 1
PVA mixing macromolecule degradable material.Setting initial strength is 56MPa, and initial molecular weight Mn0 is 9.26 × 104G/mol,
Material nodeless mesh when initial.Degradation Probability p is set as 0.1, and the calculating parameter α and β in the combined type of multinomial strength model are set respectively
For 1.23 and 1.2, φ,For 15 and 1.6, λ,For 23 and 2.6, state of strength figure is expanded 1.3 times, is chosen in phasor
A part observes the variation of its state of strength in different time period, such as Fig. 6.By molecular weight, crystallinity, the analogue value of intensity
It is compared with experiment value, fitting effect is preferable, such as Fig. 8.
Specific step is as follows:
(1) different according to the phasor under different moments, read phasor；
(2) hole, slight crack identification are carried out, after identification, hole, slight crack is separately added into hole, in slight crack set, determine hole
Hole, slight crack region area；
(3) closeness division is carried out, so that it is determined that influencing intensity factor；
(4) according to state of strength, intensity factor and area are substituted into matched strength calculation formula；
(5) by the ratio of the cellular of the state of strength and total cellular number, i.e. accounting, using accounting and the product of intensity value as
Final strength value；
(6) intensity value under state of strength in three is added, obtains material total intensity value；
(7) analog intensity value is compared with laboratory strength value, constantly adjusts ginseng optimization；Obtain optimal models；
In the present embodiment 2, it is that Duek teaches done experiment, uses original state for the PLA material of nodeless mesh, setting is just
Beginning intensity is 190.4MPa, and initial molecular weight Mn0 is 1.61 × 105g/mol.Degradation Probability p is set as 0.05, multinomial strength model
Combined type in calculating parameter α and β be set to 1.35 and 3.6, φ,For 95 and 3, λ,For 51 and 2.6, by intensity shape
State figure expands 5 times, chooses a part in phasor, observes the variation of its state of strength in different time period, such as Fig. 8.It will divide
Sub- amount, crystallinity, the analogue value of intensity are compared with experiment value, and fitting effect is preferable, such as Figure 10.Step as example 1,
It is not repeating.
In the present embodiment 1 and embodiment 2, calculated value and experimental data for sample are compared as shown in Fig. 8, Figure 10, right
Than the correctness of the degradation intensity modeling and simulation method of high molecular polymer described in result verification the present embodiment.Before this
Research, the cellular for defaulting hole state always have no effect to the intensity of material, intensity 0, and the intensity value of crystalline state and crystallization
The relationship of degree and molecular weight is unknown, and the simulation to degradation process is all the formula created according to result, it is understood that there may be uncomfortable
The property used.By read phasor in varying strength state value, by degradation process molecular weight, percent crystallization in massecuite, intensity the analogue value with
Experiment value is compared, and fitting effect is preferable, illustrates that model has good applicability and accuracy.
It as shown in figure 11, is a kind of computing device of the polymer degradation intensity based on intensity phasor of the present invention,
Include:
Reading unit 31 reads the intensity phasor of polymer to be calculated；
Judging unit 32 begins stepping through from the coordinate origin of the intensity phasor, judges the cellular on each pixel
State of strength；
Acquiring unit 33, according to neighbour-border extension algorithm, to the pixel institute of the same intensity state in the phasor
The region of composition carries out hole, slight crack identification, obtains the region inside aperture of varying strength state, the region area of slight crack and intensive
Degree；
First computing unit 34, according to the region inside aperture of the varying strength state, the region area of slight crack and intensive
Degree, and the strength model built up in advance calculate the intensity of the cellular of varying strength state；
Second computing unit 35 calculates the polymer according to the intensity of the cellular of the varying strength state
Overall strength.
Those of ordinary skill in the art will appreciate that realizing all or part of the process in above-described embodiment method, being can be with
Relevant hardware is instructed to complete by computer program, the program can be stored in a computer-readable storage medium
In, the program is when being executed, it may include such as the process of the embodiment of above-mentioned each method.Wherein, the storage medium can be magnetic
Dish, CD, read-only memory (Read-Only Memory, ROM) or random access memory (Random Access
Memory, RAM) etc..
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
In the technical scope disclosed by the present invention, any changes or substitutions that can be easily thought of by those familiar with the art, all answers
It is included within the scope of the present invention.Therefore, protection scope of the present invention should be subject to the protection scope in claims.
Claims (7)
1. a kind of calculation method of the polymer degradation intensity based on intensity phasor characterized by comprising
Step 1, the intensity phasor of polymer to be calculated is read；
Step 2, it is begun stepping through from the coordinate origin of the intensity phasor, judges the state of strength of the cellular on each pixel；
Step 3, according to neighbour-border extension algorithm, area composed by the pixel to the same intensity state in the phasor
Domain carries out hole, slight crack identification, obtains the region inside aperture of varying strength state, the region area and concentration of slight crack；
Step 4, according to the region inside aperture of the varying strength state, the region area and concentration of slight crack, and in advance
The strength model built up calculates the intensity of the cellular of varying strength state；
Step 5, according to the intensity of the cellular of the varying strength state, the overall strength of the polymer is calculated.
2. the method according to claim 1, wherein the step 3 includes:
S301 is successively read the phasor according to time sequencing, carries out hole, slight crack identification；
S302 is begun stepping through from the coordinate origin of cellular, and the region composed by the object pixel is horizontal and vertical to be both greater than equal to N
When, it is judged as hole；The region composed by the object pixel only meet longitudinal length be greater than/or be equal to N, or only meet laterally
When length is greater than or equal to N, it is judged as slight crack；
S303 first counts the number of hole and area；If target pixel points belong to the hole, traversal mark is accessed
Will is set as 1；
S304, then slight crack is traversed, specifically: longitudinal traversal is carried out to slight crack first, if pixel access traversal mark
Knowing is 0, and belongs to slight crack, its access flag is arranged to 1, slight crack is added in the set of storage slight crack；Longitudinal traversal terminates
Afterwards, then to slight crack lateral traversal is carried out, it is identical as longitudinal method, and then respectively obtain the set of slight crack and hole；
S305, according to phasor picture or state found out come slight crack, hole set, loop through；Again by current slight crack
Or the coordinate of hole, minimum and maximum coordinate value is found out, determines the range of current region；
S306 judges whether the coordinate in current slight crack is labeled；If marked, perform the next step rapid；Do not mark,
Then carry out S304；
S307 carries out outside rectangle expansion to the range of current slight crack, increases N number of pixel for the first time, increase for second 2*N according to
It is secondary to analogize, current extent is traversed, checks in slight crack or hole set whether there is the model for being present in and currently being traversed
In enclosing；If it does, adjacent slight crack or hole set is added, and current adjacent slight crack number adds 1；If split before except when
Trace or hole are adjacent without other regions, or when constituting the gross area in region and accounting for less than expand area 1/3, then stop expanding
It fills, above-mentioned adjacent slight crack is stored in set, circulation terminates；
S308 successively traverses each of slight crack set slight crack or hole, until terminating, obtains the region of varying strength state
Inside aperture, the region area of slight crack and concentration.
3. the method according to claim 1, wherein the step 4 includes:
When the state of strength of the cellular is crystalline phase, the intensity of the cellular is calculated specifically:
In formula, σZFor crystalline phase intensity value, αzhFor the number of slight crack in region, αzcFor the number of region inside aperture, SzhIt is region
The area summation of interior slight crack, SzcIt is the area summation of region inside aperture, φ,It is predefined parameter related with material properties, l ×
L is cellular number；
When the state of strength of the cellular is amorphous phase, the intensity of the cellular is calculated specifically:
In formula, σaFor crystalline phase intensity value, αahFor the number of slight crack in region, αacFor the number of region inside aperture, SahIt is region
The area summation of interior slight crack, SacIt is the area summation of the region inside aperture, γ, η are predefined parameter related with material properties, l
× l is cellular number；
When the state of strength of the cellular is Kongxiang, the intensity of the cellular is calculated specifically:
In formula, σeFor Kongxiang intensity value, αehFor the number of slight crack in region, αecFor the number of region inside aperture, SehIt is in region
The area summation of slight crack, SecIt is the area summation of region inside aperture, μ, ν are predefined parameters related with material properties, and l × l is
Cellular number.
4. according to the method described in claim 3, it is characterized in that, the step 5 includes:
In formula,For overall strength, XA(t)、XZ(t)、XE(t) cellular under respectively amorphous, crystallization, three kinds of Kongxiang state
Account for the ratio of cellular sum；σZFor crystalline phase intensity value, σaFor crystalline phase intensity value, σeFor Kongxiang intensity value.
5. the method according to claim 1, wherein the method for pre-establishing model includes:
Step 21, the intensity phasor of the polymer of modeling is read；
Step 22, it is begun stepping through from the coordinate origin of the intensity phasor of modeling, judges the cellular on each pixel
State of strength；
Step 23, according to neighbour-border extension algorithm, to the pixel institute group of same intensity state of the modeling in the phasor
At region carry out hole, slight crack identification, obtain the region inside aperture of varying strength state, the region area of slight crack and intensive journey
Degree；
Step 24, the cellular of varying strength state is calculated according to the initial value of model according to the state of strength of the cellular
Intensity；
Step 25, according to the intensity of the cellular of the varying strength state, the modeling overall strength of polymer is calculated；
Step 26, modeling is compared with overall strength with the experiment value of the overall strength of the polymer of modeling, continuous mode transfer type
Parameter, optimal models are obtained, as the model built up.
6. the method according to claim 1, wherein the step 1 specifically:
The intensity phasor of polymer to be calculated is read, the pixel of the intensity phasor is L*L；L is greater than predetermined value.
7. a kind of computing device of the polymer degradation intensity based on intensity phasor characterized by comprising
Reading unit reads the intensity phasor of polymer to be calculated；
Judging unit is begun stepping through from the coordinate origin of the intensity phasor, judges the intensity shape of the cellular on each pixel
State；
Acquiring unit, according to neighbour-border extension algorithm, composed by the pixel to the same intensity state in the phasor
Region carries out hole, slight crack identification, obtains the region inside aperture of varying strength state, the region area and concentration of slight crack；
First computing unit, according to the region inside aperture of the varying strength state, the region area and concentration of slight crack, with
And the strength model built up in advance, calculate the intensity of the cellular of varying strength state；
Second computing unit calculates the overall strength of the polymer according to the intensity of the cellular of the varying strength state.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910525798.4A CN110334414B (en) | 2019-06-18 | 2019-06-18 | Method and device for calculating degradation strength of polymer based on strength phase diagram |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910525798.4A CN110334414B (en) | 2019-06-18 | 2019-06-18 | Method and device for calculating degradation strength of polymer based on strength phase diagram |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110334414A true CN110334414A (en) | 2019-10-15 |
| CN110334414B CN110334414B (en) | 2020-11-13 |
Family
ID=68142461
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910525798.4A Active CN110334414B (en) | 2019-06-18 | 2019-06-18 | Method and device for calculating degradation strength of polymer based on strength phase diagram |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110334414B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113450879A (en) * | 2021-07-07 | 2021-09-28 | 北京科技大学 | Simulation and strength calculation method for polymer degradation process |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6993446B2 (en) * | 2003-03-17 | 2006-01-31 | Schlumberger Technology Corporation | Method and apparatus for predicting the time to failure of electronic devices at high temperatures |
| CN105868563A (en) * | 2016-04-01 | 2016-08-17 | 北京科技大学 | Modeling and simulating method for degradation process of degradable high polymer material |
| CN107908921A (en) * | 2017-11-08 | 2018-04-13 | 北京科技大学 | A kind of method for simulating not homogeneity phase polymer degradation Strength Changes |
| CN108334734A (en) * | 2017-09-04 | 2018-07-27 | 北京科技大学 | A kind of high-molecular copolymer degradation modeling and simulation method |
-
2019
- 2019-06-18 CN CN201910525798.4A patent/CN110334414B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6993446B2 (en) * | 2003-03-17 | 2006-01-31 | Schlumberger Technology Corporation | Method and apparatus for predicting the time to failure of electronic devices at high temperatures |
| CN105868563A (en) * | 2016-04-01 | 2016-08-17 | 北京科技大学 | Modeling and simulating method for degradation process of degradable high polymer material |
| CN108334734A (en) * | 2017-09-04 | 2018-07-27 | 北京科技大学 | A kind of high-molecular copolymer degradation modeling and simulation method |
| CN107908921A (en) * | 2017-11-08 | 2018-04-13 | 北京科技大学 | A kind of method for simulating not homogeneity phase polymer degradation Strength Changes |
Non-Patent Citations (1)
| Title |
|---|
| 张桃红 等: "《基于多尺度的生物医用高聚物降解强度模型》", 《高分子材料科学与工程》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113450879A (en) * | 2021-07-07 | 2021-09-28 | 北京科技大学 | Simulation and strength calculation method for polymer degradation process |
| CN113450879B (en) * | 2021-07-07 | 2024-02-09 | 北京科技大学 | Polymer degradation process simulation and strength calculation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110334414B (en) | 2020-11-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Schnipke et al. | Modeling item response times with a two‐state mixture model: A new method of measuring speededness | |
| Hallett et al. | Assessing evidence for behaviour change affecting the course of HIV epidemics: a new mathematical modelling approach and application to data from Zimbabwe | |
| Burrus et al. | Person-to-person interventions targeted to parents and other caregivers to improve adolescent health: a community guide systematic review | |
| Bechtold et al. | Does mother know best? Adolescent and mother reports of impulsivity and subsequent delinquency | |
| CN105868563B (en) | A kind of modeling and simulation method of degradable high polymer material degradation process | |
| Thomas | A neurocomputational model of developmental trajectories of gifted children under a polygenic model: When are gifted children held back by poor environments? | |
| CN110334414A (en) | The calculation method and device of polymer degradation intensity based on intensity phasor | |
| CN114240874A (en) | Bone age assessment method and device based on deep convolutional neural network and feature fusion and computer readable storage medium | |
| CN110176295A (en) | A kind of real-time detecting method and its detection device of Gastrointestinal Endoscopes lower portion and lesion | |
| CN107832890A (en) | A kind of congestion Forecasting Methodology and device | |
| CN107038736A (en) | A kind of animated show method and terminal device based on frame per second | |
| CN104794321B (en) | Detection device for being detected to preceding morbid state | |
| Preston et al. | Detecting faulty within-item category functioning with the nominal response model | |
| CN110378353A (en) | A kind of tongue picture feature extracting method, system and computer readable storage medium | |
| CN114445666A (en) | Deep learning-based method and system for classifying left eye, right eye and visual field positions of fundus images | |
| Bench | A literature review evaluating parental tendencies in prior adolescent substance users | |
| CN109124669A (en) | CT data measuring method before a kind of shaping | |
| Rolka et al. | Using simulation to assess the sensitivity and specificity of a signal detection tool for multidimensional public health surveillance data | |
| Buehl et al. | Statistical problems caused by missing data resulting from neodymium: YAG laser capsulotomies in long-term posterior capsule opacification studies: problem identification and possible solutions | |
| Aurora | Academic outcomes in school-age children with stunting and non-stunting | |
| Kolodzynski et al. | Abnormal soft palate movements in patients with microtia | |
| Han et al. | Strategies for using antigen rapid diagnostic tests to reduce transmission of SARS-CoV-2 in low-and middle-income countries: a mathematical modelling study | |
| Elliott | Descriptive analysis of audiometric and psychometric scores of students at a school for the deaf | |
| CN113450879B (en) | Polymer degradation process simulation and strength calculation method | |
| Kwaadsteniet | . Ridderprint |
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 |