1

我在 Abaqus 中构建模型时的一个常见做法是拟合材料属性。例如,我尝试了所有可能的材料属性,并查看模型给出的表面变形,然后找出最符合我们实验观察的那个。实际上,我将该值与模型输出和实验数据进行比较,得到一个 R 平方值并尝试最小化 -1 的值。* R2。

我一直在使用 scipy 优化工具箱在 Abaqus 中执行此操作。但是,有一个问题:在某些情况下,模型不会与优化器尝试的某些给定参数收敛。在这些情况下,我应该为 R2 设置什么值?我应该将其设置为-1。* numpy.inf 或 -1。* numpy.nan(假设将 numpy 导入为 np)?

此外,在某些情况下,我使用不支持一般约束的优化函数,例如模数_1 > 模数_2;如果它要求我提交模数_1 <=模数_2的工作,我可以只返回-1。* np.inf 或 -1。* np.nan 作为惩罚?

问题的发生是因为没有办法知道模型在参数空间中无法收敛的先验位置。任何帮助将不胜感激。太感谢了!

4

0 回答 0