CN108828132B

CN108828132B - Method for analyzing content of main components in cigarette smoke

Info

Publication number: CN108828132B
Application number: CN201810653885.3A
Authority: CN
Inventors: 吴晓炯; 罗佳; 汪旭
Original assignee: Shanghai Tobacco Group Co Ltd
Current assignee: Shanghai Tobacco Group Co Ltd
Priority date: 2018-06-22
Filing date: 2018-06-22
Publication date: 2021-02-26
Anticipated expiration: 2038-06-22
Also published as: CN108828132A

Abstract

The invention provides a method for analyzing the content of main components in cigarette smoke, which comprises the following steps: step S1, determining basic parameters of the cigarette; step S2, calculating the original parameters of the cigarette according to the basic parameters of the cigarette in the step S1; step S3, calculating the relative parameters of the cigarette according to the basic parameters of the cigarette in the step S1 and the original parameters of the cigarette in the step S2; and step S4, calculating to obtain the main component content value in the cigarette smoke according to the original parameters of the cigarette in the step S2 and the relative parameters of the cigarette in the step S3. The method for analyzing the content of the main components in the cigarette smoke can be applied to a computer, can calculate the influence of the change of the formula, the auxiliary materials and the cigarette structure on the main stream smoke, can accurately calculate the content of the main components in the cigarette smoke, does not need to modify the existing cigarette production and processing equipment, does not need to carry out a large number of repeated early-stage pre-tests, and has the advantages of convenience and quickness.

Description

Method for analyzing content of main components in cigarette smoke

Technical Field

The invention belongs to the technical field of cigarettes, and relates to a method for analyzing the content of main components in cigarette smoke.

Background

Since the 50 s of the last century, a great deal of related research has been carried out at home and abroad in order to improve the quality of cigarette products, accelerate the design speed of the products and optimize the experience of the products. Foreign cigarette enterprises mostly start from the first physical and chemical rules, and carry out related researches in subdivided fields such as cigarette paper ventilation rules, filter stick filtering capacity, cigarette combustion mechanism and the like, so that the overall understanding of the cigarette combustion and filtration mechanism is formed; domestic research is generally based on statistical tests, and the combination of some specific auxiliary materials and rolling parameters is analyzed to achieve a certain prediction function of the mainstream smoke. However, because of too many factors influencing the smoke of the cigarettes, the two factors do not combine the combustion rules of the whole cigarettes into a whole, and an analysis method based on the whole parameters of the cigarettes is not formed.

Specifically, the composition and content of cigarette smoke are influenced by tobacco shred components, cigarette structure, cigarette material parameters, combustion conditions and other factors, and relate to a plurality of links such as tobacco leaf production, tobacco shred processing, rolling equipment, product design and the like. In order to accurately predict the mainstream smoke of cigarettes and achieve the expected target of cigarette design, a great deal of research work is done in the industry. CN105628646 discloses an on-line cigarette tar prediction and early warning method, which utilizes infrared light to collect tobacco shred components, constructs a prediction model of tar content, and forms an early warning threshold value by calculating upper and lower limit values through factor analysis. However, the method is significantly influenced by the environment, and the measurement precision deviation is large in the data process. CN106108107 discloses a modeling design method of cigarette total ventilation rate, filter ventilation rate and draw resistance based on cigarette structure and physical parameters, which takes cigarette structure characteristic parameters and cigarette physical characteristic parameters as variables, takes cigarette total ventilation rate, filter ventilation rate and draw resistance as objective functions, and obtains an optimization method meeting the cigarette design requirements by adjusting the cigarette structure and physical parameters. However, the method only obtains the ventilation rate and the suction resistance of the cigarette, and cannot deeply research the content of the main components of the cigarette smoke.

Therefore, the method integrates indexes used in cigarette design such as a leaf group formula, a rolling auxiliary material and a rolling structure, forms a physical method based on the first property and has great value for the cigarette design, and the method has practical significance for reducing the enterprise cost. In addition, the cigarette smoke composition and content are researched by adopting a single-factor experimental design method, and the design requirements of the cigarette smoke composition and content cannot be well met generally, because the tobacco leaf composition formula of the cigarette and the cigarette structure parameters are interacted, the cigarette product is determined together. Therefore, it is necessary to construct an analysis method for the smoke composition and content of cigarettes based on the two aspects of the leaf group formula and the rolling structure.

Disclosure of Invention

In view of the above disadvantages of the prior art, the present invention provides an analysis method for the content of main components in cigarette smoke, which can effectively and accurately analyze the content of main components in cigarette smoke based on an algorithm model of basic parameters and related parameters of cigarettes, so as to further adjust the formula of tobacco leaf group and the rolling structure, and improve the design efficiency of cigarettes.

In order to achieve the above and other related objects, the present invention provides a method for analyzing the content of main components in cigarette smoke, comprising the steps of:

step S1, determining basic parameters of the cigarette;

preferably, in step S1, the basic parameters of the cigarette include, but are not limited to, tobacco type, tobacco weight, tobacco segment length, cigarette circumference, filter length, filter pressure drop, tipping paper air permeability, total tipping paper pore band width, tipping paper outer pore band spacing, forming paper air permeability, cigarette paper combustion improver content, cigarette paper gram weight, nicotine filtration coefficient, and tar filtration coefficient. The basic parameters of the cigarette are conventional data of the cigarette and can be measured and obtained. Step S2, calculating the original parameters of the cigarette according to the basic parameters of the cigarette in the step S1 by adopting a smoldering time recurrence method;

preferably, in step S2, the original parameters of the cigarette include, but are not limited to, original smoldering time, original number of mouths, original length of burn per mouth, original weight of burn per mouth, and original remaining length of burn per mouth.

Preferably, in step S2, the calculation of the original parameters of the cigarette includes the following steps:

step S21, after determining the type of the tobacco shreds, calculating the original smoldering time according to a formula (1), wherein the formula (1) is as follows:

t1＝73.3×l1+279×c-33.1×a-19.5×mc+6.86×ps+11.7×pe+6.49×pp+10.4×pl，

wherein t1 is the original smoldering time, s; l1 is the length of cigarette, cm; c is the cigarette circumference, cm; a is the content of the combustion improver of the cigarette paper in percent; mc is the gram weight of cigarette paper, g/m²(ii) a ps is the ratio of cut stems,%; pe is the proportion of the expanded filaments,%; pp is flake ratio,%; pl is the leaf-silk ratio,%.

The cut stem proportion, the cut expansion proportion, the slice proportion and the cut leaf proportion are set numerical values of cigarettes and all belong to fixed data parameter values which can be obtained by designers.

Step S22, calculating the original opening number according to the formula (2) by the original smoldering time, wherein the formula (2) is as follows:

n＝0.0121×t1+0.487，

wherein n is the original number of ports; t1 is the original smoldering time, s.

Step S23, calculating the original length of each port combustion according to the formulas (3), (4) and (5) through the original smoldering time and the original port number, wherein the formula (3) is as follows:

v＝l1/t1，

wherein v is smoldering speed, cm/s; l1 is the length of cigarette, cm; t1 is the original smoldering time, s;

the formula (4) is:

t2＝58×n，

wherein t2 is the actual smoldering time, s; n is the original number of ports;

the formula (5) is:

l2＝v×(t1-t2)/n，

wherein l2 is the original length of each port combustion, cm; v is smoldering speed, cm/s; t1 is the original smoldering time, s; t2 is the actual smoldering time, s; n is the original number of ports.

Step S24, calculating the original weight of each port of combustion according to a formula (6) through the original length of each port of combustion, wherein the formula (6) is as follows: m2 ═ m1 xl 2/l1,

wherein m2 is the original weight per burner, g; m1 is the weight of tobacco shred, g; l2 is the original length of each port burning, cm; l1 is the length of the tobacco rod in cm.

More preferably, the initial weight of each combustion is calculated according to the formulas (a) and (b) respectively as the initial calculated weight of nicotine and the initial calculated weight of tar,

the formula (a) is: m3 is m2 xnic,

wherein m3 is the original calculated weight of nicotine, g; m2 is the original weight per port, g; nic is the nicotine experience coefficient of the tobacco shred and is set to be 1.11. The nicotine experience coefficient nic of the tobacco shreds is a coefficient for adjusting the nicotine release amount of different leaf groups;

the formula (b) is: m4 ═ m2 × tar,

wherein m4 is the original calculated weight of tar, g; m2 is the original weight per port, g; tar is the empirical coefficient of tar in tobacco shreds, and is set to 1.22. The tar empirical coefficient tar is a coefficient for adjusting the tar release amount of different leaf groups.

Step S25, calculating the original residual length of each port of combustion according to formulas (7) and (8) through the original length of each port of combustion, wherein the formula (7) is as follows: l3[1] ═ l1-l2,

wherein l3[1] is the original residual length of combustion of each port 1, cm; l1 is the length of cigarette, cm; l2 is the original length of each port burning, cm;

the formula (8) is: l3[ i ] ═ l3[ i-1] -58v-l2,

wherein l3[ i ] is the original residual length of combustion of each port of the ith port, and i is 2, …, n, cm; l3[ i-1] is the original residual length of combustion of each port i-1, i is 2, …, n, cm; v is smoldering speed, cm/s; l2 is the original length of each port burned in cm.

The formulas (7) and (8) are iterative algorithms.

Preferably, in step S2, the smoldering time recurrence method is to calculate the original number of burned cigarettes through the original smoldering time of the cigarettes (only related to the type of tobacco, the rolling structure, and the cigarette paper), further calculate the original burning length and the original burning weight of each cigarette during burning, further obtain the original calculated weights of nicotine and tar of each cigarette (assuming that the tobacco of the cigarettes is uniformly distributed), and calculate to obtain the original remaining burning length of each cigarette.

Preferably, in step S2, the original parameters of the cigarette are parameters calculated by ignoring the perforation condition of the tipping paper. The tipping paper punching condition refers to the condition that the holes of the tipping paper are blocked by adhesive tapes.

Step S3, calculating relative parameters of the cigarette according to the basic parameters of the cigarette in the step S1 and the original parameters of the cigarette in the step S2 by adopting a pressure drop matching regression algorithm;

preferably, in step S3, the relative parameters of the cigarette include, but are not limited to, relative number of ports, relative remaining length of combustion per port, and relative weight of combustion per port.

Preferably, in step S3, the relative parameters of the cigarette are parameters calculated under the conditions of perforation of the tipping paper.

Preferably, in step S3, the calculation of the relative parameters of the cigarette includes the following steps:

step S31, calculating the relative port number according to the formula (9), wherein the formula (9) is as follows:

rn＝l1/(l6avg+58v)，

wherein rn is the relative number of ports; l1 is the length of cigarette, cm; l6avg is the average relative burn length, cm; v is smoldering speed, cm/s.

More preferably, the average relative burn length is calculated according to formula (c) which is:

l6avg＝Simga(l6[i])/n，

wherein l6avg is the average relative burning length, cm; simga is the running sign Sigma, which is the sum of the logarithmic group l6[ i ]; l6[ i ] is the relative combustion length after ventilation of the ith port, i is 1, …, n, cm; n is the original number of ports.

Further preferably, the relative combustion length after the ith port ventilation is calculated according to the formula (d): l6[ i ] ═ qi [ i ]/17.5 × l2,

wherein l6[ i ] is the relative combustion length after the ith port is ventilated, i is 1, …, n, cm; qi [ i ] is the cigarette combustion inlet flow rate qi value at the ith port, i is 1, …, n, ml/s; l2 is the original length of each port burned in cm.

Still more preferably, said cigarette combustion inlet flow rate qi value at the ith port is calculated according to the formula (e): qi [ i ] ═ 0.9 x (17.5-q [ i ]),

wherein qi is the cigarette burning inlet flow rate qi value of the ith port, i is 1, …, n, ml/s; q [ i ] is the ventilation flow rate q value of the filter at the length l3[ i ] of the ith port, i is 1, …, n, ml/s.

Still further preferably, the cigarette combustion inlet flow rate of the ith port conforms to the flow rate identity in the formula (f): qi [ i ] + qc [ i ] + q [ i ] ═ 17.5,

wherein qi is the cigarette burning inlet flow rate qi value of the ith port, i is 1, …, n, ml/s; qc [ i ] is the flow rate qc of the wrapping paper at the ith port, i is 1, …, n, ml/s; q [ i ] is the ventilation flow rate q value of the filter at the length l3[ i ] of the ith port, i is 1, …, n, ml/s.

Still further preferably, the flow rate of the filter is calculated according to a regression pressure drop matching equation in formula (g), wherein the formula (g) is:

q＝FindRoot[ptpw[q]＝pr[q]+pre[q]+pf[q]，q->0]，

wherein q is the flow rate through the filter tip, ml/s; findroot is the initial value (q->0) Searching the numerical solution of an equation, ml/s; q->0 means that q tends towards 0; ptpw [ q ]]Pressure drop function, cmH, for a combination of tipping paper and forming paper at a flow rate q₂O；pr[q]Considering for a segment the pressure drop function, cmH, of the air permeability of the cigarette paper at a flow rate q₂O；pre[q]Is a function of the pressure drop, cmH, of the cigarette in the cross section at a flow rate of q₂O；pf[q]Is a function of pressure drop of the filter rod at a flow rate of q, cmH₂O。

Particularly preferably, the pressure drop function of the combination of the tipping paper and the forming paper at the flow velocity q is calculated according to a formula (h), the pressure drop function of the cigarette section under the ventilation condition of the cigarette paper is calculated according to a formula (i) by considering the pressure drop function of the air permeability of the cigarette paper at the flow velocity q, the pressure drop function of the cigarette at the cross section with the flow velocity q is calculated according to a formula (j) by calculating the pressure drop function of the cigarette under the ventilation condition of no cigarette paper, the pressure drop function of the filter stick with the flow velocity q is calculated according to a formula (k),

the formula (h) is: ptpw [ q ]]＝(q/c)²×(2.56×10⁸/a2²+3.6×10⁴/a3²)，

Wherein, ptpw [ q ]]Pressure drop function, cmH, for a combination of tipping paper and forming paper at a flow rate q₂O; q is the flow rate through the filter tip, ml/s; c is the cigarette circumference, cm; a2 is the forming paper air permeability, CU; a3 is the tipping paper air permeability, CU;

the formula (i) is:

pr[q]＝pre0×Tanh(Sqrt(a1×c×l3[i]×pre0×(7×10^-5+2.9×10^-6×(17.5-q)))/Sqrt(a1×c×l3[i]×pre0×(7×10^-5+2.9×10^-6×(17.5-q)))，

wherein, pr [ q ]]For cigarettes with considerationPressure drop function of paper permeability at flow rate q, cmH₂O; pre0 is the pressure drop value, cmH, of cigarette segment neglecting the air permeability of cigarette paper₂O; tanh is a hyperbolic tangent function; sqrt is a root function; a1 is the air permeability of cigarette paper, CU; c is the cigarette circumference, cm; l3[ i]The original residual length of combustion of each port of the ith port, i is 2, …, n, cm; q is the flow rate through the filter tip, ml/s;

the formula (j) is:

pre[q]＝0.143×(1-E)²×l5×(17.5-q)/E³/s，

wherein, pre [ q ]]Is a function of the pressure drop, cmH, of the cigarette in the cross section at a flow rate of q₂O; e is the porosity of the tobacco shreds, and is set to be 0.76; l5 is the cross length, cm, of the tipping paper and the cigarette paper; s is the sectional area of cigarette in cm²；

The formula (k) is: pf [ q ] ═ kx (17.5-q) × l4,

wherein, pf [ q ]]Is a function of pressure drop of the filter rod at a flow rate of q, cmH₂O; k is the pressure drop coefficient of the filter rod, cmH2O s/cm³(ii) a q is the flow rate through the filter tip, ml/s; l4 is the length of the filter stick before the tipping paper is punched, cm. The filter stick pressure drop coefficient is obtained by calculation according to the filter stick pressure drop and the filter stick length, and each filter stick unit is provided with the corresponding filter stick pressure drop and the corresponding filter stick length.

Most preferably, the cigarette branch is calculated according to the formula (l) ignoring the pressure drop value of the cigarette paper air permeability, said formula (l) being: 2.5X (1-E) of pre0²×l3[i]/E³/s，

Wherein pre0 is the pressure drop value, cmH, of cigarette section neglecting the air permeability of cigarette paper₂O; e is the porosity of the tobacco shreds, and is set to be 0.76; l3[ i]The original residual length of combustion of each port of the ith port, i is 2, …, n, cm; s is the sectional area of cigarette in cm²。

Most preferably, the cigarette sectional area is calculated according to the formula (m), wherein the formula (m) is as follows:

s＝c×c/4/Pi，

wherein s is the sectional area of the cigarette in cm²(ii) a c is the cigarette circumference, cm; pi is the circumference ratio.

Most preferably, the length of the intersection of the tipping paper and the cigarette paper is calculated according to formula (n) which is:

l5＝lp-lf，

wherein l5 is the crossing length of the tipping paper and the cigarette paper, cm; lp is the width of the tipping paper, cm; lf is the length of the filter stick, cm.

Most preferably, the tipping paper perforated leading segment filter rod length is calculated according to formula (o):

l4＝lf-lt，

wherein l4 is the length of the filter stick before the tipping paper is punched, and is cm; lf is the length of the filter stick, cm; and lt is the distance between the outer hole belt of the tipping paper and the edge, cm.

Step S32, calculating the relative residual length of each combustion according to the formulas (10) and (11), wherein the formula (10) is as follows:

l7[1]＝l1-l6[1]，

wherein l7[1] is the relative residual length of combustion of each port 1, cm; l1 is the length of cigarette, cm; l6[1] is the relative burning length after the 1 st vent ventilation, cm;

the formula (11) is:

l7[i]＝l7[i-1]-58v-l6[i]，

wherein l7[ i ] is the relative residual length of combustion per port of the ith port, i is 2, …, n, cm; l7[ i-1] is the relative residual length of combustion per port of the i-1 st port, i is 2, …, n, cm; v is smoldering speed, cm/s; l6[ i ] is the relative combustion length after venting at the ith port, i being 1, …, n, cm.

Step S33, calculating the relative weight of each port of combustion according to the formula (12), wherein the formula (12) is as follows:

rm2[i]＝m1×l6[i]/l1，

wherein rm2[ i ] is the relative weight of each port, i.e. the combustion relative weight of the ith port, i is 1, …, n, g; m1 is the weight of tobacco shred, g; l6[ i ] is the relative combustion length after ventilation of the ith port, i is 1, …, n, cm; l1 is the length of the tobacco rod in cm.

More preferably, the relative weight per combustion is calculated by the formulas (p), (q) respectively as the relative weight per combustion of nicotine and the relative weight per combustion of tar,

the formula (p) is:

rm3[i]＝rm2[i]×nic，

wherein rm3[ i ] is the relative burning weight of nicotine per mouth, namely the relative burning weight of nicotine at the ith mouth, i is 1, …, n, g; rm2[ i ] is the relative weight of combustion at the ith port, i is 1, …, n, g; nic is the nicotine experience coefficient of the tobacco shred and is set as 1.11;

the formula (q) is:

rm4[i]＝rm2[i]×tar，

wherein rm4[ i ] is the relative weight of tar per combustion, i.e. the relative weight of tar at the ith combustion, i is 1, …, n, g; rm2[ i ] is the relative weight of combustion at the ith port, i is 1, …, n, g; tar is the empirical coefficient of tar in tobacco shreds, and is set to 1.22.

Preferably, in step S3, the pressure drop matching regression algorithm is that, based on the original combustion length of each cigarette cut obtained in S2, the pressure drop of the cigarette is matched with the pressure drop of the tipping paper, the forming paper and the filter rod system by a formula, the filter ventilation flow rate and the combustion cone inlet flow rate corresponding to the ith cut are obtained by each matching, and the combustion cone inlet flow rate determines the relative cut number, the relative residual length of each combustion cut and the relative weight of each combustion cut (including the weight of the cut tobacco relative to the combustion nicotine and tar).

And step S4, calculating to obtain the main component content value in the cigarette smoke by adopting a tobacco shred mouth-to-mouth release-filtration equation according to the original parameters of the cigarette in the step S2 and the relative parameters of the cigarette in the step S3.

Preferably, in step S4, the main components in the cigarette smoke include, but are not limited to, tar and nicotine.

Preferably, in step S4, the calculation of the main component content value (i.e. the release amount) in the cigarette smoke includes the following steps:

step S41, calculating the absolute release amount of nicotine in the whole cigarette according to a formula (13), wherein the formula (13) is as follows: f0tot ═ Sigma [ f0[ i ] ],

wherein f0tot is the absolute release amount of nicotine in the whole cigarette, mg; simga is the running sign Sigma, which is the sum of the logarithmic group f0[ i ]; f0[ i ] is the absolute release of nicotine from the combustion of the ith cigarette, i is 1, …, n, mg.

More preferably, the absolute release amount of nicotine during combustion of the ith cigarette is calculated according to the formula (r):

f0[i]＝g×qi[i]×(m3+f1[i-1,1]+f1[i-2,2]+…+f1[i,j])×e^(-a4×l2)，i＝1,…,n，

wherein f0[ i ] is the absolute release amount of nicotine burnt by the ith cigarette, mg; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; m3 is the original calculated weight of nicotine, g; a4 is an empirical coefficient, and the set value is 0.22; e is the base number of the natural logarithm; l2 is the original length of each port burning, cm; f1[ i, j ] is the absolute nicotine function, mg, retained by the cigarette in the combustion section of the ith port.

Further preferably, the absolute nicotine function retained by the cigarette in the ith combustion section in the jth combustion section is calculated according to the formula(s):

f1[i,j]＝g×qi[i]×(m3+f1[i-1,1]+f1[i-2,2]+…+f1[1,j])×e^{(-a4×l3[j+i-1])}×(1-e^(-a4×l2))，

i<j＝1,…,n；i+j-1＝1,…,n，

wherein f1[ i, j ] is the absolute nicotine function, mg, retained by the cigarette in the combustion section of the ith port in the jth combustion section; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; m3 is the original calculated weight of nicotine, g; a4 is an empirical coefficient, and the set value is 0.22; l3[ j + i-1] is the original remaining length of each port of j + i-1 in terms of cm; l2 is the original length of each port burning, cm; e is the base of the natural logarithm.

Step S42, calculating the absolute release amount of tar in the whole cigarette according to a formula (14), wherein the formula (14) is as follows: f02tot Sigma [ f02[ i ] ],

wherein, f02tot is the absolute release amount of tar in the whole cigarette, mg; simga is the running sign Sigma, which is the sum of the logarithmic group f02[ i ]; f02[ i ] is the absolute release of tar from the ith cigarette, i is 1, …, n, mg.

More preferably, the absolute release amount of tar during the combustion of the ith cigarette is calculated according to a formula (t), wherein the formula (t) is as follows:

f02[i]＝g×qi[i]×(m4+f2[i-1,1]+f2[i-2,2]+…+f2[i,j])×e^(-a4×l2)，i＝1,…,n，

wherein f02[ i ] is the absolute release amount of tar burnt by the ith cigarette, mg; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; m4 is the original calculated weight of tar, g; a4 is an empirical coefficient, and the set value is 0.22; e is the base number of the natural logarithm; l2 is the original length of each port burning, cm; f2[ i, j ] is the absolute tar function of the i th burning section of the cigarette, mg, of the j th burning section.

Further preferably, the absolute tar function retained in the j combustion section of the ith combustion cigarette is calculated according to a formula (u), wherein the formula (u) is as follows:

f2[i,j]＝g×qi[i]×(m4+f2[i-1,1]+f2[i-2,2]+…+f2[1,j])×e^{(-a4×l3[j+i-1])}×(1-e^(-a4×l2))，

i<j＝1,…,n；i+j-1＝1,…,n，

wherein f2[ i, j ] is the absolute tar function of the cigarette burning in the ith burning section, mg; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; m4 is the original calculated weight of tar, g; a4 is an empirical coefficient, and the set value is 0.22; l3[ j + i-1] is the original remaining length of each port of j + i-1 in terms of cm; l2 is the original length of each port burning, cm; e is the base of the natural logarithm.

Step S43, calculating the relative release amount of nicotine in the whole cigarette according to a formula (15), wherein the formula (15) is as follows: rf0tot ═ Sigma [ rf0[ i ] ] × (1-nicheff),

wherein, the rf0tot is the relative release amount of nicotine in the whole cigarette, mg; simga is a consecutive sign Sigma, which is the sum of logarithmic groups rf0[ i ]; rf0[ i ] is the relative release of nicotine from mouth i, i is 1, …, n, mg; the niceff is the nicotine filtration coefficient of the filter stick.

More preferably, the relative release of nicotine from the ith port is calculated according to formula (v):

rf0[i]＝g×qi[i]×(rm3[i]+rf1[i-1,1]+rf1[i-2,2]+…+rf1[i,j])×e^(-a4×l6[i])，i＝1,…,n，

wherein rf0[ i ] is the relative release amount of nicotine in mouth i, mg; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; rm3[ i ] is the relative weight of nicotine burning at mouth i, g; rf1[ i, j ] is a relative nicotine function, mg, of the cigarette burned in the ith burning section in the jth burning section; a4 is an empirical coefficient, and the set value is 0.22; e is the base number of the natural logarithm; l6[ i ] is the relative combustion length, cm, after venting at the ith port.

Further preferably, the relative nicotine function retained by the cigarette in the ith combustion section in the jth combustion section is calculated according to the formula (w):

rf1[i,j]＝g×qi[i]×(rm3[i]+rf1[i-1,1]+rf1[i-2,2]+…+rf1[1,j])×e^{(-a4×l7[j+i-1])}×(1-e^{(-a4×l6[i+j-1])})，

i<j＝1,…,n；i+j-1＝1,…,n，

wherein rf1[ i, j ] is a relative nicotine function, mg, intercepted by the cigarette in the combustion section of the ith combustion port; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; rm3[ i ] is the relative weight of nicotine in mouth i combustion, i is 1, …, n, g; a4 is an empirical coefficient, and the set value is 0.22; l7[ j + i-1] is the relative residual length of combustion of each port of j + i-1, cm; l6[ j + i-1] is the relative combustion length, cm, of the vent at the j + i-1; e is the base of the natural logarithm.

Step S44, calculating the relative release amount of tar in the whole cigarette according to a formula (16), wherein the formula (16) is as follows: rf02tot ═ Sigma [ rf02[ i ] ] × (1-tareff),

wherein, the rf02tot is the relative release amount of tar in the whole cigarette, mg; simga is a consecutive sign Sigma, which is the sum of logarithmic groups rf02[ i ]; rf02[ i ] is the relative release amount of tar in mouth i, mg; and tareff is the tar filtering coefficient of the filter stick.

More preferably, the relative amount of tar released from the ith port is calculated according to the formula (x):

rf02[i]＝g×qi[i]×(rm4[i]+rf2[i-1,1]+rf2[i-2,2]+…+rf2[i,j])×e^(-a4×l6[i])，i＝1,…,n，

wherein rf02[ i ] is the relative release amount of tar of the ith mouth in mg; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; rm4[ i ] is the relative weight of the i-th port combustion of tar, i is 1, …, n, g; a4 is an empirical coefficient, and the set value is 0.22; rf2[ i, j ] is a relative tar function, mg, of the cigarette burned in the ith burning section in the jth burning section; e is the base number of the natural logarithm; l6[ i ] is the relative combustion length, cm, after venting at the ith port.

Further preferably, the relative tar function of the cigarette burned at the ith burning section is calculated according to a formula (y):

rf2[i,j]＝g×qi[i]×(rm4[i]+rf2[i-1,1]+rf2[i-2,2]+…+rf2[1,j])×e^{(-a4×l7[j+i-1])}×(1-e^{(-a4×l6[i+j-1])})，

i<j＝1,…,n；i+j-1＝1,…,n，

wherein rf2[ i, j ] is a relative tar function, mg, intercepted by the cigarette burning in the ith burning section; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; rm4[ i ] is the relative weight of the i-th port combustion of tar, i is 1, …, n, g; a4 is an empirical coefficient, and the set value is 0.22; e is the base number of the natural logarithm; l7[ j + i-1] is the relative residual length of combustion of each port of j + i-1, cm; l6[ j + i-1] is the relative combustion length, cm, after ventilation at port j + i-1.

Preferably, in step S4, the tobacco shred mouth-by-mouth release-filter equation means that the original smoldering time, the original mouth number, the original length of combustion per mouth, the original weight of combustion per mouth, and the original remaining length of combustion per mouth are obtained in step S2, the relative mouth number, the relative remaining length of combustion per mouth, and the relative weight of combustion per mouth are obtained in step S3, and then the absolute release amount of tar and nicotine (original state) and the relative release amount of tar and nicotine, that is, f0tot, f02tot, rf0tot, and rf02tot are calculated mouth by assisting the tobacco shred mouth-by-mouth release-filter equation with the information of S2 and S3. The absolute release amount of tar and the relative release amount of tar jointly determine the content of tar in the cigarette smoke, and the absolute release amount of nicotine and the relative release amount of nicotine jointly determine the content of nicotine in the cigarette smoke.

The invention also provides a computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the analysis method described above.

The invention further provides a computer processing device comprising a processor and a computer readable storage medium, the processor executing a computer program on the computer readable storage medium to implement the steps of the analysis method described above.

As mentioned above, the method for analyzing the content of the main components in the cigarette smoke provided by the invention has the following beneficial effects:

(1) the method for analyzing the content of the main components in the cigarette smoke is based on the first physical law, assists the algorithm of the semi-empirical semi-theory, can be applied to a computer, calculates the influence of the change of a formula, auxiliary materials and a cigarette structure on the main stream smoke, and has certain universality.

(2) The method for analyzing the content of the main components in the cigarette smoke can accurately predict the content of the main components in the cigarette smoke, does not need a large number of repeated preliminary tests, and has the advantages of convenience and quickness.

(3) According to the analysis method for the content of the main components in the cigarette smoke, the actual cigarette design target can be met only by reasonably adjusting the correlation among the parameters and properly and reasonably adjusting the leaf group formula, the auxiliary materials and the rolling structure parameters according to the prediction result, and the cost reduction and efficiency improvement effects can be well achieved without modifying the existing cigarette production and processing equipment.

Drawings

FIG. 1 is a schematic flow chart of the method for analyzing the content of the main components in the cigarette smoke according to the present invention.

Detailed Description

The present invention is further illustrated below with reference to specific examples, which are intended to be illustrative only and not to limit the scope of the invention.

The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention.

It should be understood that the processing equipment or devices not specifically mentioned in the following examples are conventional in the art; all pressure values and ranges refer to relative pressures.

Furthermore, it is to be understood that one or more method steps mentioned in the present invention does not exclude that other method steps may also be present before or after the combined steps or that other method steps may also be inserted between these explicitly mentioned steps, unless otherwise indicated; it is also to be understood that a combined connection between one or more devices/apparatus as referred to in the present application does not exclude that further devices/apparatus may be present before or after the combined device/apparatus or that further devices/apparatus may be interposed between two devices/apparatus explicitly referred to, unless otherwise indicated. Moreover, unless otherwise indicated, the numbering of the various method steps is merely a convenient tool for identifying the various method steps, and is not intended to limit the order in which the method steps are arranged or the scope of the invention in which the invention may be practiced, and changes or modifications in the relative relationship may be made without substantially changing the technical content.

As shown in figure 1, the invention provides a method for analyzing the content of main components in cigarette smoke, which comprises the following steps:

step S1, determining basic parameters of the cigarette;

step S2, calculating the original parameters of the cigarette according to the basic parameters of the cigarette in the step S1 by adopting a smoldering time recurrence method;

In step S1, the basic parameters of the cigarette include, but are not limited to, tobacco type, tobacco weight, tobacco length, cigarette circumference, filter length, filter pressure drop, tipping paper air permeability, total tipping paper pore band width, tipping paper outer pore band spacing, forming paper air permeability, cigarette paper combustion improver content, cigarette paper gram weight, nicotine filtration coefficient, and tar filtration coefficient.

In step S2, the original parameters of the cigarette include, but are not limited to, original smoldering time, original number of mouths, original length of burn per mouth, original weight of burn per mouth, and original remaining length of burn per mouth.

In step S3, the relative parameters of the cigarette include, but are not limited to, relative number of ports, relative remaining length of combustion per port, and relative weight of combustion per port.

In step S4, the main components in the cigarette smoke include, but are not limited to, tar and nicotine.

Example 1

In a computer, according to the figure 1, the main components in the cigarette smoke are analyzed, and the steps are as follows:

1. selecting cigarette experimental samples, determining basic parameters of the cigarettes, wherein specific data of the basic parameters of the cigarettes are shown in table 1.

TABLE 1 basic parameters of cigarette experimental samples

Cigarette basic parameters	Numerical value
		Kind of cut tobacco	100% of shredded tobacco leaves
Weight of tobacco shred	0.7g
		Length of cigarette segment	5.1cm
Cigarette circumference	2.45cm
		Length of filter stick	3cm
Pressure drop of filter stick	8cmH2
		Coefficient of nicotine filtration	0.4
Coefficient of tar filtration	0.54
		Air permeability of tipping paper	Electrostatic drilling 150CU
Air permeability of formed paper	10000CU
		Outer hole strip distance edge of tipping paper	1.45cm
Width of tipping paper	3cm
		Air permeability of cigarette paper	60CU
Combustion improver content of cigarette paper	1.8% pure potassium salt
		Gram weight of cigarette paper	27g

2. And (4) calculating the original parameters of the cigarettes according to the basic parameters of the cigarettes in the step S1 by adopting a smoldering time recurrence method, wherein the specific data of the original parameters of the cigarettes are shown in the table 2.

TABLE 2 original parameters of cigarette experimental samples

Raw parameters of cigarette	Numerical value
		Original smoldering time t1	570s
Number of original ports n	7.4 mouth
		Original weight of each combustion m2(g)	0.023、0.023、0.023、0.023、0.023、0.023、0.023、0.009
Original residual length l3[ i ] of each port combustion](cm)	5.2、4.6、3.9、3.3、2.6、2.0、1.3、0.8

The specific calculation process of the original parameters of the cigarette is as follows:

step 1, after determining that the type of cut tobacco is 100% cut tobacco, calculating the original smoldering time t1 according to a formula (1), wherein the formula (1) is as follows: t1 ═ 73.3 × l1+279 × c-33.1 × a-19.5 × mc +6.86 × ps +11.7 × pe +6.49 × pp +10.4 × pl, giving t1 of 570 s.

Step 2, calculating an original opening number n according to a formula (2) through the original smoldering time t1, wherein the formula (2) is as follows: n is 0.0121 × t1+0.487, and 7.4 mouths as n is obtained.

Step 3, calculating the original combustion length l2 of each port according to formulas (3), (4) and (5) through the original smoldering time t1 and the original port number n, wherein the formula (3) is as follows: l1/t1, the formula (4) is: t2 is 58 × n, and the formula (5) is: l2 ═ v × (t1-t2)/n, giving l2 of 0.18 cm.

Step 4, calculating the original weight m2 of each port of combustion according to a formula (6) through the original length l2 of each port of combustion, wherein the formula (6) is as follows: m2 ═ m1 × l2/l1, giving m2 of 0.023. Wherein, the original weight m2 of each combustion is respectively calculated according to the formulas (a) and (b) to obtain the original calculated weight m3 of nicotine and the original calculated weight m4 of tar, and the formula (a) is as follows: m3 ═ m2 xnic, and the formula (b) is: m4 ═ m2 × tar, yielding m3 of 0.026 and m4 of 0.028.

Step 5, calculating the original residual length l3[ i ] of each port of combustion according to formulas (7) and (8) through the original length l2 of each port of combustion, wherein the formula (7) is as follows: l3[1] ═ l1-l2, and the formula (8) is: l3[ i ] ═ l3[ i-1] -58v-l2, giving l3[ i ] of 5.2, 4.6, 3.9, 3.3, 2.6, 2.0, 1.3, 0.8.

3. Adopting a pressure drop matching regression algorithm according to the basic parameters of the cigarettes in the step S1 and the original parameters of the cigarettes in the step S2

The relative parameters of the cigarettes are calculated, and the specific data of the relative parameters of the cigarettes are shown in the table 3.

TABLE 3 relative parameters of cigarette experimental samples

Relative parameters of cigarettes	Numerical value
		Relative number of openings rn	7.7 mouth
Relative residual length of combustion per port l7[ i ]](mm)	5.2、4.6、4.0、3.3、2.7、2.0、1.4、0.8
		Relative weight of combustion per port rm2[ i ]](g)	0.022、0.022、0.022、0.022、0.022、0.022、0.022、0.016
Mouth i combustion relative weight rm3[ i ] of nicotine](g)	0.024，0.024，0.024，0.024，0.024，0.024，0.024，0.017
		The ith port burn of Tar relative weight rm4[ i](g)	0.027，0.027，0.027，0.027，0.027，0.027，0.027，0.019

The specific calculation process of the relative parameters of the cigarettes is as follows:

step 1, calculating a pressure drop-flow function pf [ q ] of the filter stick according to a formula (k), wherein the formula (k) is as follows: pf [ q ] ═ k × (17.5-q) × l4, resulting in pf [ q ] ═ 0.214 (17.5-q).

Step 2, calculating a pressure drop-flow function pre [ q ] of the cigarette under the ventilation condition without the cigarette paper according to a formula (j)]The formula (j) is: pre [ q ]]＝0.143×(1-E)²×l5×(17.5-q)/E³S, obtaining pre [ q ]]0.016 (17.5-q). The cigarette sectional area s is calculated according to a formula (m), wherein the formula (m) is as follows: s ═ cxc/4/Pi, and s was obtained as 0.478.

Step 3, calculating a pressure drop-flow function pr [ q ] under the ventilation condition of the cigarette paper according to a formula (i), wherein the formula (i) is as follows:

obtaining pr [ q ]]＝0.07Tanh[52.7×Sqrt[7×10^-5+2.9×10^-6(17.5-q)]]/Sqrt[7×10^-5+2.9×10^-6(17.5-q)]. Wherein, the pressure drop value pre0 of the cigarette section neglecting the air permeability of the cigarette paper is calculated according to the formula (l): 2.5X (1-E) of pre0²×l3[i]/E³And/s, a pre0 value of (3.64, 3.17, 2.74, 2.28, 1.83, 1.38, 0.93, 0.56) is obtained.

Step 4, calculating a tipping paper-forming paper pressure drop-flow function value ptpw [ q ] according to a formula (h)]The formula (h) is: ptpw [ q ]]＝(q/c)²×(2.56×10⁸/a2²+3.6×10⁴/a3²) Obtaining ptpw [ q ]]Is 0.36q²。

Step 5, calculating a q value according to a regression pressure drop matching equation in a formula (g), wherein the formula (g) is as follows: q ═ FindRoot [ ptpw [ q ] ═ pr [ q ] + pre [ q ] + pf [ q ], q- >0],

a q of (2.15,2.09,2.02,1.94,1.85,1.76,1.67,1.58) was obtained.

Step 6, calculating the flow velocity value of the cigarette combustion inlet according to a formula (e), wherein the formula (e) is as follows: qi [ i ] ═ 0.9 × (17.5-q [ i ]), obtaining qi [ i ] of (13.8, 13.87, 13.94, 14, 14.1, 14.2, 14.25, 14.3).

And 7, calculating the relative combustion length of the ventilation of the ith port according to a formula (d) by using the flow rate value of the combustion inlet of the cigarette, wherein the formula (d) is that l6[ i ] = qi [ i ]/17.5 × l2, and l6[ i ] is obtained to be (0.168, 0.169, 0.17, 0.171, 0.171, 0.173, 0.174, 0.174).

Step 8, calculating the average relative burning length according to a formula (c) through the relative burning length after the ventilation of the ith opening, wherein the formula (c) is as follows: l6avg ═ Simga (l6[ i ])/n, giving l6avg of 0.17.

Step 9, calculating a relative port number rn according to a formula (9), wherein the formula (9) is as follows: rn ═ l1/(l6avg +58v), giving an rn of 7.7.

Step 10, calculating the relative residual length l7[ i ] of each port according to the formulas (10) and (11), wherein the formula (10) is as follows: l7[1] ═ l1-l6[1], the formula (11) being: l7[ i ] ═ l7[ i-1] -58v-l6[ i ], yielding l7[ i ] of 5.2, 4.6, 4.0, 3.3, 2.7, 2.0, 1.4, 0.8.

Step S11, calculating the combustion relative weight rm2[ i ] of the ith port according to the formula (12), wherein the formula (12) is as follows: rm2[ i ] ═ m1 × l6[ i ]/l1, yielding rm2[ i ] of (0.022, 0.016). Then, the ith combustion relative weight of nicotine and the ith combustion relative weight of tar are respectively calculated according to formulas (p) and (q), wherein the formula (p) is as follows: rm3[ i ] ═ rm2[ i ] × nic, said formula (q) is: rm4[ i ] ═ rm2[ i ] × tar, rm3[ i ] was obtained as (0.024, 0.024, 0.024, 0.024, 0.024, 0.024, 0.024, 0.017), rm4[ i ] was obtained as (0.027, 0.027, 0.027, 0.027, 0.027, 0.027, 0.019).

4. And (4) calculating to obtain the content value of the main components in the cigarette smoke by adopting a tobacco shred mouth-to-mouth release-filtration equation according to the original parameters of the cigarette in the step S2 and the relative parameters of the cigarette in the step S3, wherein the specific data are shown in a table 4.

TABLE 4 content of main constituents in smoke of cigarette experimental samples

Main component of cigarette smoke	Numerical value
		Absolute release of primary nicotine f0tot	1.7mg
Initial nicotine test values	1.6mg
		Relative deviation of original nicotine	6.25％
Relative release amount of nicotine rf0tot	0.92mg
		Experimental value of nicotine	0.97mg
Relative deviation of nicotine	5.15％
		Absolute amount of release f02tot of original tar	23.7mg
Original tar experimental value	23.8mg
		Relative deviation of original tar	0.4％
Relative release amount of tar rf02tot	13mg
		Experimental value of tar	12mg
Relative deviation of tar	8.33％

The specific calculation process of the main component content value in the cigarette smoke is as follows:

step 1, calculating an absolute nicotine function retained by a cigarette burning at the ith burning section according to a formula(s), wherein the formula(s) is as follows: f1[ i, j ]]＝g×qi[i]×(m3+f1[i-1,1]+f1[i-2,2]+…+f1[1,j])×e^{(-a4×l3[j+i-1])}×(1-e^(-a4×l2))，i<j is 1, …, n; i + j-1 ═ 1, …, n, yielding f1[1,1]＝0.0041,f1[1,2]＝0.0038,…f[6,1]＝0.0018。

Step 2, calculating the absolute release amount f0[ i ] of the nicotine burnt by the ith cigarette according to a formula (r), wherein the formula (r) is as follows:

f0[ i ] (0.16, 0.18, 0.2, 0.22, 0.25, 0.27, 0.29, 0.09) was obtained.

Step 3, calculating the absolute release amount of nicotine (in the original state) in the whole cigarette according to a formula (13), wherein the formula (13) is as follows: f0tot ═ Sigma [ f0[ i ] ], giving an f0tot of 1.7 mg.

Step 4, calculating an absolute tar function of the cigarette burned in the ith port at the jth burning section according to a formula (u), wherein the formula (u) is as follows: f2[ i, j ]]＝g×qi[i]×(m4+f2[i-1,1]+f2[i-2,2]+…+f2[1,j])×e^{(-a4×l3[j+i-1])}×(1-e^(-a4×l2))，i<j is 1, …, n; i + j-1 ═ 1, …, n, yielding f2[1,1]＝0.285,f2[1,2]＝0.027,…,f2[6,1]＝0.035。

Step 5, calculating the absolute release amount of the burnt tar of the ith cigarette according to a formula (t), wherein the formula (t) is as follows:

f02[ i ] (1.24, 1.74, 2.4, 3.3, 4.3, 5.7, 4.3, 0.73) was obtained.

Step 6, calculating the absolute release amount of tar (in the original state) in the whole cigarette according to a formula (14), wherein the formula (14) is as follows: f02tot ═ Sigma [ f02[ i ] ], giving f02tot of 23.7 mg.

Step 7, calculating a relative nicotine function of the cigarette burned in the ith burning section in the jth burning section according to a formula (w), wherein the formula (w) is as follows:

i < j ═ 1, …, n; i + j-1 is 1, …, n, yielding rf1[1,1] 0.004, rf1[1,2] 0.0038, …, rf1[6,1] 0.0018.

Step 8, calculating the relative release amount of nicotine of the ith mouth according to a formula (v), wherein the formula (v) is as follows: rf0[ i]＝g×qi[i]×(rm3[i]+rf1[i-1,1]+rf1[i-2,2]+…+rf1[i,j])×e^(-a4×l6[i])I 1, …, n, to obtain rf0[ i]＝(0.08，0.095，0.11，0.12，0.13，0.15，0.14，0.09)。

Step 9, calculating the relative release amount rf0tot of nicotine in the whole cigarette according to a formula (15), wherein the formula (15) is as follows: rf0tot ═ Sigma [ rf0[ i ] ] × (1-nicheff), and rf0tot of 0.92mg was obtained.

Step S10, calculating the relative tar function of the cigarette burned in the ith burning section in the jth burning section according to a formula (y), wherein the formula (y) is as follows:

i < j ═ 1, …, n; i + j-1 ═ 1, …, n, yielding rf2[1,1] ═ 0.028, rf2[1,2] ═ 0.027, …, rf2[6,1] ═ 0.035.

Step S11, calculating the relative release amount of tar in the ith port according to a formula (x), wherein the formula (x) is as follows:

obtain rf02[ i ] (0.63, 0.9, 1.3, 1.7, 2.3, 3.1, 2.3, 0.73).

Step 12, calculating the relative release amount of tar in the whole cigarette according to a formula (16), wherein the formula (16) is as follows: rf02tot ═ Sigma [ rf02[ i ] ] × (1-tareff), and rf02tot of 13mg was obtained.

In the steps 2 and 4, detecting tar in the cigarette smoke according to the regulation of the national standard GB/T19609-once 2004 to obtain an original tar experimental value and a tar experimental value; detecting nicotine in cigarette smoke according to the regulation of the national standard GB/T23355-2009 to obtain an original nicotine experimental value and a nicotine experimental value. The specific detection equipment adopts the regulations in the national standard GB/T16450-2004 to set the parameters of the experimental equipment: the suction volume was 35mL, the suction frequency was 60 s/port, and the suction time was 2 s. The relative deviation of the content of the main components in the cigarette smoke is calculated according to a formula (z), wherein the formula (z) is w ═ Abs ((w1-w2)/w2), w is the relative deviation, w1 is a predicted value of a certain component in the smoke, and w2 is an experimental value of a corresponding component in the smoke.

As can be seen from Table 4, when the method provided by the invention is used for analyzing the contents of tar and nicotine in cigarette smoke, the experimental value and the predicted value are very close, the relative deviation of the contents of tar and nicotine is 8.33% and 5.15%, respectively, the relative deviation of the predicted value of the content of the main components in the obtained cigarette smoke and the experimental value of the corresponding components is less than 10%, and the design requirement of the cigarette is completely met.

In conclusion, the method for analyzing the content of the main components in the cigarette smoke can be applied to a computer, can calculate the influence of the change of the formula, the auxiliary materials and the cigarette structure on the main stream smoke, can accurately predict the content of the main components in the cigarette smoke, does not need to modify the existing cigarette production and processing equipment, does not need to perform a large number of repeated early-stage pre-tests, and has the advantages of convenience and quickness. Therefore, the invention effectively overcomes various defects in the prior art and has high industrial utilization value.

The foregoing embodiments are merely illustrative of the principles and utilities of the present invention and are not intended to limit the invention. Any person skilled in the art can modify or change the above-mentioned embodiments without departing from the spirit and scope of the present invention. Accordingly, it is intended that all equivalent modifications or changes which can be made by those skilled in the art without departing from the spirit and technical spirit of the present invention be covered by the claims of the present invention.

Claims

1. A method for analyzing the content of main components in cigarette smoke comprises the following steps:

step S1, determining basic parameters of the cigarette;

step S4, calculating to obtain the main component content value in the cigarette smoke by adopting a tobacco shred mouth-by-mouth release-filtration equation according to the original parameters of the cigarette in the step S2 and the relative parameters of the cigarette in the step S3;

the basic parameters of the cigarette include but are not limited to the type of cut tobacco, the weight of the cut tobacco, the length of a cigarette segment, the circumference of the cigarette, the length of a filter stick, the pressure drop of the filter stick, the air permeability of tipping paper, the total width of a pore zone of the tipping paper, the width of the tipping paper, the distance between outer pore zones of the tipping paper, the air permeability of forming paper, the air permeability of cigarette paper, the content of a combustion improver of the cigarette paper, the gram weight of the cigarette paper, the nicotine filtration coefficient;

in step S2, the original parameters of the cigarette include, but are not limited to, original smoldering time, original number of cigarettes, original length of burning of each cigarette, original weight of burning of each cigarette, and original remaining length of burning of each cigarette; the method for calculating the original parameters of the cigarettes comprises the following steps:

t1＝73.3×l1+279×c-33.1×a-19.5×mc+6.86×ps+11.7×pe+6.49×pp+10.4×pl，

wherein t1 is the original smoldering time, s; l1 is the length of cigarette, cm; c is the cigarette circumference, cm; a is the content of the combustion improver of the cigarette paper in percent; mc is the gram weight of cigarette paper, g/m²(ii) a ps is the ratio of cut stems,%; pe is the proportion of the expanded filaments,%; pp is flake ratio,%; pl is leaf-silk ratio,%;

n is 0.0121 × t1+0.487, wherein n is the number of original ports; t1 is the original smoldering time, s;

step S23, calculating the original length of each port combustion according to the formulas (3), (4) and (5) through the original smoldering time and the original port number, wherein the formula (3) is as follows: l1/t1, wherein v is smoldering speed, cm/s; l1 is the length of cigarette, cm; t1 is the original smoldering time, s;

the formula (4) is: t2 ═ 58 × n, where t2 is the actual smoldering time, s; n is the original number of ports;

the formula (5) is: l2 ═ v × (t1-t2)/n, where l2 is the original length of each port burned, cm; v is smoldering speed, cm/s; t1 is the original smoldering time, s; t2 is the actual smoldering time, s; n is the original number of ports;

step S24, calculating the original weight of each port of combustion according to a formula (6) through the original length of each port of combustion, wherein the formula (6) is as follows: m2 ═ m1 × l2/l1, where m2 is the original weight per port burned, g; m1 is the weight of tobacco shred, g; l2 is the original length of each port burning, cm; l1 is the length of cigarette, cm;

step S25, calculating the original residual length of each port of combustion according to formulas (7) and (8) through the original length of each port of combustion, wherein the formula (7) is as follows: l3[1] ═ l1-l2, where l3[1] is the original residual length of combustion per port 1, cm; l1 is the length of cigarette, cm; l2 is the original length of each port burning, cm;

the formula (8) is: l3[ i ] ═ l3[ i-1] -58v-l2, where l3[ i ] is the original residual length of combustion per port of the ith port, i ═ 2, …, n, cm; l3[ i-1] is the original residual length of combustion of each port i-1, i is 2, …, n, cm; v is smoldering speed, cm/s; l2 is the original length of each port burning, cm;

in step S3, the relative parameters of the cigarette include, but are not limited to, relative number of cigarettes, relative remaining length of each cigarette, and relative weight of each cigarette;

in the step S2, the smoldering time recurrence method is to calculate the original number of burning mouths of the cigarette through the original smoldering time of the cigarette, further calculate the original burning length and the original burning weight of each mouth of the cigarette during burning, further obtain the original calculated weight of nicotine and tar of each mouth of the cigarette, and calculate to obtain the original burning remaining length of each mouth of the cigarette;

in the step S3, the pressure drop matching regression algorithm is that the cigarette pressure drop is matched with the pressure drop of tipping paper, forming paper and a filter stick system through a formula on the basis of the original burning length of each port of the cigarette obtained through S2, the filter tip ventilation flow rate and the burning cone inlet flow rate corresponding to the ith port are obtained through each matching, and the relative port number, the relative residual length of each burning port and the relative burning weight of each burning port are determined through the burning cone inlet flow rate;

the tobacco shred mouth-by-mouth release-filtration equation is that the original smoldering time, the original mouth number, the original burning length of each mouth, the original burning weight of each mouth and the original remaining length of each mouth are obtained through the step S2, the relative mouth number, the relative remaining length of each mouth and the relative burning weight of each mouth are obtained through the step S3, and the absolute release amount of tar and nicotine in the original state and the relative release amount of tar and nicotine are obtained through mouth-by-mouth calculation through the information of S2 and S3 to assist the tobacco shred mouth-by-mouth release-filtration equation; the absolute release amount of tar and the relative release amount of tar jointly determine the content of tar in the cigarette smoke, and the absolute release amount of nicotine and the relative release amount of nicotine jointly determine the content of nicotine in the cigarette smoke.

2. The method according to claim 1, wherein in step S24, the initial weight of each combustion is calculated according to the formula (a) and (b) respectively as the initial calculated weight of nicotine and the initial calculated weight of tar,

the formula (a) is: m3 ═ m2 xnic, where m3 is the original calculated weight of nicotine, g; m2 is the original weight per port, g; nic is the nicotine experience coefficient of the tobacco shreds, which is set as 1.11, and the nicotine experience coefficient nic of the tobacco shreds is a coefficient for adjusting the nicotine release amount of different leaf groups;

the formula (b) is: m4 ═ m2 × tar, where m4 is the original calculated weight of tar, g; m2 is the original weight per port, g; tar is the empirical coefficient of tar of the tobacco shreds, and is set to be 1.22; the tar empirical coefficient tar is a coefficient for adjusting the tar release amount of different leaf groups.

3. The method for analyzing the content of the main components in the cigarette smoke according to the claim 1, wherein in the step S3, the calculation of the relative parameters of the cigarette comprises the following steps:

step S31, calculating the relative port number according to the formula (9), wherein the formula (9) is as follows: rn l1/(l6avg +58v),

wherein rn is the relative number of ports; l1 is the length of cigarette, cm; l6avg is the average relative burn length, cm; v is smoldering speed, cm/s;

step S32, calculating the relative residual length of each combustion according to the formulas (10) and (11), wherein the formula (10) is as follows: l7[1] ═ l1-l6[1],

the formula (11) is: l7[ i ] ═ l7[ i-1] -58v-l6[ i ],

wherein l7[ i ] is the relative residual length of combustion per port of the ith port, i is 2, …, n, cm; l7[ i-1] is the relative residual length of combustion per port of the i-1 st port, i is 2, …, n, cm; v is smoldering speed, cm/s; l6[ i ] is the relative combustion length after ventilation of the ith port, i is 1, …, n, cm;

step S33, calculating the relative weight of each port of combustion according to the formula (12), wherein the formula (12) is as follows: rm2[ i ] ═ m1 × l6[ i ]/l1,

4. The method for analyzing the content of the main component in the cigarette smoke according to claim 3, wherein in the step S31,

the average relative burn length is calculated according to equation (c) which is: l6avg ═ Simga (l6[ i ])/n,

wherein l6avg is the average relative burning length, cm; simga is the running sign Sigma, which is the sum of the logarithmic group l6[ i ]; l6[ i ] is the relative combustion length after ventilation of the ith port, i is 1, …, n, cm; n is the original number of ports;

the relative combustion length after the ventilation of the ith port is calculated according to a formula (d), wherein the formula (d) is as follows: l6[ i ] ═ qi [ i ]/17.5 × l2,

wherein l6[ i ] is the relative combustion length after the ith port is ventilated, i is 1, …, n, cm; qi [ i ] is the cigarette combustion inlet flow rate qi value at the ith port, i is 1, …, n, ml/s; l2 is the original length of each port burning, cm;

the value of the cigarette combustion inlet flow rate qi at the ith port is calculated according to a formula (e), wherein the formula (e) is as follows: qi [ i ] ═ 0.9 x (17.5-q [ i ]),

wherein qi is the cigarette burning inlet flow rate qi value of the ith port, i is 1, …, n, ml/s; q [ i ] is the ventilation flow rate q value of the filter at the length l3[ i ] of the ith port, i is 1, …, n, ml/s;

the flow rate of the combustion inlet of the cigarette at the ith port conforms to a flow rate constant equation in a formula (f), wherein the formula (f) is as follows: qi [ i ] + qc [ i ] + q [ i ] ═ 17.5,

wherein qi is the cigarette burning inlet flow rate qi value of the ith port, i is 1, …, n, ml/s; qc [ i ] is the flow rate qc of the wrapping paper at the ith port, i is 1, …, n, ml/s; q [ i ] is the ventilation flow rate q value of the filter at the length l3[ i ] of the ith port, i is 1, …, n, ml/s;

the flow rate of the filter tip is calculated according to a regression pressure drop matching equation in a formula (g), wherein the formula (g) is as follows: q ═ FindRoot [ ptpw [ q ] ═ pr [ q ] + pre [ q ] + pf [ q ], q- >0],

5. The method for analyzing the content of the main components in the cigarette smoke according to claim 4, wherein the pressure drop function of the combination of the tipping paper and the forming paper at the flow velocity q is a pressure drop-flow function of the tipping paper-forming paper calculated according to a formula (h), the pressure drop-flow function of the cigarette containing paper under the ventilation condition is calculated according to a formula (i) by considering the pressure drop function of the air permeability of the cigarette paper at the flow velocity q, the pressure drop-flow function of the cigarette without the ventilation condition of the cigarette paper is calculated according to a formula (j) by considering the pressure drop function of the air permeability of the cigarette paper at the flow velocity q, the pressure drop-flow function of the filter stick at the flow velocity q is calculated according to a formula (k),

the formula (i) is:

wherein, pr [ q ]]Considering for a segment the pressure drop function, cmH, of the air permeability of the cigarette paper at a flow rate q₂O; pre0 is the pressure drop value of cigarette section neglecting the air permeability of cigarette paper，cmH₂O; tanh is a hyperbolic tangent function; sqrt is a root function; a1 is the air permeability of cigarette paper, CU; c is the cigarette circumference, cm; l3[ i]The original residual length of combustion of each port of the ith port, i is 2, …, n, cm; q is the flow rate through the filter tip, ml/s;

the formula (j) is: pre [ q ]]＝0.143×(1-E)²×l5×(17.5-q)/E³/s，

The formula (k) is: pf [ q ] ═ kx (17.5-q) × l4,

wherein, pf [ q ]]Is a function of pressure drop of the filter rod at a flow rate of q, cmH₂O; k is the pressure drop coefficient of the filter rod, cmH2O s/cm³(ii) a q is the flow rate through the filter tip, ml/s; l4 is the length of the filter stick before the tipping paper is punched, and is cm; the filter stick pressure drop coefficient is obtained by calculation according to the filter stick pressure drop and the filter stick length, and each filter stick unit is provided with the corresponding filter stick pressure drop and the corresponding filter stick length.

6. The method for analyzing the content of main components in cigarette smoke according to claim 5,

the pressure drop value of the cigarette section neglecting the air permeability of the cigarette paper is calculated according to a formula (l), wherein the formula (l) is as follows: 2.5X (1-E) of pre0²×l3[i]/E³/s，

Wherein pre0 is the pressure drop value, cmH, of cigarette section neglecting the air permeability of cigarette paper₂O; e is the porosity of the tobacco shreds, and is set to be 0.76; l3[ i]The original residual length of combustion of each port of the ith port, i is 2, …, n, cm; s is the sectional area of cigarette in cm²；

The cigarette sectional area is calculated according to a formula (m), wherein the formula (m) is as follows: s-c × c/4/Pi,

wherein s is the sectional area of the cigarette in cm²(ii) a c is the cigarette circumference, cm; pi is the circumference ratio;

the crossing length of the tipping paper and the cigarette paper is calculated according to a formula (n), wherein the formula (n) is as follows: l5 ═ lp-lf,

wherein l5 is the crossing length of the tipping paper and the cigarette paper, cm; lp is the width of the tipping paper, cm; lf is the length of the filter stick, cm;

the length of the filter stick at the front section of the tipping paper perforation is calculated according to a formula (o), wherein the formula (o) is as follows: l4 ═ lf-lt,

7. The method for analyzing the content of main components in cigarette smoke according to claim 3, wherein in step S33, the relative weight of burning of each mouth is calculated according to the formulas (p), (q) to calculate the relative weight of burning of each mouth of nicotine and the relative weight of burning of each mouth of tar,

the formula (p) is: rm3[ i ] ═ rm2[ i ] × nic,

the formula (q) is: rm4[ i ] ═ rm2[ i ] × tar,

8. The method for analyzing the content of the main component in the cigarette smoke according to claim 1, wherein in step S4, the main component in the cigarette smoke includes but is not limited to tar, nicotine; the calculation of the main component content value, namely the release amount in the cigarette smoke comprises the following steps:

wherein f0tot is the absolute release amount of nicotine in the whole cigarette, mg; simga is the running sign Sigma, which is the sum of the logarithmic group f0[ i ]; f0[ i ] is the absolute release amount of nicotine burnt by the ith cigarette, i is 1, …, n, mg;

wherein, f02tot is the absolute release amount of tar in the whole cigarette, mg; simga is the running sign Sigma, which is the sum of the logarithmic group f02[ i ]; f02[ i ] is the absolute release amount of tar burned by the ith cigarette, i is 1, …, n, mg;

wherein, the rf0tot is the relative release amount of nicotine in the whole cigarette, mg; simga is a consecutive sign Sigma, which is the sum of logarithmic groups rf0[ i ]; rf0[ i ] is the relative release of nicotine from mouth i, i is 1, …, n, mg; the niceff is the nicotine filtration coefficient of the filter stick;

9. The method for analyzing the content of the main components in the cigarette smoke according to claim 8, wherein in step S41, the absolute release amount of nicotine during combustion of the ith cigarette is calculated according to a formula (r), and the formula (r) is as follows:

wherein f0[ i ] is the absolute release amount of nicotine burnt by the ith cigarette, mg; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; m3 is the original calculated weight of nicotine, g; a4 is an empirical coefficient, and the set value is 0.22; e is the base number of the natural logarithm; l2 is the original length of each port burning, cm; f1[ i, j ] is the absolute nicotine function, mg, retained by the cigarette in the combustion section of the ith port in the jth combustion section;

the absolute nicotine function of the cigarette burning at the ith burning section is calculated according to a formula(s):

i<j＝1,…,n；i+j-1＝1,…,n，

10. The method for analyzing the content of the main components in the cigarette smoke according to claim 8, wherein in the step S42, the absolute release amount of the tar during the combustion of the ith cigarette is calculated according to a formula (t), and the formula (t) is as follows:

wherein f02[ i ] is the absolute release amount of tar burnt by the ith cigarette, mg; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; m4 is the original calculated weight of tar, g; a4 is an empirical coefficient, and the set value is 0.22; e is the base number of the natural logarithm; l2 is the original length of each port burning, cm; f2[ i, j ] is the absolute tar function of the i th burning cigarette retained in the j th burning section, mg;

the absolute tar function of the cigarette burned in the ith port intercepted in the jth combustion section is calculated according to a formula (u), wherein the formula (u) is as follows:

i<j＝1,…,n；i+j-1＝1,…,n，

11. The method for analyzing the content of the main components in the cigarette smoke according to claim 8, wherein in step S43, the relative release amount of nicotine in the ith port is calculated according to the formula (v):

wherein rf0[ i ] is the relative release amount of nicotine in mouth i, mg; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; rm3[ i ] is the relative weight of nicotine burning at mouth i, g; rf1[ i, j ] is a relative nicotine function, mg, of the cigarette burned in the ith burning section in the jth burning section; a4 is an empirical coefficient, and the set value is 0.22; e is the base number of the natural logarithm; l6[ i ] is the relative combustion length, cm, after ventilation of the ith port;

the relative nicotine function of the cigarette burning at the ith burning section is calculated according to a formula (w), wherein the formula (w) is as follows:

i<j＝1,…,n；i+j-1＝1,…,n，

12. The method for analyzing the content of the main components in the cigarette smoke according to claim 8, wherein in the step S44, the relative release amount of the tar in the ith mouth is calculated according to a formula (x), and the formula (x) is as follows:

wherein rf02[ i ] is the relative release amount of tar of the ith mouth in mg; g is an empirical coefficient, and the set value is 0.297; qi [ i ] is the value of the flow rate qi of the cigarette burning inlet at the ith port, ml/s; rm4[ i ] is the relative weight of the i-th port combustion of tar, i is 1, …, n, g; a4 is an empirical coefficient, and the set value is 0.22; rf2[ i, j ] is a relative tar function, mg, of the cigarette burned in the ith burning section in the jth burning section; e is the base number of the natural logarithm; l6[ i ] is the relative combustion length, cm, after ventilation of the ith port;

the relative tar function of the cigarette burned in the ith burning section is calculated according to a formula (y), wherein the formula (y) is as follows:

i<j＝1,…,n；i+j-1＝1,…,n，

13. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the analysis method according to any one of claims 1 to 12.

14. A computer processing apparatus comprising a processor and a computer readable storage medium according to claim 13, wherein the processor executes a computer program on the computer readable storage medium to perform the steps of the analysis method according to any one of claims 1 to 12.