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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 激光反射吸收热源模型！！！！代码重量爆料！！！！！求解答！！！！！！！！！
1171/3返回列表123下一页
查看: 5778  |  回复: 116
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

焊接网格划分

新虫 (初入文坛)

[交流] 激光反射吸收热源模型！！！！代码重量爆料！！！！！求解答！！！！！！！！！

代码有点错误，我现在是想用最小二乘法拟合自由表面的曲线，首先轮循所有自由表面网格，计算自由表面网格数量N，开辟数值存储自由表面上网格的重心（优良拟合三次曲线方程），gt,GS,max三个函数是用来求三次曲线方程y=a+bx+cx2+dx*x*x四个系数。然后把求出的四个系数保存到C_UDMI中，用于与离散的激光光线求交。代码编译之后有错误，可能思路上都已经与udf编程思想向偏离了，求各路大神指点，重谢！！！！！！！！！！！至少400金币！！！！呵呵
#include "udf.h"

#include "sg.h"

#include "sg_mphase.h"

#include "flow.h"

#include <math.h>

double gt(double x[],double y[]);
double GS(double a[4][4],double b[4],double x[4]);
double max(double array[4]);

DEFINE_ADJUST(store_gradient,domain)
{
  int phase_domain_index = 0;
  real xc[ND_ND],xd[ND_ND],x[4],n1,n2,thita_l,thita_g;
  Thread *t;
  Thread **pt;
  cell_t c;

  Domain *pDomain = DOMAIN_SUB_DOMAIN(domain,phase_domain_index);
  {
   Alloc_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
   Scalar_Reconstruction(pDomain, SV_VOF,-1,SV_VOF_RG,NULL);
   Scalar_Derivatives(pDomain,SV_VOF,-1,SV_VOF_G,SV_VOF_RG,Vof_Deriv_Accumulate);
  }
  mp_thread_loop_c(t,domain,pt)
{
  if (FLUID_THREAD_P(t))
  {
   int N=0;
   Thread *ppt = pt[phase_domain_index];

   begin_c_loop (c,t)
   {
          C_CENTROID(xd,c,t);
      if(C_VOF(c, ppt)>0.1&&C_VOF(c, ppt)<0.9&&fabs(xd[0]-0)>0.0005)
             N++;
   }
   end_c_loop (c,t)

   double a[N],b[N],c[N],d[N],e[N],f[N],g[N];

   int i;
   for(i=0; i<N; i++)
   {
       a=1;
   }

   int j=0;
   begin_c_loop (c,t)
   {
     C_CENTROID(xc,c,t);
     if(C_VOF(c, ppt)>0.1&&C_VOF(c, ppt)<0.9&&fabs(xc[0]-0)>0.0005)
     {
        b[j]=xc[0];
        e[j]=xc[1];
        c[j]=b[j]*b[j];
        d[j]=b[j]*b[j]*b[j];
                j++;  
     }
   }
   end_c_loop (c,t)
   f[0]=gt(a,a);f[1]=gt(a,b);f[2]=gt(a,c);f[3]=gt(a,d);f[4]=gt(b,b);f[5]=gt(b,c);f[6]=gt(b,d);
   f[7]=gt(c,c);f[8]=gt(c,d);f[9]=gt(d,d);g[0]=gt(a,e);g[1]=gt(b,e);g[2]=gt(c,e);g[3]=gt(d,e);
   double array[4][4]={{f[0],f[1],f[2],f[3]},{f[1],f[4],f[5],f[6]},{f[2],f[5],f[7],f[8]},{f[3],f[6],f[8],f[9]}};
   GS(array,g,x);
   C_UDMI(c,t,0) = x[0];
   C_UDMI(c,t,1) = x[1];
   C_UDMI(c,t,2) = x[2];
   C_UDMI(c,t,3) = x[3];
  }
}
  Free_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
}

double gt(double x[],double y[])
{   
    double sum=0;
        int i;
    for(i=0;i<N;i++)
       sum=x*y+sum;
    return sum;
}

double GS(double a[4][4],double b[4],double x[4])
{
    double c[4]={0};
    double x0[4]={0};
    int i,j,k;
    double sum=0;
    for(k=1;;k++)
    {
        for(i=0;i<4;i++)
        {
            for(j=0;j<4;j++)
            {
                sum=a[j]*x0[j]+sum;
            }
            x=x0+(b-sum)/a;
            c=fabs(x-x0);
            x0=x;
            sum=0;
        }
        r=max(c);
        if(r<0.0001)
        {
            break;
         }
    }
    return 0;
}

double max(double array[4])
{
    double a=array[0];
        int i;
    for(i=0;i<4;i++)
   {
        if(a<array)
           a=array;
    }
    return a;
}
回复此楼

» 猜你喜欢

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

查看全部散金贴
1楼 2013-12-02 10:11:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qingxiang827

铁虫 (知名作家)
★
焊接网格划分(金币+3): 谢谢参与
愿楼主心想事成
回复此楼
11楼2013-12-02 11:06:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gitme3388

铁虫 (知名作家)
★
焊接网格划分(金币+3): 谢谢参与
Good luck
回复此楼
12楼2013-12-02 11:06:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我爱喝奶

银虫 (小有名气)
★
焊接网格划分(金币+3): 谢谢参与
不太了解  祝福楼主吉星高照
赞一下 回复此楼
20楼2013-12-02 11:26:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crysmat

铁杆木虫 (知名作家)
★
焊接网格划分(金币+3): 谢谢参与
顶一个，祝楼主一切顺利
赞一下 回复此楼
22楼2013-12-02 11:37:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

焊接网格划分

新虫 (初入文坛)
虫友们！给大家提个方向去修改，可以用FORTRAN编写函数gt,GS,max,先预编译成.obj文件然后在ADJUST宏中调用。希望高手们多多指教！！！！！！！！
赞一下 回复此楼
37楼2013-12-03 10:27:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ttmany7788

铁虫 (知名作家)
★
焊接网格划分(金币+1): 谢谢参与
愿楼主心想事成
回复此楼
42楼2013-12-03 11:03:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freedom767

铁虫 (知名作家)
★
焊接网格划分(金币+1): 谢谢参与
愿楼主心想事成
回复此楼
43楼2013-12-03 11:04:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangwei19

荣誉版主 (文学泰斗)
★
焊接网格划分(金币+1): 谢谢参与
不会
回复此楼
48楼2013-12-03 11:23:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
湮没黑洞2楼
2013-12-02 10:21   回复  
焊接网格划分(金币+3): 谢谢参与
051031113楼
2013-12-02 10:25   回复  
焊接网格划分(金币+3): 谢谢参与
风雨飘然4楼
2013-12-02 10:29   回复  
焊接网格划分(金币+3): 谢谢参与
virologik5楼
2013-12-02 10:30   回复  
焊接网格划分(金币+3): 谢谢参与
haixiawu6楼
2013-12-02 10:30   回复  
焊接网格划分(金币+3): 谢谢参与
tzynew7楼
2013-12-02 10:33   回复  
焊接网格划分(金币+3): 谢谢参与
slz20088楼
2013-12-02 10:34   回复  
焊接网格划分(金币+3): 谢谢参与
yaoyao1459楼
2013-12-02 10:48   回复  
焊接网格划分(金币+3): 谢谢参与
wlxiao10楼
2013-12-02 10:57   回复  
焊接网格划分(金币+3): 谢谢参与
峰锋风逢烽13楼
2013-12-02 11:08   回复  
焊接网格划分(金币+3): 谢谢参与
doctorbang14楼
2013-12-02 11:14   回复  
焊接网格划分(金币+3): 谢谢参与
wullww15楼
2013-12-02 11:18   回复  
焊接网格划分(金币+3): 谢谢参与
aboluo19833116楼
2013-12-02 11:23   回复  
焊接网格划分(金币+3): 谢谢参与
pzchen17楼
2013-12-02 11:24   回复  
焊接网格划分(金币+3): 谢谢参与
xianzhi201018楼
2013-12-02 11:24   回复  
焊接网格划分(金币+3): 谢谢参与
fujiating19楼
2013-12-02 11:24   回复  
焊接网格划分(金币+3): 谢谢参与
wendytxy21楼
2013-12-02 11:33   回复  
焊接网格划分(金币+3): 谢谢参与
abczhang23楼
2013-12-02 11:41   回复  
焊接网格划分(金币+3): 谢谢参与
Spinage24楼
2013-12-02 11:41   回复  
焊接网格划分(金币+3): 谢谢参与
顶
chlq314101a25楼
2013-12-02 11:43   回复  
焊接网格划分(金币+3): 谢谢参与
宋全伟26楼
2013-12-02 11:44   回复  
焊接网格划分(金币+3): 谢谢参与
， [ 发自手机版 http://muchong.com/3g ]
会飞的乌龟--27楼
2013-12-02 11:44   回复  
焊接网格划分(金币+3): 谢谢参与
wander051528楼
2013-12-02 11:45   回复  
焊接网格划分(金币+3): 谢谢参与
yuzhichao29楼
2013-12-02 11:50   回复  
焊接网格划分(金币+3): 谢谢参与
祝福
Sammy230楼
2013-12-02 11:52   回复  
焊接网格划分(金币+3): 谢谢参与
doublehaohao31楼
2013-12-02 11:58   回复  
焊接网格划分(金币+3): 谢谢参与
尘轩无梦32楼
2013-12-02 11:59   回复  
焊接网格划分(金币+3): 谢谢参与
、 [ 发自手机版 http://muchong.com/3g ]
lgward33楼
2013-12-02 12:01   回复  
焊接网格划分(金币+3): 谢谢参与
huiminchlq34楼
2013-12-02 12:01   回复  
焊接网格划分(金币+3): 谢谢参与
xiejf35楼
2013-12-02 12:09   回复  
焊接网格划分(金币+1): 谢谢参与
silence136楼
2013-12-02 12:25   回复  
ts_chenchao38楼
2013-12-03 11:00   回复  
焊接网格划分(金币+1): 谢谢参与
天涯201239楼
2013-12-03 11:02   回复  
焊接网格划分(金币+1): 谢谢参与
neu23440楼
2013-12-03 11:02   回复  
焊接网格划分(金币+1): 谢谢参与
miaojiabing41楼
2013-12-03 11:03   回复  
焊接网格划分(金币+1): 谢谢参与
gyfwon44楼
2013-12-03 11:09   回复  
焊接网格划分(金币+1): 谢谢参与
gander45楼
2013-12-03 11:09   回复  
焊接网格划分(金币+1): 谢谢参与
v
liu12333846楼
2013-12-03 11:11   回复  
焊接网格划分(金币+1): 谢谢参与
lefantian12247楼
2013-12-03 11:18   回复  
焊接网格划分(金币+1): 谢谢参与
zhoupeng8749楼
2013-12-03 11:32   回复  
焊接网格划分(金币+1): 谢谢参与
·
base0803030850楼
2013-12-03 11:34   回复  
焊接网格划分(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 焊接网格划分 的主题更新
1171/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +8 QUjhNVAcOSff 2026-09-13 9/450 2026-09-17 22:05 by cuPZTDXS3VOd
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 5BDX0d0WFp7t 2026-09-15 4/200 2026-09-17 16:29 by 0ngsZVSCTS7z
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-17 16:03 by 0ngsZVSCTS7z
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-17 13:29 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 yaCh4X3wd045 2026-09-14 4/200 2026-09-17 13:17 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 13:15 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-17 12:39 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 mibUvS8DDCwf 2026-09-15 4/200 2026-09-17 06:26 by k06BKNblrGRH
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-17 05:14 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-17 04:16 by k06BKNblrGRH
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-17 03:52 by k06BKNblrGRH
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-13 6/300 2026-09-17 02:14 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-16 19:15 by rRvxEj3zorCy
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +7 QUjhNVAcOSff 2026-09-13 7/350 2026-09-16 14:36 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
信息提示
关闭
请填处理意见
关闭