我已经安装了 Optim1.6.1 包。我想在某些节点处使用约束为 u=1 来求解非线性有限元分析中的方程组。u=0 在某些节点。通常我发现 kmat 和 res 如下
[kmat,res]=systemmatrices(u,load,boundarycondition,material,constraints)
octave中的优化语法如下
[x, obj, info, iter, nf, lambda]=sqp(x0, @g, @h, upperbound, lowerbound, iter, tol)
- 在这里，“g”似乎是作为变量和方程的函数给出的目标函数。同样，“h”被认为是符号形式的方程组，例如 x1、x2 等。我在这个模板中看到了几个例子。我不清楚如何输入刚度矩阵、残差矩阵，因为这些矩阵是在单独的脚本中作为多个变量的函数计算的。
- 一般来说，机械载荷非线性 FEA 问题的目标函数是什么？它的能量残留？我能找到一个简单网格的简单示例吗？
我觉得很长一段时间都被这个问题困住了。如果这些太基础，请给我推荐一些阅读材料。