ANSYSFLUENT计算pmv自定义函数代码.docVIP

ANSYSFLUENT计算pmv自定义函数代码 /* 使用时将此文档的扩展名改为.c 然后编译到fluent当中每个参数都有说明用户可以根据需要调整调用UDM需要首先设置Define User-Defined Memory... 设置为1即可此处只存储了一组数据*/ #include udf.h /* Author:wanghao SOUTH CHINA UNIVERSITY of TECHNOLOGY BEEL 2014*/ DEFINE_ON_DEMAND(on_demand_PMVcalc) { /*定义real数据类型fluent会转译为float数据类型*/ real t_wh; /*空气温度*/ real speed_u_wh; /* u方向速度分量*/ real speed_v_wh; /* v方向速度分量*/ real speed_w_wh; /* w方向速度分量*/ real RH_wh=40; /* 输入相对湿度水蒸气饱和时输入100 */ real pa_wh; /*物质质量分数C_YI(c,t,0) 后面通过相对湿度计算水蒸气分压力*/ real icl_wh=0.155; /*服装热阻在计算时单位不是clo 此处将0.3clo换算成0.05㎡K/W*/ real fcl_wh; /*穿衣服人体外表面与裸体人体表面积之比*/ real tcl_wh; /*服装外表面温度*/ real tcl1_wh; real tcl2_wh; real temperary1_wh; real temperary2_wh; real hc_wh; /*表面换热系数*/ real tr_wh; /*平均辐射温度可以后期自己定义*/ real M_wh=69.78; /*新陈代谢量*/ real W_wh=0.0; /*机械做功*/ real a_wh; /*heat loss by radiation refer to */ real b_wh; /*heat loss by convection refer to */ real c_wh; /*heat loss diff.through skin refer to */ real d_wh; /*heat loss by sweating(comfort) refer to */ real e_wh; /*latent respiration heat loss refer to */ real f_wh; /*dry respiration heat loss refer to */ real L_wh; real speed_wh; real PMV_wh; Domain *d; /* declare domain pointer since it is not passed as an argument to the DEFINE macro */ Thread *t; /* 声明Thread指针*/ cell_t c; d = Get_Domain(1); /*!!!!!这个要置于所有定义的变量之后Get th e domain using Fluent utility*/ thread_loop_c(t,d) /* Loop over all cell threads in the domain 因为loop是针对每个cell的所以针对cell的计算都要包含到loop里面*/ { begin_c_loop(c,t) { t_wh=C_T(c,t)-273.15; /*获取空气温度℃*/ speed_u_wh=C_U(c,t); /*获取u方向速度分量m/s */ speed_v_wh=C_V(c,t); /*获取v方向速度分量m/s */ speed_w_wh=C_W(c,t); /*获取w方向速度分量m/s */ tr_wh=t_wh; /*此处把平均辐射温度近似处理为空气温度℃可以后期自己定义*/ pa_wh=RH_wh*10*exp(16.6536-4030.183/(t_wh+235)); /*计算水蒸气分压力*/ speed_wh=sqrt(pow(speed_u_wh,2.0)+pow(speed_v_wh,2.0)+pow(speed_w_wh,2.0)); /*计算绝对空气相对速度*/ if (icl_wh0.078) fcl_wh=1.00+1.290*icl_wh; else fcl_wh=1.05+0.645*icl_wh; /*判断穿衣服人体外表面与裸体人体表面积之比*/ temperary1_wh=2.38*pow(fabs(tcl1_wh-t_wh-273),0.25); /*自然换热系数*/ temperary2_wh=12.