🙋
「Evans图」是什么？图中有两条线，为什么交点那么重要？
🎓
简单说，Evans图是在同一张图上看金属溶解（阳极·氧化）反应和还原反应（阴极·还原）的图。横轴是电位，纵轴是电流密度的对数（Tafel直线），两个反应平衡的地方就是「腐蚀电位Ecorr」，那里的电流密度就是「腐蚀电流密度icorr」。这个值越大，腐蚀越快。试试在模拟器里把阴极E₀滑块向右拖。你会看到交点向右上方移动，腐蚀速度上升。
🙋
真的！阴极电位升高腐蚀速度就上升。现实中什么情况下会这样？
🎓
比如铁管外面接触溶解氧多的水（开放水路），氧还原反应的平衡电位会升高（阴极直线向上移）。所以开放水路的腐蚀往往比密闭的快。相反，如果用脱氧水（比如锅炉），腐蚀就能减慢——这就是这个原理。你在「灵敏度分析」标签页试试把E₀c降低，能看到腐蚀速度下降。
🙋
「Tafel常数βc」增大时腐蚀速度怎么样？我在「Tafel灵敏度」标签页试了一下……腐蚀速度竟然下降了！
🎓
眼光敏锐！βc大的意思是阴极反应对电位变化「不敏感」。直线斜率变缓，交点向左下移。现实里βc大（≈120 mV/dec）是产氢反应的典型值，比氧还原（≈60 mV/dec）更难促进腐蚀。这个参数在腐蚀抑制剂设计里很重要——抑制剂会附着在表面改变反应机制，故意把β变大，从而抑制腐蚀。
🙋
腐蚀抑制剂能改变β，有意思！把金属换成「锌」时，腐蚀电位变得特别低。
🎓
这正是锌皮屋顶和镀锌的原理。锌的腐蚀电位（标准电位）比铁低，接触铁时锌会优先溶解、牺牲自己来保护铁——这叫「牺牲阳极」。船底和海底管道现在还装着锌块来防护。在模拟器里切换铁和锌，对比腐蚀电位的差（原电池电势差），直观感受就强多了。
电位与电流密度的关系用Tafel式表示。阳极（氧化）和阴极（还原）的过电压η分别为：
$$\eta_a = \beta_a \log\!\left(\frac{i}{i_{0,a}}\right), \quad \eta_c = -\beta_c \log\!\left(\frac{i}{i_{0,c}}\right)$$
$\beta_a, \beta_c$: Tafel常数 [V/decade]（通常 0.04〜0.12 V/decade）、$i_{0,a}, i_{0,c}$: 交换电流密度 [A/cm²]。腐蚀电位$E_{corr}$是阳极和阴极电流密度相等的点。
从腐蚀电流密度$i_{corr}$用Faraday定律计算腐蚀速度：
$$CR\, [\text{mm/year}] = \frac{i_{corr} \cdot M}{n \cdot F \cdot \rho} \times 3.156 \times 10^{10}$$
$M$: 金属原子量 [g/mol]、$n$: 溶解价数（例：Fe→Fe²⁺时n=2）、$F$: Faraday常数（96485 C/mol）、$\rho$: 密度 [g/cm³]。
电化学腐蚀工程的实际应用
管道·结构物使用寿命预测：海水中的钢铁结构和埋地管道的腐蚀速度用这个模型推算，用于耐用年数设计和检查周期决定。与实测icorr对比，制作腐蚀分布图。
防护技术设计：外部电流（阴极防护）或牺牲阳极（锌·铝）来降低Ecorr、把腐蚀速度降近零，这是设计的基础理论。在模拟器里把E₀c降低使Ecorr低于阳极平衡电位，能看到腐蚀停止的原理。
腐蚀抑制剂评估：定量化抑制剂对β值和交换电流密度i₀的影响，决定最有效的添加量。与电化学阻抗谱（EIS）等实验数据结合使用。
常见问题
在电化学池中连续改变试样金属电位、测量电流的「极化曲线测定」求得。在Ecorr±50〜300 mV范围内log(i) vs E作图基本成直线（Tafel区域）时，从斜率读取βa、βc。用自动极化装置（电势计）通常测定时间15〜30分钟。
一般腐蚀速度评价标准为：＜0.1 mm/year: 优秀、0.1〜0.5 mm/year: 良好、0.5〜1.0 mm/year: 可接受（薄壁部件要注意）、＞1.0 mm/year: 危险。不过判断还要看结构物厚度和要求使用年限。如海洋结构物板厚20mm，0.1 mm/year经过20年也会薄2mm，设计时要考虑。
涂层防护是物理隔绝金属和腐蚀环境，在Evans图上让两个电流密度（i₀）实质变近零。阴极防护是外部电流或牺牲阳极把Ecorr降至阳极平衡电位以下，热力学上禁止阳极反应。实务中两者常结合使用（如埋地管道涂层+外部电源阴极防护），涂层破损时阴极防护能补偿。
本工具是均匀腐蚀（全面腐蚀）的一维模型。实际腐蚀还有：①孔蚀（点蚀）：氯化物导致局部腐蚀，Evans图看不到、②应力腐蚀开裂（SCC）：应力和腐蚀相乘、③原电池腐蚀：异种金属接触、④缝隙腐蚀：氧浓差电池、⑤微生物诱发腐蚀（MIC）等。Evans图用于均匀腐蚀趋势把握和防护设计基础，局部腐蚀风险需要另外评估。
在Ecorr附近施加很小的ΔE（±10〜20 mV左右）测电流的非破坏性方法。用线性极化电阻Rp和腐蚀电流密度的关系式（Stern-Geary式）$i_{corr} = B/R_p$（B = βaβc/(2.303(βa+βc))）能迅速推算icorr而基本不伤害试样。本模拟器的Tafel常数βa、βc也能用来计算B。
不锈钢和铝在升高电位时会先增电流（活性溶解区），超过某电位后突然电流大幅下降——这是「钝化」。Evans图上阳极直线超过活性区后进入「钝化区」电流下降的形态。本模拟器假设简单Tafel直线不能表现钝化，但不锈钢腐蚀评价中比较钝化维持电位和腐蚀电位很重要。
电化学腐蚀模拟器简介
本模拟器的物理模型中，阳极反应和阴极反应的极化行为用Tafel式描述。阳极电流密度$i_a$和阴极电流密度$i_c$分别用平衡电位$E_{eq,a}$、$E_{eq,c}$和Tafel斜率$\beta_a$、$\beta_c$表示为：
$$
i_a = i_{0,a} \exp\left( \frac{E - E_{eq,a}}{\beta_a} \right), \quad i_c = i_{0,c} \exp\left( -\frac{E - E_{eq,c}}{\beta_c} \right)
$$
其中$i_{0,a}$、$i_{0,c}$是交换电流密度。腐蚀电位$E_{corr}$是满足$i_a = i_c$的电位，那时的电流密度称为腐蚀电流密度$i_{corr}$。腐蚀速度$v$（mm/year）用Faraday定律，基于金属原子量$M$、价数$n$、密度$\rho$、Faraday常数$F$算出：
$$
v = \frac{M}{\rho n F} \cdot i_{corr} \cdot 3.15 \times 10^7
$$
本模型让用户输入Tafel常数后求Evans图交点，实时评估腐蚀行为。
现实应用
产业实际案例：汽车行业中，车身镀锌钢板和铝合金部件的耐腐蚀性评价用本模拟器。比如汽车生产厂选排气系（消音器）材料时，用本模拟器提前预测异种金属接触腐蚀（原电池腐蚀）风险，反映到防护设计（涂层厚度、牺牲阳极配置）。海洋结构（海上风电基础部件）和化工厂管道材料选型时，也考虑实际环境盐分浓度和pH，用本工具推算腐蚀速度。
研究·教育应用：大学材料工学系或腐蚀防护工程实验课，学生拿实验得到的极化曲线数据输入本工具，对比理论值同时学习腐蚀机制。研究者开发新的耐腐蚀合金或腐蚀抑制剂时，把Tafel常数作参数，瞬间确认腐蚀电流密度变化，提高实验计划效率。
CAE分析的结合和实务定位：本模拟器不只单独做腐蚀诊断，还与有限元（FEM）结构CAE联动。比如汽车车身骨架模型中反映腐蚀造成的板厚减少率，预测长期使用后的强度下降——这是多物理场分析的一部分。实务上在设计阶段材料选择和防护规格的妥当性确认，以及既有设备剩余寿命评估（风险管理保养）的基础数据提供中发挥作用。
常见误解与注意事项
容易误认为「Tafel常数输入正确就能精确计算腐蚀速度」，但实际上Evans图是假设阳极·阴极反应为简单活化支配的理想系统，没考虑扩散支配、钝化皮膜影响、溶液电阻的IRドロップ等现实的复杂因素。所以计算值只是「推定值」，与实环境腐蚀速度可能有偏差。
「腐蚀电位贵（高）的话腐蚀速度小」容易误解，但实际腐蚀电位和腐蚀电流密度不是独立参数。比如阴极反应Tafel斜率陡峭的话，同一腐蚀电位腐蚀电流密度会大幅变化，单靠电位判断耐腐蚀性危险。
「mm/year数值小所以安全」的误解，局部腐蚀（孔蚀或晶界腐蚀）会在均匀腐蚀速度小的情况下突然穿透。本模拟器假设全面均匀腐蚀，局部腐蚀风险需另外关注。
使用指南
- 输入阳极反应的标准电极电位E0a（例：钢的氧化反应-0.61 V vs SHE）和交换电流密度i0a（例：1×10⁻⁶ A/cm²）
- 输入阴极反应的标准电极电位E0c（例：氧还原+0.40 V vs SHE）和交换电流密度i0c（例：1×10⁻⁸ A/cm²）
- 模拟器自动描画Evans图，计算显示腐蚀电位Ecorr·腐蚀电流密度icorr·腐蚀速度（mm/year）
具体计算例
软钢在海水中腐蚀时：阳极E0a=-0.61 V、i0a=2×10⁻⁶ A/cm²、Tafel常数βa=50 mV；阴极E0c=+0.10 V、i0c=1×10⁻⁸ A/cm²、βc=100 mV设定，得到腐蚀电位Ecorr≈-0.52 V、腐蚀电流密度icorr≈5×10⁻⁵ A/cm²，腐蚀速度约0.23 mm/year（使用原子量55.85 g/mol、密度7.87 g/cm³）。
实务注意点
- 电极电位因参照电极（SHE·CSE·Ag/AgCl）不同变化±0.2～0.4 V，输入前统一基准
- 交换电流密度i0对温度·pH·化学成分强烈依存（例：pH降低i0增加10倍），优先使用运行环境数据
- Tafel常数βa·βc文献值（βa通常40～60 mV、βc60～120 mV）变动±20%，要做灵敏度分析
- 离子强度高的环境（海水vs淡水）里i0急变，icorr可能数倍变化