一、經(jīng)典力學(xué)的"近視"問題：把材料當(dāng)成無限可分的點(diǎn)

經(jīng)典的固體力學(xué)建立在一個(gè)看似合理的假設(shè)上：材料是連續(xù)的，可以被無限分割成沒有內(nèi)部結(jié)構(gòu)的"材料點(diǎn)"。

這個(gè)假設(shè)在宏觀世界非常成功——計(jì)算大橋變形、飛機(jī)機(jī)翼應(yīng)力都很準(zhǔn)確。但當(dāng)我們把目光投向微納米尺度（MEMS傳感器、微納電子器件）或應(yīng)變集中問題時(shí)，奇怪的事情發(fā)生了：

• 微懸臂梁：厚度從8μm減到2μm，測(cè)得的彈性模量從115 GPa飆升到175 GPa（變硬了50%）帶小孔的板：孔徑從5mm減到0.4mm，應(yīng)變集中系數(shù)從3.0降到1.2（按理說應(yīng)該不變）微壓痕測(cè)試：壓得越淺，算出來的"彈性模量"越大（著名的壓痕尺寸效應(yīng)）

經(jīng)典理論完全無法解釋這些現(xiàn)象，因?yàn)樗暮诵募僭O(shè)出了問題：真實(shí)的RVE（代表性體積單元）從來不是無限小的。

二、均質(zhì)化誤差的本質(zhì)：平均值≠中心值

關(guān)鍵洞察：當(dāng)物理場(chǎng)在RVE內(nèi)非線性分布時(shí)，體積平均值不等于幾何中心處的真實(shí)值。

直觀理解

想象你測(cè)量一個(gè)房間的溫度：

• 經(jīng)典做法：假設(shè)溫度在房間內(nèi)均勻分布，用房間中心的溫度代表整個(gè)房間真實(shí)情況：如果暖氣片在一側(cè)，溫度呈梯度分布，平均溫度≠中心溫度

數(shù)學(xué)上，這可以通過泰勒展開描述。對(duì)于任意物理量 p （可以是應(yīng)變、應(yīng)力或應(yīng)變能密度）：

其中 h 是RVE的尺寸（比如金屬的晶粒尺寸、聚合物的高分子鏈回轉(zhuǎn)半徑），Δp 是拉普拉斯算子（描述場(chǎng)的"彎曲程度"）。

關(guān)鍵發(fā)現(xiàn)：經(jīng)典理論只保留了第一項(xiàng)，忽略了  和  項(xiàng)——這就是"均質(zhì)化誤差"的來源。

三、高階應(yīng)變能：一個(gè)參數(shù)解釋兩種矛盾現(xiàn)象

解決方法：重新定義均質(zhì)化應(yīng)變能密度： 

其中  是經(jīng)典應(yīng)變能，  是高階修正項(xiàng)，包含應(yīng)變與其梯度的交叉乘積項(xiàng)，如  。

為什么能解釋"變硬"和"變軟"？

這些交叉項(xiàng)的符號(hào)取決于應(yīng)變的空間分布：

變形模式交叉項(xiàng)符號(hào)材料表現(xiàn)典型案例微梁彎曲	正（高階項(xiàng)為正）	硬化	厚度↓ → 剛度↑
小孔附近拉伸	正	硬化	孔徑↓ → 應(yīng)力集中↓
微壓痕中心	負(fù)（高階項(xiàng)為負(fù)）	軟化	壓深↑ → 表觀模量↓

物理圖像：

• 硬化：高階項(xiàng)"幫正忙"，讓材料"更難變形"，需要更多能量軟化：高階項(xiàng)"幫倒忙"，讓變形"更容易發(fā)生"，儲(chǔ)存能量更少

這與傳統(tǒng)應(yīng)變梯度理論（如Mindlin理論）有本質(zhì)區(qū)別——傳統(tǒng)理論需要引入多個(gè)本構(gòu)參數(shù)（18個(gè)分量的六階張量、54個(gè)分量的八階張量），而新理論只需要一個(gè)參數(shù) h （RVE尺寸），且有明確的物理意義。

四、數(shù)值實(shí)現(xiàn)：增廣拉格朗日法的巧妙之處

高階理論的最大障礙是求解困難。控制方程包含四階甚至六階微分算子，傳統(tǒng)有限元需要  或  連續(xù)性的形函數(shù)（極其復(fù)雜）。

論文采用了增廣拉格朗日方法（ALM）的巧妙策略：

1. 引入輔助變量：將應(yīng)變  作為獨(dú)立變量，不再通過位移求導(dǎo)得到約束條件：通過拉格朗日乘子 λ 和懲罰項(xiàng) r ，強(qiáng)制滿足 降階優(yōu)勢(shì)：最高階導(dǎo)數(shù)降為二階，普通二次單元即可求解

最終形成15×15的分塊剛度矩陣（每個(gè)節(jié)點(diǎn)3個(gè)位移+6個(gè)應(yīng)變+6個(gè)拉格朗日乘子），雖然自由度增加，但避免了復(fù)雜的 C1 單元和迭代處理高階邊界條件的麻煩。

五、實(shí)驗(yàn)驗(yàn)證：三個(gè)典型案例

案例1：微懸臂梁的尺寸效應(yīng)

• 預(yù)測(cè)：厚度從8μm→2μm，名義楊氏模量從115 GPa→175 GPa實(shí)驗(yàn)：Choi等人的純銅微梁實(shí)驗(yàn)數(shù)據(jù)吻合良好

名義彈性模量隨梁厚度的變化

• 機(jī)制：高階應(yīng)變能占比從6.5%增至15.7%，整體表現(xiàn)為硬化中性面不“中性”

中性面區(qū)域的高階應(yīng)變能占比

案例2：帶孔板的應(yīng)變集中

• 經(jīng)典理論：應(yīng)力集中系數(shù)恒為3，與孔徑無關(guān)實(shí)驗(yàn)觀測(cè)：孔徑越小，應(yīng)變集中系數(shù)越小

實(shí)驗(yàn)觀測(cè)與新理論預(yù)測(cè)的應(yīng)變集中因子（倒數(shù)）

• 新理論預(yù)測(cè)：與DIC（數(shù)字圖像相關(guān)）實(shí)驗(yàn)測(cè)量吻合

應(yīng)變分布隨孔徑減小逐漸趨向于平緩

案例3：微壓痕的"軟化"現(xiàn)象

• 反常現(xiàn)象：壓痕越深，反算出的彈性模量越小機(jī)制解釋：壓痕下方剪切應(yīng)變的拉普拉斯為負(fù)，導(dǎo)致交叉項(xiàng)為負(fù)，高階能量"減負(fù)"，整體表現(xiàn)為軟化

六、理論意義與工程價(jià)值

與傳統(tǒng)高階理論的對(duì)比

特征Mindlin應(yīng)變梯度理論新理論（均質(zhì)化應(yīng)變能）尺度參數(shù)	多個(gè)	1個(gè)（RVE尺寸 h）
物理意義	擬合參數(shù)	明確的微觀結(jié)構(gòu)尺度
硬化/軟化	僅能硬化	同時(shí)描述兩種行為
本構(gòu)關(guān)系	廣義應(yīng)力-應(yīng)變梯度	能量密度直接定義
邊界條件	高階項(xiàng)復(fù)雜難處理	通過ALM自然處理

工程應(yīng)用前景

1. MEMS器件設(shè)計(jì)：準(zhǔn)確預(yù)測(cè)微結(jié)構(gòu)的剛度，避免過度保守設(shè)計(jì)增材制造：考慮打印微觀結(jié)構(gòu)尺度的力學(xué)行為優(yōu)化生物力學(xué)：骨骼、軟組織等具有特征微結(jié)構(gòu)的材料建模壓痕測(cè)試：正確解讀微納壓痕數(shù)據(jù)，區(qū)分真實(shí)材料參數(shù)與尺寸效

七、從"數(shù)學(xué)點(diǎn)"回到"物理實(shí)體"

回歸物理本質(zhì)：經(jīng)典力學(xué)的失效不是因?yàn)閿?shù)學(xué)不夠復(fù)雜，而是因?yàn)楹雎粤?材料點(diǎn)"實(shí)際上有有限尺寸這一基本事實(shí)。

通過泰勒展開揭示的均質(zhì)化誤差，我們理解了為什么不均勻變形會(huì)觸發(fā)尺寸效應(yīng)——當(dāng)應(yīng)變梯度與RVE尺寸相當(dāng)時(shí)，"平均"和"真實(shí)"的差異變得不可忽略。

更重要的是，新理論用一個(gè)參數(shù)統(tǒng)一了看似矛盾的實(shí)驗(yàn)現(xiàn)象，這符合物理學(xué)的美學(xué)追求。正如論文所指出的，這種方法可以進(jìn)一步擴(kuò)展到彈塑性和損傷問題，為計(jì)算材料力學(xué)開辟新的道路。

參考文獻(xiàn)：Yuheng Cao, Chunyu Zhang*, Biao Wang. A New High-order Deformation Theory and Solution Procedure Based on Homogenized Strain Energy Density. International Journal of Engineering Science, 195 (2024) 103990.