| 查看: 796 | 回复: 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) } |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有6人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有6人回复
-
现代”学阀”该如何界定
已经有14人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有4人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
我的奶奶
已经有3人回复