| 5 | 1/1 | 返回列表 |
| 查看: 1067 | 回復: 5 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
求助幫忙改正一個FLUENT的UDF!急急急 已有1人參與
|
||
|
寫了一個UDF,描述閥芯在彈簧力和流體壓力作用下的運動,可是執(zhí)行兩次后壓力總是變成很大的負值,進而速度、位移都變得很大,網(wǎng)格更新就不成功了。能否有人幫忙看一下? 附UDF源碼: #include <stdio.h> #include "udf.h" #define UDF_FILENAME "udf_loc_velo.txt" # define K_SPRING 31400 /* read current location and velocity from file */ DEFINE_CG_MOTION(spool, dt, vel, omega, time, dtime) { Thread *t = DT_THREAD (dt); face_t f; static real force; static real loc; static real velo; static real init_disp; static real s_force ; static real dv; static real total_force; static FILE *fp; /* Open data files */ 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; } /* compute force on piston wall */ force = 0.0; begin_f_loop (f, t) { real *AA; AA = F_AREA_CACHE (f, t); force += F_P (f, t) * AA[2]; } end_f_loop (f, t) /* add in spring force */ init_disp = 0.01; s_force = K_SPRING * (loc + init_disp); total_force = force - s_force-346-10*velo; /* compute change in velocity */ dv = dtime * total_force / 0.554; velo += dv; loc += velo * dtime; vel[2] = velo; write_loc_velo_file(loc,velo); Message ("\nUDF spool: time = %f, z_vel = %f, force = %f,total_f = %f, loc(m)= %f\n", time, velo, force,total_force, loc); } /* write current location and velocity in file*/ 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); } |
新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 283求調劑(080500) +7 | A child 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 311求調劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
[考研] 【求調劑】085601材料工程專碩 | 總分272 | +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 一志愿南昌大學324求調劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
[考研] 277跪求調劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|
[考研] 求調劑 +9 | 張zz111 2026-03-27 | 10/500 |
|
|
[考研] 297求調劑 +11 | 田洪有 2026-03-26 | 11/550 |
|
|
[考研] 375求調劑 +4 | 雨夏整夜 2026-03-29 | 4/200 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 0703本科鄭州大學求調劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
[考研] 086000調劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 求調劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
[考研] 336材料求調劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 321求調劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 中國科學院深圳先進技術研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 環(huán)境專碩324分求調劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 材料專碩331求調劑 +4 | 鮮當牛 2026-03-24 | 4/200 |
|
|
[考研] 277分求調劑,跨調材料 +3 | 考研調劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 接收2026碩士調劑(學碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|