| 查看: 2842 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
pmsl666至尊木虫 (著名写手)
|
[求助]
FreeFem++ 怎么将所有单元的所有节点的坐标及位移值输出 已有1人参与
|
|
|
我使用FreeFem++ 求解一个偏微分方程，采用的是P2单元，P2单元按理说应该是6节点三角形单元，我想按单元顺序，将所有单元的所有节点的坐标及位移值输出， 我看了帮助文档，貌似里面只有三角形3个顶点坐标及位移输出，怎么样才能将三角形顶点以及中点坐标位移都输出到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 ] |
回复此楼» 猜你喜欢
某些机构，以效率低为荣，以效率低作为存在感
已经有9人回复
-
今天务委会开完了，明天出结果吗
已经有3人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有6人回复
-
放榜前的不淡定
已经有18人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
-
人气不行了
已经有11人回复
-
2026国自然函评费到账
已经有22人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在Ansys时间后处理的载荷位移曲线中怎样提取最大载荷对应的横坐标值
已经有6人回复
- FreeFem++ 中知道网格信息，如何得到每个顶点相邻的顶点编号
已经有8人回复
5楼2016-03-27 09:43:45
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
pmsl666: 金币+20, ★★★★★最佳答案 2015-06-15 17:57:15
感谢参与，应助指数 +1
pmsl666: 金币+20, ★★★★★最佳答案 2015-06-15 17:57:15
{ 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]点的位移，边界中点坐标可以根据三角形顶点求得 |
2楼2015-06-15 11:31:32
3楼2015-11-05 12:27:42
young_sir
银虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 3838.2
- 红花: 2
- 帖子: 526
- 在线: 297.2小时
- 虫号: 3579112
- 注册: 2014-12-05
- 性别: GG
- 专业: 声学
4楼2015-12-22 15:41:59
