分类目录归档:技术博客

针对一些技术的展示

3min速解COMSOL 中湍流模型分类

在 COMSOL 6.3 的 CFD 模块中,湍流模拟提供了多种建模方法,包括 RANS-EVM(涡黏模型)、RANS-RSM(雷诺应力模型)、大涡模拟(LES) 和 分离涡模拟(DES)。这些方法在计算精度、资源消耗和适用场景上各有侧重。下面我们将逐一解析它们的原理、优缺点及典型应用,帮助你根据实际问题做出明智选择。

🔧 一、RANS-EVM:涡黏模型 —— 快速高效的工程利器

✅ 原理: RANS(Reynolds-Averaged Navier-Stokes)通过对 Navier-Stokes 方程进行时间平均,将湍流效应建模为附加的“湍流黏度”。

EVM(Eddy Viscosity Models) 假设雷诺应力与平均应变率成正比,用湍流黏度来封闭方程。

📌 常见模型

  • k-ε(标准、Realizable)
  • k-ω(包括 SST)
  • Spalart-Allmaras
  • v2-f、低雷诺数 k-ε 等

✅ 优点:

  • 计算速度快,资源消耗低
  • 对大多数稳态或弱非稳态工程流动足够准确
  • 易于设置,适合初步设计或参数研究

❌ 缺点:

  • 无法准确模拟强旋流、分离流、各向异性湍流
  • 对复杂几何或曲率大的流动预测能力有限

📍 适用场景

  • 管道流、换热器、风力机叶片、建筑风载等常规工程问题

🔍 二、RANS-RSM:雷诺应力模型 —— 更精细的湍流描述

✅ 原理: RSM(Reynolds Stress Models) 不再假设雷诺应力与应变率成正比,而是直接求解雷诺应力张量的输运方程,能更好地反映湍流的各向异性。

📌 COMSOL 6.3 新增模型:

  • Wilcox R-ω
  • SSG–LRR(Speziale–Sarkar–Gatski / Launder–Reece–Rodi)

✅ 优点:

  • 能准确模拟强旋涡、二次流、曲率影响大的流动
  • 比 EVM 更适合复杂三维湍流结构

❌ 缺点:

  • 计算量大,收敛性差
  • 对网格质量和初值敏感

📍 适用场景:

  • 涡轮机械、旋转管道、旋风分离器、强旋燃烧室等

🌀 三、LES:大涡模拟 —— 追求瞬态结构的精度之选

✅ 原理: LES(Large Eddy Simulation) 使用空间滤波将湍流分为大尺度(解析)和小尺度(建模)。大涡直接模拟,小涡用亚网格模型(如Smagorinsky)建模。

✅ 优点:

  • 能捕捉瞬态大尺度涡结构,精度远高于 RANS
  • 适用于研究湍流机理、噪声、混合等问题

❌ 缺点:

  • 网格要求极高,计算量大
  • 时间步长小,模拟时间长
  • 对初边值条件敏感

📍 适用场景:

  • 湍流噪声预测、燃烧模拟、污染物扩散、流体-结构耦合等

🔁 四、DES:分离涡模拟 —— RANS 与 LES 的混合策略

✅ 原理:

DES(Detached Eddy Simulation) 是一种混合模型:

  • 在边界层使用 RANS(节省资源)
  • 在分离区自动切换为 LES(提高精度)

✅ 优点:

  • 在分离流动、非定常尾涡等区域比 RANS 更准确
  • 相比纯 LES,计算量更小,适合工程应用

❌ 缺点:

  • 网格设计复杂,需人工经验
  • 模型切换区域可能引入误差(如“灰区”问题)

📍 适用场景:

  • 汽车绕流、飞行器尾涡、建筑物风载、流体诱导振动等

🧭 总结:如何选择湍流模型?

模型类型 精度 计算成本 适用流动类型
RANS-EVM ★★☆☆☆ ★☆☆☆☆ 稳态/弱非稳态工程流动
RANS-RSM ★★★☆☆ ★★☆☆☆ 强旋、各向异性湍流
LES ★★★★☆ ★★★★☆ 瞬态大涡结构研究
DES ★★★☆☆ ★★★☆☆ 分离流、非定常尾涡

✅ 建议:

  • 做初步设计或工业仿真?用 RANS-EVM
  • 有强旋涡或复杂曲率?选 RANS-RSM
  • 想研究瞬态涡结构或噪声?上 LES
  • 想兼顾精度与效率,模拟分离流?用 DES

📚 参考资料 : COMSOL 中国. (2025). CFD 模块更新 – COMSOL® 6.3 发布亮点. https://cn.comsol.com/release/6.3/cfd-module

: COMSOL 中国. CFD 建模和仿真软件. https://cn.comsol.com/cfd-module

: 知乎专栏. (2025). RANS vs DES vs LES vs DNS:湍流模拟的区别及选择指南. https://zhuanlan.zhihu.com/p/18028581363

如需进一步了解如何在 COMSOL 中设置这些模型,或获取教学案例,欢迎继续提问!

吃透 COMSOL 后处理(五):破解“不连续”之谜——薄层结构后处理的终极指南

开篇引言:

“明明物理上连续,云图却出现‘断崖式’跳变;薄壁两侧温度本该平滑,颜色条却硬生生劈成两半?薄层、壳、膜、内部壁……这些降维建模的‘黑科技’,后处理时却常让人怀疑人生。今天,我们迎战后处理领域的终极挑战:如何正确解读薄层结构里的‘假不连续’,让 up、down、side 三大算子替你拆穿一切陷阱。”

正文核心:

1. 为什么‘薄层’会出现‘假不连续’?——先救认知
  • 建模真相
    • 3D 实体壁厚 = 几何真实网格;
    • 薄层 = 2D 边界 + 数学厚度(无体积网格)。
  • 后处理副作用
    • 边界只存一份解,却与两侧域相邻;
    • 表面图默认“平均”两侧结果 → 相邻边界间无梯度 → 视觉跳变;
    • 域图(切面、体图)只能选单侧 → 切换边界时数值突变 ≠ 错误。
2. up & down 算子:只看‘上/下’侧的真相
  • 语法
    • up(solid.mises) ← 外表面(+z 法向)
    • down(solid.mises) ← 内表面(-z 法向)

一、实操案例:

1.层压复合壳(laminated shell thermal expansion)

模型亮点

  • 6 层不同纤维角 → 弯曲-热耦合
  • 壳物理场 + 多层材料 → 数学厚度 2 mm(无实体网格)
  • 已求解 → 直接后处理

实操

结果 > 三维绘图组 > 表面 → 表达式 lshell.misesGp

外观:波浪形剖面连续云图

陷阱:壳默认只绘“参考面”,上下表面被隐藏!

  • 复制表面图 → 重命名 Outer → 表达式改为 up(solid.mises) ← 外表面(+z 法向)
  • 再复制 → Inner → 表达式 down(solid.mises) ← 内表面(-z 法向)
  • 色标同步 0 – 30 MPa

一句话写进论文
“图中 可见层间弯曲导致外纤维应力较内纤维高 ,壳模型一次性捕获梯度。”

3. side 算子:多域交汇点的‘定向显微镜’
  • 语法side(dom, 表达式) → 在 n-1 维实体上,显式提取与指定域 dom 相邻侧的结果。
  • 优势
    • 可跨维度:边界、边、点皆可用;
    • 避开‘up/down 仅升/降维’限制。

多层剖面:厚度方向“假几何”放大 20×

操作:结果 > 三维绘图组 > 多层材料切面 → 选择“最外层” & “最内层”

  • 高度表达式:缩放因子 20 → 2 mm 变 40 mm 可视
  • 颜色:应力 solid.mises

输出:波浪形剖面,外红内蓝 → 审稿人一眼 get 梯度方向。

2.壳-梁连接(shell_beam_connection)

路径:App 库 > 结构力学 > Shell > 壳和梁连接

1. 模型亮点

  • 3 mm 薄壁壳 + 梁腿 → 角焊缝弯曲载荷
  • 已求解 → 直接后处理

2. Step-1 上下表面应力“台阶”线

操作:

  • 三维截线 → 沿壳壁母线 (0,0,0) → (0,0,0.3) m
  • 线图: – 数据1 up(solid.mises) 管半径 0.25 mm 红色 – 数据2 down(solid.mises) 管半径 0.25 mm 蓝色
  • 结果:一条“红蓝双管”→ 台阶高 15 MPa → 弯曲主导

3. Step-2 批量筛选“壳面朝向”

需求:只显示与载荷方向夹角 < 90° 的壳面

  • 定义 > 选择 > 壳边界 → 角度选择:与 (0,1,0) 点积 > 0
  • 表面图 → 表达式 up(solid.mises) → 选择上述集合

二、对比验证:壳 vs 实体(三分钟搞定)

  • 同一几何 → 复制研究 → 改用“固体力学” + 3 层实体网格(厚度方向)
  • 截面母线应力:实体最大 76 MPa 壳 up 74 MPa
  • 误差 2.7 % → 计算时间从 8 min → 1 min

结论句: “在 2.7 % 误差范围内,壳模型节省 87 % 计算时间,适用于初始设计迭代。”


三、Shell 后处理 Checklist(官方模型版)

步骤自查问题在层压壳/壳-梁模型中的操作
想看哪一侧?up/down 显式指定外/内表面
双侧差多少?up(σ) – down(σ) 直接量化梯度
厚度方向可视化?多层材料切面 + 高度表达式放大
批量筛选朝向?角度选择 + side(dom,expr) 一键过滤
壳结果可信吗?与实体网格截面应力对比误差

结尾总结:

“薄层结构的后处理,陷阱不在物理,而在‘维度降维’后的解读方式。掌握 up、down、side 三大算子,你就拥有了‘定向透镜’:想看哪一侧,就看哪一侧;想证明连续,就能拿出单侧数据;想展示跳变,也能一条线画出台阶。至此,五篇系列全部完结——从数据体检、图像美颜、深度剖析、专业制图,到破解不连续,你已集齐 COMSOL 后处理的完整技能树。愿你下一次汇报,不再只说‘看起来差不多’,而是指着图自信地说:‘这里,就是设计优化的关键点!’”

“吃透 COMSOL 后处理”系列:从图像到价值的有限元仿真进阶之路

COMSOL Multiphysics 的“后处理”往往被当成锦上添花的一步,实则却是决定仿真价值能否真正落地的“最后一公里”。初学COMSOL的用户往往在“后处理”环节陷入“数据堆砌”的困境:颜色泛滥、图窗杂乱、关键物理量彼此掩盖,最终难以用一张图讲清故事。针对这一痛点,为了让“数据”顺利升格为“信息”,再沉淀为“洞察”,这篇博文将给出一条可复制的“零代码”后处理流水线,详尽地介绍 COMSOL 中的各种后处理形式,包括:

  • 基本绘图与表达式: 涵盖创建自定义表达式、使用派生值(如全局计算、平均值)验证模型、绘制网格图来评估网格质量等。
  • 高级可视化技巧: 包括使用变形、高度表达式和缩放来增强可视化效果,以及利用截面、截线、合并数据集和过滤器来聚焦和分析特定区域的结果。
  • 专用绘图类型: 介绍了针对特定物理场的绘图,如射频(RF)分析的远场图和粒子追踪图。
  • 降维后处理: 针对薄层结构,详细讲解了如何处理和解释在边界(二维)上绘制体(三维)结果时可能出现的“不连续性”,并介绍了updownside 算子的使用。
  • 视图、相机与灯光控制: 说明了如何通过调整视角、相机位置和多种光源(定向光、点光源、聚光灯)来创建更具表现力和清晰度的图像。
  • 动画与结果导出: 简要提及了如何创建动画并导出结果。

分步规划:

  • 第一篇:基础篇——《数据“体检”与表达:如何确保你的仿真结果值得信赖?》
    • 核心内容:聚焦仿真结果的基础验证与核心绘图表达。
    • 关键技巧:自定义表达式、派生值(全局计算)、网格质量评估。
  • 第二篇:技巧篇——《让图像“说话”:提升仿真图表现力的核心技巧》
    • 核心内容:介绍增强可视化效果的通用技巧。
    • 关键技巧:变形、高度表达式、缩放、视图与相机控制、灯光与场景设置。
  • 第三篇:高级篇——《洞察与聚焦:如何深度剖析你的仿真数据?》
    • 核心内容:讲解如何对数据进行降维、切片和聚焦分析。
    • 关键技巧:截面与截线、合并数据集、过滤器。
  • 第四篇:专业篇——《专业图表与特殊应用:满足你的“小众”需求》
    • 核心内容:介绍针对特定物理场的专用绘图类型。
    • 关键技巧:远场图(RF/声学)、粒子追踪图(庞加莱图、相图)。
  • 第五篇:挑战篇——《破解“不连续”之谜:薄层结构后处理的终极指南》
    • 核心内容:深度解析薄层结构中的降维后处理问题。
    • 关键技巧:理解边界与域图差异、掌握updownside 算子。

第一篇:数据“体检”与表达

标题:

吃透 COMSOL 后处理(一):数据“体检”与表达——如何确保你的仿真结果值得信赖?

开篇引言:

“当你的仿真结果新鲜出炉,你是否也曾感到无从下手?面对五彩斑斓的云图,除了‘看起来不错’,我们如何确信它真的反映了物理世界的真相?今天,我们将从最基础也是最核心的步骤开始,学习如何像医生一样,为你的仿真数据做一次全面的‘体检’,并掌握表达核心信息的语言。”

正文核心:

1. 自定义表达式:让软件绘制你真正关心的物理量
  • 痛点:默认的物理量列表(如温度、应力)有时无法满足我们的分析需求。例如,在电磁-热耦合分析中,你可能想直接观察“焦耳热功率密度”的分布,以确定电路板上的热点。
  • 解决方案:COMSOL 允许你像写数学公式一样,手动输入任意表达式。
  • 实操案例
  1. 打开 COMSOL 案例库中的“母线板(Busbar)”模型。
  2. 添加一个“三维绘图组”,在“表面”图中,将表达式从默认的ec.normJ(电流密度模)改为自定义表达式ec.Jx*ec.Ex + ec.Jy*ec.Ey + ec.Jz*ec.Ez

3.对比:你会发现,这个手动输入的表达式与软件预定义的ec.Qrh(电阻损耗)结果完全一致。这验证了我们表达式的正确性。

  • 进阶技巧
    • 解构分析:我们可以进一步将焦耳热公式拆解,分别绘制ec.Jx*ec.Exec.Jy*ec.Ey。通过对比,你会发现热量主要集中在 X 方向的电流路径上,而 Y 和 Z 方向的贡献微乎其微。这种“解构”能力是自定义表达式最强大的地方,它能帮你洞察现象背后的主导因素。

ec.Jx*ec.Ex

ec.Jy*ec.Ey

2. 派生值:用几个关键数字,为模型做一次“体检”
  • 痛点:云图只能提供定性的视觉感受。如何快速、定量地验证模型的正确性?答案是通过计算一些关键的“派生值”。
  • 核心思想:利用全局计算(Global Evaluation)来检查最基本的物理守恒定律是否满足。
  • 实操案例
    1. 打开“散热器(Heat Sink)”模型。能量守恒检查:在“结果”下的“派生值”中,添加一个“全局计算”。在表达式中输入ht.ntefluxInt(总净能量速率)。计算结果应非常接近于你施加在模型上的总热功率(例如1W)。如果计算结果是0.99996 W,那就证明能量在系统内是守恒的,你的模型在能量层面是可信的。
    操作路径:模拟开发器窗口下—结果节点下—右击“派生值”选择并单击“全局计算”—在设置窗口找到“表达式”栏在输入框里输入ht.ntefluxInt
  • 质量守恒检查:再添加一个“全局计算”,输入表达式spf.out1.Mflow(质量流)。将这个值与入口速度、流体密度和入口面积手动计算出的理论质量流量进行对比。如果两者非常接近,则证明质量也是守恒的。操作路径:模拟开发器窗口下—结果节点下—右击“派生值”选择并单击“全局计算”—在设置窗口找到“表达式”栏在输入框里输入spf.out1.Mflow
  • 价值:派生值是验证模型正确性的“试金石”。只有当这些基本物理量都守恒时,我们才能对后续的复杂分析有信心。
3. 网格图:为你的“画布”做一次质量评估
  • 痛点:再漂亮的云图,如果建立在一张粗糙、扭曲的“画布”(网格)上,其结论也是不可靠的。因此,后处理的第一步,往往是从检查网格本身开始。
  • 解决方案:创建一个“网格图”,无需任何求解结果,即可直观地查看网格单元的质量和尺寸分布。
  • 实操案例
    1. 在“散热器”模型的“网格”节点上,直接点击“绘制”。解读:在生成的网格图中,颜色通常代表单元质量(Quality),范围为0到1。红色区域表示质量差(细长、扭曲)的单元。如果高应力、高梯度区域恰好位于这些红色单元上,那么你需要重新细化网格,否则结果可能失真。

价值:网格图是“事前诸葛”,它能帮助你在进行复杂分析之前,就发现并修正潜在的数值误差源头,避免“垃圾进,垃圾出”的悲剧。

结尾总结:

“至此,我们已经完成了对仿真数据的基础‘体检’。我们学会了如何‘说’出自己关心的物理量,用几个关键数字验证了模型的基本功,并检查了我们所依赖的‘画布’——网格——的质量。这是所有后续高级分析的地基。下一篇,我们将进入‘美工’阶段,学习如何让这些可信的数据,变成一张震撼、清晰、能打动人心的图像。敬请期待!”

接触类仿真 Contact Simulation

接触类仿真:从生活卡扣到工程应用的深度探索

引子:生活中的卡扣

在日常生活中,我们随处可见各种卡扣结构——从手机壳的卡扣到家具的拼接部件,从汽车内饰的固定装置到电子设备的连接接口。这些看似简单的结构背后,大家知道吗,其实他们背后蕴含着复杂的接触力学原理。作为一名仿真工程师,我常常思考:如何准确地模拟这些接触行为?不同的仿真软件又是如何处理接触问题的?

基于本人以往仿真经验,咱们讨论一下主流仿真软件(Abaqus和Comsol)两款软件中的接触算法

一、Abaqus中的接触算法

Abaqus作为通用有限元分析软件的标杆,提供了丰富的接触算法选项,其中接触类型有:

1.通用接触算法:适用于复杂接触场景,能自动处理接触对的识别

2.面-面接触:精度较高但计算成本也较高

3.点-面接触:计算效率较高但精度相对较低

其中接触方法包含硬接触和软接触,硬接触选择有:

1. 直接强制方式:使用Lagrange multiplier,严格满足无穿透条件

2. 接触罚函数法:通过”弹簧”概念处理穿透问题

3. 增广拉格朗日方式:增加叠代数,使用Penalty的方式,逼近真实行为。

软接触选择则包含有指数形式,线性形式和表格三种形式定义接触压力和过盈的关系。

Abaqus特别擅长处理大变形、大滑移的非线性接触问题,其接触算法在汽车碰撞、橡胶密封等工业应用中表现优异。

二、COMSOL中的接触算法

COMSOL作为多物理场耦合仿真专家,其接触算法具有自身特点:

1.基于物理场的接触模型:可以自然耦合到其他物理场(如热、电等)

2.惩罚接触法:默认采用的方法,计算效率高

3.增广拉格朗日法:提供更严格的接触约束

4.粘附接触:可模拟粘接-脱粘过程

5.参数化曲面接触:对特定几何形状有优化

COMSOL的突出优势在于接触与其他物理场的无缝耦合,比如可以直接研究接触压力对接触电阻的影响。

三、其中两种软件核心差异对比如下:

软件AbaqusCOMSOL
算法侧重点力学精度和稳定性多物理场耦合
大变形处理优异良好
多物理场耦合需要用户设置内置自然耦合
接触对自动识别强大相对简单
计算效率取决于算法选择通常较高
学习曲线较陡峭相对平缓

四、接触仿真的重要拓展

1.接触力学与接触电阻

在电气连接器、开关设备等应用中,接触压力会显著影响接触电阻:

1.接触面积理论:实际接触面积远小于表观面积

2.膜电阻与收缩电阻:表面氧化膜和电流线收缩效应

3.Holm接触理论:经典的电接触理论框架

通过耦合力学接触和电接触模型,可以预测连接器的接触电阻随插拔次数、振动环境的变化。

2.接触传热分析

接触热阻是许多散热设计中的关键参数:

1.表面粗糙度影响:微凸体决定了实际热传导路径

2.界面材料作用:导热硅脂、相变材料等的影响

3.接触压力-热阻关系:通常呈非线性关系

典型应用包括芯片封装散热、航天器热控制等。

五、科研与工程应用案例

案例1:连接器卡扣耐久性分析

使用Abaqus模拟插针卡扣在10万次开闭循环中的应力变化,预测塑性累积和失效位置,指导设计改进。

案例2:连接器接触电阻优化

在comsol中建立多物理场模型,分析不同镀层材料、接触正压力对接触电阻的影响,找到最佳参数组合。

案例3:动力电池模组膨胀力分析

研究锂离子电池充放电过程中的体积膨胀对电池模组接触压力的影响,优化约束设计。

结语

接触仿真作为连接虚拟与现实的关键技术,其精确度直接影响产品设计的可靠性。理解不同软件接触算法的特点,根据具体应用场景选择合适的工具和方法,是每位仿真工程师的必修课。随着计算技术的发展,接触仿真正从宏观向微观尺度延伸,从单物理场向多物理场耦合深化,为产品创新提供越来越强大的支持。

【Python】Streamlit 应用 “COMSOL 颜色表预览”

Streamlit APP “COMSOL 颜色表预览”

有些客户在拿到基于COMSOL的分析结果后,会提出“当前云图颜色配色不满意,需要拿到COMSOL文件自己进行进行调色”。但是大部分时候我是不提供COMSOL文件的。不提供文件的原因有如下几点:

  • 客户其实不会用这个软件,自己操作半天还是不会弄,最后还是返回来要我处理。
  • 对文件内些有独有技术的保护,不想外流。
  • 有些客户不懂有限元,看到文件内的操作不理解,会出现反复问询的情况,避免陷入反复自证吃了几碗粉的情况。
  • 有些客户拿到文件后自己的保密意识不到位,有外流风险。(以前出现客户拿到模型后毕业礼把模型文件卖给模型收集者,然后模型收集者又将文件挂到进行兜售交易。)

为了回绝客户的 “当前云图颜色配色不满意,需要拿到COMSOL文件自己进行进行调色”, 我开发一个在线的颜色表预览App,客户可以预览所有的COMSOL预置颜色表,选择好意向颜色表后进行反馈。

“COMSOL 颜色表预览” 应用入口

目前此APP仅对合作伙伴开放使用,如遇APP使用请及时联系我。欢迎各位给出宝贵意见~!

纯Python开发,使用了streamlit,pandas,pyvista等模块。

更新记录

2024年11月21日 V2.1

  • 增加了2D、3D默认视图的选择,方便对不同的维度的结果进行快速预览
  • 增加了颜色表的顺滑

2024年11月15日 V2.0上线

  • 增加COMSOL格式的colorbar渲染
  • 增加单个订单的多个场景预览选择
  • 在无查询订单状态下,增加了几个官方案例库里的后处理场景
  • 布局调整

COMSOL Matlab livelink多参数研究仿真 2

model = mphopen('1113.mph'); %打开文件

param_name = ["d1" "a3" "a2" "a1"]; % 参数名称
param_unit = ["um" "um" "um" "um"]; % 参数单位
param_Value = [0.1 0.3 0.5;...  %参数值
                50 75 100;...
                1 3 5;...
                0.1 1 2];

% For循环运行所有参数组
para_num=size(param_Value,2);
T_max=zeros(para_num^4,1);
Stress1=zeros(para_num^4,1);
Stress2=zeros(para_num^4,1);
para=ones(para_num^4,4);
i=1;
for para1=1:para_num % 循环研究参数名
    for para2=1:para_num
        for  para3=1:para_num
            for para4=1:para_num
                %记录输出
                tic;
                fprintf('Start simulation %d: %s=%.1f | %s=%.1f | %s=%.1f | %s=%.1f \n',i,...
                    param_name(1),param_Value(1,para1),...
                    param_name(2),param_Value(2,para2),...
                    param_name(3),param_Value(3,para3),...
                    param_name(4),param_Value(4,para4));
                %设定参数
                model.param.set(param_name(1),num2str(param_Value(1,para1)),param_unit(1));
                model.param.set(param_name(2),num2str(param_Value(2,para2)),param_unit(2));
                model.param.set(param_name(3),num2str(param_Value(3,para3)),param_unit(3));
                model.param.set(param_name(4),num2str(param_Value(4,para4)),param_unit(4));
                %运行计算
                model.sol('sol1').runAll;
                model.sol('sol2').runAll;
                model.sol('sol3').runAll;
                %记录数据
                para(i,:)=[param_Value(1,para1) param_Value(2,para2) param_Value(3,para3) param_Value(4,para4)];
                T_max(i) = mphglobal(model,'dom1','dataset','dset1','unit','degC'); %导出指定数据集的探针数据
                Stress1(i) = mphglobal(model,'dom2','dataset','dset3','unit','MPa');
                Stress2(i) = mphglobal(model,'dom6','dataset','dset4','unit','MPa');
                fprintf('Done simulation %d: T_max = %.2f | Stress1 = %.2f | Stress2 = %.2f\n',i,T_max(i),Stress1(i),Stress2(i));
                i=i+1;
                toc
            end
        end
    end
end
save('data.mat','para','T_max','Stress1','Stress2');

COMSOL Matlablive 多参数研究,批量化运行+导出+自动整理

COMSOL软件的研究求解模块其实提供了比较丰富的参数化运行结构,但是对不同参数的研究的结果导出方面,还是需要很多手动设置操作的才能导出批量化数据。

比如以下情况:研究对线是一个催化反应器,需要入口温度、入口流速、反应器尺寸、反应物初始孔隙率等四个不同的输入参数进行单变量的不同输入值的影响研究。那么如果直接使用COMSOL来做的话,就需要使用参数化扫描对单个参数进行扫描,然后导出动画进行分解,并保存四个不同的mph文件,用于分别研究四个参数。

针对这种情况,其实可以用Matlab livelink来进行批量的单个参数的修改,并运行求解导出相关结果(mph、图、报告、表数据)。

% 1、定义要运行的参数组
model = mphopen('param_run.mph');  % 打开模型
% 设定参数组
param_name = ["T_in" "U_in" "d_ball" "eps_init"]; % 参数名称
param_unit = ["degC" "m/s" "um" ""]; % 参数单位
param_Value = [600 650 700 750 800;...  %参数值
                0.1 0.2 0.3 0.4 0.5;...
                75 150 200 300 500;...
                0.03 0.05 0.1 0.15 0.2];

% 2、 For循环运行所有参数组
for name=1:length(param_name) % 循环研究参数名
   for value=1:size(param_Value,2) % 参数水平
       % 记录输出
       tic;
       fprintf('Start simulation: %s=%s [%s]\n',param_name(name),num2str(param_Value(name,value)),param_unit(name));
       % 设定参数
       model.param.set(param_name(name),num2str(param_Value(name,value)),param_unit(name));
       % 运行计算
       model.sol('sol1').runAll;
       % 导出数据
       model.result.report('rpt1').run; % 导出报告
       model.result.export('anim1').run; % 导出动画
       model.result.export('anim2').run;
       model.result.export('anim3').run;
       model.result.export('anim4').run;
       model.result.export('anim5').run;
       model.result.export('anim6').run;
       model.result.export('anim7').run;
       model.result.export('plot1').run; % 导出出图数据表
       model.result.export('plot2').run;
       model.result.export('plot3').run;
       % 整理导出数据
       movefile("export_dir", strcat(param_name(name),"_", num2str(param_Value(name,value))));
       copyfile("export_file",strcat(param_name(name),"_", num2str(param_Value(name,value))));
       % 完成
       toc
   end
   % 设定回原参数
   model.param.set(param_name(name),num2str(param_Value(name,1)),param_unit(name));
end

【Python】GifSpliter “GIF分割器”

“GifSpliter” GIF分割器

用PyQt5重写了一个“GIF分割器”的小应用:

功能:

  • 主要就是为了能快速分解GIF文件,将所有帧保存成PNG格式。
  • 支持文件拖拽进窗口,也是唯一的文件输入方式,就是简单暴力,直接方便,xD。
  • 可选择的输出位置:可以指定位置,也可以在原文件所在位置给你输出分割帧。
  • 可以拖拽文件夹输入,只要你文件夹里有gif,都给你处理了。
  • 支持混合输入模式。可以同时选中多个gif,或者gif+文件夹的多个文件输入。

现成的软件 GifSpliter_v3.2:

链接:https://pan.baidu.com/s/1JJLOF1vDYpd7pRzrtr93Xg?pwd=ucuo
提取码:ucuo
–来自百度网盘超级会员V6的分享

为什么做这个?

COMSOL仿真导出GIF图很方便,但是想要导出不同时刻参数下的多个图像很不方便。那……不如直接导出GIF后再分割?

但是网上没有这种轻量级的,简单暴力,合我需求的应用,那直接自己撸一个小工具!

现在,直接选中文件拖进去,等几秒就好了,多好?

Update History

v3.2 231031

  • 丰富了控制台打印信息,增加了彩色的控制台输出,并添加了处理进度条;

v3.1 230831

  • 保留了控制台窗口,并优化了控制台输出,现在可以通过控制台输出信息查看处理

    情况

v3.0 230631

  • PyQt5设计了全新的窗口;
  • 增加了文件夹模式;
  • 增加了多种输出模式的选择;

v2.0

  • 增加了多文件处理模式。

v1.0

  • Tinkter界面,单gif文件的处理方式,默认输出到exe所在文件夹。

SourceCode

import sys
import os
from PyQt5.QtWidgets import *
from PyQt5.QtGui import *
from ui_GifSpliter import *
from spliter import GifSplitter
import res_rc

class MyGifSpliter(QWidget):
    def __init__(self, parent=None):
        super().__init__(parent)  # 调用父类构造函数
        self.__ui = Ui_GifSpliter()
        self.__ui.setupUi(self)
        self.appDir = sys.path[0]
        print(self.appDir)
        self.gifDir = self.appDir
        self.gifDir = self.appDir
        self.specDir = self.appDir
        self.outDir = self.appDir
        self.__ui.textEdit_DirSpec.setText(self.appDir)
        # 调用Drops方法
        # self.drawn()
        self.__ui.radioButton_DirSpec.clicked.connect(self.do_setSpecDir)


    # 鼠标拖入事件
    def dragEnterEvent(self, event):
        if event.mimeData().hasUrls:
            event.accept()
        else:
            event.ignore()

    # 鼠标放开执行
    def dropEvent(self, event):
        """获取拖拽文件的路径"""
        files = []
        for url in event.mimeData().urls():
            path = url.toLocalFile()
            files.append(path)
        # 打印拖住进入窗口的文件路径
        self.gifDir = os.path.split(files[0])[0]
        # print(f'拖拽文件路径为:"{self.gifDir}"')

        # 判断输出模式
        if self.__ui.radioButton_DirGif.isEnabled():
            self.outDir = self.gifDir + '/分割后PNG'
            if not os.path.exists(self.outDir):
                os.mkdir(self.outDir)
        elif self.__ui.radioButton_DirApp.isEnabled():
            self.outDir = self.appDir
        elif self.__ui.radioButton_DirSpec.isEnabled():
            self.outDir = self.specDir
        # print(f'输出文件夹为:"{self.outDir}"')

        # 处理获取的文件
        for file in files:
            # 文件夹
            if os.path.isdir(file):
                # print(f'{file}是文件夹')
                self.gifDir = file
                self.outDir = self.gifDir + '/分割后PNG'
                if not os.path.exists(self.outDir):
                    os.mkdir(self.outDir)
                GifSplitter.split_gifs_in_folder(file, save_folder=self.outDir)
                print(f'\033[92m完成对文件夹:"{file}"的处理~~~\033[0m')
            # 文件
            else:
                # 判断扩展名是否为.gif
                if file.lower().endswith('.gif'):
                    # print(f'{file}是gif文件')
                    GifSplitter(file).export_frames(save_folder=self.outDir)
                else:
                    print(f'\033[91m! "{os.path.basename(file)}" 不是GIF文件~~~!\033[0m')
        if self.__ui.checkBox_openDirSpec.isChecked():
            os.startfile(self.outDir)

    #设定背景图片
    # def drawn(self):
    #     self.__ui.palette = QPalette()
    #     self.__ui.palette.setBrush(QPalette.Background, QBrush(QPixmap("BackGround1.png")))
    #     self.setPalette(self.__ui.palette)

    # 根据radioButton设定文件路径
    def do_setSpecDir(self):
        if self.__ui.radioButton_DirSpec.isChecked():
            self.specDir = QFileDialog.getExistingDirectory(None,"选取PNG输出文件夹","")
            self.__ui.textEdit_DirSpec.setText(self.specDir)
            print(f'设定指定输出路径为:{self.specDir}')


if __name__ == "__main__":
    app = QApplication(sys.argv)
    e = MyGifSpliter()
    e.show()
    sys.exit(app.exec_())
# -*- coding: utf-8 -*-

# Form implementation generated from reading ui file 'ui_GifSpliter.ui'
#
# Created by: PyQt5 UI code generator 5.15.9
#
# WARNING: Any manual changes made to this file will be lost when pyuic5 is
# run again.  Do not edit this file unless you know what you are doing.


from PyQt5 import QtCore, QtGui, QtWidgets


class Ui_GifSpliter(object):
    def setupUi(self, GifSpliter):
        GifSpliter.setObjectName("GifSpliter")
        GifSpliter.resize(640, 480)
        GifSpliter.setMinimumSize(QtCore.QSize(640, 480))
        GifSpliter.setMaximumSize(QtCore.QSize(640, 480))
        font = QtGui.QFont()
        font.setFamily("微软雅黑")
        font.setPointSize(11)
        font.setBold(True)
        font.setUnderline(False)
        font.setWeight(75)
        GifSpliter.setFont(font)
        GifSpliter.setAcceptDrops(True)
        icon = QtGui.QIcon()
        icon.addPixmap(QtGui.QPixmap(":/icon/C cyn.png"), QtGui.QIcon.Normal, QtGui.QIcon.Off)
        GifSpliter.setWindowIcon(icon)
        GifSpliter.setStyleSheet("")
        self.verticalLayout_2 = QtWidgets.QVBoxLayout(GifSpliter)
        self.verticalLayout_2.setSizeConstraint(QtWidgets.QLayout.SetMinimumSize)
        self.verticalLayout_2.setObjectName("verticalLayout_2")
        self.label = QtWidgets.QLabel(GifSpliter)
        self.label.setMinimumSize(QtCore.QSize(0, 0))
        font = QtGui.QFont()
        font.setFamily("微软雅黑")
        font.setBold(False)
        font.setWeight(50)
        self.label.setFont(font)
        self.label.setAcceptDrops(True)
        self.label.setFrameShape(QtWidgets.QFrame.NoFrame)
        self.label.setFrameShadow(QtWidgets.QFrame.Plain)
        self.label.setLineWidth(1)
        self.label.setTextFormat(QtCore.Qt.AutoText)
        self.label.setAlignment(QtCore.Qt.AlignCenter)
        self.label.setOpenExternalLinks(True)
        self.label.setObjectName("label")
        self.verticalLayout_2.addWidget(self.label)
        self.horizontalLayout = QtWidgets.QHBoxLayout()
        self.horizontalLayout.setSizeConstraint(QtWidgets.QLayout.SetMinimumSize)
        self.horizontalLayout.setObjectName("horizontalLayout")
        self.radioButton_DirGif = QtWidgets.QRadioButton(GifSpliter)
        self.radioButton_DirGif.setChecked(True)
        self.radioButton_DirGif.setObjectName("radioButton_DirGif")
        self.horizontalLayout.addWidget(self.radioButton_DirGif)
        self.radioButton_DirSpec = QtWidgets.QRadioButton(GifSpliter)
        self.radioButton_DirSpec.setChecked(False)
        self.radioButton_DirSpec.setObjectName("radioButton_DirSpec")
        self.horizontalLayout.addWidget(self.radioButton_DirSpec)
        self.radioButton_DirApp = QtWidgets.QRadioButton(GifSpliter)
        self.radioButton_DirApp.setChecked(False)
        self.radioButton_DirApp.setObjectName("radioButton_DirApp")
        self.horizontalLayout.addWidget(self.radioButton_DirApp)
        self.verticalLayout_2.addLayout(self.horizontalLayout)
        self.horizontalLayout_2 = QtWidgets.QHBoxLayout()
        self.horizontalLayout_2.setSizeConstraint(QtWidgets.QLayout.SetMinimumSize)
        self.horizontalLayout_2.setObjectName("horizontalLayout_2")
        self.label_2 = QtWidgets.QLabel(GifSpliter)
        self.label_2.setObjectName("label_2")
        self.horizontalLayout_2.addWidget(self.label_2)
        self.textEdit_DirSpec = QtWidgets.QTextEdit(GifSpliter)
        sizePolicy = QtWidgets.QSizePolicy(QtWidgets.QSizePolicy.Expanding, QtWidgets.QSizePolicy.Ignored)
        sizePolicy.setHorizontalStretch(0)
        sizePolicy.setVerticalStretch(0)
        sizePolicy.setHeightForWidth(self.textEdit_DirSpec.sizePolicy().hasHeightForWidth())
        self.textEdit_DirSpec.setSizePolicy(sizePolicy)
        font = QtGui.QFont()
        font.setPointSize(9)
        font.setBold(False)
        font.setWeight(50)
        self.textEdit_DirSpec.setFont(font)
        self.textEdit_DirSpec.setObjectName("textEdit_DirSpec")
        self.horizontalLayout_2.addWidget(self.textEdit_DirSpec)
        self.checkBox_openDirSpec = QtWidgets.QCheckBox(GifSpliter)
        self.checkBox_openDirSpec.setObjectName("checkBox_openDirSpec")
        self.horizontalLayout_2.addWidget(self.checkBox_openDirSpec)
        self.verticalLayout_2.addLayout(self.horizontalLayout_2)
        self.verticalLayout_2.setStretch(0, 1)

        self.retranslateUi(GifSpliter)
        QtCore.QMetaObject.connectSlotsByName(GifSpliter)

    def retranslateUi(self, GifSpliter):
        _translate = QtCore.QCoreApplication.translate
        GifSpliter.setWindowTitle(_translate("GifSpliter", "GIF文件分割器 v3.0"))
        self.label.setText(_translate("GifSpliter", "<html><head/><body><p align=\"center\">GIF文件分割器:   可以快速将GIF文件分割成逐帧PNG文件</p><p align=\"center\">-------------------------------------------------------------</p><p align=\"center\">使用方法:</p><p align=\"center\">将GIF文件或含有多个GIF的文件夹拖拽进此窗口即可</p><p align=\"center\">程序会自动分解所有拖拽进来的GIF并分解成单张PNG文件</p><p align=\"center\">可以在下面单选按钮中选择PNG文件输出位置</p><p align=\"center\">如果需要选择生成到特定文件夹,需要指定输出文件夹</p><p align=\"center\"><br/>-------------------------------------------------------------</p><p align=\"center\"><a href=\"https://guagefangzhen.cn/\"><span style=\" text-decoration: underline; color:#0000ff;\">更多分享:https://guagefangzhen.cn</span></a><br/></p></body></html>"))
        self.radioButton_DirGif.setText(_translate("GifSpliter", "生成到GIF所在目录"))
        self.radioButton_DirSpec.setText(_translate("GifSpliter", "生成到指定目录"))
        self.radioButton_DirApp.setText(_translate("GifSpliter", "生成到程序所在目录"))
        self.label_2.setText(_translate("GifSpliter", "指定输出目录:"))
        self.textEdit_DirSpec.setHtml(_translate("GifSpliter", "<!DOCTYPE HTML PUBLIC \"-//W3C//DTD HTML 4.0//EN\" \"http://www.w3.org/TR/REC-html40/strict.dtd\">\n"
"<html><head><meta name=\"qrichtext\" content=\"1\" /><style type=\"text/css\">\n"
"p, li { white-space: pre-wrap; }\n"
"</style></head><body style=\" font-family:\'SimSun\'; font-size:9pt; font-weight:400; font-style:normal;\">\n"
"<p style=\" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px;\">指定分割帧文件输出路径</p></body></html>"))
        self.checkBox_openDirSpec.setText(_translate("GifSpliter", "分割后打开目录"))
import res_rc
# 用于定义gif分割程序
import os
import imageio


class GifSplitter:
    def __init__(self, gif_file):
        self.gif_file = gif_file
        self.gif = imageio.get_reader(gif_file)

    def get_gif_info(self):
        """获取gif信息"""
        info = {}
        info['file_name'] = os.path.basename(self.gif_file)
        info['folder_path'] = os.path.dirname(self.gif_file)
        info['total_frames'] = len(self.gif)
        return info

    def export_frames(self, save_folder=None):
        """导出gif所有帧为png图片"""
        if save_folder is None:
            save_folder = os.path.join(self.get_gif_info()['folder_path'],
                                       'splitted')
            if not os.path.exists(save_folder):
                os.mkdir(save_folder)

        for i, frame in enumerate(self.gif):
            imageio.imsave(os.path.join(save_folder, '{}-{:03d}.png'.format(
                self.get_gif_info()['file_name'], i)), frame)
        print(f'完成分割: "{(os.path.basename(self.gif_file))}" | 输出文件夹-->"{save_folder}"')


    @classmethod
    def split_gifs_in_folder(cls, folder_path, save_folder=None):
        """分割文件夹内所有的gif文件"""
        print(f'拖拽文件夹为:"{folder_path}"')
        print(f'输出文件夹为:"{save_folder}"')
        if save_folder is None:
            save_folder = os.path.join(folder_path, 'splitted')
            if not os.path.exists(save_folder):
                os.mkdir(save_folder)

        gif_files = [file for file in os.listdir(folder_path)
                     if file.endswith('.gif')]

        i = 1
        for gif_file in gif_files:
            gif_splitter = cls(os.path.join(folder_path, gif_file))
            gif_splitter.export_frames(save_folder)
            print(f'文件夹处理进度: {i}/{len(gif_files)}')
            i = i + 1

if __name__ == '__main__':
    # 分割单个gif文件
    gif_splitter = GifSplitter('test.gif')
    gif_splitter.export_frames(save_folder="gif")

    # 分割文件夹内所有gif文件
    GifSplitter.split_gifs_in_folder('gif_folder')

【Python】Streamlit应用 “这个班上的值不值”

“这个班上的值不值” Web版

做了一个Streamlit的练手应用,用于测算“工作的性价比”。

非严谨计算,请勿认真对待,xD。

挂在服务器上了,比较轻量,应该不会下线,可以分享给朋友玩一下http://175.24.226.62:8501/

缘由

去年在某位群友的分享下,碰到一个很有意思的一个excel小应用。

一个挺有意思的小应用,一点职场人的茶余饭后的小话题。玩了两下后思考了一下,这个做成excel属实是不方便好友分享,如果做成一个web版的或者小程序版本,不是更好?这个想法就这样一直留在todo list里了。

最近工作太忙了,但是为了不被眼下的工作完全淹没自己,还是逼迫自己慢慢的完成一些todo list上的事。因为,我坚持认为,人除了要赚眼前的填饱肚子的月供之外,也要持续的坚持学习新东西,做一些对长远规划有用的东西。

五月份某天在学习调研dashboard和streamlit的时候,想起这个“这个班上的值不值”的应用,就把它当做一个学习练手的小task做了,并且部署到了我的轻量级服务器上。


做法分享:

简单来讲,就是分为两步:

1、基于python streamlit写脚本,开发一个应用。

2、将调试好的streamlit应用部署到服务器上。

主要要点:

1、我的开发环境是Anaconda+Pycharm,这也是主流的比较推荐的python科学研究的开发环境;

2、Streamlit相关的资料还是比较少,主要建议上官网看他的一些API说明:

Streamlit官方Documentation(推荐):https://docs.streamlit.io/

也有个国人做的学习手册(施工中):http://cw.hubwiz.com/card/c/streamlit-manual/

然后有了一些基础概念后,就可以面对AI编程了,哈哈。

3、我是梳理了一下应用,基本很快就确定了一个框架:siderbar作为参数输入,然后右边做输出和说明展示的。进而开始coding。

import streamlit as st

with st.sidebar:
    st.title('输入参数')
    salary = st.number_input('平均日薪酬(RMB)', value=100, min_value=50, step=25)
    work_t = st.number_input('工作时长(单位:小时,下班时间-上班时间)', value=8.0, min_value=1.0, step=0.5)
    trans_t = st.number_input('通勤时长(单位:小时,上下班花在路上的通勤时间)',value=1.0,step=0.25)
    slack_t = st.number_input('摸鱼时间(单位:小时,吃饭+不干活+午休)',value=1.0,min_value=0.0,step=0.25)

    st.title('相关系数')
    edu = st.number_input('学历系数',value=1.0,step=0.2,min_value=0.8,max_value=2.0)
    with st.expander('系数说明'):
        st.text('专科及以下 0.8\n普通本科  1.0\n高级本科  1.2\n普通硕士  1.4\n高级硕士  1.6\n普通博士  1.8\n高级博士  2.0')
    env_working = st.number_input('工作环境系数',value=1.0,step=0.2,min_value=0.8,max_value=1.1)
    with st.expander('系数说明'):
        st.text('偏僻地区  0.8\n工厂户外  0.9\n普通     1.0\n体制内   1.1')
    env_female = st.number_input('异性环境指数',value=1.0,min_value=0.9,max_value=1.1,step=0.1)
    with st.expander('异性指数说明'):
        st.text('没有      0.9\n不多不少  1.0\n很多     1.1')
    env_coleg = st.number_input('同事环境系数', value=1.0, min_value=0.95, max_value=1.05, step=0.05)
    with st.expander('同事环境系数说明'):
        st.text('SB很多   0.95\n普通很多   1.0\n优秀很多  1.05')
    env_career = st.number_input('职业资格系数', value=1.0, min_value=1.0, max_value=1.15, step=0.05)
    with st.expander('职业资格系数说明'):
        st.text('无要求、二级    1.0\n建造造价监理    1.05\n建筑岩土结构    1.1\n主任医师、教授  1.15')
    early_working = st.number_input('是否在8:30前上班', value=1.0, min_value=0.95, max_value=1.0, step=0.05)
    st.text('是否8:30前上班? 是:0.95, 否:1.0')
    final_ratio = st.number_input('综合环境系数(不要改)', value=1.0, max_value=1.1, min_value=0.9, step=0.05)

point = salary*final_ratio/(35*(work_t+trans_t-0.5*slack_t)*edu*env_career)*early_working*env_coleg*env_female

st.title(':male-technologist:这个班上的值不值?')


st.caption(f"您的工作性价比指数为: _%.2f_" % point)
if point<=0.8:
    text_out = '参考评价:您的工作:green[很惨]~~~  :weary:'
elif point<=1.5:
    text_out = '参考评价:您的工作:blue[一般]. :expressionless:'
elif point<=2.0:
    text_out = "参考评价:您的工作:red[很爽]~~!  :stuck_out_tongue_winking_eye:"
else:
    text_out = "参考评价:您的工作:red[爽到爆炸]~~~!!!  :sunglasses:"
st.caption(text_out)
st.caption('_这是一个娱乐性APP,通过对在左侧侧边栏<-输入的几个参数进行计算,评估自己这个班上的值不值~_')
st.caption('_仅为博君一笑,请勿认真对待~!_')

with st.expander("计算公式说明:"):
    st.image('img.png')
st.write('[Powered by 瓜哥 using "streamlit"](https://guagefangzhen.cn)')

4、Streamlit的发布到服务器,我百度了看到的基本都是用Docker环境安装streamlit环境。我的想法更简单: 服务器上装了miniconda –> 在miniconda里装Streamlit环境 –> 在Streamlit虚拟环境中启动上传的 JobPerformanceCalculator.py

5、Streamlit中的启动,需要添加nohup命令。这样你就可以在运行后关闭这个命令窗口,而应用一直会在服务器后台运行了。

nohup streamlit run JobPerformanceCalculator.py