24小时热门版块排行榜    

查看: 2653  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 招08考数学 +4 laoshidan 2026-03-20 11/550 2026-03-23 08:47 by laoshidan
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
信息提示
请填处理意见