| 查看: 2635 | 回复: 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 ] |
回复此楼» 猜你喜欢
材料调剂
已经有4人回复
-
面上模板改不了页边距吧？
已经有6人回复
-
307求调剂
已经有6人回复
-
304求调剂
已经有5人回复
-
317一志愿华南理工电气工程求调剂
已经有8人回复
-
272求调剂
已经有3人回复
-
化工专硕348，一志愿985求调剂
已经有6人回复
-
292求调剂
已经有3人回复
-
290求调剂
已经有6人回复
-
295求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在Ansys时间后处理的载荷位移曲线中怎样提取最大载荷对应的横坐标值
已经有6人回复
- FreeFem++ 中知道网格信息，如何得到每个顶点相邻的顶点编号
已经有8人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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 (幼儿园)
- 金币: 3828.2
- 红花: 2
- 帖子: 524
- 在线: 296.4小时
- 虫号: 3579112
- 注册: 2014-12-05
- 性别: GG
- 专业: 声学
4楼2015-12-22 15:41:59
5楼2016-03-27 09:43:45
dhxyadeng
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 15.5
- 帖子: 6
- 在线: 1.3小时
- 虫号: 7444396
- 注册: 2017-11-08
- 性别: GG
- 专业: 力学中的基本问题和方法
6楼2020-07-13 15:25:43
