




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
结晶器的选择与设计在现代结晶器的设计中考虑结晶动力学设计产品的尺寸,尺寸产品的质量。结晶器工艺设计
1.确定操作模型
(1).连续型
(2).间歇型
2.产生过饱和度的方法冷却、蒸发、反应、其它
3。结晶过程质量衡算质量、热量、粒数衡算(略)颗粒衡算与质量、热量衡算的关系
结晶器的选择与设计在现代结晶器的设计中3.设计流程操作模型的选择连续:大型生产一般产生的晶体的尺寸分布较宽可用一些控制手段来完成大晶体窄分布的产品。间歇:相对小的生产能力(50T/天)一般可产生分布较窄的晶体也需要适当的控制方法。操作过程的选择冷却:溶解度随温度变化较大的体系,适宜用冷却的方法。蒸发:如果溶解度随温度变化不大的体系可用蒸发法。在溶解度高时,为提高回收率用蒸发的方法。反应:能产生沉淀的反应体系。其它:溶剂萃取,高压结晶等。3.设计流程结晶器类型的选择DTB结晶器(drafttubeandbaffle)图5.7(导流筒、挡板式)特征:固体悬浮较好,可带有结晶排除系统和产品颗粒分级系统,实现晶体粒度的控制。结晶器类型的选择DTB结晶器(drafttubeandFC(forcedcirculation)强制循环图5.115.175.18强制循环结晶器可用于很多过程,可根据不同的过程来调节,其也可设置粒度控制系统。FC(forcedcirculation)强制循环第六章-结晶器设计课件结晶器的计算生产能力和产品质量,和粒度要求计算要求的基础数据与模型溶解度数据溶液的基本性质和颗粒的基本性质如密度、形状系数等结晶成长速率模型与成核动力学模型溶液的初始浓度C0结晶器的计算生产能力和产品质量,和粒度要求计算:
MT—悬浮密度
C1—溶液的终止浓度,根据最后排出体系的温度而定,即C1=f(Tf)
VW—蒸发水量
V0—溶液的进料体积流率,此参数对多组分系统,可根据母液中杂质浓度而定,对单组分系统MT是一个操作参数,根据结晶器的流动状态而定。因此,可用上式计算体积、流率。当结晶器操作于MSMPR模型结晶器体积:进料流率*停留时间这样设计的结晶器比较简单,但因为是MSMPR的假设下,而以其结果很难适用于大型结晶器。计算:假设在稳定连续操作的结晶过程中,其产品的颗粒的累积尺寸分布可用Rosion-Rammler方程表示:
设计图及设计程序r—相对产品颗粒个数的累积分布,当所有的晶体考虑在内r=1m—改良的均一系数,不同的m值表示不同的分布函数,假设在一个过程中的产品其m值不变,也不随颗粒的尺寸而变。
L—晶体的尺寸。
L*—参考尺寸,其定义为在此尺寸下,r值为0.3679与m无关。
x—无因次晶体尺寸,x=L/L*假设在稳定连续操作的结晶过程中,其产品的颗粒的累积尺寸分布可在用方程(1)表示积累产品晶体分布时,颗粒的粒数密度(即在一定的颗粒尺寸范围内的晶体个数)为
(2)
在用方程(1)表示积累产品晶体分布时,颗粒的粒数密度(即在一因此,在此尺寸范围内的晶体质量为(3)
—产品晶体的质量分布kg/m·h—晶体密度kg/m3—晶体的体积形状系数
—以个数为基准的晶体生产速率,它等于晶体的成核速率#/h因此,在此尺寸范围内的晶体质量为如果以产品中最大质量的晶体作为晶体尺寸来表示,这种晶体尺寸可以从方程(3)得出即(4)
注意:当m=1时,xd=3,其分布式和MSMPR的结果一致。当m=0时,产品晶体为同一尺寸,表示一个连续操作的带有理想产品分级的操作过程,对方程(3)进行积分,得到以质量为基础的生产速率。
(5)如果以产品中最大质量的晶体作为晶体尺寸来表示,这种晶体尺寸可当m=1是全部结晶器内成核速率,如果用V表示结晶器体积,成核速率表示为B,即(6)方程(6)变成这一结果和MSMPR的结果相同当m=1是全部结晶器内成核速率,如果用V表示结晶器体积以上分析是从产品晶体部分来考虑,下面再从结晶器内部的晶体变化来考虑。在稳定操作状态下,细小的晶体或者从外部加入(加晶种)或者来源于系统内的成核,假设在操作过程中成长的晶体没有破碎或聚并。因此,在结晶器内一定尺寸的晶体的晶数密度与产品中大于此晶体尺寸的积累尺寸分布成正比,也就是说,在这种假设下,在结晶器内的晶体密度分布正比于方程(1)的r。如果把晶种(成核)在dx范围内的个数表示为n0dx(注意:n0为核的粒数密度n0,在无因次晶体尺寸x的颗粒数为以上分析是从产品晶体部分来考虑,下面再从结晶器内部的晶体变化因此,在结晶器内全部晶体重量为
从另一方面讲,在结晶器内晶体的全部个数也可以表示为(8)Θ—晶体的平均停留时间从晶体的平均停留时间的定义
Θ—平均停留时间一定尺寸的晶体个数一定尺寸晶体的停留时间(7)因此,在结晶器内全部晶体重量为(7)等号两边同时等于结晶器内晶体的总个数(相对值,因为r为相对分布),从晶体生长的角度我们有
或(9)
由方程(8)和(9)
粒数密度(10)
等号两边同时等于结晶器内晶体的总个数(相对值,因为r为相对分把方程(10)代入晶体总质量方程(7)
由生产速率方程(方程(5))和方程(11)我们有
(11)(12)把方程(10)代入晶体总质量方程(7)(11)(12)从结晶器内的悬浮密度考虑,如果我们用ε表示溶液的含量分数,结晶器内晶体的质量可表示为:
最后得到(13)(14)从结晶器内的悬浮密度考虑,如果我们用ε表示溶液的含量分数,结分析方程(14)此方程给出生产强度(P/ρCV)悬浮密度(1-ε),晶体尺寸L(或X),成长速率,成核速率方程(6),从产品角度分析而得到的粒数,质量衡算。方程(14)是从结晶器内的结晶过程而得到的粒数,质量衡算方程,这两个方程必须建立一定的联系,从此而得到一个结晶器的设计图。设(15)
(16)(6)分析方程(14)(6)
用Ld代替L*方程(6)可表示为方程(14)可表示为
(17)(18)用Ld代替L*(17)(18)作图把方程(17)取对数在一定的Ld下,与为一条直线,斜率为1这条直线取决于要求的产品粒度Ld和相应m值为特征。同一直线下,表示在一定的m值下的同一粒径,这条直线方程做特征线,就是说,这些平行线表示不同的粒径。另一方面,在一定的生产能力下,由方程(17),Ld和Bkv立方根成正比,因此,Bkv的轴可以和Ld轴平行而置。作图把方程(17)取对数第六章-结晶器设计课件作一个表示Ld,P/ρCV,Bkv的关系图,m=1作出特征线并确定Ld在Ld轴上的值,当m=1时,I3可被计算为0.225
当P/ρCV=0.01,Bkv=1×108,Ld=7.631×10-4P/ρCV=0.01,Bkv=1×109,Ld=3.542×10-4图8用纵坐标表示P/ρCV(对数值)用一横坐标表示Bkv(对数值)另一横坐标表示Ld轴(对数值)在Bkv轴上找到1×108和P/ρCV=0.01作一直线通过此点,但斜率等于1标此直线与Ld轴交点为7.631×10-4与上相同做出,P/ρCV=0.01,Bkv=1×109不同的m值做图。作一个表示Ld,P/ρCV,Bkv的关系图,m=1作出特征线第六章-结晶器设计课件
按道理讲,不同的m值,其P/ρCV、Bkv和Ld应遵循不同的规律,可以完全做成另一个图,然而,要是用一张图表示即存在的的信息总量较为简单。在特征方程中,如果我们把P/ρCV考虑成一个常数,Bkv=k(m)1/Ld3这个方程表明,当P/ρCV不变时,Bkv与Ld3成反比,如不同的m值,它们的关系只改变其系数但比例关系不变,也就是说这条直线的位置有些变化,但变化趋势是一样的。这样我们就可以用一个Ld与Bkv之间的关系。按道理讲,不同的m值,其P/ρCV、Bkv和Ld应遵循不同例如:m=1,P/ρCV=0.01,Bkv=1×108,Ld=7.631×10-4m=2,P/ρCV=0.01,Bkv=1×108,Ld=6.013×10-4在例1图上作出m=2的Bkv轴和Ld轴做图:选画一条辅助线画一条横线为当m=2时的Bkv线从m=1的Bkv轴上找到Bkv=1×108通过此点作垂直线交m=2于H点作平行于辅助线的直线通过H点交P/ρCV=0.01的横线于G
例如:m=1,P/ρCV=0.01,Bkv=1×108,Ld通过G作斜率等于1的特征线,这条特征线的Ld必须等于6.013×10-4从m=1的轴上找到6.013×10-4,做垂直线交于此条直线位F通过F点作出直线即为m=2得Ld轴这里强调一下,在这样的图中,如果m值不同BkV和Ld得值是不同的。到此,我们在一张图上把Ld,BkV,m,P/ρV的关系表示出来。通过G作斜率等于1的特征线,这条特征线的Ld必须等于6.01第六章-结晶器设计课件结晶器内的变化规律图由方程从产品方程可得代入如果假设为常数则
即从结晶器内的粒数与质量衡算结晶器内的变化规律图
与为直线斜率斜率为,这条线与表示与G乘积相同,所以做坐标()轴和G轴时,这个方程也一定要满足,要等于方程左面的乘积。即()或可以任选一个值而另一个值必须满足上面的方程,把这些关系标在同一张图上,成为图7与为直第六章-结晶器设计课件成核速率对晶体尺寸的影响分析在生产强度一定的条件下,即,P/ρCV为常数控制结晶器中的成核速率成核速率增加,操作点向右移动,晶体尺寸减小。反之增加。另一方面,从晶体的成长与质量横算考虑,在图中有要满足操作线,即悬浮密度与成长速率的关系。成核速率对晶体尺寸的影响分析在生产强度一定的条件下,即,P/第六章-结晶器设计课件生产速率增加生产速率的增加,在设计图上操作点上移产品的晶体尺寸可能增加,已可能降低。要看其造作是怎样提高生产能力的在原操作点左方,晶体尺寸增加。在原操作点右方,晶体尺寸减小。生产速率增加生产速率的增加,在设计图上操作点上移第六章-结晶器设计课件保持产品粒度,而增加(或降低)生产能力保持产品粒度,在设计图上是早同一条特征线,但其操作线在变化。同一条特征线,在不同的生产能力下,要求有不同的操作线。其悬浮密度必须保持恒定,以确保产品的粒度不变,因此,成长速率必须成比例增加。保持产品粒度,而增加(或降低)生产能力保持产品粒度,在设计图第六章-结晶器设计课件设计图的应用
-用试验数据估计结晶动力学数据试验可测得粒度分布,生产速率,结晶器有效体积,悬浮密度根据粒度分布,可计算粒数分布获得平均尺寸Lw(即积累尺寸分布为36.8%的晶体尺寸)和在参数n,n为(Rosin-Rammlardiagram的斜率)由n可以得到m,Lw可转化为Ld由测得的晶体尺寸分布可得到Ld,和m从图中根据悬浮密度Ld,m和生产能力确定成长速率和成核速率(如果实验中通过改变不同的参数和测量起过饱和度,可得成核速率模型)设计图的应用
-用试验数据估计结晶动力学数据试验可测得m与nde关系,和LD与LW,m的关系m与nde关系,和LD与LW,m的关系第六章-结晶器设计课件改善晶体尺寸改善晶体尺寸
已知:m,Lw,或m,Ld(产品要求)
生产能力p,生产速率计算:悬浮密度(质量衡算)设计方法:由悬浮密度和成长速率确定操作线由Ld确定特征线由两条线的交点即为系统点有特征线上读出P/ρV得到V同时得到BKV值,即如果要达到产品的粒度要求成核速率必须控制,在此成核速率不能达到时,系统操作必须调节悬浮密度或成长速率。结晶器体积的设计已知:m,Lw,或m,Ld(产品要求)结晶器体积的第六章-结晶器设计课件结晶过程的模拟(基于平仓贤理论)缺欠动力学数据(成核速率在生产上的要求与其规律没有联系)图的准确程度用完全计算的方法实现结晶器的设计,考虑到体系的成核与成长速率设计方程可以改写为可用试差求解此两个方程
结晶过程的模拟(基于平仓贤理论)缺欠计算方法已知:p,ρc,m,Ld,MT(其中可给范围值)
给一个成长速率G
计算V从方程a
计算成核速率从方程b
用成核速率模型计算成核速率比较如不匹配,重设成长速率进行系统模拟。计算方法已知:p,ρc,m,Ld,MT(其中可给范围值)
结晶器的选择与设计在现代结晶器的设计中考虑结晶动力学设计产品的尺寸,尺寸产品的质量。结晶器工艺设计
1.确定操作模型
(1).连续型
(2).间歇型
2.产生过饱和度的方法冷却、蒸发、反应、其它
3。结晶过程质量衡算质量、热量、粒数衡算(略)颗粒衡算与质量、热量衡算的关系
结晶器的选择与设计在现代结晶器的设计中3.设计流程操作模型的选择连续:大型生产一般产生的晶体的尺寸分布较宽可用一些控制手段来完成大晶体窄分布的产品。间歇:相对小的生产能力(50T/天)一般可产生分布较窄的晶体也需要适当的控制方法。操作过程的选择冷却:溶解度随温度变化较大的体系,适宜用冷却的方法。蒸发:如果溶解度随温度变化不大的体系可用蒸发法。在溶解度高时,为提高回收率用蒸发的方法。反应:能产生沉淀的反应体系。其它:溶剂萃取,高压结晶等。3.设计流程结晶器类型的选择DTB结晶器(drafttubeandbaffle)图5.7(导流筒、挡板式)特征:固体悬浮较好,可带有结晶排除系统和产品颗粒分级系统,实现晶体粒度的控制。结晶器类型的选择DTB结晶器(drafttubeandFC(forcedcirculation)强制循环图5.115.175.18强制循环结晶器可用于很多过程,可根据不同的过程来调节,其也可设置粒度控制系统。FC(forcedcirculation)强制循环第六章-结晶器设计课件结晶器的计算生产能力和产品质量,和粒度要求计算要求的基础数据与模型溶解度数据溶液的基本性质和颗粒的基本性质如密度、形状系数等结晶成长速率模型与成核动力学模型溶液的初始浓度C0结晶器的计算生产能力和产品质量,和粒度要求计算:
MT—悬浮密度
C1—溶液的终止浓度,根据最后排出体系的温度而定,即C1=f(Tf)
VW—蒸发水量
V0—溶液的进料体积流率,此参数对多组分系统,可根据母液中杂质浓度而定,对单组分系统MT是一个操作参数,根据结晶器的流动状态而定。因此,可用上式计算体积、流率。当结晶器操作于MSMPR模型结晶器体积:进料流率*停留时间这样设计的结晶器比较简单,但因为是MSMPR的假设下,而以其结果很难适用于大型结晶器。计算:假设在稳定连续操作的结晶过程中,其产品的颗粒的累积尺寸分布可用Rosion-Rammler方程表示:
设计图及设计程序r—相对产品颗粒个数的累积分布,当所有的晶体考虑在内r=1m—改良的均一系数,不同的m值表示不同的分布函数,假设在一个过程中的产品其m值不变,也不随颗粒的尺寸而变。
L—晶体的尺寸。
L*—参考尺寸,其定义为在此尺寸下,r值为0.3679与m无关。
x—无因次晶体尺寸,x=L/L*假设在稳定连续操作的结晶过程中,其产品的颗粒的累积尺寸分布可在用方程(1)表示积累产品晶体分布时,颗粒的粒数密度(即在一定的颗粒尺寸范围内的晶体个数)为
(2)
在用方程(1)表示积累产品晶体分布时,颗粒的粒数密度(即在一因此,在此尺寸范围内的晶体质量为(3)
—产品晶体的质量分布kg/m·h—晶体密度kg/m3—晶体的体积形状系数
—以个数为基准的晶体生产速率,它等于晶体的成核速率#/h因此,在此尺寸范围内的晶体质量为如果以产品中最大质量的晶体作为晶体尺寸来表示,这种晶体尺寸可以从方程(3)得出即(4)
注意:当m=1时,xd=3,其分布式和MSMPR的结果一致。当m=0时,产品晶体为同一尺寸,表示一个连续操作的带有理想产品分级的操作过程,对方程(3)进行积分,得到以质量为基础的生产速率。
(5)如果以产品中最大质量的晶体作为晶体尺寸来表示,这种晶体尺寸可当m=1是全部结晶器内成核速率,如果用V表示结晶器体积,成核速率表示为B,即(6)方程(6)变成这一结果和MSMPR的结果相同当m=1是全部结晶器内成核速率,如果用V表示结晶器体积以上分析是从产品晶体部分来考虑,下面再从结晶器内部的晶体变化来考虑。在稳定操作状态下,细小的晶体或者从外部加入(加晶种)或者来源于系统内的成核,假设在操作过程中成长的晶体没有破碎或聚并。因此,在结晶器内一定尺寸的晶体的晶数密度与产品中大于此晶体尺寸的积累尺寸分布成正比,也就是说,在这种假设下,在结晶器内的晶体密度分布正比于方程(1)的r。如果把晶种(成核)在dx范围内的个数表示为n0dx(注意:n0为核的粒数密度n0,在无因次晶体尺寸x的颗粒数为以上分析是从产品晶体部分来考虑,下面再从结晶器内部的晶体变化因此,在结晶器内全部晶体重量为
从另一方面讲,在结晶器内晶体的全部个数也可以表示为(8)Θ—晶体的平均停留时间从晶体的平均停留时间的定义
Θ—平均停留时间一定尺寸的晶体个数一定尺寸晶体的停留时间(7)因此,在结晶器内全部晶体重量为(7)等号两边同时等于结晶器内晶体的总个数(相对值,因为r为相对分布),从晶体生长的角度我们有
或(9)
由方程(8)和(9)
粒数密度(10)
等号两边同时等于结晶器内晶体的总个数(相对值,因为r为相对分把方程(10)代入晶体总质量方程(7)
由生产速率方程(方程(5))和方程(11)我们有
(11)(12)把方程(10)代入晶体总质量方程(7)(11)(12)从结晶器内的悬浮密度考虑,如果我们用ε表示溶液的含量分数,结晶器内晶体的质量可表示为:
最后得到(13)(14)从结晶器内的悬浮密度考虑,如果我们用ε表示溶液的含量分数,结分析方程(14)此方程给出生产强度(P/ρCV)悬浮密度(1-ε),晶体尺寸L(或X),成长速率,成核速率方程(6),从产品角度分析而得到的粒数,质量衡算。方程(14)是从结晶器内的结晶过程而得到的粒数,质量衡算方程,这两个方程必须建立一定的联系,从此而得到一个结晶器的设计图。设(15)
(16)(6)分析方程(14)(6)
用Ld代替L*方程(6)可表示为方程(14)可表示为
(17)(18)用Ld代替L*(17)(18)作图把方程(17)取对数在一定的Ld下,与为一条直线,斜率为1这条直线取决于要求的产品粒度Ld和相应m值为特征。同一直线下,表示在一定的m值下的同一粒径,这条直线方程做特征线,就是说,这些平行线表示不同的粒径。另一方面,在一定的生产能力下,由方程(17),Ld和Bkv立方根成正比,因此,Bkv的轴可以和Ld轴平行而置。作图把方程(17)取对数第六章-结晶器设计课件作一个表示Ld,P/ρCV,Bkv的关系图,m=1作出特征线并确定Ld在Ld轴上的值,当m=1时,I3可被计算为0.225
当P/ρCV=0.01,Bkv=1×108,Ld=7.631×10-4P/ρCV=0.01,Bkv=1×109,Ld=3.542×10-4图8用纵坐标表示P/ρCV(对数值)用一横坐标表示Bkv(对数值)另一横坐标表示Ld轴(对数值)在Bkv轴上找到1×108和P/ρCV=0.01作一直线通过此点,但斜率等于1标此直线与Ld轴交点为7.631×10-4与上相同做出,P/ρCV=0.01,Bkv=1×109不同的m值做图。作一个表示Ld,P/ρCV,Bkv的关系图,m=1作出特征线第六章-结晶器设计课件
按道理讲,不同的m值,其P/ρCV、Bkv和Ld应遵循不同的规律,可以完全做成另一个图,然而,要是用一张图表示即存在的的信息总量较为简单。在特征方程中,如果我们把P/ρCV考虑成一个常数,Bkv=k(m)1/Ld3这个方程表明,当P/ρCV不变时,Bkv与Ld3成反比,如不同的m值,它们的关系只改变其系数但比例关系不变,也就是说这条直线的位置有些变化,但变化趋势是一样的。这样我们就可以用一个Ld与Bkv之间的关系。按道理讲,不同的m值,其P/ρCV、Bkv和Ld应遵循不同例如:m=1,P/ρCV=0.01,Bkv=1×108,Ld=7.631×10-4m=2,P/ρCV=0.01,Bkv=1×108,Ld=6.013×10-4在例1图上作出m=2的Bkv轴和Ld轴做图:选画一条辅助线画一条横线为当m=2时的Bkv线从m=1的Bkv轴上找到Bkv=1×108通过此点作垂直线交m=2于H点作平行于辅助线的直线通过H点交P/ρCV=0.01的横线于G
例如:m=1,P/ρCV=0.01,Bkv=1×108,Ld通过G作斜率等于1的特征线,这条特征线的Ld必须等于6.013×10-4从m=1的轴上找到6.013×10-4,做垂直线交于此条直线位F通过F点作出直线即为m=2得Ld轴这里强调一下,在这样的图中,如果m值不同BkV和Ld得值是不同的。到此,我们在一张图上把Ld,BkV,m,P/ρV的关系表示出来。通过G作斜率等于1的特征线,这条特征线的Ld必须等于6.01第六章-结晶器设计课件结晶器内的变化规律图由方程从产品方程可得代入如果假设为常数则
即从结晶器内的粒数与质量衡算结晶器内的变化规律图
与为直线斜率斜率为,这条线与表示与G乘积相同,所以做坐标()轴和G轴时,这个方程也一定要满足,要等于方程左面的乘积。即()或可以任选一个值而另一个值必须满足上面的方
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 温州理工学院《音乐》2023-2024学年第一学期期末试卷
- 山东中医药高等专科学校《比较宪法》2023-2024学年第二学期期末试卷
- 上海体育大学《流体机械CAD》2023-2024学年第二学期期末试卷
- 武汉电力职业技术学院《食用菌栽培与加工技术》2023-2024学年第二学期期末试卷
- 商场铺位租赁合同书二零二五年
- 荒山荒地承包合同书范例
- 二零二五离婚抚养费给付标准
- 质押借款合同模板二零二五年
- 二零二五社保补偿金协议
- 私人房屋装修安全协议书
- 2025年会计政策更新试题及答案
- 2025分布式光伏工程验收标准规范
- 2024年初级药师考试历年真题试题及答案
- 浙江省杭州市萧山区高桥初中教育集团2024学年第二学期4月份素养调研九年级语文试卷题卷
- 二级造价师水利工程考试真题卷(2025年)
- 2024年云南省气象部门事业单位招聘考试真题
- 9.2法律保障生活 课件 -2024-2025学年统编版道德与法治七年级下册
- 幼儿园牛奶知识普及课件
- 公司泥头车管理制度
- 2025中美关税大战“对等关税”政策解读课件
- 《九州通医药公司应收账款管理现状、问题及对策》13000字(论文)
评论
0/150
提交评论