劉宇飛 李忠海 黎繼烈 陸 俊

（中南林業(yè)科技大學食品科學與工程學院，湖南 長沙 410004）

塑料以其優(yōu)良的特性成為當今使用最為普遍的食品包裝材料，但包裝材料本身的安全性，包裝材料對食品品質的影響，越來越多的受到人們的關注。塑料不是一種純物質，它是由多種添加助劑配制加工而成，其中高分子聚合物（合成樹脂）是其主要成分，抗氧化劑、增塑劑、抗靜電劑等作為添加劑加入其中用于改善塑料加工性能或彌補合成樹脂本身性能的不足［1］。塑料中加入抗氧化劑的目的在于降低塑料在加工過程中受熱氧化，減緩塑料在儲存使用過程中由于曝露在光線下而進行的自發(fā)氧化，防止塑料變黃、變碎以延長塑料商品壽命。但是由于成品塑料中的抗氧化劑自加入之時起就在不停地反應以保證塑料品質，所以塑料食品包裝材料中會殘留有抗氧化劑反應產(chǎn)生的裂解物質以及抗氧化劑分子單體。尤其，對于聚乙烯、聚丙烯等包裝材料而言，抗氧化劑及其裂解物是食品包裝材料中遷移進入食品的主要物質［2］。這些物質進入食品必定會影響食品品質，甚至危害消費者的安全。食品塑料包裝中的這些抗氧化劑及其反應產(chǎn)物是否會進入食品，是如何進入食品，進入食品的量是多少等關于食品包裝安全性的問題正引起越來越多的國內外學者的關注和研究。文章主要綜述國內外學者對于食品接觸塑料包裝材料中常用4抗氧化劑（BHT、Irganox1010、Irganox1076、Irgafos168）向食品（模擬物）的遷移研究現(xiàn)狀。

1 食品塑料包裝材料分類

用于包裝食品的塑料，根據(jù)用途主要有塑料包裝薄膜、塑料瓶、塑料盒和塑料袋等；根據(jù)材料組成主要有聚乙烯（polyethylene，PE）、高 密 度 聚 乙 烯（high density polyethylene，HDPE）、低密度聚乙烯（low density polyethylene，LDPE）、線性低密度聚乙烯（linear low density polyethylene，LLDPE）、聚丙烯（polypropylene，PP）、聚苯乙烯（polystyrene，PS）、聚氯乙烯（polyvinyl chloride，PVC）、聚偏二氯乙烯（polyvinyl dichloride，PVDC）、聚苯二甲酸乙二醇酯（polyethylene terephthalate，PET）等共聚塑料［3］。

2 塑料包裝材料中常用抗氧劑

塑料中常加入的抗氧化劑有取代酚類、芳香族胺類、亞磷酸酯類、含硫酯類等幾種［4］，使用最廣泛的有4種：①Irganox1076（β－（3，5－二叔丁基－4－羥基苯基）丙酸正十八碳醇酯）CAS編號：0002082－79－3，是白色粉末物質，易溶于丙酮、環(huán)己烷等，不溶于水，是一種眾所周知的聚合物材料抗氧化劑，特別適用于聚烯烴材料［5］；②Irganox1010（四［β－（3，5－二叔丁基－4－羥基苯基）丙酸］季戊四醇酯）CAS編號：0006683－19－8，也是一種結晶白色粉末，由于其高性能和低擴散系數(shù)（相較于其它抗阻酚類），是大多數(shù)商品聚丙烯的熱穩(wěn)定性添加 劑［6］；③BHT（2，6－二 叔 丁 基－4－甲基苯酚）CAS 編 號：0000128－37－0，作為抗氧化劑應用比較廣泛，可直接添加在食品中作為食品添加劑能有效延遲食物酸敗的發(fā)生；④Irgafos168（三（2，4－二 叔 丁 基 苯 基）亞磷酸酯）CAS 編 號：0031570－04－4，是一種白色粉末物質，有低揮發(fā)性，不溶于水但對水解十分敏感［7］。BHT、Irganox1010、Irganox1076都屬于取代酚類抗氧化劑一般會與亞磷酸酯抗氧化劑Irgafos168共同加入塑料之中產(chǎn)生協(xié)同作用以達到最好的抗氧化效果?？傮w來看，大多數(shù)抗氧化劑穩(wěn)定作用良好而且無毒，然而，由于苯基的存在，一些苯基取代的亞磷酸酯類被認為具有一定毒性［8］；此外，一些抗氧化劑的分解產(chǎn)物，如Irgafos168的主要氧化裂解產(chǎn)物鄰苯二甲酸丁二酯，毒理學試驗［9］證明其有潛在的致癌作用，被美國和歐盟列為在部分直接與皮膚接觸化妝品禁止檢出和在兒童保育品中僅允許微量檢出。因此，包裝塑料抗氧化劑雖在許多國家都被允許使用，但其添加量和添加純度都有相當嚴格的要求。

3 抗氧化劑遷移試驗研究進展

3.1 模擬物選擇

食品種類繁多，形狀差別也很大，而且從成分本就復雜的食品里檢測分析抗氧化劑相對困難，因此為了簡便規(guī)范遷移試驗，通常會選擇能表征某一類食品的食品模擬物來替代食品進行遷移試驗。

對于遷移試驗的規(guī)定，國外法規(guī)制定較早，歐共體先后頒布過指令2004／1／EC，2004／19／EC，2005／79／EC 和2007／19／EC，這些指令經(jīng)多次修改完善，目前統(tǒng)一由（EU）No.10／2011取代；美國FDA 相關法令也對遷移試驗的條件及方法，甚至遷移單元的選取給出了較為詳盡的說明［10］；中國的GB 9687、GB 9688標準所制定的塑料遷移試驗方法主要是用于總遷移量的測定，參照歐盟標準也制定了具體的浸提試驗標準如GB／T 5009.156——2003，但總的來說中國相關的標準比較模糊，方法老舊且針對性不強。

美國FDA 依照食品特性推薦將食品分為水性、酸性、酒精類和脂肪類食品4類［10］。根據(jù)標準，這4類食品分別對應模擬物A、B、C、D，歐盟［11］對于遷移用模擬物的規(guī)定與其相近但劃分的更為詳細。特別的，對應水性、酸性和含乙醇食品可供選擇的分別是純水、低濃度的乙醇和乙酸水溶液作為模擬物，而對于油脂食物的模擬物的選取，普遍認為直接用植物油特別是橄欖油做遷移試驗最符合實際，有些研究［4，5，7］也提到純的異辛烷、正己烷和98%乙醇水溶液作模擬物來替代脂肪食物。

3.2 常用檢測方法

目前，高效反相液相色譜（HPLC）分析技術是中國對塑料包裝中抗氧劑檢測的最主要方法，有報道［8］也提到使用高效液相色譜－質譜聯(lián)用法（HPLC－MS）做檢測更為精確，這是由于質譜儀的使用能在較大質量范圍內為液相色譜提供高靈敏的物質鑒別能力。

3.3 抗氧化劑遷移試驗研究

遷移試驗研究的目的就是在試驗條件下，能近似模擬出包裝材料在實際使用中的環(huán)境條件，同時利用現(xiàn)代分析儀器，測出一個時間段的遷移變化量，分析遷移規(guī)律。遷移研究主要有總遷移（overall migration，OM）和特定遷移（specific migration，SM）兩個方面，它們分別定義食品接觸材料向食品中遷移的所有物質和某一種或一類物質。對于抗氧化劑的遷移研究主要集中于特定遷移的研究。

溫度和時間是影響抗氧化劑遷移的主要因素。J.M.López－Vilario等［12］研究了LDFE 中 的Irganox1010、Irganox1076、BHT、Irgafos168，40 ℃條件下在水性食品模擬物中遷移10d，發(fā)現(xiàn)只有Irgafos168 被檢出，且含量低于歐盟的特定遷移量。

包裝材料差異對抗氧化劑遷移也會有較大影響。Monica Bertoldo等［13］研究了水浸提兩種不同結構的PP 材料中Irganox1010，研究發(fā)現(xiàn)與PE共聚的PP材料比PP均聚物材料在相同條件下更易遷移出抗氧化劑，雖然在模擬物中使用高效液相色譜未直接檢出Irganox1010，但Irganox1010的氧化降解產(chǎn)物卻能一直檢出；M.J.Galotto等［14］對不同厚度的LDPE膜中的Irganox1076，在40 ℃，20d，脂肪和水兩類食品模擬物種的遷移試驗進行研究，結果表明材料厚度越薄，遷移出抗氧化劑所需時間越短，相同條件下脂肪模擬物中遷移出的抗氧化劑量多于水性食品模擬物。

進入20世紀，微波、輻照技術被引入食品加工領域并得到了廣泛的應用，經(jīng)過輻照、微波處理后的塑料包裝安全性成為關注的熱點。針對輻照、微波技術對于食品包裝中主要抗氧化劑的影響，Dae Hoon Jeon等［15］研究了輻照線性低密度聚乙烯（LLDPE）膜，發(fā)現(xiàn)樣品在5kGy輻照下Irgafos168并未減少，而在10kGy輻照下Irganox1076比初始濃度減少了34.9%。200kGy輻照量下，在低密度聚乙烯包裝袋中的Irgafos168并未在模擬物中檢測出，而Irganox1076 卻測出隨著輻照劑量的增加而檢測量下降的趨勢；Jonas Alin等［16］研究了微波條件下PP 中Irganox1076和Irgafos168在4種食品模擬物中遷移過程，結果顯示在研究的材料中PP均聚物是最難遷移的，在水和酸性食物模擬物中的遷移量甚微，幾乎檢測不到。

此外，李美紅［17］還對HDPE 膜中抗氧劑（Irganox1076、Irgafos168、BHT）在6種食用油中的遷移動力學過程進行了研究。

3.4 抗氧劑在遷移過程中的不穩(wěn)定性研究

由于成品塑料中抗氧劑在塑料加工使用時已經(jīng)開始與空氣接觸并發(fā)生反應，同時遷移試驗萃取液（模擬液）中的抗氧劑在試驗過程中也會反應分解影響試驗結果，因此全面掌握各種抗氧劑在食品模擬液中遷移過程的不穩(wěn)定性尤為重要。劉志剛 等［18］研究了BHT、Irganox1076 和 在20，40，60 ℃條件下，在不同食品模擬物內的不穩(wěn)定性，結果表明，除Irgafos168外的其它兩種抗氧劑在60 ℃條件下95%的乙醇水溶液中的比較穩(wěn)定性，而在50%乙醇水溶液中不穩(wěn)定性顯著增加；M.S.Dopico－García等［19］研究了Irganox1076、Irganox1010、Irgafos168、BHT 在5，40，70 ℃條件下，在4種模擬物遷移20d的變化情況，研究發(fā)現(xiàn)酚類抗氧劑在水性食品模擬物中的穩(wěn)定性要比在脂肪食品模擬物中差很多。對于水性模擬物來說，低溫狀態(tài)下水和10%濃度乙醇溶液中的抗氧劑較為穩(wěn)定；C.S imoneau等［20］研究了Irganox1076在水和油性食物及其模擬物中的穩(wěn)定性，發(fā)現(xiàn)抗氧劑在油性模擬物中的穩(wěn)定性要好于在其它模擬物，此外，溫度變化特別是高溫對于抗氧劑的不穩(wěn)定性也有相當大的影響。

4 遷移理論數(shù)學模型研究進展

遷移試驗費時費力，而且對于含量極微的遷移物質檢測相對困難，因此在試驗基礎上建立能夠廣泛應用的數(shù)學模型便成為探討物質遷移的新方向。

4.1 Fick定律的遷移數(shù)學模型

當遷移物滿足已知其在包裝材料內的擴散系數(shù)及其在包裝和食品間的分配系數(shù)時，可根據(jù)Fick第二定律的二階偏微分方程預測遷移過程。

Fick第二定律

式中：

C—— 擴散物質的體積濃度，kg／m3；

t—— 擴散時間，s；

x—— 距離，m；

D—— 擴散系數(shù)，m2／s。

假定，①初始時刻遷移物均勻地分布在聚合物中；②食品模擬物混合充分，遷移物在模擬物中沒有濃度梯度；③在遷移過程中，擴散系數(shù)D 和分配系數(shù)KP，F(xiàn)是常數(shù)；④在聚合物和溶劑界面上無傳質阻力；⑤聚合物基體不會被改變（無溶脹）［17］。

初始時，食品模擬物中遷移物含量為零，經(jīng)過某一段時間（t）后，食品模擬物中的遷移量趨于平衡值MF，t。

同時滿足以上兩個條件后，將方程（1）利用分離變量（傅里葉級數(shù)）法可求得t時刻食品中遷移量MF，t解的方程（2），此方程常用來作為單層膜遷移的模擬函數(shù)。

式中：

MF，t—— 食品模擬物中經(jīng)過時間t的遷移量，kg；

MP，0—— 遷移物在聚合物材料中的初始量，kg；

D—— 擴散系數(shù)，m2／s；

L—— 聚合物厚度，m；

VP，VF—— 塑料和食品的體積，m3；

KP，F(xiàn)—— 分配系數(shù)，即平衡時塑料與食品中的遷移物質量濃度的比值；

qn—— 方程tan（α）＝qnα的非零正根，其中n 為指數(shù)變量；

α—— 平衡時食品模擬物中的遷移物與材料中的遷移物的質量比，α＝VF／（VP·KP，F(xiàn)）。

多層材料遷移預測要考慮到多層膜的加工工藝、結構以及膜與膜之間的可能相互作用，簡單的多層材料預測模型假設大多是以單層膜假設為基礎建立的。Begley－Hollifield［22］模型，該模型是直接將2 層包裝視為一層膜遷移；Laoubi－Vergnaud［21］模型，其在2層同種材料的前提下，討論了雙層薄膜層內物質向層間遷移的過程，缺點是未考慮雙層共擠加工時抗氧化劑可能發(fā)生的遷移；Franz［23］模型針對由阻隔層和再生層組成的成分不同的雙層材料。Laoubi－Vergnaud［24］研究了3層薄膜模型，但是該模型只討論了膜間遷移卻忽略抗氧化劑在膜與食品間的作用過程。

4.2 非Fick擴散

Fick定律是用來描述物質內部擴散現(xiàn)象的，相較于Fick第一定律擴散方程只適應于擴散質不隨時間變化的遷移過程，F(xiàn)ick第二定律擴散方程則是考慮到了擴散物質流量會隨著時間和距離變化而推導產(chǎn)生的，其作為遷移預測模型比較合理。但是以上介紹的模型，均是以理想的假設條件作為前提建立的，在實際情形下，包裝與食品的相互作用遠沒有理想中的規(guī)律。實際的遷移過程還會受到食品溶液溶脹、邊界條件、食品粘稠性［25］等諸多因素的影響，在這些因素的單獨或共同作用下，物質極有可能產(chǎn)生非Fick擴散或者無規(guī)律遷移，這就使僅能符合Fick定律的簡單遷移模型顯得不再適用和有效。非Fick擴散特指一種與Fick擴散定律的線性方程結構相近，卻無法描述此種定律擴散現(xiàn)象中出現(xiàn)的不連續(xù)的一個遠離平衡的熱力學過程［26］。由于這些造成無規(guī)則擴散的影響因素眾多，人們也不了解這些因素相互之間的影響過程，所以目前還沒有適用定量分析計算方法。

4.3 影響遷移數(shù)學模型的兩個參數(shù)

由于遷移包括擴散和分配兩個過程，所以可以用擴散系數(shù)（diffusion coefficient）和分配系數(shù)（partition coefficient）來表征遷移［27］。

作為建立遷移數(shù)學模型的一個重要參數(shù)，擴散系數(shù)很難通過直接試驗測定其數(shù)值。目前國外報道較多的擴散系數(shù)預測模型有Baner－Piringer［28］模型、Limm－Hollifield［29］模型、Brandsch［30］模型和Helmroth［31］模型。在中國，擴散系 數(shù)的研究主要以介紹分析國外研究進展為主。

分配系數(shù)KP，F(xiàn)主要呈現(xiàn)形式是其定義式，其定義式：

式中：

KP，F(xiàn)—— 分配系數(shù)；

CP，∞—— 遷移達到平衡時聚合物中遷移物濃度，kg／m3；

CF，∞—— 遷移達到平衡時食品（食品模擬物）中遷移物濃度，kg／m3。

可以假設，包裝材料在與食品接觸并發(fā)生遷移的過程中，初始時食品內不含有遷移物，于是將式（3）結合質量守恒定律，得到：

式中：

CF，0——t＝0 時食品（食品模擬物）中遷移物濃度，kg／m3；

mF—— 遷移單元中食品樣品質量，kg；

mP—— 遷移單元中聚合物樣品質量，kg。

目前，國內外沒有針對食品包裝中抗氧化劑遷移的專用數(shù)學模型的報道，一般研究都是以單層數(shù)學遷移模型作為基礎，基于不同的塑料包裝材料特性建立簡單模型。Jong－Koo Han等［33］分別以3層的HDPE和LDPE包裝材料中的抗氧劑BHT 在食品模擬物中的遷移過程，利用有限元方法量化多層結構遷移，用計算機建立了一個數(shù)值處理，該模型比較成功的預測了遷移模擬結果與試驗數(shù)據(jù)。I.E.Helmroth等［34］研究了PE包裝材料中抗氧化劑Irganox1076在不同濃度乙醇水溶液中的遷移過程并建立模型，確定了擴散系數(shù)為1.1×10－13m2／s的一種Fick擴散方程對于不同遷移時間和溶液濃度條件下的遷移過程均能準確描述出來。劉志剛等［35］使用計算機軟件研究了抗氧劑BHT、Irganox1076 和Irgafos168在不同的試驗溫度下向純乙醇的遷移過程，并將計算機模擬得出的數(shù)值與試驗所得結果進行對比，分析顯示模擬計算得到的擴散系數(shù)能預測遷移過程。

5 問題及展望

遷移試驗的條件需要盡量模擬包裝實際使用的環(huán)境，試驗條件相對平和，試驗中所遷移得到的抗氧化劑極其微量，如何保存并純化所遷移得到的抗氧化劑做到定性定量分析，這就要求在試驗中必須十分注意試驗方案的制定和檢測儀器的選擇；同時，由于遷移試驗費時費力，研究遷移數(shù)學模型輔助研究的就顯得大有裨益，但是目前國內外所研究的數(shù)學模型大部分是Fick擴散定理假設基礎下的理想模型，嘗試建立適用于抗氧化劑遷移的非Fick的擴散模型可以作為研究模型建立的研究新重點?？寡趸瘎┳鳛樗芰习b材料中主要的添加助劑，其在包裝與食品間的傳質過程，遷入的抗氧化劑在食品中的變化過程以及含有遷出抗氧化劑的安全性等相關研究還有很多工作需要去完成。

1 陳宇，王朝暉，鄭德.實用塑料助劑手冊［M］.北京：化學工業(yè)出版社，2007：128～248.

2 朱文亮.食品塑料包裝材料污染物遷移的研究進展［J］.食品與機械，2010，26（6）：89～93.

3 楊潔彬，王晶，王柏琴.食品安全性［M］.北京：中國輕工業(yè)出版社，2004：179～191.

4 劉志剛，胡長鷹，王志偉.3種聚烯烴抗氧劑遷移的試驗分析及數(shù)值模擬［J］.包裝科學與工程，2007，28（1）：1～3.

5 Sanches Silva，R Sendón García，I Cooper，et al.Compilation of analytical methods and guidelines for the determination of selected model migrants from plastic packaging A［J］.Food Science＆Technology，2009，17（2 006）：535～546.

6 Comyn J.Polymer permeability［M］.New York：Elsevier，1985：723～731.

7 Demertzis P G，F(xiàn)ranz R.Systematic study on the stability of selected polymer antioxidants in EU official aqueous and alternative food simulants using HPLC［J］.Food Additives and Contaminants，1998，15（1）：93～99.

8 陳志鋒，潘健偉.塑料食品包裝材料中有毒有害化學殘留物及分析方法［J］.食品與機械，2006，22（2）：3～7.

9 Haider N，Karlsson S.Loss and transformation products of the aromatic antioxidants in MDPE film under long－term exposure to biotic and abiotic conditions［J］.Journal of Applied Polymer Science，2002，85（5）：974～988.

10 US FDA／CFSAN.Guidance for industry preparation of premarket notifications for food contact substances：Chemistry recommendations［DB／OL］.Code of Federal Regulations，Title 21 Food and Drugs.2011［2012－9－15］.http：／／www.fda.gov／Food／GuidanceComplianceRegulatoryInformation／Guidance－Documents／FoodIngredientsandPackaging／ucm081818.htm.

11 EU.No.10／2011On plastic materials and articles intended to come into contact with food［S］.Brussels：Official Journal of the European Union，2011.

12 Dopico－García M S，López－Vilario J M，González－Rodríguez M V.Determination of antioxidant migration levels from low－density polyethylene films into food simulants［J］.Journal of Chromatography A，2003，10（18）：53～62.

13 Monica Bertoldoa，F(xiàn)rancesco Ciardellib.Water extraction and degradation of a sterically hindered phenolic antioxidant in polypropylene films［J］.Polymer，2004，45（36）：8 751～8 759.

14 Galotto M J，Torres A，Guarda A，et al.Experimental and theoretical study of LDPE versus different concentrations of Irganox1076and different thickness［J］.Food Research International，2010，44（11）：566～574.

15 Dae Hoon Jeon，Gun Young Park，In Shin Kwak.Antioxidants and their migration into food simulants on irradiated LLDPE film［J］.LWT，2007，40（12）：151～156.

16 Jonas Alin，Minna Hakkarainen.Type of polypropylene material significantly influences the migration of antioxidants from polymer packaging to food simulants during microwave heating［J］.Journal of Applied Polymer Science，2010，118（19）：1 084～1 093.

17 李美紅.微波作用下塑料包裝材料中化學物遷移動力學研究［D］.廣州：暨南大學，2008.

18 劉志剛，王志偉.3種聚烯烴抗氧劑在不同食品模擬物內的不穩(wěn)定性試驗研究［J］.包裝工程，2008，29（4）：1～3.

19 Dopico－GarcáM S，López－Vilarino J M，González－Rodíguez M V.Effect of temperature and type of food simulant on antioxidant stability［J］.Journal of Applied Polymer Science，2006，100（59）：656～663.

20 Simoneau C，Hannaert P.Stability testing of selected plastics additives for food contact in EU aqueous，fatty and alternative simulants［J］.Food Additives and Contaminants，1999，16（5）：197～206.

21 王平利，王志偉，胡長鷹.遷移預測模型中擴散系數(shù)的研究［J］.包裝工程，2009，30（1）：11～14.

22 Crank J.The mathematics of diffusion［M］.2nd ed.Oxford：Clarendon Press，1975：129～135.

23 Lauo，Wong S.Contamination in food from packaging material［J］.Journal of Chromatography A，2000，882（1～2）：255～270.

24 Laoubis，Vergnaud J M.Food sandwich packaging with a recycled polymer between two functional barriers of different thicknesses［J］.Polymer Testing，1996，15（3）：269～279.

25 Grassi，Colombo，Lapasin.Experimental determination of the theophylline diffusion coefficient in swollen sodium－alginate membranes［J］.Journal of Controlled Release，2001，76（8）：93～105.

26 黃肖紅.HDPE膜中抗氧化劑1076向食品或食品模擬物的遷移研究［D］.廣州：暨南大學，2008.

27 Begl Eyth，Brandschj，Limmw，et al.Diffusion behaviour of additives in polypropylene in correlation with polymer properties［J］.Food Additives and Contaminants，2008，25（11）：1 409～1 415.

28 Baner A L，Brandsch J，F(xiàn)ranz R，et al.The estimation of migration potential of additives from polymeric materials into foodstuffs［C］／／In Proceedings of the 17th Annual International Conference in Stabilization and Degradation of Polymers.Luzern：Fraunhofer Publica，1995：11～30.

29 Limmw，Hollifieldhc.Modelling of additive diffusion in polymers［J］.Food Additives and Contaminants，1996（13）：949～967.

30 Masaro L，Zhu X X.Physical models of diffusion for polymer solutions，gels and solids［J］.Progress Polymer Science，1999，24（12）：731～775.

31 Helmrothe，Varekampc，Dekkerm.Stochastic modeling of migration from polyolefins［J］.Journal of the Science of Food and Agriculture，2005，85（18）：909～916.

32 王志偉，孫彬青，劉志剛.包裝材料化學物遷移研究［J］.包裝工程，2004，25（5）：1～4.

33 Jong－Koo Han，Susan E Selke，Theron W Downes，et al.Application of a computer model to evaluate the ability of plastics to act as functional barriers［J］.Packing Technology and Science，1999，36（9）：231～238.

34 Helmroth I E，Bekhuis H A M，Linssen J P H，et al.Direct measurement of additive migration from low－density polyethylene as a function of space and time［J］.Journal of Applied Polymer Science，2002，86（10）：3 185～3 190.

35 劉志剛，胡長鷹，王志偉.3種聚烯烴抗氧劑遷移的試驗分析及數(shù)值模擬［J］.包裝工程，2007，28（1）：26～33.