梁衛(wèi)華，吳大勇，舒均國

（1中國科學(xué)院理化技術(shù)研究所，北京100190；2湖南中鋰新材料有限公司，湖南 常德415001）

20 年來，鋰離子電池逐步成為應(yīng)用最廣泛的儲能器件，為消費(fèi)電子產(chǎn)品和電動汽車提供動力。中國已建成完整的鋰離子電池產(chǎn)業(yè)鏈，根據(jù)中國電子顯示：2019 年中國鋰離子電池產(chǎn)量為157.2 億只，同比增長12.4%；按容量計(jì)算，鋰離子電池產(chǎn)量148.0 GW·h，同比增長19.2%。其中，動力電池產(chǎn)量為85.4 GW·h，同比增長21%；消費(fèi)型電池產(chǎn)量56.2 GW·h，同比增長15.4%[1]。中國已經(jīng)成為世界上最大的鋰離子電池生產(chǎn)國，且市場需求仍在不斷擴(kuò)大[2]。雖然我國已經(jīng)在鋰電池制造工藝和裝備制造業(yè)方面都取得了長足的進(jìn)步，但是在涉及卷對卷制造工藝領(lǐng)域的基礎(chǔ)研究幾乎是停滯的。在鋰電池極片生產(chǎn)過程中應(yīng)用卷對卷工藝包括涂布、輥壓、分切、卷繞等工序。其中涂布工藝，涉及涂布漿料的均質(zhì)與輸運(yùn)、流體力學(xué)和表面張力控制下的潤濕與鋪展、漿料的薄膜干燥與固化等化工單元過程和多臺電機(jī)電氣控制和基膜張力控制等單元控制過程。因此涂布工序被認(rèn)為是鋰電池生產(chǎn)中最為復(fù)雜和重要的工序。作為生產(chǎn)極片的工程技術(shù)人員，即使并不需要深入了解各種涂布設(shè)備制造的技術(shù)細(xì)節(jié)，至少也需要掌握各種涂布參數(shù)的意義和設(shè)置方法。但是，各種關(guān)于涂布技術(shù)的專業(yè)書籍和文章，由于涉及流體力學(xué)，閱讀起來深奧難懂，致使我國產(chǎn)業(yè)界理論與技術(shù)水平和美國、日本等存在巨大差距。

涂布有多種方式，如順輥涂布、逆輥涂布、刮刀涂布、擠壓涂布、坡流涂布、凹版涂布等?；谕坎紳{料的性質(zhì)和涂布厚度精度的要求，往往可以選出2～3 種涂布方法。例如，適用于電池極片涂布的方法有刮刀涂布、擠壓涂布等。條縫擠壓涂布是一種重要的涂布方式，其涂布量由計(jì)量泵設(shè)定，對涂布厚度/面密度的控制比較精確[3]，生產(chǎn)效率也很高。因此，條縫擠壓涂布正逐漸替代逗號刮刀轉(zhuǎn)移涂布，成為鋰離子電池極片制造的主流方法。但是，條縫涂布也有不足，即高速工作時(shí)不能進(jìn)行間歇涂布留出焊接極耳的空白。逗號刮刀轉(zhuǎn)移涂布是另一種應(yīng)用廣泛的極片生產(chǎn)方式[4]（圖1），它的結(jié)構(gòu)簡單，制造成本較低，清理方便，特別適用于極片的分段涂布，因此在消費(fèi)電子型鋰電池極片制造和超級電容的極片制造中，仍被普遍采用。但是在鋰電池極片涂布生產(chǎn)過程中，操作工人通常需要憑借經(jīng)驗(yàn)多次調(diào)整涂布縫隙H0來獲得需要的涂布厚度和面密度，難以避免的試涂工作浪費(fèi)了大量的物料并降低了生產(chǎn)效率。

圖1 (a)逗號刮刀轉(zhuǎn)移涂布；(b)涂布縫隙放大圖及坐標(biāo)系Fig.1 (a)Schematic diagrams of comma-roll coating;(b)enlarged view of coating gap and coordinate system

從20世紀(jì)80—90年代起，國際上涂布技術(shù)的研究主要集中在順流對輥涂布和逆流對輥涂布方式上[5-6]，而關(guān)于逗號刮刀涂布的理論分析和實(shí)驗(yàn)研究較少。近年來，包能勝等[7-8]從雷諾潤滑方程出發(fā)，結(jié)合涂膜的成型機(jī)理，建立雙輥涂布的理論模型求解涂布厚度并與流體仿真結(jié)果進(jìn)行了對比得到比較吻合的結(jié)果，并預(yù)測了雙輥涂布的涂布厚度。劉正士等[9]從材料力學(xué)的角度出發(fā)計(jì)算了流場壓力引起刮刀變形對涂布橫向精度的影響。不過，該研究沒有求解關(guān)于涂布最關(guān)心的問題——涂布厚度/面密度。本文作者認(rèn)為，非常有必要針對逗號刮刀涂布建立模型，開展理論分析并結(jié)合流體力學(xué)專業(yè)軟件進(jìn)行數(shù)值模擬；探究刮刀輥、涂布輥徑、縫隙高度等上游形狀因素對流場和涂布最終厚度影響，導(dǎo)出涂布濕厚度與涂布縫隙的關(guān)系式，將模擬數(shù)值與理論值、實(shí)際涂布結(jié)果相對照，從而揭示涂布濕厚度與涂布刮刀縫隙之間的關(guān)系，提高涂布工藝的效率和水平。

根據(jù)表觀黏度與剪切速率的關(guān)系，常把涂布漿料分成3種：剪切變稀型流體、剪切變稠型流體和牛頓型流體。鋰離子電池極片涂布漿料是復(fù)雜的混合體系，包括溶劑、黏合劑、活性材料和導(dǎo)電添加劑等。它通常表現(xiàn)為非牛頓剪切變稀材料，其流變特性包括黏度剪切變稀、黏彈性、屈服應(yīng)力和時(shí)間依賴性[10]。涂布過程是一個連續(xù)穩(wěn)態(tài)過程，漿料的黏彈性、屈服應(yīng)力等對涂布的影響較小，可忽略不計(jì)。漿料加工設(shè)備和工藝條件會影響漿料的黏度。當(dāng)漿料所含黏結(jié)劑較少時(shí)，呈現(xiàn)牛頓型流體特性；由球磨機(jī)和3D 混合器制備的漿料也呈現(xiàn)牛頓型流體特性[10-11]。由于漿料體系實(shí)際情況比較復(fù)雜，所以在涂布研究和涂布裝置的設(shè)計(jì)中往往將漿料視為均質(zhì)流體。對于非牛頓流體的漿料也可以在一定的剪切率范圍內(nèi)將其近似為牛頓型流體，取特征剪切率下的黏度進(jìn)行分析和計(jì)算。這種處理方法在文獻(xiàn)中已有報(bào)道[3,8,12-13]。

1 刮刀涂布理論分析與求解

1.1 刮刀涂布數(shù)學(xué)模型的推導(dǎo)

逗號刮刀轉(zhuǎn)移涂布如圖1(a)所示，常用于鋰電池極片涂布。涂布輥以速度U＊順時(shí)針轉(zhuǎn)動將涂布漿料帶出刮刀縫隙H0在涂布輥上形成穩(wěn)定涂布濕厚度Hw；橡膠背輥順時(shí)針轉(zhuǎn)動并與涂布輥負(fù)間隙接觸，將漿料完全轉(zhuǎn)移到鋁箔或銅箔上。上下移動刮刀可以改變涂布縫隙H0，從而改變涂布最終濕厚度Hw。如何確定H0與Hw的關(guān)系，是技術(shù)人員十分關(guān)心的重要問題。

特定漿料涂布的特性是流體力學(xué)研究的范疇。納維-斯托克斯方程是流體力學(xué)的基本方程

式中，u為流體在x方向的速度；v為y方向的速度；p為流體靜壓力；μ和ρ分別為流體的黏度系數(shù)和密度。

根據(jù)遠(yuǎn)場理論，遠(yuǎn)離刮刀縫隙處的流場流動對刮刀縫隙處的流動影響很小，可以忽略不計(jì)。因此將研究重點(diǎn)放在逗號刮刀涂布縫隙附近的區(qū)域，建立以逗號刮刀圓弧與垂直線交點(diǎn)為坐標(biāo)原點(diǎn)、坐標(biāo)體系如圖1(b)所示的研究區(qū)域。在該區(qū)域內(nèi)，即刮刀-涂布輥縫隙范圍內(nèi)，按雷諾潤滑理論Y 方向的速度很小可以近似v=0，從式(3)可得?u/?x=0、?2u/?x2=0，代入式(1)、式(2)中，得到

式(4)和式(5)即為一維形式的雷諾潤滑方程，由此兩式可知流場壓力只關(guān)于x 的一維函數(shù)p(x)，速度u是只關(guān)于y的一維函數(shù)u(y)。

對式(4)進(jìn)行y坐標(biāo)方向上的一次積分

二次積分為

式中C1和C2未知。這時(shí)需要引入邊界條件：刮刀圓弧處速度為零即y=0 時(shí)u=0；涂布輥速度為U＊，即y=H(x)時(shí)u=U＊，U＊為涂布速度。

式中H(x)為刮刀與涂布輥之間的縫隙高度，是關(guān)于x 的形狀函數(shù)。在圖1(b)坐標(biāo)體系下，定義刮刀縫隙最小高度為H0，則刮刀與涂布輥之間的縫隙高度H(x)為

式中當(dāng)量半徑R=RcRr/(Rc+Rr)，式(7)是關(guān)于y的顯函數(shù)，同時(shí)也是關(guān)于x 的隱函數(shù)，因?yàn)槭街衐p/dx和H(x)都是關(guān)于x的函數(shù)。目前來看，dp/dx關(guān)于x的函數(shù)是未知的。此時(shí)需要引入系統(tǒng)的質(zhì)量守恒定律，假設(shè)Q是涂布漿料在單位寬度上的體積流量，那么沿x方向的流量是恒定的都為Q。

式中，Hw為最終涂布濕厚度，未知。將式(7)代入式(10)中

至此，式(12)推導(dǎo)出了基于雷諾潤滑理論的刮刀縫隙涂布理論模型。

1.2 理論模型的求解

結(jié)合邊界條件例如壓力導(dǎo)數(shù)邊界條件或壓力邊界條件和幾何形狀函數(shù)H(x)可以對上述模型進(jìn)行求解，最終求出最終涂布厚度Hw。如圖1所示，刮刀和涂布輥的左側(cè)安裝有涂布料槽。在進(jìn)行刮刀涂布的過程中，料槽內(nèi)充滿涂布漿料，并且人工或液面跟蹤系統(tǒng)保持涂布液面穩(wěn)定。又因?yàn)樵谑?1)～式(3)中忽略重力因素對流場的影響，可以假設(shè)在-∞≤x≤0區(qū)域內(nèi)充滿涂布漿料。那么在x=-∞處，壓力為外界空氣壓力，即

在x=0處，壓力為外界空氣壓力。

將式(9)代入式(12)整理得

對式(16)進(jìn)行積分，

在上述邊界條件中，假設(shè)p(-∞)=0，顯然這個假設(shè)與涂布料槽的實(shí)際情況不符。那么提出另外一種假設(shè)，邊界條件為x=-R 當(dāng)量半徑處的壓力為零，即

其中x=-R=-RcRr/(Rc+Rr)，該點(diǎn)所處流場的位置見附錄1。

為簡化起見，令R/2H0=K，這是一個關(guān)于流場輪廓幾何形狀的參數(shù)，得

至此推導(dǎo)出了基于雷諾潤滑方程理論的逗號刮刀幾何形狀參數(shù)對涂布厚度影響的關(guān)系式(23)。在已知刮刀輥、涂布輥當(dāng)量半徑R和刮刀縫隙H0時(shí)，可以通過該式計(jì)算最終涂布濕厚度Hw。另外在計(jì)算出Hw/H0后，可以得到流場壓力隨位置θ變化的函數(shù)

其中式(23)與文獻(xiàn)[14]中的刮刀涂布厚度公式特定條件下（刮刀傾斜角度為0）得到的公式相同，只是我們的推導(dǎo)過程主要針對逗號刮刀，而不是文獻(xiàn)[14]中的通用型刮刀。

2 基于Fluent 軟件的逗號刮刀流場模擬

上述的刮刀涂布理論分析與求解是建立在刮刀縫隙附近小區(qū)域內(nèi)，流場符合潤滑理論假設(shè)的基礎(chǔ)上，并在設(shè)定p(0)=0和p(-R)=0的邊界條件下求解得出的。進(jìn)而，我們希望利用專業(yè)軟件建立流場范圍更廣、更接近逗號刮刀涂布流動實(shí)際情況的2D模型；通過設(shè)定適當(dāng)?shù)倪吔鐥l件，用CFD 方法求解流場模型。可以通過兩者結(jié)果的比較來判斷刮刀涂布理論邊界條件是否合適，并分析前者的適用范圍。因此，采用Fluent軟件進(jìn)行了逗號刮刀流場模擬。對逗號刮刀流場進(jìn)行流體力學(xué)數(shù)值模擬計(jì)算，首先使用Fluent 軟件的網(wǎng)格模塊Gambit 建立漿料料槽和刮刀涂布輥組成的2D 模型，并定義邊界條件。如圖2涂布輥邊界條件定義為moving wall，逗號刮刀邊界條件定義為wall，料槽頂端邊界條件設(shè)為Pressure-inlet，刮刀縫隙出口的邊界條件設(shè)為Pressure-outlet，其他的邊界設(shè)為wall 條件。因?yàn)椴豢紤]重力對涂布厚度的影響，所以在操作條件中未激活重力因素選項(xiàng)。

圖2 刮刀涂布流場邊界條件Fig.2 Physical model of 2D and boundary conditions

其次，是2D 模型網(wǎng)格的劃分。以涂布輥半徑Rr=20 mm，刮刀半徑Rc=20 mm，涂布縫隙H0=0.2 mm 為例，采用Submap 方式，單元為Quad型，對2D 區(qū)域進(jìn)行網(wǎng)格劃分。其中，在壓力出口處進(jìn)行了細(xì)化，劃分?jǐn)?shù)為40，結(jié)果如圖3所示。其中Cell 數(shù)量22400 個，F(xiàn)ace 45200 個，Node 數(shù)量22841 個。最小網(wǎng)格面積5.0×10-6m2，最大網(wǎng)格面積5.45×10-4m2。網(wǎng)格驗(yàn)無關(guān)性驗(yàn)證的工作證實(shí)在壓力出口處劃分大于等于40，整體面網(wǎng)格大于等于45200個后壓力出口處的流量已經(jīng)不隨網(wǎng)格數(shù)量的變化而發(fā)生變化了。

再者，設(shè)定涂布漿料密度ρ=1000 kg/m3，漿料黏度μ=5 Pa·s，涂布速度U＊=0.1 m/s，料槽頂端邊界條件設(shè)為Pressure-inlet，刮刀縫隙的邊界條件設(shè)為Pressure-outlet，兩者的數(shù)值都設(shè)為0。計(jì)算該區(qū)域的雷諾數(shù)Re=H0U＊ρ/μ=0.004，因此在Fluent V6.3.26 軟件中，采用層流模型，忽略重力因素，用有限體積法求解控制方程，選擇采用壓力隱式穩(wěn)態(tài)定常求解器?？刂品匠虊毫﹄x散采用標(biāo)準(zhǔn)式，動量離散采用一階迎風(fēng)格式，壓力速度耦合方程組采用Simple 算法及欠松弛因子迭代求解。

迭代次數(shù)及收斂標(biāo)準(zhǔn)：本算例中采用，連續(xù)方程和速度方程的殘差要求控制在小于10-6和10-7迭代次數(shù)為6613 次。計(jì)算區(qū)域和刮刀縫隙區(qū)域流函數(shù)分布如圖4所示，全流場內(nèi)的壓力云圖和刮刀縫隙區(qū)域壓力云圖如圖5所示。在遠(yuǎn)離狹縫的大部分區(qū)域內(nèi)壓力分布在-1.20×10-3～1.14×10-3Pa之間，近似為0。僅在距狹縫出口5 mm 范圍內(nèi)，壓力先顯著上升再急劇下降，表現(xiàn)為只與x相關(guān)的一維函數(shù)p(x)幾乎不隨y 變化。這與前文雷諾潤滑方程的推論相同。需要指出的是，在圖5中負(fù)壓出現(xiàn)在涂布輥剛剛進(jìn)入料槽處，液體速度急劇上升引起壓力下降，成為全流場靜壓最低點(diǎn)。壓力出口處的速度沿y軸分布如圖6所示。

圖3 刮刀涂布流場的網(wǎng)格劃分Fig.3 Grid of comma-roll coating model

圖4 逗號刮刀流場的流函數(shù)Fig.4 Stream function distribution

圖5 刮刀縫隙處的壓力分布Fig.5 Pressure distribution at coating gap

圖6 縫隙出口截面（x=0）速度曲線Fig.6 Velocity distribution curve at x=0

最后，根據(jù)計(jì)算結(jié)果中涂布縫隙H0處輸出的質(zhì)量流量m=0.013280711 kg/s，計(jì)算涂布濕厚度與涂布縫隙的比值Hw/H0=m/ρU＊H0=0.6640。

采用與以上算例類似的方法和過程，分別用軟件模擬了不同輥徑、縫隙高度、黏度和速度下流場情況，并計(jì)算Hw/H0數(shù)值，結(jié)果列入表1、表2和表3。

3 結(jié)果與討論

3.1 刮刀縫隙帶出率Hw/H0與幾何形狀參數(shù)K 的關(guān)系

通過對逗號刮刀涂布流場雷諾潤滑理論建模和求解，得到了關(guān)于刮刀輥、涂布輥輥徑、縫隙高度等上游形狀因素對流場和涂布最終厚度濕厚度的關(guān)系式(23)。求解公式(23)，得到了K與Hw/H0的函數(shù)關(guān)系（圖7），并將表1 中流體模擬結(jié)果Hw/H0引入圖7。結(jié)果表明，逗號刮刀涂布的濕厚度Hw與刮刀縫隙高度H0呈正比關(guān)系，同時(shí)受到上游當(dāng)量半徑與縫隙高度比值R/2H0=K 影響。K 越大，Hw/H0的比值越大。當(dāng)K 趨向于+∞，即涂布縫隙趨近于0 時(shí)，Hw/H0→2/3。這意味著牛頓流體的逗號刮刀流場流出率的極限值Hw/H0=2/3。

表1 輥徑大小和縫隙高度對涂布厚度影響Table 1 Simulation results of different geometric parameters

表2 漿料黏度對涂布厚度的影響Table 2 Simulation results of different viscosity

表3 速度對涂布厚度的影響Table 3 Simulation results of different velocity

圖7 Hw/H0關(guān)于幾何參數(shù)K的函數(shù)圖Fig.7 Comma roll coater flow rate predicted by lubrication theory,with comparison to numerical model

由于在實(shí)際應(yīng)用中K=R/2H0的取值范圍一般在10～2000范圍內(nèi)，所以Hw/H0的取值范圍，多在略小于2/3 的0.66 附近。圖7 所示，潤滑理論計(jì)算結(jié)果和表1 中2D 流體力學(xué)軟件模擬計(jì)算的結(jié)果在K值（10～2000）范圍內(nèi)高度吻合，說明逗號刮刀涂布模型是合理的。涂布工程師和操作技術(shù)人員可以使用該公式來設(shè)定合適涂布間隙H0來進(jìn)行第一次試涂。

3.2 流場中的壓力分布

圖8 不同K值時(shí)流場壓力與位置關(guān)系圖Fig.8 Relationship between flow field pressure and position at different K values

理論推導(dǎo)過程中邊界條件p(-R)=0的假設(shè)是否成立，可以通過軟件計(jì)算來驗(yàn)證。在表1不同幾何參數(shù)的軟件模擬結(jié)果壓力分布云圖（圖5）中取坐標(biāo)(-R,0)點(diǎn)的壓力p(-R,0)，代入式(25)，考察此處的無因次壓力P＊(-R,0)，結(jié)果列于表4。可見，除了K=5.55的無因次壓力P＊(-R,0)=0.0154外，在其他K值條件下P＊(-R,0)與P＊max相比都可以近似為零。因此，p(-R)=0的邊界條件假設(shè)基本合理。

表4 不同幾何參數(shù)K條件下的-R處的無因次壓力Table 4 Dimensionless pressure at point of(-R,0)

3.3 漿料黏度和涂布速度對涂布厚度的影響

由式(23)可知，作為牛頓流體涂布漿料的黏度對涂布厚度沒有影響。這一點(diǎn)通過Fluent軟件模擬也得出相同的結(jié)論，表2 中漿料黏度在0.1～10 Pa·s 的變化，引起的厚度變化不到千分之一，說明漿料黏度變化對厚度幾乎沒有影響。同樣地，在Fluent模擬結(jié)果中（表3）涂布速度從0.005 m/s變化到0.4 m/s，涂布厚度只有萬分之一的變化。這表明速度對涂布厚度沒有明顯的影響。在式(11)中涂布流量與涂布速度U＊和dp/dx有關(guān)

把式(16)代入(11)中，則可以消除黏度項(xiàng)μ。因此，涂布厚度與漿料黏度并不相關(guān)。又因?yàn)镠w=Q/U＊，在計(jì)算過程中又消除了涂布速度U＊。因此，涂布厚度也不受涂布速度大小的影響。

3.4 重力作用產(chǎn)生的水力壓頭對涂布厚度的影響

在理論模型的建立和求解中忽略了重力作用對涂布厚度的影響。這是因?yàn)楣蔚锻坎挤较驗(yàn)樗椒较颍c重力方向不一致，其水力壓頭沿水平方向x的偏微分?pg/?x不能求出解析解。為了進(jìn)一步分析重力作用對涂布厚度的影響，將上述公式中的靜態(tài)壓力p 由其靜壓和水頭壓力組成的等效壓力pe=p+pg代替，其中在刮刀縫隙處x 軸的水頭壓力pg=ρgRc，Rc為刮刀輥半徑?？紤]重力因素后，式(17)的邊界條件應(yīng)該改為p(θk)=ρgRc和p(0)=0。求解該方程可得

式(27)中第1 項(xiàng)與式(23)相同，是未考慮重力時(shí)的涂布帶出率；第2 項(xiàng)為狹縫進(jìn)出口水頭壓力差對涂布厚度的影響，定義為Δ Hg。因?yàn)棣萲=arctan( - R 2H0)∈（-π/2,0)，所以第2 項(xiàng)分母為負(fù)值。涂布厚度在考慮水力壓頭時(shí)會比不計(jì)水力壓頭時(shí)有所增大。根據(jù)式(27)計(jì)算如表1 的各參數(shù)時(shí)水力壓頭對涂布厚度的影響ΔHg，見表5。

同理也可以計(jì)算出表3 和表4 條件下水力壓頭對涂布厚度的影響，見表6。

表5 不同幾何參數(shù)K條件下的水頭壓力產(chǎn)生對涂布厚度影響Table 5 Influence of hydraulic head on coating thickness at different K

表6 不同黏度或涂布速度時(shí)，水力壓頭對涂布厚度的影響Table 6 Influence of hydraulic head on coating thickness at different viscosity and speed

觀察表5和表6發(fā)現(xiàn)只有在幾何參數(shù)K=11.11、5.555 即刮刀間隙很大時(shí)（H0=1 mm、2 mm）水力壓頭才會對涂布厚度產(chǎn)生0.017和0.047的增量，其他縫隙小于等于0.5 mm 情況下，黏度μ∈[0.1,10] Pa·s 和涂布速度U＊∈[0.005, 0.4] m/s 范圍內(nèi)，水力壓頭對涂布厚度的影響可以忽略不計(jì)。通過改變Fluent軟件操作條件增加重力因素的方法也可以模擬水力壓頭對涂布厚度的影響，得到的結(jié)果也演證了式(27)和表5、表6 的結(jié)果。由于篇幅所限，不再贅述。分析表明，常用鋰電池極片刮刀涂布工況下重力產(chǎn)生的水力壓頭對涂布厚度的影響可以忽略不計(jì)。

3.5 表面張力對逗號刮刀涂布厚度的影響

圖9 刮刀縫隙出口自由液面形狀Fig.9 Free surface after coating gap

在式(14)中，假設(shè)壓力出口（x=0）處涂布漿料的壓力p(0)等于外界大氣的壓力，即表壓p(0)=0，是雷諾潤滑方程邊界條件的常用設(shè)置。在使用Fluent 軟件設(shè)置邊界條件時(shí)也可以把出口壓力設(shè)為0。但在實(shí)際涂布工作中，漿料流出壓力出口后形成了由液體表面張力、黏性力和液體壓力的合力控制自由表面流動。Coyne經(jīng)過研究得出了涂布自由液面高度h(x)與濕厚度hw之比的關(guān)系式h(x)/hw=1+kexp(-λx)[11-12]。如圖9所示，在x=+∞處h(x)=hw，p(+∞)=0；在x=0 處p(0)=-σ/r＜0。其中，σ 為涂布漿料的表面張力，r 為x=0 處的液面彎曲半徑。雖然σ 可以測量，但r 是未知量，所以p(0)也是未知量。因此，能否簡單地把p(0)設(shè)置為0來應(yīng)用，需要按不同的條件進(jìn)行討論。Kistler 等[14]認(rèn)為，在刮刀（剛性刮刀和片式刮刀）涂布流場中設(shè)置入口和出口壓力為零的潤滑理論，可以準(zhǔn)確地預(yù)測涂布厚度和液體壓力。Kawabata[15]測量了高黏度漿料（μ=4～6 Pa·s）刮刀流場壓力，證實(shí)刮刀縫隙出口處壓力為0。Sullivan等[16]通過有限元法數(shù)值計(jì)算和實(shí)驗(yàn)證實(shí)，刮刀縫隙出口處p(0)的數(shù)值與峰值壓力相比小至可以忽略不計(jì)；在毛細(xì)準(zhǔn)數(shù)不是很大（Ca=μU/σ=3.7和0.9）時(shí)，表面張力的影響也可以忽略不計(jì)。通常，鋰電池極片涂布的黏度在1～10 Pa·s范圍，水性負(fù)極漿料表面張力0.04～0.05 N/m[3]，NMP 漿料體系的表面張力約0.04 N/m（20 ℃）。逗號刮刀涂布速度為0.05～0.2 m/s，由此可以計(jì)算出電池漿料的毛細(xì)準(zhǔn)數(shù)Ca＞1。毛細(xì)準(zhǔn)數(shù)Ca體現(xiàn)了液體黏性力與表面張力之比，Ca 越大，表面張力對流動的影響也越小。在這種情況下，刮刀縫隙出口處的壓力為零的假設(shè)成立。在鋰電池極片涂布速度較高的工況下，表面張力對涂布厚度的影響可以忽略不計(jì)。

4 結(jié) 論

本文通過對N-S方程的合理簡化，得出刮刀縫隙附近流場的雷諾潤滑方程。應(yīng)用該方程在合理邊界條件下求解逗號刮刀涂布流場適用于牛頓流體的力學(xué)行為發(fā)現(xiàn)，盡管逗號刮刀上游幾何尺寸差距明顯，但是它的無因次壓力-位置曲線幾乎重合；涂布濕厚度Hw與刮刀縫隙高度H0的比值在K=R/2H0的取值范圍（10～2000）內(nèi)約等于2/3，而且不受涂布漿料黏度變化和涂布速度變化的影響。通過刮刀涂布理論模型的得出的涂布濕厚度Hw以及壓力-位置關(guān)系等計(jì)算結(jié)果與Fluent軟件模擬的二維涂布流場的CFD結(jié)果高度吻合。

圖A1 (-R,0)處于流場中的位置
Fig.A1 Geometrical parameters of comma-roll coating system