| 查看: 809 | 回复: 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文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有11人回复
