24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2524  |  回复: 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的回帖

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的回帖
查看全部 6 个回答

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的回帖
信息提示
请填处理意见