实际生产中存在不少卧式圆筒形的容器,液位计算也不是个简单的线性关系,实际上还相当复杂。
已知下图高度x求取面积要相对简单些,而反过来根据面积求x却异常麻烦。
Maple解:
当x=0时,解值为-Pi*r^2/2 ,定积分则为 A+ Pi*r^2/2 。
可以使用逼近法来根据已知面积来解x 。
At := Pi*(r**2) ;
x := r ;
loop
As := ( (r**2)*arctan((-r+x)/((2.0*r*x-x**2)**0.5))-(r-x)*((2.0*r*x-x**2)**0.5) ) + 0.5*Pi*(r**2) ;
若As < At ,
K := (At-As)/At ;
x := x + x * K ;
若As>=At ,计算结束。