24小时热门版块排行榜    

查看: 703  |  回复: 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)
}
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 晓嘉嘉 的主题更新
信息提示
请填处理意见