GeoBC / research log
← 任务目录

R02 Connecting rod 仿真代码

已完成
01

验收 / 总结 / 介绍

完整FEn完成20个加载卸载状态,与同网格细Cal满载结果的位移相对L2差为1.06e-7。

02

问题 / review记录

细Cal卸载未完成;当前本构交叉验证只覆盖未屈服分支,不代表三方全路径验收通过。

评论与审阅

昵称:
自己在本浏览器提交的评论可再次编辑。

展开详情:参数、方法、交付链接与完整review记录
R02 / SIMULATION

Connecting rod 仿真代码

固定物理定义、明确本构范围,保留可复核的原始场。粗Cal、细Cal与GPU FEn三方消融已经完成。

阶段第1周
开始日期2026-10-22
截止日期2026-10-23
状态已完成
前置任务R01

单位与材料

长度mm、力N、应力MPa;E=115100MPa,ν=0.28,密度4430kg/m³。CalculiX使用有限应变J2;材料摘要点不外推,不将其当作损伤、断裂或疲劳认证。

五个目标工况

工况曲柄角净压力 [MPa]转速 [rpm]宽度偏载κ
压缩30°430000
拉伸360°060000
组合弯曲60°1.545000.3
超压15°645000
超速360°090000

边界与求解

  • 大孔参考径向MPC,加切向位移、轴向位移的面积加权积分均值约束。
  • 载荷按原生网格重新组装,包含承压面分布力与完整XY惯性体载荷。
  • Cal:C3D10、NLGEOM、严格Newton收敛控制,逐加载/卸载状态输出。
  • FEn:同一曲边P2网格、P2位移和4点积分,仅实现与Cal源码匹配的J2初始未屈服分支;任一积分点达到781.6MPa则停止。
R02验收完成 · GPU迁移可行

20个加载卸载状态已完成,耗时22分51.911秒;满载复核通过,细Cal卸载仍不完整。GPU承担PETSc线性代数,CPU承担单元装配;不声称整个有限元过程都在GPU上。

脚本与阶段结果

源码为review快照;大体积原始求解输入未随展示站发布,不宣称本站本身是独立完整复现包。

问题 / review记录

问题已确认事实后续
细Cal完整协议未完成达到满载约77分33秒;运行103分26秒后触发资源保护,仅完成两个卸载增量保留满载快照,完整卸载结果不得冒充完成
资源保护原因日志未记录当时具体越线资源,无法断言是RAM还是磁盘新FEn守护记录实际资源值
跨求解器验证小型曲边基准及完整连杆满载逐点复核通过;Cal全路径对照不完整保留满载复核;完整Cal卸载仍待完成

结果与交付链接

每日更新与review

追加更新 →

正在读取记录…

10月24日更新 · 粗Cal / 细Cal / GPU FEn消融完成

四栏交互三维对比 ↗
指标细Cal同网格FEn
满载最大位移4.426359184 mm4.426358799 mm
原始GP应力峰值304.040976 MPa304.040989 MPa
达到满载累计时间4652.822秒707.173秒
整个加载卸载路径未完成,6206.06秒时资源保护停止20状态完成,1371.911秒
逐节点位移相对L2差1.06144e-7;最大绝对差5.11268e-7 mm
同物理GP张量体积加权相对差1.58026e-7;GP坐标匹配误差2.28737e-13 mm

FEn采用CPU装配与CUDA线性代数,全部20状态未屈服、PEEQ为0,卸载后最大位移1.11076e-8 mm。满载约束平衡及三个既定应力ROI检查通过。消融结论:粗Cal/细Cal量化网格敏感性;细Cal/同网格FEn隔离求解器与CUDA后端差异;跨单元平均只改变显示。GPU改造在0006未屈服分支上验证可行。后续正式求解统一使用通过质量门槛的细网格,粗网格只保留为本次消融对照。细Cal卸载不完整,因此不外推为通用塑性或三方全路径等价;这也不等于细网格的真实离散误差为零。