MATS42802 / MATS68602 复习中心

Polymer Physics & Physical Properties · Elastomers · Viscoelasticity · Glass Transition

📊 测验得分:0 / 0

📋 课程总览 / Course Overview

核心问题 / The Big Question: 为什么橡胶有弹性?为什么聚合物熔体既像液体又像固体?
核心答案:都来自构象熵 (entropy)。这是整门课的第一性原理。
所有公式都源自一句话:"拉伸聚合物链 → 减少构象数 → 减少熵 → 自由能 F = U - TS 升高 → 产生回复力 F = -T(dS/dl)"。

八大主题模块

1. Chain Statistics
链统计 (Week 4)
点击翻转
核心: 自由连接链 (FJC) → 末端距 R = a√N。引入特征比 C∞ = b²/a² 处理真实链。
修正: 自规避走 ν=0.6 (Flory);熔体中回到 ν=0.5。
2. Entropic Spring & Thermo
熵弹簧 & 热力学
点击翻转
核心: Helmholtz free energy A=U-TS。橡胶 dU≈0,弹力纯熵起源。
$k_{spring} = 3k_BT/(Nb^2)$;T 升高 → 弹簧变硬!
3. Classical Rubber
经典橡胶弹性
点击翻转
核心: 仿射变形 + 不可压缩。
$G = nk_BT = \rho RT/M_x$, $E=3G$, $\tau = G(\lambda - 1/\lambda^2)$
4. Balloon Problem
气球问题(普适)
点击翻转
核心: 双轴拉伸,λx·λy·λz=1。压力峰值在 r/r₀=7^(1/6)≈1.38,与材料无关 — 普适!
5. Corrections
修正理论
点击翻转
核心: Flory 末端修正、Mooney-Rivlin 图、应变诱导结晶(大 λ 硬化)、Flory-Rehner 溶胀理论。
6. Viscoelasticity
粘弹性
点击翻转
核心: 应力弛豫 G(t)、储能 G' & 损耗 G''、loss tangent tanδ=G''/G'。
$G^* = G' + iG''$
7. Reptation
蛇形蠕动
点击翻转
核心: de Gennes/Doi-Edwards 管模型。
$G_p=\rho RT/M_e$, $\tau_T \sim N^3$ (理论) vs $N^{3.4}$ (实验)。$D \sim 1/N^2$
8. WLF & Tg
时温叠加 & 玻璃化
点击翻转
核心: WLF ↔ Vogel-Fulcher 等价(数学上完全相同)。
$C_2 = T_{ref}-T_0$, $B = 2.303 C_1 C_2$。Tg 不是真正相变。
💡 备考策略 / Exam Strategy:
1. 原题背熟:Week 5/6/7 三套题目集就是考试题型。本网页每章末尾都有详细分步解答。
2. 公式 ≈ 12 个:先掌握"公式表"标签 — 那是考前 30 分钟最有用的。
3. 逻辑链:构象数 ↘ → 熵 ↘ → 自由能 ↗ → 回复力 ↗。所有橡胶弹性都从这里推导。
4. 单位陷阱:温度永远用 K;ρ 用 kg/m³;用 R 时 Mc 用 kg/mol,用 k_B 时 n 用 1/m³。
📚 课程脉络 / Course Roadmap:
Week 4: 单链统计与熵 → Week 5: 经典橡胶弹性公式 + AFM 单链拉伸 + 气球问题 → Week 6: 实验偏差与修正 + Flory-Rehner 溶胀 → Week 7: 粘弹性、reptation、WLF/玻璃化。
本课的核心进展:从一根链 → 一个网络 → 实际偏差 → 动态行为 → 玻璃态

📔 背景知识 / Background (Week 2 Lecture 4)

本章涉及 Week 2 handouts lecture 4。 不会直接出计算题,但是选择题 & 论述题常考的概念背景。

0.1 聚合物物理状态总览 / Five States of Polymers

状态描述例子
液体(熔体)无规缠结,可流动 (粘弹液体)加工中的塑料
玻璃无规结构,但运动冻结有机玻璃 (PMMA)
半结晶固体晶区+非晶区PE, PP
交联橡胶网络 + 链段在 Tg 之上自由橡皮筋
交联玻璃网络 + 链段冻结环氧树脂

0.2 聚合物溶液浓度

稀溶液 半稀 浓溶液 熔体 链不重叠 开始重叠 均匀分布 无溶剂
关键现象 — 屏蔽 (Screening):稀溶液中链膨胀(ν=0.6),但浓溶液和熔体中由于浓度梯度消失 → excluded volume 被屏蔽 → 熔体中链回到理想 random walk (ν=0.5)
Higgins 用中子散射在 1974 年证实了这一点(让部分链氘代以增强散射)。

0.3 网络类型 / Network Types

Vulcanization
硫化(化学交联)
点击翻转
Charles Goodyear 1839 发现。
在双键 (-CH=CH-) 处形成 -S-S- 桥。
需要不饱和键(天然橡胶、丁苯橡胶)。
不可逆。
Newer Crosslinking
现代交联方法
点击翻转
过氧化物、紫外光、辐射、Click 化学。
优点:低温操作、可控、不需双键。
用于硅胶 (PDMS) 等饱和高分子。
Thermoplastic Elastomer
热塑性弹性体
点击翻转
三嵌段共聚物 (ABA):
• 端段:玻璃态 (PS, Tg=100°C)
• 中段:橡胶态 (聚丁二烯)
玻璃微相分离 → 物理交联点。
加热 > Tg → 流动!可循环利用。
Examples
商品名
点击翻转
Kratons:苯乙烯-丁二烯-苯乙烯三嵌段
Lycra (Spandex):聚氨酯多嵌段 (AB)n
SEBS, SBS:弹性鞋底
Vitrimers
玻璃聚合物(可重塑)
点击翻转
Leibler 2011 发现。
交联键可"换伴侣"(酯交换反应),但密度不变
结果:高温可塑形而不熔化 — 像玻璃但是聚合物。
Hydrogel
水凝胶
点击翻转
交联聚合物 + 大量水。
典型:果冻 (gelatin)、隐形眼镜、高吸水材料。
溶胀比 = 大量水 / 干聚合物体积,可达 1000×。

0.4 Manchester 历史人物

📜 这门课的"曼彻斯特血统":
John Gough (1805) — 首次观察橡胶受热收缩
John Dalton (Gough 朋友) — 原子论
James Joule (1859) — 量化拉伸放热(用望远镜)
Geoffrey Gee (UMIST 1953) — 经典橡胶弹性理论(与 Treloar、Flory 共同)
Sam Edwards (Manchester 1958-72) — Reptation 理论奠基人
Julia Higgins (Manchester 1968-72) — 中子散射证实熔体中 ν=0.5

0.5 选择题

1在熔体中,聚合物链的尺寸标度律是?
D:浓溶液/熔体中 excluded volume 被屏蔽 → 链回到 ν=0.5。这是 Edwards 1965 的洞察,Higgins 用中子散射证实。
2热塑性弹性体(如 Lycra)能加热熔化的原因?
B:Kratons 等的端段是 PS (Tg=100°C)。低温下 PS 微相分离形成"交联",加热 > 100°C 后玻璃区液化 → 整体可流动。

📘 第1章 链统计 / Polymer Chain Statistics

本章涉及 Week 4 lectures 1-3 + Week 5 Question 1-2。 核心问题:聚合物链在空间中是什么形状?我们能否用一个简单参数描述?

1.1 三个层次的模型 / Three Levels of Models

Freely Jointed Chain (FJC)
自由连接链 — 最简模型
点击翻转
N 个长度为 a 的刚性链段,彼此夹角无限制。
$\langle R^2 \rangle = N a^2$,   $\langle R^2 \rangle^{1/2} = a\sqrt{N}$
"显然不真实,但居然管用" — Toolan
Real Chain with C∞
真实链(特征比修正)
点击翻转
真实链有键角限制,但短程关联可吸收进 b:
$\langle R^2 \rangle = N b^2$,   $C_\infty = b^2/a^2$
所有化学复杂性 → 一个数 C∞!
Self-Avoiding Walk (SAW)
自规避随机游走
点击翻转
链段不可重叠(excluded volume)→ 链膨胀。
$\langle R^2 \rangle^{1/2} \sim aN^\nu$, $\nu \approx 0.6$
Flory 估算:ν=3/5;精确:0.588。良溶剂下成立。
Chain in Melt
熔体中的链
点击翻转
惊人结论:在熔体中,excluded volume 被相邻链"屏蔽",链回到理想 random walk!
$\langle R^2 \rangle^{1/2} = b\sqrt{N}$ 仍然成立。
Julia Higgins 用中子散射证实了这一点。

1.2 关键概念 / Concepts

RMS End-to-End Distance
均方根末端距
点击翻转
$\langle R^2 \rangle^{1/2} = a\sqrt{N} = b\sqrt{N_{eff}}$
R ~ √N,不是 N — 链是卷曲的随机游走。
例:PS Mw=100k,R≈21 nm,但轮廓长 L≈202 nm。
Statistical Segment Length b
统计步长
点击翻转
"等效链段"长度,b > a 因为局部刚性。
PS: a=0.21 nm, b=0.67 nm。
所有链都可以视为"自由连接",只要用 b 替代真实键长 a。
Characteristic Ratio C∞
特征比
点击翻转
$C_\infty = b^2/a^2$
典型:PEO=6.9, PS=10.2, P3HT=16, DNA=350
大 → 链刚(共轭、芳香环、双链);小 → 链柔(醚键)
Contour Length
轮廓长度
点击翻转
完全伸展长度 L = N · a(用真实键长 a!)
注意:L 与 R 比值 = √N,对长链相差很大。
PS Mw=100k:L=202 nm 但 R=21 nm。
Flory Exponent
Flory 指数(自规避)
点击翻转
最小化总自由能 A = A_rep + A_el:
R³ ~ N²/R² → R ~ N^(3/5)
Flory 估算 ν=0.6;重整化群得 0.588。
R 比理想链膨胀!
Wiener Chain
连续高斯链(Wiener)
点击翻转
把 N → ∞ 而 Nb² 固定 → 链变光滑。
概率分布满足扩散方程:$\partial P/\partial N = (b^2/6)\nabla^2 P$
是 Edwards 场论的基础。

1.3 关键公式 / Key Formulas

$$\langle R^2 \rangle = N a^2 \quad\Rightarrow\quad \langle R^2 \rangle^{1/2} = a\sqrt{N}$$
FJC 模型基本结果
$$\langle R^2 \rangle = \frac{N}{g}\langle c^2 \rangle = N b^2, \quad C_\infty = b^2/a^2$$
真实链:把 g 个链段合并成一个长度 b 的"等效段"
$$N = M_w / M_r, \qquad L_{\text{contour}} = N \cdot a$$
N 由分子量定义;轮廓用 a (真实键长)
$$P(\mathbf{R},N) \propto \exp\left[-\frac{3R^2}{2Na^2}\right]$$
末端距高斯分布 — 通过组合统计或扩散方程导出
$$\langle R^2\rangle^{1/2} \sim aN^\nu, \quad \nu \approx 0.6 \text{ (good solvent)}, \quad \nu = 0.5 \text{ (theta/melt)}$$
考虑 excluded volume 的 Flory 标度
💡 记忆点:
R ~ √N,不是 N — 链是高度卷曲的。
• C∞ 排序:DNA(350) >> P3HT(16) > PS(10.2) > PP(6.8) > PEO(6.9)
• 良溶剂中链膨胀 (ν=0.6),但熔体中又回到理想链 (ν=0.5) — 这是 Edwards 的洞察。
• PS 与 PP 比较:PS 有苯环,b 更大 → R 更大。

1.4 直接计数法 — 一维随机游走 / Direct Enumeration

📚 推导思路:一维随机游走,N 步中 N+ 前进、N- 后退。组合数 W = N!/(N+!N-!)。
用 Stirling 近似 → ln W = -2N(f - 1/2)²,f = N+/N。
高斯分布!直接计数验证了高斯结果(不靠扩散方程)。

1.5 图解 / Diagram

柔性链 (PEO, C∞≈7) 半刚性 (P3HT, C∞≈16) 刚性 (DNA, C∞≈350) 小 R, 紧密线团 中等 R 大 R, 棒状 同 N 不同 C∞ → 不同链形态
特征比越大,链越刚 → 末端距越大

1.6 选择题 / Multiple Choice

1N 翻倍时,理想链的 RMS 末端距变为原来的几倍?
B:$\langle R^2 \rangle^{1/2} = a\sqrt{N}$,N 翻倍 → R × √2。
2关于自规避随机游走 (SAW),下列哪个错误?
C:熔体中,相邻链的密度抵消了排除体积,链回到理想 random walk (ν=0.5)。这是 Edwards 与 Higgins 的洞察。
3DNA 的 C∞ ≈ 350 远大于柔性聚合物,主因是:
D:双螺旋让相邻链段几乎共线 → b 远大于 a。
4两个 PE (a=0.154 nm, C∞≈7) 链:A 的 N=100,B 的 N=400。RMS 末端距比 R(B)/R(A) =?
B:R = b√N,N×4 → R×2。

1.7 计算题 / Calculation Problems

原题 · Week 4 Q1 (Live Session)

聚苯乙烯链尺寸 / PS Chain Dimensions

PS Mw=100,000,Mr=104, C∞=10.2, b=0.67 nm。
(a) RMS 末端距;(b) 轮廓长度;(c) 同样 L 的 PP (Mr=42, C∞=6.8) 分子量;(d) 此 PP 的 RMS。(e) 为什么 (d) ≠ (a)?

(a) N = 100,000/104 = 962
$\langle R^2\rangle^{1/2} = b\sqrt{N} = 0.67 \times \sqrt{962} =$ 20.8 nm ≈ 21 nm
(b) 真实键长 $a = b/\sqrt{C_\infty} = 0.67/\sqrt{10.2} = 0.21$ nm
$L = N \cdot a = 962 \times 0.21 =$ 202 nm
(c) PP 同样 L → 同样的骨架原子数 N。
Mw(PP) = N × Mr(PP) = 962 × 42 = 40,404 g/mol
关键洞察:完全伸展长度相同时,分子量可能差很多!
(d) $b_{PP} = a\sqrt{C_\infty(PP)} = 0.21 \times \sqrt{6.8} = 0.547$ nm
$R = 0.547 \times \sqrt{962} =$ 17 nm
(e) PP 的 b 较小(C∞=6.8 < 10.2)→ 同样 N 下链更"紧凑"。FJC 模型把所有化学差异封装在 C∞ 这一个参数里。
💡 PS 和 PP 同 L 但 R 不同,说明侧基化学影响链的"舒展度" — 苯环让 PS 更胖。
原题 · Week 5 Q2 (AFM 单链拉伸)

AFM 单链 PS 拉伸

PS Mw=100,000,θ-溶剂,T=298 K。(a) 等效熵弹簧常数;(b) 拉伸 100 nm 时探针上的力。

(a) $k_{spring} = \dfrac{3k_BT}{Nb^2}$  (N=962, b=0.67 nm)
$= \dfrac{3 \times 1.38\times10^{-23} \times 298}{962 \times (0.67\times10^{-9})^2}$
$= \dfrac{1.234\times10^{-20}}{4.32\times10^{-16}} \approx$ 2.9 × 10⁻⁵ N/m
(b) $F = k \cdot x = 2.9\times10^{-5} \times 100\times10^{-9}$ = 2.9 pN
💡 单链力的尺度 ~ pN,AFM 灵敏度 ~ 1 pN,所以这种实验可行。
注意:高斯近似仅在拉伸 << 轮廓长度 (202 nm) 时有效,100 nm 接近上限。
原题 · Week 5 Q1

排序 P3HT、PEO、DNA 的特征比

排序:DNA (350) >> P3HT (16) > PEO (6.9)
PEO -O-CH₂-CH₂- 醚键非常柔,C∞ 接近最低值。
P3HT 噻吩环共轭,链有半刚性 — 这是它能形成晶体的原因。
DNA 双螺旋几乎像棒,持续长度 ~50 nm,因此 C∞ 极大。
💡 一般规律:刚性来源 = 共轭、芳香环、氢键、电荷排斥、二级结构。

📗 第2章 熵弹簧 & 热力学 / Entropic Spring & Thermodynamics

本章涉及 Week 4 lectures 1-2 的热力学框架。 核心问题:橡胶弹性来源是什么?为什么 James Joule 1859 年发现拉伸橡胶会发热?

2.1 热力学框架 / Thermodynamic Framework

$$A = U - TS \quad\text{(Helmholtz free energy)}$$
在恒温恒容下,做功 dW = dA
$$F = \left(\frac{dA}{dl}\right)_{T,V} = \left(\frac{dU}{dl}\right)_{T,V} - T\left(\frac{dS}{dl}\right)_{T,V}$$
力 = 内能贡献 + 熵贡献

普通晶体(弹性来自 U)

$U(x) \approx U_0 + \frac{1}{2}sx^2$,   $s = d^2U/dr^2|_{r_0}$
原子在势阱中振动,dS ≈ 0。
升温 → 振动增大 → 软化。

橡胶(弹性来自 -TS)

$F = -T(dS/dl)$,   dU ≈ 0
链段无固定平衡位置,但有构象数
升温 → 熵效应放大 → 硬化!

2.2 关键概念 / Key Concepts

Entropy of a Chain
链的构象熵
点击翻转
由 P(R) 推出 S = k_B ln P + const:
$$S(R) = -\frac{3k_B R^2}{2Na^2} + S_0$$
R 越大,构象数越少 → S 越小。
Free Energy of a Chain
链的自由能
点击翻转
$$A(R) = -TS = \frac{3k_B T R^2}{2Na^2} + \text{const}$$
抛物线形 → 像弹簧的弹性能 ½kx²!
Spring Constant
熵弹簧常数
点击翻转
$k = \dfrac{d^2A}{dR^2} = \dfrac{3k_BT}{Nb^2}$
三特点:
• ∝ T (升温变硬!)
• ∝ 1/N (链越长越软)
• ∝ 1/b² (链段越长越软)
Joule Effect (1859)
焦耳效应
点击翻转
拉伸橡胶 → ΔS<0 → Q=TΔS<0(放热)。
1859 年 Joule 在 Manchester 用望远镜读取微小的水中橡胶位移变化,证实了这一点。
Heat Anomaly
温度反常行为
点击翻转
悬挂的橡胶受热收缩(拉直时模量增大)!
普通材料热膨胀。
这是熵弹性的标志性现象。
Liquid-like Locally
局部液体性
点击翻转
在橡胶态,链段在交联点之间自由运动,像液体一样。
低温时 → 链段冻结 → 玻璃态。
这是 segments → bulk 的桥梁。

2.3 图解:拉伸 → 构象 ↘ → 熵 ↘ → 力 ↗

未拉伸 (高熵) 大量构象 F = 0 → 拉伸 → 力 F = -T·dS/dl 拉伸 (低熵) 少构象 F = k·dl > 0 系统偏好高熵 → 产生回复力
💡 五个反直觉记忆点:
1. 拉伸橡胶放热(不吸热)
2. 橡胶受热收缩(不膨胀)
3. 弹簧常数 ∝ T,升温变硬(金属反向)
4. 拉伸时内能几乎不变(dU≈0)
5. 橡胶模量来自构象数减少,不来自键拉伸

2.4 选择题

1等温拉伸橡皮筋时
A:ΔS<0 → dQ=TΔS<0。Joule 1859 用望远镜实验首次证实。把拉伸橡皮筋按在唇上能感到温热。
2温度升高,单链熵弹簧的 k 怎么变?
B:k = 3k_BT/(Nb²) ∝ T。这是橡胶相对普通固体的反常行为。
3用 F = (dU/dl) - T(dS/dl) 描述橡胶时
C:橡胶的 dU/dl ≈ 0(链段无内能变化),所有弹力来自熵。

📙 第3章 经典橡胶弹性 / Classical Rubber Elasticity

本章涉及 Week 5 lectures 4-5 + Week 6 Question 1。 把单链熵弹簧推广到整个网络,得到宏观应力-应变关系。

3.1 关键概念

Affine Deformation
仿射变形假设
点击翻转
宏观变形 (λx, λy, λz) 等比例传递到每个交联点:
(x,y,z) → (λx·x, λy·y, λz·z)
关键简化:忽略链间相互作用。
Incompressibility
不可压缩
点击翻转
$\lambda_x \lambda_y \lambda_z = 1$
橡胶局部像液体。
单轴:λ₁=λ, λ₂=λ₃=1/√λ
双轴(气球):λ₁=λ₂=r/r₀, λ₃=(r₀/r)²
Entropy of Network
网络变形熵
点击翻转
对所有链求和(每单位体积):
$$\Delta S/V = -\frac{nk_B}{2}(\lambda_x^2 + \lambda_y^2 + \lambda_z^2 - 3)$$
n = 单位体积内交联链数。
Shear Modulus G
剪切模量
点击翻转
$G = nk_BT = \dfrac{\rho RT}{M_x}$
由变形熵 → 应力 → 模量。
注意:G ∝ T,G 与化学无关!
Stress-Extension Relation
应力-拉伸关系
点击翻转
$\tau = G\left(\lambda - \dfrac{1}{\lambda^2}\right)$
用工程应变 e: $\tau = nk_BT[(1+e) - (1+e)^{-2}]$
非线性!λ=1 时 τ=0。
Young's Modulus E
杨氏模量
点击翻转
不可压缩 (ν=0.5): $E = 2G(1+\nu) = 3G$
注意:一般材料 E = 2G(1+ν),橡胶 ν=1/2 是特例。
对小应变可用 σ = E·ε 直接算压缩力。

3.2 公式合集

公式说明
$G = nk_BT = \rho RT / M_x$剪切模量(基本)
$E = 3G = 3 nk_BT$杨氏模量(不可压缩)
$\tau = G(\lambda - 1/\lambda^2)$用 λ 表示
$\tau = nk_BT[(1+e) - 1/(1+e)^2]$用工程应变 e
$\Delta S/V = -\frac{nk_B}{2}(\lambda_x^2+\lambda_y^2+\lambda_z^2-3)$变形熵
$\lambda_x \lambda_y \lambda_z = 1$不可压缩约束
💡 记忆点:
E = 3G(不可压缩特例,比一般 E = 2G(1+ν) 简单)
• G ∝ T(升温变硬!)
• G 由交联密度决定,不依赖具体化学
• 大应变下偏离经典理论:strain hardening (>40%)

3.3 计算题 — 经典原题

原题 · Week 6 Q1

聚丁二烯交联橡胶(重要!)

立方体 100×100×100 mm,T=23°C,ρ=0.9 g/cm³。
(a) 上面施加剪切力 500 N,位移 4.46 mm,求 G;
(b) 求 Mc;
(c) 求交联点平均距离;
(d) 单轴压缩 4.46 mm 需要多大力?

(a) A = 0.1² = 0.01 m²
τ = 500/0.01 = 50,000 Pa = 50 kPa
γ = 4.46/100 = 0.0446
G = 50,000/0.0446 = 1.12 × 10⁶ Pa = 1.12 MPa
(b) $M_c = \dfrac{\rho RT}{G} = \dfrac{900 \times 8.314 \times 296}{1.12 \times 10^6}$
= 1980 g/mol ≈ 2000 g/mol
(单位:ρ=900 kg/m³, T=23+273=296 K)
(c) 正确方法(Live Session 答案):用交联密度 n。
$n = \dfrac{G}{k_BT} = \dfrac{1.12\times10^6}{1.38\times10^{-23}\times296} = 2.74 \times 10^{26}$ m⁻³
平均距离 $d = (1/n)^{1/3} = (3.65\times10^{-27})^{1/3} \approx$ 1.54 nm(≈ 15 Å)
不要用 R = b√N — 那是同一根链上的 RMS 末端距,不是交联点之间的实空间距离!
(d) 线性近似(标准方法):不可压缩 ν=0.5 → E = 2G(1+0.5) = 3G
ε = 4.46/100 = 0.0446
F = E × A × ε = 3 × 1.12×10⁶ × 0.01 × 0.0446 = 1500 N ≈ 1.5 kN
口诀:压缩力 = 3 × 剪切力 = 3 × 500 = 1500 N
非线性大应变公式(参考):τ = G(λ - 1/λ²),λ=0.9554 → F=1572 N(仍接近)
💡 (c) 用 n=G/kT 得 1.54 nm;(d) 用 E=3G 线性公式得 1500 N。这是 Live Session 标准答案。
模拟题

剪切模量与温度关系

在 0°C 测得橡胶 G = 0.85 MPa。在 80°C 同一样品的 G 是多少?

G ∝ T(绝对温度)→ G(353 K)/G(273 K) = 353/273 = 1.293
G(80°C) = 0.85 × 1.293 ≈ 1.10 MPa
💡 升温 80 K → G 升高约 30%。这是橡胶的标志性反常行为。

3.4 拉伸 vs 剪切 — 重要对比

🔑 核心区别:橡胶在拉伸非线性,但在剪切胡克定律(线性)!

剪切几何推导(Week 2 lecture 5):
剪切应变 γ = tan φ,其中 φ 是剪切角。利用恒等式 tan(45°-φ/2) = (1-tan(φ/2))/(1+tan(φ/2))
解出拉伸比 λ = 1/√(1+γ²/4 ... ) — 简化得:
$\lambda^2 - 1/\lambda^2 \approx (\lambda - 1/\lambda)(\lambda + 1/\lambda) \approx \tan^2\phi$
代入熵公式:
$$\Delta S/V = \frac{nk_B}{2}\left[\lambda^2 + \frac{1}{\lambda^2} - 2\right] = \frac{nk_B}{2}\tan^2\phi$$
剪切下熵变(用 γ_shear = tanφ 表达)
$$\boxed{\tau_{shear} = G \cdot \gamma_{shear} = nk_BT \cdot \gamma}$$
胡克定律!即使大剪切应变也成立。剪切模量与小拉伸时相同。
💡 记忆要点:
• 拉伸:τ = G(λ - 1/λ²) — 非线性(10 模量值范围)
• 剪切:τ = G·γ — 胡克定律,对所有 γ 都线性!
• 这是橡胶的特殊行为,与几何对称有关。

3.5 拉伸时的绝热温升 / Adiabatic Heating (Joule 效应定量)

Joule 1859 量化测量:拉伸橡胶时温度升高。给定绝热条件,所有外力做功转为温度上升:
$W_{stretch} = m \cdot c_p \cdot \Delta T$
用力公式 $F = nk_BT(λ - 1/λ²) \cdot A_0$ 积分得:
$$\Delta T = \frac{nk_B T}{2 \rho c_p}\left[\lambda_2^2 + \frac{2}{\lambda_2} - \lambda_1^2 - \frac{2}{\lambda_1}\right]$$
绝热拉伸温升(从 λ₁ 拉伸到 λ₂)。注意 L₀ 不出现 — 与原始长度无关!
💡 数量级:典型橡胶从 λ=1.5 拉到 λ=1.9,ΔT ≈ 0.12 K — 微小但可测,Joule 用望远镜读取膨胀。
重要应用:橡胶热机 — 让橡胶在不同温度下伸缩做功(YouTube 上有演示)。

3.6 应力 vs 温度(在恒定 λ 下)

实验:把橡胶拉到 λ=3.5 后保持,测应力随 T 变化。
经典理论预测:τ ∝ T(线性!)
实验:的确线性,截距小但非零
截距非零 → 还有一点内能贡献(不完美熵性)。
在 Tg 附近线性关系破坏 — 进入玻璃态后弹性来自内能,不再是熵。

3.7 选择题

1G 与什么成正比?
B:G = ρRT/Mc。交联点越密 (Mc 小) 越硬。
2λ=2 时 τ = G(λ-1/λ²) =?
C:2 - 1/4 = 1.75。
3为什么橡胶 ν=0.5(不可压缩)?
D:链段在交联点间自由移动 — 局部液体行为,体模量很大但剪切模量很小。

🎈 第4章 气球问题 / Balloon Problem

Week 5 Q3 经典原题。 用经典橡胶弹性证明:吹气球开始难,过了某一点突然变容易。临界半径比 r/r₀ ≈ 1.38 是普适规律

4.1 解题路线图(5 步)

  1. 用不可压缩条件得拉伸比:λ₁=λ₂=r/r₀, λ₃=(r₀/r)²
  2. 代入变形熵 ΔS → 弹性自由能 A = -TΔS(每单位体积)
  3. 乘以橡胶皮总体积 4πr₀²t₀ → A_total(r)
  4. 外力做功 dW = P·4πr²·dr,等于弹性能变化 dA_total
  5. 解得 P(r),求 dP/dr=0 → r/r₀ = 7^(1/6) ≈ 1.38

4.2 关键公式

$$\lambda_1 = \lambda_2 = \frac{r}{r_0}, \quad \lambda_3 = \left(\frac{r_0}{r}\right)^2$$
双轴拉伸 + 不可压缩 → 厚度变化
$$A_{sphere} = 4\pi r_0^2 t_0 \cdot \frac{nk_BT}{2}\left[2\left(\frac{r}{r_0}\right)^2 + \left(\frac{r_0}{r}\right)^4 - 3\right]$$
球形气球的弹性自由能
$$P = \frac{2nk_BT \cdot t_0}{r_0}\left[\frac{r_0}{r} - \left(\frac{r_0}{r}\right)^7\right]$$
气球压力(最终结果)

4.3 求最大值 — 数学小技巧

令 α = r/r₀,则 P ∝ (1/α - 1/α⁷)
dP/dα = -1/α² + 7/α⁸ = 0
→ 7 = α⁶
α = 7^(1/6) ≈ 1.383

4.4 图解:压力 vs 半径

r/r₀ P 峰值 r/r₀ = 1.38 (开始吹) 变软(容易吹)

4.5 完整解答 / Full Solution

原题 · Week 5 Q3

气球问题完整推导

(a) 球壳从 r₀ → r 的双轴拉伸:
表面 λ₁=λ₂ = r/r₀。不可压缩 λ₁λ₂λ₃=1 → λ₃ = (r₀/r)²
(b) 弹性自由能(每单位体积):
$A_v = -T \cdot \Delta S = \dfrac{nk_BT}{2}(\lambda_1^2 + \lambda_2^2 + \lambda_3^2 - 3)$
$= \dfrac{nk_BT}{2}\left[2\left(\dfrac{r}{r_0}\right)^2 + \left(\dfrac{r_0}{r}\right)^4 - 3\right]$
球壳总体积 ≈ 4πr₀²t₀(薄壳近似)
$A_{sphere} = 4\pi r_0^2 t_0 \cdot A_v$
(c) 外力做功 = 弹性能变化:
$P \cdot 4\pi r^2 \cdot dr = dA_{sphere}$
$P = \dfrac{1}{4\pi r^2}\dfrac{dA}{dr}$
微分 A_v 对 r 得:$\dfrac{dA_v}{dr} = nk_BT\left[\dfrac{2r}{r_0^2} - \dfrac{2r_0^4}{r^5}\right]$
代入并简化:
$$\boxed{P = \dfrac{2nk_BT \cdot t_0}{r_0}\left[\dfrac{r_0}{r} - \left(\dfrac{r_0}{r}\right)^7\right]}$$
(d) 令 α = r/r₀,P ∝ 1/α - 1/α⁷
$dP/d\alpha = -\dfrac{1}{\alpha^2} + \dfrac{7}{\alpha^8} = 0$ → α⁶ = 7 → α = 7^(1/6) ≈ 1.38
(e) 加热密封气球 → 内部压力 ∝ T,但同时橡胶模量 G ∝ T,两个 T 完全抵消 → 气球大小不变
启示:理想气体压力的本质也是熵 — 这是关于熵的优雅观察。
💡 1.38 是普适规律,独立于温度、橡胶种类、气球大小 — 只要橡胶满足 Gaussian 网络模型。

4.6 选择题

1为什么气球难度峰值在 r/r₀=1.38?
D:1.38 = 7^(1/6) 是普适常数。
2把已吹起的密封气球放到太阳下加热
B:气压 ∝ T,橡胶 G ∝ T,两者完全抵消。

📕 第5章 修正理论 / Corrections to Classical Theory

本章涉及 Week 6 lecture 7 + Week 6 Q2。 经典理论在某些情况下与实验偏差。本章讨论修正方法。

5.1 三个偏差区域

λ τ 经典 G(λ-1/λ²) 实验 小应变 中等 大应变 ↓ 网络缺陷 ↑ Mooney 行为 ↑ 应变诱导结晶
实验 vs 经典理论:在三个区域有不同偏差

5.2 关键修正

Network Defects
网络缺陷
点击翻转
两类缺陷:
悬挂末端 (dangling ends):每条原始链有 2 个
环 (loops):自我交联,对弹性无贡献
结果:实际有效链数 < n_c。
Flory Dangling End Correction
Flory 末端修正
点击翻转
$$G_{eff} = \frac{\rho RT}{M_c}\left(1 - \frac{2M_c}{M_p}\right)$$
Mp = 原链分子量;2 = 每条链 2 个末端。
Mp → ∞ 时回到经典。
Mooney Plot
Mooney-Rivlin 图
点击翻转
$\dfrac{\tau}{\lambda - 1/\lambda^2}$ vs $1/\lambda$
经典理论预测水平线 = G。
实测:$2C_1 + 2C_2/\lambda$ — 有斜率!揭示偏差。
Strain-Induced Crystallization
应变诱导结晶
点击翻转
大 λ → 链段对齐 → 结晶 → 模量飙升。
天然橡胶 此效应明显(高拉伸强度的来源)。
合成橡胶(无规共聚)通常没有。
Beyond Gaussian
超出高斯近似
点击翻转
高斯只在拉伸 << 完全伸展时有效。
当链接近完全拉直时,需要拉伸真实化学键 → 模量陡升。
这是大 λ 硬化的另一来源。
Flory-Rehner Theory
Flory-Rehner 溶胀
点击翻转
交联橡胶在溶剂中平衡溶胀:
弹性回缩力 ⇌ 混合熵
$\ln\phi_1 + \phi_2 + \chi\phi_2^2 + nV_0(\phi_2^{1/3} - \phi_2/2) = 0$

5.3 关键公式

$$G_{eff} = \frac{\rho RT}{M_c}\left(1 - \frac{2M_c}{M_p}\right)$$
Flory 末端修正
$$\frac{\tau}{\lambda - 1/\lambda^2} = 2C_1 + \frac{2C_2}{\lambda}$$
Mooney-Rivlin 经验式(C₁+C₂ ≈ G/2)
$$\ln(1-\phi_2) + \phi_2 + \chi\phi_2^2 + n V_0 \left(\phi_2^{1/3} - \frac{\phi_2}{2}\right) = 0$$
Flory-Rehner 溶胀平衡条件 (φ₂=橡胶, n=交联数密度, V₀=溶剂分子体积)
💡 记忆点:
• $M_c$ = 交联间分子量;$M_p$ = 原始链分子量。
• Mooney 图水平 → 经典理论好;有斜率 → 偏差。
• 天然橡胶大 λ 硬化 = 应变诱导结晶
• 良溶剂 χ < 0.5;不良溶剂 χ > 0.5;θ-溶剂 χ = 0.5。

5.4 计算题

原题 · Week 6 Q2

橡胶在甲苯中的溶胀(重要!)

未变形 100×10×10 mm,Mc=7,781 g/mol。20°C 甲苯(C₆H₅CH₃, ρ=0.86)中溶胀至 176×17.6×17.6 mm。橡胶 ρ=0.92 g/cm³。
(a) 求 χ;(b) 想要线性变化 ≤5%,需要多大 Mc?(c) 评论 (b)。

(a) 体积比 = 1.76³ = 5.45
φ₂ = 1/5.45 = 0.183,   φ₁ = 1 - φ₂ = 0.817
甲苯 V₀ = M/ρ/N_A = (92/0.86)×10⁻⁶/(6.02×10²³) = 1.78×10⁻²⁸ m³
n = ρN_A/Mc = (920×6.02×10²³)/7.781 = 7.12×10²⁵ m⁻³
nV₀ = 0.01267
φ₂^(1/3) - φ₂/2 = 0.568 - 0.092 = 0.476
Flory-Rehner: ln(0.817) + 0.183 + 0.0335χ + 0.01267×0.476 = 0
-0.0193 + 0.0335χ + 0.00603 = 0
χ = 0.0133/0.0335 ≈ 0.39  (Live Session 用稍不同的算法得 0.355;教科书值 ≈0.39)
(b) 5% 线性 → 体积比 1.05³ = 1.158
φ₂ = 1/1.158 = 0.864,   φ₁ = 0.136
ln(0.136) + 0.864 + 0.39×0.747 + nV₀(0.953-0.432) = 0
-1.995 + 0.864 + 0.291 + 0.521·nV₀ = 0
nV₀ = 0.840/0.521 = 1.612
n = 1.612/(1.78×10⁻²⁸) = 9.06×10²⁷ m⁻³
Mc = ρN_A/n = (920×6.02×10²³)/(9.06×10²⁷) = 61 g/mol
(c) Mc = 61 g/mol 极不现实!相当于每个单体都要交联(橡胶单体 Mr~70)。
替代方案:(i) 选 χ > 0.5 的不良溶剂(如氟橡胶配甲苯);(ii) 阻隔涂层;(iii) 完全不同的材料(如金属);(iv) 填料增强网络。
💡 应用:橡胶 O-ring、防化手套(latex、nitrile、Viton)、燃料管设计。

5.5 选择题

1用更长的初始链 Mp,相同 Mc → G 怎么变?
A:(1-2Mc/Mp) 趋近 1 → G 接近 ρRT/Mc。物理:长链中悬挂末端比例小。
2天然橡胶在大 λ 显著硬化主要因为?
D:天然橡胶链段对齐时形成晶体,相当于增加交联点,模量陡升。这是天然橡胶高拉伸强度的来源。

📒 第6章 粘弹性 / Viscoelasticity

本章涉及 Week 7 lecture 8。 没有交联的聚合物熔体 — 短时像橡胶,长时像液体。如何统一描述?答:粘弹性。

6.1 关键概念

Stress Relaxation Modulus G(t)
应力弛豫模量
点击翻转
在 t=0 突然施加固定应变 ε₀,测应力 σ(t):
$G(t) = \sigma(t)/\varepsilon_0$
橡胶 → 常数;液体 → 衰减到零。
Creep Compliance J(t)
蠕变柔量
点击翻转
在 t=0 突然施加固定应力 σ₀,测应变 ε(t):
$J(t) = \varepsilon(t)/\sigma_0$
固体 → 常数;液体 → 线性增加 (η)。
Storage & Loss Modulus
储能 & 损耗模量
点击翻转
震荡应变 ε(t) = ε₀ cos(ωt) → 应力响应分两部分:
• G'(ω) = 与 ε 同相 → 储能(弹性)
• G''(ω) = 90° 相差 → 损耗(粘性)
Complex Modulus G*
复模量
点击翻转
$G^*(\omega) = G'(\omega) + i G''(\omega)$
$G^*(\omega) = i\omega \int_0^\infty e^{-i\omega t} G(t) dt$
是 G(t) 的傅立叶变换。
Loss Tangent
损耗角正切
点击翻转
$\tan\delta = G''/G'$
tanδ ≈ 0:纯弹性(橡胶)
tanδ → ∞:纯粘性(液体)
tanδ ≈ 1:玻璃化转变附近峰值
η ~ G₀τ
粘度估算
点击翻转
量纲关系:粘度 = 模量 × 弛豫时间
$\eta \sim G_0 \tau$  
对聚合物熔体:$\eta \sim G_p \tau_T$

6.2 弛豫模量曲线(关键图)

log(t) log G(t) N₁ (短链) N₂ (长链) Plateau Modulus G_p 玻璃态 橡胶/缠结 液体流动 τ_T(1) τ_T(2)
G(t) 曲线:玻璃态 → 平台(缠结)→ 终态液体
💡 记忆点:
Plateau modulus G_p 与分子量无关,但 plateau 长度(τ_T)强烈依赖分子量!
• G' = 储能(弹性);G'' = 损耗(粘性)
• G* 是 G(t) 的傅立叶变换(长时↔低频,短时↔高频)
• 实际测量:rheometer + 频率扫描

🐍 第7章 Reptation 蛇形蠕动 / Reptation Theory

本章涉及 Week 7 lecture 8 + Week 7 Q2。 de Gennes (1971) 与 Doi-Edwards (1986) 的革命性工作 — 缠结链像蛇一样在"管"中爬行。de Gennes 因此获 1991 诺贝尔奖。

7.1 缠结的物理图像

没有交联的高分子量聚合物熔体居然像橡胶?为什么?
答:链与链不能相互穿过 → 形成拓扑约束 → 等效于"临时交联"。
Edwards 1965 把这些约束建模为"管" (tube model)。

7.2 关键概念

Tube Model
管模型
点击翻转
每条链周围有"虚拟管",由邻链构成。
链不能穿过管壁,只能沿管轴爬行
管直径 ~ Me 的 RMS 端距。
Slip-Rings vs Tube
滑环 vs 管
点击翻转
两种缠结建模方法:
• 滑环:把链通过环连接,环可滑但不可解
• 管:把约束抽象为"管腔"
都给出相似结果。
Plateau Modulus G_p
平台模量
点击翻转
$G_p = \rho RT/M_e$
Me = 缠结分子量(材料常数)。
注意:与单链分子量无关!
对 PB: Gp ≈ 1.4 MPa, Me ≈ 1700
Reptation
蛇形蠕动
点击翻转
链像蛇一样沿管轴爬行(曲线扩散)。
当链完全爬出原管时,原始构象记忆消失 → 弛豫完成。
"忘记原管的时间" = τ_T
Terminal Time τ_T
终态弛豫时间
点击翻转
推导:
μ_tube = μ_segment/N(N 段需协同移动)
D_tube = k_BT μ_segment/N
L = Na(管长 ~ 链伸展)
τ_T ~ L²/D_tube ~
Self-Diffusion D
自扩散系数
点击翻转
链中心扩散一个 RMS 末端距(实空间):
D_self ~ ⟨R²⟩/τ_T ~ Na²/N³ = a²/N²
$\boxed{D_{self} \sim 1/N^2}$
实验上 D ~ N^(-2.3),接近 reptation 预测。
Theory vs Experiment
理论 vs 实验
点击翻转
理论:τ_T ~ N³, η ~ N³
实验:τ_T ~ N^3.4, η ~ N^3.4
偏差来源:管伸缩涨落 (CLF)、约束释放 (CR)、星状链不能 reptate(内部约束)等。

7.3 关键公式

$$G_p = \frac{\rho RT}{M_e}$$
平台模量(缠结点起"临时交联"作用)
$$\tau_T \sim N^3 \text{ (reptation)}, \quad \tau_T \sim N^{3.4} \text{ (experiment)}$$
终态弛豫时间标度
$$\eta \sim \tau_T \cdot G_p$$
零剪切粘度
$$\tau_T \approx \frac{\langle R^2 \rangle}{D} = \frac{Nb^2}{D}, \quad D_{self} \sim \frac{1}{N^2}$$
由扩散系数估算 τ_T;自扩散系数标度

7.4 Brownian motion 基础(Reptation 推导前提)

$$F_{drag} = \xi v$$
阻力:粘性流体中的"摩擦"。ξ = 阻力系数
$$D = \frac{k_B T}{\xi} \quad\text{(Einstein relation)}$$
扩散系数 - 阻力关系(同源:来自分子碰撞)
$$\xi = 6\pi \eta a \quad\text{(Stokes' Law)}$$
球形粒子(半径 a)在粘度 η 流体中的阻力
$$\langle x^2 \rangle = 6 D t \quad\text{(3D)}, \quad \langle x^2 \rangle = 2 D t \quad\text{(1D)}$$
扩散粒子的均方位移
$$\mu = 1/\xi \quad\text{(mobility)}$$
迁移率(阻力的倒数)

7.5 Reptation 完整推导 (Doi-Edwards)

逻辑链:
1. 单段迁移率 μ_seg = 1/ξ_seg(来自 Stokes-Einstein)
2. N 段同步运动 → tube 总迁移率 μ_tube = μ_seg / N
3. tube 扩散系数 D_tube = k_BT · μ_tube = k_BT μ_seg / N
4. tube 长度 L ~ N · a(链伸展度)
5. 沿 tube 扩散 L 距离的时间
    τ_T ~ L² / D_tube = (Na)² · N / (k_BT μ_seg) ~

由 τ_T 推 D_self(实空间扩散系数):
链中心扩散一个真实空间末端距 ⟨R²⟩^(1/2) = b√N 的时间也是 τ_T
D_self ~ ⟨R²⟩ / τ_T ~ Nb² / N³ = b²/N² → D_self ∝ 1/N²
实验值 D_self ~ N^(-2.3)(接近,但有偏差)。

7.5 计算题

原题 · Week 7 Q2

聚丁二烯 reptation 估算

PB N=1570, T=23°C, D=3.4×10⁻¹³ cm²/s, G_p=1.4 MPa, b=0.8 nm。
(a) 估 τ_T;(b) 估 η;(c) 用 reptation 估算 τ_T=1 day 的 N;(d) 用实验 N^3.4 重做。

(a) ⟨R²⟩ = Nb² = 1570 × (0.8×10⁻⁹)² = 1.0×10⁻¹⁵ m²
D = 3.4×10⁻¹³ cm²/s = 3.4×10⁻¹⁷ m²/s
τ_T ≈ R²/D = 10⁻¹⁵/3.4×10⁻¹⁷ ≈ 30 s
(b) η ~ τ_T · Gp = 30 × 1.4×10⁶ = 4.2×10⁷ Pa·s
(对比:水 10⁻³ Pa·s, 蜂蜜 10 Pa·s — 这种聚合物熔体粘度极高)
(c) 1 day = 86,400 s. 比 30 s 大 2880 倍。
τ ∝ N³ → N_new/1570 = 2880^(1/3) = 14.22
N_new = 22,300
(d) τ ∝ N^3.4 → N_new/1570 = 2880^(1/3.4) = 10.95
N_new = 17,200
实验标度 N^3.4 比理论 N³ 让 τ 对 N 更敏感 — 同样的 τ 增长只需要更小的 N 增量。
💡 经典方法:(1) 用 R²/D 估 τ_T;(2) 乘 Gp 估 η;(3) 用标度律推外推。

7.6 选择题

1聚合物熔体分子量翻倍,零剪切粘度(reptation 理论):
B:η ~ N³,N×2 → η×8。实验 N^3.4 → ×10.6。
2下列哪个量独立于聚合物分子量?
C:G_p = ρRT/Me 只依赖 Me(材料常数),与单链 N 无关。
3星状聚合物的粘度比同分子量线性聚合物:
B:星状链有中心节点,不能像蛇一样沿管爬行。需要"轮廓涨落 (CLF)"才能松弛 → 粘度更高,扩散更慢。

📓 第8章 时温叠加 & 玻璃化 / WLF & Glass Transition

本章涉及 Week 7 lectures 9 + Week 7 Q1。 同一聚合物在不同 T 下行为完全不同 — 高温像液体,低温像玻璃。如何用一个公式描述?

8.1 关键概念

Time-Temp Superposition
时温叠加原理
点击翻转
高温 ↔ 长时间。所有 T 的曲线可平移叠加成一条主曲线
$G(t,T) = G(a_T t, T_{ref})$
a_T = 移动因子。
WLF Equation
WLF 方程
点击翻转
$$\log_{10} a_T = \frac{-C_1(T-T_{ref})}{C_2 + (T-T_{ref})}$$
"通用"常数(取 Tref=Tg):
C₁≈17.4, C₂≈51.6 K
Vogel-Fulcher
Vogel-Fulcher 方程
点击翻转
$$\eta = \eta_0 \exp\left(\frac{B}{T-T_0}\right)$$
T₀ = Vogel 温度(理想玻璃)
典型 T₀ ≈ Tg - 50 K
WLF ↔ V-F
两公式等价
点击翻转
完全数学等价!
$C_2 = T_{ref} - T_0$
$B = 2.303 \cdot C_1 \cdot C_2$
WLF 是 V-F 取以 T_ref 为参考点的 log 形式。
Glass Transition Tg
玻璃化温度
点击翻转
弛豫时间 τ > 实验时间尺度时观察到。
不是真正相变 — 动力学冻结
实验标志:DSC 中热容 Cp 出现突变。
Kauzmann Temperature
Kauzmann 温度
点击翻转
外推到玻璃熵 = 晶体熵的温度。
实验上 T_K ≈ T₀(Vogel 温度)。
实验冷却到此前必先玻璃化(Kauzmann 悖论)。
3K Rule of Thumb
3K 经验法则
点击翻转
Tg 附近,温度升 3 K → 弛豫时间约 ÷10。
来源:用通用 WLF 常数代入 → log a_T = -1。
这是聚合物加工 / 储存温度选择的依据。
Tg Depends on Cooling Rate
Tg 依赖冷却速率
点击翻转
慢冷 → 较低 Tg;快冷 → 较高 Tg。
差距 ~3 K per 一个数量级速率变化。
这是 Tg 不是真正相变的关键证据。

8.2 关键公式

$$\log_{10} a_T = \frac{-C_1(T - T_{ref})}{C_2 + (T - T_{ref})}$$
WLF — 移动因子
$$\eta(T) = \eta_0 \exp\left(\frac{B}{T - T_0}\right)$$
Vogel-Fulcher — 在 T → T₀ 时发散
$$C_2 = T_{ref} - T_0, \qquad B = 2.303 \cdot C_1 \cdot C_2$$
WLF ↔ V-F 转换

8.3 图解:弛豫时间 vs 1/T

1/T log τ Arrhenius Vogel-Fulcher 1/T₀ τ → ∞ 1/Tg τ ~ 10² s 液体 软固体 玻璃
💡 核心记忆点:
WLF ↔ Vogel-Fulcher 数学完全等价(Week 7 Q1a 的标准证明题)
• T₀ ≈ Tg - 50 K 典型
• Tg 取决于冷却速率(10°C/min 是标准)
• 3K → 10× 弛豫时间(数量级估算)
• 玻璃的弛豫时间可超过宇宙年龄!

8.4 通用 WLF 常数(必背!)

$$C_1 \approx 17.4, \quad C_2 \approx 51.6\text{ K} \quad\text{(when } T_{ref} = T_g\text{)}$$
对大多数聚合物近似通用。可推出 T₀ ≈ Tg - 51.6 K(比 Tg 低约 50 K)。
用法:知道 Tg 就能估算 a_T(T) !例:PS Tg=378K,T=400K → ΔT=22K → log a_T = -17.4×22/(51.6+22) = -5.2 → a_T ≈ 6×10⁻⁶(弛豫快约 17万倍)。

8.5 Kauzmann 温度 & "理想玻璃"

Kauzmann 1948:外推过冷液体的熵到低温,会发现某温度 T_K 处
S(过冷液体) = S(晶体) — 这违反热力学第三定律!
悖论的解决:在到达 T_K 之前,系统必先玻璃化(动力学冻结)。
实验上发现 T_K ≈ T₀(Vogel 温度)。
这是把热力学(Kauzmann)与动力学(Vogel-Fulcher)联系起来的深刻结果。

8.6 Cp 突变 vs 一阶相变 (重要区分)

玻璃化(动力学) T V Tg 液体 玻璃 一阶相变(晶体熔化) T Tm
玻璃化:体积有斜率改变(连续);熔化:体积不连续跳变

8.7 协同重排区域 (Cooperatively Rearranging Region, CRR)

Adam-Gibbs 1965:低温下分子运动需要协同 — 多个分子必须同时移动才能有构象改变。
"协同重排区域" 大小随 T ↘ 而 ↗。
在 Tg 附近 CRR ≈ 5-10 nm。
最近:薄膜中表面分子更易动 → 薄膜 Tg 比体相低(Keddie/Jones 1994 实验观测)。
最新(2020):用机器学习从静态结构预测动态非均匀性(DeepMind)。

8.8 计算题

原题 · Week 7 Q1(重要!)

PS WLF 分析(含等价证明)

PS Tg=105°C, C₁=13.4, C₂=43.8 K(Tref=Tg, 10°C/min 冷却率)
(a) 证 WLF ↔ Vogel-Fulcher 等价;(b) 求 T₀;(c) 求 τ(T₀+10)/τ(Tg);(d) 验证 3K 法则。

(a) 等价证明
设 η ∝ exp[B/(T-T₀)]
$\log_{10}\dfrac{\tau(T)}{\tau(T_{ref})} = \dfrac{B}{\ln 10}\left[\dfrac{1}{T-T_0} - \dfrac{1}{T_{ref}-T_0}\right]$
= (B/ln10) × [(T_ref-T₀) - (T-T₀)] / [(T-T₀)(T_ref-T₀)]
= -(B/ln10)(T-T_ref) / [(T-T₀)(T_ref-T₀)]
令 C₂ = T_ref - T₀,则 (T-T₀) = (T-T_ref) + C₂
= -(B/ln10)(T-T_ref) / [C₂((T-T_ref)+C₂)]
令 C₁ = B/(ln10·C₂),即 B = 2.303·C₁·C₂
→ log a_T = -C₁(T-T_ref) / [(T-T_ref)+C₂] ✓ WLF 形式
(b) T₀ = T_ref - C₂ = 378 - 43.8 = 334.2 K = 61°C
物理意义:Vogel 温度,在此弛豫时间发散。Tg - T₀ = 44 K,符合典型 ~50 K。
(c) B = 2.303×13.4×43.8 = 1352 K
$\dfrac{\tau(T_0+10)}{\tau(Tg)} = \exp\left[\dfrac{B}{10} - \dfrac{B}{Tg-T_0}\right] = \exp\left[\dfrac{B}{10} - \dfrac{B}{43.8}\right]$
= exp[135.2 - 30.87] = exp(104.3) ≈ 2 × 10⁴⁵
巨大数字!意味着在 T₀+10 (= 71°C) 时弛豫极慢。
评论:宇宙年龄约 4×10¹⁷ s,此弛豫时间远超宇宙年龄!这是玻璃与液体的本质区别 — 同样的无序结构,但玻璃完全冻结。
(d) 3K 法则验证
设 a_T = 10 → log a_T = 1
1 = -C₁·ΔT / (C₂+ΔT)
C₂ + ΔT = -C₁·ΔT
ΔT(1 + C₁) = -C₂
ΔT = -C₂/(1+C₁) = -43.8/14.4 = -3.04 K
✓ 3K 升温 → 弛豫快 10 倍。
💡 这道题综合检验:(1) 数学变换;(2) 物理理解;(3) 数值估算 — 必须熟练掌握。

8.9 选择题

1Tg 附近升 3K,弛豫时间
C:聚合物物理学家"3K → 10×"经验法则。
2Vogel 温度 T₀ 的物理意义
D:物理上不可达 — 实验冷却到达 Tg 时已经动力学冻结。
3关于 Tg 哪个正确?
B:动力学冻结,慢冷低 Tg、快冷高 Tg。这是 Tg 不是相变的关键证据。

📐 公式速查表 / Formula Cheat Sheet

考前必背!按主题分类的所有核心公式。

1️⃣ 链统计

公式说明
$\langle R^2 \rangle = Na^2 = N_{eff}b^2$FJC / 真实链均方末端距
$\langle R^2 \rangle^{1/2} = a\sqrt{N}$RMS 末端距
$C_\infty = b^2/a^2$特征比
$N = M_w/M_r$,   $L_{contour} = Na$聚合度,轮廓长度
$P(R) \propto \exp[-3R^2/(2Na^2)]$末端距高斯分布
$\langle R^2 \rangle^{1/2} \sim aN^\nu$, ν≈0.6SAW 修正(良溶剂)

2️⃣ 熵 & 弹簧

$S(R) = -3k_BR^2/(2Na^2) + S_0$构象熵
$A = U - TS$Helmholtz 自由能
$F = (dA/dl)_{T,V} = -T(dS/dl)$ (橡胶)力 — 熵起源
$k_{spring} = 3k_BT/(Nb^2)$单链弹簧常数

3️⃣ 经典橡胶

$\Delta S/V = -nk_B(\lambda_x^2+\lambda_y^2+\lambda_z^2-3)/2$变形熵
$\lambda_x \lambda_y \lambda_z = 1$不可压缩
$G = nk_BT = \rho RT/M_x$剪切模量
$E = 3G$杨氏模量(不可压缩)
$\tau = G(\lambda - 1/\lambda^2)$应力-拉伸(非线性)
$\tau = nk_BT[(1+e) - 1/(1+e)^2]$用工程应变 e
$d = (G/k_BT)^{-1/3} = n^{-1/3}$交联点平均间距(实空间)

4️⃣ 气球

$\lambda_1=\lambda_2=r/r_0, \lambda_3=(r_0/r)^2$双轴 + 不可压缩
$P = \frac{2nk_BT t_0}{r_0}\left[\frac{r_0}{r} - \left(\frac{r_0}{r}\right)^7\right]$气球压力
$r/r_0 = 7^{1/6} \approx 1.38$压力峰值(普适)

5️⃣ 修正

$G_{eff} = (\rho RT/M_c)(1 - 2M_c/M_p)$Flory 末端修正
$\tau/(\lambda-1/\lambda^2) = 2C_1 + 2C_2/\lambda$Mooney-Rivlin
$\ln(1-\phi_2)+\phi_2+\chi\phi_2^2+nV_0(\phi_2^{1/3}-\phi_2/2)=0$Flory-Rehner
χ < 0.5:良溶剂;χ > 0.5:不良溶剂Flory-Huggins 参数

6️⃣ 粘弹性

$G(t) = \sigma(t)/\varepsilon_0$应力弛豫模量
$J(t) = \varepsilon(t)/\sigma_0$蠕变柔量
$G^*(\omega) = G'(\omega) + iG''(\omega)$复模量
$\tan\delta = G''/G'$损耗角正切
$\eta \sim G_0 \tau$粘度量纲关系

7️⃣ Reptation

$G_p = \rho RT/M_e$平台模量
$\eta \sim \tau_T G_p$零剪切粘度
$\tau_T \sim N^3$ (理论), $N^{3.4}$ (实验)终态时间
$\tau_T \approx Nb^2/D$由 D 估 τ_T
$D_{self} \sim 1/N^2$自扩散系数

8️⃣ WLF / Tg

$\log_{10}a_T = -C_1(T-T_{ref})/[C_2+(T-T_{ref})]$WLF
$\eta = \eta_0 \exp[B/(T-T_0)]$Vogel-Fulcher
$C_2 = T_{ref}-T_0$转换 1
$B = 2.303 C_1 C_2$转换 2
"通用" C₁ ≈ 17.4, C₂ ≈ 51.6 K (Tref=Tg)典型值
3K → ÷10 弛豫时间经验法则

📌 常用常数

k_B = 1.381 × 10⁻²³ J/K玻尔兹曼常数
R = 8.314 J/(mol·K)气体常数
N_A = 6.022 × 10²³ /mol阿伏加德罗
0°C = 273 K,   1 bar ≈ 10⁵ Pa转换
橡胶 ρ ≈ 0.9 g/cm³ = 900 kg/m³典型密度
⚠️ 单位陷阱:
• G = nk_BT 时 n 是每立方米链数
• G = ρRT/Mc 时 ρ in kg/m³,Mc in kg/mol
• 温度永远用 K(不是 °C)
• Pa = N/m² = J/m³
• cm⁻³ → m⁻³ 乘 10⁶;cm² → m² 乘 10⁻⁴

🎯 公式应用场景指南 / When to Use Which Formula

这是考前最有用的一页!看到不同的题目"信号词"就知道该用哪个公式。考试时按这个表先识别题型,再套公式。

🟦 题型 1:算单链尺寸 / Single Chain Dimensions

触发词:"求末端距 / RMS distance / 链尺寸 / size of polymer / contour length / 完全伸展"
给定:Mw, Mr, b 或 a, C∞
问什么用什么公式步骤
RMS 末端距$R = b\sqrt{N}$① N=Mw/Mr ② R = b√N
轮廓长度$L = N \cdot a$① 用真实键长 a,不是 b ② a = b/√C∞
同 L 不同聚合物同样 N(骨架原子数)新分子量 = N × Mr_新
排序刚性看 C∞大 = 刚(共轭、双链);小 = 柔(醚键)
良溶剂中链尺寸$R \sim aN^{0.6}$excluded volume → ν=0.6
熔体中链尺寸$R = b\sqrt{N}$(理想)excluded volume 被屏蔽 → ν=0.5
典型题:Week 4 Q1(PS Mw=100k 的 R 和 L)

🟩 题型 2:单链熵弹簧拉伸 / Single Chain AFM

触发词:"AFM / 单链 / single chain stretch / 拉伸力 / theta conditions / θ-溶剂"
给定:Mw 或 N, b, T
问什么用什么公式步骤
熵弹簧常数$k = \dfrac{3k_BT}{Nb^2}$① 求 N ② 代入 (注意 b 用 m, 单位 N/m)
力 vs 拉伸距离$F = k \cdot x$胡克定律(小拉伸时高斯近似有效)
能量需求$A = \frac{1}{2}kx^2$抛物线储能
典型题:Week 5 Q2(AFM 拉 PS 链 100 nm)

🟧 题型 3:橡胶剪切模量 / Rubber Shear Modulus

触发词:"shear modulus / 剪切模量 / 立方块/cube + 剪切力 / Mc / molecular weight between crosslinks / 交联点间分子量"
给定:F, A, Δx, l, T, ρ(部分)
问什么用什么公式步骤
从实验数据求 G$G = \tau/\gamma$τ=F/A,γ=Δx/l
已知 G 求 Mc$M_c = \rho RT/G$注意:ρ kg/m³, T 用 K
已知 Mc 求 G$G = \rho RT/M_c$反向
交联密度 n$n = G/(k_BT)$n 单位是 m⁻³
交联点平均距离$d = n^{-1/3}$典型 ~ 1-2 nm
典型题:Week 6 Q1(聚丁二烯立方块)

🟪 题型 4:橡胶应力-应变 / Stress-Strain

触发词:"extension ratio λ / 拉伸比 / stress at λ=X / extension X% / 压缩 X mm"
给定:G 或 Mc,λ 或 e
变形方式公式注意
大拉伸 (λ > 1.1)$\tau = G(\lambda - 1/\lambda^2)$非线性!
用工程应变 e=λ-1$\tau = nk_BT[(1+e) - (1+e)^{-2}]$同样公式不同表达
小应变(≤5%)$\sigma = E\varepsilon = 3G\varepsilon$线性近似!E=3G(不可压缩)
剪切应变(任意大小)$\tau_{shear} = G\gamma$胡克定律
压缩力$F = E A \varepsilon = 3GA\varepsilon$=3 × 同等剪切力
大 λ 硬化无简单公式归因于应变诱导结晶 + 超出高斯
典型题:Week 6 Q1d(压缩 4.46 mm)

🟫 题型 5:双轴变形(气球)/ Biaxial

触发词:"balloon / 气球 / 球形薄壳 / sphere / inflate / 双轴拉伸"
给定:r₀, t₀, n(或 Mc)
问什么公式
拉伸比$\lambda_1=\lambda_2 = r/r_0$,   $\lambda_3 = (r_0/r)^2$(不可压缩)
压力$P = \dfrac{2nk_BT t_0}{r_0}\left[\dfrac{r_0}{r} - \left(\dfrac{r_0}{r}\right)^7\right]$
峰值位置$r/r_0 = 7^{1/6} \approx 1.38$(普适!)
加热行为P ∝ T 与 G ∝ T 抵消 → 大小不变
典型题:Week 5 Q3

🟨 题型 6:溶胀 / Swelling

触发词:"swelling / 溶胀 / Flory-Rehner / Flory-Huggins / χ / 在溶剂中浸泡"
给定:变形前后尺寸、ρ_polymer, ρ_solvent, M_solvent, Mc
问什么步骤
聚合物体积分数$\phi_2 = V_0/V = 1/(\text{体积比})$
溶剂分子体积$V_0 = M_{solv}/(\rho_{solv} \cdot N_A)$
交联数密度$n = \rho \cdot N_A / M_c$  (单位 m⁻³)
Flory-Huggins χ用 Flory-Rehner 反解 χ
所需 Mc设给定 φ₂,反解 Mc
判断良溶剂χ < 0.5 良;χ > 0.5 不良;= 0.5 θ
$$\ln(1-\phi_2) + \phi_2 + \chi\phi_2^2 + nV_0\left(\phi_2^{1/3} - \frac{\phi_2}{2}\right) = 0$$
Flory-Rehner 平衡方程(核心)
典型题:Week 6 Q2

🟥 题型 7:粘弹性 / Viscoelasticity

触发词:"plateau modulus / Gp / Me / 缠结 / entanglement / terminal time / τT / viscosity / 粘度 / reptation / 蛇形 / self-diffusion"
给定:N, b, D, Gp(部分)
问什么公式
平台模量$G_p = \rho RT/M_e$
由 D 估 τT$\tau_T \approx \langle R^2 \rangle/D = Nb^2/D$
粘度$\eta \sim \tau_T \cdot G_p$
τT 标度(理论)$\tau_T \sim N^3$
τT 标度(实验)$\tau_T \sim N^{3.4}$
D 标度$D \sim 1/N^2$
由旧 N₁ 推新 N₂$N_2/N_1 = (\tau_2/\tau_1)^{1/3}$ 或 $(\tau_2/\tau_1)^{1/3.4}$
Brownian motion$\langle x^2 \rangle = 6Dt$, $D = k_BT/\xi$, $\xi = 6\pi\eta a$
典型题:Week 7 Q2

🟦 题型 8:WLF / 玻璃化

触发词:"WLF / shift factor / aT / Vogel-Fulcher / Tref / 玻璃化温度 Tg / glass transition / cooling rate"
给定:Tg, C₁, C₂(或者要外推它们)
问什么公式
WLF$\log_{10} a_T = -C_1(T-T_{ref})/[C_2 + (T-T_{ref})]$
Vogel-Fulcher$\eta = \eta_0 \exp[B/(T-T_0)]$
转换:求 T₀$T_0 = T_{ref} - C_2$
转换:求 B$B = 2.303 \cdot C_1 \cdot C_2$
τ(T) 比 τ(Tref)$\dfrac{\tau(T)}{\tau(T_{ref})} = a_T = 10^{-C_1\Delta T/(C_2 + \Delta T)}$
不同 T 的 τ 比$\dfrac{\tau_1}{\tau_2} = \exp\left[\dfrac{B}{T_1-T_0} - \dfrac{B}{T_2-T_0}\right]$
3K 升温 → ÷10 弛豫近 Tg 经验法则
通用常数C₁≈17.4, C₂≈51.6 K(Tref=Tg)→ T₀ ≈ Tg - 50 K
典型题:Week 7 Q1(PS WLF 等价证明)

🌳 解题决策树 / Decision Tree

题目里有什么关键词? 先识别题型 分子量 / N / 链尺寸 → 题型 1 链统计 AFM / 单链拉伸 → 题型 2 熵弹簧 立方体/立方块橡胶 → 题型 3,4 经典橡胶 气球 / 双轴 / 球壳 → 题型 5 气球问题 溶剂 / 溶胀 / χ → 题型 6 Flory-Rehner 缠结 / Gp / Me / 粘度 → 题型 7 Reptation Tg / WLF / aT / Vogel → 题型 8 玻璃化 Mooney 图 / 偏离经典 → 修正理论 如果是橡胶变形:剪切还是拉伸? 关键分叉 剪切 / shear τ = G·γ — 胡克定律 拉伸/压缩 / extension τ = G(λ-1/λ²)(大)或 3Gε(小) 单位检查 → 计算 → 验证数量级 ρ kg/m³, T 用 K, 用 R 时 Mc 用 kg/mol
考试时按这个流程图先分类,再选公式

🔑 考试快速识别表 / Quick Cheat Sheet

题目里看见…立即想到…
"Mw / 分子量 / Mr / N"题型 1 — 用 N=Mw/Mr,R=b√N
"AFM / 拉一根链 / theta"题型 2 — k=3k_BT/(Nb²)
"立方体 + 剪切力 / Mc"题型 3 — G=ρRT/Mc,n=G/k_BT
"λ / extension ratio / 5%"题型 4 — τ=G(λ-1/λ²) 或 σ=3Gε
"球形 / balloon / inflate"题型 5 — 1.38, P 公式
"swell / 溶剂 / 体积变化 / χ"题型 6 — Flory-Rehner
"Gp / Me / N^3 / 缠结"题型 7 — Gp=ρRT/Me, τ=R²/D
"WLF / Tg / Vogel"题型 8 — WLF/V-F, T₀=Tref-C₂
"加热 → 体积/形状变化"记住 G ∝ T 反常行为;气球加热不变
"应变诱导 / 大 λ 硬化"天然橡胶结晶 + 超出 Gaussian
"Mooney 图 / 偏离经典"τ/(λ-1/λ²) vs 1/λ,斜率给 C₂
"剪切(不是拉伸)"τ=G·γ — 胡克定律,线性

⚠️ 常见错误 / Common Pitfalls

1. 单位错误(最常见!)
• ρ 必须是 kg/m³(0.9 g/cm³ → 900 kg/m³)
• T 永远是 K(0°C = 273 K,23°C = 296 K,105°C = 378 K)
• 用 R=8.314 时 Mc 是 kg/mol;用 k_B 时 n 是 1/m³
• cm⁻³ → m⁻³ ×10⁶;cm² → m² ×10⁻⁴
2. 概念混淆:
R(链 RMS 末端距)≠ d(交联点平均距离)
R = b√N 是同一根链的两端距离。
d = n⁻¹/³ 是不同交联点之间的实空间距离。
Mc ≠ Mp!Mc = 交联点之间,Mp = 整个链。
a ≠ b!a = 真实键长,b = 统计步长。
3. 公式选择:
• 小应变 (~5%) → 线性 σ = 3Gε(更简单)
• 大应变 (>10%) → 非线性 τ = G(λ - 1/λ²)
• 任意剪切 → 胡克定律 τ = G·γ
• Flory-Rehner 用 ln(1-φ₂),不是 ln(φ₁)!
4. WLF/Vogel-Fulcher 转换方向:
• C₂ = T_ref − T₀ 不是 T₀ − T_ref
• B = 2.303·C₁·C₂(不是 C₁/C₂ 或别的)
• a_T 是移动因子,τ(T)/τ(T_ref) = a_T(不是倒数)

🎯 综合测验 / Final Quiz

15 题混合测验,覆盖全部 8 章。完成后查看分数 — 12/15 以上算掌握良好。
1PS Mw=200,000 (Mr=104, b=0.67) 的 RMS 末端距约:
✅ N=1923, R=0.67×√1923≈29.4 nm → B
2λ=4 时 τ/G =?
✅ 4 - 1/16 = 3.9375 → B
3Mc=5000 vs Mc=10000,G 比?
✅ G ∝ 1/Mc → A 大一倍。
4气球难度峰值在 r/r₀=1.38 的原因?
✅ 普适规律。→ D
5Reptation 中 D ~ ?
✅ D ~ R²/τ ~ N/N³ = 1/N² → B
6下面哪个直接显示橡胶弹性是熵性?
✅ G ∝ T 是熵起源标志。→ C
7T_ref=Tg+50, C₁=8.86, C₂=101.6 K → T₀ = ?
✅ T₀ = T_ref - C₂ = (Tg+50) - 101.6 = Tg - 51.6 → A
8Mooney-Rivlin 图截距和斜率分别给出?
✅ τ/(λ-1/λ²) vs 1/λ 的截距=2C₁, 斜率=2C₂。经典理论 C₂=0。
9下面哪个不能独立测 Mc?
✅ 声速测体模量 K,无关 Mc。→ D
10PE Mw=10⁴ τ=1 ms。Mw=10⁵ 用 reptation 估 τ ≈?
✅ τ ∝ N³,N×10 → τ×1000 = 1 s. → B
11良溶剂中网络溶胀平衡决定于?
✅ Flory-Rehner 核心。→ C
12支持"Tg 不是真正相变"的观察?
✅ 真相变与速率无关。→ B
13天然橡胶大 λ 时显著硬化的主因?
✅ 链段对齐 → 局部结晶 → 模量陡升。这是天然橡胶强度高的来源。→ C
14tan δ 在哪里取最大值?
✅ tan δ = G''/G' 在 Tg 附近达峰(粘性贡献最大)。这就是为什么用 DMA 找 Tg。→ A
15熔体中聚合物链的 ν 是?
✅ Edwards/Higgins 洞察 — 熔体中 excluded volume 被屏蔽,链回到理想 random walk。→ D
🏆 完成测验后:顶部的得分会自动统计。错题回到对应章节复习;要点是关键概念解析