我来发个问题

rp_no

问题是这样的，用weak来作一个普通的扩散方程： n*test+D*(test(ux)*ux+test(uy)*uy)

这个简单

然后给一个函数max(x,y),希望在n<max时候按照扩散方程计算，但n不能超过max的值。

按原来的思路，需要在公式中加入判断，n是否大于max，但这个表达式设计变量n本身，无法求解。

有什么好的思路么

sky_walker

可以用step()函数进行判断

rp_no

sky_walker 发表于 2014-10-13 19:38
可以用step()函数进行判断

能具体说说么

rp_no

sky_walker 发表于 2014-10-13 19:38
可以用step()函数进行判断

感觉应该是不能在方程中使用dependent variable本身。。所以貌似这问题无解。

sky_walker

step()函数的路径为Definition>Function<Step, step()函数中可以使用dependent variable本身啊，我经常这样做，用step()来定义分段函数。

rp_no

sky_walker 发表于 2014-10-14 10:09
step()函数的路径为Definition>Function

是不是STEP函数的算法，可以保证比较好的连续性和高阶导数连续性？

sky_walker

STEP函数可以设置光滑度(smoothing),从而保证较好的连续性

rp_no

sky_walker 发表于 2014-10-17 16:52
STEP函数可以设置光滑度(smoothing),从而保证较好的连续性

昨天试了一下，确实。但是如果过度区间大的话，阶跃性就不好了。我想用这个来做一个状态判断，接触的，所以要求阶跃很大。

后来也查了一下CONTACT的WEAK，居然直接就是用的逻辑判断来做的接触状态判断。
当然我也怀疑是否COMSOL解这类问题比较弱，又或者是我对SOLVER的设置不对

之前做abaqus分析复杂的大面积接触都是很稳定的，而COMSOL稍微复杂一点的接触收敛性感觉就很不好了

rp_no

sky_walker 发表于 2014-10-17 16:52
STEP函数可以设置光滑度(smoothing),从而保证较好的连续性

不好意思，看了下1楼，这是2个问题。。

但是都设计这种跳变函数的问题。。。顶楼的问题还设计一个是否可以在公式中使用dependent variable自身做判断的问题。

sky_walker

做接触还是用abaqus吧，comsol只能做一些简单的接触问题

rp_no

sky_walker 发表于 2014-10-18 13:34
做接触还是用abaqus吧，comsol只能做一些简单的接触问题

1楼的问题基本解决了
不用稳态，用瞬态，就比较好的能计算step了

接触那个其实就是一个状态非线性。我有个计算必须用COMSOL来做。略有点头疼。别的都还好。。。