ImageVerifierCode 换一换
格式:DOCX , 页数:33 ,大小:254.64KB ,
资源ID:9800694      下载积分:3 金币
快捷下载
登录下载
邮箱/手机:
温馨提示:
快捷下载时,用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)。 如填写123,账号就是123,密码也是123。
特别说明:
请自助下载,系统不会自动发送文件的哦; 如果您已付费,想二次下载,请登录后访问:我的下载记录
支付方式: 支付宝    微信支付   
验证码:   换一换

加入VIP,免费下载
 

温馨提示:由于个人手机设置不同,如果发现不能下载,请复制以下地址【https://www.bdocx.com/down/9800694.html】到电脑端继续下载(重复下载不扣费)。

已注册用户请登录:
账号:
密码:
验证码:   换一换
  忘记密码?
三方登录: 微信登录   QQ登录  

下载须知

1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。
2: 试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。
3: 文件的所有权益归上传用户所有。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 本站仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

版权提示 | 免责声明

本文(Get格雅温度控制系统设计计算机控制技术课程设计.docx)为本站会员(b****7)主动上传,冰豆网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知冰豆网(发送邮件至service@bdocx.com或直接QQ联系客服),我们立即给予删除!

Get格雅温度控制系统设计计算机控制技术课程设计.docx

1、Get格雅温度控制系统设计计算机控制技术课程设计温度控制系统设计(计算机控制技术课程设计)学 号: 计算机控制技术课程设计题 目温度控制系统设计学 院专 业班 级姓 名指导教师2021年月日课程设计任务书学生姓名:_ 专业班级:_指导教师: 周申培 工作单位:_题 目: 温度控制系统设计初始条件:被控对象为电炉,采用热阻丝加热,利用大功率可控硅控制器控制热阻丝两端所加的电压大小,来改变流经热阻丝的电流,从而改变电炉炉内的温度。可控硅控制器输入为05伏时对应电炉温度0300,温度传感器测量值对应也为05伏,对象的特性为积分加惯性系统,惯性时间常数为T140秒。要求完成的主要任务: 包括课程设计工

2、作量及其技术要求,以及说明书撰写等具体要求1设计温度控制系统的计算机硬件系统,画出框图;2编写积分别离PID算法程序,从键盘接受Kp、Ti、Td、T及的值;3计算机仿真被控对象,编写仿真程序;4通过数据分析Td改变时对系统超调量的影响。5. 撰写设计说明书。课程设计说明书应包括:设计任务及要求;方案比较及认证;系统滤波原理、硬件原理,电路图,采用器件的功能说明;软件思想,流程,源程序;调试记录及结果分析;参考资料;附录:芯片资料,程序清单;总结。时间安排:2021年5月23日2021年6月3日摘要比例-积分-微分控制简称PID控制,是控制系统中应用最为广泛的一种控制规律。实际运行的经验和理论的

3、分析都说明,这种控制规律对许多工业过程进行控制时,都能得到满意的效果。利用计算机可以很好地使用PID算法对控制对象进行控制,具有较高的精度,并且可以很方便的改变PID参数,以到达不同的控制效果。本设计的控制对象为电炉,控制量为电炉温度,利用单片机对大功率可控硅导通角的控制,可以很方便地改变电热丝两端的电压,从而起到调节温度的作用。而热电偶配合单片机编程,能够较精确地得到炉温,使单片机能够实时发出控制信号,快速将炉温调节为给定值。当外界出现干扰使炉温发生变化时,单片机能够通过PID算法快速使炉温回到给定值。为了使PID控制更加稳定可靠,本设计参加了积分别离的改进措施,当偏差较大时取消积分作用,利

4、用PD控制快速使系统趋于稳定;当偏差小于某一个值时,参加积分作用,以消除静差。利用Matlab软件,可以通过仿真得到Td改变对系统超调量的影响。关键词:PID控制;Matlab;系统超调量1 设计任务及要求1.1 设计任务要求被控对象为电炉,采用热阻丝加热,利用大功率可控硅控制器控制热阻丝两端所加的电压大小,来改变流经热阻丝的电流,从而改变电炉炉内的温度。可控硅控制器输入为05伏时对应电炉温度0300,温度传感器测量值对应也为05伏,对象的特性为积分加惯性系统,惯性时间常数为T140秒。要求完成的主要任务:1设计温度控制系统的计算机硬件系统,画出框图;2编写积分别离PID算法程序,从键盘接受K

5、p、Ti、Td、T及的值;3计算机仿真被控对象,编写仿真程序;4通过数据分析Td改变时对系统超调量的影响。5. 撰写设计说明书。1.2 任务要求分析本系统的控制对象为电炉,被控量为温度,利用温度传感器实时检测电炉温度;将测得的数据经过A/D转换后送入计算机;计算机系统将检测得到的温度与炉温给定值进行比较,并计算偏差;按照预置的控制算法,对可控硅控制器的导通角进行调节,从而可以控制热阻丝两端的电压,起到温度调节的作用。为了实现数据采集、转换、处理以及PID算法控制,并通过键盘对温度设定值和PID控制参数进行修正,因此可以使用单片机或PLC。由于本次控制对象为电炉,其时间常数较大,因此采用周期不宜

6、过小,防止系统响应过于频繁,降低计算机系统的效率并使控制品质变坏;但也不能太大,否那么会使误差不能及时消除。2 方案比较及认证2.1 方案设计用温度传感器来检测炉的温度,将炉温转变成毫伏级的电压信号,经温度变送器放大并转换成电流信号。由电阻网络将电流信号变成电压信号,送入A/D转换器,通过采样和模数转换,所检测到的电压信号和炉温给定值的电压信号送入计算机程序中作比较,得出给定值与实际值之间的偏差,并与进行比较,从而确定算法。计算得到的控制量输出给可控硅控制器,改变可控硅的导通角,到达调压的目的,是电阻丝两端的电压增大或较小,进而实现对炉温的控制。方案一:使用PLC直接对铁块温度进行PID调节。

7、其根本思想为,触摸屏要设定温度给PLC,PLC控制光耦电路,然后对电阻丝进行控制加热。具体流程为,由触摸屏设定温度给PLC,由PLC通过电路控制铁块的加热或不加热,热电偶把温度以电压的方式经过AD转换反响给PLC,PLC再经过PID来控制温度。其结构框图如图2-1所示。图2-1 方案一设计结构图方案二:单片机热电偶温度自动控制。主要的控制芯片采用AT89C51,要求传感器测量的电压范围和可控硅控制器的电压在0-5,所以A/D与D/A转换芯片采用ADC0809和DAC0832。炉温控制在0300内,因此采用镍铬-铜镍热电偶,同时选用运算放大器将信号放大。图2-2 方案二设计结构图2.2 方案认证

8、根据设计要求,综合考虑选择方案二。在系统中,利用热电偶测得电阻炉实际温度并转换成毫伏级电压信号。该电压信号经过温度检测电路转换成与炉温相对应的数字信号进入单片机,单片机进行数据处理后,通过液晶显示器显示温度,同时将温度与设定温度比较,根据设定计算出控制量,根据控制量通过控制继电器的导通和关闭从而控制电阻丝的导通时间,以实现对炉温的控制。程序流程图如图2-3所示。图2-3 程序流程图3 系统软件设计3.1 PID控制算法模拟PID控制器的控制规律为3-1在PID调节中,比例控制能迅速反响误差,从而减小误差,但比例控制不能消除稳态误差,KP的加大,会引起系统的不稳定;积分控制的作用是:只要系统存在

9、误差,积分控制作用就不断地积累,输出控制量以消除误差,因而,只要有足够的时间,积分控制将能完全消除误差,积分作用太强就会使系统超调增大,甚至使系统出现振荡;微分控制可以减小超调量,克服振荡,使系统的稳定性提高,同时加快系统的动态响应速度,减小调节时间,从而改善系统的动态性能。为了便于计算机实现PID控制算法,必须将式3-1变换成差分方程,以得到数字PID位置型控制算式3-2根据式3-2可写出u(k-1)的表达式3-3将式3-2与式3-3相减,可以得到数字PID增量型控制算式为3-4式中,Kp为比例增益;Ki=Kp*T/Ti为积分系数;Kd=Kd*Td/T为微分系数。相对于位置型算法,增量型算法

10、不需要做累加,计算误差或计算精度对控制量的计算影响较小,而位置型算法要用到过去的累加值,容易产生较大的累加误差。位置型算法不仅要占用较多的内存单元,而且不便于编写程序,并且逐渐增大的累加误差可能引起系统冲击,严重影响系统的稳定性。综合考虑,应该使用增量型数字PID控制算法来增加系统的稳定性以及控制精度。3.2 积分别离的PID控制控制算法在一般的PID控制中,当有较大的扰动或大幅度改变给定值时,由于此时有较大的偏 差,以及系统有惯性和滞后,故在积分项的作用下,往往会产生较大的超调和长时间的波动。特别对于温度等变化缓慢的过程,这一现象更为严重,为此,可采用积分别离措施,即偏差e(k)较大时,取消

11、积分作用;当偏差较小时才将积分作用投入。亦即当时,采用PD控制;当时,采用PID控制。积分别离阈值应根据具体对象及控制要求。假设值过大时,那么达不到积分别离的目的;假设值过小,那么一旦被控量y(t)无法跳出个积分别离区,只进行PD控制,将会出现残差,为了实现积分别离,编写程序时必须从数字PID差分方程式中别离出积分项,进行特殊处理。积分别离PID控制算法流程图如图3-1所示。图3-1 积分别离PID控制算法流程图4 系统硬件设计4.1 系统滤波原理一般微机应用系统的模拟输入信号中,均含有种种噪音和干扰,它们来自信号源本身、传感器、外界干扰等。噪音有两大类:一类为周期性的,另一类为不规那么的。前

12、者可采用双积分A/D转换器,有效地消除其影响。后者为随机信号,可用数字滤波方法予以消除。 算术平均值法式要按输入的N个采样为周期ixi=1N,寻找这样一个y,使y与各采样值间的偏差的平方和为最小,使 由一元函数求值原理可得4.2 硬件设计原理该温度控制硬件设计采用了单片机AT89C52,A/D转换器ADC0809与D/A转换器DAC0832。其设计思想为:用热电偶来检测炉的温度,将炉温转变成毫伏级的电压信号,经温度变送器放大并转换成电流信号。由电阻网络讲电流信号变成电压信号,送入A/D转换器ADC0809,通过采样和模数转换,所检测到的电压信号和炉温给定值的电压信号都转换成数字量送入单片机AT

13、80C52进行比较,其差值即为实际炉温和给定炉温的偏差,以单片机为核心的数字PID控制器对偏差按照给定的方法运算,运算结果送入D/A转换器DAC0832转换成模拟电压,经功率放大器放大后送入晶闸管调压器,触发晶闸管并改变其导通角的大小,从而控制电阻炉的加温电压,起到炉温调节的作用。其电路图如图4-1所示。图4-1 硬件设计电路图5 系统仿真5.1 仿真程序及图形被控对象为采用simulink仿真,通过simulink模块实现积分别离PID控制算法。 选择适宜的Kp,Ki,Kd是系统的仿真效果趋于理想状态。MATLAB程序如下所示。clear all;close all;ts=2; %采样时间2

14、ssys=tf(1,40,1,0); %令sys为系统传递函数 dsys=c2d(sys,ts,zoh); %将sys离散化并加零阶保持器num,den=tfdata(dsys,v); %求sys多项式模型参数kp=1;ti=50;td=1;beta=0.1;ki=kp*ts/ti;kd=kp*td/ts;Simulink仿真图如图5-1所示。图5-1 Simulink仿真图5.2 仿真结果(1)当Td=0.1时,仿真波形图如图5-2所示。图5-2(2)当Td=1时,仿真波形图如图5-3所示。图5-3 当Td=1时的仿真波形图(3)当Td=10时,仿真波形图如图5-4所示。图5-4 当Td=1

15、0时的仿真波形图(4)当Td=15时,仿真波形图如图5-5所示。图5-5 当Td=15时的仿真波形图(5)当Td=20时,仿真波形图如图5-6所示。图5-6 当Td=20时的仿真波形图(6)当Td=30时,仿真波形图如图5-6所示。图5-7 当Td=30时的仿真波形图(7)当Td=40时,仿真波形图如图5-8所示。此时系统出现振荡,系统不稳定。图5-7 当Td=40时的仿真波形图5.3 结果分析增大微分时间Td,有利于加快系统响应,使超调量减小,稳定性增加,但系统对扰动的抑制能力减弱,对扰动有较敏感的响应。通过仿真可以看出,当K,Td,T,取了适宜的值后,改变Td既会改变系统的超调量,又会使系

16、统的稳定时间发生较大的改变。在一定范围内Td增大,而相应的系统超调量减小,说明此时微分作用会使系统超调量减小。但超过该范围,随Td增大超调量增大。仿真结果中,当Td=40时,系统甚至出现不稳定的情况。因此,积分时间常数Td必须取恰当值,微分作用过强,会使系统超调过大,甚至不稳定。因此,Td增大即微分作用的增强有助于增加系统的稳定性,还可以改善系统的动态特性,如可以明显减少超调量,缩短调节时间等,提高控制精度。但Td值偏大都会适得其反。6 心得体会通过本次温度控制系统设计,个人有了很多获得:一,了解了自身各种理论知识的缺乏之处,强化了自身知识水平;二,知道了理论必须与实践结合,各种技术才能真正到

17、达熟练运用的程度;三,科学技术是严密谨慎的,不能有一丝马虎,否那么可能会出现较大误差;四,通过本次温度控制设计,进一步学习了单片机设计与Matlab仿真设计。也学习了数字PID控制,比较了数字PID位置型控制算法与数字PID增量型控制算法。同时了解了数字PID控制器改进的方法,如积分别离。为了实现积分别离,编写程序时必须从数字PID差分方程式中别离出积分项,进行特殊处理。参考文献1 于海生. 计算机控制技术M. 北京:机械工业出版社,2021:101-111.2 张毅刚,彭喜元. 单片机原理及应用M. 北京:高等教育出版社,2021:43-82.3 刘教瑜. PLC应用技术M. 北京:人民邮电

18、出版社,2021:276-278.4 刘红丽. 传感与检测技术M. 北京:国防工业出版社,2021:73-81.5 王孝武,方敏,葛锁良. 自动控制理论M. 北京:机械工业出版社,2021:117-118.附录A 芯片资料附A1 ADC0809芯片功能A/D转换器是将模拟电压或电流转换成数字量的期间或装置,它是一种模拟系统和计算机之间的接口,它在数据采集和控制系统中,得到了广泛的应用,常用的A/D转换器有ADC0809. 它是一种带有8通道模拟开关的8位逐次逼近式A/D转换器,转换时间为100us左右,线性误差为1/2LSB,采用28脚双立直插式封装,ADC0809由8通道模拟开关、通道选择逻

19、辑、8位A/D转换器及三态输出锁存缓冲器组成。(1)8通道模拟开关及通道选择逻辑该局部的功能是实现8选1操作,通道选择信号C、B、A与所选通道的关系如下:表A-1 ADC0809模拟开关及通道选择逻辑 地址锁存允许信号ALE、正脉冲用于通道选择信号C、B、A的锁存。加至C、B、A上的通道选择信号在ALE的作用下送入通道选择逻辑后,通道i上的模拟输入被送至A/D转换器转换。(2)8位A/D转换器图A-1 ADC0809引脚图如图A-1所示。IN0IN7为模拟信号的8个输入通道。VREF,VREF为基准电压的正极和负极。 ADDA 、ADDB和ADDC为模拟信号输入通道的地址选择线。ALE为地址锁

20、存信号,由低电平到高电平正跳变时讲地址选择线的状态锁存,一选通相应的输入通道。 START为启动信号,正脉冲的上升沿使内部存放器清零,从下降沿开始进行A/D转换。 OEC为转换结束信号,在START信号之后变低,转换结束为高电平,用来申请中断。OE为输出允许信号,有效时将输出存放器中的数据放到数据总线上。2-82-1为数码输出端,2-8为最低有效位,2-1为最高有效位。附A2 DAC0832芯片功能D/A转换器的功能是把二进制数字量电信号转换为与其数值成正比的模拟量电信号。常用D/A转换器为DAC0832芯片。DAC0832工作在单缓冲存放器方式,即当CS信号来时,D0D7数据线送来的数据直通

21、进行D/A转换,当IOW变高时,那么此数据便被锁存在存放器中,因此D/A转换的输出也保持不变。DAC0832讲输入的数字量转换成差动的电流输出Iout1和Iout2,为了将其编程电压输出,须经过运算放大器,使其输出05VVref为-5V或010VVref为-10V,假设要形成负电压输出,那么Vref需接正的基准电压。图A-2 DAC0832引脚图如图A-2 所示,DAC0832是双列直插式8位D/A转换器。能完成数字量输入到模拟量(电流)输出的转换。图4为DAC0832的引脚图。其主要参数如下:分辨率为8位,转换时间为1s,满量程误差为1LSB,参考电压为(+10/span-10)V,供电电源

22、为(+5+15)V,逻辑电平输入与TTL兼容。在DAC0832中有两级锁存器,第一级锁存器称为输入存放器,它的允许锁存信号为ILE,第二级锁存器称为DAC存放器,它的锁存信号也称为通道控制信号XFER。当ILE为高电平,片选信号 CS和写信号 WR1为低电平时,输入存放器控制信号为1,这种情况下,输入存放器的输出随输入而变化。此后,当WR1由低电平变高时,控制信号成为低电平,此时,数据被锁存到输入存放器中,这样输入存放器的输出端不再随外部数据DB的变化而变化。对第二级锁存来说,传送控制信号XFER和写信号WR2同时为低电平时,二级锁存控制信号为高电平,8位的DAC存放器的输出随输入而变化,此后

23、,当WR2由低电平变高时,控制信号变为低电平,于是将输入存放器的信息锁存到DAC存放器中。其余各引脚的功能定义如下: DI7DI0:8位的数据输入端,DI7为最高位。IOUT1:模拟电流输出端1,当DAC存放器中数据全为1时,输出电流最大,当 DAC存放器中数据全为0时,输出电流为0。IOUT2:模拟电流输出端2,IOUT2与IOUT1的和为一个常数。RFB:反响电阻引出端,DAC0832内部已有反响电阻,所以 RFB端可以直接接到外部运算放大器的输出端,相当于将一个反响电阻接在运算放大器的输出端和输入端之间。 VREF:参考电压输入端,此端可接一个正电压,也可接一个负电压,它决定0至255的

24、数字量转化出来的模拟量电压值的幅度,VREF范围为+10。VREF端与D/A内部T形电阻网络相连。VCC:芯片供电电压,范围为5V15V。GND :模拟量地/数字量地,即模拟/数字电路接地端。附A3 AT89C51单片机AT89C51提供以下的功能标准:4K字节闪烁存储器,128字节随机存储器,32个I/O口,2个16位定时/计数器,1个5向量两级中断结构,1个串行通信口,片内震荡器和时钟电路。另外,AT89C51还可以进行0HZ的静态逻辑操作,并支持两种软件的节电模式。8051单片机的内部结构十分复杂,但封装之后,只有引脚是面向用户的,所以使用者需要熟悉各引脚的用途。常用的8051芯片是用双

25、列直插40脚封装。如图A-3所示。图A-3 AT89C51引脚图其各引脚功能如下所示:VCC:供电电压。 GND:接地。 P0口:P0口为一个8位漏级开路双向I/O口,每脚可吸收8TTL门电流。当P0口的管脚第一次写1时,被定义为高阻输入。P0能够用于外部程序数据存储器,它可以被定义为数据/地址的第八位。在FIASH编程时,P0口作为原码输入口,当FIASH进行校验时,P0输出原码,此时P0口外部必须被拉高。P1口:P1口是一个内部提供上拉电阻的8位双向I/O口,P1口缓冲器能接收输出4TTL门电流。P1口管脚写入1后,被内部上拉为高,可用作输入,P1口被外部下拉为低电平时,将输出电流,这是由

26、于内部上拉的缘故。在FLASH编程和校验时,P1口作为第八位地址接收。P2口:P2口为一个内部上拉电阻的8位双向I/O口,P2口缓冲器可接收,输出4个TTL门电流,当P2口被写“1时,其管脚被内部上拉电阻拉高,且作为输入。并作为输入时,P2口的管脚被外部拉低,将输出电流,这是由于内部上拉的缘故。P2口当用于外部程序存储器或16位地址外部数据存储器进行存取时,P2口输出地址的高八位。在给出地址“1时,利用内部上拉优势,当对外部八位地址数据存储器进行读写时,P2口输出其特殊功能存放器的内容。P2口在FLASH编程和校验时接收高八位地址信号和控制信号。P3口:P3口管脚是8个带内部上拉电阻的双向I/

27、O口,可接收输出4个TTL门电流。当P3口写入“1后,它们被内部上拉为高电平,并用作输入。P3口也可作为AT89C51的一些特殊功能口,P3口同时为闪烁编程和编程校验接收一些控制信号。RST:复位输入。当振荡器复位器件时,要保持RST脚两个机器周期的高电平时间。ALE/PROG:当访问外部存储器时,地址锁存允许的输出电平用于锁存地址的地位字节。在FLASH编程期间,此引脚用于输入编程脉冲。在平时,ALE端以不变的频率周期输出正脉冲信号,此频率为振荡器频率的1/6,可用作对外部输出的脉冲或用于定时目的。 PSEN:外部程序存储器的选通信号。在由外部程序存储器取指期间,每个机器周期两次PSEN有效

28、。但在访问外部数据存储器时,这两次有效的PSEN信号将不出现。 EA/VPP:当EA保持低电平时,那么在此期间外部程序存储器0000H-FFFFH,不管是否有内部程序存储器。注意加密方式1时,EA将内部锁定为RESET;当EA端保持高电平时,此间内部程序存储器。在FLASH编程期间,此引脚也用于施加12V编程电源VPP。XTAL1:反向振荡放大器的输入及内部时钟工作电路的输入。XTAL2:来自反向振荡器的输出。附录B 程序清单附B1 单片机程序代码单片机程序代码如下所示:积分别离PID控制算法子程序:START:MOV 68H,KP ;分别将KP ,TI ,TD, T,送入指定的存储单元MOV

29、 54H,TI MOV 55H,TDMOV 56H,TMOV 57H,MOV A,68H;计算KI=KP*T/TI MOV B,56HMUL ABMOV B,54HDIV ABMOV 69H,AMOV A,68H;计算KD=KP*TD/TMOV B,55HMUL ABMOV B,56HDIV ABMOV 6AH,ALOOP0: MOV DPTR #7FF0H;读取预定温度值,送入ADC0809的IN0口地址 MOV DPTR,A ;启动A/D转换LOOP1 JB P3.3,LOOP1;等待转换数据MOVX A,DPTR;读取ADC0809的IN0口转换后的数据MOV 5CH,A;将预定值数据放

30、入指定的存储单元MOV DPTR,#7FF1H;读取采样温度值,送ADC0809的IN1口地址MOVX DPTR,A;启动A/D转换LOOP2: JB P3.3,LOOP2;等待转换数据 MOVX A,DPTR;读取ADC0809的IN1口转换后的数据MOV 49H,A ;将采样值数据放入指定的存储单元MOV A,5CH ;计算e(i),先取温度给定值CLR CSUBB A,50H ;温度给定值-采样值JNC AA0;判断e(i)的正负,如果为正,跳至AA0CPL A;e(i)为负,下两条指令求补ADD A,#01H AA0MOV R0,57HSUBB A,57H;|e(i)|-JNC AA1;|e(i)|跳至AA1SJMP AA2;|e(i)|跳至AA2 AA1: LCALL PD;调用PD算法 AA2: LCALL PID;调用PID算法MOV A,7CH;将Ui通过DAC0832输出MOV DPTR,#7FF2HMOVX DPTR,AINC DPTRMOVX DPTR,ALCALL DELAY;调用延时子程序,等待下一次采样计算SJMP LOOP0;进入下一次控制计算DELAY PROC NEARDL0: MOV R6,#FFHDL1: MOV R7,#FFHDL2: MOV R5,#FFHDLS: DJNZ R5,DLSDJNZ R7,DL2DJNZ R6

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1