去年帮东莞的一个加工厂调铝合金微弧氧化的膜层附着力，一开始死磕单因素试验，试了三周，电解液温度从20℃调到40℃，电流密度从5A/dm²爬到15A/dm²，每次都觉得摸到临界点了，一到产线放大就掉膜。
真的烦。
之前总觉得单因素把每个变量的影响摸透就行，直到产线干了二十年的老师傅提醒，温度高了之后电解液里的硅酸钠会析出来，电流的实际分布根本不是你实验室小槽子测的那个数，两个变量根本不是独立的。
别把它当成”高级版正交试验”
很多人刚接触响应曲面法的时候，觉得这玩意儿不就是多做几轮试验，把回归方程算出来，找个最优解？本质上和正交试验没差。
错。
当时一开始设计试验的时候，选了三个变量：电解液温度、电流密度、处理时间，响应值是膜层的划痕临界载荷。按教材里的Box-Behnken设计选了17个试验点，实验室小槽做出来的拟合R²有0.92，算出来的最优参数是32℃，10A/dm²，12min。
结果上产线就崩了。
小槽的电解液循环靠磁力搅拌，产线的是大流量泵，槽边和槽中心的温度差能到6℃，你实验室测的是均匀体系的响应，到了产线根本没有那个均匀度。说起来都好笑，当时拿着拟合的曲面图给产线主管看，人家扫了一眼就问，你这图里的温度，是槽子里哪个点的温度？我当场卡壳。
拟合出来的最优解，要先卡边界
很多人拿到软件算出来的最优值，第一反应是照着试，从来不想这个解落在什么限制条件里。当时那个拟合方程算出来的最优划痕载荷是32N，但是你算一下对应的电解液寿命：32℃的时候，硅酸钠的分解速率比25℃高40%（估算值，没做长期跟踪），每两周就要补一次主盐，产线的运维成本直接翻三倍。
响应曲面的边界外推是绝对不能碰的。比如当时为了找更高的附着力，有人说把电流密度再往上调2A，看看能不能到40N的载荷——用现有方程外推的结果看起来很漂亮，实际测出来膜层直接被电流击穿，表面满是灼蚀坑，根本没法用。
真正能用的最优解，从来不是曲面上那个数学最大值，是把所有工艺限制条件叠上去之后的那个可行域交点。我们最后选的参数是28℃，9A/dm²，10min，算出来的划痕载荷是27N，比理论最优低一点，但是电解液的温度差在产线的允许波动范围内，主盐的补加周期能到一个月，工人操作的时候不用频繁调整冷却系统，良率能稳在93%以上。
别迷信软件跑出来的结果
软件不知道你实验室的小槽子和产线的槽子差多少，不知道你的搅拌系统能不能做到±1℃的温度均匀性，不知道工人换挂具的时候会不会把工件贴到槽壁上。
当时做试验的时候，有一个中心点的平行样，测出来的载荷比其他平行样低了8N，一开始软件自动把这个点归为异常值删掉了，我当时觉得不对，翻了试验记录，那天刚好实验室的冷却水坏了，槽液温度比设定值高了3℃——这个点根本不是异常，是真实的温度波动带来的响应变化，要是直接删掉，最后拟合出来的曲面在温度轴上的坡度就会偏缓，到了产线温度飘一点，结果就差很多。
这个坑我踩过。
现在很多人做响应曲面，恨不得把所有试验点都按软件提示的顺序排，连试验顺序的随机化都懒得做，前面的槽液还没稳定就开始测，后面的结果飘得没边，最后还怪响应曲面法不准。
说实话，响应曲面法本质上就是帮你把几个变量之间的交叉影响理清楚，不是什么能直接给出黄金参数的魔法。它比单因素、正交试验多走的那一步，就是告诉你变量之间不是各管各的——温度变了，电流密度的最优值也会跟着变，不是你固定了一个参数就能一直用。
到现在我们也没完全摸透那个大槽的液流分布规律，只是每次做试验之前，先在槽里五个固定点测一遍温度，把偏差超过2℃的点记下来，算进响应的修正项里。
上次去那个加工厂，他们现在调工艺的时候，还是会拿响应曲面先做实验室的小试，但是每拿到一个拟合结果，先对着产线的设备参数列个三页的限制条件，再也不会拿着实验室的最优值直接往产线上搬了。
