摘要
裝着三公斤物品的紙袋，提起後最初似乎能保持原狀；走一段路，手挽卻緩慢變長，令人懷疑物件是否在偷偷增加重量。本文把此日常錯覺化為一項有明確物理邊界的模型研究：若外力恆定，具有黏彈性的纖維網絡是否仍可累積變形？以四段有效受拉帶、負載角度與橫截面描述力學幾何；比較 Kelvin–Voigt、標準線性固體與引入永久流動的 Burgers 類模型，利用階躍負載的蠕變柔量及卸載疊加原理，推導提袋期間的伸長、放下即時收縮、延遲恢復及假設性永久殘餘。設定三公斤、帶寬十二毫米、厚度零點三毫米、有效長度八十毫米、與垂直方向三十度、瞬時模量六百 MPa、延遲模量二百五十 MPa、遲滯時間三百秒、持續提拿十分鐘，標準線性固體預測卸載前伸長 0.96763 毫米，其中即時彈性 0.31466 毫米、延遲部分 0.65297 毫米；卸載十五分鐘後延遲殘餘約 0.03251 毫米，理論上最終完全復原。另一個加入不可恢復流動的示範模型則預測 0.22655 毫米的永久分量。四萬次合成參數試驗產生 0.218375 的「伸長超過自訂 1.2 毫米界線」機率，但此概率只對指定先驗分布成立，絕不能解作市售紙袋的不合格率。五張圖揭示提袋時間、承重、帶寬、濕度代理參數與異質性假設對視覺感受的影響。本文指出，恆定重量下的延遲伸長可以是力學上正常的黏彈響應；但若沒有卸載恢復曲線、真實幾何、含水率及損傷檢查，就不能把「長了」直接判作「快斷了」，亦不能把可恢復蠕變等同於永久塑性。
Abstract
Why might the handle of a paper shopping bag grow longer under a constant load? We formalize the question using a geometry-aware synthetic model of viscoelastic creep. A four-branch tension approximation converts the carried mass, strip width and thickness, and strap inclination into an effective stress. The Kelvin–Voigt, standard linear solid (SLS), and exploratory Burgers-like models are compared through creep compliance and unloading superposition. Under deliberately assumed parameters for a 3 kg bag, the SLS predicts 0.96763 mm of extension immediately before unloading at 600 s. This comprises 0.31466 mm of instantaneous elastic response and 0.65297 mm of delayed deformation; the latter decays to 0.03251 mm after 15 min of recovery and vanishes asymptotically. Adding an illustrative irreversible viscous branch yields a permanent component of 0.22655 mm. A seeded simulation of 40,000 synthetic variations produces a 0.218375 probability of exceeding an arbitrary 1.2 mm visibility threshold, which is conditional on assumed parameter distributions and is not an estimate of failure in commercial bags. Five research figures evaluate competing constitutive descriptions, holding time, geometry, humidity proxies and uncertainty. The study distinguishes recoverable creep, assumed irreversible strain and unmodelled fracture, and specifies measurements required for empirical validation.
關鍵詞： 紙袋；手挽；蠕變；黏彈性；Kelvin–Voigt；標準線性固體；Burgers 模型；纖維網絡；濕度；參數不確定性；合成蒙地卡羅
1. 引言：重量沒有變，長度卻是時間的函數
假設某人把三公斤購物品放在紙袋裏，從商店走回家，沿途並沒有加入物品。從牛頓力學看，袋內重物對手挽的垂直合力約固定於二十九點四三牛頓；若手挽可以視為純線彈性的桿，且其截面、方向及溫度均不變，那麼應變會在受力瞬間建立，之後保持近乎定值。然而，許多由紙、纖維或聚合物構成的物件，即使沒有更多負重，其變形亦可能繼續隨時間增加。這種在持續應力下的時間依賴應變稱為蠕變（creep），是黏彈性研究的基本現象，並不違反力與變形的關係；相反，它提示讀者不能把材料的本構關係縮減成一個不含時間的彈簧常數（Ward & Sweeney, 2012）。
日常觀察有三個互不等價的解釋。首先，手挽可能只是發生可恢復的延遲拉長：放下紙袋後，長度慢慢回復，說明材料有記憶效應，但未必永久損壞。其次，接縫或材料可能發生滑移、局部塑性或纖維間不可逆重排，使卸載很久後仍留有一段殘餘長度。第三，手挽的兩端、黏合位置或孔洞可能出現局部損傷，甚至接近破斷。肉眼看到的「變長」無法單獨區分三者。研究若把伸長量直接解讀成破斷概率，便會把位移指標、疲勞／損傷指標及材料承載能力錯誤地混合。
本文的問題故意有些荒誕：紙袋究竟會不會在路上自己「增磅」？嚴肅答案涉及材料工程中最基礎的區分：外力可以固定，但材料內部的應力傳遞與微觀重排仍可能隨時間演變。紙纖維形成帶方向性、含有孔隙及纖維間接觸的網絡，受外加拉伸與水分變動時會出現複雜非線性時間響應；真正的材料紙張絕非本文的一維均質固體，這點可從 Samantray 等人（2021）的多尺度紙纖維模型得到方法論上的支持。本文因此選擇可檢查的最小模型，不把簡化當成工業認證。
核心目的有四：第一，證明在恆定載重下，時間依賴伸長與機械平衡並不矛盾；第二，區分可恢復蠕變與人們容易誤認為同一回事的永久變形；第三，評估負載時間及受力幾何如何改變同一個重量造成的位移；第四，用透明的合成參數不確定性展示單一示範曲線如何不足以提供產品族群的風險頻率。本文不是預測某個品牌手挽何時破裂，而是針對一項荒誕日常觀察建立最小充分的力學解釋與可重現計算流程。
這個研究也具有一般化價值。塑膠袋提帶、書包肩帶、紙製包裝邊緣、布料拉帶與黏着連接都可能在一定程度上存在時間效應；可是材料種類、濕度、載重模式及製作工藝各異，不能把紙袋假設參數移植到其他物件。為阻止過度外推，文中每項數值都標示為研究者選取的示例，每張圖均標記為合成計算，並把「可見伸長」與「危險載重」分成兩套概念。讀者所學應是建立材料模型的方法，而不是一張似是而非的市售紙袋安全表。
2. 文獻回顧：由抽象彈簧到吸濕的纖維網絡
Ward 與 Sweeney（2012）在聚合物力學著作討論線性黏彈性、Boltzmann 疊加原理、蠕變柔量及應力鬆弛，並介紹 Kelvin、Voigt、Maxwell 和標準線性固體等元件模型。這些理論為「固定力、長度仍變」提供數學語言。其框架討論的是材料類別的本構理論，並沒有量度本文假設的紙袋。因此我們引用的是模型形式與疊加方法，而不會引用其作為三公斤手挽應變的經驗來源。最小彈簧與黏壺模型亦會錯過實際紙纖維的多時間尺度和不可線性化行為，這些侷限應在研究起點便列明。
Lin（2024）的標準線性固體研究仔細比較 Maxwell 與 Kelvin 形式的三元件結構，說明兩者能以適當參數轉換描述等價的線性響應，而各元素的物理詮釋必須保持清晰。本文選擇「瞬時彈簧串聯 Kelvin–Voigt 元件」的 Kelvin 表示，因為其階躍蠕變柔量可以直接拆成即時伸長與延遲伸長，方便討論人眼何時發現紙袋下垂。這一選擇並不表示真實紙袋的微觀纖維可一一對應成兩條理想彈簧和一個黏壺。該文主要解釋本構模型等價性，並沒有為本文所定模量提供材料量測。
紙張與紙板文獻使抽象模型與研究題目產生真正關聯。Samantray、Massart、Peerlings 與 Geers（2021）在《Mechanics of Materials》以纖維層級速率依賴模型及網絡尺度分析研究紙的力學響應，指出外載與水分變動會令宏觀應變具時間依賴性。這項研究足以支持以蠕變為出發點的學術合理性；它沒有測試本文的紙袋手挽結構，亦不表示任何固定厚度的紙材皆可以使用同一組模量。從纖維尺度到實物手挽涉及構造、加工、黏合、折線與孔洞的多層轉譯，這正是本文無法省略的外部效度問題。
Olsson、Salmén、Eder 與 Burgert（2007）研究木纖維的機械吸濕蠕變，文獻摘要表明受拉纖維在循環相對濕度環境中的蠕變，可以大於最高恆定濕度下的蠕變；Dong、Olsson 與 Salmén（2010）則檢視纖維形態與該現象的關係。兩項研究提醒我們：就算購物重量完全不變，走過空調商場與室外潮濕空氣時的材料歷程可能不同。但這些結果並不能在沒有試驗的情況下變成一條通用的濕度修正係數；本文第四幅圖因此明確採用任意但公開的「濕度代理參數」展示敏感度，並不聲稱其曲線曾由上述論文擬合。
Gunderson 與 Tobey（1990）直接討論紙板在循環相對濕度下的拉伸蠕變及濕度變化速率，為「紙類材料 + 張力 + 水分」的研究交集提供較接近的資料來源。其摘要提及特定條件下蠕變對濕度改變速率的依賴結論，不能被簡化成「濕度總是不重要」。事實上本文正是缺乏特定產品的含水率、厚度、織向及膠合條件，才不會以固定相對濕度一個數字冒充完整紙袋材料的狀態。對紙袋手挽而言，纖維方向性與黏合點兩者都可能令一維平均應力估計過分樂觀。
Olsson 與 Salmén（2014）研究紙漿纖維和紙的機械吸濕蠕變，從標題、書目及公開摘要可確認研究主題涉及水分與應變；本文不宣稱已檢查其全部原始實驗，亦不引用未能查明的精確測試係數。各項研究合併可支持一個謹慎命題：紙纖維網絡可能對時間、外力和濕度歷程敏感；但它們不能直接決定某紙袋的承重上限、可見伸長界線或損壞機率。剩下的具體數值全是本文清楚定義的模型產物。
由此可見一個值得研究的空位。既有力學文獻可以告訴我們如何測定應變及怎樣建立本構方程，但普通使用者只得到一個不完整觀察：手挽彷彿愈提愈長。從此觀察到科學解釋，必須同時處理「重物的力如何分配到細帶」、「一維本構對長時間施載會產生甚麼函數形狀」、「卸載後會否復原」及「圖形上的跨情境差異如何影響所見」。本文不是宣稱發現一種未報道過的纖維機制，而是把這些可檢查的層次合成一份完整的荒誕研究。
3. 研究問題、操作定義及可被推翻的假設
RQ1 問：在重量、提袋角度、截面及外界條件固定時，手挽應變是否仍能增加？H1 預測，只要延遲柔量為正且尚未完全飽和，標準線性固體在持續負載期間的伸長率便為正，無須增加重量。這是一個模型內的理論可檢驗命題。未來若對真實手挽施加恒定張力並取得反向或零伸長時間序列，便需要重新判斷其本構參數或現有模型的適用性。
RQ2 問：放下紙袋後，剛才增加的長度是否會永久留下？H2 將 Kelvin–Voigt 與標準線性固體的延遲恢復和具永久黏流支路的模型分開：純線性固體應在足夠長的卸載後使所有應變趨近於零，含恆定黏性流動的模型則可能留下永久分量。此假設並不是說「真實紙袋必定復原」；它提出的是可用卸載觀察辨別模型的實驗設計。
RQ3 問：相同三公斤在不同帶寬、厚度與握提角度下會否產生同樣伸長？H3 預測，在均質軸向受力近似中，截面愈小或受力帶與垂直線夾角愈大，每個有效受拉分支的張力愈高，伸長亦愈大。然而，此預測要求有效分支數不變、各支受力均勻，亦不涵蓋局部應力集中。
RQ4 問：若外界材料參數有變動，某條固定「變長超過一點二毫米」界線的機率能否被計算？H4 指出，在清楚指定分布時，可以計算條件概率；改變分布便會改變概率。這一題的研究價值不在百分數本身，而在防止把任意合成先驗造成的頻率誤稱為消費者實際面對的失效率。紙袋破斷往往由缺陷、接縫、濕潤與局部應力共同控制，本文的「超過伸長界線」完全沒有被校準為破斷標準。
本文將「手挽長度」操作性定義為沿有效受拉路徑、參考卸載長度為八十毫米的一段帶材的伸長量；這只是觀察窗口，並非整個手提環的外輪廓。本文沒有模擬人的手指位置改變、紙帶環形曲率拉直、袋口翹曲或袋體重心晃動。若日常眼見的「變長」主要來自幾何拉直而非材料應變，本模型的機制權重就需要另外重估。把觀察定義寫清楚，是讓數學結果可被否證的前提。
4. 受力幾何：一袋重物對四段紙帶並不客氣
先從最簡單的力學平衡建立模型。令袋內質量為 m，重力加速度 g0 為 9.81 m/s²。假設兩條手挽各可分解成兩段有效直線受拉部分，共有 n=4 段，而且各段與垂直方向夾角為 θ。四段受力對稱，承受相同張力 F_t；以垂直方向合力守恆可得四段張力的垂直分量之和恰好平衡重量。這是一個簡化的四支幾何模型，不代表所有紙袋均設有相同的帶形與對稱性。
若每段有效承力截面的寬度 b、厚度 h，且應力在截面上均勻分布，工程應力 σ 等於 F_t 除以面積 bh。使用牛頓及平方毫米時，σ 的單位為 N/mm²，即 MPa，因此公式中的 g0 可保持以 N/kg 表達，而 b 與 h 用毫米。若改成 SI 米及帕斯卡而又忘記一百萬的換算，會造成嚴重單位錯誤；本研究程式全部採用 MPa–mm–s 一致制。
取 m=3 kg、b=12 mm、h=0.30 mm、θ=30°，得 σ 約為 2.35992 MPa。這個數值是可檢查的合力及面積運算結果，不是測量得到的紙袋應力。假如只有三段有效分攤重量，或當握提導致兩條手挽長短不一，便須重新解決靜力平衡甚至接觸約束，而不是繼續使用分母中的四。尤其當某段手挽因製造缺陷而局部脫膠，等分受力假設可以瞬間失效；本模型的有效應力只是基線，不應被稱為最大局部應力。
夾角會直接改變張力。θ=0° 時，每段承受重量四分之一；θ 增大到六十度，cosθ 變為二分之一，每段張力即比垂直情境增加一倍。當 θ 接近九十度，簡化式會發散，這不是說紙袋可以產生無窮張力，而是表示固定四段、純軸向拉力、無重新配置的理想幾何失去物理可行性。因此所有合成情境把角度限制於十八至四十二度，既避免奇異點，也提醒讀者不要把公式外推到近乎水平的手挽。
帶寬和厚度以乘積方式進入分母，看似容易透過增加任一個尺寸降低應力。但實際紙張厚度可能包含塗層、皺褶與內部孔隙，名義截面不等於全部有效承力截面；此外，折線與孔洞附近存在應力集中。這些情況不能僅以一個平均截面積修正完畢，可能需要平面應力分析、局部材料實測及破壞力學方法。本文選擇初等軸向應力，目的僅是把「同一袋重量」轉換成材料模型的清楚輸入。
5. 三種本構機制：即時、延遲和不能忘記的變形
把材料純化成線彈性彈簧時，應力與應變由 σ=Eε 描述，E 是楊氏模量；只要 σ 不變，ε 亦不變。這適合說明剛提起紙袋瞬間的一項拉長，卻不可能單獨解釋恆載時長度繼續增加。黏彈性加入黏壺後，材料可用內部變數描述緩慢傳遞，恆定外力遂不代表即時平衡配置已達最後狀態。值得警惕的是，這些「彈簧」、「黏壺」是數學機械類比，不是在紙纖維網絡內找到三個一模一樣的實物零件。
Kelvin–Voigt 元件讓彈簧 E1 與黏壺 η 並聯；兩支承受相同應變，應力相加。在恆定應力 σ0、初始應變為零的載入階躍下，應變滿足一階微分方程。對其解施加零初始值，就得到漸近的指數增長。此模型沒有瞬時應變突跳，因為黏壺不能在有限應力下瞬間形變；因此它可以描述純延遲蠕變，卻難以涵蓋真實手挽最初一拉就下沉的快速彈性部分。
標準線性固體（standard linear solid；SLS）的 Kelvin 表示在上述並聯元件之外再串聯一支彈簧 E0。串聯意味兩支承受相同應力、總應變相加，所以它在加載瞬間擁有 σ0/E0 的即時變形，再加上一個逐步建立、最後飽和的延遲分量。在本文公式內，E0、E1 都以 MPa 表示，τ 以秒表示，η=E1τ 的單位為 MPa·s。只要 E0、E1 和 τ 均為正，蠕變曲線從正值開始，單調增加，但在固定應力下有有限極限。
其中 J(t) 是階躍蠕變柔量，單位 MPa⁻¹；t 為加載後經過的秒數。極短時間的柔量為 1/E0，極長時間的柔量為 1/E0+1/E1。其時間導數是 exp(−t/τ)/(E1τ)，永遠非負，且隨 t 增大而減少。這證明了 H1 在 SLS 的條件下成立：重量不用增加，手挽仍可繼續變長；但每秒新增的長度逐漸下降，並非永遠按固定速度拉長。若真實觀察顯示晚期反而加速伸長，則單一延遲時間的 SLS 可能嚴重不足。
不過，完全可恢復的 SLS 不會留下真正永久應變。為使 RQ2 的比較具有辨別力，本文再設一個 Burgers 類探索模型：在 SLS 柔量上加入恆定黏性流動項 t/η_p，η_p 是另一個正的黏度參數。這使加載時產生隨時間線性增加的額外應變；負載移除後，該項不會恢復。稱它為「類」模型，是因為我們直接透過串聯柔量建構此響應，不把它誤解成任何已被證實的紙袋塑性機理。巨觀永久殘餘也可能來自接縫滑移、纖維損傷或局部變形，本構中的流動黏壺不能識別其微觀起因。
為甚麼不把實物紙袋直接假設成 Burgers？因為一條長時間的上升曲線往往可以由多個延遲時間組成的可恢復模型近似，也可以由不可恢復流動近似；單看受載階段，兩者可能不易區分。只有在充分控制應力、環境與幾何之後觀察卸載恢復，才能對永久分量是否存在提出更有力的判斷。即使放下一小時仍有伸長，也可能是極慢的恢復時間，而非真正無限時間的殘餘。這個識別限制是本研究選擇三個模型而不是只畫一條好看曲線的原因。
6. 卸載疊加：放下紙袋後的那一下與之後的十五分鐘
研究若只計算 t=600 秒當下有多長，便無法把黏彈性與損傷區分。對具線性時間平移不變性的本構系統，在 t=0 施加 σ0，並於 t=T 取消負載，可把歷程寫成一個正階躍與一個在 T 時刻開始的負階躍。Boltzmann 疊加下，任一時刻的應變可以用兩份柔量之差表達。這個推導要求應變小、材料未因負載而改變參數、環境恒定且沒有非線性損傷；實際紙袋若已明顯撕裂，這些條件顯然失效。
式中 H 是階躍函數，在 t<T 等於零，在 t≥T 等於一。對 SLS，卸載瞬間的純彈性項會立即消失，因此剛卸載的殘餘應變只等於延遲元件已建立的部分。卸載後 s 秒，剩餘延遲應變以 exp(−s/τ) 衰減；它永遠非負且趨向零。用 L0 表示有效參考長度，就可以把無量綱應變乘成可觀測位移 δ，單位毫米。
因此加載時間愈長，剛卸載後仍未回復的長度愈大；但一旦卸載，在這個理想模型內不再持續向下伸長。本文不把卸載後殘餘位移稱為「永久變形」，除非模型具有不能隨 s 增加消失的項。對 Burgers 類流動分支，卸載後永久應變為 σ0T/η_p，因為加載 T 秒期間流動所累積的應變在撤除應力後固定留下。其餘延遲分量仍會照 SLS 形式衰減。
這些式子同時提供模型驗證的方法。在 t=T 前一瞬，SLS 總伸長是 L0σ0/E0 加上延遲項；在卸載後一瞬，則只剩延遲項；兩者的差必須恰好等於 L0σ0/E0。假如程式計算所得不是這個數字，便代表時間階躍實作或單位換算出錯。進一步地，在遠大於 τ 的恢復時間，SLS 的殘餘必須極接近零，Burgers 類模型則趨向正的永久分量。本文的程式把兩種檢查分別寫入數值稽核輸出。
需強調，「瞬時」在本研究是數學階躍的極限；人手提起或放下袋子需要有限時間，提起與放下亦伴隨慣性、摩擦、接觸重新排列及手挽曲率變化。若以高速相機或位移計記錄真實紙袋，應以有限升載時間代替不連續階躍，並量度握持造成的額外幾何位移。否則就算觀測曲線符合指數，也不能唯一證明紙纖維的遲滯時間等於本文的 τ。
7. 計算設計、單位守恆及可重現性
本文定義一個可重跑的基準情境。袋內質量三公斤，四段等分受力的手挽，單段帶寬十二毫米、厚度零點三毫米，與垂直方向成三十度，單段有效參考長度八十毫米。材料採用 E0=600 MPa、E1=250 MPa、τ=300 秒的示意 SLS，因此並聯黏壺等效黏度 η=E1τ=75,000 MPa·s。提拿持續時間 T=600 秒，隨後卸載一千二百秒以觀察恢復。另為 Burgers 類探索支路選擇 η_p=500,000 MPa·s。所有數字只作量級清楚、算術容易核對的展示，未從文獻反推或由實物校準。
| 參數 | 基準設定 | 單位與詮釋 |
|---|---|---|
| m | 3.0 | kg；袋內固定質量 |
| b、h | 12.0、0.30 | mm；每段帶寬、有效厚度 |
| θ、L0 | 30、80 | 度、mm；受拉角與參考長度 |
| E0、E1 | 600、250 | MPa；假設性瞬時與延遲模量 |
| τ、T | 300、600 | s；單一遲滯時間、提拿時間 |
| η_p | 500,000 | MPa·s；探索性不可恢復流動項 |
| δ* | 1.20 | mm；研究者自定觀測界線，非安全標準 |
數值運算採用 Python、NumPy 與 Matplotlib。analysis.py 從頭產生五幅圖，並把摘要輸出至 data/summary.json；data/sampled_cases.csv 保留合成試驗前一百五十筆輸入與輸出供肉眼檢查。所有圖皆以英文座標及 SI 相容單位標示，避免中文字形在繪圖環境缺字。繪圖不需任何外部私人數據。隨機試驗固定種子 2026102701，樣本數四萬；單次執行會覆寫所生成的圖和資料檔，以便核對數值一致性。原始論文、程式、參考文獻核對紀錄與模擬參數都保存在同一論文資料夾。
為了測試不確定性而又不冒充實測資料，本文自行指定七項互相獨立的合成變動。質量、帶寬、厚度、E0、E1 與 τ 分別服從以各自基準值為算術平均數的對數常態分布；其變異係數依次為 12%、10%、12%、15%、15% 與 20%。角度則獨立服從十八至四十二度的均勻分布。對數常態的優點是取值自然為正，避免在普通常態抽樣下產生負厚度與負模量；缺點是高值尾部也屬研究者假定，沒有材料統計證據。角度均勻分布同樣是演算方便的選擇，絕非人類提袋姿勢的調查結果。
如某參數 X 的指定平均數為 μ，變異係數為 c，程式採用 log(X) 服從常態分布，其對數變異為 ln(1+c²)，對數平均為 ln(μ) 減去該變異的一半。如此生成的 X 才會以 μ 為算術平均數。若直接用 μ 當對數平均數，所生成參數的平均值便會偏高，導致伸長分布在不知不覺間偏移。這是合成研究亦需要嚴守的統計定義，否則重現同一個隨機種子仍不能保證模型意義相同。
每次合成試驗先計算該次的 σ，再由該次的 E0、E1、τ 和 T 求出載重末期伸長；最後計算高於自訂界線 δ*=1.2 毫米的指示值。研究不引入無資料根據的「紙袋不合格」分類，不計算破斷機率，也不使用不存在的真人樣本作統計顯著性檢定。蒙地卡羅估計的二項標準誤只衡量這個指定生成機制下有限次抽樣的隨機誤差；即使其數值很小，仍不能掩蓋材料分布、形狀及本構模型的結構性不確定性。
8. 研究結果一：固定負重下為甚麼會愈提愈長？
根據第四節求得的應力約 2.359919 MPa，基準 SLS 的零時刻即時應變是 σ/E0=0.003933，乘上八十毫米參考長度，得到 0.314656 毫米的即時伸長。經過十分鐘，延遲柔量已建立 1−exp(−600/300)=1−exp(−2) 約 86.47% 的長期可恢復份額，額外增加約 0.652972 毫米。兩項合計，負載撤去之前的預測伸長為 0.967628 毫米，精確程式輸出見附錄資料。這三項數字並不需要說「重物偷偷變重」；只需要承認材料對相同 σ 的應變會隨時間緩慢追上另一個漸近狀態。
圖一把 Kelvin–Voigt、SLS 與探索性 Burgers 類模型置於同一坐標。Kelvin–Voigt 從零開始，負載時間增長時逐漸伸長；SLS 則一開始已包含瞬時彈性分量。兩個完全可恢復模型在卸載後均朝零收斂，但卸載當刻 SLS 因移除串聯彈簧的張力而出現一個即時下降。Burgers 類曲線在負載期間多一項線性時間依賴伸長，卸載後還保留非零殘餘。三者同樣可在重量固定時出現「愈提愈長」，證明單一生活觀察不足以識別材料支路的種類。
圖 1： 三個模型於固定三公斤示例載荷下的載入與卸載位移。第十分鐘的虛線標示卸載時刻。圖形來源：本文 analysis.py 的純理論數值計算；不含實測紙袋資料。SLS 與 Kelvin–Voigt 將長期回復至零；額外黏流項只屬探索性永久殘餘機制。
從力學上解釋，固定應力並不等於固定應變速率，更不等於固定應變。SLS 在 t 時刻的伸長速率為 L0σ exp(−t/τ)/(E1τ)，逐步下降到零。於載入初期，應變速率較快，之後逐漸慢下來；如果觀察者只在開始與十分鐘後拍兩張照片，會忽略中間的曲率，錯把差值當作一段恆定速度的拉長。相反，永久黏流支路會使長時間斜率不歸零，此處正好提出實驗上應延長觀察窗的理由。只把一條短時曲線擬合得很漂亮，並不等於已識別長期本構。
這一結果支持 H1 的有限範圍版本：在所選 SLS 且等分受力的假設下，負載末期總位移大於剛開始的位移，且外力沒有變大。然而，「伸長必然安全」不是此式可推論的結果。不同的紙袋可能在可見一毫米變形前已從黏合處破壞；也有紙帶可能伸長數毫米後仍遠未達破斷極限。若要作產品安全建議，必須另有強度、缺陷分布、疲勞歷程及法規測試資料。本文的位移只是一個行為變數，不是可靠度函數。
9. 研究結果二：卸載之後，甚麼才叫真正的永久變長？
在 T=600 秒卸載的前一瞬，SLS 總位移為 0.967628 毫米；撤去外力後，瞬時彈簧釋放的位移為 0.314656 毫米，因而剛卸載時還留 0.652972 毫米。再經過十五分鐘，也就是九百秒或三倍 τ，延遲殘餘下降到 0.032510 毫米。假如延長恢復到無窮時間，這項殘餘必定漸近於零。數值差值亦核對：卸載瞬間下降的位移與獨立計得的 L0σ/E0 相符至浮點精度。這是模型自身的重要一致性檢查。
引入 η_p=500,000 MPa·s 的探索性額外黏流支路後，六百秒內累積的不可恢復位移為 L0σT/η_p=0.226552 毫米。注意，這並非「我們量到紙袋永久變長 0.23 毫米」。它只表示：如果一個特定材料可合理用這條額外線性流動支路描述，而且這項黏度參數真的是該值，那麼在此載重歷程之後會產生這個永久分量。若 η_p 趨於無窮大，該項消失，模型退回純 SLS；若長期外力很大，線性流動式又可能不符合真實材料的損傷或應變硬化機制。
圖二把承重時間 T 當作自變數，顯示卸載前總位移、剛卸載後殘餘以及純即時彈性部分。即時彈性在固定負重與幾何下不隨 T 改變；可恢復的延遲部分則在初期增加較快，逐漸趨近飽和值。圖上即使顯示提三十分鐘仍比五分鐘更長，也不能把增長當成永久損傷的證據。最直接的實驗反證是卸載後再測同一標記距離，尤其在多個相隔 τ 的時間點收集恢復軌跡。
圖 2： 持續負載時間改變總位移及卸載立即留下的可恢復延遲位移；純即時彈性分量在固定載重下保持不變。來源：假設性 SLS 的解析公式，並非長期觀察資料。圖示的「卸載後殘餘」不能稱為永久塑性，因為足夠長時間後會消失。
這項分辨有一個實務上的思想實驗。設有人從店舖出來提了一小時，放下紙袋後馬上發現手挽仍較長，便說手挽「已經壞了」。假如只等三秒，這個判斷忽略了材料的延遲恢復；假如只等三十分鐘，亦仍未能排除比三百秒更慢的支路。要證明不可恢復的部分，需要長時間恢復測試或多時間尺度辨識，並檢查物理損傷，不能以一個任意等待期直接定義無限時間的塑性。SLS 之所以有教學價值，是它給出清楚的零殘餘對照，讓人理解真正永久變形需要額外假設。
另一個重要含意是受載終點與卸載終點之間沒有同義關係。圖一中 SLS 在第十分鐘有約 0.97 毫米伸長，圖二中它卸載後卻只留下約 0.65 毫米；兩個數字對應不同時刻和邊界條件。任何把「負載伸長」誤抄為「永久伸長」的報告，都會高估本模型內的不可恢復變形。研究者即使沒有實際量測設備，至少應在符號與圖例中分清楚「負載中」、「剛卸載」與「恢復後」三種狀態。
10. 研究結果三：相同重量之下的幾何放大
固定材料參數與提袋時間，基準位移與應力成正比，而應力又與 m/(bh cosθ) 成正比。質量增加一倍，理想模型位移增加一倍；寬度加倍，位移減半；厚度增加一倍，位移亦減半；角度則依 cosθ 的倒數調節。這種比例律源自靜力平衡與線性本構，是確定的數學推論，並非統計迴歸關係。現實紙袋隨質量增加可能令角度、有效分支數、局部折線與截面同時變化，所以比例律不能在大載重下任意外推。
圖三固定 h=0.30 mm、θ=30°、T=600 秒，將帶寬分成八、十二與十六毫米三條示例曲線。橫向質量範圍由半公斤到六公斤，只為展示各假設式如何改變預測伸長，絕不表示六公斤位於任何紙袋的許可負載內。對十二毫米的基準帶寬，伸長達 1.20 毫米所需質量約為 3.72044 公斤；這個臨界質量完全由自訂位移界線決定，與破斷強度無關。若把界線改為 0.60 毫米，臨界質量在其他條件相同時便減半。此處的「臨界」不代表材料性能突然轉變。
圖 3： 假設不同有效帶寬時，載重質量與第十分鐘伸長的線性關係。水平虛線為自訂 1.20 mm 可見伸長界線；沒有破斷、認證或產品安全資料參與計算。來源：本文靜力平衡及 SLS 理論計算。
把 δ*=L0σJ(T) 與第四節的幾何式聯立，可反求某個位移門檻下的質量 m*。它跟 b、h、cosθ 成正比，跟 L0 及 J(T) 成反比。由於 J(T) 會隨材料軟化與負載時間改變，同一條線若沒有標明環境與時間，便不適合作工程圖。更重要的是，m* 只是位移門檻的反向運算結果，絕非材料強度極限。
若真實紙袋的兩條手挽長度不同，重量會不等分；若一條局部脫膠，另一條的實際張力也可能上升。手挽與袋口的接駁位置還會發生彎曲、剝離與剪切，不符合純軸向膜應力。這些局部現象可能主導實際失效，與可見平均伸長沒有一對一映射。因此圖三的價值在於揭示幾何如何影響位移，而不是替任何紙袋設定「最多可放幾公斤」。
11. 濕度代理：水分歷程不等於當下的空氣濕度
紙纖維受水分影響令本研究不只是理想彈簧的遊戲。Olsson 等人（2007）、Dong 等人（2010）及 Samantray 等人（2021）均為水分、纖維與時間依賴力學的關係提供研究背景。但是相對濕度 RH 是外部環境的水蒸氣狀態，並不等於材料內部即時含水率。紙袋從乾燥商店走到潮濕街道時，材料吸濕需要時間，局部層狀結構亦可能有水分梯度。因此，只用「RH=80%」替代整個歷程，最多只能作參數敏感度示範。
圖四採用完全公開而未經實測校準的代理映射。以 RH=50% 作基準，令兩個模量隨 RH 指數下降，τ 亦隨 RH 指數下降，再把變化後的參數代回 SLS。此設定刻意產生「濕度增加時，材料較軟且較快形成延遲應變」的研究情境；它不是從文獻取出的普遍定律，更不保證真實紙袋在所有 RH 區間均單調變形。若不同紙種的黏合劑遇水改變載荷路徑，或紙張出現不可逆纖維滑移，簡單 RH 代理便可能錯判位移與失效形式。
圖 4： 指定指數軟化與時間參數映射下，相對濕度由 30% 改至 90% 的示意伸長曲線。這是假設性代理計算，沒有用真實紙袋濕度資料校準；曲線不具獨立的經驗預測資格。
濕度敏感度實驗反而提醒研究者，不應把「濕度影響」當成可任意調高所有模型參數的理由。若未來取得實物資料，應量度紙張吸濕的時間常數及纖維方向性，然後估計 E0、E1、τ 如何隨含水率變化；若循環濕度時的蠕變高於恆定濕度情境，或許需要額外的歷程項，而非只改變當下模量。Gunderson 與 Tobey（1990）討論循環濕度改變速率的工作，正好提示環境歷程與單一 RH 數字有距離。
本文刻意不把 RH 加入蒙地卡羅的隨機變數，避免把尚未校準的代理假設與既定的六項對數常態分布混在一起。圖五的合成機率條件在參考環境下成立，而圖四純粹檢視另一組環境情境。若要聯合計算溫濕波動與負重變化，應先描述參數相關性、材料環境滯後及動態應力歷程，再作較大的合成試驗。
12. 研究結果四：四萬袋「虛構紙袋」的統計意義
以第七節分布產生四萬組參數，每組都計算質量、截面、受力角、模量與遲滯時間，再求第十分鐘位移。程式輸出平均為 1.025251 毫米，中位數 0.996148 毫米，樣本標準差約 0.248260 毫米；第五及第九十五百分位分別為 0.671331 和 1.475488 毫米。這些數字全部是指定假設生成機制下的分布摘要。其均值高於以平均參數計算的 0.967628 毫米，部分因截面積位於分母而出現非線性效應，說明「把平均參數代入函數」不等於「對函數的輸出取平均」。
圖五以歸一化直方圖顯示四萬次合成位移。對研究者自訂的 δ*=1.20 毫米界線，模擬超出比例為 0.218375，名義二項抽樣標準誤約 0.002066。在該獨立抽樣生成模型下，粗略正態近似 95% 蒙地卡羅誤差區間約為 0.21433 至 0.22242。這個區間只衡量有限模擬抽樣誤差，不是市面紙袋伸長比例的可信區間，也沒有涵蓋本構選錯、忽略黏合點或材料分布不明的結構性不確定性。
圖 5： 固定種子 2026102701、四萬次合成參數下的十分鐘預測伸長。虛線代表自訂 1.20 mm 位移界線；直方圖是合成數據，並非市場抽樣、事故紀錄或破斷試驗。抽樣規則與代碼見 analysis.py。
兩種相反解讀都會犯錯。一種說四萬次試驗證明有兩成紙袋太脆弱，忽略全部樣本都是由研究者假設的先驗抽取，連強度測試也沒有做。另一種說既然是合成數據便毫無價值，忽略它確實檢查了方程對參數變動的反應，揭示平均數代入偏差，並提供日後規劃試驗的可用代碼。正確表述必須同時肯定模型內推導的效度，及承認對現實世界外推的限度。
模擬沒有人工挑選符合預期的樣本，也沒有丟棄高伸長尾部。所有輸入值為正且角度遠離幾何奇異點，但對數常態分布的尾部形狀仍是人為指定。如果目的是分析實際破斷風險，應先建立經由實物校準的材料強度與缺陷分布，並可能需要極值統計或重要性抽樣。本文的位移超界百分比不具備那項證據基礎。
13. 模型辨識、內部效度與反例檢查
有人可能反駁：「手挽看似變長，只是紙帶在手指上滑動。」此反駁十分合理。肉眼看到的輪廓長度可能同時包含材料軸向應變、手提環彎曲部分被拉直、手指壓痕、接縫滑移及攝影視角。本文的定義只涵蓋八十毫米有效帶段的軸向伸長，刻意把其餘機制固定；因此就算數學完全正確，也只是其中一種解釋。若沒有在帶段兩端畫上標記並沿相同路徑量度長度，就不能從外觀確定軸向蠕變的實際比例。
另一個反例是材料應力與應變關係可能高度非線性。紙纖維在吸濕過程中可能有微觀重排、微裂紋或損傷累積。若模量隨應變下降，便可能在恒定重量下加速伸長，而本文的 SLS 只會預測伸長速度下降。若見到加速曲線，仍不能單憑該曲線判定是應力分配改變、黏合逐步失效還是材料本構變軟。需要用截面測量、施載控制及局部影像確認。本文故意保留容易反證的指數模型，並公開其適用邊界，不會把一切不確定性藏入任意黏度。
時間常數亦不是材料的唯一物理指紋。在有限時間窗內，兩個不同指數衰減支路可能近似一個指數；有限速率的提袋動作亦可能被誤擬合成較大的 τ。如果只觀察載重後十秒，E0 與 E1 可能互相補償，令參數估計不穩定。若量度器有零點偏移，即時彈性位移更容易被當作儀器偏差。合理的實驗應包括多個負載時間、多個卸載時段、不同小應力水平及重複測量。本文沒有虛構這些實驗或聲稱已完成真實模型辨識。
模型內部有三項可被程式驗證的必要條件。首先 SLS 載入時的 J(t) 對所有 t≥0 都非負，伸長隨時間不減；其次卸載瞬間下降的位移應等於 L0σ/E0；最後在遠超 τ 的恢復時間，SLS 殘餘必須接近零。summary.json 保存這些數值檢查。這些測試不能證明模型符合現實，但能防止程式與模型自身矛盾；漂亮的圖片不能替代最低限度的理論及數值自洽驗證。
從可否證性看，最有用的下一步不是單純把合成試驗次數加到一百萬，而是尋找可將模型分開的實驗。若卸載後位移立即回縮，之後按多個時間尺度恢復，便支持存在可恢復部分；若充分控制環境、長期等待後仍有殘餘，則有理由加入永久支路；若卸載後反而繼續伸長，便顯示本文的單調恢復本構不適合該試件。這些是對未來數據的條件性檢驗準則，不是現成的實驗發現。
14. 討論：重量、受力、時間及安全性不能混作一談
第一個常見陷阱是把外力不變理解成物件不能再形變。機械平衡描述各分支怎樣分擔重物，材料本構則描述給定應力歷程如何變成應變歷程；兩者屬不同層次。SLS 可以在每個時刻滿足力平衡，同時因內部黏壺逐步變形令手挽變長。就像帶阻尼的彈簧系統，外力停留在固定值後，內部配置仍可緩慢演化。數學沒有讓重物增加質量，而是推翻了「力必須立刻對應固定長度」的過度簡化直覺。
第二個陷阱是把全部時間效應稱為塑性。本文基準 SLS 的延遲位移即使在卸載後仍短暫存在，亦會逐漸回復。真正不可恢復的變形需要額外假設或損傷機制；單靠放下一小會兒的照片不能區分。反過來，局部黏合點滑動可能尚未令整條手挽變長很多便造成危險。以定力試驗、定長試驗和長期卸載追蹤分別量度材料響應，比單看最高伸長更有辨識價值。
第三個陷阱是把「袋內重量不變」誤寫成「每一條手挽受力不變」。質量固定只決定總重力，無法決定各分支張力；一條手挽滑離手指或局部脫膠後，其他分支可能承受更高負荷。即使各手挽完整，行走時的加速度與提拿角度變化也會改變瞬時張力。本文固定幾何只是為隔離蠕變效應；若要預測實際走路情境，應同步記錄手挽受力、角度、加速度與位移，再以動態載重取代恒定階躍。
第四個陷阱是把合成分布當作產品監管證據。圖五的超界比例是研究者指定模量、厚度與角度分布所得，沒有任何市場隨機樣本，更沒有測試破斷。界線 1.20 毫米也不是材料失效點。研究可以用曲線啟發實驗，例如選幾種帶寬、保持時間及控制濕度情況分組試驗；但不能據此宣稱某種紙袋合格、某品牌不安全或紙袋可負重至三點七公斤。
比較不同手挽時還需要留意長度標準化。相同應變 ε 乘以不同 L0 會產生不同絕對伸長，故一毫米對八十毫米長的材料與一百六十毫米長的材料並不等價。實驗報告至少應提供無量綱應變、標距及整個手提環幾何；若只把攝影上的下垂距離稱為材料蠕變，就容易把環形彎曲、手部變位與纖維伸長混在一起。將觀測量的意義說清楚，往往比再增加一個模型參數更能提升論文品質。
荒誕題目令人發笑，是因為購物者通常只問袋子放了多少公斤，卻少問已經拿了多久、紙帶有多薄、兩手挽是否等分受力、剛才有沒有淋雨。這些因素即使不改變重物質量，也可能改變看見的長度。但「有可行機制」與「某個實際案例確由該機制造成」有重要分別。本文達成前者，並明確保留後者所需的實驗工作；嚴謹的荒誕研究不需要藉誇張宣稱才能得到有趣結論。
15. 研究限制與有效範圍
本文最重要的限制是完全沒有實測紙袋。模量、時間常數、厚度、截面均勻性與黏流黏度全部由研究者假設，不能稱為實際材料常數。雖然引用了真實紙纖維和紙板文獻，但沒有把某文獻試件的數值直接移植到商業手挽。紙材纖維方向、回收比例、加工與黏合方式、孔洞、裂紋、含水率等都可能令模型失效。尤其實際手挽可以是扭紙繩、多層貼合帶或複合材料，未必具備本文假設的平面矩形截面。
其次，四段等張力模型忽略袋口柔性、手指支點、提袋轉角和行走時的動態負荷。局部剝離、剪切、孔洞尖端及折線應力集中可能主導真實破斷，需另行作結構試驗或有限元素分析。本文沒有破斷準則、材料強度分布或疲勞損傷方程，所以任何伸長預測都不能轉為安全荷載、事故概率或法規合規結論。
再次，SLS 是線性時間不變模型，面對大應變、多重遲滯時間、濕度循環及不可逆滑移可能過度簡化。探索性黏流支路只是一種可比較的數學形式，不能唯一識別永久變形的物理機制。圖四的濕度映射完全未校準，亦忽略材料內部水分傳輸。若要建立工程級預測，必須先以特定紙材的蠕變、卸載恢復及環境控制試驗校準。
蒙地卡羅還假設多項參數互相獨立，且以對數常態和均勻分布代表不明的材料差異；實際同廠紙袋的厚度、強度、纖維取向可能彼此相關。四萬次模擬抽樣可以把條件概率的計算誤差減小，但無法解決先驗分布錯誤。部分學術來源已核實作者、標題及 DOI，個別論點依據公開摘要而未逐篇全文審讀；因此文內刻意沒有從無法確認的文獻抄取手挽強度、實測伸長或隨機參數範圍。
16. 結論及未來研究
本研究回答了一項聽似荒誕的問題：紙袋沒有變重，手挽為何仍可愈提愈長？在可檢查的四段受拉近似及線性黏彈性模型下，恒定應力仍可產生延遲應變。示例 SLS 於提拿十分鐘後伸長 0.967628 毫米，卸載後立即回復 0.314656 毫米，十五分鐘後只剩約 0.032510 毫米可恢復殘餘；加入另一個假設性黏流分支，才產生 0.226552 毫米的永久分量。所有數字只屬示例模型計算。
相同質量可以因帶寬、厚度、提拿角度及材料參數改變而有不同位移。四萬次合成情境下，跨越自訂一點二毫米界線的比例為 0.218375，但絕不代表實際紙袋的破斷風險。後續應用標準化標距、定力裝置與多種保持時間量度真實位移曲線，並在固定及循環濕度下重複卸載恢復試驗。若要研究安全，須另測接縫剝離、局部缺陷、疲勞及最終破壞；不能把一條可恢復蠕變曲線當作安全認證。
參考文獻
Dong, F., Olsson, A.-M., & Salmén, L. (2010). Fibre morphological effects on mechano-sorptive creep. Wood Science and Technology, 44, 475–483. https://doi.org/10.1007/s00226-009-0300-3
Gunderson, D. E., & Tobey, W. E. (1990). Tensile creep of paperboard—Effect of humidity change rates. Materials Research Society Symposium Proceedings, 197, 213. https://doi.org/10.1557/PROC-197-213
Lin, C.-Y. (2024). Rethinking and researching the physical meaning of the standard linear solid model in viscoelasticity. Mechanics of Advanced Materials and Structures, 31(11), 2370–2385. https://doi.org/10.1080/15376494.2022.2156638
Olsson, A.-M., Salmén, L., Eder, M., & Burgert, I. (2007). Mechano-sorptive creep in wood fibres. Wood Science and Technology, 41, 59–67. https://doi.org/10.1007/s00226-006-0086-5
Olsson, A.-M., & Salmén, L. (2014). Mechano-sorptive creep in pulp fibres and paper. Wood Science and Technology, 48, 569–580. https://doi.org/10.1007/s00226-014-0624-5
Samantray, P., Massart, T. J., Peerlings, R. H. J., & Geers, M. G. D. (2021). Modeling the effect of creep in paper fibres under the influence of external loading and changes in moisture. Mechanics of Materials, 163, 104075. https://doi.org/10.1016/j.mechmat.2021.104075
Ward, I. M., & Sweeney, J. (2012). Linear viscoelastic behaviour. In Mechanical Properties of Solid Polymers (Chapter 5). Wiley. https://doi.org/10.1002/9781119967125.ch5