发明名称 一种基于模糊滑模的智能汽车横向控制方法
摘要 本发明公开了一种基于模糊滑模的智能汽车横向运动控制方法,包括以下步骤:智能汽车横向滑模控制设计和模糊滑模控制设计。首先,将横向偏差y<sub>L</sub>和方位偏差ε<sub>L</sub>融合为集成偏差e<sub>L</sub>,用e<sub>L</sub>设计滑模切换函数s。其次,选择s及其微分<img file="DDA00003539375100011.GIF" wi="29" he="47" />作为模糊滑模控制的输入变量,前轮转角δ<sub>fc</sub>作为输出变量。设计模糊滑模控制方法。本发明将模糊控制与滑模控制相结合,不同于常规模糊控制方法,本发明模糊控制的输入变量是滑模切换函数及其微分,这种改变既克服了滑模控制所容易产生的抖振现象,又增加了模糊控制不依靠系统精准模型的优点,可有效地克服汽车的非线性特性、参数不确定性以及外界的干扰,提高了控制器的精度、可靠性和稳定性。
申请公布号 CN103439884B 申请公布日期 2015.12.23
申请号 CN201310306007.1 申请日期 2013.07.19
申请人 大连理工大学 发明人 李琳辉;连静;王宏旭;孙永朝;郭烈;王蒙蒙;丁新立;李明
分类号 G05B13/04(2006.01)I;B60W10/20(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 大连东方专利代理有限责任公司 21212 代理人 李洪福
主权项 一种基于模糊滑模的智能汽车横向控制方法,其特征在于:包括以下步骤:A、智能汽车横向滑模控制设计A1、以横向偏差y<sub>L</sub>和方位偏差ε<sub>L</sub>作为控制模型的状态变量,以前轮转角δ<sub>fc</sub>作为控制模型的输入变量,建立含不确定因素和干扰的智能汽车横向运动控制模型;具体模型如式(1)所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>q</mi><mrow><mo>.</mo><mo>.</mo></mrow></mover><mo>=</mo><mrow><mo>(</mo><mi>A</mi><mo>+</mo><mi>&Delta;A</mi><mo>)</mo></mrow><mover><mi>q</mi><mo>.</mo></mover><mo>+</mo><mrow><mo>(</mo><mi>E</mi><mo>+</mo><mi>&Delta;E</mi><mo>)</mo></mrow><mi>q</mi><mo>+</mo><mrow><mo>(</mo><mi>B</mi><mo>+</mo><mi>&Delta;B</mi><mo>)</mo></mrow><mi>u</mi><mo>+</mo><mi>d</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000801224250000011.GIF" wi="1488" he="84" /></maths>式中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>A</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>A</mi><mn>4</mn></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000801224250000012.GIF" wi="288" he="164" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>E</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>E</mi><mn>3</mn></msub></mtd><mtd><msub><mi>E</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>E</mi><mn>6</mn></msub></mtd><mtd><msub><mi>E</mi><mn>5</mn></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000801224250000013.GIF" wi="329" he="159" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>B</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>B</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>B</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000801224250000014.GIF" wi="219" he="159" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>d</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>d</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>d</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000801224250000015.GIF" wi="194" he="159" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>A</mi><mn>1</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>v</mi><mi>x</mi></msub><mo>-</mo><msub><mi>a</mi><mn>11</mn></msub><mi>L</mi><mo>+</mo><msub><mi>a</mi><mn>12</mn></msub><mo>-</mo><msub><mi>a</mi><mn>21</mn></msub><msup><mi>L</mi><mn>2</mn></msup><mo>+</mo><msub><mi>a</mi><mn>22</mn></msub><mi>L</mi><mo>)</mo></mrow><mo>;</mo><msub><mi>E</mi><mn>2</mn></msub><mo>=</mo><mrow><mo>(</mo><mo>-</mo><msub><mi>a</mi><mn>11</mn></msub><msub><mi>v</mi><mi>x</mi></msub><mo>-</mo><msub><mi>a</mi><mn>21</mn></msub><msub><mi>v</mi><mi>x</mi></msub><mi>L</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000801224250000016.GIF" wi="1167" he="80" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>E</mi><mn>3</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>a</mi><mn>11</mn></msub><mo>+</mo><msub><mi>a</mi><mn>21</mn></msub><mi>L</mi><mo>)</mo></mrow><mo>;</mo><msub><mi>A</mi><mn>4</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>a</mi><mn>21</mn></msub><mi>L</mi><mo>-</mo><msub><mi>a</mi><mn>22</mn></msub><mo>)</mo></mrow><mo>;</mo><msub><mi>E</mi><mn>5</mn></msub><mo>=</mo><msub><mi>a</mi><mn>21</mn></msub><msub><mi>v</mi><mi>x</mi></msub><mo>;</mo><msub><mi>E</mi><mn>6</mn></msub><mo>=</mo><mo>-</mo><msub><mi>a</mi><mn>21</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000801224250000017.GIF" wi="1131" he="74" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>B</mi><mn>1</mn></msub><mo>=</mo><mrow><mo>(</mo><mo>-</mo><msub><mi>b</mi><mn>1</mn></msub><mo>-</mo><msub><mi>b</mi><mn>2</mn></msub><mi>L</mi><mo>)</mo></mrow><mo>;</mo><msub><mi>B</mi><mn>2</mn></msub><mo>=</mo><msub><mi>b</mi><mn>2</mn></msub><mo>;</mo><msub><mi>d</mi><mn>1</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>a</mi><mn>11</mn></msub><mi>L</mi><mo>-</mo><msub><mi>a</mi><mn>12</mn></msub><mo>+</mo><msub><mi>a</mi><mn>21</mn></msub><msup><mi>L</mi><mn>2</mn></msup><mo>-</mo><msub><mi>a</mi><mn>22</mn></msub><mi>L</mi><mo>)</mo></mrow><msub><mi>v</mi><mi>x</mi></msub><mi>&rho;</mi><mo>;</mo></mrow>]]></math><img file="FDA0000801224250000018.GIF" wi="1197" he="75" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>d</mi><mn>2</mn></msub><mo>=</mo><mrow><mo>(</mo><mo>-</mo><msub><mi>a</mi><mn>21</mn></msub><mi>L</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>v</mi><mi>x</mi></msub><mover><mi>&rho;</mi><mo>.</mo></mover><mo>;</mo></mrow>]]></math><img file="FDA0000801224250000019.GIF" wi="415" he="68" /></maths>其中,<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>a</mi><mn>11</mn></msub><mo>=</mo><mo>-</mo><mrow><mo>(</mo><msub><mi>C</mi><mi>r</mi></msub><mo>+</mo><msub><mi>C</mi><mi>f</mi></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>mv</mi><mi>x</mi></msub><mo>;</mo><msub><mi>a</mi><mn>12</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>L</mi><mi>r</mi></msub><msub><mi>C</mi><mi>r</mi></msub><mo>-</mo><msub><mi>C</mi><mi>f</mi></msub><msub><mi>L</mi><mi>f</mi></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>mv</mi><mi>x</mi></msub><mo>-</mo><msub><mi>v</mi><mi>x</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA00008012242500000110.GIF" wi="987" he="78" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msub><mi>a</mi><mn>21</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>L</mi><mi>r</mi></msub><msub><mi>C</mi><mi>r</mi></msub><mo>-</mo><msub><mi>L</mi><mi>f</mi></msub><msub><mi>C</mi><mi>f</mi></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>I</mi><mi>z</mi></msub><msub><mi>v</mi><mi>x</mi></msub><mo>;</mo><msub><mi>a</mi><mn>22</mn></msub><mo>=</mo><mo>-</mo><mrow><mo>(</mo><msubsup><mi>L</mi><mi>f</mi><mn>2</mn></msubsup><msub><mi>C</mi><mi>f</mi></msub><mo>+</mo><msubsup><mi>L</mi><mi>r</mi><mn>2</mn></msubsup><msub><mi>C</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>I</mi><mi>z</mi></msub><msub><mi>v</mi><mi>x</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA00008012242500000111.GIF" wi="1010" he="79" /></maths><maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msub><mi>b</mi><mn>1</mn></msub><mo>=</mo><msub><mi>C</mi><mi>f</mi></msub><mo>/</mo><mi>m</mi><mo>;</mo><msub><mi>b</mi><mn>2</mn></msub><mo>=</mo><msub><mi>L</mi><mi>f</mi></msub><msub><mi>C</mi><mi>f</mi></msub><mo>/</mo><msub><mi>I</mi><mi>z</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA00008012242500000112.GIF" wi="523" he="73" /></maths><maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><mi>q</mi><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>y</mi><mi>L</mi></msub></mtd><mtd><msub><mi>&epsiv;</mi><mi>L</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>]]></math><img file="FDA00008012242500000113.GIF" wi="302" he="91" /></maths>为控制模型的状态变量,<img file="FDA00008012242500000114.GIF" wi="44" he="84" />为控制模型状态变量的一阶微分,<img file="FDA00008012242500000115.GIF" wi="49" he="89" />为控制模型状态变量的二阶微分,u为控制模型的输入变量即前轮转角δ<sub>fc</sub>,L<sub>f</sub>、L<sub>r</sub>分别为质心到前后轴的距离(m),I<sub>z</sub>为智能汽车转动惯量(kg·m<sup>2</sup>),C<sub>f</sub>、C<sub>r</sub>分别为前后轮胎的侧偏刚度(N/rad),m为智能汽车质量(kg),L为智能汽车的预瞄距离,v<sub>x</sub>为智能汽车横向速度(m/s),ρ为路径曲率,d为干扰项,ΔA、ΔB、ΔE为不确定项;A2、计算横向偏差和方位偏差;经图像处理算法得到智能汽车的参考路径曲线后,在参考路径曲线与过预瞄点且平行于图像横向中心线的直线相交处作切线,获得智能汽车运动控制所需的路径信息,即预瞄点处智能汽车与参考路径的横向偏差和方位偏差信息;具体计算方法如下:定义图像的中心点位置O代表智能汽车的纵向中心线和横向中心线的交点,O<sub>1</sub>为预瞄点在图像中的位置,直线O<sub>1</sub>X<sub>1</sub>为过预瞄点O<sub>1</sub>且平行于图像横向中心线OX的直线;A21、横向偏差的计算定义参考路径曲线和直线O<sub>1</sub>X<sub>1</sub>的交点P(x<sub>0</sub>,y<sub>0</sub>)处的横坐标与图像预瞄点O<sub>1</sub>的横坐标之差为智能汽车在图像中的横向偏差:<maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><msub><mi>y</mi><mi>L</mi></msub><mo>=</mo><msub><mi>x</mi><mn>0</mn></msub><mo>-</mo><mfrac><msub><mi>w</mi><mn>1</mn></msub><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000801224250000021.GIF" wi="1239" he="130" /></maths>其中,w<sub>1</sub>为图像的宽度,横向偏差y<sub>L</sub>为像素值,需转换成实际横向偏差:<maths num="0015" id="cmaths0015"><math><![CDATA[<mrow><msub><mi>y</mi><mi>L</mi></msub><mo>=</mo><mi>&gamma;</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>-</mo><mfrac><msub><mi>w</mi><mn>1</mn></msub><mn>2</mn></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000801224250000022.GIF" wi="1207" he="141" /></maths>式中,γ为像素点与实际距离比例系数;A22、方位偏差的计算得到参考路径曲线之后,在其与直线O<sub>1</sub>X<sub>1</sub>的交点处P(x<sub>0</sub>,y<sub>0</sub>)作切线,设切线方程为:<maths num="0016" id="cmaths0016"><math><![CDATA[<mrow><mi>y</mi><mo>=</mo><msub><mi>K</mi><mi>d</mi></msub><mi>x</mi><mo>+</mo><msub><mi>B</mi><mi>d</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000801224250000023.GIF" wi="1217" he="99" /></maths>其中,K<sub>d</sub>为切线斜率,B<sub>d</sub>为切线截距,智能汽车在图像中的方位偏差定义为:<maths num="0017" id="cmaths0017"><math><![CDATA[<mrow><msub><mi>&epsiv;</mi><mi>L</mi></msub><mo>=</mo><mi>arctan</mi><mrow><mo>(</mo><msub><mi>K</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000801224250000024.GIF" wi="1245" he="89" /></maths>A3、针对式(1)所示的横向控制模型,将横向偏差和方位偏差按照无量纲化加权规则进行融合,使用融合后的集成偏差来设计滑模控制的切换函数s;基于数据处理的需要,将横向偏差和方位偏差分别无量纲化:<maths num="0018" id="cmaths0018"><math><![CDATA[<mrow><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>L</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><mn>2</mn><mfrac><mrow><msub><mi>y</mi><mrow><mi>L</mi><mi>max</mi></mrow></msub><mo>-</mo><msub><mi>y</mi><mi>L</mi></msub></mrow><mrow><msub><mi>y</mi><mrow><mi>L</mi><mi>max</mi></mrow></msub><mo>-</mo><msub><mi>y</mi><mrow><mi>L</mi><mi>min</mi></mrow></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000801224250000025.GIF" wi="1229" he="137" /></maths><maths num="0019" id="cmaths0019"><math><![CDATA[<mrow><msub><mover><mi>&epsiv;</mi><mo>&OverBar;</mo></mover><mi>L</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><mn>2</mn><mfrac><mrow><msub><mi>&epsiv;</mi><mrow><mi>L</mi><mi>max</mi></mrow></msub><mo>-</mo><msub><mi>&epsiv;</mi><mi>L</mi></msub></mrow><mrow><msub><mi>&epsiv;</mi><mrow><mi>L</mi><mi>max</mi></mrow></msub><mo>-</mo><msub><mi>&epsiv;</mi><mrow><mi>L</mi><mi>min</mi></mrow></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000801224250000026.GIF" wi="1215" he="147" /></maths>y<sub>Lmax</sub>和y<sub>Lmin</sub>分别为横向偏差的最大值和最小值,ε<sub>Lmax</sub>和ε<sub>Lmin</sub>分别为方位偏差的最大值和最小值;对式(6)和(7)处理后的横向偏差和方位偏差进行融合时加入权重系数,因此集成偏差表达为:<maths num="0020" id="cmaths0020"><math><![CDATA[<mrow><msub><mi>e</mi><mi>L</mi></msub><mo>=</mo><msub><mi>&eta;</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>L</mi></msub><mo>+</mo><msub><mi>&eta;</mi><mn>2</mn></msub><mo>&CenterDot;</mo><msub><mover><mi>&epsiv;</mi><mo>&OverBar;</mo></mover><mi>L</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000801224250000027.GIF" wi="1177" he="91" /></maths>其中η<sub>1</sub>、η<sub>2</sub>为权重系数,均大于0且η<sub>1</sub>+η<sub>2</sub>=1,e<sub>L</sub>为集成偏差;将集成偏差e<sub>L</sub>作为滑模切换函数的参数,设计滑模切换函数s,其表达式为:<maths num="0021" id="cmaths0021"><math><![CDATA[<mrow><mi>s</mi><mo>=</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mi>e</mi><mi>L</mi></msub><mo>+</mo><msub><mover><mi>e</mi><mo>.</mo></mover><mi>L</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000801224250000031.GIF" wi="1139" he="86" /></maths>c<sub>1</sub>为常数,<img file="FDA0000801224250000032.GIF" wi="60" he="81" />为集成偏差变化率;B、模糊滑模控制设计B1、选择切换函数s及其微分<img file="FDA0000801224250000033.GIF" wi="39" he="69" />作为模糊滑模控制的输入变量,前轮转角δ<sub>fc</sub>作为模糊滑模控制的输出变量;输入变量和输出变量对应模糊子集语言变量均设置为{NB,NM,NS,Z,PS,PM,PB},其中NB、NM、NS、Z、PS、PM、PB分别称为“负大”、“负中”、“负小”、“零”、“正小”、“正中”、“正大”;输入变量s、<img file="FDA00008012242500000310.GIF" wi="39" he="65" />和输出变量δ<sub>fc</sub>的模糊子集隶属度函数相同,对NM、NS、Z、PS和PM这五个变量级采用三角形函数,对NB采用Z型函数,PB采用S型函数;B2、采用专家经验法确定模糊滑模的控制规则,控制规则如表1所示:表1 控制规则表<img file="FDA0000801224250000034.GIF" wi="1986" he="1221" />每一条模糊控制规则由下面形式的“IF‑THEN”模糊语句构成:R<sup>(i)</sup>:如果s是<img file="FDA0000801224250000035.GIF" wi="80" he="88" />且<img file="FDA0000801224250000036.GIF" wi="44" he="75" />是<img file="FDA0000801224250000037.GIF" wi="88" he="88" />则δ<sub>fc</sub>是B<sup>i</sup>;其中,<img file="FDA0000801224250000038.GIF" wi="71" he="82" />和<img file="FDA0000801224250000039.GIF" wi="90" he="87" />为输入变量模糊子集的语言变量,B<sup>i</sup>是输出变量模糊子集的语言变量;i=1,2,...,49代表模糊控制规则的数目;举例说明,其中一条模糊控制规则可表示为:如果s是PB且<img file="FDA0000801224250000045.GIF" wi="44" he="70" />是PB则δ<sub>fc</sub>是NB该控制规则理解为s的趋势是增加的,为了减小偏差,需给车辆前轮一个负大转角来纠正车辆偏差;模糊推理采用曼达尼的max‑min合成法,采用重心法进行解模糊运算,确定出输出变量δ<sub>fc</sub>:<maths num="0022" id="cmaths0022"><math><![CDATA[<mrow><msub><mi>&delta;</mi><mi>fc</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>49</mn></munderover><msub><mi>&mu;</mi><msub><mi>C</mi><mi>i</mi></msub></msub><mo>&times;</mo><msub><mi>&delta;</mi><msub><mi>f</mi><mi>i</mi></msub></msub></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>49</mn></munderover><msub><mi>&mu;</mi><msub><mi>C</mi><mi>i</mi></msub></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000801224250000041.GIF" wi="1291" he="341" /></maths>其中<img file="FDA0000801224250000042.GIF" wi="101" he="92" />为输出变量<img file="FDA0000801224250000043.GIF" wi="78" he="102" />处的隶属度,即输出变量的隶属度函数在<img file="FDA0000801224250000044.GIF" wi="78" he="107" />处的取值。
地址 116024 辽宁省大连市高新园区凌工路2号