304不锈钢贮箱箱底封闭液胀 航天级精简复刻方案（HyperMesh+LS-DYNA CEL）
一、软件链路（一套直接复刻，无多余工具）
1. 几何：CATIA/SolidWorks 出STEP 2. 前处理网格：HyperMesh（LS-DYNA Profile） 3. 求解内核：LS-DYNA R14（CEL流固耦合唯一选择） 4. 监控后处理：LS-PrePost 5. 单位统一：mm、N、ton、MPa、s
二、工况速记（对应你的结构）
上下曲底+环形片一体封闭坯料（304）；模具235B仅轴向端面摩擦限位、无径向包覆；内腔预灌满清水无气体；缓慢密闭增压液胀，同步整形上下曲率；全程封闭无溢流口。
三、上午/下午分步复刻操作（8步，3小时可搭完模型）
步骤1：几何清理（20min）
1. 导入STEP，分3个组件Collector
- BLANK：304箱底整体坯料
- DIE_UP/DIE_DOWN：上下235B轴向限位凹模
- FLUID：内腔欧拉流体包裹体 2. 清理微小特征：删除≤1mm工艺倒角、小孔；保留曲底圆角、环形过渡圆弧、坯料轴向摩擦端面 3. 流体域生成：提取坯料内壁封闭面，整体向外偏移2mm做封闭六面体包裹腔，充液小孔封闭，无外流通道 4. 干涉检查：凹模内径＞坯料最大胀形外径，径向全程不接触
步骤2：网格划分（40min，航天网格标准）
1）BLANK 304坯料（SHELL163，5点厚度积分）
- 加密区（2.5–3mm）：曲底过渡圆角、环形连接段、环-底衔接圆弧（起皱/减薄危险区）
- 平缓区（4–6mm）：球面中心
- 质控：全四边形、长宽比≤3、雅克比≥0.7，无自由边
- 输入实测板厚t
2）DIE上下凹模（SHELL163 刚体壳）
网格6–10mm，仅保证环形接触面完整，无需加密
3）FLUID欧拉域（SOLID164 六面体，禁用四面体）
流体网格=坯料网格×1.2；环形、圆角内侧局部加密；全域连续无空洞，初始100%液相填充。
步骤3：材料+状态方程赋值（15min，参数直接抄）
1）BLANK 304不锈钢
plaintext
*MAT_PIECEWISE_LINEAR_PLASTICITY 1,7.93e-9,193000,0.3,0.001 真实应力-塑性应变曲线（0~0.75应变） *MAT_ANISOTROPIC_HILL48 轧制板材各向异性，抑制环向异常流动 *FLD 导入304成形极限，判破裂
2）DIE 235B刚体
plaintext
*MAT_RIGID 2,7.85e-9,206000,0.3
刚体约束：仅开放Z轴向平移，X/Y平动+三轴旋转全锁死
3）FLUID 清水可压缩CEL（封闭腔核心）
plaintext
*MAT_NULL 3,1.0e-9,0,0 *EOS_GRUNEISEN 3,1480,0.19,0,0,0
初始欧拉域气相分数0，无空气。
步骤4：接触搭建（15min，匹配轴向摩擦限位）
1. 凹模-坯料端面摩擦接触 *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE 主面：凹模环形端面；从面：坯料上下端面 μs=0.15，μd=0.09；不加固支约束，允许坯料轴向滑动补料；刚度缩放0.8防穿透 2. 坯料自接触： *CONTACT_AUTOMATIC_SINGLE_SURFACE 防止环形片褶皱嵌套 3. CEL自动流固耦合，无需手动建流体-固体接触
步骤5：工艺三阶段加载（20min，缓慢液胀整形）
总计算时长0.6s，分3连续步，封闭腔无流体进出
步1 平衡步（0~0.01s）
凹模Z向全固定；流体全域预压2MPa，消除网格间隙震荡
步2 缓慢增压主胀形（0.01~0.4s，核心）
采用 *LOAD_CEL_VOLUME 线性增压，升压速率0.5MPa/s；无主动推模，仅靠液压胀形，坯料端面沿凹模滑移补料；同步拉伸整形上下曲底
步3 保压整形（0.4~0.6s）
压力维持峰值恒定，静压力校平曲面回弹、修正曲率凹陷
步骤6：求解控制关键字（复制粘贴直接用，解决CEL报错）
plaintext
*CONTROL_CEL 2,0,0,0,0,0 !二阶界面、关闭欧拉流出（封闭腔） *CONTROL_HOURGLASS 1,0.05,0.08,0,0,0 !全局0.05，环形加密区局部0.08防假皱 *CONTROL_TIMESTEP 0,0,0,0,0,0 !关闭质量缩放，航天高精度强制要求 *CONTROL_ENERGY 1,1,1,1,1 !监控沙漏能/内能，合格标准＜5% *DAMPING_GLOBAL 0.02,0,0,0,0,0 !弱阻尼消除环形小幅振动 *CONTROL_SHELL 20,0,0,0,0,0 !全积分壳抑制剪切自锁
步骤7：输出设置（5min）
plaintext
*DATABASE_BINARY_D3PLOT 0.005,0,0,0,0,0 *DATABASE_GLSTAT 0.001 *DATABASE_RCFORC 0.001
输出：壁厚、等效应变、环/径向分应变、内腔全场压力、凹模支反力、端面滑移距离、曲面径向位移（曲率偏差）
步骤8：导出K文件提交计算（5min）
1. Export Solver Deck输出.k文件 2. 单机32核并行提交；实时监控LS-PrePost：沙漏能、压力曲线、无负体积、无网格畸变
四、仿真运行物理过程（快速理解，方便调参）
1. 平衡：低压消除初始接触间隙 2. 缓慢密闭增压：液体可压缩使内腔压力匀速上升；坯料径向自由胀形，环形片环向受压、曲底双向拉伸；坯料两端沿235B凹模滑动补料，摩擦仅限位不锁死 3. 保压：恒定压力释放局部回弹，校正上下椭球曲率尺寸
五、后处理快速评判标准（仿真跑完立刻出结论）
1. 减薄：曲底圆角最大减薄率≤15%合格 2. 起皱：环形段环向压应变＞-0.06存在失皱风险 3. 曲率整形：曲面母线径向位移偏差≤0.2mm 4. 数值可信：沙漏能/总内能＜5%，无负体积报错 5. 压力对标：压力-时间曲线线性平稳，无剧烈冲击震荡
六、快速排错清单（建模极易踩坑，下午直接对照）
1. 计算负体积：流体网格加密、降低升压速率、打开二阶CEL界面 2. 环形片大量假皱：提高局部沙漏阻尼、加大全局弱阻尼、降低动摩擦系数提升补料 3. 压力曲线剧烈跳变：欧拉域不允许存在气相、细化圆角流体网格 4. 坯料无法滑移、过度减薄：删除端面固支约束，检查接触摩擦设置 5. 计算速度过慢：适度粗化欧拉流体网格，仅危险区加密
七、两套使用模式（按需选择）
1. 工程快速迭代（3–6h跑完）：适度粗化流体网格，省略平衡步，仅输出壁厚、峰值压力、起皱风险，批量调试压力、摩擦参数 2. 航天鉴定高精度版（论文/工艺评审）：完整三步+二阶CEL+Hill48+FLD失效、无质量缩放，输出全域压力梯度、屈曲演化、曲率回弹全套数据，可直接用于仿真鉴定报告