版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

pmsl666

至尊木虫 (著名写手)

应助: 2 (幼儿园)
金币: 24199.5
散金: 199
红花: 56
帖子: 1179
在线: 684.5小时
虫号: 945274
注册: 2010-01-20
性别: GG
专业: 金属材料的合金相、相变及

[求助] 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 ]

回复此楼

» 猜你喜欢

请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

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

在Ansys时间后处理的载荷位移曲线中怎样提取最大载荷对应的横坐标值已经有6人回复
FreeFem++ 中知道网格信息，如何得到每个顶点相邻的顶点编号已经有8人回复

独学则无友，孤陋必寡闻。

1楼2015-06-12 15:23:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yangwww123

新虫 (初入文坛)

应助: 3 (幼儿园)
金币: 99
帖子: 3
在线: 10.7小时
虫号: 1683283
注册: 2012-03-11
性别: GG
专业: 偏微分方程

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 41.8
帖子: 13
在线: 48小时
虫号: 1837103
注册: 2012-05-28
专业: 应用数学方法

请问啊，我安装了FreeFem++,不知道怎么使用它来编辑我的数据？一点击快捷图标，先弹出一个黑框框，然后就又出来一个选择算例的窗口，这是怎么回事呢？我自己要是编辑一个算例，在哪里编辑数据呢

赞一下

回复此楼

3楼2015-11-05 12:27:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

young_sir

银虫 (正式写手)

应助: 3 (幼儿园)
金币: 3823.2
红花: 2
帖子: 523
在线: 296.3小时
虫号: 3579112
注册: 2014-12-05
性别: GG
专业: 声学

引用回帖:

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

请问您解决这个问题了吗？我也刚接触，完全不知道怎么入手…

赞一下

回复此楼

4楼2015-12-22 15:41:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖