| 查看: 663 | 回复: 1 | |||
LG_Maxwell金虫 (著名写手)
|
[交流]
关于matlab计算精度的问题
|
|
我想由已知的点应力矩阵，计算出三个主应力和主方向，理论很成熟就是代公式即可，但是结果上出现了问题，怀疑为计算精度问题，具体如下： 计算主应力： for i=1:1:r T{i}=[xx(i),xy(i),zx(i);xy(i),yy(i),yz(i);zx(i),yz(i),zz(i)]; [~,principalstress{i}]=eig(T{i}); p1(i)=principalstress{i}(1,1); p2(i)=principalstress{i}(2,2); p3(i)=principalstress{i}(3,3); end 计算主方向：（以主应力1为例） l1=1; m1=l1*(zx*xy-yz*(xx-p1))/(xy*yz-zx*(yy-p1)); n1=l1*((xx-p1)*(yy-p1)-xy*xy)/(xy*yz-zx*(yy-p1)); abcd1=sqrt(l1^2+m1^2+n1^2);%归一化 l1=l1/abcd1; m1=m1/abcd1; n1=n1/abcd1; 这套方法计算文献算例也好，自己随便输入的点应力矩阵（类似[1，2，3；2，4，5；3，5，6]这种）也好，都没问题；但是计算我的实验数据时得到的主应力还算合理，得到的主方向问题很大，即三个方向矢量点乘的结果完全不能认为是0，这不符合物理规律。 怀疑是程序计算精度不足导致的，因为实验测得的点应力，六个分量之间数量级从e-4到e-12不等，差距很大，公式又有乘除运算，可能导致浮点计算误差很大。 所以想请教各位大佬，有什么办法可以提高运算精度，或者修改算式回避这个问题，再或者有没有擅长应力理论的大佬，发现是其他问题导致的结果错误。 谢谢谢谢！  |
回复此楼» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有119人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
LG_Maxwell
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1230.6
- 散金: 704
- 沙发: 1
- 帖子: 1014
- 在线: 37小时
- 虫号: 22140562
- 注册: 2020-05-15
- 性别: GG
- 专业: 光学
2楼2020-06-08 11:30:33