首页> 中国专利> 用于预测轮胎和道路表面之间的物理相互作用效应的方法

用于预测轮胎和道路表面之间的物理相互作用效应的方法

摘要

本发明涉及一种预测轮胎和道路表面(12)之间的至少一个物理相互作用效应的方法,包括:绘制所述道路表面(12)的实轮廓(16);将所述实轮廓(16)分割成多个层,每个层对应于所述实轮廓(16)的一个剖面;按照所述实轮廓(16)的剖面,从所述实轮廓(16)中确定出所述道路表面(12)的至少一个描述符的至少一个变化规律;对于每个变化规律来说,在剖面的所谓的相互作用值上确定了每个描述符的所谓的相互作用值;并且基于至少一个描述符的相互作用值来预测物理相互作用效应的值。

著录项

  • 公开/公告号CN102472690A

    专利类型发明专利

  • 公开/公告日2012-05-23

    原文格式PDF

  • 申请/专利号CN201080033120.1

  • 发明设计人 D·罗舒;

    申请日2010-07-19

  • 分类号G01M17/02(20060101);G01N33/42(20060101);

  • 代理机构11314 北京戈程知识产权代理有限公司;

  • 代理人程伟;王锦阳

  • 地址 法国克莱蒙-费朗

  • 入库时间 2023-12-18 05:21:27

法律信息

  • 法律状态公告日

    法律状态信息

    法律状态

  • 2017-10-27

    专利权的转移 IPC(主分类):G01M17/02 登记生效日:20171010 变更前: 变更后: 变更前:

    专利申请权、专利权的转移

  • 2015-07-29

    授权

    授权

  • 2012-12-05

    专利申请权的转移 IPC(主分类):G01M17/02 变更前: 变更后: 登记生效日:20121101 申请日:20100719

    专利申请权、专利权的转移

  • 2012-07-04

    实质审查的生效 IPC(主分类):G01M17/02 申请日:20100719

    实质审查的生效

  • 2012-05-23

    公开

    公开

说明书

技术领域

本发明涉及一种用于预测轮胎和道路表面之间的物理相互作用效 应的方法。它适用于对轮胎在路面上的滚动噪音进行预测,但是并不 限制于此。

背景技术

从现有技术中已知了对道路路面进行模拟的方法。在该方法中, 测量了路面的三维实轮廓。所述实轮廓包括由掺入基体(例如柏油) 的鹅卵石或砾石碎屑所形成的实压头。

为了简化所获得的实轮廓的使用,确定实压头的实几何描述符包 括,例如,压头的平均高度、压头的平均密度等。

每个实压头通过虚压头来进行模拟,例如基于实描述符的球体。 通过虚描述符来描述球体,例如每单位面积中的球体数量、每个球体 的半径、每个球体的位置等。因此,虚描述符可以获得虚轮廓,在该 虚轮廓中,具有大尺寸的实压头是通过具有大半径的球体来表示的, 而具有小尺寸的实压头是通过具有较小半径的球体来表示的。因此, 虚描述符在几何形状方面描述了虚压头,该虚压头对于其而言模拟了 实压头。因此,虚压头简化了实压头的几何描述符。

然而,因此而确定的虚描述符并不考虑路面和轮胎之间的相互作 用的物理实际情况来描述给定的物理效应,例如滚动噪音。

因此,有必要借助修正因子来修正虚描述符,其物理意义是忽略 了大部分时间。例如,计算了PSI(当前效用指数)描述符,或者DI (恶化指数)描述符以复杂的方式描述了路面的恶化。并且,这些修 正因子无法以量化且精确的方式预测所考虑的物理效应。

发明内容

本发明的目的是提供一种预测方法,其可以正确地描述轮胎和路 面之间的物理相互作用。

为此,本发明的主题是一种预测轮胎和路面之间的至少一个物理 相互作用效应的方法,其特征在于:

-绘制所述路面的实轮廓,

-将所述实轮廓分割成若干层,每个层对应于所述实轮廓的高度,

-基于所述实轮廓来确定出所述路面的至少一个描述符的至少一 个变化规律,所述变化规律是所述实轮廓的高度的函数,

-对于每个变化规律来说,在所述高度的所谓的相互作用值上确定 了每个描述符的所谓的相互作用值,并且

-基于至少一个描述符的相互作用值来预测所述物理相互作用效 应的值。

根据本发明的方法可以正确地描述轮胎和路面之间的物理相互作 用。通过确定变化规律并找到每个描述符的相互作用值,因此可以正 确地描述轮胎和路面之间的物理相互作用效应。因此,每个虚描述符 是路面的物理描述符而并不是每个实压头的简单的几何描述符,这与 现有技术中的方法并不相同。所述高度的相互作用值对应于轮廓高度, 在此高度下轮胎的压扁真实地体现了轮胎在路面上的实压扁。

根据所述方法的可选特征:

-所述描述符或多个描述符包括虚压头的至少一个所谓的虚描述 符,以用来模拟所述实压头;

-每个虚描述符的相互作用值一方面基于所述轮胎和所述路面之 间的相互作用规律而被确定,另一方面基于所述虚压头的每个变化规 律而被确定,并且

-基于至少一个虚描述符的相互作用值来预测所述物理相互作用 效应的值。

待被操作的数据量相对于相同的物理相互作用效应的数字处理操 作而言更少,后者可以导致过高的数据量并且导致计算时间太长。这 是因为对于参数的每个值而言,虚压头总是相同的,即使它们所模拟 的实压头全都不相同。并且对于参数的每个值来说,每个实压头通过 虚压头来进行模拟,该虚压头传递了轮胎和路面之间的相互作用的物 理实际情况,而不是像现有技术那样传递几何实际情况。

根据所述方法的可选特征:

-至少一个所谓的次要虚描述符的至少一个变化规律被确定为至 少一个所谓的主要虚描述符的函数。

-所述实压头的至少一个所谓的实描述符的至少一个变化规律被 确定为参数的函数,并且对于参数的每个值而言,每个主要虚描述符 的值是基于每个实描述符的每个变化规律并且基于每个次要虚描述符 的每个变化规律而被确定。

-所述实描述符和次要虚描述符分别描述了所述轮胎与实压头之 间和轮胎与虚压头之间的相互作用的类似的物理量。

有利地,所述实描述符包括:

ο通过所述实压头产生的轮胎的实变形,和

ο在轮胎中通过所述实压头而压入的实体积,

所述次要虚描述符包括:

ο通过所述虚压头产生的轮胎的虚变形,和

ο在轮胎中通过所述虚压头而压入的虚体积。

有利地,对于可变参数的每个值而言,通过相同的球体来表示虚 压头。

该球体虚压头可以真实地模拟实压头所形成的约束场,该实压头 在交通引起的刨平效应下一般呈圆形。此外,球体虚压头的虚描述符 很容易进行计算。这些描述符特别包括球体或透镜体的体积、球体半 径等。

根据本发明的方法的可选特征,基于作为至少一个描述符的相互 作用值的函数的用于预测滚动噪音的规律来预测轮胎在路面上的滚动 噪音。

由于每个实描述符的每个相互作用值正确地模拟了轮胎和路面之 间的相互作用,因此基于预测噪音的规律而能够精确地并且以量化的 方式确定出轮胎在路面上的滚动噪音。

这种预测规律对于路面设计特别有利。实际上,不同于在铺设路 面的道路上例如通过声音传感器实验测量滚动噪音,可以基于仅一米 或数米的道路样品来预测噪音。这提供了道路建设成本的节约和实验 成本的节约。

在给定条件下滚动的给定轮胎的滚动噪音特别取决于路面的粗糙 度。该粗糙度导致轮胎的振动,所述轮胎的振动产生噪音。路面也展 现出吸收和反射特性。滚动噪音随这些特性而变化。在已考虑针对物 理相互作用效应的物理方法下,根据本发明的方法可以分开这些特性 的每一个的贡献,且道路设计者可以确定导致有或多或少的噪音的路 面的一个或多个特性。相反,噪音测量是综合性的,不可能区分噪音 起源的物理机理。单独的测量不提供特性和滚动噪音之间的清楚关系, 且不允许控制路面设计。

有利地,所述描述符是路面的每单位面积中的实压头的实数量, 并且相互作用值是路面的每单位面积中的实压头的最大数量值。

已经令人惊讶的发现,可以借助每单位面积的实压头的最大数量 值来正确地预测轮胎的滚动噪音。此外,与现有技术的噪音预测方法 相对照而言,路面的每单位面积的压头数量很容易得到计算,在现有 技术的噪音预测方法中,必须不可避免地计算数十个描述符。

有利地,所述描述符是球体的虚半径并且相互作用值是球体的半 径的相互作用值。

已经发现,可以通过增加球体的虚半径的相互作用值来提高预测 噪音规律的精确度,该球体的虚半径的相互作用值是用来精确预测滚 动噪音的相关值。

本发明的主题还涉及一种计算机程序,其包括代码指令,当该计 算机程序在计算机上执行时,所述代码指令能够控制例如以上限定的 方法的步骤的执行。

本发明的主题还涉及一种用于记录数据的媒介,其包括记录形式 的诸如以上限定的程序。

本发明的主题还涉及使诸如以上限定的程序能够在远程通信网络 上获得,其目的是为了它的下载。

本发明的主题最后还涉及一种制造路面的方法,其包括例如以上 限定的方法的预测步骤。

附图说明

通过阅读如下说明书并同时参照附图可更好地理解本发明,所述 说明书仅以非限制性的实例给出,在附图中:

-图1表示道路路面的实轮廓和已过滤的路面的实轮廓;

-图2显示了将实轮廓划分为多个层的步骤;

-图3表示作为高度z的函数的实轮廓的实压头数目的变化规律;

-图4显示了通过实压头确定轮胎的实变形的步骤;

-图5显示了确定轮胎中的由实压头压入的实体积的步骤;

-图6显示了通过根据本发明的方法获得的虚轮廓;

-图7表示轮胎和路面之间的模拟的相互作用;

-图8显示了用于确定轮胎的虚变形和由轮胎中的虚压头压入的 虚体积所参照的球体的透镜部分;

-图9和10表示作为高度z的函数的虚描述符的变化规律。

在图中,已表示出对应于路面的常规纵向(X)、横向(Y)和竖 直(Z)取向的互相正交的轴X、Y、Z。

具体实施方式

已在图1和2中显示了旨在与轮胎相互作用的道路路面12。路面 12包括其中掺入不同尺寸的砾石碎屑的基体,所述基体包含柏油、沥 青和/或水泥。这些砾石碎屑(也称为粗砂)在路面12的表面上形成压 头14。压头14旨在与轮胎接触。

路面12根据下文所述的预测方法而制得。

优选利用激光轮廓仪绘制路面12的样品的二维的所谓实轮廓16。 在此情况中,实轮廓16的深度在长度L的路面上进行测量。通过获得 二维轮廓16,对于绘制三维轮廓,收集时间和待实施的硬件减少。实 轮廓16在平面X,Z中在图1中由虚线表示,并显示出波峰18和波谷 20。认为每个实轮廓16表示路面的厚度dy的切片。

作为变型,绘制数个二维轮廓16以构建三维样品。

轮廓16通过截止滤波器22进行过滤,所述截止滤波器22由图1 中的折线表示,并称为波峰线。波峰线22选择并连接显示高度z的波 峰18,使得在轮廓16上滚动的轮胎变得仅与所选的波峰22接触。最 后,通过将所选的波峰22的高度重新初始化至零高度而调整轮廓16。 然后获得经调整的轮廓24,其由图1中的实线表示。

在图2中已经显示了轮廓24。将轮廓24离散化为轮廓24的可变 参数的函数。在此情况中,可变参数为轮廓24的高度z。然后沿Z轴 将轮廓24划分为数个层26,例如一百个层,每个层对应于一个高度z。 所述层基本上为平面,并平行于X,Y平面。

定义了压头14的所谓的实描述符NR、γR和VR。下标R表明这些 描述符基于实轮廓24测得。也定义了高度z的剖面,其对应于垂直于 Z轴的水平面并在高度z延伸。对于每个实描述符NR、γR和VR,现在 将计算作为可变参数z的函数的变化规律。

参照图2和3,NR表示每单位面积的实压头14的数目。对于每个 高度z,计算由高度z的每个剖面所截取的每单位面积的轮廓24的压 头14的数目。然后变化规律NR(z)确定为可变参数z的函数。接着, 基于变化规律NR(z)确定所谓的相互作用值NR,E。相互作用值NR,E为路 面12的每单位面积的实压头14的数目的最大值。在图3中,对于zR,E=17 mm,每平方米为NR,E=115个压头。

在三维轮廓的情况中,对于每个高度z确定与轮胎接触的实压头 14的数目。对于路面的同一个切片dy而言,如果波峰18在所考虑的 切片dy中不共面,则与轮胎接触的实压头14的数目变化。

在二维轮廓的情况中,假设对于每个高度z,实压头14的数目的 最大值为实压头14的数目的实值。实际上,局部极大值通常表明在所 研究的切片dy中沿横向方向存在波峰18。

参照图4,γR表示由每个压头14所导致的轮胎的实局部变形。对 于每个高度z,计算由γz,i=cotan(αi)定义的每个变形γz,i,其中αi为在高 度z存在的每个压头14i的顶点的半角。接着,对于每个高度z,计算 变形γz,i的平均值。因此变化规律γR(z)确定为可变参数z的函数。

参照图5,VR表示在轮胎中由实压头14压入的每单位面积的实体 积。对于每个高度z,计算由定义的在高度z的由每个 压头14i压入的每个体积ΔVz,i,其中Lcz,i为在高度z的每个压头i的弦 长,hm,i为在长度Lcz,i上每个压头14i的平均高度。接着通过加入体积 ΔVz,i(z)来计算每个高度z的VR(z)。因此确定变化规律VR(z)为可变参 数z的函数,并由所定义,其中L表示轮廓24的长度。

路面12通过在图6中表示的虚轮廓28进行模拟。虚轮廓28包括 用于模拟实压头14的虚压头30。定义虚压头30的所谓虚描述符。每 个虚描述符根据变化规律而作为可变参数z的函数进行变化。

这些虚描述符包括主要虚描述符av、Rv、hv和次要虚描述符γV、 VV。下标V表明所述描述符为虚的。次要虚描述符γV、VV和主要描 述符γR、VR描述了轮胎分别与虚压头30和实压头14之间的相互作用 的类似物理量。以类似于实轮廓24的方式,定义虚轮廓28的高度z 的剖面,所述高度z的剖面对应于垂直于Z轴的水平面并在高度z延 伸。

虚压头30由球体32表示。作为变型,虚压头30可由适合于待模 拟的路面的类型和适合于物理相互作用的椎体、圆柱条或其他形状表 示。

轮胎和半径Rv的虚球体32之间的相互作用已在图7中表示。对 于施加于轮胎的负载Z,球体32产生深度为2hv的压痕,并包括与具 有最大半径av和高度hv的轮胎接触的透镜状部分34。轮胎和球体32 之间的相互作用遵循Hertz相互作用定律:

av=3(1-v2).Z.Rv4.E

其中v和E分别为轮胎橡胶的泊松比和刚度模量。

半径为Rv的球体32的透镜状部分34在图8中表示。对于每个次 级虚描述符γv、Vv,现在将计算每个变化规律作为主要虚描述符av、 Rv、hv的至少一个的函数,在本文作为av和Rv的函数。

γv表示由每个虚压头30(在本文由半球部分34)导致的轮胎的虚 变形。实变形基于在两个维度上的实压头14的轮廓16进行计算。虚 变形基于在三个维度上的球体32进行计算。因此,根据穿过每个球体 32的X,Z平面组(就像球体32沿虚轮廓28的随机设置),认为虚变 形对应于由每个球体32所导致的变形的平均值。因此,基于二维实轮 廓,所述方法有可能确定虚轮廓的三维描述符。

在每个剖面X,Z中,平均变形对应于由中间接触点施加的变形, 所述中间接触点即坐标为x=(av.cosθ)/2,y=(av.sinθ)/2的点。在中间接 触点的斜率由如下关系给出:

dzdx=-av2.cosθ.[Rv2-av24.(1-3.sin2θ]-12

为了获得在剖面X,Z组上的轮胎的虚变形γV,我们计算:

γm,v=2π0π/2-av2.cosθ.[Rv2-av24.(1-3.sin2θ]-12

因此,变化规律γv(z)确定为av(z)和Rv(z)的函数,因此作为z的函 数:

γm,v(z)=avπ(Rv2-av24)12+18π[av(Rv2-av24)12]3

VV表示由轮胎中的虚压头30(在本文由半球部分34)压入的虚体 积。

因此,变化规律VV(z)以如下方式确定为av(z)和Rv(z)的函数,因 此为z的函数:

VV(z)=π.Rv3.(23+cos3β3-cosβ)其中cosβ=Rv-hvav

接着,对于可变参数z的每个值,基于每个实描述符γR(z)、VR(z) 的每个变化规律和每个次要虚描述符γv(z)、Vv(z)的每个变化规律,确 定每个主要虚描述符av(z)、Rv(z)的值。在此情况中,对于可变参数z 的每个值,在一方面在实压入体积VR(z)和虚压入体积Vv(z)之间建立 等式,在另一方面在实变形γR(z)和虚变形γv(z)之间建立等式。因此每 个主要虚描述符av(z)、Rv(z)的变化规律确定为可变参数z的函数。然 后可谈及具有施加的变形和施加的体积的模型。

将会注意到,对于可变参数z的每个值,虚球体32是相同的,并 均由相同的主要描述符av(z)、Rv(z)进行描述。

具有施加的变形和施加的体积的主要虚描述符av(z)和Rv(z)的变化 规律分别在图9A和9B中表示。

然后,在一方面基于上文指出的轮胎和路面之间的相互作用的 Hertz定律,在另一方面基于主要虚描述符av(z)、Rv(z)的每个变化规律, 确定可变参数的所谓的相互作用值(本文为高度z的相互作用值zE)。 在示出的实例中,数字求解出由这三个规律形成的系统,获得zE=0.19 mm。对于每个变化规律,也确定在相互作用值zE下的每个描述符av(z)、 Rv(z)的相互作用值。获得av(zE)=1.85mm和Rv(zE)=5.14mm。

接着,预测在轮胎和路面之间的物理相互作用。在此情况中,基 于用于预测B=g(V1,...,Vn)类型的噪音B的规律,预测轮胎在路面12上的 滚动噪音B,其中V1,...,Vn为路面的描述符D1,...,Dn的值。

在第一实施方案中,n=2,D1为路面12的压头的特征数目,且D2为路面12的压头的特征尺寸。在示出的实例中,压头的特征数目为每 单位面积路面12的实压头12的数目的相互作用值NR,E,压头的特征 尺寸为虚球体32的半径的相互作用值Rv(zE)。

在该实施方案中,基于如下定义的滚动噪音B预测规律来预测噪 音B:

B=20log[(RRo)15(NoN)35(CCo)k(0.751-v2EEo)15]+Bo(等式1)

其中C为轮胎的滚动速度,v为橡胶的泊松比,且E为橡胶的刚度模 量,Bo为在具有参照虚描述符Ro和No的参照路面上以参照速度Co 滚动的具有刚度模量Eo的参照轮胎的参照滚动噪音。k属于区间 [1.2-2]。优选地,k=1.6。作为变型,k=1.2。在另一变型中,k=2。Bo 取决于进行测量的位置:在带有轮胎的车辆内部,在铺设路面的轨道 边缘等。

在第二实施方案中,n、D1和D2与第一实施方案相同。基于如下 定义的滚动噪音B预测规律来预测噪音B:

B=20log[(Rv(zE)Ro)15(NoNR,E)35(CCo)k]+Bo,(等式2)

当使用参照轮胎比较数个不同的路面时,橡胶的刚度模量E和泊 松比v是已知的,并包含于参照常数Bo′中。k属于区间[1.2-2]。优选 地,k=1.6。作为变型,k=1.2。在另一变型中,k=2。

在第三实施方案中,在已考虑对于某些类型的路面描述符Rv相对 于其他描述符具有可忽略的影响下,基于如下定义的滚动噪音B预测 规律来预测给定轮胎的滚动噪音B:

B=20log[(NoNR,E)35(CCo)k]+Bo,,(等式3)

在该实施方案中,n=1,路面压头的特征数目D1为仅有的所用的 描述符。k属于区间[1.2-2]。优选地,k=1.6。作为变型,k=1.2。在另 一变型中,k=2。通常观察到在通过等式2和3所预测的噪音B之间的 小差异。

在第四实施方案中,基于如下定义的滚动噪音B预测规律来预测 任意轮胎的滚动噪音B:

B=20log[(NoNR,E)35(CCo)k(0.751-v2EEo)15]+Bo,,,(等式4)

在该实施方案中,n=1,路面压头的特征数目D1为仅有的所用的 描述符。k属于区间[1.2-2]。优选地,k=1.6。作为变型,k=1.2。在另 一变型中,k=2。

在第五实施方案中,基于如下定义的滚动噪音B预测规律来预测 给定轮胎的滚动噪音B:

B=20log[(Rv(zE)Ro)15(NoNR,E)35]+Bo,,,,(等式5)

仅考虑路面的特性的该关系可以获得良好的预测。实际上,所述 规律考虑Rv(zE),尽管其显示比NR,E更低的数学权数,但其可以间接描 述用于制造路面的方法的特性,如岩石性质、破碎模式的影响。该关 系也可以了解用于制造路面的方法的特性的更改如何更改了滚动噪 音。

预测给定雕刻式轮胎的滚动噪音B的实例。

利用在各种路面T1-T4上的雕刻式轮胎进行试验。已根据本发明 的方法确定每个这些路面的描述符NR,E和Rv(zE)。利用激光轮廓仪获得 每个路面T1-T4的1米样品的实轮廓,所述激光轮廓仪的空间分辨率 沿X轴为dx=0.1mm,沿Z轴为dz=0.01mm。使用第二实施方案的 预测规律(等式2),其中,值Bo′=45dB。该Bo′值对应于对客运车辆 轮胎在车辆外部通常测得的试验平均值。

给定车辆在每个路面T1-T4上以60km/h行驶。然后在三个不同位 置测量噪音:接近轮胎、在车辆内部和当车辆经过时在车辆外部(也 称为滑行噪音)。结果在下文表1中给出。

表1

对于接近轮胎测得的噪音,获得如下关系B测量=0.73×B预测+52, 相关系数R2=0.8。

对于在车辆外部轮胎测得的噪音,获得如下关系B测量=1.26×B预测- 21.40,相关系数R2=0.96。

对于在车辆内部轮胎测得的噪音,获得如下关系B测量=1.23×B预测- 22.10,相关系数R2=0.97。

注意到,可通过根据测量位置校正Bo′的值而更准确地接近实验值 B测量

无论测量位置是哪里,预测模型和进行的测量之间的关系是非常 令人满意的,并有可能基于一个或多个路面的样品比较并预测由给定 轮胎的滚动所产生的噪音B。

预测给定无花纹轮胎的滚动噪音B的实例。

为了表征道路的粗糙度而不与轮胎的刻纹相互作用,使用无花纹 轮胎。在显示道路网的路面的多样性的各种路面R1-R6上进行试验。 已根据本发明的方法确定每个这些路面的描述符NR,E和Rv(zE)。利用激 光轮廓仪获得给定路面的1.7m样品的实轮廓,所述激光轮廓仪的空间 分辨率沿X轴为dx=0.4mm,沿Z轴为dz=0.01mm。使用第二和第 三实施方案的预测规律(等式2和3),值Bo′Bo″=45dB。

给定车辆在每个路面R1-R6上以60km/h行驶。进行接近轮胎的 数个噪音实验测量,并平均每个路面R1-R6的实验测量值。结果在下 文表2中给出。

表2

在一方面,注意到路面组R1-R6的滚动噪音B的预测是优良的, 由此确定所选择描述符的相关性以及用于非常不同的路面的滚动噪音 B预测规律的适合性。

此外,注意到,利用根据第三实施方案的规律所预测的噪音B的 值大于利用根据第二实施方案的规律所预测的噪音B的值。在该实施 例的架构内,相比于根据第二实施方案的规律,根据第三实施方案的 规律可以获得更接近于实验值的噪音B的预测值。

最后,注意到预测的噪音B的精确度大于之前实施例的预测的噪 音B的精确度。

根据本发明的方法的全部或部分可通过代码指令实施,当所述方 法的步骤在计算机上执行时,所述代码指令能够控制所述方法的步骤 的执行。所述指令可以从用于记录数据的介质(例如硬盘或闪存、CD 或DVD类型)上所记录的计算机程序发出。可以为这种程序做好准备 使其可以利用,其目的是为了在远程通信网络(如英特网或无线网) 上进行下载。因此,可以经由该网络将程序的升级发送至与所述网络 连接的计算机。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号