24小时热门版块排行榜    

查看: 2654  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 北科281学硕材料求调剂 +8 tcxiaoxx 2026-03-20 8/400 2026-03-23 10:16 by 冠c哥
[考研] 284求调剂 +6 Zhao anqi 2026-03-22 6/300 2026-03-23 09:23 by king123!
[考研] 310求调剂 +4 baibai1314 2026-03-16 4/200 2026-03-22 20:19 by edmund7
[考研] 287求调剂 +8 晨昏线与星海 2026-03-19 9/450 2026-03-22 17:01 by i_cooler
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
信息提示
请填处理意见