24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2667  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本科新能源科学与工程,一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 一志愿厦门大学化学学硕307求调剂 +10 y7czhao 2026-03-26 10/500 2026-03-28 14:23 by 唐沐儿
[考研] 求调剂 +6 争取九点睡 2026-03-28 6/300 2026-03-28 14:14 by Ainin_
[考研] 286求调剂 +4 丢掉懒惰 2026-03-27 7/350 2026-03-28 08:07 by baoball
[考研] 材料求调剂 一志愿哈工大总分298分,前三科223分 +5 dongfang59 2026-03-27 5/250 2026-03-28 04:53 by wxiongid
[考研] 340求调剂 +5 jhx777 2026-03-27 5/250 2026-03-28 04:18 by fmesaito
[考研] 求调剂推荐 材料 304 +15 荷包蛋hyj 2026-03-26 15/750 2026-03-28 04:13 by fmesaito
[考研] 一志愿南师大0703化学 275求调剂 +4 Ripcord上岸 2026-03-27 4/200 2026-03-27 17:00 by zhyzzh
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +9 vv迷 2026-03-22 9/450 2026-03-27 15:59 by 不吃魚的貓
[考研] 307求调剂 +8 超级伊昂大王 2026-03-24 9/450 2026-03-27 15:34 by 超级伊昂大王
[考研] 化学308分求调剂 +8 你好明天你好 2026-03-23 9/450 2026-03-27 14:01 by 杨光于青云
[考研] 一志愿北化求调剂 +3 Jsman 2026-03-22 3/150 2026-03-26 21:06 by ajpv风雷
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +5 星沉uu 2026-03-26 6/300 2026-03-26 19:02 by macy2011
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 302求调剂 +4 锦衣卫藤椒 2026-03-25 4/200 2026-03-25 16:29 by 功夫疯狂
[考研] 求调剂一志愿武汉理工大学材料工程(085601) +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 277分求调剂,跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 求调剂一志愿海大,0703化学学硕304分,有大创项目,四级已过 +6 幸运哩哩 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
信息提示
请填处理意见