需将有限元模型写为 Python文件,并基于 Python 的二次开发提取模型计算的最大主应力分布和接触应力写入文本文档,后将整个有限元计算流程写入批处理(bat)文件中,再结合优化设计软件 ISIGHT的 SIMCODE 模块,进行自动化提交运算并进行优化设计。其次,利用 MIGA 算法进行全局优化,以最小化 SM 和 SC 为优化目标进行多目标优化设计,寻找到最优的应力过渡层的厚度、层数和幂函数的分布指数,软件流程图如图 5-3(b)所示,值得注意的是,应力过渡层的厚度在取值在 0-4 mm,层数在 10 层以下,幂函数的分布指数取值为0-10,综上,上述优化设计流程的如式 5-3 所示。
现在已经有初始遗传算法可以使用,只需要更改球的尺寸,模型的整体尺寸,球应该不是刚体了,应该是牙釉质的参数,输出球和梯度材料两个部件上的接触应力。



