| 6 | 1/1 | 返回列表 |
| 查看: 1060 | 回復(fù): 5 | |||
卑良新蟲 (小有名氣)
|
[求助]
求助幫忙改正一個FLUENT的UDF!急急急 已有1人參與
|
|
寫了一個UDF,描述閥芯在彈簧力和流體壓力作用下的運(yùn)動,可是執(zhí)行兩次后壓力總是變成很大的負(fù)值,進(jìn)而速度、位移都變得很大,網(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); } |
|
代碼沒問題,你需要檢查相關(guān)物性參數(shù)設(shè)置是否合理,比如彈簧的倔強(qiáng)系數(shù)是否合理,物體質(zhì)量是否合理等等 fluent培訓(xùn) |
新蟲 (小有名氣)
新蟲 (小有名氣)
新蟲 (小有名氣)
新蟲 (小有名氣)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項目,四級已過 +5 | 幸運(yùn)哩哩 2026-03-22 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +3 | Zhao anqi 2026-03-22 | 3/150 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +3 | vv迷 2026-03-22 | 3/150 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +3 | niko- 2026-03-22 | 3/150 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|