| 6 | 1/1 | 返回列表 |
| 查看: 2657 | 回復(fù): 5 | |||
pmsl666至尊木蟲 (著名寫手)
|
[求助]
FreeFem++ 怎么將所有單元的所有節(jié)點(diǎn)的坐標(biāo)及位移值輸出 已有1人參與
|
|
我使用FreeFem++ 求解一個偏微分方程,采用的是P2單元,P2單元按理說應(yīng)該是6節(jié)點(diǎn)三角形單元,我想按單元順序,將所有單元的所有節(jié)點(diǎn)的坐標(biāo)及位移值輸出, 我看了幫助文檔,貌似里面只有三角形3個頂點(diǎn)坐標(biāo)及位移輸出,怎么樣才能將三角形頂點(diǎn)以及中點(diǎn)坐標(biāo)位移都輸出到txt文件中。特此求助。 { ofstream ff("grapht.txt" ); for (int i=0;i<Th.nt;i++) { for (int j=0; j <3; j++) ff<<Th[j].x << " "<< Th[j].y<< " "<<u[][Vh(i,j)]<<endl; ff<<Th[0].x << " "<< Th[0].y<< " "<<u[][Vh(i,0)]<<endl <<endl<<endl; } } [ Last edited by feixiaolin on 2015-6-12 at 16:00 ] |

新蟲 (初入文壇)
{ ofstream ff("graph1.txt" ;for (int i=0;i<Th.nt;i++) {ff<<Th[0].x << " "<< Th[0].y<< " "<<u(Th[0].x,Th[0].y)<<endl; ff<<(Th[0].x+Th[1].x)/2 << " "<< (Th[0].y+Th[1].y)/2<< " "<<u((Th[0].x+Th[1].x)/2,(Th[0].y+Th[1].y)/2)<<endl; ff<<Th[1].x << " "<< Th[1].y<< " "<<u(Th[1].x,Th[1].y)<<endl; ff<<(Th[1].x+Th[2].x)/2 << " "<< (Th[1].y+Th[2].y)/2<< " "<<u((Th[1].x+Th[2].x)/2,(Th[1].y+Th[2].y)/2)<<endl; ff<<Th[2].x << " "<< Th[2].y<< " "<<u(Th[2].x,Th[2].y)<<endl; ff<<(Th[0].x+Th[2].x)/2 << " "<< (Th[0].y+Th[2].y)/2<< " "<<u((Th[0].x+Th[2].x)/2,(Th[0].y+Th[2].y)/2)<<endl; ff<<Th[0].x << " "<< Th[0].y<< " "<<u(Th[0].x,Th[0].y)<<endl <<endl<<endl; } } u(x,y)就能直接輸出[x,y]點(diǎn)的位移,邊界中點(diǎn)坐標(biāo)可以根據(jù)三角形頂點(diǎn)求得 |

新蟲 (初入文壇)
銀蟲 (正式寫手)
新蟲 (初入文壇)
鐵蟲 (初入文壇)

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +4 | jiajunX 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 工科0856求調(diào)劑 +4 | 沐析汀汀 2026-03-21 | 4/200 |
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +8 | 材料學(xué)257求調(diào)劑 2026-03-23 | 9/450 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +9 | niko- 2026-03-19 | 10/500 |
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|