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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » COMSOL Multiphysics & HYPERWORKS » Comsol自定义场变量，只要曾经温度T>Ts，对应场变量元素为1，一旦为1，则永久保持
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 3102  |  回复: 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

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
5楼: Originally posted by bfamq at 2018-11-13 21:27:34
遇到这个问题的人不少，comsol中的变量都是基于全部节点的，条件改变，所有节点值都会变。除非可以调用每个节点才行。

Comsol调用每个节点不是问题，主要问题在于无法获知节点编号。就像我前面说的，我尝试过一种思路，就是在matlab中定义与comsol网格节点对应的场变量，但是comsol调用matlab函数的时候，传给matlab的节点是分好几拨的，而且无法获取某一波网格节点编号。
赞一下 回复此楼
6楼2018-11-15 11:10:40
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 271求调剂 +8 勒布朗@ 2026-03-31 10/500 2026-03-31 22:53 by wwytracy
[考研] 339求调剂 +3 zjjkt 2026-03-31 3/150 2026-03-31 20:34 by 83503孙老师
[考研] 一志愿北京科技大学材料学硕328分求调剂 +4 1段时间 2026-03-31 5/250 2026-03-31 19:56 by wxiongid
[考研] 375求调剂 +7 雨夏整夜 2026-03-29 7/350 2026-03-31 18:52 by xhai2011
[考研] 材料调剂 +9 Eujd1 2026-03-31 10/500 2026-03-31 18:41 by JourneyLucky
[考研] 一志愿南昌大学324求调剂 +6 hanamiko 2026-03-29 6/300 2026-03-31 16:35 by hypershenger
[考研] 354求调剂 +3 lxb598 2026-03-31 4/200 2026-03-31 13:42 by sophie2180
[考研] 调剂求院校招收 +7 鹤鲸鸽 2026-03-28 7/350 2026-03-31 11:21 by oooqiao
[考研] 吉大生物学326分求调剂 +3 sunnyupup 2026-03-31 3/150 2026-03-31 09:28 by longlotian
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
[考研] 11408总分309，一志愿东南大学求调剂，不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
[考研] 292求调剂 +13 是妍子也是研子 2026-03-30 13/650 2026-03-30 18:01 by 小徐0109
[考研] 一志愿中南大学化学0703总分337求调剂 +6 niko- 2026-03-27 6/300 2026-03-30 10:25 by herarysara
[考研] 070305高分子化学与物理 304分求调剂 +12 c297914 2026-03-28 12/600 2026-03-29 16:21 by Serene1974
[考研] 321求调剂 +7 璞玉~~ 2026-03-25 8/400 2026-03-29 06:41 by 544594351
[考研] 316求调剂 +7 江辞666 2026-03-26 7/350 2026-03-28 21:28 by sanrepian
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 348求调剂 +4 小懒虫不懒了 2026-03-27 5/250 2026-03-27 12:47 by 果果妈咪
[考研] 321求调剂 +6 wasdssaa 2026-03-26 6/300 2026-03-26 20:57 by sanrepian
信息提示
关闭
请填处理意见
关闭