CN113532711B - Printing pressure identification method for central embossing cylinder of satellite type flexographic printing machine - Google Patents
Printing pressure identification method for central embossing cylinder of satellite type flexographic printing machine Download PDFInfo
- Publication number
- CN113532711B CN113532711B CN202110620780.XA CN202110620780A CN113532711B CN 113532711 B CN113532711 B CN 113532711B CN 202110620780 A CN202110620780 A CN 202110620780A CN 113532711 B CN113532711 B CN 113532711B
- Authority
- CN
- China
- Prior art keywords
- printing
- vibration
- printing pressure
- lambda
- identification
- 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.)
- Active
Links
- 238000007639 printing Methods 0.000 title claims abstract description 166
- 238000004049 embossing Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000004044 response Effects 0.000 claims abstract description 48
- 238000013499 data model Methods 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims description 29
- 230000001133 acceleration Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000001914 filtration Methods 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 16
- 238000002790 cross-validation Methods 0.000 claims description 10
- 238000010276 construction Methods 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 239000000758 substrate Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007647 flexography Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/0028—Force sensors associated with force applying means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/10—Measuring force or stress, in general by measuring variations of frequency of stressed vibrating elements, e.g. of stressed strings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/0061—Force sensors associated with industrial machines or actuators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Rotary Presses (AREA)
- Printing Methods (AREA)
Abstract
The invention provides a satellite type flexographic printing machine central embossing cylinder printing pressure identification method, firstly, discretizing modeling of a flexographic printing machine embossing system; secondly, establishing a data model of vibration response containing noise interference; thirdly, establishing a regularized printing pressure identification basic model, and realizing the identification of the coefficient load by using the l1 norm approximation, and realizing the identification of the zero load by using the l2 norm approximation; finally, the proposed method is used to identify the central embossing cylinder printing pressure of a star flexo printing machine. The invention accurately identifies the printing pressure of the printing roller, is convenient for proper setting of the printing pressure in the printing process, prevents the phenomena of printing deformation caused by overlarge printing pressure or printing blurring caused by insufficient printing pressure, and the like, and ensures the printing quality of the satellite type soft plate printer.
Description
Technical Field
The invention belongs to the technical field of load identification, and relates to a printing pressure identification method for a central embossing cylinder of a satellite type flexographic printing machine.
Background
With the full implementation of the basic national policy of saving resources and protecting environment, green printing becomes the main stream of industry development. Flexography has been strongly developed as one of green printing techniques. However, the flexibility of the printing plate makes the printing pressure critical to the quality of the printing process, and improper setting of the printing pressure will cause distortion, blurring and even color level variation of the printed pattern. The detection and control of printing pressure becomes a key indicator for evaluating printing performance.
In the printing process, the flexible printing plate and the central embossing roller form a printing pressure pair, and the printing stock is used for printing ink on the flexible printing plate on the printing stock under the action of double-sided pressure of the flexible printing plate and the central embossing roller. In this process, the flexographic plate and the central embossing cylinder are both in a rotating state, and the printing pressure varies with the variation of the substrate, making it impossible to directly measure the printing pressure between the flexographic plate and the central embossing cylinder. Therefore, how to measure or identify the printing pressure becomes a key technology for improving the quality of the flexible printing.
It is noted that this section is intended to provide a background or context for the embodiments of the invention that are recited in the claims. The description herein is not admitted to be prior art by inclusion in this section.
Disclosure of Invention
The invention aims to provide a printing pressure identification method for a central embossing cylinder of a satellite type flexographic printing machine, which can accurately identify the printing pressure of the printing cylinder, is convenient for proper setting of the printing pressure in the printing process, prevents the phenomena of printing deformation caused by overlarge printing pressure or printing blurring caused by insufficient printing pressure, and ensures the printing quality of the satellite type flexographic printing machine.
The invention adopts the following technical scheme to realize the purposes:
the method for identifying the printing pressure of the central embossing cylinder of the satellite type flexographic printing machine comprises the following steps:
s1: collecting vibration signals of a printing plate roller and a central embossing roller;
s2: constructing a vibration signal data model;
s3: carrying out convolution filtering pretreatment on the maximum correlation kurtosis of the imprinting vibration signal;
s4: building a printing pressure identification model;
s5: the printing pressure of the satellite type flexographic printing machine is identified.
Further, the step S1 specifically includes:
and acquiring vibration acceleration response data and axle center tracks of two ends of the printing plate roller and the central embossing roller in the printing process.
Further, the vibration acceleration response data of the printing plate roller means that each measuring shaft end comprises three groups of measuring data including horizontal, vertical and axial; the vibration acceleration response data of the central embossing cylinder means that each measuring shaft end comprises three groups of horizontal, vertical and axial measuring data; the axis track of the printing plate roller at least comprises two groups of axis track data in three positions of two ends and the middle of the roller; the central embossing cylinder axis track refers to two sets of axis track data comprising at least three of the cylinder ends and the cylinder middle.
Further, the step S2 specifically includes:
the vibration signal data model construction comprises discretization modeling of a printing system of a flexographic printing machine and vibration response vibration signal model construction of the printing system with noise interference;
s201: discrete modeling of the flexographic printing press impression system; in a finite time, the product of the number of discrete micro-elements and the pulse width deltat is a constant value, the number of discrete micro-elements is controlled by the size of the pulse width deltat, and the discrete form can be expressed as:
wherein R is 1…n Is the vibration response of the discretization structure of the imprinting system; h is a 1…n A vibration response transfer function of a discretized structure of the imprinting system; f (F) 0…n-1 The imprinting system discretizes the structure driving load; Δt is the pulse width;
s202, based on the vibration data measurement, constructing a data structure, and establishing a data model of vibration response containing noise interference, wherein the response basic model is as follows:
Rδ＝R+δ＝HF+δ (2)
wherein: rδ is the vibrational response of noise-containing disturbances; r is a data model of vibration response without noise interference; delta represents a measurement error in the vibrational response; h is an imprinting system characteristic matrix; f is the printing pressure.
Further, the step S3 specifically includes:
the maximum correlation kurtosis convolution filtering pretreatment of the vibration signal of the imprinting system comprises the following calculation formula:
wherein: CK (CK) M (T) is kurtosis;is a filter coefficient, T is the period of the vibration signal, and y is the signal component; m is the number of displacements.
Further, the step S4 specifically includes:
s401: for the influence of noise signals on the identification load, a regularized printing pressure identification basic model is established, and the basic model can be expressed as:
F λ,δ ＝(H T H+λI) -1 H T R δ (4)
wherein: f (F) λ,δ Printing pressures identified for regularization; h T The transpose matrix is the characteristic matrix of the imprinting system; lambda is a regularization parameter; i is an identity matrix; r is R δ A vibrational response including noise disturbance;
s402: tikhonov generalized cross-validation criterion function definition based on impression system feature matrix, and the calculation formula is as follows:
wherein: v (V) G (lambda) is a Tikhonov generalized cross-validation criterion function;tr () represents the product of the matrices;
s403: the optimal regular parameter of the Tikhonov generalized cross-validation criterion function of the imprinting system is calculated, and a calculation formula of the regular parameter can be obtained through a formula (4):
wherein: lambda (lambda) * Is the optimal regularization parameter
S404: the Tikhonov generalized cross regularized solution identifies smaller load and has poorer sparse load identification capability for more zero solutions such as impact load, pulse load and the like; therefore, on the basis of tikgonov generalized cross regularization, adopting l 1 Norm approximation identifies sparse load, l 2 The norm approximation realizes zero load identification, and the Tikhonov generalized cross regularization model can be written as:
wherein:is l 1 +l 2 Norm approximation Tikhonov generalized cross regularization parameter lambda 2 Is l 2 Norms regular parameters; lambda (lambda) 1 Is l 1 Normative parameters->Penalty term for l2, F 1 Penalty term for l 1;
s405: solving for l 1 +l 2 Norm regularization problem parameter setting, i.e. l as shown in formula (6) 1 +l 2 Norm regularization form, determining λ from a priori information about sparse loading 1 And lambda (lambda) 2 ，λ 1 Can take lambda 1 ∈[0.001λ m ，0.1λ m ]，λ m ＝||2H T R||∞；λ 2 ∈[0.1λ * ,λ * ]。
Further, the step S5 specifically includes:
will l 1 +l 2 And (4) introducing a norm regularization parameter and Tikhonov generalized cross regularization into the formula (4) to obtain the printing pressure of the satellite type flexographic printing machine.
The invention has the beneficial effects that:
(1) The invention discretizes vibration response and transfer function aiming at discretization modeling of the imprinting system of the flexographic printing machine, refines system structure and comprehensively characterizes dynamic characteristics of the imprinting system of the printing machine;
(2) The invention establishes a data model of vibration response containing noise interference, fuses the multiple data such as radial vibration, axial vibration track and the like of the printing plate roller and the embossing roller, and forms a vibration data structure of the embossing system of the full-scale flexographic printing machine;
(3) The invention adopts the convolution filtering pretreatment filtering with the maximum correlation kurtosis, eliminates the background noise of the machine operation in the vibration signal acquisition process, reserves the vibration response of each driving load of the system operation, and prevents the influence of noise signals on the printing pressure identification structure;
(4) Aiming at the regularized printing pressure identification basic model, the invention solves the influence of the singular value of the reverse problem on the identification result, converts the random unfixed problem in the identification process into the proper problem, and avoids the identification error;
(5) The invention realizes the identification of the coefficient load by using the l1 norm approximation, and realizes the identification of the zero load by using the l2 norm approximation.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a flow chart of the present invention;
FIG. 3 is a flow chart of the filtering of the vibration signal of the roller according to the present invention.
Detailed Description
Example embodiments will now be described more fully with reference to the accompanying drawings. However, the exemplary embodiments may be embodied in many forms and should not be construed as limited to the examples set forth herein; rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of the example embodiments to those skilled in the art. The described features or characteristics may be combined in any suitable manner in one or more embodiments.
The invention will be described in detail with reference to the accompanying drawings, wherein the specific steps are as follows:
as shown in fig. 1 and 2, the printing pressure identification method comprises five major parts of vibration signal acquisition, vibration signal data model construction, maximum correlation kurtosis convolution filtering treatment, printing pressure identification model construction and printing pressure calculation of a printing plate roller and a central embossing roller.
Collecting vibration signals of the printing plate roller and the central embossing roller, and collecting vibration response of the printing plate roller and the central embossing roller under the driving of printing pressure; the vibration signal data model construction is to discretize modeling of the imprinting system structure of the flexographic printing machine, and the imprinting system vibration response vibration signal model construction of the model containing noise interference is realized according to the discretization modeling; the maximum correlation kurtosis convolution filtering processing is the filtering processing of noise signals and vibration signals driven by non-printing pressure in the vibration signals; the printing pressure identification model is constructed by adopting an imprinting characteristic matrix Tikhonov generalized cross-validation criterion function and l on the basis of a regularized printing pressure identification basic model 1 +l 2 Establishing a printing pressure recognition theoretical model by approximating a Tikhonov generalized cross regularization parameter by a norm; and calculating the printing pressure, and identifying the printing pressure by adopting the processing method.
The method comprises the steps of collecting vibration signals of a printing plate roller and a central embossing roller, and carrying out live measurement on test data of at least three printing standard working conditions; the vibration acceleration response measurement of the printing plate roller comprises vibration acceleration response data of two ends of the printing plate roller and the central embossing roller; the experimental measurement of the axle center track comprises the axle center track measurement of the two ends of the printing plate roller and the central embossing roller; the vibration acceleration response measurement of the printing plate roller comprises 3 groups of horizontal, vertical and axial measurement data at each measurement shaft end; a central embossing cylinder vibration acceleration response measurement, each measurement shaft end containing 3 sets of horizontal, vertical and axial measurement data; the axle center track measurement of the printing plate roller at least comprises two groups of axle center track data at the two ends of the roller and the middle 3 of the roller; the central embossing cylinder axis locus measurement comprises at least two sets of axis locus data at the two ends of the cylinder and at the middle 3 of the cylinder.
As shown in FIG. 3, the maximum relevant kurtosis convolution filtering process of the vibration signal of the flexographic printing imprint system comprises convolution calculation initialization, vibration signal convolution calculation, relevant kurtosis value calculation, optimal filter vector calculation and deconvolution calculation.
Wherein, the convolution calculation is initialized, and the length of the filter is set to be 350; and judging the kurtosis value, wherein the threshold value is 3.
The invention provides a specific implementation condition of a printing pressure identification method of a printing cylinder of a flexographic printing machine, which comprises the following steps:
during printing in a flexographic printing machine, the plate cylinder and the central embossing cylinder are in contact with the inner and outer sides of the substrate, respectively, and there is a printing pressure between the plate cylinder, the substrate, and the central embossing cylinder. As the printing amount increases, the printing pressure may vary with the fluctuation of the printing amount, and the printing pressure may be excessively large to cause printing deformation, or insufficient to cause printing blurring, or the like. The printing pressure identification method for the embossing cylinder of the flexographic printing machine detects the shaft end vibration acceleration and the shaft center track of the printing plate cylinder and the central embossing cylinder in real time, models the embossing system of the flexographic printing machine in a discretization mode, discretizes the vibration response and the transfer function, and comprehensively characterizes the dynamic characteristics of the embossing system of the printing machine; further, the multi-element data such as radial vibration, axial vibration track and the like of the printing plate roller and the embossing roller are fused to form a vibration data structure of the embossing system of the full-scale flexographic printing machine, and the maximum relevant kurtosis convolution filtering preprocessing filtering is adopted to eliminate the background noise of the machine operation in the vibration signal acquisition process, retain the vibration response of each driving load of the system operation and prevent the influence of noise signals on the printing pressure identification structure; meanwhile, a regularized printing pressure identification basic model is established, the influence of a reverse problem singular value on an identification result is solved, the random unfixed problem in the identification process is converted into the unfixed problem, the identification error is avoided, the identification of the coefficient load is realized by utilizing the l1 norm approximation, and the zero load is identified by utilizing the l2 norm approximation; and finally, the accurate identification of the printing pressure of the central embossing cylinder of the satellite type flexographic printing machine is realized.
The identification method of the printing pressure of the embossing roller of the satellite type flexographic printing machine is mainly realized by the following steps:
and step 1, measuring corresponding axle center tracks by vibration acceleration of a central embossing cylinder and a printing plate cylinder in the printing process.
The vibration acceleration response measurement of the printing plate roller comprises vibration acceleration response data of two ends of the printing plate roller and the central embossing roller; the experimental measurement of the axle center track comprises the axle center track measurement of the two ends of the printing plate roller and the central embossing roller; the vibration acceleration response measurement of the printing plate roller comprises 3 groups of horizontal, vertical and axial measurement data at each measurement shaft end; the vibration acceleration response measurement of the central embossing cylinder comprises three groups of horizontal, vertical and axial measurement data at each measurement shaft end; the axle center track measurement of the printing plate roller at least comprises two groups of axle center track data at the two ends of the roller and the middle 3 of the roller; the central embossing cylinder axis track measurement at least comprises two groups of axis track data at the two ends of the cylinder and at the middle 3 of the cylinder;
step 2, discretization modeling of the imprinting system of the flexographic printing machine, wherein in a finite time, the product of the number of discrete microelements and the pulse width delta t is a constant value, the number of the discrete microelements can be controlled by the size of the pulse width delta t, and the discrete form can be expressed as:
wherein R is 1…n Is the vibration response of the discretization structure of the imprinting system; h is a 1…n A vibration response transfer function of a discretized structure of the imprinting system; f (F) 0…n-1 The imprinting system discretizes the structure driving load; Δt is the pulse width.
Step 3, based on the vibration data measurement, constructing a data structure, and establishing a data model of vibration response containing noise interference, wherein the response basic model is as follows:
R δ ＝R+δ＝HF+δ (2)
wherein: r is R δ A vibrational response that is noisy; r is a data model of vibration response without noise interference; delta represents a measurement error in the vibrational response; h is an imprinting system characteristic matrix; f is the printing pressure;
and 4, carrying out convolution filtering pretreatment on the maximum correlation kurtosis of the vibration signal of the imprinting system, wherein the calculation formula of the filtering pretreatment is as follows:
wherein: CK (CK) M (T) is kurtosis;is a filter coefficient, T is the period of the vibration signal, and y is the signal component; m is the displacement number;
step 5, establishing a regularized printing pressure identification basic model aiming at the influence of noise signals on identification load, wherein the basic model can be expressed as:
F λ,δ ＝(H T H+λI) -1 H T R δ (4)
wherein: f (F) λ,δ Printing pressures identified for regularization; h T The transpose matrix is the characteristic matrix of the imprinting system; lambda is a regularization parameter; i is an identity matrix; r is R δ A vibrational response including noise disturbance;
step 6, tikhonov generalized cross-validation criterion function definition based on an imprinting system feature matrix, wherein a calculation formula is as follows:
wherein: v (V) G (lambda) is a Tikhonov generalized cross-validation criterion function;tr () represents the product of the matrices;
step 7, calculating optimal regular parameters of Tikhonov generalized cross-validation criterion functions of the imprinting system, wherein a calculation formula of the regular parameters can be obtained through a formula (4):
wherein: lambda (lambda) * Is the optimal regularization parameter
The solution of step 8, tikhonov generalized cross regularization is more prone to identify smaller loads, and the sparse load identification capability of more zero solutions such as impact and pulse loads is poorer. Therefore, on the basis of tikgonov generalized cross regularization, adopting l 1 Norm approximation identifies sparse load, l 2 The norm approximation realizes zero load identification, and the Tikhonov generalized cross regularization model can be written as:
wherein:is l 1 +l 2 Norm approximation Tikhonov generalized cross regularization parameter lambda 2 Is l 2 Norms regular parameters; lambda (lambda) 1 Is l 1 Normative parameters->Penalty term for l2, F 1 Penalty term for l 1;
step 9, solving for l 1 +l 2 Norm regularization problem parameter setting, i.e. l as shown in formula (6) 1 +l 2 Norm regularization form, determining λ from a priori information about sparse loading 1 And lambda (lambda) 2 ，λ 1 Can take lambda 1 ∈[0.001λ m ，0.1λ m ]，λ m ＝||2H T R||∞；λ 2 ∈[0.1λ * ,λ * ]。
And step 10, introducing the l1+l2 norm regularization parameter and Tikhonov generalized cross regularization into a formula (4) to obtain the printing pressure of the satellite type flexographic printing machine.
The satellite type flexographic printing machine embossing cylinder printing pressure identification can be obtained through the steps (1) to (10). If there is a deviation, the filtering performance can be modified by changing the filter length in the step (4); adjusting the identification of the l1 norm approximation modification coefficient load; the adjustment of the l2 norm approximation modifies the identification of the zero load. Likewise, the adjustment of the recognition performance can be completed by using the relation between the recognition parameters and the recognition accuracy in the steps (2) to (9).
Other embodiments of the invention will be apparent to those skilled in the art from consideration of the specification and practice of the invention disclosed herein. This application is intended to cover any variations, uses, or adaptations of the invention following, in general, the principles of the invention and including such departures from the present disclosure as come within known or customary practice within the art to which the invention pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the invention being indicated by the following claims.
Claims (4)
1. A method for identifying the printing pressure of a central embossing cylinder of a satellite-type flexographic printing machine, comprising the steps of:
s1: collecting vibration signals of a printing plate roller and a central embossing roller;
s2: constructing a vibration signal data model;
the step S2 specifically comprises the following steps:
the vibration signal data model construction comprises discretization modeling of a printing system of a flexographic printing machine and vibration response vibration signal model construction of the printing system with noise interference;
s201: discrete modeling of the flexographic printing press impression system; in a finite time, the product of the number of discrete micro-elements and the pulse width deltat is a constant value, the number of discrete micro-elements is controlled by the size of the pulse width deltat, and the discrete form can be expressed as:
wherein R is 1…n Is the vibration response of the discretization structure of the imprinting system; h is a 1…n Is pressed intoA printing system discretization structure vibration response transfer function; f (F) 0…n-1 The imprinting system discretizes the structure driving load; Δt is the pulse width;
s202, based on the vibration data measurement, constructing a data structure, and establishing a data model of vibration response containing noise interference, wherein the response basic model is as follows:
Rδ＝R+δ＝HF+δ (2)
wherein: rδ is the vibrational response of noise-containing disturbances; r is a data model of vibration response without noise interference; delta represents a measurement error in the vibrational response; h is an imprinting system characteristic matrix; f is the printing pressure;
s3: carrying out convolution filtering pretreatment on the maximum correlation kurtosis of the imprinting vibration signal;
the step S3 specifically comprises the following steps:
the maximum correlation kurtosis convolution filtering pretreatment of the vibration signal of the imprinting system comprises the following calculation formula:
wherein: CK (CK) M (T) is kurtosis;is a filter coefficient, T is the period of the vibration signal, and y is the signal component; m is the displacement number;
s4: building a printing pressure identification model;
the step S4 specifically includes:
s401: for the influence of noise signals on the identification load, a regularized printing pressure identification basic model is established, and the basic model can be expressed as:
F λ,δ ＝(H T H+λI) -1 H T R δ (4)
wherein: f (F) λ,δ Printing pressures identified for regularization; h T The transpose matrix is the characteristic matrix of the imprinting system; lambda is a regularization parameter; i is an identity matrix; r is R δ Vibration response with noise disturbance；
S402: tikhonov generalized cross-validation criterion function definition based on impression system feature matrix, and the calculation formula is as follows:
wherein: v (V) G (lambda) is a Tikhonov generalized cross-validation criterion function;
tr () represents the product of the matrices;
s403: the optimal regular parameter of the Tikhonov generalized cross-validation criterion function of the imprinting system is calculated, and a calculation formula of the regular parameter can be obtained through a formula (4):
wherein: lambda (lambda) * Is the optimal regularization parameter
S404: the Tikhonov generalized cross regularized solution identifies smaller load and has poorer sparse load identification capability for more zero solutions such as impact load, pulse load and the like; therefore, on the basis of tikgonov generalized cross regularization, adopting l 1 Norm approximation identifies sparse load, l 2 The norm approximation realizes zero load identification, and the Tikhonov generalized cross regularization model can be written as:
wherein:is l 1 +l 2 Norm approximation Tikhonov generalized cross regularization parameter lambda 2 Is l 2 Norms regular parameters; lambda (lambda) 1 Is l 1 Normative parameters->Penalty term for l2, F 1 Penalty term for l 1;
s405: solving for l 1 +l 2 Norm regularization problem parameter setting, i.e. l as shown in formula (6) 1 +l 2 Norm regularization form, determining λ from a priori information about sparse loading 1 And lambda (lambda) 2 ，λ 1 Can take lambda 1 ∈[0.001λ m ，0.1λ m ]，λ m ＝||2H T R||∞；λ 2 ∈[0.1λ * ,λ * ]；
S5: the printing pressure of the satellite type flexographic printing machine is identified.
2. The method for identifying the printing pressure of a central embossing cylinder of a satellite flexographic printing machine according to claim 1, characterized in that said step S1 is in particular:
and acquiring vibration acceleration response data and axle center tracks of two ends of the printing plate roller and the central embossing roller in the printing process.
3. A satellite flexographic printing press central embossing cylinder printing pressure identification method according to claim 2, wherein:
the vibration acceleration response data of the printing plate roller means that each measuring shaft end comprises three groups of horizontal, vertical and axial measuring data; the vibration acceleration response data of the central embossing cylinder means that each measuring shaft end comprises three groups of horizontal, vertical and axial measuring data; the axis track of the printing plate roller at least comprises two groups of axis track data in three positions of two ends and the middle of the roller; the central embossing cylinder axis track refers to two sets of axis track data comprising at least three of the cylinder ends and the cylinder middle.
4. A satellite flexographic printing press central embossing cylinder printing pressure identification method according to claim 3, characterized in that said step S5 is in particular:
will l 1 +l 2 And (4) introducing a norm regularization parameter and Tikhonov generalized cross regularization into the formula (4) to obtain the printing pressure of the satellite type flexographic printing machine.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110620780.XA CN113532711B (en) | 2021-06-03 | 2021-06-03 | Printing pressure identification method for central embossing cylinder of satellite type flexographic printing machine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110620780.XA CN113532711B (en) | 2021-06-03 | 2021-06-03 | Printing pressure identification method for central embossing cylinder of satellite type flexographic printing machine |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113532711A CN113532711A (en) | 2021-10-22 |
| CN113532711B true CN113532711B (en) | 2024-04-05 |
Family
ID=78095599
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110620780.XA Active CN113532711B (en) | 2021-06-03 | 2021-06-03 | Printing pressure identification method for central embossing cylinder of satellite type flexographic printing machine |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113532711B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004306378A (en) * | 2003-04-04 | 2004-11-04 | Seiko Epson Corp | Control method for ink-jet printer, and ink-jet printer |
| CN104866671A (en) * | 2015-05-26 | 2015-08-26 | 山东大学 | Vibration speed signal-based method for recognizing maximal pressure rise rate |
| CN105403242A (en) * | 2015-05-22 | 2016-03-16 | 华东交通大学 | Locomotive pantograph-catenary hard point photoelectric vibration comprehensive detection and GPS positioning method and system |
| CN105844250A (en) * | 2016-03-31 | 2016-08-10 | 山东大学 | Method for identifying maximum pressure rising rate based on vibration acceleration signal |
| CN106951197A (en) * | 2017-01-19 | 2017-07-14 | 北京印刷学院 | A kind of printing equipment remote fault diagnosis method and system based on signal testing |
| CN107073923A (en) * | 2014-10-16 | 2017-08-18 | 温德默勒&霍乐沙两合公司 | Body for receiving printing picture structure |
| CN107729612A (en) * | 2017-09-18 | 2018-02-23 | 北京印刷学院 | A kind of hidden danger appraisal procedure for printing equipment critical component |
| CN107748891A (en) * | 2017-09-18 | 2018-03-02 | 北京印刷学院 | A kind of printing equipment critical component state identification method based on normal operation in normal domain |
| CN108920850A (en) * | 2018-07-09 | 2018-11-30 | 西安理工大学 | A kind of flexo pressure prediction method based on convolutional neural networks |
-
2021
- 2021-06-03 CN CN202110620780.XA patent/CN113532711B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004306378A (en) * | 2003-04-04 | 2004-11-04 | Seiko Epson Corp | Control method for ink-jet printer, and ink-jet printer |
| CN107073923A (en) * | 2014-10-16 | 2017-08-18 | 温德默勒&霍乐沙两合公司 | Body for receiving printing picture structure |
| CN105403242A (en) * | 2015-05-22 | 2016-03-16 | 华东交通大学 | Locomotive pantograph-catenary hard point photoelectric vibration comprehensive detection and GPS positioning method and system |
| CN104866671A (en) * | 2015-05-26 | 2015-08-26 | 山东大学 | Vibration speed signal-based method for recognizing maximal pressure rise rate |
| CN105844250A (en) * | 2016-03-31 | 2016-08-10 | 山东大学 | Method for identifying maximum pressure rising rate based on vibration acceleration signal |
| CN106951197A (en) * | 2017-01-19 | 2017-07-14 | 北京印刷学院 | A kind of printing equipment remote fault diagnosis method and system based on signal testing |
| CN107729612A (en) * | 2017-09-18 | 2018-02-23 | 北京印刷学院 | A kind of hidden danger appraisal procedure for printing equipment critical component |
| CN107748891A (en) * | 2017-09-18 | 2018-03-02 | 北京印刷学院 | A kind of printing equipment critical component state identification method based on normal operation in normal domain |
| CN108920850A (en) * | 2018-07-09 | 2018-11-30 | 西安理工大学 | A kind of flexo pressure prediction method based on convolutional neural networks |
Non-Patent Citations (2)
| Title |
|---|
| 基于AR模型的印刷机滚筒扭矩及其振动试验研究;王仪明;赵明明;武淑琴;许文才;;振动与冲击;20160215;35(3);全文 * |
| 基于多元图像特征分析的印刷单元故障诊断方法;徐卓飞;张海燕;宋希龙;翟一红;吴欣阳;;包装工程;20150410(07);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113532711A (en) | 2021-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5395242B2 (en) | Method for detecting occurrence of a printing error on a printing medium during processing of the printing medium on a printing press | |
| CN109791094B (en) | Method and system for efficient load identification | |
| GB2108264A (en) | Register control in the production of multi-colour printed sheets | |
| CN113532711B (en) | Printing pressure identification method for central embossing cylinder of satellite type flexographic printing machine | |
| EP1958775B1 (en) | Ink splitting correction method | |
| CN106739485A (en) | A kind of printing machine longitudinal direction alignment on-line checking and method for diagnosing faults and device | |
| CN112765881B (en) | Dynamic weighing method and system capable of being expanded to multiple bridges based on neural network | |
| US20020050219A1 (en) | Method and apparatus for setting register on a multicolor printing machine | |
| CN109476151B (en) | Production facility with inspection system and method for adjusting the sensitivity of an inspection system | |
| CN113128338A (en) | Intelligent diagnosis method for printing machine roller fault under small sample | |
| US20070051265A1 (en) | Method for printing correction | |
| DE10139310A1 (en) | Method and imaging device for register setting | |
| CN101332698B (en) | Metering zone for determining lipid-sensing limitation during printing | |
| CN110667258B (en) | Method for analyzing printing quality by means of neural network | |
| Durnyak et al. | Determination of the Optical Density of Two-Parameter Tone Transfer for a Short Printing System of the Sixth Dimension. | |
| CN113156832A (en) | Method and system for rapidly analyzing performance parameters of ink path system of offset press | |
| DE102018200994B4 (en) | Puncture register | |
| US6318165B1 (en) | Device for measuring the stickiness of a flowable medium | |
| CN207942800U (en) | A kind of printing plate cylinder mechanism of printing machine | |
| JP2002137529A (en) | Method for suppressing overfilling or underfilling in production of printing original image to minimum | |
| CN113320286A (en) | Optimizing detection mode of waste page | |
| AT524837B1 (en) | SYSTEM OF MEASUREMENT SYSTEM AND RIM | |
| JPH029636A (en) | Flexographic rotary press | |
| DE102013017055A1 (en) | Printing plate for offset printing unit, has surface with first surface area, which is color-friendly, second surface area, which is ink-repellent and third surface area, which is color-friendly | |
| EP3984749A1 (en) | Method for producing printed products |
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 |