| 查看: 780 | 回复: 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) } |
回复此楼» 猜你喜欢
梦想
已经有9人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
20