| 查看: 761 | 回复: 0 | ||
[求助]
udf中正弦变化的位移边界加载问题
|
|
在udf中像位移边界为x(r,t)=0.2*e^（-r^6)*sin(2*pi*f*t）这种该怎么编写程序，其中x为垂直方向的分量，r为水平方向的分量，r属于[0,6]。自己尝试了一下感觉不太对，程序如下，就是这个r变量不知道怎么处理。还望老师和各位大牛指导一下。 #include"udf.h" #difine pi 3.14159265 DEFINE_CG_MOTION（zhengxian,dt,vel,omega,time,dtime) //定义边界的位移 { face_t f; //扫描的面 real_t=CURRENT_TIME; //定义当前时间 real f=33000 //振动频率 real dr=0.0001; //横向坐标距离微分 real r_intial=0.0; //横向初始位移 real r; //横向位移变量 if(r_intial<=6) { r_intial += dr; r=pow(r_intial,6); } if(r>6） return; begin_f_loop(f,thread,t) //循环开始 { vel(0)=0; vel(1)=0.2e(-10)*exp((200*r)*cos(2*pi*f*time)*2*pi*f) //垂直方向的速度 vel(2)=0; omega(0)=0.0; omega(1)=0.0; omega(2)=0.0; } end_f_loop（f,thread,f) } |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
100