24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

whcjordan123

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 2100.1
散金: 20
帖子: 75
在线: 35.7小时
虫号: 1406140
注册: 2011-09-18
性别: GG
专业: 机械动力学

[交流] ICEM网格划分问题已有9人参与

小弟最近在做高速列车的外流场仿真，用的ICEM划分网格，现在有点迷茫的是，我应该对计算域（包括列车在内）整体进行块划分，不断切割得到结构六面体网格，还是对计算域和列车分别划分不同类型的网格，但又怕最后两者的兼容问题。求高手指点。
另外，如果分开划分，是不是需要设定不同的part或体啊？

红色为高速列车

回复此楼

» 猜你喜欢

请问社科规划办官网上的“项目查询模块”与“项目申报模块”的账号密码一样吗？已经有3人回复
用NBS溴化时溶剂除了CCl4外还有什么已经有9人回复
40+中年人失业后，还能找到好工作吗？已经有15人回复
二审4个多月了还是Awaiting Reviewer Scores 已经有5人回复
申请科研项目纯靠实力能成功吗？已经有12人回复
极限拉扯，如何抉择已经有15人回复

有容乃大！

1楼 2012-04-22 20:09:57

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

whcjordan123

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 2100.1
散金: 20
帖子: 75
在线: 35.7小时
虫号: 1406140
注册: 2011-09-18
性别: GG
专业: 机械动力学

引用回帖:

5楼: Originally posted by zxbok1979 at 2012-04-23 11:44:38
主要看你想做到什么程度，是否一定要结构网格？是否考虑转向架、连接部位、受电弓等？还是简化处理，不考虑上述部位，除了车头其他地方很规则？如果考虑上述三个部位，那么很不幸的告诉你，很难；单一个转向架，一个 ...

大侠，能不能帮我看看我的动网格udf哪里出错了，就是高速列车驶过声屏障.万分感谢！
#include “udf.h”
#if !RP_NODE
#define UDF_FILENAME “udf_loc_velo”
/* read current location and velocity from file */
Static void
Read_loc_velo_file(real*loc,real*velo)
{
  FILE*fp=fopen(UDF_FILENAME,”r”);
  if(fp !=NULL)
{
  float read_loc,read_velo;
  fscanf (fp,”%e %e”,&read_loc,&read_velo);
fclose(fp);
*loc=(real) read_loc;
*velo=(real) read_velo;
}
else
{
  *loc=0.0;
  *velo=0.0;
}
} /* write current location and velocity in file*/
static void
write_loc_velo_file(real loc,real velo)
{
FILE*fp=fopen(UDF_FILENAME,”w”);
  if(fp!=NULL)
{
  fprintf(fp,”%e %e”,loc,velo);
fclose(fp);
}
else
Message(“\n Warning:cannot write %s file”,UDF_FILENAME);
}
#endif /*!RP_NODE*/

DEFINE_CG_MOTION(mvz250z,dt,cg_vel,cg_omega,time,dtime)
{
real f;//=RP_Get_Real(“motion-par/freq_s”);
real a;//=RP_Get_Real(“motion-par/amplitude”);
real omega; real vel;
/*real amplitude=.20/* in inches */
/* real conv =0.0254;*/
omega=2.0*M_PI*f;
/*vel=amplitude*conv*omega*sin(omega*time);*/
Vel=250.0/3.6;
cg_vel[0]=0.0;/*x-velocity*/
cg_vel[1]=0.0;
cg_vel[2]=vel;
NV_S(cg_omega,=,0.0);/* no angular motion*/
}
DEFINE_CG_MOTION(muf250z,dt,cg_vel,cg_omega,time,dtime)
{
real f;//=RP_Get_Real(“motion-par/freq_s”);
real a;//=RP_Get_Real(“motion-par/amplitude”);
real omega;
real vel;
/*real amplitude =.20;/* in inches*/
/*real conv = 0.0254;*/
omega=2.0*M_PI*f;
/*vel=amplitude*conv*omega*sin(omega*time);*/
Vel=-250.0/3.6;
cg_vel[0]=0.0;/*x_velocity*/
cg_vel[1]=0.0;
cg_vel[2]=vel;
NV_S(cg_omege,=,0.0);/* no angular motion*/
}