坐标变换和动态数学数学模型的简化.ppt

上传人:wuy****n92 文档编号:77664367 上传时间:2023-03-16 格式:PPT 页数:49 大小:592KB
返回 下载 相关 举报
坐标变换和动态数学数学模型的简化.ppt_第1页
第1页 / 共49页
坐标变换和动态数学数学模型的简化.ppt_第2页
第2页 / 共49页
点击查看更多>>
资源描述

《坐标变换和动态数学数学模型的简化.ppt》由会员分享,可在线阅读,更多相关《坐标变换和动态数学数学模型的简化.ppt(49页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、第三章第三章 高动态性能变频调速系统高动态性能变频调速系统山东大学山东大学v问题的提出问题的提出v三相异步电动机的动态数学模型三相异步电动机的动态数学模型v坐标变换和动态数学模型的简化坐标变换和动态数学模型的简化v矢量控制的变频调速系统矢量控制的变频调速系统v直接转矩控制变频调速系统直接转矩控制变频调速系统v无速度传感器变频调速系统无速度传感器变频调速系统本章提要3.2 坐标变换和动态数学模型的简化 上节中虽已推导出异步电机的动态数学模上节中虽已推导出异步电机的动态数学模型,但是,要分析和求解这组非线性方程显然型,但是,要分析和求解这组非线性方程显然是十分困难的。在实际应用中必须设法予以简是十

2、分困难的。在实际应用中必须设法予以简化,简化的基本方法是化,简化的基本方法是坐标变换坐标变换。直流电动机调速性能优异且便于控制是因为其具备以下几个条件:直流电动机调速性能优异且便于控制是因为其具备以下几个条件:(1)直流电动机的主磁场由直流励磁电流产生,补偿绕组基本上克)直流电动机的主磁场由直流励磁电流产生,补偿绕组基本上克服了电枢反应,所以一般认为其服了电枢反应,所以一般认为其主磁场是一个稳定的直流磁场主磁场是一个稳定的直流磁场。(2)当电刷位于几何中性线上时,电枢磁场与主磁场在空间是垂直)当电刷位于几何中性线上时,电枢磁场与主磁场在空间是垂直的(互差的(互差90 电角度),电角度),是自然

3、解耦的是自然解耦的。(3)励磁电流和电枢电流)励磁电流和电枢电流互相独立互相独立,各自在不同的回路中,控制简,各自在不同的回路中,控制简单,易于实现。单,易于实现。(4)直流电动机的动态数学模型只有一个输入)直流电动机的动态数学模型只有一个输入/输出变量输出变量电枢电电枢电压压/转速,在工程允许的一些假定条件下,直流电动机可以描述成转速,在工程允许的一些假定条件下,直流电动机可以描述成单单输入单输出的二阶线性系统输入单输出的二阶线性系统。异步电动机和直流电动机相比有着本质上的区别:异步电动机和直流电动机相比有着本质上的区别:(1)三相异步电动机的定子通以三相平衡正弦交流电流,产生一个)三相异步

4、电动机的定子通以三相平衡正弦交流电流,产生一个随时间和空间都在变化的旋转磁场随时间和空间都在变化的旋转磁场。(2)转子电流也产生旋转磁场,它和定子旋转磁场相位不同,但稳)转子电流也产生旋转磁场,它和定子旋转磁场相位不同,但稳态时都是态时都是同步旋转的,在空间上不存在垂直关系同步旋转的,在空间上不存在垂直关系。(3)三相异步电动机(鼠笼式)的)三相异步电动机(鼠笼式)的转子是短路的转子是短路的,只能调节定子电只能调节定子电流。流。(4)异步电动机的数学模型至少是一个七阶的模型异步电动机的数学模型至少是一个七阶的模型,其输入量为电,其输入量为电压(电流)和频率,输出量为磁链和转速,是一个多变量系统

5、。在压(电流)和频率,输出量为磁链和转速,是一个多变量系统。在静止的静止的A、B、C坐标系中,异步电动机的数学模型为时变方程组。坐标系中,异步电动机的数学模型为时变方程组。如果能够简化异步电动机的动态数学模型,从如果能够简化异步电动机的动态数学模型,从而像直流电动机那样分别独立控制励磁电流和而像直流电动机那样分别独立控制励磁电流和转矩电流,并使它们的磁场在空间位置上也互转矩电流,并使它们的磁场在空间位置上也互差差90 电角度,就可以获得像直流电动机那样电角度,就可以获得像直流电动机那样优异的调速性能优异的调速性能。由以上分析可以推想:由以上分析可以推想:坐标变换坐标变换 交流电机的物理模型二、

6、坐标变换二、坐标变换坐标变换:从一种坐标轴系转换到另一种坐标轴系的变换坐标变换:从一种坐标轴系转换到另一种坐标轴系的变换矢量控制相关的三种坐标系:静止的三相矢量控制相关的三种坐标系:静止的三相ABC坐标系、静坐标系、静止的二相止的二相、坐标系和旋转的二相坐标系和旋转的二相d、q坐标系;坐标系;由机电能量转换的基本原理可知,电动机内气隙磁场是进行由机电能量转换的基本原理可知,电动机内气隙磁场是进行能量转换的媒介,由定子侧输入的能量正是通过气隙磁场传能量转换的媒介,由定子侧输入的能量正是通过气隙磁场传递到转子的。递到转子的。在进行坐标变换时,只要能使变换前后产生的气隙基波合成在进行坐标变换时,只要

7、能使变换前后产生的气隙基波合成磁势不变(幅值和空间相位相同),两者就是等效的。磁势不变(幅值和空间相位相同),两者就是等效的。因此,因此,磁势不变是不同坐标系间进行变换的一项基本原则。磁势不变是不同坐标系间进行变换的一项基本原则。(一)坐标变换的原则(一)坐标变换的原则设在某坐标系下的电路或系统的电压和电流向量分别为设在某坐标系下的电路或系统的电压和电流向量分别为u和和i,在新的坐标系下,电压和电流向量变成,在新的坐标系下,电压和电流向量变成u和和i,定义新向,定义新向量与原向量的坐标变换关系为量与原向量的坐标变换关系为 Cu u=u (3-33)Ci i=i (3-34)其中其中Cu和和Ci

8、分别为电压和电流变换阵。分别为电压和电流变换阵。当满足功率不变的约束条件时,当满足功率不变的约束条件时,Cu和和Ci的关系为的关系为 CuT Ci=I (3-35)式中式中 I为单位矩阵为单位矩阵。这里对这里对Cu和和Ci的选择并没有加任何约束,它们可以是任意的。的选择并没有加任何约束,它们可以是任意的。在一般情况下,为了使变换阵简单易记,令在一般情况下,为了使变换阵简单易记,令Cu=Ci=C即把电压和电流变换阵取为同一矩阵,则式(即把电压和电流变换阵取为同一矩阵,则式(3-35)变成)变成 CT C=I 或或 CT=C-1(3-36)式(式(3-36)就是)就是坐标变换满足功率不变的约束条件

9、坐标变换满足功率不变的约束条件,且取,且取电压和电流变换阵相同时对变换矩阵的要求,这样的坐标电压和电流变换阵相同时对变换矩阵的要求,这样的坐标变换属于正交变换。变换属于正交变换。众所周知,交流电机三相对称的静止绕组众所周知,交流电机三相对称的静止绕组 A、B、C,通,通以三相平衡的正弦电流时,所产生的合成磁动势是旋转磁以三相平衡的正弦电流时,所产生的合成磁动势是旋转磁动势动势F,它在空间呈正弦分布,以同步转速,它在空间呈正弦分布,以同步转速 1(即电流的(即电流的角频率)顺着角频率)顺着 A-B-C 的相序旋转。的相序旋转。ABCABCiAiBiCF1图图a 三相交流绕组三相交流绕组(二)(二

10、)3s/2s变换变换 旋转磁动势的产生旋转磁动势的产生然而,旋转磁动势并不一定非要三相不可,除单相以外,二然而,旋转磁动势并不一定非要三相不可,除单相以外,二相、三相、四相等任意对称的多相绕组,通以平衡的多相电相、三相、四相等任意对称的多相绕组,通以平衡的多相电流,都能产生旋转磁动势,当然以两相最为简单。流,都能产生旋转磁动势,当然以两相最为简单。不同电机模型彼此等效的原则是:不同电机模型彼此等效的原则是:在不同坐标下所产生的磁在不同坐标下所产生的磁动势完全一致。动势完全一致。根据电机学原理,异步电动机三相绕组的作用,完全可以根据电机学原理,异步电动机三相绕组的作用,完全可以用在空间上互相垂直

11、的两个静止的用在空间上互相垂直的两个静止的、绕组来代替,如图绕组来代替,如图3-6所示。由三相所示。由三相ABC轴系变换到两相轴系变换到两相轴系轴系以产生同样以产生同样的旋转磁势为准则,并需要满足功率不变的约束条件的旋转磁势为准则,并需要满足功率不变的约束条件。将交流电机的物理模型等效地变换成类似直流电机的模式,将交流电机的物理模型等效地变换成类似直流电机的模式,分析和控制就可以大大简化。分析和控制就可以大大简化。坐标变换正是按照这条思路坐标变换正是按照这条思路进行的。进行的。uCuAuCi600600AuuBBiAiiBiC图3-6 定子ABC轴系到静止轴系的变换(2)等效的两相交流电机绕组

12、)等效的两相交流电机绕组图B 两相交流绕组 两相静止绕组两相静止绕组 和和 ,它们在空间互差它们在空间互差90,通以时,通以时间上互差间上互差90的两相平衡交流的两相平衡交流电流,也产生旋转磁动势电流,也产生旋转磁动势 F。当两个旋转磁动势大小当两个旋转磁动势大小和转速都相等时,即认为图和转速都相等时,即认为图b的两相绕组与图的两相绕组与图a的三相绕组的三相绕组等效。等效。Fii1(3)旋转的直流绕组与等效直流电机模型)旋转的直流绕组与等效直流电机模型1FMTimitMT图c 旋转的直流绕组 再看图再看图c中的两个匝数相等且互相垂直的绕组中的两个匝数相等且互相垂直的绕组 M 和和 T,其中分别

13、通以直流电流,其中分别通以直流电流 im 和和it,产生合成磁动,产生合成磁动势势 F,其位置相对于绕组来说是固定的。,其位置相对于绕组来说是固定的。如果让包含两个绕组在内的整个铁心以同步转速如果让包含两个绕组在内的整个铁心以同步转速旋转,则磁动势旋转,则磁动势 F 自然也随之旋转起来,成为旋转磁自然也随之旋转起来,成为旋转磁动势。动势。把这个旋转磁动势的大小和转速也控制成与图把这个旋转磁动势的大小和转速也控制成与图 a 和图和图 b 中的磁动势一样,那么这套旋转的直流绕组中的磁动势一样,那么这套旋转的直流绕组也就和前面两套固定的交流绕组都等效了。当观察者也就和前面两套固定的交流绕组都等效了。

14、当观察者也站到铁心上和绕组一起旋转时,在他看来,也站到铁心上和绕组一起旋转时,在他看来,M 和和 T 是两个通以直流而相互垂直的静止绕组。是两个通以直流而相互垂直的静止绕组。如果控制磁通的位置在如果控制磁通的位置在 M 轴上,就和直流电轴上,就和直流电机物理模型没有本质上的区别了。这时,绕组机物理模型没有本质上的区别了。这时,绕组M相当相当于励磁绕组,于励磁绕组,T 相当于伪静止的电枢绕组。相当于伪静止的电枢绕组。等效的概念等效的概念 由此可见,以产生同样的旋转磁动势为准则,图由此可见,以产生同样的旋转磁动势为准则,图a的三相交流绕组、图的三相交流绕组、图b的两相交流绕组和图的两相交流绕组和图

15、c中整体旋转中整体旋转的直流绕组彼此等效。或者说,在三相坐标系下的的直流绕组彼此等效。或者说,在三相坐标系下的 iA、iB、iC,在两相坐标系下的,在两相坐标系下的 i、i 和在旋转两相坐标系和在旋转两相坐标系下的直流下的直流 im、it 是等效的,它们能产生相同的旋转磁动是等效的,它们能产生相同的旋转磁动势。势。有意思的是:就图有意思的是:就图c的的M、T两个绕组而言,当两个绕组而言,当观观察者站在地面看上去,察者站在地面看上去,它们是与三相交流绕组等效的它们是与三相交流绕组等效的旋转直流绕组;如果旋转直流绕组;如果跳到旋转着的铁心上看跳到旋转着的铁心上看,它们就,它们就的的确确是一个直流电

16、机模型了。这样,通过坐标系的的确确是一个直流电机模型了。这样,通过坐标系的变换,可以找到与交流三相绕组等效的直流电机模的变换,可以找到与交流三相绕组等效的直流电机模型。型。现在的问题是,如何求出现在的问题是,如何求出iA、iB、iC 与与 i、i 和和 im、it 之间准确的等效关系,这就是之间准确的等效关系,这就是坐标变换坐标变换的任务的任务。注意注意:在这里,不同电机模型彼此在这里,不同电机模型彼此等效的原则等效的原则是:是:在在不同坐标下所产生的磁动势完全一致。不同坐标下所产生的磁动势完全一致。2.三相三相-两相变换(两相变换(3/2变换)变换)现在先考虑上述的第一种坐标变换现在先考虑上

17、述的第一种坐标变换-在三相静在三相静止绕组止绕组A A、B B、C C和两相静止绕组和两相静止绕组、之间的变换,或之间的变换,或称三相静止坐标系和两相静止坐标系间的变换,称三相静止坐标系和两相静止坐标系间的变换,简称简称 3/2 3/2 变换变换。下图中绘出了下图中绘出了 A、B、C 和和 、两个坐标系,为两个坐标系,为方便起见,取方便起见,取 A 轴和轴和 轴重合。设三相绕组每相有效轴重合。设三相绕组每相有效匝数为匝数为N3,两相绕组每相有效匝数为,两相绕组每相有效匝数为N2,各相磁动势,各相磁动势为有效匝数与电流的乘积,其空间矢量均位于有关相为有效匝数与电流的乘积,其空间矢量均位于有关相的

18、坐标轴上。由于交流磁动势的大小随时间在变化着,的坐标轴上。由于交流磁动势的大小随时间在变化着,图中磁动势矢量的长度是随意的。图中磁动势矢量的长度是随意的。三相和两相坐标系与绕组磁动势的空间矢量三相和两相坐标系与绕组磁动势的空间矢量 AN2iN3iAN3iCN3iBN2i60o60oCB 设磁动势波形是正弦分布的,当三相总磁动势与二相设磁动势波形是正弦分布的,当三相总磁动势与二相总磁动势相等时,两套绕组瞬时磁动势在总磁动势相等时,两套绕组瞬时磁动势在 、轴上的投轴上的投影都应相等,影都应相等,写成矩阵形式,得写成矩阵形式,得(3-37)匝数比应为匝数比应为 当三相总磁动势与二相总磁动势相等时,并

19、考虑变换前当三相总磁动势与二相总磁动势相等时,并考虑变换前后总功率不变,得后总功率不变,得 代入式(代入式(3-373-37),得),得(3-37-1)令令 C3s/2s 表示从三相坐标系变换到两相坐标系的变换矩阵,表示从三相坐标系变换到两相坐标系的变换矩阵,则则(3-38)三相三相两相坐标系的变换矩阵两相坐标系的变换矩阵 如果三相绕组是如果三相绕组是Y形联结不带零线,则有形联结不带零线,则有 iA+iB+iC=0,或或 iC=iA iB。代入式(。代入式(3-38)并整理后得)并整理后得 按照所采用的条件,电流变换阵也就是电压变换阵,按照所采用的条件,电流变换阵也就是电压变换阵,同时还可证明

20、,它们也是磁链的变换阵。同时还可证明,它们也是磁链的变换阵。3.两相两相两相旋转变换(两相旋转变换(2s/2r变换)变换)从上图等效的交流电机绕组和直流电机绕组物理模从上图等效的交流电机绕组和直流电机绕组物理模型的图型的图 b 和图和图 c 中从两相静止坐标系到两相旋转坐标系中从两相静止坐标系到两相旋转坐标系 d、q 变换称作两相变换称作两相两相旋转变换,简称两相旋转变换,简称 2s/2r 变换变换 其其中中 s 表示静止,表示静止,r 表示旋转。表示旋转。把两个坐标系画在一起,即得下图。把两个坐标系画在一起,即得下图。dqi1 cosdiiqisiid1jsindi 图中,两相交流电流图中,

21、两相交流电流 i、i 和两个直流电流和两个直流电流 id、iq 产生同样的以同步转速产生同样的以同步转速 1旋转的合成磁动势旋转的合成磁动势 Fs 。由于。由于各绕组匝数都相等,可以消去磁动势中的匝数,直接用各绕组匝数都相等,可以消去磁动势中的匝数,直接用电流表示,例如电流表示,例如 Fs 可以直接标成可以直接标成 is。但必须注意,这。但必须注意,这里的电流都是空间矢量,而不是时间相量。里的电流都是空间矢量,而不是时间相量。M,T 轴和矢量轴和矢量 Fs(is)都以转速)都以转速 1旋转,分量旋转,分量 id、iq 的长短不变,相当于的长短不变,相当于d,q绕组的直流磁动势。绕组的直流磁动势

22、。但但 、轴是静止的,轴是静止的,轴与轴与 d轴的夹角轴的夹角 随时随时间而变化,因此间而变化,因此 is 在在 、轴上的分量的长短也随时轴上的分量的长短也随时间变化,相当于绕组交流磁动势的瞬时值。间变化,相当于绕组交流磁动势的瞬时值。2s/2r变换变换公式公式由图可见,由图可见,i、i 和和 id、iq 之间存在下列关系之间存在下列关系写成矩阵形式,得写成矩阵形式,得 是两相旋转坐标系变换到两相静止坐标系的变换阵。是两相旋转坐标系变换到两相静止坐标系的变换阵。式中式中 两相旋两相旋转转两相静止坐两相静止坐标标系的系的变换变换矩矩阵阵(3-41)对两边都左乘以变换阵的逆矩阵,即得对两边都左乘以

23、变换阵的逆矩阵,即得(3-40)则两相静止坐标系变换到两相旋转坐标系的变换阵是则两相静止坐标系变换到两相旋转坐标系的变换阵是 电压和磁链的旋转变换阵也与电流(磁动势)旋转电压和磁链的旋转变换阵也与电流(磁动势)旋转变换阵相同。变换阵相同。两相静止两相静止两相旋两相旋转转坐坐标标系的系的变换变换矩矩阵阵v电压和磁链的旋转变换与电流的旋转变换相同。电压和磁链的旋转变换与电流的旋转变换相同。vdq坐标系产生的气隙磁势同坐标系产生的气隙磁势同坐标系一样,也正是坐标系一样,也正是ABC坐标系中三相绕组产生的气隙旋转磁势。但与坐标系中三相绕组产生的气隙旋转磁势。但与坐标系相坐标系相比,比,dq坐标系产生该

24、旋转磁势的方法不同:坐标系产生该旋转磁势的方法不同:它是在同步旋转的它是在同步旋转的dq线圈中通入两直流量线圈中通入两直流量id和和iq,合成磁势,合成磁势F1相对相对dq轴系是静止的,依靠轴系是静止的,依靠dq轴系本身的同步旋转,使轴系本身的同步旋转,使F1成为同步旋成为同步旋转的圆形磁势转的圆形磁势。v正是通过正是通过坐标系到坐标系到dq坐标系的变换,最终将三相正弦交坐标系的变换,最终将三相正弦交流电流变换为两相直流量。流电流变换为两相直流量。令矢量is和M轴夹角为s,已知id、iq,求is和1,就是直角坐标/极坐标变换,简称K/P变换。显然,其变换式应为 4.直角坐标直角坐标/极坐标变换

25、极坐标变换(K/P)(3-42)变换过程变换过程ABC坐标系坐标系 坐标系坐标系dq坐标系坐标系3/2变换变换C2s/2r三、异步电动机在三、异步电动机在、静止静止坐标系上的数学模型坐标系上的数学模型 把把异异步步电电机机在在三三相相静静止止ABC坐坐标标系系上上的的数数学学模模型型变变换换到到两两相相坐坐标标系系上上,由由于于两两相相坐坐标标轴轴互互相相垂垂直直,两两相相绕绕组组之之间间没没有有磁磁的的耦耦合合,仅仅此此一一点点,就就会会使使数数学学模模型简单了许多。型简单了许多。图图3-9 用用两两相相静静止止坐坐标标系系表表示示的的异异步步机机等等效电路效电路RsLmLmLsLsRsLm

26、LrLmLrRrRr1.电压方程电压方程 式中,下标式中,下标s和和r分别表示定子和转子变量;下标分别表示定子和转子变量;下标 和和 分别表示分别表示 轴和轴和 轴变量轴变量.坐标系定子等效两相绕组的互感;坐标系定子等效两相绕组的互感;(3-43)2.磁链方程磁链方程 ABC三相坐标系的磁链方程经坐标变换简化为以下三相坐标系的磁链方程经坐标变换简化为以下坐标系磁链方程:坐标系磁链方程:在在两两相相坐坐标标系系中中,定定子子和和转转子子的的等等效效绕绕组组落落在在互互相相垂垂直直的的两两根根轴轴上上,它它们们之之间间没没有有耦耦合合关关系系,互互感感磁磁链链只只在在同同轴绕组之间存在,所以式中的

27、每个磁链分量只剩下两项。轴绕组之间存在,所以式中的每个磁链分量只剩下两项。(3-44)3.3.电磁转矩方程电磁转矩方程 以上电压方程、磁链方程和电磁转矩方程再加上式以上电压方程、磁链方程和电磁转矩方程再加上式(3-1)运动方程和式()运动方程和式(3-2)转角微分方程构成了)转角微分方程构成了静静止坐标系上的异步电动机数学模型。止坐标系上的异步电动机数学模型。这种在两相静止坐这种在两相静止坐标系上的数学模型又称作标系上的数学模型又称作Kron异步电机方程式或双轴异步电机方程式或双轴原型电机(原型电机(Two Axis Primitive Machine)基本方程式。)基本方程式。(3-45)4

28、.在两相同步旋转坐标系(dq坐标系)上的数学模型 两相同步旋两相同步旋转转dq坐标系的旋坐标系的旋转速度等于定子转速度等于定子电源的同步角速电源的同步角速度度 1。用。用dq坐标坐标系表示的异步电系表示的异步电动机等效电路如动机等效电路如图图3-10所示。所示。图图3-10 异异步步电电动动机机在在同同步步旋旋转转dq坐坐标标系系的的等效电路等效电路idridsudsuqsiqsiqr1LmLmLmLsLsLrLrRsRsrdquqr=0udr=0LmRrRr1.电压方程v dq坐标系相对于转子的旋转角速度为坐标系相对于转子的旋转角速度为 1-s,即,即转差角速度。转差角速度。式(式(3-46

29、)的电压方程右边系数矩阵的每)的电压方程右边系数矩阵的每一项都是非零的,这说明异步机在二相同步旋转坐标系一项都是非零的,这说明异步机在二相同步旋转坐标系下的数学模型仍是强耦合的。下的数学模型仍是强耦合的。(3-46)2.磁链方程3.3.电磁转矩方程电磁转矩方程 由由于于dq坐坐标标系系与与电电动动机机气气隙隙磁磁场场同同步步旋旋转转,彼彼此此之之间间无无相相对对运运动动,当当A、B、C坐坐标标系系中中的的变变量量为为正正弦弦函函数数时时,dq坐坐标标系系中中的的变变量量将将是是直直流流量量,已已经经非非常常接接近近直直流流电电动动机机了了。但但是是,直直流流电电动动机机的的电电枢枢回回路路和和励励磁磁回回路路是是解解耦耦的的,而而异异步步机机在在二二相相同同步步旋旋转转坐坐标标系系下下的的数数学学模模型型仍仍是是强耦合的。强耦合的。(3-47)(3-48)v1.异步异步电机的坐机的坐标变换是在哪三个坐是在哪三个坐标系下系下?各有各有何特点何特点?v2.异步异步电机矢量控制方程包括哪些机矢量控制方程包括哪些?是如何得到的是如何得到的?

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 教育专区 > 初中资料

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁