⚗️ 线性粘弹性理论 LVE
交互式公式推导演示

安东帕流变仪技术支持 · 公式推导 · 动画演示 · 实测案例

📐 Boltzmann叠加 🌀 振荡测试 🌡️ 时温叠加TTS ⏱️ 松弛&蠕变
📐 Boltzmann叠加原理与时域表示
应力响应卷积积分
核心卷积方程 · 线性粘弹性基础
核心公式
σ(t)=∫−∞tG(t−t')·γ̇(t')dt'
1
线性叠加假设
各历史时刻t'的应变增量dγ=γ̇(t')dt'对当前应力贡献独立,可线性叠加。
2
单次阶跃贡献
dσ(t)=G(t−t')·γ̇(t')dt'
3
对全历史积分叠加
σ(t)=∫−∞tG(t−t')γ̇(t')dt'
4
因果性:G(τ)=0 当τ<0
材料不预知未来,积分上限取当前时刻t。
G(t)
松弛模量 Pa,随时间衰减
γ̇(t')
应变速率 s⁻¹,历史变形
γ̇(t') 应变速率
G(t−t')·γ̇ 贡献
σ(t) 累积应力
🏭 PP熔体注塑残余应力预测
螺杆产生复杂剪切历史γ̇(t'),结合MCR 702测得的Prony级数G(t),通过卷积积分预测模腔内残余应力分布,指导保压参数优化。
T=230°Cγ̇=100~1000 s⁻¹MCR 702e
⚠️测试前需预剪切+等待≥5λmax,确保无应变历史初始状态。
G(t) 与 J(t) 互倒关系
时域卷积 · Laplace域 · 数值转换
时域卷积互倒
0tG(t−t')·(t')dt'=1
Laplace域
(s)·(s)=1/s²
1
阶跃应变定义G(t)
γ(t)=γ₀H(t) → σ(t)=G(t)·γ₀
2
同一应力的柔量响应
γ(t)=∫₀ᵗJ(t−t')·σ̇(t')dt'
3
代入σ̇=G'(t')·γ₀后化简
H(t)=∫₀ᵗJ(t−t')·G'(t')dt'
4
Laplace变换+卷积定理
G̃(s)·J̃(s)=1/s²
⚠️G(t)·J(t)≠1,是卷积意义互倒,非逐点乘积!
G(t) 衰减
J(t) 增长
卷积结果=t
🔧 RheoCompass数值互转操作
Schwarzl近似:G(t)≈1/[J(t)−J²·(dJ/dt)·t/2]
路径:分析→转换→G(t)↔J(t),基于Honerkamp-Weese正则化反演,误差<5%。
𝓕
G*(ω) ↔ G(t) Fourier变换
频域时域互转 · 储能损耗模量推导
复数模量
G*(ω)=iω∫0G(t)e−iωtdt
G'(ω) 正弦变换
G'(ω)=ω∫0G(t)sin(ωt)dt
G''(ω) 余弦变换
G''(ω)=ω∫0G(t)cos(ωt)dt
1
代入振荡应变
γ=γ₀e^(iωt), γ̇=iωγ₀e^(iωt)
2
换元τ=t−t'代入卷积积分
σ=iωγ₀e^(iωt)·∫₀^∞G(τ)e^(−iωτ)dτ
3
定义G*(ω)=σ/γ
G*(ω)=iω∫₀^∞G(τ)e^(−iωτ)dτ
4
展开e^(−iωτ)分离实虚部
G'=ω∫G·sindt, G''=ω∫G·cosdt
2.0
G(t)
G·sin(ωt)→G'
G·cos(ωt)→G''
🔬 松弛数据→频率扫描预测
对G(t)做Fourier变换,将时间窗0.01~10000s扩展为等效频率10⁻⁴~100 rad/s,弥补低频测试耗时之不足。RheoCompass中:分析→转换→G(t)→G*(ω)。
时域→频域Prony级数
λ
连续松弛时间谱 H(λ)
谱表示 · 由谱计算G'/G'' · MWD表征
G(t)谱展开
G(t)=Ge+∫H(λ)e−t/λd(lnλ)
由谱计算G'
G'=Ge+∫H(λ)(λω)²/[1+(λω)²]d(lnλ)
由谱计算G''
G''=∫H(λ)λω/[1+(λω)²]d(lnλ)
1
广义Maxwell离散谱
G(t)=Σᵢ Gᵢ·exp(−t/λᵢ)
2
连续极限:Gᵢ→H(λ)d(lnλ)
G(t)=Gₑ+∫H(λ)e^(−t/λ)d(lnλ)
3
代入Fourier变换,交换积分次序
G'=Gₑ+∫H(λ)·(λω)²/[1+(λω)²]d(lnλ)
4
单Maxwell元权重函数
w'=(λω)²/[1+(λω)²], w''=λω/[1+(λω)²]
1.0
H(λ) 松弛谱
w'·H→G'贡献
w''·H→G''贡献
谱形状材料类型
单峰窄谱单分散聚合物标样
单峰宽谱宽MWD工业聚合物
双峰谱双峰HDPE/嵌段共聚物
宽平台交联网络/凝胶
🧬 MWD流变学表征
H(λ)对高分子量尾端极敏感(η₀∝M³·⁴),可检测GPC无法分辨的痕量HMW组分。RheoCompass→分析→松弛谱→Tikhonov正则化反演。
η
η₀、Je⁰ 与 λz
积分流变参数 · 分子量与分布双重表征
零剪切粘度三种等价定义
η₀=∫0G(t)dt=limω→0G''/ω=limω→0|η*|
稳态可恢复柔量
Je=(1/η₀²)∫0G(t)dt
特征松弛时间
λz=η₀·Je⁰=limω→0G'/(ω·G'')
1
η₀=G(t)曲线下面积
η₀=∫₀^∞G(t)dt (稳态流动极限)
2
频域等价
η'(ω)=G''(ω)/ω → η₀ 当ω→0
3
Je⁰弹性记忆度量
Je⁰=(1/η₀²)∫₀^∞t·G(t)dt
4
特征时间尺度
λz=η₀·Je⁰,ω<1/λz时进入终端区
G(t)
∫G dt=η₀
t·G(t)→Je
M³·⁴
η₀∝Mw³·⁴
PDI²
Je⁰∝PDI²
η₀Je
λz=特征时间
🏗️ 分子量与分布双重表征
η₀∝M³·⁴→测分子量;Je⁰∝PDI²→测分散度;λz=η₀Je⁰→判断终端流动区入口频率。
⚠️测η₀需确认G'∝ω²且G''∝ω,即ω<1/λz已进入终端流动区。
🌀 振荡测试核心公式
~
线性区振荡响应
G' G'' 定义 · 相位角 · 物理意义
应变输入
γ(t)=γ₀·sin(ωt)
线性区应力响应
σ(t)=γ₀[G'sin(ωt)+G''cos(ωt)]
模量定义
G'=(σ₀/γ₀)cosδ  |  G''=(σ₀/γ₀)sinδ
1
线性区输出同频有相位差δ
σ(t)=σ₀·sin(ωt+δ)
2
三角展开
σ=σ₀cosδ·sinωt + σ₀sinδ·cosωt
3
对比定义G' G''
G'=(σ₀/γ₀)cosδ, G''=(σ₀/γ₀)sinδ
4
物理意义
G':同相分量,弹性储存能;G'':异相分量,粘性耗散能。
δ=0°
纯弹性固体(胡克体)
δ=90°
纯粘性液体(牛顿流体)
35°
γ(t) 应变
σ(t) 应力
相位差δ
🧪 聚合物溶液凝胶化判断
1 Hz下测G'和G'':稀溶液G''≫G'(δ→90°);凝胶点G'=G''(δ=45°,Winter-Chambon判据);强凝胶G'≫G''(δ→0°)。
f=1Hzγ₀=1%(LVE内)MCR 102
G*
复数模量 |G*| 与相位角 δ
复平面表示 · tanδ损耗因子 · 能量
复数模量幅值
|G*|=σ₀/γ₀=√(G'²+G''²)
相位角与损耗因子
δ=arctan(G''/G')  |  tanδ=G''/G'
每周期能量收支
EG'γ₀²  |  EG''γ₀²
1
复数表示
G*(ω)=G'(ω)+i·G''(ω)
2
模量幅值
|G*|=√(G'²+G''²)=σ₀/γ₀
3
每周期能量
E储=π·G'·γ₀², E耗=π·G''·γ₀²
4
tanδ物理含义
每周期耗散能/储存能之比,tanδ=1为液固转变临界点。
状态tanδ范围δ
理想弹性固体0
粘弹性固体0.1~16°~45°
粘弹性液体1~∞45°~90°
理想粘性液体90°
1.00
🔊 汽车降噪材料筛选
tanδ峰值温度=阻尼最优温度。要求-10°C~60°C内tanδ>0.3,通过温度扫描优化共聚物组成,MCR 702 TwinDrive+Peltier温控精度±0.05°C。
f(ω)
频率扫描特征点分析
终端区 · 交叉点 · Cox-Merz · G_N⁰
复数粘度
η*(ω)=G*(ω)/(iω)=η'−iη''
Cox-Merz经验规则
|η*(ω)|≈η(γ̇)   当ω=γ̇ (线性缠结聚合物)
终端区斜率律
G'∝ω²  |  G''∝ω   (ω→0)
1
终端流动区(ω→0)
G'∝ω², G''∝ω, tanδ→∞(液态特征)
2
交叉点ωc:G'=G''
λc≈1/ωc,特征松弛时间的简便估算
3
橡胶平台区
G'≈G_N⁰(缠结模量),G''极小值
4
缠结分子量Me
Me=ρRT/G_N⁰,PE约1000~2000 g/mol
频率区间G'斜率G''斜率
终端流动区2(斜率2)1(斜率1)
橡胶平台区≈0(平台)极小值
玻璃化转变区急速上升峰值
3.0
G'(ω)
G''(ω)
交叉点ωc
🏭 PE熔体加工窗口确定
通过频率扫描确定ωc和G_N⁰,计算Me=ρRT/G_N⁰。Cox-Merz规则将振荡数据转换为稳态粘度曲线,指导挤出螺杆转速-温度匹配。
T=190°CPP25平行板MCR 502
J*
复数柔量 J* — 应力控制响应
J' J'' 定义 · 与G*互倒 · 低模量样品
复数柔量(注意负号)
J*=J'−iJ''=1/G*(ω)
J' J'' 与 G' G'' 转换
J'=G'/|G*|²  |  J''=G''/|G*|²
应力控制振荡应变响应
γ(t)=σ₀[J'sinωt−J''cosωt]
1
应力控制输入
σ(t)=σ₀·e^(iωt),J*=γ/σ=1/G*
2
分离实虚部
1/(G'+iG'')=(G'−iG'')/|G*|²
3
得J' J''
J'=G'/|G*|², J''=G''/|G*|²
4
符号约定(重要!)
J*=J'−iJ''(负号),G*=G'+iG''(正号),初学者常混淆。
J' 储存柔量
=G'/|G*|²,Pa⁻¹
J'' 损耗柔量
=G''/|G*|²,Pa⁻¹
45°
G* (第一象限)
J*=1/G* (第四象限)
🩺 生物凝胶低模量表征
透明质酸溶液(G*<1 Pa)用CSR应力控制模式测J*(ω),避免DSO模式惯性误差。Je⁰=limω→0J'(ω)直接给出稳态弹性柔量。
CSR模式MCR 302CP25-1锥板
⚠️G*<1 Pa时优先用应力控制(CSR),应变控制(DSO)可能因仪器惯性引入系统误差。
LVE
线性区判据 — 应变幅值扫描
5%判据 · Payne效应 · 工作点选择
工程判据(Anton Paar默认)
[G'(γ₀)−G'LVE]/G'LVE < 5%
工作点选择建议
γwork = 0.1×γLVE(保留10倍安全裕量)
1
应变扫描必须先于频率扫描
在目标频率(通常1 rad/s)下,γ₀从小到大扫描,确定G'偏离平台5%的临界点。
2
G'平台区 = LVE线性区
低应变下G'为常数;超过γ_LVE后G'下降(结构破坏)或先升后降(应变诱导结构形成)。
3
Payne效应(炭黑/纳米填料体系)
填料颗粒网络在极低应变(<0.01%)即开始破坏,G'急剧下降,G''出现峰值,需格外注意。
4
G''峰判据(更灵敏)
G''出现极大值对应的γ₀往往比G'偏离5%更早,可作为LVE边界的保守估计。
G'(γ₀)
G''(γ₀)
γ_LVE边界
材料类型典型γ_LVE注意事项
聚合物熔体1%~100%通常宽LVE区
浓聚合物溶液0.1%~10%浓度越高越窄
弱凝胶/乳液0.1%~1%结构敏感
强凝胶/颗粒体0.01%~0.1%脆性破坏
炭黑复合物<0.01%Payne效应显著
🌡️ 时温叠加原理 TTS
WLF
WLF方程 — 近Tg区域
自由体积理论 · 水平位移因子aT
WLF方程
log(aT)=−C₁(T−Tref)/[C₂+(T−Tref)]
通用常数(Tref=Tg
C₁=17.44  |  C₂=51.6 K
主曲线位移操作
G'(ω,T)=G'(aT·ω, Tref)
1
自由体积理论基础
松弛时间τ∝exp[B/f(T)],f(T)为自由体积分数,温度升高自由体积增大。
2
自由体积线性展开
f(T)=f₀+αf(T−Tg),αf≈4.8×10⁻⁴K⁻¹
3
位移因子定义
log(aT)=log[τ(T)/τ(Tref)]
4
代入自由体积表达式推导WLF
log(aT)=−C₁(T−Tref)/[C₂+(T−Tref)]
适用范围
Tg ~ Tg+100°C
失效条件
结晶/降解/非热流变简单
150°C
T=120°C
T=150°C(ref)
T=180°C
T=210°C
🔬 PS主曲线构建(Tref=150°C)
在120°C~200°C间每10°C测频率扫描,RheoCompass TTS模块自动拟合WLF参数C₁=8.86、C₂=101.6K,构建跨越15个量级频率的主曲线。
PS: C₁=8.86C₂=101.6KTref=Tg+50°C
⚠️TTS适用前提:热流变简单(thermorheologically simple),即仅时间尺度变化而谱形不变。验证:van Gurp-Palmen图(δ vs |G*|)各温度曲线重叠。
Ea
Arrhenius方程 — 远离Tg区域
流动活化能 · 高温熔体 · 1/T线性化
Arrhenius位移因子
ln(aT)=Ea/R·(1/T−1/Tref)
粘度温度依赖
η(T)=A·exp(Ea/RT)
流动活化能估算
Ea=R·Δlnη/Δ(1/T)
1
Eyring过渡态理论
流动需克服能垒Ea,流动速率∝exp(−Ea/RT),松弛时间τ∝exp(+Ea/RT)。
2
位移因子
aT=τ(T)/τ(Tref)=exp[Ea/R·(1/T−1/Tref)]
3
取对数线性化
ln(aT)=(Ea/R)·(1/T−1/Tref)
4
ln(aT) vs 1/T → 直线斜率=Ea/R
适用于T≫Tg的非极性聚合物熔体(PE、PP等)。
PE典型Ea
25~30 kJ/mol
PS典型Ea
100~120 kJ/mol
50
ln(aT) vs 1/T
斜率=Ea/R
🏗️ PP/PE挤出温度曲线优化
测3个温度(190/210/230°C)下η₀,Arrhenius拟合得Ea,外推预测任意加工温度下粘度,指导螺杆转速-温度匹配,优化能耗与产量平衡。
PP: Ea≈40 kJ/molPP25平行板MCR 302
MC
主曲线构建与垂直位移因子 bT
水平+垂直双重位移 · 密度修正 · 叠加验证方法
完整TTS位移(含垂直修正)
bT·G(ω,T)=G(aT·ω, Tref)  |  bT=ρrefTref/(ρT)
位移操作
ωshifted=aT·ω  |  G'shifted=G'(ω)/bT
1
选定参考温度Tref
通常选中间温度或Tg+50°C,aT(Tref)=1,bT(Tref)=1。
2
水平位移:频率轴平移
ω_reduced=aT·ω(低T→高频移,高T→低频移)
3
垂直位移:模量轴修正(橡胶弹性熵效应)
bT=ρrefTref/(ρT),对聚合物熔体通常<5%修正
4
叠加质量验证
各温度曲线平滑拼接;van Gurp-Palmen图(δ vs |G*|)温度无关;aT温度依赖满足WLF或Arrhenius。
验证方法判断标准适用
G' G''主曲线目视各温度曲线平滑拼接所有体系
van Gurp-Palmen图δ vs |G*|温度无关非热流变简单筛查
aT(T)拟合优度WLF/Arrhenius R²>0.999定量验证
T=120°C
T=150°C(ref)
T=180°C
T=210°C
主曲线
🔭 极低频等效预测(延伸12个量级)
实际测试最低频率约0.01 rad/s(每点需100s),通过TTS主曲线可延伸至10⁻⁶ rad/s,等效预测100天时间尺度的蠕变行为,用于食品、涂料、化妆品长期储存稳定性评估。
实测: 0.01~100 rad/s等效: 10⁻⁶~10⁶ rad/s延伸12个量级
⚠️ TTS失效案例
半结晶聚合物(PP在结晶温度附近):结晶程度随温度变化,材料结构本身改变,不满足热流变简单假设,TTS不可用。此时van Gurp-Palmen图各温度曲线明显不重叠。
⏱️ 应力松弛与蠕变
G(t)
应力松弛 — G(t)测量
阶跃应变 · Maxwell模型 · Prony级数拟合
松弛模量定义
G(t)=σ(t)/γ₀  (施加阶跃应变γ₀后)
广义Maxwell—Prony级数
G(t)=GeGᵢ·exp(−t/λᵢ)
单Maxwell元极限
G(t)=G₀·exp(−t/λ)  |  λ=η/G₀
1
Maxwell模型本构方程
dσ/dt + σ/λ = G₀·dγ/dt
2
阶跃应变:γ=γ₀H(t),dγ/dt=γ₀δ(t)
t>0时:dσ/dt + σ/λ = 0
3
常微分方程解
σ(t)=G₀γ₀·exp(−t/λ) → G(t)=G₀exp(−t/λ)
4
广义Maxwell:N个元并联
G(t)=Gₑ+Σᵢ Gᵢ·exp(−t/λᵢ)
3.0 s
阶跃应变γ(t)
松弛应力σ(t)
G(t)=σ/γ₀
🔧 RheoCompass Prony级数拟合流程
① 施加小应变γ₀(LVE内,通常1%)
② 保持5~10个最长松弛时间
③ 分析→拟合→广义Maxwell→自动确定N个{Gᵢ,λᵢ}
④ 导出Prony参数用于FEM(Abaqus/ANSYS格式)
γ₀=1%PP25平行板DSO快速施加
⚠️应变施加时间需<0.1×λ_min,否则早期数据被扭曲。建议使用DSO(直接应变振荡)模式,响应时间<10ms。
J(t)
蠕变与蠕变恢复 — J(t)测量
阶跃应力 · Voigt模型 · 稳态柔量
蠕变柔量定义
J(t)=γ(t)/σ₀  (施加阶跃应力σ₀后)
广义Voigt蠕变柔量
J(t)=J₀+ΣJₖ[1−exp(−t/τₖ)]+t/η₀
长时极限(稳态流动)
J(t)→Je+t/η₀  (t→∞,斜率=1/η₀)
1
Voigt元本构
σ=Jk⁻¹·γ+ηk·γ̇, τk=ηk·Jk
2
阶跃应力σ₀,求γ(t)
γ=σ₀·Jk·[1−exp(−t/τk)]
3
加入瞬时弹性J₀和粘性流动项
J(t)=J₀+ΣJk[1−exp(−t/τk)]+t/η₀
4
稳态极限
斜率1/η₀→η₀;截距→Je
1000
σ(t) 施加/撤除
γ(t) 蠕变
γrec 恢复
J(t)=γ/σ₀
🧴 乳液/膏体长期稳定性预测
低应力(LVE内)蠕变测试300s+恢复300s:斜率1/η₀→静置粘度;恢复率=(γmax−γ∞)/γmax×100%→弹性恢复能力。据此预测货架期沉降与分层风险。
σ₀=0.1~1 Pat=300s蠕变+300s恢复CP25-1
⚠️蠕变应力必须在LVE区内,需预先由振荡应力扫描确定σ_LVE,通常σ₀<0.1×σ_LVE。
⚙️
力学模型全景对比
Maxwell · Voigt · SLS · 广义Maxwell
模型G(t)G(t→∞)适用
MaxwellG₀exp(−t/λ)0聚合物熔体松弛
Voigt-KelvinG₀+ηδ(t)G₀固体蠕变
SLS (Zener)G∞+G₁exp(−t/λ)G∞>0软物质全程
广义MaxwellGₑ+ΣGᵢexp(−t/λᵢ)Gₑ≥0实际聚合物
SLS标准线性固体本构方程
σ+λ₁σ̇=(G+G1)γ+Gλ₁γ̇
SLS G(t) 与 G*(ω)
G(t)=G+G1e−t/λ  |  G'=G+(λω)²G1/[1+(λω)²]
Maxwell
SLS (Zener)
广义Maxwell (N=3)
1
聚合物熔体 → 广义Maxwell,Gₑ=0
G(t→∞)→0,无平衡弹性,适合未交联聚合物。
2
交联弹性体 → 广义Maxwell + Gₑ>0
Gₑ=νkT,ν为交联密度,橡胶弹性平台模量。
3
软物质/水凝胶 → SLS或广义Voigt
参数少物理意义明确,适合拟合有限测试窗口数据。
4
Prony级数N的选取原则
每十倍时间1个模式,总体跨越时间范围决定N,RheoCompass可自动优化N避免过拟合。
📋
安东帕 MCR 系列完整LVE测试流程
RheoCompass操作路径 · 常见错误排查 · 参数速查表
1
样品制备与装样
聚合物熔体:热压成型→修边成圆形→装入平行板(PP25或PP8);溶液/软物质:注射器加样,锥板(CP25-1或CP50-1)。装样后在测试温度下等待惰性气氛完全置换(N₂保护熔体)。
2
温度平衡
Peltier系统:±0.1°C精度,平衡5 min;对流烘箱(CTD):±0.5°C,平衡15~20 min。热电偶确认温度稳定后方可开始测试。
3
应变扫描 → 确定γ_LVE和γ_work
f=1 Hz,γ₀从0.01%扫到1000%(对数间隔,每decade 5~10点),以G'偏离平台值5%的γ₀为γ_LVE,取γ_work=0.1×γ_LVE用于后续所有振荡测试。
4
频率扫描 → G'(ω) G''(ω) η*(ω)
γ₀=γ_work,ω从100→0.01 rad/s(从高频到低频,避免低频长时间劣化);每十倍频5~10个点;稳态判据:连续3个采样点偏差<1%。
5
数据分析(RheoCompass路径)
拟合→广义Maxwell→Prony级数{Gᵢ,λᵢ};分析→松弛谱→H(λ)(Tikhonov正则化);TTS→WLF拟合→主曲线;自动计算→η₀、Je⁰、λz、G_N⁰;导出→FEM格式。
异常现象可能原因解决方案
G'在所有频率<G''应变超出LVE区重做应变扫描,降低γ₀
低频G'翘起样品降解/溶剂蒸发加溶剂陷阱或缩短测试时间,N₂保护
G' G''高频噪声大扭矩<0.5μNm(量程下限)升温或换小锥板(CP8)
TTS主曲线不叠加结晶/相变/多组分分相检查van Gurp-Palmen图,排查热流变简单性
η₀测不到终端区频率范围不足(λz过大)补充低频蠕变测试,间接求η₀=Je⁰/λz
松弛早期数据异常应变施加速率过慢启用DSO快速模式,或检查液压伺服响应时间
蠕变数据散乱应力过低(<σ_LVE/10)适当提高σ₀,但确保仍在LVE内
参数典型范围获取方式物理/工程意义
η₀ (Pa·s)10²~10⁷低频G''/ω;蠕变斜率∝Mw³·⁴,分子量表征
Je⁰ (Pa⁻¹)10⁻⁵~10⁻³蠕变恢复;频率扫描低频J'∝PDI²,分散度表征
G_N⁰ (Pa)10⁴~10⁶G''最小值处的G'缠结平台模量→Me
λz (s)0.01~1000η₀·Je⁰;低频G'/ωG''特征松弛时间,加工时间尺度
Ea (kJ/mol)20~150Arrhenius拟合aT(T)流动活化能,温度敏感性
tanδ峰温 (°C)≈Tg温度扫描1 Hz玻璃化转变温度,阻尼特性
γ_LVE (%)0.001~100应变扫描G'偏离5%点线性区边界,测试工作点基准