 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 623 | 回復(fù): 1 | |||
LG_Maxwell金蟲 (著名寫手)
|
[交流]
關(guān)于matlab計(jì)算精度的問(wèn)題
|
|
我想由已知的點(diǎn)應(yīng)力矩陣，計(jì)算出三個(gè)主應(yīng)力和主方向，理論很成熟就是代公式即可，但是結(jié)果上出現(xiàn)了問(wèn)題，懷疑為計(jì)算精度問(wèn)題，具體如下： 計(jì)算主應(yīng)力： 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 計(jì)算主方向：（以主應(yīng)力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; 這套方法計(jì)算文獻(xiàn)算例也好，自己隨便輸入的點(diǎn)應(yīng)力矩陣（類似[1，2，3；2，4，5；3，5，6]這種）也好，都沒(méi)問(wèn)題；但是計(jì)算我的實(shí)驗(yàn)數(shù)據(jù)時(shí)得到的主應(yīng)力還算合理，得到的主方向問(wèn)題很大，即三個(gè)方向矢量點(diǎn)乘的結(jié)果完全不能認(rèn)為是0，這不符合物理規(guī)律。 懷疑是程序計(jì)算精度不足導(dǎo)致的，因?yàn)閷?shí)驗(yàn)測(cè)得的點(diǎn)應(yīng)力，六個(gè)分量之間數(shù)量級(jí)從e-4到e-12不等，差距很大，公式又有乘除運(yùn)算，可能導(dǎo)致浮點(diǎn)計(jì)算誤差很大。 所以想請(qǐng)教各位大佬，有什么辦法可以提高運(yùn)算精度，或者修改算式回避這個(gè)問(wèn)題，再或者有沒(méi)有擅長(zhǎng)應(yīng)力理論的大佬，發(fā)現(xiàn)是其他問(wèn)題導(dǎo)致的結(jié)果錯(cuò)誤。 謝謝謝謝！  |
回復(fù)此樓» 猜你喜歡
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生
已經(jīng)有4人回復(fù)
-
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有246人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有1人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
LG_Maxwell
金蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1230.6
- 散金: 704
- 沙發(fā): 1
- 帖子: 1014
- 在線: 37小時(shí)
- 蟲號(hào): 22140562
- 注冊(cè): 2020-05-15
- 性別: GG
- 專業(yè): 光學(xué)
2樓2020-06-08 11:30:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 LG_Maxwell 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 哈工大材料324求調(diào)劑 +3 | 閆旭東 2026-03-10 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +4 | 小木蟲085600 2026-03-09 | 7/350 |
|
|
|
[考研] 復(fù)試調(diào)劑 +6 | 呼呼？~+123456 2026-03-08 | 7/350 |
|
|
|
[考研] 考研材料與化工，求調(diào)劑 +3 | 戲精丹丹丹 2026-03-09 | 3/150 |
|
|
|
[考研] 0703求調(diào)劑 +10 | c297914 2026-03-06 | 11/550 |
|
|
|
[考研] 求調(diào)劑 +3 | 鶴遨予卿 2026-03-09 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +4 | xxxcm 2026-03-08 | 7/350 |
|
|
|
[考研] 0856材料與化工290求調(diào)劑 +7 | Nebulala 2026-03-08 | 8/400 |
|
|
|
[考研] 286求調(diào)劑 +12 | Faune 2026-03-06 | 14/700 |
|
|
|
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過(guò)六級(jí)有競(jìng)賽求調(diào)劑 +5 | galaxary 2026-03-07 | 7/350 |
|
|
|
[考研] 一志愿211 化學(xué)305分求調(diào)劑 +3 | 0703楊悅305分 2026-03-05 | 3/150 |
|
|
|
[考研] 一志愿211 085600 280數(shù)二英二求調(diào)劑 +3 | 月山斜 2026-03-06 | 3/150 |
|
|
|
[考研] 267化工調(diào)劑求助 +7 | 聰少OZ 2026-03-04 | 7/350 |
|
|
|
[考研] 復(fù)試調(diào)劑 +7 | 呼呼？~+123456 2026-03-05 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +5 | danyyyy 2026-03-04 | 5/250 |
|
|
|
[考研] 材料調(diào)劑 +4 | L9370 2026-03-05 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +9 | 土土小蟲 2026-03-03 | 11/550 |
|
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|