一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法

未命名 07-19 阅读:239 评论:0


1.本发明涉及可信性评估领域,更具体的说是涉及一种飞机结构数字孪生的可信性评估方法。


背景技术:

2.由于飞机运行环境复杂多变,且飞机不同结构层次之间关系繁杂,基于数字孪生的结构预测、优化决策过程中存在各种不确定性,这些不确定性会在飞机结构数字孪生系统中传递、迭代,最终导致仿真和预测结果的不确定性。现有的数字孪生预测过程中很少考虑不确定性因素对预测结果的影响,而不确定性因素深刻地影响着数字孪生系统的预测结果。因此为了提高数字孪生预测的精度与准确性,需要对数字孪生系统进行可信性评估,可信性评估可以识别数字孪生系统存在的潜在问题和不确定性,从而提高数字孪生的精度和准确性。
3.数字孪生可信性评估是提高数字孪生技术应用和发展的关键工具。现有技术对数字孪生技术的评估方法,有根据数字孪生系统的性能特点的8项模型评价准则和对应的29个量化指标,最终通过加权的方式获得综合评估指标的评估方法,但该方法指标繁杂,且容易引入评价人员的主观性。还有针对航空航天设备数字孪生系统的可信性评估的10维成熟度模型。尽管对数字孪生的评估方法已有多种,但是针对飞机结构数字孪生系统进行的可信性评估的方法还比较欠缺,因此,如何提供一种能够有效开展飞机结构数字孪生可信性评估的方法是本领域技术人员亟需解决的问题。


技术实现要素:

4.本发明的目的在于提供一种针对飞机典型结构数字孪生的可信性评估方法。可客观高效地对飞机结构数字孪生系统的可信性展开评估,进而降低基于数字孪生的飞机结构安全性设计和健康维护的决策风险,建立数字孪生系统的仿真可信性。
5.本发明的技术方案如下:步骤1,对要评估飞机结构数字孪生系统开展层次化解析,建立该数字孪生系统的参数关联模型;步骤2,分析所述参数关联模型的可信性关键要素,定义可信性关键要素成熟度等级划分标准,基于该标准得到该参数关联模型每个单模型的可信性关键要素的成熟度等级;步骤3,根据步骤2得到的可信性关键要素的成熟度等级,构建单模型的预测能力成熟度决策模型,依据可信性关键要素的权重系数,确定单模型的预测能力成熟度等级;步骤4,基于预测能力成熟度决策模型与贝叶斯网络,建立由可信性关键要素-模型预测能力成熟度-数字孪生系统可信性传递的贝叶斯网络,最终得到该飞机结构数字孪生系统的可信等级,可信等级满足要求则完成评估,不满足要求则需要对飞机结构数字孪生系统进行修正与更新。
6.作为优选,参数关联模型为该飞机结构数字孪生系统各个模型之间的联系而构建的参数模型,包括该数字孪生系统的每个单模型及其之间的关联。
7.作为优选,步骤2具体为:定义每个单模型的可信性关键要素,并将其划分为k个成熟度等级,定义每个可信性关键要素的成熟度等级划分标准;分析要评估飞机结构的每个单模型的可信性关键要素,获得对应的成熟度等级。
8.作为优选,步骤3中,基于下式计算单模型的预测能力成熟度等级:,式中,表示模型的预测能力成熟度等级,范围在0~k-1之间,表示该模型下每种关键要素的成熟度等级,表示权重系数,范围在0~1之间,n为每个单模型的关键要素数量。
9.作为优选,步骤4中,基于单模型的预测能力成熟度等级以及对模型进行灵敏度分析获取的灵敏度指标来确定贝叶斯网络内部的条件传递概率,层层传递,最终得到该数字孪生系统可信等级。
10.作为优选,灵敏度分析为:根据参数关联模型分别研究每个单模型的不确定性对数字孪生系统响应的敏感性和对可靠度的灵敏度,辨识出对数字孪生系统响应或可靠度影响较大的单模型。
11.作为优选,通过下式确定贝叶斯网络的传递概率,可得到该飞机结构数字孪生系统的可信等级:,式中,为可信等级,范围与的范围一致,在0~k-1之间,m为单模型的数量,下标i表示第i个单模型,表示第i个单模型的预测能力成熟度等级,表示通过灵敏度分析得到的第i个单模型的sobol灵敏度指标,取值范围在0~1之间,取值越大表示该模型发生变化或受不确定性影响时对数字孪生系统输出响应程度越大,取值越小则对输出响应程度越弱。
12.作为优选,建立的可信性关键要素-模型预测能力成熟度-数字孪生系统可信性传递的贝叶斯网络为三层网络结构,第一层为要素层,包含每个单模型的可信性关键要素的成熟度等级,第二层为模型层,包含每个单模型的预测能力成熟度等级,第三层为系统层,包含该数字孪生系统的可信等级;要素层到模型层的传递是基于预测能力成熟度决策模型,利用可信性关键要素的成熟度等级以及权重系数确定条件传递概率;模型层到系统层基于灵敏度分析确定条件传递概率层层递进。
13.作为优选,若步骤4得到的可信等级不满足用户要求,则需要对飞机结构数字孪生系统的模型进行修正或更新,如按照真实的机翼结构对模型的几何细节进行丰富或对模型进行不确定性量化和灵敏度分析等工作。
14.本发明的有益效果为:
1、本发明提供的评估方法实现了对飞机结构数字孪生系统可信性的量化,克服了现有数字孪生系统评估方法中指标繁杂、主观因素较强且没有考虑多源不确定性的不足。本发明将飞机结构数字孪生系统中存在的多源不确定性考虑在内,以尽量客观的方法验证飞机结构数字孪生模型与物理实体的符合性,减少了主观因素的影响,且对不同的数字孪生对象可以根据其特性调整可信性关键要素及贝叶斯网络结构,具有较好的灵活性,能够适用于一般飞机结构数字孪生系统的可信性评估。
15.2、本发明采用全新的数字孪生系统可信性评估方法,利用预测能力成熟度决策模型与贝叶斯网络,通过对该飞机结构数字孪生系统下每个单模型的可信性关键要素进行分析,得到每个单模型下可信性关键要素的成熟度等级,建立由可信性关键要素-模型预测能力成熟度-数字孪生系统可信性传递的贝叶斯网络,构建了一套飞机结构数字孪生可信性评估框架。
附图说明
16.图1为本发明以飞机柔性翼结构数字孪生系统为例的流程图;图2为本发明一个实施例的参数关联模型;图3为本发明一个实施例的贝叶斯网络传递示意图;图4为本发明一个实施例可信性评估的示意图。
具体实施方式
17.下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
18.如图1所示,本发明采用一种基于预测能力成熟度决策模型与贝叶斯网络的方法来评估飞机结构数字孪生系统的可信性,通过建立由可信性关键要素-模型预测能力成熟度-数字孪生系统可信性传递的贝叶斯网络,得到一套飞机结构数字孪生系统可信性评估框架及综合评估结果。
19.以飞机柔性翼结构数字孪生系统为例:步骤1:对飞机柔性翼结构数字孪生系统的模型结构开展层次化解析,如图2所示,飞机柔性翼结构数字孪生系统可认为由气动力学仿真模型和结构有限元仿真模型2个单模型组成。
20.步骤2:对空气动力学仿真模型和结构有限元仿真模型2个单模型的几何模型保真度(rgf)、物理及材料模型保真度(pmmf)、验证与确认(v&v)、灵敏度分析(ssa)、不确定性量化(uq)、数据适用性(da)6项可信性关键要素进行分析。
21.在飞机柔性翼结构数字孪生系统中,几何模型保真度(rgf)指仿真模型的详细表征水平或几何特征规范性;物理及材料模型保真度(pmff)主要指模型被校准的程度或模型基于物理的程度;验证与确认(v&v)分为模型验证和模型确认,模型验证是分析计算模型能否准确地代表基本数学模型及其求解过程,模型确认是确定一个模型能够准确描述真实物理世界的程度;灵敏度分析(ssa)是研究模型输入的不确定性对其输出影响的过程;不确定
性量化(uq)是识别、描述并对影响仿真模型精度的因素进行定量的过程;数据适用性(da)指获取的信息及数据适用于模型更新或进化的程度。
22.对几何模型保真度(rgf)、物理及材料模型保真度(pmmf)、验证与确认(v&v)、灵敏度分析(ssa)、不确定性量化(uq)、数据适用性(da)6项可信性关键要素划分为0级、1级、2级和3级4种成熟度等级,划分标准如下:几何模型保真度(rgf)的成熟度等级按以下描述划分:0级:模型的建立依赖判断和经验或者直接从模型库导出,没有针对真实的飞机结构进行细节化处理,很少或没有对模型的几何保证真度进行量化。
23.1级:对模型的主要部件进行了几何保真度量化,但是对真实飞机结构的细节处理仍然不够细致,具体的,对于飞机柔性翼结构,如铆接处的细节处理、表面光滑度、部件与部件之间的间隙等。
24.2级:对模型的绝大部分部件进行几何保真度量化,并且部分几何保真度量化工作由外部独立审查小组完成,对飞机结构的细节进行了真实化处理,具体的,对于飞机柔性翼结构,对铆钉连接处的细节处理,以及装配、磨损和疲劳引起的部件缺陷等都考虑在内。
25.3级:对模型的所有部件都进行了几何保真度量化,并且全部由外部独立审查小组完成。具体的,对于飞机柔性翼结构,完全依照真实柔性翼结构进行建模,没有对细节处进行近似处理或简化。
26.物理及材料模型保真度(pmff)成熟度等级按以下描述划分:0级:模型的建立基于经验,或者已知信息并不完整,模型中不同功能的原件大多为松耦合或者没有耦合,没有对模型的重要特征及参数进行校准。
27.1级:模型的建立并不完全基于经验,部分建模基于物理原理,模型整体的耦合性较差,对于模型的重要特征及参数进行了简易的校准。
28.2级:所有重要的物理过程模型和材料模型都是基于物理的。使用分离效应测试(set)和积分效应测试(iet)对模型的重要参数进行校准。重要物理过程之间耦合性良好,并且已经对物理和材料模型进行了一定程度的评审,例如非正式评审或内部评审。
29.3级:所有物理模型和材料模型都是基于物理的,使用分离效应测试(set)和积分效应测试(iet)进行校准的需求很小。所有的物理过程之间为紧耦合且是双向耦合。对物理和材料模型进行了独立的同行评审,例如:由建模与仿真工作的专家进行评审或外部专业审查小组进行的正式评审。
30.验证与确认(v&v)成熟度等级按以下描述划分:0级:模型的准确性评估几乎完全基于判断和经验。
31.1级:利用有限的计算结果和实验结果对模型的准确性进行了评估,进行评估的系统响应量(srq)与模型领域不直接相关或实验条件与所感兴趣应用不直接相关,系统响应量(srq)或实验条件表征中的实验不确定性在很大程度上是未确定的或基于经验的。
32.2级:对分离效应测试(set)实验和积分效应测试(iet)实验的一些关键系统响应量(srq)的计算结果和实验结果进行定量分析,并且对系统响应量(srq)和实验条件中存在的不确定性进行了表征。对模型验证和模型确认进行了评审,例如非正式评审或内部评审。
33.3级:对模型验证和模型确认进行专业评审,由外部专业审查小组完成。
34.灵敏度分析(ssa)成熟度等级按以下描述划分:
0级:仅依靠判断和经验对模型的灵敏度进行了评估,没有使用适于该模型的灵敏度评估指标。
35.1级:对仿真模型中的灵敏度分析是非正式或较为片面的,而不是系统的,使用了适用于该模型的灵敏度评估指标。
36.2级:对模型中大多数参数进行了灵敏度分析,使用了适用于该模型的灵敏度评估指标,对灵敏度分析进行了评审,如非正式的评审或内部评审。
37.3级:对参数不确定性和模型形式不确定性进行了全面且系统的灵敏度分析,使用高效且准确的灵敏度指标进行分析,并且数值逼近或采样误差对模型的灵敏度影响较小,对灵敏度分析进行了专业评审,由外部专业审查小组完成。
38.不确定性量化(uq)成熟度等级按以下描述划分:0级:模型的不确定性量化基于判断和经验或仅对模型进行了确定性分析。
39.1级:对模型中感兴趣领域中的不确定性进行了表征,但是没有对偶然不确定性与认知不确定性进行区分。对不确定性量化做出了一些假设,如:大多数概率密度函数被表征为高斯分布、不确定参数被认为与所有其他参数无关。
40.2级:对模型中存在的绝大多数不确定性进行了量化,并且对偶然不确定性和认知不确定性进行了分离,对不确定性量化进行评审,如非正式的评审或内部评审。
41.3级:对认知不确定性和偶然不确定性进行全面量化,对不确定性量化进行了专业评审,由外部专业审查小组完成。
42.数据适用性(da)成熟度等级按以下描述划分:0级:采集得到的信息数据不够完整和准确,难以用于模型的更新或进化,有用的数据占比较少,且飞机结构(实体)与数字孪生模型(虚体)之间数据交互的时延较长。
43.1级:采集得到的部分数据满足模型的更新或进化的要求,但是数据完整性及准确性有待完善。飞机结构(实体)与数字孪生模型(虚体)之间数据交互的时延较短。
44.2级:采集得到数据可以较为完整和准确,可以较好的满足模型更新或进化,对数据适用性进行了评审,例如非正式评审或内部评审。且飞机结构(实体)与数字孪生模型(虚体)之间的数据交互的时延短,可以满足飞机结构数字孪生系统大部分性能的实时性要求。
45.3级:采集得到数据完整、准确,且可靠性高,可以很好的满足模型更新或进化,并且对数据适用性进行了评审,由外部独立审查小组完成。且飞机结构(实体)与数字孪生模型(虚体)之间的数据交互时延极短,可以满足飞机结构数字孪生全部性能的实时性要求。
46.利用以上的划分标准,对飞机柔性翼结构数字孪生系统的空气动力学仿真模型和结构有限元仿真模型的可信性关键要素的成熟度等级进行划分,共2个单模型,每个单模型有6项关键要素,可得到12项关键要素的成熟度等级,得到的各项成熟度等级如下表1所示。
47.表1 关键要素的成熟度等级

48.步骤3:依据步骤2得到的可信性关键要素成熟度等级,通过下式计算单模型的预测能力成熟度等级:,式中,表示模型的预测能力成熟度等级,表示该模型下每种关键要素的成熟度等级,表示权重系数,范围在0~1之间,用户可以根据自己的需求以及模型本身的特性,设置每个可信性关键要素的权重系数。
49.本实施例中认为2个单模型的每个可信性关键要素同等重要,权重系数均取0.5,利用上式可得到空气动力学仿真模型的预测能力成熟度等级为2.167,2.167表示该模型成熟度等级在2和3之间,小数部分为0.167,即认为16.7%的概率为3级成熟度,83.3%的概率为2级成熟度;结构有限元仿真模型的成熟度等级为2,即认为100%的概率为2级成熟度。
50.预测能力成熟度模型是用来评估模型建模完整性与输出准确性的一种手段,模型预测能力成熟度等级有0级、1级、2级和3级4种成熟度等级,等级越高代表该模型完整性以及准确性越高。
51.模型的预测能力成熟度等级具有以下含义:0级:模型可信性较差,在这种情况下需要对模型进行修正与更新。
52.1级:模型可信性一般,基本能够满足仿真的需求,但是模型保真度一般,模型仿真得到的数据与真实数据之间存在一定差异。
53.2级:模型可信性良好,基本满足构建飞机结构数字孪生系统的要求。
54.3级:模型可信性极高,模型仿真得到数据与真实数据差异极小,满足构建飞机结构数字孪生系统的要求。
55.步骤4:现在得到了空气动力学仿真模型和结构有限元仿真模型两个单模型的预测能力成熟度等级。建立由可信性关键要素-模型预测能力成熟度-数字孪生系统可信性传递的贝叶斯网络。
56.根据参数关联模型分别研究两个单模型的不确定性对柔性翼结构数字孪生系统响应的敏感性和对可靠度的灵敏度,辨识出对数字孪生系统响应或可靠度影响较大的单模型(灵敏度分析)。最终通过下式确定贝叶斯网络的传递概率,可得到该飞机柔性翼结构数
字孪生系统的可信等级:,式中:表示通过灵敏度分析得到该模型的sobol灵敏度指标,取值范围在0~1之间,取值越大表示该模型发生变化或受不确定性影响时对数字孪生系统输出响应程度越大,取值越小则对输出响应程度越弱。此处两个单模型的均取0.5,最终可得出该飞机柔性翼结构数字孪生系统可信等级为2.084, 2.084表示该模型可信等级在2和3之间,小数部分为0.084,即认为8.4%的概率为3级可信等级,91.6%的概率为2级可信等级。
57.飞机结构数字孪生系统可信等级有0级、1级、2级和3级4种等级,具体含义如下:0级:飞机结构数字孪生系统可信性较差,在这种情况下可以认为数字孪生系统难以满足我们的需求,需要进行校准和更新。
58.1级:飞机结构数字孪生系统可信性一般,基本可以满足数字孪生的要求,但是无法实现高精度的预测。
59.2级:飞机结构数字孪生系统可信性较高,可以满足大部分对数字孪生系统的要求。
[0060] 3级:飞机结构数字孪生系统可信性极高,完全满足对数字孪生系统的要求,可以实现高精度的预测,达到实时预测和优化决策的目标。
[0061]
建立的可信性关键要素-模型预测能力成熟度-数字孪生系统可信性传递的贝叶斯网络为三层网络结构,第一层为要素层,包含每个单模型的可信性关键要素的成熟度等级,第二层为模型层,包含每个单模型的预测能力成熟度等级,第三层为系统层,包含该数字孪生系统的可信等级。要素层到模型层的传递是基于预测能力成熟度决策模型,利用可信性关键要素的成熟度等级以及权重系数确定条件传递概率;模型层到系统层的传递是基于步骤1得到的参数关联模型,通过灵敏度分析确定条件传递概率。如图3所示,三层贝叶斯网络层层递进,最终确定飞机柔性翼结构数字孪生系统的可信等级。依据飞机柔性翼结构数字孪生系统建立的贝叶斯网络如表2所示:
[0062]
表2 贝叶斯网络的内容。
[0063]
如图4所示,最终通过分析得到该飞机柔性翼结构数字孪生系统可信等级大概率为2级,表明该飞机柔性翼结构数字孪生系统可信性较高,如果不是对数字孪生系统要求特
别高,则认为该数字孪生系统是比较可靠的。
[0064]
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

技术特征:
1.一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法,其特征在于,包括以下步骤:步骤1,对要评估飞机结构数字孪生系统开展层次化解析,建立该数字孪生系统的参数关联模型,包括若干个单模型及单模型之间的联系;步骤2,分析所述参数关联模型的可信性关键要素,定义可信性关键要素成熟度等级划分标准,基于该标准得到每个单模型的可信性关键要素的成熟度等级;步骤3,根据步骤2得到的可信性关键要素的成熟度等级,构建单模型的预测能力成熟度决策模型,并依据可信性关键要素的权重系数,评估单模型的预测能力成熟度等级;步骤4,基于预测能力成熟度决策模型与贝叶斯网络,建立由可信性关键要素-模型预测能力成熟度-数字孪生系统可信性传递的贝叶斯网络,最终得到该飞机结构数字孪生系统的可信等级。2.根据权利要求1所述的一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法,其特征在于,所述参数关联模型为该飞机结构数字孪生系统各个模型之间的联系而构建的参数模型,包括该数字孪生系统的每个单模型及其之间的关联。3.根据权利要求2所述的一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法,其特征在于,所述步骤2具体为:定义每个单模型的可信性关键要素,并将其划分为k个成熟度等级,定义每个可信性关键要素的成熟度等级划分标准;分析要评估飞机结构的每个单模型的可信性关键要素,获得对应的成熟度等级。4.根据权利要求3所述的一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法,其特征在于,步骤3中,基于下式计算单模型的预测能力成熟度等级:,式中,表示模型的预测能力成熟度等级,范围在0~k-1之间,表示该模型下每种关键要素的成熟度等级,表示权重系数,范围在0~1之间,n为每个单模型的关键要素数量。5.根据权利要求4所述的一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法,其特征在于,步骤4中,基于单模型的预测能力成熟度等级以及对模型进行灵敏度分析获取的灵敏度指标来确定贝叶斯网络内部的条件传递概率,层层传递,最终得到该数字孪生系统可信等级。6.根据权利要求5所述的一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法,其特征在于,所述灵敏度分析为:根据参数关联模型分别研究每个单模型的不确定性对数字孪生系统响应的敏感性和对可靠度的灵敏度,将单模型按照对数字孪生系统响应或可靠度的影响程度排序。7.根据权利要求6所述的一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法,其特征在于,通过下式确定贝叶斯网络的传递概率,得到该飞机结构数字孪生系统的可信
等级:,式中,为可信等级,范围在0~k-1之间,m为单模型的数量,下标i表示第i个单模型,表示第i个单模型的预测能力成熟度等级,表示通过灵敏度分析得到的第i个单模型的sobol灵敏度指标,取值范围在0~1之间,取值越大表示该模型发生变化或受不确定性影响时对数字孪生系统输出响应程度越大,取值越小则对输出响应程度越弱。8.根据权利要求6所述的一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法,其特征在于,建立的可信性关键要素-模型预测能力成熟度-数字孪生系统可信性传递的贝叶斯网络为三层网络结构,第一层为要素层,包含每个单模型的可信性关键要素的成熟度等级,第二层为模型层,包含每个单模型的预测能力成熟度等级,第三层为系统层,包含该数字孪生系统的可信等级;要素层到模型层的传递是基于预测能力成熟度决策模型,利用可信性关键要素的成熟度等级以及权重系数确定条件传递概率;模型层到系统层基于灵敏度分析确定条件传递概率层层递进。9.根据权利要求8所述的一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法,其特征在于,若步骤4得到的可信等级不满足用户要求,则对飞机结构数字孪生系统进行修正与更新,直到可信等级满足用户要求时完成评估。

技术总结
本发明公开了一种基于贝叶斯网络的飞机结构数字孪生可信性评估方法。首先依据飞机结构数字孪生系统的模型结构开展模型层次解析并划分,建立数字孪生系统的参数关联模型。并通过分析单模型的几何模型保真度、物理及材料模型保真度、验证与确认、灵敏度分析、不确定性量化、数据适用性等6项可信性关键要素,构建单模型的预测能力成熟度决策模型并确定其预测能力成熟度等级。然后,建立由可信性关键要素-模型预测能力成熟度-数字孪生系统可信性传递的贝叶斯网络。最后,依据飞机结构数字孪生的参数关联模型,根据灵敏度分析获取的灵敏度指标确定贝叶斯网络的条件传递概率,量化飞机结构数字孪生系统的可信性。构数字孪生系统的可信性。构数字孪生系统的可信性。


技术研发人员:张浩宇 李维 李洪双
受保护的技术使用者:南京航空航天大学
技术研发日:2023.06.14
技术公布日:2023/7/18
版权声明

本文仅代表作者观点,不代表航家之家立场。
本文系作者授权航家号发表,未经原创作者书面授权,任何单位或个人不得引用、复制、转载、摘编、链接或以其他任何方式复制发表。任何单位或个人在获得书面授权使用航空之家内容时,须注明作者及来源 “航空之家”。如非法使用航空之家的部分或全部内容的,航空之家将依法追究其法律责任。(航空之家官方QQ:2926969996)

航空之家 https://www.aerohome.com.cn/

飞机超市 https://mall.aerohome.com.cn/

航空资讯 https://news.aerohome.com.cn/

分享:

扫一扫在手机阅读、分享本文

相关推荐