Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
|
![]() |
#1 |
高级会员
注册日期: 2008-01-02
帖子: 247
声望力: 22 ![]() |
![]()
表达式
D=solve('S^3-I1*S^2+I2*S-I3') 中,原I1,I2,I3的值并没有传入方程 S^3-I1*S^2+I2*S-I3=0 因为solve是符号求解器。 破解之法是求数值解: clear all clc s=eye(3); %计算应力不变量 I1=trace(s);%计算第一应力不变量I1 I2=det(s(2:3,2:3))+det(s(1:2:3,1:2:3))+det(s(1:2,1:2));%计算I2 I3=det(s);%计算I3 %解三次方程,求主应力。 c = [1 -I1 I2 -I3];%这就是三次多项式 S^3-I1*S^2+I2*S-I3 D=roots(c); |
![]() |
![]() |
![]() |
#2 |
初级会员
注册日期: 2009-02-28
年龄: 37
帖子: 4
声望力: 0 ![]() |
![]()
非常感谢mathjiang的解答,醍醐灌顶!对有些东西还是认识太浅,还需加强学习。
|
![]() |
![]() |