 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 621 | 回復: 1 | |||
LG_Maxwell金蟲 (著名寫手)
|
[交流]
關于matlab計算精度的問題
|
|
我想由已知的點應力矩陣，計算出三個主應力和主方向，理論很成熟就是代公式即可，但是結果上出現(xiàn)了問題，懷疑為計算精度問題，具體如下： 計算主應力： 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]這種）也好，都沒問題；但是計算我的實驗數(shù)據(jù)時得到的主應力還算合理，得到的主方向問題很大，即三個方向矢量點乘的結果完全不能認為是0，這不符合物理規(guī)律。 懷疑是程序計算精度不足導致的，因為實驗測得的點應力，六個分量之間數(shù)量級從e-4到e-12不等，差距很大，公式又有乘除運算，可能導致浮點計算誤差很大。 所以想請教各位大佬，有什么辦法可以提高運算精度，或者修改算式回避這個問題，再或者有沒有擅長應力理論的大佬，發(fā)現(xiàn)是其他問題導致的結果錯誤。 謝謝謝謝！  |
回復此樓» 猜你喜歡
重慶交通大學光子學微結構與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復
-
一志愿鄭大材料學碩298分，求調(diào)劑
已經(jīng)有6人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有165人回復
-
尋合作：應力腐蝕多尺度模擬
已經(jīng)有5人回復
-
考研交流
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
幫我的英語口語老師找學生
已經(jīng)有2人回復
LG_Maxwell
金蟲 (著名寫手)
- 應助: 0 (幼兒園)
- 金幣: 1230.6
- 散金: 704
- 沙發(fā): 1
- 帖子: 1014
- 在線: 37小時
- 蟲號: 22140562
- 注冊: 2020-05-15
- 性別: GG
- 專業(yè): 光學
2樓2020-06-08 11:30:33
相關版塊跳轉
我要訂閱樓主 LG_Maxwell 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] A區(qū)一本交叉課題組，低分調(diào)劑，招收機械電子信息通信等交叉方向 +4 | lisimayy 2026-03-04 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 282求調(diào)劑 +3 | 夕～日 2026-03-05 | 3/150 |
|
|
|
[考研] 0703 東華大學 理學化學均過 a 線，總分 281 求調(diào)劑 +4 | 香香 Lu 2026-03-03 | 7/350 |
|
|
|
[考研] 0856材料求調(diào)劑 +17 | hyf hyf hyf 2026-02-28 | 18/900 |
|
|
[考研] 一志愿西交大材料學碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
|
[考研] 一志愿東北大學化學314分求調(diào)劑 +5 | lr1212.. 2026-03-02 | 5/250 |
|
|
|
[考研]
材料工程專碩283求調(diào)劑
 5+8
|
，�。� 2026-03-02 | 10/500 |
|
|
|
[考研] 11408，學碩276求調(diào)劑 +3 | 崔wj 2026-03-02 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | repeatt?t 2026-02-28 | 7/350 |
|
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
|
[考研] 一志愿華南理工大學材料與化工326分，求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
|
[考研] 281求調(diào)劑 +5 | 2026計算機_誠心 2026-03-01 | 8/400 |
|
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
|
[考研] 化工299分求調(diào)劑 一志愿985落榜 +5 | 嘻嘻(*^ω^*) 2026-03-01 | 5/250 |
|
|
|
[基金申請]
剛錄用，沒有期刊號，但是在線可看的論文可以放為代表作嗎
10+3
|
arang1 2026-03-01 | 3/150 |
|
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 3/150 |
|
|
|
[考研]
 尋找調(diào)劑
+4
|
LYidhsjabdj 2026-02-28 | 4/200 |
|