用python自动计算拟合材料曲线

        之前我曾经发表过如何用JohnsonCook模型延申应力应变曲线到1以及有效应力应变曲线是怎么来的这两篇文章，后来处理的材料曲线实在是太多了。索性就写了个自动处理材料曲线的软件。

        写着写着把，突然想把材料拟合这部分功能也加进去，然后就发现了python拟合曲线竟然如此简单高效，比我之前在Excel用平方差之和趋近于零，用规划求解器拟合快多了。今天就介绍这个吧。

这是一条材料曲线，黄色部分的b,n值就是我们要拟合计算得到的johnsoncook公式的关键的两个值。这里我已经在Excel里用规划求解器得出了

接下来我们用Python来进行拟合

第一步导入我们要用到的包

打开处理材料曲线的Excel表格，导入要拟合的有效应力应变曲线。

定义johnsoncook公式的函数

用popt,pcov=curve_fit()函数进行拟合

这里用matplotlib把原始的有效应力应变曲线和拟合后的曲线以散点图的方式表现出来，用作对比。

拟合的还可以，就是最后那段有点问题

这里是拟合后得到输出的b,n值。

对比一下，没什么差别。