版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1/1气液分离装置的建模与仿真第一部分气液分离装置的建模方法和基础 2第二部分传质模型的建立和求解技术 4第三部分流动模型的构建和求解方法 7第四部分状态方程和物性参数的选取 10第五部分分离效率评估与影响因素分析 11第六部分优化模型的建立与求解策略 14第七部分仿真结果验证与模型精度评价 18第八部分气液分离装置仿真在设计和优化中的应用 20
第一部分气液分离装置的建模方法和基础关键词关键要点主题名称:流体动力学模型
1.气液分离装置的核心在于流体动力学模型,它描述了流体在装置中的流动特性。
2.模型考虑流体粘度、密度、速度分布、湍流效应和压力降等因素。
3.模型可以预测气液两相流在分离器内的分布、流动模式和分离效率。
主题名称:传质模型
气液分离装置的建模方法和基础
导言
气液分离装置广泛应用于石油、化工、制药等行业,其设计和优化对工艺效率和经济效益至关重要。气液分离过程的建模与仿真是实现装置设计和优化不可或缺的手段,为预测和分析装置性能提供理论依据。
气液分离装置的建模方法
气液分离装置的建模方法主要分为以下几类:
1.粒子方法
*离散相模型(DEM):将气液两相视为由离散粒子组成的集合,通过求解粒子运动方程来模拟流场。
*格子玻尔兹曼方法(LBM):将流体视为由若干个小速度格子组成的连续介质,通过求解格子碰撞-传播方程来模拟流场。
2.连续相方法
*欧拉-拉格朗日法:将气相视为连续介质,将液相视为离散粒子,通过求解连续相动量方程和离散相运动方程来模拟流场。
*欧拉-欧拉法:将气液两相都视为连续介质,通过求解两相连续性方程、动量方程和能量方程来模拟流场。
3.混合方法
*混合欧拉-拉格朗日法:将部分相视为连续介质,部分相视为离散粒子,结合欧拉-拉格朗日法和欧拉-欧拉法进行建模。
*混合格子玻尔兹曼方法:将部分流体区域视为离散粒子,部分区域视为连续介质,结合DEM和LBM进行建模。
气液分离装置的建模基础
1.流场求解
*质量守恒方程:表示流体质量的守恒。
*动量守恒方程:描述流体受到力的影响产生的运动。
*能量守恒方程:描述流体能量的变化。
2.相互作用模型
*曳力模型:描述气液两相之间的相互作用引起的阻力。
*质量传递模型:描述气液两相之间的物质传递。
*热传递模型:描述气液两相之间的热量传递。
3.边界条件
*入口边界条件:指定流体在入口处的速度、温度等参数。
*出口边界条件:指定流体在出口处的压力、温度等参数。
*壁面边界条件:指定流体与壁面之间的相互作用,如无滑移条件或滑移条件。
4.数值方法
*有限差分法(FDM):将流场划分为离散网格,并在网格上求解控制方程。
*有限体积法(FVM):将流场划分为一系列控制体,并在控制体上求解控制方程。
*有限元法(FEM):将流场划分为一系列单元,并在单元上求解控制方程。
5.模型验证与校正
*模型验证:通过与实验数据对比验证模型的有效性。
*模型校正:根据验证结果调整模型参数,以提高模型精度。
结论
气液分离装置的建模与仿真是一项复杂的工程技术,需要综合考虑流场求解、相互作用模型、边界条件、数值方法和模型验证等因素。通过合理选择建模方法和掌握其基础,可以建立准确可靠的模型,为气液分离装置的设计和优化提供有力支撑。第二部分传质模型的建立和求解技术传质模型的建立和求解技术
1.传质模型的建立
传质模型是描述物质在气液分离装置中传质过程的数学模型。根据系统特性和传质机理,可以建立不同的传质模型:
-平衡模型:假设气液间存在瞬时平衡,忽略传质速率的限制。
-非平衡模型:考虑传质速率的限制,将传质分为各个独立的步骤。
-准平衡模型:兼顾平衡模型和非平衡模型的优点,假定在局部区域内存在平衡,但整体传质过程存在速率限制。
2.传质模型的求解技术
传质模型的求解技术主要包括:
2.1解析法
对于一些简单的传质模型,可以通过解析方法直接求解。例如,对于平衡模型中的简单塔板,可以使用解析公式计算塔板效率。
2.2数值法
对于复杂传质模型,通常需要使用数值方法求解。常用的数值方法包括:
-网格法:将传质过程离散化为一系列单元格,每个单元格内求解传质方程。
-有限差分法:将传质方程离散化为差分方程,逐点迭代求解。
-有限元法:将传质区域离散化为单元网格,使用加权余量法求解单元内传质方程。
2.2.1隐式方法
隐式方法在求解时,将未知数移到下一个时间步或网格点,通过求解联立方程组得到未知数。隐式方法的优点是无条件稳定,可以采用较大时间步长或网格尺寸。
2.2.2显式方法
显式方法在求解时,将未知数留在当前时间步或网格点,通过直接代入方式得到未知数。显式方法的优点是计算简单,但不稳定,需要采用较小时间步长或网格尺寸。
3.常见传质模型
3.1平衡模型
平衡模型适用于气液传质速率很大的情况,其传质方程为:
```
Y=Hx
```
其中:
-Y:液相组分摩尔分数
-x:气相组分摩尔分数
-H:平衡常数
3.2非平衡模型
非平衡模型适用于气液传质速率较小的情况,其传质方程为:
```
G(Y,x)=K(Y-Hx)
```
其中:
-G(Y,x):气相传质通量
-K:传质系数
-H:平衡常数
3.3准平衡模型
准平衡模型使用平衡常数定义局部平衡,并使用传质系数描述整体传质速率,其传质方程为:
```
d(Y-Hx)/dt=K(Y-Hx)
```
其中:
-Y:液相组分摩尔分数
-x:气相组分摩尔分数
-H:平衡常数
-K:传质系数
-t:时间
4.传质模型的验证
传质模型建立后,需要对其精度进行验证。常用的验证方法包括:
-与实验数据进行比较
-使用已知解对模型进行仿真
-分析模型的敏感性第三部分流动模型的构建和求解方法关键词关键要点【连续相流动的建模】:
1.采用欧拉-拉格朗日方法,分离相流体和颗粒相流体。
2.利用连续介质模型描述分离相流体的宏观行为,应用动量、能量和质量守恒方程。
3.考虑分离相的湍流特性,采用雷诺平均纳维-斯托克斯方程(RANS)模型计算湍流流场。
【分散相流动的建模】:
流动模型的构建和求解方法
1.模型构建
气液分离器内的流动模型基于质量守恒和动量守恒定律构建,考虑流体不可压缩、恒温条件。模型中包含以下守恒方程:
*连续性方程:描述流体质量守恒。
*动量方程:描述流体动量守恒,包括对流项、扩散项、压力梯度项和重力项。
2.求解方法
模型求解采用有限体积法,将计算域离散为有限体积,并在每个体积上应用守恒方程。求解过程涉及以下步骤:
2.1压力的计算
*SIMPLE算法:压力校正方程(PPE)求解压力,该方程从动量方程中导出,使用速度和压力间的关系。
*PIMPLE算法:半隐式压力校正与隐式求解耦合,提高了求解稳定性。
2.2速度的计算
*MomentumInterpolation(SIMPLEC)算法:使用压力校正结果更新速度。
*修正对流项(Rhie-Chow插值):消除压力-速度耦合引起的奇偶数网格交错问题。
2.3湍流模型的求解
*k-ε湍流模型:广泛用于模拟湍流,基于湍动能(k)和耗散率(ε)的概念。
*k-ω湍流模型:可更好地模拟近壁面湍流,基于湍动能(k)和比耗散率(ω)的概念。
2.4界面捕捉方法
*LevelSet方法:将界面表示为标量函数ϕ,ϕ>0表示液体区域,ϕ<0表示气体区域。
*VolumeofFluid(VOF)方法:将流体表示为占有计算体积的体积分数,0表示气体,1表示液体。
3.特殊求解技术
*BoundaryConditionExtrapolation(BCE)算法:用于处理出口边界,避免错误的压力梯度干扰计算。
*CoordinatedAdaptiveRefinement(CAR)技术:自适应网格细化,在高梯度区域局部加密网格,提高精度。
*Multigrid方法:多尺度求解技术,将不同频率分量的问题逐级分解求解,提高效率。
4.模型验证和精度分析
*网格无关性研究:验证解决方案随网格细化趋于稳定。
*实验或数值基准比较:将模拟结果与实验数据或其他数值模型的结果进行比较,评估模型精度。
*敏感性分析:研究模型对输入参数变化的敏感性,确定关键参数。第四部分状态方程和物性参数的选取关键词关键要点【状态方程的选择】:
1.理想气体状态方程适用于低压、低温气体,但对于高压、高密度气体或液相,其精度较低。
2.范德瓦尔斯方程考虑了分子间相互作用,但对于极性分子或非理想气体,其精度有限。
3.彭-罗宾逊方程是相对准确的状态方程,适用于广泛的压力、温度和组分范围,但在临界点附近精度有所下降。
【物性参数的选取】:
状态方程和物性参数的选取
选择合适的热力学状态方程和物性参数对于气液分离装置的建模和仿真至关重要。
状态方程
状态方程描述了流体在给定的温度、压力和体积下的热力学行为。对于气液分离应用,常用的状态方程包括:
*理想气体方程:适用于低压、低温条件。
*Peng-Robinson(PR)方程:适用于中压、中温条件,能够预测液相和气相的相行为。
*Soave-Redlich-Kwong(SRK)方程:与PR方程类似,但对非极性流体更加准确。
*CubicPlusAssociation(CPA)方程:考虑了流体之间的缔合效应,适用于处理极性流体。
状态方程的选择取决于流体的物性、操作条件和建模精度要求。
物性参数
物性参数是状态方程中用于表征流体热力学性质的常数。对于气液分离装置,需要考虑以下主要物性参数:
*临界压力和温度:流体从气相转变为液相的条件。
*分子量:流体的质量与一个质子质量之比。
*极性因子:表征流体分子之间的偶极相互作用。
*缔合因子:表征流体分子之间的缔合能力。
*二氧化碳因子:用于表示二氧化碳在流体混合物中的比例。
物性参数的选择
物性参数可以通过实验测量获得,也可以从热力学数据库或文献中获取。对于非常见的流体,可能需要使用估算方法来确定物性参数。
选择物性参数时,需要考虑以下因素:
*流体的纯物质性质:可以从热力学数据库中获得。
*混合物性质:可以从混合规则或经验方程中估算,例如Kay's方法或Peng-Robinson混合规则。
*温度和压力效应:物性参数通常会随着温度和压力的变化而变化,需要考虑这些效应。
*建模目的和精度:对于概念性建模,可以使用简化的物性参数;对于高精度建模,需要使用更准确的物性参数。
通过仔细选择状态方程和物性参数,可以确保气液分离装置的建模和仿真结果的准确性和可靠性。第五部分分离效率评估与影响因素分析关键词关键要点分离效率评估
1.分离效率是衡量气液分离装置性能的重要指标,通常用分离系数表示,反映了装置将气体和液体有效分离的能力。
2.分离系数受多种因素影响,包括装置设计、操作条件和进料特征,如流速、流型、液滴尺寸和气泡分布。
3.评估分离效率常用的方法包括分离系数法、流体动力学模型和实验测量,其中实验测量提供最直接和准确的结果。
影响分离效率的因素
1.气液比:高气液比会导致气泡破裂和液滴夹带,降低分离效率。
2.液滴尺寸:大液滴容易被气流带走,小液滴则更容易被捕集,因此液滴尺寸对分离效率至关重要。
3.气泡分布:均匀的气泡分布有利于气液接触,提高分离效率。反之,不均匀分布会导致死区和分离不充分。分离效率评估与影响因素分析
分离效率评价指标
气液分离装置的分离效率通常通过以下指标评价:
*分离因子(SF):分离器入口气相组分浓度与出口气相组分浓度的比值,衡量气相中目标组分的富集程度。
*净化因子(PF):分离器入口液相组分浓度与出口液相组分浓度的比值,衡量液相中目标组分的去除程度。
*分离度(SR):分离因子与净化因子的几何平均值,综合衡量分离装置的气液分离能力。
影响分离效率的因素
影响气液分离装置分离效率的因素主要包括:
1.分离器类型
不同类型的气液分离器具有不同的分离机理,因此分离效率也不同。常见的分离器类型有旋流分离器、湍流沉降分离器、纤维束分离器等。
2.流动参数
气液两相的流速、流型和湍流度对分离效率有较大影响。一般来说,气液流速越高,分离效率越低;湍流度越高,分离效率越高。
3.物理性质
气液两相的密度、黏度和表面张力等物理性质也会影响分离效率。密度差越大,分离效率越高;黏度越大,分离效率越低;表面张力越大,分离效率越高。
4.雾滴尺寸
雾滴尺寸是影响分离效率的关键因素。雾滴尺寸越小,分离难度越大。因此,减小雾滴尺寸有利于提高分离效率。
5.操作条件
分离器的工作压力、温度和进料组分组成等操作条件也会影响分离效率。通常情况下,压力越高,分离效率越高;温度越高,分离效率越低;进料组分中目标组分浓度越高,分离效率越低。
6.设备结构参数
分离器的尺寸、形状和内部结构等结构参数也会影响分离效率。例如,分离器高度增加,分离效率提高;分离器截面积减小,分离效率提高。
影响因素分析
1.分离因子与气液流速的关系
分离因子随气液流速的变化而变化。在低流速下,分离因子较低,随着流速增加,分离因子升高,达到一定流速后,分离因子不再变化。
2.分离因子与雾滴尺寸的关系
分离因子与雾滴尺寸呈负相关关系。雾滴尺寸越大,分离因子越低;雾滴尺寸越小,分离因子越高。
3.分离因子与密度差的关系
分离因子与气液两相的密度差呈正相关关系。密度差越大,分离因子越高。
4.分离度与气液比的关系
分离度随气液比的变化而变化。在低气液比下,分离度较高,随着气液比增加,分离度下降,达到一定气液比后,分离度不再变化。
5.分离效率与设备高度的关系
分离效率随设备高度的增加而提高。这是因为设备高度增加,气液两相停留时间延长,分离更加充分。
6.分离效率与设备直径的关系
分离效率随设备直径的减小而提高。这是因为直径减小,气液两相流动速度增加,湍流度增加,分离效果增强。第六部分优化模型的建立与求解策略关键词关键要点目标函数的制定
1.确定优化目标,如最大化分离效率、最小化能耗或投资成本。
2.基于目标函数构建数学模型,考虑气液分离装置的实际运行条件和过程参数。
3.采用合适的优化算法,如线性规划、非线性规划或混合整数规划,以求解目标函数。
变量的选取与约束条件的设定
1.识别影响气液分离性能的关键变量,包括气液比、流量、压力和温度。
2.设置合理的约束条件,如气液两相的物性参数、管道规格和设备容量。
3.考虑工艺限制和安全要求,以确保优化方案的可行性和安全性。
优化算法的选择
1.根据目标函数和变量约束条件,选择合适的优化算法。
2.考虑算法的收敛速度、计算复杂度和鲁棒性。
3.评估不同算法的性能,并选择最能满足优化需求的算法。
参数估计与敏感性分析
1.利用实验数据或理论模型估计优化模型中的参数。
2.进行敏感性分析,了解参数变化对优化结果的影响。
3.识别关键参数,并优化这些参数以获得最佳的分离性能。
多目标优化
1.考虑气液分离装置的多个优化目标,如分离效率、能耗和投资成本。
2.构建多目标优化模型,考虑不同目标之间的权重和权衡。
3.采用多目标优化算法,如NSGA-II或MOEA/D,以求解多目标问题。
鲁棒优化
1.考虑到气液分离装置在实际运行中可能存在的参数不确定性和扰动。
2.构建鲁棒优化模型,使优化方案对不确定性具有鲁棒性。
3.采用鲁棒优化算法,如鲁棒线性规划或鲁棒非线性规划,以求解鲁棒优化问题。气液分离装置优化模型的建立与求解策略
#优化模型的建立
优化模型的建立涉及以下关键步骤:
1.模型选择:选择适合气液分离装置的物理模型,如平衡模型、非平衡模型或CFD模型。
2.目标函数定义:确定优化目标,如最大化分离效率、最小化能耗或其他性能指标。
3.决策变量识别:确定影响目标函数的决策变量,如进料流量、进料温度、塔径和塔板数。
4.约束条件定义:建立约束条件以确保模型的物理可行性和操作安全性,如平衡关系、压力降限制和设备尺寸限制。
#求解策略
优化模型建立后,需要采用合适的求解策略求解最佳决策变量。常见求解策略包括:
1.线性规划(LP):适合处理线性目标函数和决策变量的优化问题。用于优化物料平衡和能量平衡等线性问题。
2.非线性规划(NLP):适用于含非线性目标函数或决策变量的优化问题。用于优化具有复杂平衡关系或非线性热力学性质的气液分离装置。
3.混合整数规划(MILP):当决策变量包含整数时使用。例如,塔板数或设备类型。
4.启发式算法:当模型过于复杂或难以求解时使用,如遗传算法、模拟退火或粒子群优化。
#模型求解过程
模型求解过程通常包括以下步骤:
1.模型参数估计:根据实验数据或工艺知识估计模型参数,如进料组分、相平衡常数和热传递系数。
2.模型验证:将模型预测结果与实际操作数据进行比较以验证模型的准确性。
3.求解优化模型:使用所选求解器求解优化模型以获得最佳决策变量值。
4.敏感性分析:研究决策变量对目标函数和约束条件的影响,以了解模型的鲁棒性和关键参数。
5.模型优化:根据敏感性分析结果,进一步优化模型结构或求解策略以提高模型的精度和效率。
#实施和监控
优化模型求解后,需要将其实施到工艺设计和控制中。监控和评估模型性能至关重要,以确保其持续有效性和准确性。
1.模型实施:将优化结果整合到气液分离装置的设计、操作或控制系统中。
2.模型监控:定期监测模型预测与实际操作数据的比较,以确保模型准确性。
3.模型更新:随着工艺条件或设备性能的变化,更新模型参数或结构以保持其相关性和准确性。
#结论
优化模型的建立与求解策略对于气液分离装置的优化设计和操作至关重要。通过选择合适的模型、求解器和实施策略,工程师可以获得最佳决策变量值,从而改善分离效率、降低能耗并增强工艺可靠性。持续监控和更新模型对于确保其有效性并满足不断变化的工艺需求至关重要。第七部分仿真结果验证与模型精度评价关键词关键要点主题名称:模型精度评价方法
1.采用统计指标,如均方根误差(RMSE)、平均绝对误差(MAE)和决定系数(R^2),来量化模型预测结果与实际观测值之间的差异。
2.利用敏感性分析,考察模型输入参数的扰动对预测结果的影响,评估模型的鲁棒性。
3.结合定性和定量评价,通过交叉验证、拟合优度检验和专家判断,综合考虑模型的准确性和适用性。
主题名称:仿真结果与实验数据的比较
仿真结果验证与模型精度评价
一、仿真结果验证
仿真结果验证通过对比仿真结果与实际测量数据或公认模型的结果来进行,目的是确保模型捕获了系统关键特性并能够准确预测实际操作。
1.与实际测量数据对比
*定性验证:检查仿真结果与测量数据之间的总体趋势和形状相似性。
*定量验证:使用统计指标(如均方根误差、相关系数等)量化仿真结果与测量数据的差异。
2.与公认模型对比
*定性验证:比较不同模型的预测结果,寻找关键趋势和行为的一致性。
*定量验证:通过统计指标评估不同模型预测结果之间的差异和相似性。
二、模型精度评价
1.误差分析
*绝对误差:仿真结果与实际测量值或公认模型预测值之间的绝对差异。
*相对误差:绝对误差与实际测量值或公认模型预测值之比。
*最大相对误差:仿真结果与实际测量值或公认模型预测值之间的最大相对误差。
2.统计指标
*均方根误差(RMSE):仿真结果与实际测量值或公认模型预测值之间平方差的平均值开平方根。
*相关系数(R):仿真结果与实际测量值或公认模型预测值之间的相关程度,介于-1到1之间。
*均方误差(MSE):仿真结果与实际测量值或公认模型预测值之间平方差的平均值。
3.敏感性分析
*通过改变模型输入参数,评估模型输出对这些参数的敏感性。
*有助于识别模型中关键不确定性并指导进一步的模型优化。
4.建模不确定性量化
*利用概率方法,考虑模型输入参数和假设的不确定性,对模型输出的置信区间进行量化。
*有助于评估预测结果的鲁棒性和可靠性。
三、提高模型精度
为了提高模型精度,可以采取以下措施:
*收集高精度测量数据:作为验证和校准模型的基础。
*选择合适的建模工具:匹配系统的复杂性和建模目标。
*仔细校准模型参数:根据测量数据或公认模型调整模型参数。
*进行全面验证:使用多种验证方法确保模型精度。
*持续优化模型:通过敏感性分析和不确定性量化不断提高模型的预测能力。第八部分气液分离装置仿真在设计和优化中的应用气液分离装置仿真在设计和优化中的应用
气液分离装置仿真技术在气液分离过程的设计和优化中发挥着至关重要的作用。通过构建和求解装置模型,仿真可以预测其行为,识别潜在问题并制定优化策略。
设计阶段的应用
在设计阶段,仿真可用于:
*评估不同方案的可行性:比较不同设计方案的性能,确定最佳选择。例如,可以仿真不同塔盘数和直径的蒸馏塔,以确定实现特定分离目标的最佳配置。
*优化设备参数:确定最佳操作参数,如塔板间距、再沸器热负荷、回流比等,以最大化装置效率。通过仿真,可以预测这些参数的变化对分离性能的影响,并确定其最佳值。
*识别潜在问题:识别和解决设计中的潜在问题,如液泛、压降过大、塔内分布不均等。仿真可以预测这些问题的发生,并帮助设计人员采取措施加以避免或减轻。
*减少试车时间和成本:通过仿真,可以提前预测装置的启动和操作行为。这有助于优化试车计划,减少试车时间和成本。
优化阶段的应用
在优化阶段,仿真可用于:
*提高装置效率:识别和纠正影响装置效率的因素,如塔内结垢、催化剂失活、设备故障等。通过仿真,可以量化这些
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 仓库消防方案
- 电影院行业发展规模预测与投资潜力研究报告
- 修复磨损或部分损坏的发动机行业三年发展预测分析报告
- 关于专利地图的法律咨询行业发展趋势研判及战略投资深度研究报告
- fy2023初中英语个性化初一九月考卷简单教师版
- 英语文摘中英对照File
- 2024年仁爱版三年级数学(上册)模拟试卷及答案(各版本)
- 2024年非磁介质复制品合作协议书
- 有关知识产权的顾问服务行业发展趋势研判及战略投资深度研究报告
- 2024年模糊控制变频空调器项目合作计划书
- 2024年绿化工程合同参考样本(二篇)
- 2.2 做更好的自己 课件-2024-2025学年统编版道德与法治七年级上册
- 地球在宇宙中的位置 课件-2024-2025学年七年级地理上学期(2024)人教版
- 2024年白山长白朝鲜族自治县公开招聘“乡村振兴专员”77名工作人员(高频重点复习提升训练)共500题附带答案详解
- 广告行业市场调研合同
- 2024陕西西安市轨道交通集团限公司运营分公司招聘易考易错模拟试题(共500题)试卷后附参考答案
- 初中语文 八下第五单元整体公开课一等奖创新教学设计
- 2024年译林英语三年级上册Unit1 第3课时 (教学课件)Wrap-up time
- 2024秋季开学第一课主题班会教案设计巴黎奥运精神3篇
- 非暴力沟通改善亲子关系
- 汽修合伙人合同协议书
评论
0/150
提交评论