中性筆出墨量與測試條件參數之間的數學模型構建

吳瑛戟

上海文采實業有限公司，上海 200072

自從中性筆發展以來，在其各項檢測指標中出墨量一直是技術人員關注的焦點，因為出墨量的大小直接影響到書寫的體驗。出墨量太大或太小對書寫感覺都有不利的影響。出墨量大，書寫手感比較順滑，但筆頭容易積墨，書寫線跡易有積墨點，字跡干燥時間相對變長，墨水消耗快。對于有些滲透性強的墨水，紙張背面易滲透墨跡，嚴重的會致使背面無法書寫利用。出墨量太小，書寫手感較澀，易造成字跡不飽滿、空心，嚴重的甚至會斷線，另外書寫過程對筆頭球座的磨損量也會相對較大，總體書寫的體驗較差。

然而要準確測試中性筆的出墨量并不是件容易的事，因為中性筆出墨量測試過程中，影響出墨量大小的因素較多。這些因素包括：書寫角度、書寫載荷、書寫間距、書寫紙性能、書寫速度等。所以如何準確測定中性筆的出墨量一直是行業內專家討論和研究的熱門課題之一。從以往的研究和發表的文獻來看，姚鶴忠[1-2]從圓珠筆筆頭結構參數和墨水的性能參數與出墨量關系進行了研究。陳基銘[3]考慮了墨水、筆頭、紙張、球珠表面的墨水層等因素，構建了墨水與筆頭匹配的出墨量仿真模型，利用計算機軟件進行模擬分析。王通虎等[4]研究了測試時劃線間距變化對中性筆出墨量的影響。吳瑛戟[5]等從中性筆出墨量測量不確定度評定的角度研究了測試中的主要影響因素。從中性筆出墨量的研究來分析，主要集中在兩個方面：一是從筆頭墨水本身固有的性能出發來研究影響出墨量的規律；二是從出墨量測試方法的角度來研究影響出墨量的規律。從已發表的文獻來看，從筆頭、墨水本身固有的性能出發來研究的論文較多、較完整和全面，而從測試方法的角度出發來研究并不夠全面和充分。

為了系統地研究測試過程參數對中性筆出墨量測試的影響，我們對各因素進行了分析。為了使問題得到一定的簡化，我們假定書寫測試在恒溫恒濕的狀況下進行，測試環境條件：（23±2）℃，（50±5）%RH，這樣我們可以忽略測試環境對測試結果的影響。另外，書寫紙采用滿足ISO 14145要求的書寫紙。

1 中性筆出墨量與測試條件參數的關系分析

1.1 書寫間距與出墨量的關系分析

出墨量測試中，書寫間距指的是書寫測試時兩條螺旋線之間的距離。在文獻[6]中規定書寫間距一般在1～5 mm之間選擇。從我們平時測試中積累的經驗來看，在其他條件不變的情況下，書寫間距越大，測試的出墨量越大。關于出墨量與書寫間距的關系，在文獻[4]中有比較好的闡述。

1.2 書寫紙與出墨量的關系分析

在中性筆出墨量的測試中，標準中一般推薦使用70 g/m2的書寫紙。目前國際上通常使用滿足ISO 14145標準的紙張，我們國內書寫紙的標準為GB/T 12654，這一標準的實際指標同ISO 14145標準紙張的指標并不能很好地對應。出墨量測試中書寫紙對出墨量的影響也較大，通常國產書寫紙測試的出墨量較ISO 14145的紙張出墨量要大一些。部分具體的數據分析可以參閱文獻[5]。

1.3 書寫角度和書寫載荷與出墨量的關系分析

圖1是劃圓書寫機上筆芯和夾具的示意圖。從圖中來看，書寫紙和筆頭成θ時，實際的書寫壓力F，為書寫載荷G（夾具和專用砝碼的總重量）的一個軸向分力，大小為G。在文獻[6]中，書寫測試時的書寫角度一般在60°～70°之間選擇，推薦使用70°。隨著θ角的變大，F力也逐漸增大。所以，我們在研究書寫角度和書寫載荷與出墨量的關系時要同時結合研究。

圖1 書寫機上筆芯和夾具的示意圖

以上是我們對書寫測試中的主要參數與出墨量測試中的一些定性分析，下面我們從書寫角度和書寫載荷兩個因素開始進行基礎建模。至于其他的如書寫速度、書寫紙和書寫間距等對出墨量的影響，從以往測試的經驗來看，我們可以通過試驗，求出一個相應的系數，逐一加入到模型中。

2 中性筆出墨量與書寫角度及書寫載荷之間的數學模型建立

2.1 數學建模的意義

一般來說，數學建?？梢悦枋鰹椋簩τ诂F實世界的一個特定對象，為了一個特定目的，根據特有的內在規律，作出一些必要的簡化假設，運用適當的數學工具得到的一個數學結構[7]。在社會實際中的問題復雜多變，量與量之間的關系并不明顯，并不是套用某個數學公式或只w用某個學科、某個領域的知識就可以圓滿解決的。數學建模就是要從眾多的事物和現象中找到共同的、本質的東西，抓住問題的主要矛盾，從大量數據和定量分析中尋找并發現規律。然而，數學建模本身就是一個創造性的思維過程，需要綜合運用所掌握的知識和方法，創造性地分析解決實際中的問題。數學建模有許多方法，諸如線性規劃法、回歸分析法、層次分析法、微分方程法、量綱分析法和集合分析法等，這里我們采用量綱分析法。

2.2 量綱分析方法

所有的物理量都是由自身的物理屬性，以及由此而規定的量度單位這兩個要素構成的[8]。所謂量綱就是撇開物理量量度單位的具體名稱，而著重于該物理量的類別（屬性）[8]。顯然，量綱是物理量的實質[8]。我們采用dimA表示物理量A的量綱，通常以L代表長度量綱，M代表質量量綱，T代表時間量綱，Θ代表溫度量綱。這幾個為基本量綱，我們這里只需選取M、L、T為3個基本量綱，其他物理量的量綱為導出量綱。有些量是無量綱量，如圓周率π為周長與直徑的比值，角度α為弧長與曲率半徑的比值，這類同類量的比值都是無量綱量。任一物理量A的量綱都可用3個基本量綱的冪次形式來表示：

當用數學公式表示一個物理公式時，等號兩端必須保持量綱的一致性。這種性質稱為量綱齊次性原理。當方程中的各項具有相同的量綱時，這個方程被稱為是量綱齊次的。也只有具有相同量綱的量才能作比較或相加、減，這是量綱分析的基礎。量綱分析法是將物理量之間復雜的函數關系式，轉換成一些無量綱量之間關系式的一種方法[8]。在量綱齊次原理上發展起來的量綱分析法有兩種：一種稱為瑞利（Rayleigh）法；一種稱為布金漢（Buckingham）定理（又稱π定理）[8]。我們這里選用瑞利法。

瑞利法的基本原理是，倘若描述某一個物理現象同n個物理量有關[8]：

而其中的某一個物理量Bi可表示為其他物理量的指數乘積[8]：

寫成量綱形式[8]：

將式(4)中各項物理量的量綱按式(1)表示為基本量綱的冪次乘積形式，并根據量綱齊次原理，確定a，b，…，p,就可得出表達該物理現象的方程式。

2.3 出墨量與書寫角度及書寫載荷之間的數學模型建立

從中性筆出墨量與測試條件參數之間的關系分析和經驗判斷來看，書寫間距、書寫紙、書寫速度對出墨量的影響因素單一，但由于書寫角度和書寫載荷對出墨量測試的影響關系較為復雜，所以我們先固定上述其他因素，來探討書寫角度和書寫載荷對出墨量的影響關系，其他因素可以在此基礎上乘上對應的系數來解決。

2.3.1 條件假設和模型建立

在出墨量的測試中，我們一般以書寫長度100 m的出墨量為基準，用書寫前的質量同書寫100 m后的質量差值作為出墨量的數值。這里用m表示百米出墨量，其單位為mg。從圖1的分析可以得到，在上述其他條件固定的情況下，經分析出墨量m同書寫角度θ、書寫壓力F(為書寫載荷G的軸向分力，大小為G）、重力加速度g等相關：

我們按照瑞利法構造一個公式：

式中：

m—百米出墨量，單位mg；

m0—書寫角度為θ0（書寫參考角度，可以為60°～70°間的某一角度）時的百米出墨量，單位mg；

k1(θ)—書寫角度為θ時的角度系數，無單位，當θ為θ0時，k1(θ0)=0；

F—書寫壓力，單位N；

g—重力加速度，單位m/s2。

根據圖1中的分析，式（6）進一步進行演變，得到：

式中：

G—書寫載荷（夾具的總重量），單位N；

m1—夾具和專用砝碼的總質量，單位g；

k2(θ) —書寫角度為θ時的角度系數，等于k1(θ)(sinθ)a，無單位。

下面以式（7）為研究對象，來求解式中的指數a和b。式（7）左邊出墨量m的量綱dimm=M

將式（7）寫成量綱形式得到：

用基本量綱的（M，L，T）表示各物理量量綱得到：

根據量綱齊次原理，得：

M：a=1

L：a+b=0

T：-2（a+b）=0

解得a=1，b=-1

所以式（7）整理后，得到：

為了進一步簡化模型，根據量綱一致原則可以假設：

m0=k3m1

則可以得到：

2.3.2 模型初步分析

從得到的模型式(11)可以看出，夾具和專用砝碼的總質量m1和出墨量m之間有一定的比例關系。也就是夾具和專用砝碼的總質量m1固定的情況下，出墨量m會隨著書寫角度的變化而變化，這樣的變化可以用m和m1之間的比例系數的變化來表示。再深入分析，夾具和專用砝碼的總質量m1取一定的值，假定為m10，通過試驗，在不同的書寫角度θ，可以得到對應的出墨量m，這樣我們可以得到一組k(θ)的值，記為k(θ)10。如果m1取值改變為m11，通過試驗可以得到另一組k(θ)的值，記為k(θ)11。這樣成功地找到了一個無量綱量k(θ)，我們可以在研究k(θ)的變化規律中得到出墨量m的變化規律，使問題研究更加簡單，更加有針對性。

2.4 考慮書寫間距的模型

從書寫間距對出墨量的影響分析中可以得出，其影響同其他因素的關聯性不大，為此我們可以在式(11)基礎模型的基礎上乘上一個書寫間距系數Cdx來修正模型。由于基礎模型是在書寫間距為1 mm的狀況下建立的，其書寫間距系數可以記為Cd1，其值為1。我們可以利用實驗求出書寫間距為2 mm下出墨量（記為m2）同書寫間距1 mm時的出墨量（記為m1）進行比較，得到Cd2=m2/m1。同理可以得到Cd3、Cd4和Cd5。為此，考慮書寫間距的模型可以得到：

2.5 考慮書寫紙的模型

從書寫紙對出墨量的影響分析中可以得出，其影響同其他因素的關聯性不大，為此我們可以在式(12)的基礎上再乘上一個書寫紙系數CpX來修正模型。由于前面的模型都是在符合ISO 14145的書寫紙上建立的，其書寫紙系數可以記為Cp0，其值為1。我們可以利用實驗求出其他書寫紙，書寫間距1 mm時出墨量（記為mp1）同ISO 14145的書寫紙上，書寫間距1mm時的出墨量（記為m1）進行比較，得到Cp1=mp1/m1。為了求證書寫紙系數Cp對書寫間距的敏感性，也可以進一步試驗兩種書寫紙下其他書寫間距下的書寫紙系數CpX。為此，考慮書寫紙影響因素的模型可以寫為：

2.6 考慮書寫速度的模型

從書寫速度對出墨量的影響分析中可以得出，其影響同其他因素的關聯性不大，為此我們可以在式(13)的基礎上再乘上一個書寫速度系數CvX來修正模型。由于前面的模型書寫速度都是4.5 m/min，其書寫速度系數可以記為Cv0，其值為1。我們可以利用實驗求出其他書寫速度，書寫間距1 mm時的出墨量（記為mv1）同ISO 14145的書寫紙上，書寫速度是4.5 m/min、書寫間距1 mm時的出墨量（記為m1）進行比較，得到Cv1=mv1/m1。為此，考慮書寫紙影響因素的模型可以寫為：

3 實驗分析和模型參數求解

3.1 實驗方案設計

實驗以670/0.5筆頭、130 mm長度的某款筆芯進行試驗，此款筆芯按國家標準書寫長度要求大于400 m。隨著書寫的進行，由于筆頭磨損、墨水高度下降等原因，后200 m的出墨量會有明顯下降，所以這里考慮到筆芯百米出墨量的這一變化規律，實驗時只考察筆芯前200 m的出墨量，以此200 m出墨量的平均值作為考察對象，設置不同的書寫條件參數，考察筆芯的百米出墨量。

3.2 不同載荷與書寫角度的試驗

表1 不同載荷（夾具質量）與書寫角度時的出墨量和k(θ)

從表1的數據以平均值來比較，不同載荷下出墨量的數值會發生改變，呈現了中間高、兩頭低的現象；在載荷一定的情形下，隨著書寫角度的減小，出墨量有逐漸增大的趨勢；但是k(θ)1X卻從大到小，跟夾具質量形成負相關的線性關系，具體見圖2。

圖2 不同夾具質量時的出墨量k(θ)出墨量k (θ)1X

3.3 不同書寫間距的試驗

為進行不同書寫間距的出墨量比較，在載荷為0.993 N（夾具質量101.34 g）、書寫角度為65°、書寫速度為4.5 m/min的條件下進行書寫間距1 mm、2 mm、3 mm、4 mm、5 mm分別進行同款筆芯前200 m出墨量的測試。具體數據見表2。

表2 不同書寫間距時的出墨量和Cdx

從表2可以看出,隨著書寫間距的增加，出墨量有較大幅度的增大。同1 mm的書寫間距比較，2 mm時出墨量增加較多;隨著書寫間距的進一步增加，出墨量增加幅度減小,增加的幅度在9%～12%之間。

3.4 不同書寫紙的試驗

為進行不同書寫紙間的出墨量比較，在載荷為0.993 N（夾具質量101.34 g）、書寫角度為65°、書寫速度為4.5 m/min、書寫間距1 mm的條件下進行不同書寫紙同款筆芯前200 m出墨量的測試。具體數據見表3。從表3可以得出，兩種書寫紙間的出墨量差異較小，Cpx的平均值也可以取為1。

表3 不同書寫紙間出墨量的比較

表3 不同書寫紙間出墨量的比較

3.5 不同書寫速度的試驗

為進行不同書寫速度的出墨量比較，在載荷為0.993 N（夾具質量101.34 g）、書寫角度為65°、書寫間距1 mm的條件下進行不同書寫速度下的出墨量比較試驗，書寫速度為3.5 m/min、4.5 m/min、5.5 m/min、6.5 m/min、7.5 m/min分別進行同款筆芯前200 m出墨量的測試。具體數據見表4，書寫速度與Cvx的曲線圖見圖3。從圖3 看出，書寫速度與Cvx之間呈線性負相關。經回歸計算，相關系數R的絕對值為0.92，書寫速度與Cvx線性相關性良好。一次線性回歸方程為：

圖3 不同書寫速度下的Cvx

4 模型應用和討論

4.1 模型各系數的確定

從試驗的數據分析和數學模型的各系數求解過程來看，書寫角度和書寫載荷系數k(θ)可以從表1中查詢，或者表1中未能直接查得的，可以結合夾具質量m1通過線性插值求得。書寫間距系數Cdx可以通過表2直接查得。從分析和實際的數據來看這里的書寫紙系數Cpx取為1。書寫速度系數Cvx可以通過式（15）進行回歸計算獲得。模型計算總表見表5。

表5 模型計算總表

4.2 模型的應用

對于一款中性墨水圓珠筆筆芯，在研發階段要經過多方案、長時間的試驗才能成功。這期間筆頭的選型、墨水的確定、筆芯主要尺寸的確定等，都需要足夠的試驗。通過模型的應用，可以很好地幫助筆頭和墨水的選型，可以利用模型進行不同書寫條件下的結果比較、分析和預測，可以協調筆頭、墨水供應商與筆芯成品生產商之間檢測標準的不一致帶來的矛盾。同時模型也可以指導筆芯研發人員進行快速試驗，確定主要參數（如墨水量、出墨量）時提供指導和方便。但是模型應用中應該注意的是，文章中是以0.5 mm子彈頭和特定墨水匹配后作為研究對象，其他的筆頭和墨水的匹配模型中的系數應進行具體的試驗，本文各系數請參考使用。各種不同型號的筆頭因為結構參數的不同，與不同的墨水匹配后書寫角度和書寫載荷系數k(θ)需要重新確定，其他系數也應謹慎使用，要根據實驗室實際情況通過試驗驗證確認。

5 結論和展望

5.1 結論

對中性筆芯出墨量測試參數的分析，運用量綱分析的方法成功構建了中性筆出墨量與測試條件參數之間的數學模型，并以0.5 mm子彈頭中性筆芯為例進行了多組實際測試，求解了模型中的各個參數。從試驗和模型中可以獲得以下有用的結論：

（1）在載荷一定，也即夾具質量一定的情況下，出墨量隨著書寫角度的變化而變化，角度越小出墨量越大，這一趨勢較為明顯。

（3）隨著書寫間距的增加，出墨量增加明顯，而且書寫間距2 mm時比1 mm時出墨量增加的幅度比3 mm對比2 mm時的幅度更大，呈現了前大后小的現象。

（4）隨著書寫速度的增加，出墨量減少，兩者呈現了線性負相關的關系。

5.2 展望

通過中性筆出墨量與測試條件參數之間的數學模型構建和0.5 mm子彈頭筆芯的測試求解和驗證，基本掌握了中性筆出墨量與測試條件參數之間的變化規律，在實際工作中可以運用模型針對實際問題進行必要的預測和求解。然而，這一模型中參數的求解是基于0.5 mm子彈頭筆芯的，對于中性筆芯筆頭型號眾多，筆頭內部結構有較大區別，而且筆頭與墨水之間需要進行匹配，其匹配方式較多，因此模型使用中需要針對不同的匹配組合進行參數求解。同樣也希望針對不同的匹配進行數學模型的求解，積累數據建立模型數據庫，這是對中性筆書寫原理的進一步研究具有重要意義的一項工作。