额，大致想了一下，应该不算太难算。但是要再添加一些假定：冰块表面的消减是均一的。
否则接触空气部分和接触水部分的换热量不同，会导致上下消解的速度不同，那冰块就变形了，没电脑的情况下有点难算。将所有换热带来的质量损耗均匀的分摊在表面每个部分，这样对于分开冰块的水上水下两部分比较有利（你就当冰块一直在水面转动吧~）
接下来就简单些了：
dM/dt=A1*h1(t1-t2)+A2*h2*(t1-t3)
A1+A2=(M/p)^(1/3)
M为冰块随时间变化的函数M(t)，A1为冰块接触空气部分的面积，是时间的函数A1（t），A2为接触水部分的面积，也是时间的函数A2（t）；t1是冰的温度 t2是空气温度 t3是水的温度，按设定都是常数 p为密度 。
冰从冰柜拿出来放到室温开始算，一直到全部化成水的时间
