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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » udf中在单元格的所有面上的循环
51/1返回列表
查看: 941  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

mafeng0907

捐助贵宾 (初入文坛)

[交流] udf中在单元格的所有面上的循环

我编写的以下程序，在一个给定单元格的所有面上进行循环，得到在单元格上表面和下表面上的一个值sum3的和，但是在最后给sum3连接一个udm，该udm的值一直都是0，这是为什么呢？
real sum3= 0;
c_face_loop(cell,thread,n)
  {
      f = C_FACE(cell,thread,n);
      tf = C_FACE_THREAD(cell,thread,n);
         tc0 = THREAD_T0(tf);
         c0 = F_C0(f,tf);
         C_CENTROID(x0,c0,tc0);
         F_CENTROID(xf,c0,tc0);

if (BOUNDARY_FACE_THREAD_P(tf))
{
Tf = F_T(f,tf);
T_0 = C_T(c0,tc0);
sum3 += 0;
}
else
{
c1 = F_C1(f,tf);
tc1 = THREAD_T1(tf);

if(xf[1]=x0[1] && xf[0]>x0[0])
{
   if ((C_UDMI(c1,tc1,37)-sita_min)>0)
    {
     sum3 += (C_UDMI(c1,tc1,37)-sita_min)*M*g*gfactor/eps;
    }
   else
    {
     sum3 += 0;
    }
}
if (xf[1]=x0[1] && xf[0]      {
           if ((C_UDMI(c0,tc0,37)-sita_min)>0)
             {
                 sum3 += -(C_UDMI(c0,tc0,37)-sita_min)*M*g*gfactor/eps;
             }
          else
             {
                 sum3 += 0;
             }
     }
}
else
{
                 sum3 += 0;
}
}

if (NNULLP(THREAD_STORAGE(thread, SV_UDM_I)))
{
        C_UDMI(cell, thread, 34) = sum3;
}
回复此楼

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-06-06 21:39:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanghuandet

新虫 (初入文坛)
★
mafeng0907(金币+1): 谢谢参与
请问，我想确定非结构网格来流方向的面的有效截面积，该如何求？或者通过来流面的流量，关键是来流面的thread，如何就能使C_FACE(f,t,i)中的i是指的你想要的那个面呢？
赞一下 回复此楼
2楼2013-11-08 20:45:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vicly2013

铁虫 (小有名气)
★
mafeng0907(金币+1): 谢谢参与
进来学习一下
回复此楼
5楼2013-11-08 23:55:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
yingying15883楼
2013-11-08 23:48   回复  
mafeng0907(金币+1): 谢谢参与
yingying15884楼
2013-11-08 23:48   回复  
相关版块跳转 我要订阅楼主 mafeng0907 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭