 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 620 | 回復(fù): 1 | |||
LG_Maxwell金蟲(chóng) (著名寫(xiě)手)
|
[交流]
關(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)力矩陣（類(lèi)似[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ù)此樓» 猜你喜歡
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲(chǔ)能器件方向博士
已經(jīng)有13人回復(fù)
-
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有213人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有3人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
LG_Maxwell
金蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1230.6
- 散金: 704
- 沙發(fā): 1
- 帖子: 1014
- 在線: 37小時(shí)
- 蟲(chóng)號(hào): 22140562
- 注冊(cè): 2020-05-15
- 性別: GG
- 專(zhuān)業(yè): 光學(xué)
2樓2020-06-08 11:30:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 LG_Maxwell 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 320材料一志愿西工大專(zhuān)碩英二數(shù)二 有兩年光伏工作經(jīng)驗(yàn) +3 | 錘某人 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 材料與化工，291分，求調(diào)劑 +6 | 咕嚕咕嚕123123 2026-03-04 | 6/300 |
|
|
|
[考研] 材料復(fù)試調(diào)劑 +6 | 學(xué)材料的點(diǎn) 2026-03-01 | 7/350 |
|
|
|
[考研]
 一志愿中南大學(xué)理學(xué)化學(xué)
+6
|
15779376950 2026-03-01 | 8/400 |
|
|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過(guò) a 線，總分 281 求調(diào)劑 +3 | 香香 Lu 2026-03-03 | 5/250 |
|
|
[考研] 0857求調(diào)劑 +4 | CDDD1013 2026-03-04 | 4/200 |
|
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
 5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
|
[考研] 一志愿西工大材料學(xué)，英一數(shù)二總分321分，求調(diào)劑。 +3 | zz05zz 2026-03-02 | 5/250 |
|
|
|
[考研] 085600求調(diào)劑 +4 | LRZZZZZZ 2026-03-02 | 6/300 |
|
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
|
[考研] 一志愿東北大學(xué)化學(xué)314分求調(diào)劑 +5 | lr1212.. 2026-03-02 | 5/250 |
|
|
|
[考研] 一志愿中科大能動(dòng)297求調(diào)劑，本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
|
[考研]
計(jì)算機(jī)學(xué)碩 分?jǐn)?shù)285求調(diào)劑
+4
|
glwshine 2026-03-02 | 5/250 |
|
|
|
[考研] 267求調(diào)劑 +6 | 釣魚(yú)佬as 2026-03-02 | 6/300 |
|
|
|
[考研] 281電子信息求調(diào)劑 +3 | jhtfeybgj 2026-03-02 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
|
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5 | 半截的詩(shī)0927 2026-03-02 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +3 | cgyzqwn 2026-03-01 | 7/350 |
|
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|