我怎么用对分法求解主应力?
2视图(30天)
显示旧的评论
对于任何应力(或应变)张量为3 x3或任何给定的矩阵,利用特征方程,我制定:
I1σ=σ(1,1)+(2,2)+σ(3,3);
I2 =σ(1,1)*σ(2,2)+σ(1,1)*σ(3,3)+σ(2,2)*σ(3、3)σ(1、2)^ 2σ(1、3)^ 2σ(2、3)^ 2;
I3σ=σ(1,1)*(2,2)*σ(3,3)+ 2 *(σ(1、2)*σ(2、3)*σ(3,1))σ(1、2)^ 2 *σ(3、3)σ(2、3)σ^ 2 *(1,1)σ(3,1)σ^ 2 * (2,2);
x =根([1 i1 + I2 i3]);
我不确定如何使用二分法的时候这样做. .确定最大和最小主应力为任意应力张量
我之前有这种二分法的使用:
函数x = BisectionRoot_smith(有趣,a, b) %使用二分法查找函数的根有趣在a和b之间的间隔%。功能停止时公差小于% 10 ^ 6
tol_error = 1 e-6;如果((有趣的(a) *乐趣(b)) > 0) disp的错误:对分法需要间隔绑定一个根的其他聚合= 0;%初始化聚合变量
%如果区间端点属于可接受错误,返回
%合适作为根
如果< tol_error abs(有趣的(a))
x = a;
聚合= 1;
elseif< tol_error abs(有趣的(b))
x = b;
聚合= 1;
结束
%连续循环使用二分,直到一个可接受的公差
%获得
而~聚合
Xmid = 0.5 * (b + a);
如果abs(乐趣(Xmid)) < tol_error
%发现融合解决方案
x = Xmid;
聚合= 1;
其他的
%重复二等分的过程
如果乐趣(Xmid) *乐趣(a) > 0
一个= Xmid;
其他的
b = Xmid;
结束
结束
结束
结束
但不知道如何连接这两个
