| 查看: 2652 | 回复: 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 ] |
回复此楼» 猜你喜欢
一志愿北京化工大学 070300 学硕 336分 求调剂
已经有5人回复
-
323求调剂
已经有6人回复
-
352求调剂
已经有3人回复
-
一志愿东华大学化学070300，求调剂
已经有8人回复
-
277材料科学与工程080500求调剂
已经有7人回复
-
317求调剂
已经有18人回复
-
293求调剂
已经有5人回复
-
280分求调剂 一志愿085802
已经有7人回复
-
0854电子信息求调剂
已经有3人回复
-
263求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在Ansys时间后处理的载荷位移曲线中怎样提取最大载荷对应的横坐标值
已经有6人回复
- FreeFem++ 中知道网格信息，如何得到每个顶点相邻的顶点编号
已经有8人回复
young_sir
银虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 3833.2
- 红花: 2
- 帖子: 525
- 在线: 297.2小时
- 虫号: 3579112
- 注册: 2014-12-05
- 性别: GG
- 专业: 声学
4楼2015-12-22 15:41:59
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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
5楼2016-03-27 09:43:45
