SPS建模操作及练习课件

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,9/22/2015,#,SPS,的建模及练习,1,SPS的建模及练习1,SPS,的程序和文件构成：,INPREP,PREPR,OUTPRP,INTRAN,TRANS,OUTTRN,INGRAF,GRAFR,OUTGRF,TPORT/*,可以从,review,文件中读取冷态结果,为方便阅读，在,INPREP/INTRAN/INGRAF,中均可包含,INCLUDE,文件。,2,SPS 的程序和文件构成：2,SPS,输入文件：,命名：,aaa.inprep aaa.intran aaa.ingraf bbb.inc,编辑：推荐用,window,系统自带的写字板编写,也可以使用,model builder,3,SPS 输入文件：3,单位：,Keyword,LIQUID,GAS,Keyword,LIQUID,GAS,AFLOW,AM3/HR,AM3/HR,MASS,KG,KG,ANGLE,DEGREE,DEGREE,MFLOW,KG/S,KG/S,AREA,CM2,CM2,POWER,KW,KW,DENSITY,KG/M3,KG/M3,PRESSURE,KPAG,BARG,DIAMETER,MM,MM,ROUGHNESS,MM,MM,EFFICIENCY,FR,FR,SHEAR.RATE,1/SEC,1/SEC,ELEVATION,M,M,SHEAR.STRESS,PA,PA,FLOW,M3/HR,M3/HR,SIG.TIME,MIN,MIN,FRACTION,FR,FR,SPEED,RPM,RPM,FUEL,M3/HR,M3/HR,SWEPT.VOLUME,CM3,CM3,HEAT.RATE,KJ/KW-HR,KJ/KW-HR,TEMPERATURE,DC,DC,HEAD,M,M,TENSILE.STRESS,GPA,GPA,HEAT.CAPACITY,KJ/KG-DC,KJ/KG-DC,TIME,MIN,MIN,HEAT.COND,KJ/HR-M-DC,KJ/HR-M-DC,TFLOW,GJOULE/D,GJOULE/D,HEAT.TRANSFER,KJ/HR-M2-DC,KJ/HR-M2-DC,TORQUE,N-M,N-M,HEAT.VALUE,KJ/M3,KJ/M3,VALVE.CG,M3/HR-KPA,M3/HR-KPA,INERTIA,KG-M2,KG-M2,VALVE.COEFF,M3/HR-KPA.5,M3/HR-KPA.5,LENGTH.HEADER,KM,KM,VALVE.FLOW,M3/HR,M3/HR,LENGTH.PIPE,KM,KM,VELOCITY,M/S,M/S,LINEPACK,M3,M3,VISCOSITY,CP,CP,WALL,MM,MM,VOLUME,M3,M3,4,单位：KeywordLIQUIDGASKeywordLIQU,编写,inprep,文件,5,编写inprep文件5,inprep,文件必须编写的内容,TITLE,(,必须首行,),LIQUID,or,GAS,ENGLISH,or,METRIC,CUSTODY,PIPEPARMS,ISOTHERMAL,THERMAL,or,TRANSTHERMAL,STATE AGA,STATE BWRS,STATE CNGA,STATE SCL,VISCOSITY (non-Newtonian),WAX,STATE TABLE,),=EQUIPMENT,设备、节点、变量、曲线等,/*,字母可以大写也可小写,6,inprep文件必须编写的内容6,编写,inprep,文件：,标题：,TITLE,，第,1,行，可以是任意文字和字符,SELECT,：选择,SPS,的其他产品模块，包括仿真模块、泄漏检测模块、培训模块等，可以不输入，默认仿真模块,-SIMULATOR,。在线模拟需要输入,RTU,传输的数据。,气体或液态：,GAS,或者,LIQUID,7,编写inprep文件：7,编写,inprep,文件：,标准状态：,格式：,CUSTODY,+PRESSURE = PREF,+TEMPERATURE = TREF,举例：,CUSTODY,+PRESSURE=1.01325,+TEMPERATURE=20,注,1,：此语句的压力为绝对压力，其余没有标注的压力均为表压。标准压力,101.325kPa,（绝压）,注,2,：,API,中以,15,为标准温度，中国国标以,20,为标准温度。,注,3,：标准状态影响标态流量。,问题：表压,5bar,，绝压,=,？,8,编写inprep文件：8,编写,inprep,文件：,自定义单位：,格式：,DEFUNITS NAME = EXPRESSION,USEUNITS KEYWORD NAME,举例：,DEFUNITS NN= M3/HR*0.0001 /*,标准流量单位采用万方,/,小时,USEUNITS FLOW NN,USEUNITS FUEL NN,DEFUNITS NN1 = AM3/HR*0.0001 /*,管态流量单位采用万方,/,小时,USEUNITS AFLOW NN1,注意,1,：,KEYWORD,见,Default units,表格。,注意,2,：,NAME,可以是任何字符，但不能与,SPS,内置语句重复。,注意,3,：,EXPRESSION,使用,SPS,的内置单位，见,Default units,表格。,问题：如何用,MPa,代替,SPS,模型的内置压力单位？,9,编写inprep文件：9,编写,inprep,文件：,输入管道参数：,格式：,PIPEPARMS,+ FRICTION COLE | NIKUR | MOODY RUF | FF,+ INITIAL PINIT,+ KNOT SPAC,+ THRM.COEFF n,举例：,PIPEPARMS,+ FRICTION COLE 0.01 /*,管道粗糙度按,10,m,+ KNOT 1 /*,距离计算步长,1km,+ INITIAL 40 /*,线路最高点起始压力,40bara,（,0,流量）,+ THRM.COEFF 1.33E-6,10,编写inprep文件：10,编写,inprep,文件：,输入管道参数：,举例：,PIPEPARMS,+ FRICTION COLE 0.03 /*,管道粗糙度按,0.03mm,+ KNOT 1 /*,距离计算步长,1km,+ INITIAL 40 /*,线路最高点起始压力,40bara,（,0,流量）,+ THRM.COEFF 1.33E-6,注意,1,：距离步长越长，计算速度越快，但瞬态计算精度下降，通常,KNOT,步长不超过管道长度的,2%,。,注意,2,：起始压力只影响平衡时间，不影响计算结果。,注意,3,：原油和成品油粗糙度的选取：直缝管和无缝管取,0.05mm,；螺旋缝管取,0.125mm,（,OD=DN350,）。天然气管道粗糙度选取：无内涂层取,0.03mm,，有内涂层取,0.01mm,。,11,编写inprep文件：11,输气管道摩阻计算公式和输气效率说明,.doc,COLEBROOK,公式实际上是,Prandtle,（普朗特）水力光滑管公式和,Nikuradse,（尼古拉兹）完全粗糙管公式的数学组合。该公式为隐函数形式，计算量较大，但计算精度高。它适用于整个紊流区（，并且做为管道水力计算的基本公式被世界许多国家采用。,GB50251,中附录,A,中的气体流量计算公式由,Pandandle,公式推导出，,仅适用于手工计算，精度较差，引入了输气效率,E,作为修正。若采用,COLEBROOK,公式则需要采用电脑计算，则不需要再考虑输气效率,E,。,可以在交互状态的,spans,中修改输气效率,E,（默认值为,1,），但不建议修改。,12,12,编写,inprep,文件：,限制参数范围语句：,格式：,SET.LIMIT,+ KEY SUB ATT LL LOW HIGH HH DEF UNIT,+ KEY SUB ATT LL LOW HIGH HH DEF UNIT,+ .,举例：,SET.LIMIT,+ T * OD * * 1422 1800 * MM,/*,低于,LL,或高于,HH,的值产生,error,；,/*,在,LL,和,LOW,之间的值产生,warning,；,/*,在,HIGH,和,HH,之间的值产生,warning,；,/*,默认值用*代替，也可以用具体参数修改默认值,/*,可以在,html,文件夹中打开,limits.html,查看所有参数的限制值和默,/*,认值以及单位。,13,编写inprep文件：13,编写,inprep,文件：,关闭流体批量跟踪和混合：,格式：,NOTRACK,注意,1,：此语句只能在,STATE,和,EQUIPMENT,语句之间输入,注意,2,：如果有此语句，就不允许在,EXTERNAL,处输入流体性质。,14,编写inprep文件：14,编写,inprep,文件：,传热模型选取：,ISOTHERMAL TEMP,绝热模型，温度为恒定值。通常成品油管道采用该模型。,THERMAL TEMP,部分传热模型，管道不换热，设备换热，可以计算泵、压缩机、调节阀的前后的温度变化。,TRANSTHERMAL TEMP,设备和管道均换热,+ SP.HEAT SPH,流体比热，,BWRS,和,SCL,模型不需要输入,+ HEATCOND HCOND,流体导热系数,KJ/HR-M-DC,天然气取,0.108,原油和成品油取,0.360.58,W/M-DC=KJ/HR-M-DC*0.278,HCOND(KJ/HR-M-DC)=0.493(1-0.00054T)/DENSITY(,相对密度）,15,编写inprep文件：15,16,16,编写,inprep,文件：,+ MIN.FILM.COEF MFC,最小膜系数，,KJ/HR-M2-DC,，模拟流体和管道内壁之间的液体薄膜传热系数。默认取,0,。,+ HEAT.FRIC.EXP HFELH HFELC HFETH HFETC,热阻修正系数，默认取,0,。,+ FORCED.CONVECTION FORCED,强制对流换热基准条件，取,WALL.TEMP or MEAN.TEMP,+ FREE.CONVECTION FREE,自然对流换热基准条件，,取,WALL.TEMP or MEAN.TEMP,本页中的参数通常不需要输入，取默认值即可。对于特殊情况（特殊介质、流态、温度等），需要检索传热学资料或相关专家确定。,17,编写inprep文件：17,编写,inprep,文件：,传热模型选取：,ISOTHERMAL TEMP,绝热模型，温度为恒定值。通常成品油管道采用该模型。,THERMAL TEMP,部分传热模型，管道不换热，设备换热，可以计算泵、压缩机、调节阀的前后的温度变化。,TRANSTHERMAL TEMP,设备和管道均换热,+ SP.HEAT SPH,流体比热，,BWRS,和,SCL,模型不需要输入,+ HEATCOND HCOND,流体导热系数,KJ/HR-M-DC,天然气取,0.108,原油和成品油取,0.360.58,W/M-DC=KJ/HR-M-DC*0.278,HCOND(KJ/HR-M-DC)=0.493(1-0.00054T)/DENSITY(,相对密度）,18,编写inprep文件：18,编写,inprep,文件：,流体状态方程选取：,STATE AGA TRIVB TRIVC /*,美国煤气协会计算方法,+ GAS,SG CO2 HHV /*,相对密度，,CO2,摩尔组份含量，高位发热值,注,1,：,适用温度范围为,0C to 55C,，最高到,8.3MPa,。,注,2,：在适用范围内，此方程的计算结果与,BWRS,结果接近。,注,3,：详细方程见输气管道相关教材。,注,4,：天然气的相对密度是指在相同压力和温度条件下，天然气和干空气的密度之比。,SG=M/28.964,，,M,为天然气的平均分子量，甲烷的分子量为,16,。标准状态下空气的密度为,1.293kg/m3,，标准状态下甲烷的密度为,0.717kg/m3,。,注,5,：适用于已知密度的气体。,19,编写inprep文件：19,编写,inprep,文件：,流体状态方程选取：,STATE CNGA SG HHV /*,美国加利福尼亚天然气协会计算方法,注,1,：,适用相对密度范围为,0.55 0.7,，最高到,6.9MPa,，温度,060,。,注,2,：采用,SPS,，该方程的计算结果与,BWRS,和,AGA,差别较大。,注,3,：详细方程见输气管道相关教材。,20,编写inprep文件：20,编写,inprep,文件：,流体状态方程选取：,STATE BWRS.MOLE TRIVC TRIVB,/*.MOLE,代表摩尔组份，否则为质量组份含量,+ NAMES NAME1 NAME2 NAME3 .,+ INITIAL FR1 FR2 FR3 .,+ VISC V01 V02 V03 ./*,不同组份的粘度，默认值为,0.01cp,（动力粘度）,+ VPMI VPMI1 VPMI2 VPMI3 ./*,不同组份的粘度压力系数,+ VTMI VTMI1 VTMI2 VTMI3 ./*,不同组份的粘度温度系数,举例：,STATE BWRS.MOLE,+ NAMES C1 C2 C3 IC4 NC4 IC5,+ CO2 N2 H2S,+ INITIAL 0.925469 0.039582 0.003353 0.001158 0.000863 0.00221,+ 0.018909 0.008455 0.000001,注,1,：,BWRS,方程的适用范围最宽，计算精度最优,注,2,：除天然气外，还可以计算,CO2,和乙烯管道。,注,3,：对于以甲烷为主的天然气，不需要输入,VPMI,和,VTMI,，但对于,CO2,等特殊介质管道需要检索有关资料确定上述,2,个参数。,21,编写inprep文件：21,编写,inprep,文件：,流体状态方程选取：,STATE SCLPROP TRIVC TRIVB ASTM,/*PROP,代表模拟混油，不输入则不进行混油计算。,ASTM,表示使用,ASTM,粘度公式,22,编写inprep文件：22,编写,inprep,文件：,流体状态方程选取：,+ FLUID FNAME,+ DENSITY P0 T0 p0 PM0 TM0 PTMULT PPMULT TTMULT,/*,标准压力（绝压） 标准温度 标准密度,/*,液体的密度受温度影响较大，对压力不敏感。,/*,柴油密度,820870kg/m3,，汽油,700770kg/m3,，煤油,750830kg/m3,/*,不同原油密度差别很大，,6501060kg/m3,+ VISC (,u,0 VPMI VTMI) | (A B),/*,牛顿流体的原油管道推荐输入,A B,值。,/*,特别注意：应反复试算,A B,值，确保在输送温度范围内计算粘温曲线与实,/*,际粘温曲线尽量靠近。,/*,成品油管道和非牛顿流体输入,u0,，标准状态下的粘度,/*,对于非牛顿流体，还要额外输入,VISCOSITY,曲线数据,23,编写inprep文件：23,编写,inprep,文件：,流体状态方程选取：,+ HCAP Cv0 CvT,/*,流体热容，原油为,2.02.1kJ/kgDC,+ HCOND K0 KT,/*,流体导热系数，通常取,0.14W/mDC,+ COLOR color,/*,油品颜色,blue, red, green,等,24,编写inprep文件：24,编写,inprep,文件：,非牛顿流体的粘度输入：,对于非牛顿流体，需按下面方法输入粘度曲线值：,VISCOSITY,+ FLUID name1,+ V(T,S) s1 s2 . sm,+ t1 v11 v12 v1m,+ t2 v21 v22 v2m,+ .,+ tn vn1 vn2 vnm,/*sm,代表剪切速率，,tn,代表温度，,vnm,代表对应粘度,25,编写inprep文件：25,举例：,VISCOSITY,+ FLUID CRUDE,+ V(T,S) 10 20 30 40 50,+,0,16 16 17 18 19,+,5,13 13 14 15 16,+,10,10 11 12 13 14,+,15,8 8 9 9 10,+,20,5 5 6 6 7,+,25,4 5 5 6 6,+,30,4 4 4 5 5,26,26,编写,inprep,文件：,选择量纲类型：,METRIC,米制或公制,或者,ENGLISH,英制或美制,27,编写inprep文件：27,编写,inprep,文件：,开始输入物理模型：,= EQUIPMENT,/*,阀门、管道、泵、压缩机等物理设备在该语句以下输入,/*,详见,help,文件,28,编写inprep文件：28,编写,intran,文件,29,编写intran文件29,Intran,文件必须编写的内容,BEGIN,(,必须首行,),INTERACTIVE,(if you want to run the simulation interactively),TRENDLIST,(if you want to produce interactive and/or GRAFR time plots or run TPORT),PROFILE (INTRAN),(if you want to produce GRAFR profile plots),SHARE,(if you want to run TPORT),Note: The INTRAN file must be entered in capital letters. For more syntax rules, see,Input syntax,.,30,Intran文件必须编写的内容30,编写,intran,文件：,开始命令,BEGIN RECNO, + BEGIN.TIME = BEGIN.TIME/*,开始时间，分钟, + END.TIME = END.TIME/*,结束时间，分钟, + PRESSURE.TOLERANCE = PRES.TOLER,/*,允许压力误差，默认,2psi=13.8kPa, + TEMP.TOLERANCE = TEMP.TOLER,/*,允许温度误差，默认,0.1,华氏度,=0.37,举例：,BEGIN 0,+ BEGIN.TIME = 360,+ END.TIME = 1440,BEGIN 0,+ BEGIN.TIME=0,+ END.TIME=1E20,31,编写intran文件：31,编写,intran,文件：,执行交互命令,INTERACTIVE CRT,ROWS=ROWS, COLS=COLS/*,屏幕行列数, FONT=font, /*,字体大小,fg=fgcolor, /*,字体颜色,bg=bgcolor, /*,字体背景色,ffg=ffgcolor, /*,数字颜色,fbg=fbgcolor/*,数字背景色,32,编写intran文件：32,编写,intran,文件：,存储,timeplot,图形和参数,TRENDLIST dev_names, + KEY.LETTER = dev_kl_include, + DEVICE.EXCLUDE = dev_kl_exclude, + PEEK.MATCH = peek_names_incl, + PEEK.EXCLUDE = peek_names_excl, + UNITS.MATCH = units_include, + SUB.TYPE = subtypes_include, + ECHO = echo,/*,该指令规定,SPS,存储哪些参数的,timeplot,/*TRENDLIST *,表示存储所有参数的,timeplot,/*,存储过多的参数将导致,review,文件体积庞大,33,编写intran文件：33,编写,intran,文件：,存储,distplot,图形和参数,PROFILE INT,/*,按一定时间间隔不断存储,distplot,图形和参数,/*,交互命令状态下，不需要该语句也可以产生,distplot,/*,如果没有,GRAF,文件，则可以不输入该语句,34,编写intran文件：34,编写,intran,文件：,SHARE,SHARE dev_names, + KEY.LETTER = dev_kl_include, + DEVICE.EXCLUDE = dev_kl_exclude, + PEEK.MATCH = peek_names_incl, + PEEK.EXCLUDE = peek_names_excl, + UNITS.MATCH = units_include, + SUB.TYPE = subtypes_include,/*,用于,TPORT,从,trans,中提取参数,/*,必须输入,TRENDLIST,语句，才能使用,SHARE,35,编写intran文件：35,编写,intran,文件：,规定开始图形界面,MACRO (INIT, OVERVIEW),/*,采用,OVERVIEW.DSP,为开始画面，,OVERVIEW.DSP,是用户自定义的文件,/* OVERVIEW.DSP,需要在,trans,内建立。,/*,其余语句详见,help,文件,36,编写intran文件：36,软件操作与练习,37,软件操作与练习37,习题一：,建立简单天然气管道模型,基本参数：,38,组,分,C1,C2,C3,iC4,nC4,iC5,Mol%,93,4,1.5,0.5,0.5,0.5,管道起点压力：,4.0MPa,管道末端压力：,2.0MPa,管道起点高程：,200m,管道终点高程：,400m,管径：,D323.911mm,管道沿线地温：,10,管道起点温度：,30,管道长度：,100km,管道埋深：,1200mm,土壤导热系数,4.6,KJ/hr-m-DC,问题,1,：管道输气能力,=,？（万方,/,天）,问题,2,：管道总传热系数,=,？,问题,3,：说明管道温度降低的原因,问题,4,：管道起终点的天然气实际密度,=,？,问题,5,：管道起终点的天然气流速,=,？,问题,6,：管道摩阻和高程差产生的压降,=,？,问题,7,：输出全线压力和流量曲线,问题,8,：输气全线温度温度和地温曲线,问题,9,：起点压力提高,1MPa,后输量,=,？,问题,10,：终点压力降低,1MPa,后的输量,=,？,习题一： 建立简单天然气管道模型38组 分C1C2C3iC,习题二：,建立简单原油管道模型,基本参数：,39,管道起点压力：,4.0MPa,管道末端压力：,2.0MPa,管道起点高程：,200m,管道终点高程：,400m,管径：,D323.911mm,管道沿线地温：,10,管道起点温度：,30,管道长度：,100km,管道埋深：,1200mm,土壤导热系数,4.6KJ/KG-DC,粘度：,A=16.49 B=6.59,标准密度：,845kg/m3,问题,1,：管道输油能力,=,？（吨,/,天）,问题,2,：管道总传热系数,=,？,问题,3,：说明管道温度降低的原因,问题,4,：管道起终点的原油粘度和密度,=,？,问题,5,：管道起终点的原油流速,=,？,问题,6,：管道摩阻和高程差产生的压降,=,？,问题,7,：输出全线压力和流量曲线,问题,8,：输气全线温度粘度和地温曲线,问题,9,：起点压力提高,1MPa,后输量,=,？,问题,10,：终点压力降低,1MPa,后的输量,=,？,习题二： 建立简单原油管道模型39管道起点压力：4.0MP,习题三：复杂原油管道模拟与调试,基本参数：,40,水平管道，常温输送,首站进站压力：,500kPa,末站进站压力：,500kPa,管径：,D40611mm,管道沿线地温：,10,管道起点温度：,30,管道长度：,200km,管道埋深：,1200mm,土壤导热系数,4.6KJ/KG-DC,粘度：,A=16.49 B=6.59,标准密度：,845kg/m3,管道设计压力：,6.3MPa,设计输量：,350,万吨,/,年,问题,1,：合理布置泵站，以完成设计输量,问题,2,：设计输量下各泵站进出口压力、功率,=,？,问题,3,：将输量分别调整到,350,和,450m3/hr,，此时各泵站水头、功率,=,？,问题,4,：管道最大输送能力,=,？,问题,5,：单座泵站的最大输送能力,=,？,问题,6,：输出全线介质粘度和温度曲线,问题,7,：判断管道所处的流态区域,问题,8,：在中间泵站设置进出站压力选择控制系统，保证进站压力不低于,1000kPa,，出站压力不高于,6300kPa,。,习题三：复杂原油管道模拟与调试40水平管道，常温输送问题1：,习题四：输油管道水击模拟与分析,基本参数：,41,同习题三,出站设置单向阀,出站,ESD,阀门后设置水击泄放阀，设定值：,Po=6300kPa,，,Pc=5000kPa,。,问题,1,：在,500m3/h,输量下，第二站进站,ESD,阀门误关。分析首站出站设置和不设置水击泄放阀的差别，并输出如下数据和曲线：,1,）首站出站压力、流量变化曲线。,3,）第二站进站压力、流量变化曲线。,4,）水击过程中全线各点达到的最高压力。,5,）水击过程中,习题四：输油管道水击模拟与分析41同习题三问题1：在500m,习题五： 输气管道小时峰调节分析,基本参数：,42,组,分,C1,C2,C3,iC4,nC4,iC5,Mol%,93,4,1.5,0.5,0.5,0.5,管道起点压力：,9.8MPa,管道末端压力：,1.6MPa,水平管道，,330km,管径：,D323.911mm,管道沿线地温：,20,管道起点温度：,20,用户用气规律：,16,点：,30,万方,/,小时,718,点：,40,万方,/,小时,1924,点：,30,万方,/,小时,问题,1,：找出调节小时峰能力最大的口径,问题,2,：输出管道末端和起点的流量变化曲线,问题,3,：分析影响管道小时峰调节能力的因素,习题五： 输气管道小时峰调节分析42组 分C1C2C3iC,习题六： 输气管道压气站布站学习,基本参数：,43,组,分,C1,C2,C3,iC4,nC4,iC5,Mol%,93,4,1.5,0.5,0.5,0.5,管道起点压力：,7.0MPa,管道末端压力：,4.0MPa,管道长度：,4000km,管径：,D1016181mm,管道沿线地温：,20,管道起点温度：,20,管道设计压力：,10MPa,管道年输量：,120,亿方,压气站驱动方案：奇数站为燃驱站，其余为电驱站,管道里程和高程：,距离,高程,水平距离,km,高程,m,0,800,1000,1000,2000,2000,3000,800,4000,100,问题：按,1.35,压比布站，计算全线压缩机总轴功率和能耗。若压比增加和减小，总功率如何变化。,习题六： 输气管道压气站布站学习43组 分C1C2C3iC,习题七： 顺序输送管道建模和调试练习,基本参数：,44,以习题三为基础,柴油批量：,20000m3,90,号和,93,号汽油各,10000m3,顺序：柴,9093,柴,在,2,号泵站设置分输阀,柴油分输量,1000m3,，分输流量,100m3/h,；,90,和,93,号汽油分输量均为,500m3,；分输流量,50m3/h,问题,1,：研究柴油顶汽油通过输油泵的压力和流量变化规律，并解释变化原因。,问题,2,：研究顺序输送过程中管道总摩阻和流量的变化规律，并解释变化原因。,习题七： 顺序输送管道建模和调试练习44以习题三为基础问题1,