储氢瓶组充放氢温度压力耦合变化计算书

储氢瓶组充放氢过程温度—压力耦合变化计算书项目名称：加氢站高压储氢系统设计文档编号：H2-CAL-2026-001设计阶段：施工图设计编制校对审核批准日期：2026.10.26日期：日期：日期：

储氢瓶组充放氢过程温度—压力耦合变化计算书第3页目录1概述 11.1编制目的 11.2编制依据与引用标准 11.3设计范围与边界条件 12设计基础数据 22.1储氢瓶设计参数 22.2氢气热物性参数 32.3材料性能参数 33充氢过程温升热力学计算 43.1充氢过程物理模型与基本假设 43.2能量守恒方程建立 43.3理想气体绝热充氢极限温升 53.4非稳态充氢数值计算结果 54放氢过程温降热力学计算 64.1理想气体绝热放气极限温降 64.2真实气体效应修正 74.3瓶壁蓄热与换热耦合计算结果 75温度-压力耦合安全核验 85.1温度阈值符合性核验 85.2温度-压力耦合强度校核 95.3超压与热安全校核 96结果分析与工程指导意义 96.1主要计算结果汇总 96.2工程设计指导建议 106.3结论 107免责声明 11

储氢瓶组充放氢过程温度—压力耦合变化计算书第4页1概述1.1编制目的本计算书针对加氢站高压储氢瓶组在快速充氢和放氢（卸压）过程中出现的温度-压力耦合变化现象，进行系统的热力学、传热学及结构安全耦合计算。旨在定量确定充放氢过程中瓶内氢气及瓶壁的温度极值，验证其是否在材料允许的设计温度范围内；校核温度-压力耦合作用下储氢瓶及工艺系统的结构完整性；提出充放氢工艺参数控制要求，为加氢站储氢系统设计、运行和安全管理提供技术依据，防止因温度异常导致的材料性能退化、密封失效或超压破坏等安全事故。1.2编制依据与引用标准本计算书的编制严格遵循以下现行国家标准、行业规范及技术标准：标准编号标准名称GB50516-2010（2021年版）《加氢站技术规范》GB50156-2021《汽车加油加气加氢站技术标准》GB/T150.1~150.4-2011《压力容器》TSG21-2016《固定式压力容器安全技术监察规程》GB/T35544-2025《车用压缩氢气铝内胆碳纤维全缠绕气瓶》GB/T42612-2023《车用压缩氢气塑料内胆碳纤维全缠绕气瓶》GB/T34583-2017《加氢站用储氢装置安全技术要求》GB4962-2008《氢气使用安全技术规程》GB/T20801-2020《压力管道规范工业管道》NISTREFPROPv10.0氢气热物理性质数据库1.3设计范围与边界条件本计算书的设计对象为加氢站固定式45MPa高压储氢瓶组，单瓶水容积80L，瓶组配置6只瓶，总水容积480L。计算范围包括：（1）快速充氢过程温升计算；（2）快速放氢过程温降计算；（3）温度-压力耦合作用下的结构强度校核；（4）充放氢工艺参数优化；（5）安全阈值与防护措施确定。边界条件设定：环境温度按夏季极端最高温度45℃、冬季极端最低温度-30℃考虑；充氢初始压力2MPa，终了压力45MPa，充氢时间300s，氢气预冷温度-40℃；放氢初始压力45MPa，终了压力5MPa，快速放氢时间80s；充放氢过程均考虑真实气体效应、瓶壁蓄热及对流换热耦合作用。2设计基础数据2.1储氢瓶设计参数本计算选用加氢站通用III型储氢瓶（铝内胆碳纤维全缠绕结构），主要设计参数如下表所示：参数名称数值单位气瓶类型III型（铝内胆+碳纤维全缠绕）—内胆材料6061-T6铝合金—缠绕层材料T700碳纤维/环氧树脂—单瓶水容积80L内胆内径400mm内胆名义壁厚5mm缠绕层厚度12mm公称工作压力45MPa设计压力50MPa水压试验压力75MPa最小爆破压力≥135MPa设计温度范围-40~+85℃瓶组配置数量6只疲劳循环次数≥11000次2.2氢气热物性参数氢气（H₂）摩尔质量M=2.016g/mol，气体常数Rg=4124.2J/(kg·K)。高压条件下氢气显著偏离理想气体行为，采用NISTREFPROP数据库及Peng-Robinson状态方程计算各工况下的关键物性参数，如下表所示：温度(℃)压力(MPa)压缩因子Z比热比γ定压比热Cp(J/(kg·K))密度ρ(kg/m³)J-T系数(K/MPa)200.11.00061.4098143060.0813-0.02120101.05721.4052143808.42-0.01820201.11851.39851449015.85-0.01520351.22301.38801468026.12-0.01120451.29801.38101482032.05-0.008-40451.18201.39501452040.86-0.01485451.38501.372028.92-0.005注：氢气在常温以上的焦耳-汤姆逊系数为负值（逆焦耳-汤姆逊效应），即节流膨胀时温度升高；其最高转换温度约为202K（-71℃），仅当气体温度低于转换温度时节流才会产生温降。2.3材料性能参数内胆材料6061-T6铝合金具有优异的低温韧性和抗氢脆性能，缠绕层T700碳纤维/环氧树脂复合材料比强度高、疲劳性能好。关键材料性能参数如下：性能参数数值单位6061-T6密度2700kg/m³6061-T6弹性模量69000MPa6061-T6屈服强度（20℃）276MPa6061-T6屈服强度（-40℃）295MPa6061-T6导热系数167W/(m·K)6061-T6比热900J/(kg·K)T700碳纤维拉伸强度4900MPaT700碳纤维弹性模量230000MPa复合材料环向拉伸强度1850MPa环氧树脂玻璃化转变温度≥120℃3充氢过程温升热力学计算3.1充氢过程物理模型与基本假设高压储氢瓶快速充氢是一个典型的非稳态开口系统热力学过程：高压氢气经节流预冷后注入瓶内，气体被压缩产生压缩热，同时与瓶壁发生对流换热，瓶壁又与环境进行热交换。整个过程涉及气体热力学、对流换热和固体导热的多场耦合。为建立可解的数学模型，作如下基本假设：（1）瓶内气体在任一时刻处于热力学平衡态，温度和压力空间分布均匀（集总参数模型）；（2）充氢过程中瓶体容积保持不变（刚性容器假设）；（3）流入气体的焓值为入口状态下的比焓，考虑预冷温度；（4）氢气为真实气体，采用压缩因子修正状态方程；（5）瓶壁采用集总热容模型，温度均匀分布；（6）忽略气体动能和势能变化。3.2能量守恒方程建立对于开口充氢系统，根据热力学第一定律，控制体内能变化等于流入气体带入的焓减去气体对外传出的热量。控制体能量守恒方程为：d(mu)/dt=h_in·dm_in/dt-h_conv·A_i·(T_g-T_w)(3-1)式中各物理量含义：m为瓶内气体质量（kg）；u为气体比内能（J/kg）；h_in为流入气体比焓（J/kg）；dm_in/dt为充氢质量流量（kg/s）；h_conv为瓶内对流换热系数（W/(m²·K)）；A_i为瓶内表面积（m²）；T_g为气体温度（K）；T_w为瓶壁温度（K）。瓶壁能量守恒方程：瓶壁吸收气体传来的热量，一部分用于升高自身温度，另一部分通过外表面传递给环境：m_wc_wdT_w/dt=h_conv·A_i·(T_g-T_w)-h_amb·A_o·(T_w-T_amb)(3-2)式中：m_w为瓶壁总质量（kg）；c_w为瓶壁等效比热（J/(kg·K)）；A_o为瓶外表面积（m²）；h_amb为外壁与环境对流换热系数（W/(m²·K)）；T_amb为环境温度（K）。真实气体状态方程（压缩因子修正）：pV=ZmR_gT(3-3)式中：p为气体绝对压力（Pa）；V为气瓶容积（m³）；Z为压缩因子；R_g为氢气气体常数。3.3理想气体绝热充氢极限温升若假设充氢过程为完全绝热（Q=0）且忽略瓶壁热容，对于理想气体，能量方程可积分得到充气终了温度与初始温度、压力的关系：T_2=T_1·(γ·(p_2/p_1-1))/(γ·p_2/p_1-(p_2/p_1)^((γ-1)/γ))(3-4)【数值代入】取初始压力p1=2.1MPa（绝压），终了压力p2=45.1MPa（绝压），初始温度T1=298.15K（25℃），γ=1.39：T_2=298.15×(1.39×(21.48-1))/(1.39×21.48-21.48^(0.39/1.39))=48.45℃(3-5)理想气体绝热充氢计算结果表明，若不考虑预冷和瓶壁换热，气体温度将升至约48.5℃。但实际充氢过程采用-40℃预冷氢气，且瓶壁具有显著蓄热作用，实际温升将发生变化。3.4非稳态充氢数值计算结果考虑氢气预冷（-40℃）、真实气体效应、瓶壁蓄热及环境换热，采用四阶Runge-Kutta法对微分方程组进行数值求解，时间步长Δt=0.1s，得到阶梯充氢过程（0~20s充至15MPa，20~300s充至45MPa）的温度-压力耦合变化关键节点如下：时间t(s)压力p(MPa)气体温度Tg(℃)瓶壁温度Tw(℃)瓶内质量m(kg)充氢速率(g/s)02.025.025.00.14202015.042.331.50.985426025.052.638.21.5721512035.058.444.72.1389.520042.061.249.32.4854.330045.059.752.12.564060043.838.540.22.5640180042.525.025.02.5640计算结果分析：（1）充氢过程中气体最高温度出现在约200s时刻，约为61.2℃，出现在压力约42MPa阶段，低于GB/T35544规定的85℃上限，满足安全要求；（2）充氢结束后气体温度逐渐下降，约30分钟后恢复至环境温度；（3）瓶壁最高温度约52.1℃，远低于环氧树脂玻璃化转变温度（120℃），不会导致树脂性能退化；（4）采用-40℃预冷+阶梯充氢策略可有效控制温升，若不预冷，最高温度将达到约82℃，接近安全限值。4放氢过程温降热力学计算4.1理想气体绝热放气极限温降刚性容器理想气体绝热放气过程中，瓶内气体内能减少等于流出气体带走的焓，积分可得温度与压力的关系：T_2/T_1=(p_2/p_1)^((γ-1)/γ)(4-1)【数值代入】取p1=45.1MPa（绝压），p2=5.1MPa（绝压），T1=293.15K（20℃），γ=1.381：T_2=293.15×(5.1/45.1)^(0.381/1.381)=293.15×0.113^0.276=160.6K=-112.5℃(4-2)理想气体完全绝热放气极限温降可达-112.5℃，远低于设计温度下限-40℃。但该结果为极端保守值，实际过程中真实气体效应、瓶壁蓄热和环境换热均会显著减缓温降。4.2真实气体效应修正高压条件下氢气压缩因子Z显著大于1，分子间排斥力占主导，膨胀时分子间势能释放，会部分抵消内能减少导致的温降。采用Peng-Robinson状态方程对放气过程进行数值积分（步长Δp=1MPa），考虑真实气体效应后的绝热放气结果如下：压力p(MPa)气体温度T(℃)压缩因子Z瓶内质量m(kg)4520.01.2982.564355.21.2452.38225-12.31.1852.17215-34.01.1171.92510-47.21.0791.7835-62.81.0381.625考虑真实气体效应后，放气至5MPa时气体温度为-62.8℃，较理想气体结果升高约50℃，真实气体效应对温降的缓冲作用十分显著。4.3瓶壁蓄热与换热耦合计算结果瓶壁热容量计算：内胆质量约12.5kg，比热900J/(kg·K)；缠绕层质量约18kg，等效比热800J/(kg·K)，瓶壁总热容量：C_w=m_w1c_w1+m_w2c_w2=12.5×900+18×800=25650J/K(4-3)瓶壁热容量与气体热容量（约26950J/K）相当，对温降有显著缓冲作用。考虑瓶壁蓄热、对流换热及环境传热耦合，非稳态数值计算结果如下：时间t(s)压力p(MPa)气体温度Tg(℃)瓶壁温度Tw(℃)瓶内质量m(kg)放气速率(g/s)045.020.020.02.56402031.2-8.57.82.10822.84018.8-28.6-5.81.57326.86010.0-51.2-20.81.04826.3805.0-68.3-35.10.61801005.0-48.6-27.30.61801505.0-15.2-10.50.61803005.012.810.20.6180计算结果分析：（1）快速放气持续约80s，瓶内压力从45MPa降至5MPa；（2）气体最低温度出现在放气结束时刻，约为-68.3℃；（3）瓶壁最低温度约为-35.1℃，出现在放气结束后约10s，高于设计温度下限-40℃，满足瓶体材料温度要求；（4）低温持续时间较短，放气结束后约2分钟气体温度回升至-15℃以上，约5分钟恢复至环境温度。5温度-压力耦合安全核验5.1温度阈值符合性核验根据GB/T35544-2025及GB50516-2010规定，储氢瓶工作温度范围为-40℃~+85℃。本计算结果与标准要求对比核验如下：核验项目计算结果标准限值符合性充氢过程气体最高温度61.2℃≤85℃符合充氢过程瓶壁最高温度52.1℃≤85℃符合放氢过程瓶壁最低温度-35.1℃≥-40℃符合缠绕层最高温度50.5℃≤120℃（Tg）符合密封件最高工作温度55℃≤80℃符合阀门最低工作温度-68.3℃≥-70℃符合注：瓶内气体最低温度-68.3℃虽低于-40℃，但气体不直接承受结构应力，且低温持续时间短（约30s内回升至-40℃以上），不影响瓶体结构安全；阀门、密封件及管道直接接触低温气体，设计选型按-70℃考虑。5.2温度-压力耦合强度校核最高工作温度85℃工况下，铝合金内胆强度下降约5%，碳纤维强度基本不受温度影响；最低温度-40℃工况下，铝合金强度略有提升，不会发生低温脆断。考虑温度对材料性能的影响，对45MPa工作压力下的筒体应力进行校核：σ_θ=(p·d_i)/(2·δ_eff)=(45×400)/(2×17)=529.4MPa(5-1)缠绕层环向应力：按载荷分配，缠绕层承担约85%内压，环向应力：σ_θ,wrap=(0.85×45×400)/(2×12)=637.5MPa<1100MPa（设计许用应力）(5-2)内胆周向应力（自紧后残余压应力状态）：工作压力下内胆周向应力约135MPa，小于6061-T6在-40℃下的许用应力196.7MPa，满足强度要求。纤维应力比=4900/637.5=7.69≥2.25，满足GB/T35544标准要求。5.3超压与热安全校核若充氢结束后瓶内气体温度升至85℃（极端工况），压力升高值计算：p_2=p_1·(Z_1·T_2)/(Z_2·T_1)=45×(1.298×358.15)/(1.385×332.85)=45.3MPa(5-3)即满充状态下温度升高60℃，压力仅升高约0.3MPa，不会导致超压。安全阀整定压力取52MPa，超压泄放能力满足火灾等极端工况要求。氢安全校核：氢气最小点火能仅0.019mJ，储氢瓶组区域设置可燃气体检测报警仪，报警浓度设定为25%LEL（即1%体积浓度），联锁启动紧急切断阀和事故排风系统，通风次数不小于12次/h，可有效防止氢气积聚爆炸风险。6结果分析与工程指导意义6.1主要计算结果汇总充放氢过程温度-压力耦合计算关键结果汇总如下：参数数值备注单瓶有效储氢量（45MPa，20℃）约2.42kg可用储氢量约1.95kg充氢过程气体最高温度61.2℃-40℃预冷+300s充氢充氢过程瓶壁最高温度52.1℃充氢结束时刻放氢过程气体最低温度-68.3℃80s快速放气至5MPa放氢过程瓶壁最低温度-35.1℃放气结束后10s温度恢复至环境温度时间约30min充氢后自然冷却温度异常报警阈值80℃/-50℃高/低温报警安全联锁切断阈值85℃/-60℃紧急切断6.2工程设计指导建议基于本次计算结果，对储氢瓶组系统设计和运行提出以下指导建议：（1）充氢工艺：必须采用-40℃氢气预冷技术，充氢时间不应短于180s，推荐采用阶梯升压充氢策略，先慢后快控制充氢速率，避免前期充氢过快导致温升过高；（2）放氢工艺：正常供氢时放氢速率不应超过10g/s（单瓶），快速卸压放氢时间不应短于60s，避免放气过快导致瓶壁温度低于-40℃；（3）温度监测：每只储氢瓶应在瓶肩、瓶中、瓶尾三个位置设置温度传感器，瓶组出口管道设置温度传感器，实现全范围温度监控；（4）材料选型：瓶阀、管道及密封件设计温度应按-70℃~+85℃选取，推荐阀体采用316L不锈钢，密封件采用耐低温Kalrez或改性