应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » COMSOL Multiphysics & HYPERWORKS » Comsol自定义场变量，只要曾经温度T>Ts，对应场变量元素为1，一旦为1，则永久保持
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 3103  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

scottxhx

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

[交流] Comsol自定义场变量，只要曾经温度T>Ts，对应场变量元素为1，一旦为1，则永久保持 已有3人参与

问题：comsol能否自定义一个场变量，来标记每个网格点温度T曾经是否大于Ts，一旦某个网格点温度曾经大于Ts，该网格状态量为1，即便温度降下来了，状态量不变。

这个问题非常具有代表性，其归根结底是要定义场变量f(x,y,z,t)，这个场变量依附于每一个网格，可以初始化，可被用户自定义操作。我目前发现，comsol并不能不依附任何物理场模块自定义这样的变量，退而求其次的做法是采用comsol的数学模块里的域偏微分/域常微分方程/不可逆转变/生物传热模块，这些做法逻辑是能行得通，但是求解非常不准确。采用偏微分或者常微分模块，无论是通过系数将模块数学方程u=f，还是du/dt=f这么简单的方式，其中f通过逻辑控制，即便这样，我发现comsol求解出来的u误差很大，连u=f，这里置f=1或者0，都不能严格控制计算的u=1或者0；du/dt=f的方式求解的误差更离谱，我一直定义f为正，但是u求解出来的还有一点点负值，这样就无法通过u来判断它曾经的状态。不可逆转变/生物传热一样有这样的问题，我觉得这个东西应该很简单，而且计算也不困难，不知道大家有没有什么好的方案。

还要说的是，Matlab函数我也试过，本来想把这样的状态量在matlab里定义为全局场变量，但是发现comsol送过来的网格坐标是分批的，打个比方，模型有50万网格，comsol会分好几拨送过来，也就是Matlab函数在一次时间迭代开始前会被调用好几次，而且送过来的网格批次的先后顺序也不得而知，这就无法对齐场变量元素的位置。
回复此楼

» 猜你喜欢
1楼 2018-11-13 14:22:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scottxhx

兑换贵宾

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
7楼: Originally posted by 满手灰尘 at 2018-11-16 13:29:46
我之前也是用pde做过类似的，但没有发现你说的计算结果不准确的现象，你的结果是有多离谱，贴张图来看看

我在这里把系数PDE设置成u=f, 在计算过程中f只会等于0或者1，结果算完，u从-0.4~1.x多都有值，这里和状态转变速度有关系，comsol求解器对于转变速度很快的方程求解，会出现很大的人为光滑。
赞一下 回复此楼
8楼2018-11-21 16:46:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

landl3017

超级版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
用数学模块下的事件试试？官网博客有一个恒温箱的例子，温度降到一定值就会加热

发自小木虫Android客户端
赞一下 回复此楼
2楼2018-11-13 15:44:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scottxhx

实习版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
2楼: Originally posted by landl3017 at 2018-11-13 15:44:04
用数学模块下的事件试试？官网博客有一个恒温箱的例子，温度降到一定值就会加热

我试了一下这个事件模块怎么只能调用全局变量，比如我想写为温度的函数，运行的时候提示变量未定义
赞一下 回复此楼
3楼2018-11-13 17:56:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

landl3017

实习版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by scottxhx at 2018-11-13 17:56:14
我试了一下这个事件模块怎么只能调用全局变量，比如我想写为温度的函数，运行的时候提示变量未定义...

没用过事件，只是看过帮助。事件有显示事件和隐式事件。隐式事件还要配套另一个项目使用。

发自小木虫Android客户端
赞一下 回复此楼
4楼2018-11-13 19:35:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 335求调剂 +3 321* 2026-03-31 4/200 2026-04-01 00:00 by 321*
[考研] 080200学硕，机械工程专业277分，求带走！ +4 瓶子PZ 2026-03-31 4/200 2026-03-31 20:16 by vgtyfty
[考研] 304求调剂 +8 素年祭语 2026-03-31 11/550 2026-03-31 18:20 by 无际的草原
[考研] 生物学296求调剂 +8 汤圆包 2026-03-29 12/600 2026-03-31 17:05 by 18828373951
[考研] 329求调剂，一志愿西北工业大学，材料工程（085601） +6 小小机灵虫 2026-03-29 12/600 2026-03-31 16:58 by 记事本2026
[考研] 求调剂 +8 11ggg 2026-03-30 8/400 2026-03-31 13:56 by nanaliuyun
[考研] 270求调剂 +3 小杰pp 2026-03-31 4/200 2026-03-31 12:59 by wxiongid
[考博] 材料专业申博 +5 杜雨婷dyt 2026-03-29 5/250 2026-03-31 11:19 by oooqiao
[考研] 085600，专业课化工原理，320分求调剂 +6 大馋小子 2026-03-29 6/300 2026-03-31 10:03 by 氯化亚硝酰
[考研] 085601一志愿中山大学深圳材料工程330求调剂 +5 pipiver 2026-03-30 5/250 2026-03-31 07:37 by JourneyLucky
[考研] 281求调剂 +5 亚克西good 2026-03-26 7/350 2026-03-30 20:42 by dophin1985
[考研] 284求调剂 +14 junqihahaha 2026-03-26 15/750 2026-03-30 14:12 by 探123
[考研] 329求调剂 +8 星野? 2026-03-26 8/400 2026-03-30 13:41 by chemdavid
[考研] 考研调剂 +5 Sanmu-124 2026-03-26 5/250 2026-03-30 13:36 by chemdavid
[考研] 085600，专业课化工原理，321分求调剂 +5 大馋小子 2026-03-28 5/250 2026-03-29 08:56 by qingfeng258
[考研] 298求调剂 +4 种圣赐 2026-03-28 4/200 2026-03-29 08:42 by q1092522407
[考研] 305求调剂 +8 RuiFairyrui 2026-03-28 8/400 2026-03-29 08:22 by fmesaito
[考研] 275求调剂 +10 jjjjjjjjjjl 2026-03-27 10/500 2026-03-27 23:47 by barnett0632
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[考研] 一志愿吉大071010，316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
信息提示
关闭
请填处理意见
关闭