24小时热门版块排行榜    

查看: 961  |  回复: 4

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;
}
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注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的回帖
简单回复
2013-11-08 23:48   回复  
mafeng0907(金币+1): 谢谢参与
2013-11-08 23:48   回复  
相关版块跳转 我要订阅楼主 mafeng0907 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 一志愿华中农业071010,总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
信息提示
请填处理意见