| 查看: 720 | 回复: 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) } |
回复此楼» 猜你喜欢
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有7人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有9人回复
-
今年春晚有几个节目很不错，点赞！
已经有10人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有10人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复