传热中边界为周性性热流量输入，如何设置呢？

水加火花

         在传热模块中，某一边界条件为热源/热沉q,其中q的输入是周期性的。0～100微秒为q,100～2000微秒为0；2000～2100微秒又为q,2100～4000微秒又为0，如此循环，总共求解时间为5秒。
          试图用如下方法来实现，但是发现行不通。边界是输入q*bo(t),其中bo(t)为调用的Matlab中的波函数，bo(t)=square(2*pi*500*t,5)，即周期为2000微秒，占空比为5%，但实际求解后发现，中在0-100微秒时间段内，q*bo(t)=1*1=q; 而在5s求解时间内，均为0。
           试图利用COMSOL中函数来定义，发现有分段函数可以使用，设置xstat,xend，f; 设置了好几组时间可以求解，但是我的求解总时间为5s ,至少得输入2500组才行啊。想请问一下高手，有没有更好的解决方法呢？
                 谢谢了！

水加火花

哦，上式中bo(t)=(square(2*pi*500*t,5)+1)/2.

Vivi_Song

可以考虑用逻辑表达式来表示，定义一个范围，在这个范围内，表达式的值为1，不在这个范围内，值为零，然后表达式乘以q，就是周期性热源啦。至于这个范围怎么确定，根据你的那个周期，可以写成一个关于n的表达式。n=1、2、3...

jinxiaoli_1981

将热流写成如下形式即可
q(t)=q*(t>(0.002*floor(t/0.002)))*(t<(0.0001+0.002*floor(t/0.002)))
注意，计算中时间单位是国际单位制s