| 查看: 1397 | 回复: 0 | |||
[求助]
fluent中UDF设定边界条件，导入出错
|
|
最近在做一个fluent模拟，需要设定非稳态边界条件，写了一个UDF程序如下，导入fluent里面时，出现： Error: C:\Users\\305\313\346\257\Desktop\pexit2.c.c: line 18: parse error. Error: C:\Users\\305\313\346\257\Desktop\pexit2.c.c: line 19: invalid lvalue in assignment. 的错误，但是鄙人实在找不出错误之处，所以请教各位大神（由于是新虫，总共只有5个金币，对不起各位了！） UDF程序： #include"udf.h" #define r_1 20 #define r_2 15.4 #define a 0.77 #define w 418.879 #define H 18 #define t 3.2 #define m 1.2 #define p_s 587050 #define p_d 1980000 #define V_r pi*(pow(r_1,2)-pow(r_2,2))*H DEFINE_PROFILE(unsteady_pressure,thread,position) { face_t f; begin_f_loop(f,thread) { real t=RP_Get_Real("flow-time"  ;t=t-0.015*floor(t/0.015); if(t<0.002083 && t>0) F_PROFILE(f,thread,position)=587050; else if(t<=0.00845256 && t>=0.002083) real f_beta,h_beta,V_cbeta; f_beta=(1-pow(a,2))*418.879*t-0.5*pow((1-a),2)*sin(2.0*418.879*t)-pow(a,2)*asin((1.0/a-1)*sin(418.879*t))-a*(1-a)*sin(418.879*t)*pow((1-pow((1.0/a-1),2)*pow(sin(418.879*t),2)),0.5); h_beta=r_1*(1-(1-a)*cos(418.879*t)-sqrt(pow((1-a),2)*pow(cos(418.879*t),2)+2*a-1)); V_cbeta=V_r-0.5*H*(pow(r_1,2)*f_beta+h_beta*t); F_PROFILE(f,thread,position)=587050*pow((8932/V_cbeta),1.2); else if(t<=0.01375 && t>0.00854256) F_PROFILE(f,thread,position)=2144045.39-45844.5044*418.879*t; else if(t<=0.0141667 && t>0.01375) F_PROFILE(f,thread,position)=44548090.82-3103131857*t; else if(t<=0.015 && t>0.0146667) F_PROFILE(f,thread,position)=587050; } end_f_loop(f,thread) } DEFINE_PROFILE(unsteady_temperature,thread,position) { face_t f; begin_f_loop(f,thread) { real t=RP_Get_Real("flow-time" ;t=t-0.015*floor(t/0.015); if(t<0.002083 && t>0) F_PROFILE(f,thread,position)=308.15; else if(t<=0.00854256 && t>=0.002083) F_PROFILE(f,thread,position)=308.15*pow(587050*pow((8932/V_cbeta),1.2)/587050,0.16667); else if(t<=0.01375 && t>0.00854256) F_PROFILE(f,thread,position)=308.15*pow((2144045.39-45844.5044*418.879*t)/587050,0.16667); else if(t<=0.0141667 && t>0.01375) F_PROFILE(f,thread,position)=308.15*pow((44548090.82-3103131857*t)/587050,0.16667); else if(t<=0.015 && t>0.0141667) F_PROFILE(f,thread,position)=308.15; } end_f_loop(f,thread) } |
回复此楼» 猜你喜欢
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有10人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有6人回复
-
青基如何破局
已经有4人回复
-
基金不中，共勉
已经有8人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
Fluent 同一边界是否可以加载两个边界条件
已经有10人回复
- 求助fluent中的速度边界条件设置
已经有4人回复
- CFX或者Fluent各类出口边界条件在求解时具体是如何处理的？
已经有5人回复
- fluent 加载udf 失败 请高手不吝赐教
已经有8人回复
- fluent 加载udf出现这样错误 怎么办
已经有14人回复
- FLUENT中编写UDF时指针变量是否要初始化
已经有4人回复
- 【求助】fluent UDF设置边界条件
已经有13人回复
- fluent udf的应用
已经有20人回复
- Fluent中UDF的编译问题
已经有25人回复
- fluent 模拟结果 设为 另一个case的进口边界条件
已经有16人回复
- fluent出口边界条件和回流设置问题
已经有9人回复
- Fluent中UDF自定义了一个热流密度
已经有12人回复
- Fluent仿真电池冷却，遇到传热边界条件设置的问题，请教，附图片
已经有4人回复
- fluent中怎么能让边界默认为interior?
已经有8人回复
- udf求边界上的点出错
已经有3人回复
- fluent-udf自定义标量方程边界条件
已经有5人回复
- 边界条件的设置
已经有9人回复
- 好心人帮帮忙吧，为什么这个UDF导入fluent以后，在边界条件处没有显示~~~
已经有17人回复
- fluent中用udf编写三维速度边界条件出现错误，大家来指导一下啊
已经有9人回复
- 【求助】物料属性可以用UDF修改吗
已经有8人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
投票:
