DOI: 10.11817/j.issn.1672-7207.2019.04.017

热定型过程弹力针织物热塑性应变模型研究

周丽春，刘顺菁，金福江，楚贝贝

(华侨大学 机电及自动化学院，福建 厦门，361021)

摘要：为了得到针织物在热定型过程中的应变理论模型，首先通过分析热定型过程中弹力针织物在常温拉伸、加热伸缩、冷却收缩阶段中热塑性和应变力学特性，推导出高弹应变模型和塑性应变模型。其次通过分析针织物定张力加热过程高弹形变时应变松弛时间、弹性模量及黏度，推导出针织物定张力加热过程的应变模型；结合热塑性和应变力学特性，进而得到加热过程针织物形变模型；最后通过分析冷却过程中的收缩率模型推导出纤维收缩形变模型，联合以上2个形变模型得到弹力织物经向与纬向热塑性应变模型。结合实际热定型工艺，由热塑性应变模型得到幅宽和面密度与热定型加工的温度、加热时间以及经、纬向拉伸量的热塑性机理模型。研究结果表明：用实际生产中的结果与模型计算结果进行对比分析，证明了该热塑性机理模型的准确性和有效性。

关键词：弹力针织物；热定型；热塑性；应变

中图分类号：TP2；TS1        文献标志码：A         文章编号：1672-7207(2019)04-0892-08

Study on thermoplastic strain model of elastic knitted fabric in the heat setting process

ZHOU Lichun, LIU Shunjing, JIN Fujiang, CHU Beibei

(College of Mechanical Engineering and Automation, Huaqiao University, Xiamen 361021, China)

Abstract: In order to obtain the strain theoretical model of the knitted fabric during heat setting process, the high elastic strain model and plastic strain models were derived firstly by analyzing the thermoplastic and strain mechanical properties of elastic knitted fabrics during normal temperature stretching, heating expansion and contraction, and cooling shrinkage. Secondly, the strain model of the constant tension heating process of knitted fabrics was derived by analyzing the strain relaxation time, elastic modulus and viscosity of high elastic deformation during the heating process of knitted fabrics. Finally, the fiber shrinkage deformation model was derived by analyzing the shrinkage model in the cooling process, and the warp and weft thermoplastic strain models of the elastic fabric were obtained by combining the above two deformation models. Combined with the actual heat setting process, the thermoplastic strain model was used to obtain the thermoplastic mechanism models of the width and surface density and heat setting process temperature, heating time and warp and weft stretch. The results show that the accuracy and effectiveness of the thermoplastic mechanism model are proved by comparing the results in actual production with the results calculated by the model.

Key words: elastic knitted fabric; heat setting; thermoplastic; strain

幅宽与面密度是弹力针织物热定型加工过程中的2个重要质量指标。建立幅宽、面密度与热定型过程中的工艺参数的模型包括时间、温度、拉幅量和超喂量等是弹力针织物热定型研究的重点^[1]，为精确设计热定型工艺参数提供依据，该模型是热定型过程的关键模型，在预测产品质量和设计工艺参数方面发挥着重要作用。目前国内外对于热定型生产过程的模型研究大部分是以实验和生产数据为基础的拟合模型。齐鲁等^[2-3]利用动态黏弹谱仪通过改变温度来研究织物经向和纬向的动力学性能，在此基础上对温度与弹性针织物中纤维的松弛模量进行建模；钟丽莹等^[4]以实验为依据阐释了织物经纬向伸缩率、热定型加热时间、温度与生产效率之间的关系，建立了可以预测面密度与幅宽的模型；郜琳琳^[5]以误差最小为原则，选取多个优化目标，建立热定型温度、时间等工艺参数与织物幅宽、面密度间的模型；赵迪^[6]根据前人的研究，为了提高质量指标面密度、幅宽预测模型的准确性采取了改进工艺参数与输出变量方案的方法，取得良好效果。但是，这些模型是在小样实验数据的基础上建立的，在实验数据的采样区间发生变化或数据有误差的条件下，模型对实际生产的作用会大大降低，因此这种模型局限性较大且缺乏普遍性。对于弹力针织物在热定型过程应变的理论模型研究，生产产品的热塑性在热定型工序中会受到影响，本文作者通过对热塑性的机理分析建立弹力针织物在热定型过程中的幅宽与面密度机理模型，并运用实际生产中的实验数据验证该模型的准确性与有效性。

1  热定型弹力织物经向和纬向变化 分析

热处理过程是利用弹性纤维的热塑性，将织物加热到一定温度，消除其内应力，同时施加恒定张力将织物幅宽缓慢拉伸至生产指定规格后撤除外力并迅速冷却降温至织物的外形与尺寸处于平稳状态而不再发生改变。热定型织物经历了常温及升温拉伸和冷却收缩3个过程，其经向和纬向发生变化。

1.1  纬向变化量

上机后针织物首先固定布边在室温下进行拉伸，调整拉幅装置而产生恒力作用于织物，对织物起到纬向拉伸的作用。其次，升温过程是通过将织物放入烘箱，对烘箱设置温度，当烘箱的温度不断升高时，织物的应力会随之降低，由于织物受到张力，其纬向将被拉伸。最后，冷却过程即织物从烘箱取出后快速降温，该过程引起织物收缩。织物的纬向在这3个过程中都会发生一定的变化，其状态如图1(a)表示。由图1(a)可知：针织物纬向发生改变过程中存在4个状态，包括上机幅宽W_i、定前幅宽和定后幅宽以及下机幅宽W_o。状态间的变量包括由拉幅装置引起的拉伸形变量、在烘箱中升温时产生的热塑性形变量△W₊和冷却时的收缩变化量△W_-。故纬向变化量△W为

            (1)

式中：△W为纬向变化量，cm；为室温下由拉幅装置引起的拉伸形变量，cm。下机幅宽W_o为

                (2)

1.2  经向变化量

经向变化过程与纬向变化过程类似，其区别在于纬向伸长的初始物理作用是通过拉幅装置提供而经向变化是由超喂装置实现。具体表现为超喂装置通过对喂布量进行控制从而在经向产生恒定张力。织物的经向伸缩变化在整个过程中的情况如图1(b)所示。故经向变化量△L为

             (3)

式中：△L为经向变化量，cm；为室温下由超喂装置引起的拉伸形变量，cm；△L₊为升温产生的热塑性经向形变量，cm；△L_-为冷却时经向收缩形变量，cm。下机布长L_o为

                (4)

式中：L_i为织物的上机布长，cm。

综上所述，经向在热定型中的变化与纬向近似，其变化趋势相同，因此，为了简化计算，只需要选择经向或者纬向一种变化建立伸缩率模型即可，以经向变化量△L为例进行研究。弹力针织物是由大量弹性纤维构成的，为了使问题简单化，假设忽略纤维与纤维之间的作用力，则可以对单根纤维的伸缩率建立模型。

据上述分析可知机械拉伸量(常温)、热塑性形变(升温) △L₊和收缩形变(冷却) △L_-这3个因素是弹性纤维在热定型工序必不可少的部分。因此，弹性纤维伸缩率模型可以分解为3部分进行探讨，以下重点研究热塑性形变与收缩形变。

图1  形变量变化过程

Fig. 1  Shape variable change processes

2  弹性纤维热塑性形变模型

2.1  常温拉伸受力分析

1) 应变与应力。应变定义为在有外力作用于物体却没有产生位移时，此外力使物体表现为局部形变。应力是物体单位截面的反作用力，它是在物体发生相对变形的情况下为了抵消外力作用而产生的^[7]。应力为矢量，可分解为法线方向和切线方向2部分，分别以正应力σ和剪应力τ来表示。

2) 应力与张力。弹性纤维在热定型过程中张力的方向与横截面相垂直，大小保持不变，纤维单位横截面积的拉伸正应力σ为

                   (5)

式中：σ为拉伸正应力，N/cm²；A为张力为F时纤维的横截面积，cm²。

3) 应力与形变。纤维受到力的作用在拉伸过程中，根据(其中，A表示横截面积，L表示纤维长度)。假定纤维体积V保持不变，横截面积和长度发生改变，即面积A_I变化成A_II，长度L_I变化成L_II时，截面积与拉伸长度成反比。根据可知当应力σ由σ_I变化成σ_II时，截面积与应力成反比。综上分析，应力与拉伸长度成正比，关系可以表示为

                  (6)

在应力表达式(5)中，当A表示热定型过程中受到张力作用下的实时横截面积，此时应力σ被称作真实应力。当A_i表示弹性纤维的初始截面积时，应力σ_a则被称为许用应力^[8]。

                  (7)

由式(6)可得

                  (8)

式中：L为瞬时布长，cm；σ_a为许用应力，N/cm²。故有

        (9)

式中：△l为实时布长变化值，cm。

4) 应变与形变。假设纤维的应变不大，那么由式(10)来描述拉伸应变即拉伸后纤维的变化量与初始长度的比值^[9]。

               (10)

式中：ε为工程应变。

5) 高弹应变。其特点是弹性模量E_g小，应变大且随时间发生变化，应变相对滞后于应力^[9]。高弹应变发生的的条件是常温下实际应力足够大或正常拉伸状态时的温度高于玻璃化温度。

              (11)

式中：τ_g为高弹形变松弛时间，s。

6) 塑性应变。塑性应变的特点是应力较小的状态下产生的应变巨大，而且应变量在外力消失后无法复原^[9]。塑性应变产生的条件是常温下实际应力较大且足以克服屈服应力或者正常拉伸时的温度高于黏流 温度。

                (12)

式中：η_s为黏度，Pa·s。

2.2  定张力加热应变模型

1) 定张力加热应变。定张力热定型过程中在纤维伸长这个宏观的状态下发生新的高弹形变。因此，其应变表达式为

     (13)

式中：E_g为实时高弹模量；为原塑性应变；为原高弹应变；为松弛回复后剩余的高弹应变；σt/η_s为在应力σ作用下新塑性应变；为在σ作用下新高弹应变^[10]。

2) 恒拉力F作用下的原应变量。弹力织物在常温下其拉伸瞬间状态满足胡克定律，拉力F可由式(7)和(10)计算。在t₀时间内织物被持续拉伸，当其大于拉伸负荷时，将产生原高弹应变：

              (14)

式中：为恒拉力下提供的原应变量；为初始弹性模量；σ₀为原高弹应力，N/cm²；t₀为持续拉伸时长，s；为初始松弛时间，s。

原塑性应变为

               (15)

式中：为原塑性应变量；σ₀为原高弹应力，N/cm²； 为初始黏度，Pa·s。

3) 加热应变高弹形变松弛时间、弹性模量及黏度。式(13)中，应变ε和应力σ的计算公式已推导，见式(8)和(9)。下面计算其余3个量即松弛时间τ_g、高弹态下的弹性模量E_g以及黏度η_s。

① 高弹形变松弛时间τ_g。τ_g的一般公式为^[10]

            (16)

式中：U为摩尔高弹活化能，J；U_p为因增塑而导致高弹活化能的降低，J；a·σ_e为应力所导致高弹活化能的降低，J；R为气体常数，R=8.307 6 J/(mol·K)；T为热力学温度，K。

令为系统下降的总活化能，则式(16)变为

             (17)

当温度波动较小时，E_a可视为常数，一般对于固体高聚物有R_a=1.6×10² J/mol^[10]。

热定型过程中由于织物受到张力作用，纤维分子的运动受到限制，因此，需要克服能量增加。故此过程中的松弛时间τ_g为

            (18)

式中：△E_a为松弛阶段需增加的能量，J。

将初始温度T₀与总活化能E_a的数值代入式(18)中，可求得。因此，将理想气体常数R，总活化能E_a及增加的活化能△E_a数值代入式(18)中可得紧张热定型过程τ_g的表达式为

             (19)

② 弹性模量。高弹形变是指在高弹状态下物质发生的形变，当弹性针织物由普弹态变为高弹态时，其弹性模量受温度影响，存在松弛现象，因此，也发生相应的转变。引入弹性模量E，在室温T₀时，

               (20)

通过单位体积内网链构象熵的改变计算应力σ 有^[11]

               (21)

式中：为交联点间网链的平均相对分子质量；ρ为密度，kg/m³。

令为温模系数，与分子性质有关故b为常数，则式(20)变为

                (22)

当式(22)中温度T为室温T₀时，与式(20)相等，故可求得

              (23)

温度保持不变的条件下拉伸弹性纤维，并维持其应变不变，应力松弛现象可以表示成线性黏弹性的Maxwell模型^[11]。该模型由1个弹簧和1个液体黏壶组成的，其表达式为

                (24)

式中：τ_m为应力松弛时间，τ_m=η/E，s；η为液体黏度，Pa·s；E为弹簧弹性模量；σ₀为模型受力形变时的初始应力，这里初始应力为σ_a，N/cm²；σ(t)为t时刻模型的内应力，N/cm²。时间为t时弹性纤维的弹性模量E(t)为

          (25)

式中：E(t)为应力松弛模量；E₀为起始弹性模量；τ_m为应力松弛时间。

由于纤维满足时温等效理论^[12]，结合式(22)和(25)可得高弹模量的最终表达式为

               (26)

式中：系数c=η_g/b；η_g为弹性纤维高弹态下的黏度，Pa·s；b为弹性纤维的温模系数。

③黏度η_s。在高弹态条件下纤维由玻璃态转变成黏流态其黏弹性保持恒定，可视为常数^[13]。计算高弹黏度η_g和黏流黏度η_s可研究黏流态。此状态的黏度与时间的关系具有以下形式：

              (27)

式中：E_η为流动活化能，J。当T在一定范围内时，η_s与T成反比。

弹性纤维处于黏流态条件下，其分子网络相互交联时的黏度η_s为常数可以表示成，其中为分子从高弹态过度到黏流态的温度。且此状态的黏度η_s与高弹态黏度η_g相等^[14]。

因此，可以分析高弹态黏度η_g在平稳状态下的求解。已知黏流态下弹性纤维的黏度η_s(T)为10³~10⁴ Pa·s^[13]，由于转折温度取温度范围内的最小值，根据黏度与温度成反比的关系故其黏度应取范围的最大值，有

             (28)

2.3  加热形变模型综合

弹性织物初始应变量的累积是从上机后到室温下的拉伸直到热定型前的拉伸过程，此过程中的布长由L_i拉伸到，初始拉伸量为。

在拉伸过程中织物的高弹应变量较小，有

            (29)

由式(7)可以计算出许用应力σ_a，由式(9)可求出原高弹应力σ₀为

         (30)

令从布匹上机后到布匹置入烘箱前的距离为s₀，若车速v₁保持不变，则拉伸时间t₀为

                  (31)

由式(20)和(29)可求得弹性纤维正常状态下的弹性模量，结合式(31)可得弹性纤维的初始高弹模量E_g0为

              (32)

将式(30)~(32)代入式(21)可得初始高弹形变松弛时间为

            (33)

此时，在拉伸力F_l的作用下，实际应力小于屈服应力，因此，初始的塑性应变为零，即

                 (34)

综上所述，将式(14)，(15)，(18)，(19)，(26)，(27)，(33)和(34)代入(19)可推导出升温过程中的热塑性形变量公式：

  (35)

式中：t为纤维热塑性阶段定型时间，t=s_t/v₁；s_t为热定型过程中烘箱的长度。

2.4  冷却收缩形变模型

假设冷却过程的温度为T_c(单位为K)，使其保持恒定，故收缩过程温差△T=T-T_c，T为定型温度。弹性纤维收缩率为C=△L_-/L_i。在松弛状态下，纤维的收缩率与定型温度和收缩时间有关，温度增高或者时间增长都会使收缩率变大，因此，纤维的收缩率如下式所示^[15]：

             (36)

式中：t_c为冷却收缩时间，t_c=s_c/v₁，s；s_c为冷却过程中的部分热定型机的长度，m；K为玻尔兹曼常数；h为普朗克常数；△E_s为冷却收缩活化能，与冷却温度有关。在冷却过程中，当温度设置为70 ℃时，△E_s=39.04 kJ/mol；当温度设置为95 ℃时，△E_s=.9 kJ/mol^[16]。

综上可知冷却收缩布长△L_-为

           (37)

3  弹性纤维伸缩率模型

将式(9)，(35)和(37)代入式(3)可得弹性纤维伸缩率模型。

        (38)

其中：

            (39)

        (40)

     (41)

             (42)

4  热定型幅宽和面密度模型

弹力针织物伸缩率模型体现为△W_f=f₁(T, V_L, t, C)以及△G_f=f₂(T, V_L, t, C)，其中，△W为门幅，△G为面密度。工艺参数包括温度T，拉幅量V_L，时间t，超喂量C。忽略织物中各纤维之间的相互作用，织物的拉伸量可以用纤维的拉伸量来表示，其模型如下^[17]。

4.1  幅宽模型

由织物的经向模型公式(38)可以推出纬向的伸缩率为

      (43)

式中：X₁，X₂，X₃和X₄见式(39)~(42)。

由式(2)可得下机幅宽W_o为

     (44)

4.2  面密度模型

面密度模型中包含上机变量和下机变量，其中上机时的变量是可测的^[18]，而下机变量可由公式求得。下机幅宽W_o由式(44)计算得到。

另外，由式(4)和(44)求得下机布长L_o为

       (45)

假设织物面密度分布均匀，且进、出布时织物总重量G_A维持恒定，因此有，故可求得面密度G_o为

                (46)

将式(43)和(44)代入式(45)得出下机面密度G_o为

       (47)

5  模型实验验证

以JC14.76T+2.22T的弹力针织物为实验对象，从数据表中随机抽取10组数据，分析真实值与计算值的误差来验证该伸缩率模型的准确性。生产数据表如表1所示。

5.1  幅宽伸缩率模型验证

测得纤维的拉伸弹性系数k_w=0.417 7 cN/mm以及初始截面积A_i=1.931 tex。通过式(44)可计算出织物下机幅宽并与真实值对比，结果如表2所示。

由表2可知：幅宽计算值与真实值间相对误差最大值为0.85%，处于允许误差范围内。求得模型中残差平方和数值较小为7.583，拟合优度R²为0.984，因此，计算值与真实值的拟合效果好且精确度高。

5.2  面密度伸缩率模型验证

测得纤维拉伸弹性系数k_l=1.526 cN/mm，根据式(47)计算出织物的下机面密度，并与生产真实值进行对比，结果如表3所示。

表1  生产数据表

Table 1  Production data table

表2  织物幅宽伸缩率模型有效性验证

Table 2  Validation of fabric width expansion ratio model

由表3可知：面密度计算值与真实值间的相对误差最大值为0.89%，在允许误差范围内。求得模型中残差平方和数值较小为7.022，拟合优度R²为0.953，因此，计算值与真实值的拟合效果好且精确度高。

表3  织物面密度模型有效性验证

Table 3  Validation of fabric gram weight expansion ratio model

6  结论

1) 通过对热定型过程中弹力织物在常温拉伸、加热伸缩、冷却收缩这3个阶段经向和纬向受力分析，得到弹力织物经向和纬向综合应变模型，最终建立了烘箱温度、超喂、拉幅量等工艺参数与幅宽、面密度质量指标间的应变机理模型。

2) 该机理模型的结构虽然复杂，但可以准确描述出热定型过程弹力织物面密度、幅宽与温度、车速、超喂和拉幅量的关系，并通过实际的热定型实验数据验证了该机理模型具有较高的准确性与有效性，为解决热定型工艺参数精确、定量设计提供了可靠的理论支持，对提高染整行业热定型生产技术水平具有重要作用。

参考文献：

[1] Yang R H, Kan C W. Effect of heat setting parameters on some properties of PLA knitted fabric[J]. Fibers and Polymers, 2013, 8(14): 1347-1353.

[2] 齐鲁, 曹振林. 涤纶织物热定型后动态力学行为[J]. 高分子学报, 1987, 1(3): 177-183.

QI Lu, CAO Zhenlin. The dynamic mechanical behavior of pet fabrics after heat-setting[J]. Acta Polymerica Sinica, 1987, 1(3): 177-183.

[3] KEINATH S E, BOYER R F. Dynamic mechanical analysis of polymers in the liquid state: use of a perforated shim stock support[J]. Journal of Applied Polymer Science, 2010, 28(6): 2105-2118.

[4] 钟丽莹, 金福江. 基于LS-SVM的棉/氨纶弹力布热定形效率预测[J]. 武汉理工大学学报(信息与管理工程版), 2010, 32(1): 88-90.

ZHONG Liying, JIN Fujiang. Heat-setting efficiency model of cotton-spandex knit goods based on LS-SVM[J]. Journal of wuhan university of technology (information & management engineering), 2010, 32(1): 88-90.

[5] 郜琳琳. 棉氨针织物热定型工艺的建模与优化设计[D]. 厦门: 华侨大学信息科学与工程学院, 2013: 13-15.

GAO Linlin. Modeling and optimization design of cotton/polyurethane knitted fabric heat setting process[D]. Xiamen: Huaqiao University. College of Information Science and Engineering, 2013: 13-15.

[6] 赵迪. 染整热定型工艺模型及数字化设计系统研究[D]. 厦门: 华侨大学信息科学与工程学院, 2015: 22-26.

ZHAO Di. Research on heat setting process model and parameters digital design system[D]. Xiamen: Huaqiao University. College of Information Science and Engineering, 2015: 22-26.

[7] 王光钦, 丁桂保, 杨杰. 弹性力学[M]. 3版. 北京: 清华大学出版社, 2015: 10-15.

WANG Guangqin, DING Guibao, YANG Jie. Elastic mechanics[M]. 3rd ed. Beijing: Tsinghua University Press, 2015: 10-15.

[8] 李维特, 黄保海, 毕仲波. 热应力理论分析及应用[M]. 北京: 中国电力出版社, 2004: 4.

LI Weite, HUANG Baohai, BI Zhongbo. Analysis and application of thermal stress theory[M]. Beijing: China Electric Power Press, 2004: 4.

[9] 杨桂通. 弹性力学简明教程[M]. 2版. 北京: 清华大学出版社, 2013: 30.

YANG Guitong. A concise guide to elastic mechanics[M]. 2nd ed. Beijing: Tsinghua University Press, 2013: 30.

[10] 沈新元.高分子材料加工原理[M]. 2版. 北京: 中国纺织出版社, 2009: 188-191.

SHEN Xinyuan. Principles of polymer materials processing[M]. 2nd ed. Beijing: China Textile Press, 2009: 188-191.

[11] 代丽君, 张玉军, 姜华珺. 高分子概论[M]. 北京: 化学工业出版社, 2005: 99-102.

DAI Lijun, ZHANG Yujun, JIANG Huajun. Introduction to polymers[M]. Beijing: Chemical Industry Press, 2005: 99-102.

[12] 何曼君, 张红东, 陈维孝, 等. 高分子物理[M]. 3版. 上海: 复旦大学出版社, 2007: 55.

HE Manjun, ZHANG Hongdong, CHENG Weixiao, et al. Polymer physics[M]. 3rd ed. Shanghai: Fudan University Press, 2007: 55.

[13] FERRY J D. Viscoelastic properties of polymers[J]. Journal of the Electrochemical Society,1980, 90(1): 95-129.

[14] 张义同. 热粘弹性理论[M]. 天津: 天津大学出版社, 2002: 15-30.

ZHANG Yitong. Thermoviscoelastic theory[M]. Tianjin: Tianjin University Press, 2002: 15-30.

[15] PREVORSEK D, TOBOLSKY A V. Determination of non-flow shrinkage ratio in oriented fibers[J]. Textile Research Journal, 1963, 17: 795-802.

[16] 上海纺织工学院化纤教研组. 高收缩共聚酯纤维收缩机理研究之一[J]. 合成纤维, 1979, 4(1): 1-4.

Shanghai Textile Institute of Chemical Fiber Teaching and Research Group. Study on shrinkage mechanism of high shrinkage copolyester fiber[J]. Synthetic Fiber, 1979, 4(1): 1-4.

[17] 张翠玲, 赵国樑, 王甜甜, 等. 热定型温度对涤纶针织物性能的影响[J]. 纺织学报, 2008, 29(9): 78-81.

ZHANG Cuiling, ZHAO Guoliang, WANG Tiantian, et al. Effects of heat setting temperature on the properties of knitted polyester fabrics[J]. Journal of Textile Research, 2008, 29(9): 78-81.

[18] 陈长荣. 氨纶弹性机织物染整工艺探讨[J]. 山东纺织科技, 1996(3): 10-13.

CHEN Changrong. Dyeing and finishing process of spandex elastic fabric[J]. Shandong Textile Science & Technology, 1996(3): 10-13.

(编辑  杨幼平)

收稿日期：2018-07-29；修回日期：2018-10-22

基金项目(Foundation item)：福建省科技计划重点项目(2017H0022)；福建省科技计划高校产学合作重大项目(2016H6015)(Project(2017H0022) supported by the Science and Technology Planning Key Program of Fujian Province; Project(2016H6015) supported by the Science and Technology Planning Major Program of Industry-University Research Cooperation of Fujian Province)

通信作者：周丽春，博士，讲师，从事染整复杂系统建模与控制研究；E-mail: 13905958730@163.com