24小时热门版块排行榜    

查看: 2655  |  回复: 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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

独学则无友,孤陋必寡闻。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangwww123

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hli47

新虫 (初入文坛)

请问啊,我安装了FreeFem++,不知道怎么使用它来编辑我的数据?一点击快捷图标,先弹出一个黑框框,然后就又出来一个选择算例的窗口,这是怎么回事呢?我自己要是编辑一个算例,在哪里编辑数据呢
3楼2015-11-05 12:27:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

young_sir

银虫 (正式写手)

引用回帖:
3楼: Originally posted by hli47 at 2015-11-05 12:27:42
请问啊,我安装了FreeFem++,不知道怎么使用它来编辑我的数据?一点击快捷图标,先弹出一个黑框框,然后就又出来一个选择算例的窗口,这是怎么回事呢?我自己要是编辑一个算例,在哪里编辑数据呢

请问您解决这个问题了吗?我也刚接触,完全不知道怎么入手…
4楼2015-12-22 15:41:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绿色de清晨

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by young_sir at 2015-12-22 15:41:59
请问您解决这个问题了吗?我也刚接触,完全不知道怎么入手…...

点击这个图标,再编码
FreeFem++ 怎么将所有单元的所有节点的坐标及位移值输出
4.png

5楼2016-03-27 09:43:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dhxyadeng

铁虫 (初入文坛)

试试这几个函数
// savemesh(Th,"nameMesh.mesh";
// ffSaveVh(Th,Vh,"nameVh.txt";
// ffSaveData(u,"nameData.txt";
6楼2020-07-13 15:25:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pmsl666 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 08工科 320总分 求调剂 +11 梨花珞晚风 2026-03-17 11/550 2026-03-22 17:42 by luoyongfeng
[考研] 寻找调剂 +4 倔强芒? 2026-03-21 4/200 2026-03-22 16:14 by 木托莫露露
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
信息提示
请填处理意见