分布特性预测方法、装置和计算机可读存储介质与流程

文档序号:28264133发布日期:2021-12-31 17:45阅读:60来源:国知局
分布特性预测方法、装置和计算机可读存储介质与流程

1.本发明涉及建筑结构设计领域,尤其涉及分布特性预测方法、装置和计算机可读存储介质。


背景技术:

2.混凝土结构的质量缺陷判定,一直困扰着广大混凝土从业者。
3.目前对带缺陷钢管混凝土的承载力进行试验和非线性有限元分析研究的文献较多,但对随机生成缺陷的钢管混凝土的承载力的非线性有限元分析尚不多见。
4.现阶段建立的大多数有限元模型在考虑内部缺陷时,对于缺陷的位置、尺寸和形状等因素的设定人为干扰因素较多,建立的模型过于理想化,缺乏随机性,与实际情况有一定的出入,这也导致最终计算的承载力存在偏差。


技术实现要素:

5.本发明的主要目的在于提供一种分布特性预测方法、装置和计算机可读存储介质,旨在解决如何使得有限元建模方法对内部缺陷的布置具备随机性,进而使得有限元模型能够与实际情况更为吻合,减小最终计算的承载力存在的偏差的技术问题。
6.为实现上述目的,本发明提供一种分布特性预测方法,所述分布特性预测方法包括以下步骤:
7.建立钢管混凝土模型,基于所述钢管混凝土模型得到标准承载力;
8.选取分布模型,设置随机变量以通过所述分布模型生成随机缺陷类型;
9.在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型,基于所述随机缺陷钢管混凝土模型得到缺陷承载力;
10.将所述标准承载力和所述缺陷承载力进行对比,得到对比结果。
11.可选地,所述建立钢管混凝土模型的步骤之后包括:
12.为所述钢管混凝土模型中的材料选择本构模型和接触属性,对所述钢管混凝土模型划分网格,施加边界条件和外荷载,以得到标准钢管混凝土模型。
13.可选地,所述基于所述钢管混凝土模型得到标准承载力的步骤包括:
14.基于所述标准钢管混凝土模型通过预设算法计算得到标准承载力。
15.可选地,所述选取分布模型的步骤之前包括:
16.确定所述标准钢管混凝土模型的计算区域,在所述计算区域内,确定缺陷中心坐标、缺陷形状和缺陷体积。
17.可选地,所述确定缺陷中心坐标的步骤包括:
18.选取位置随机函数以确定位置随机的缺陷中心坐标;
19.所述确定缺陷形状的步骤包括:
20.若确定缺陷形状为多边形,通过第一预设方法确定多边形的边数,进而确定缺陷形状;
21.若确定缺陷形状为椭圆形,通过第二预设方法确定椭圆的扁平程度,进而确定缺陷形状;
22.所述确定缺陷体积的步骤包括:
23.通过随机生成的缺陷率与所述标准钢管混凝土模型的体积,确定缺陷体积。
24.可选地,所述设置随机变量以通过所述分布模型生成随机缺陷类型的步骤包括:
25.将所述缺陷体积、所述缺陷形状和所述缺陷中心坐标设置为随机变量;
26.将所述随机变量代入所述分布模型以生成随机缺陷类型。
27.可选地,所述在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型的步骤包括:
28.在所述标准钢管混凝土模型内部生成与所述随机缺陷类型相对应的内部缺陷,对所述内部缺陷选定范围内的单元自动删除或赋予预设参数,以得到随机缺陷钢管混凝土模型。
29.可选地,所述将所述标准承载力和所述缺陷承载力进行对比,得到对比结果的步骤包括:
30.将所述标准承载力和所述缺陷承载力进行对比,得到缺陷分布特性数据,将所述缺陷分布特性数据作为对比结果,保存所述缺陷分布特性数据和所述缺陷分布特性数据对应的所述随机缺陷钢管混凝土模型。
31.此外,为实现上述目的,本发明还提供一种装置,所述装置包括:存储器、处理器及存储在所述存储器上并可在所述处理器上运行的分布特性预测程序,所述分布特性预测程序被所述处理器执行时实现如上所述的分布特性预测方法的步骤。
32.此外,为实现上述目的,本发明还提供一种计算机可读存储介质,所述计算机可读存储介质上存储有分布特性预测程序,所述分布特性预测程序被处理器执行时实现如上所述的分布特性预测方法的步骤。
33.本发明提出一种分布特性预测方法、装置和计算机可读存储介质。所述分布特性预测方法通过建立钢管混凝土模型,基于所述钢管混凝土模型得到标准承载力,通过选取分布模型,设置随机变量以通过所述分布模型生成随机缺陷类型,使得有限元建模方法对内部缺陷的布置具备随机性;通过在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型,基于所述随机缺陷钢管混凝土模型得到缺陷承载力;将所述标准承载力和所述缺陷承载力进行对比,得到对比结果,进而使得有限元模型能够与实际情况更为吻合,减小最终计算的承载力存在的偏差。
附图说明
34.图1是本发明实施例方案涉及的硬件运行环境的终端结构示意图;
35.图2为本发明分布特性预测方法第一实施例的流程示意图。
36.本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
37.应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
38.本发明实施例的主要解决方案是:一种分布特性预测方法,所述分布特性预测方
法包括以下步骤:
39.建立钢管混凝土模型,基于所述钢管混凝土模型得到标准承载力;
40.选取分布模型,设置随机变量以通过所述分布模型生成随机缺陷类型;
41.在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型,基于所述随机缺陷钢管混凝土模型得到缺陷承载力;
42.将所述标准承载力和所述缺陷承载力进行对比,得到对比结果。
43.由于混凝土结构的质量缺陷判定,一直困扰着广大混凝土从业者。目前对带缺陷钢管混凝土的承载力进行试验和非线性有限元分析研究的文献较多,但对随机生成缺陷的钢管混凝土的承载力的非线性有限元分析尚不多见。现阶段建立的大多数有限元模型在考虑内部缺陷时,对于缺陷的位置、尺寸和形状等因素的设定人为干扰因素较多,建立的模型过于理想化,缺乏随机性,与实际情况有一定的出入,这也导致最终计算的承载力存在偏差。
44.本发明提供一种分布特性预测方法,通过建立钢管混凝土模型,基于所述钢管混凝土模型得到标准承载力,通过选取分布模型,设置随机变量以通过所述分布模型生成随机缺陷类型,使得有限元建模方法对内部缺陷的布置具备随机性;通过在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型,基于所述随机缺陷钢管混凝土模型得到缺陷承载力;将所述标准承载力和所述缺陷承载力进行对比,得到对比结果,进而使得有限元模型能够与实际情况更为吻合,减小最终计算的承载力存在的偏差。
45.如图1所示,图1是本发明实施例方案涉及的硬件运行环境的终端结构示意图。
46.本发明实施例终端可以是pc,也可以是平板电脑、便携计算机等能够支持abaqus有限元软件的终端设备。
47.如图1所示,该终端可以包括:处理器1001,例如cpu,网络接口1004,用户接口1003,存储器1005,通信总线1002。其中,通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏(display)、输入单元比如键盘(keyboard),可选用户接口1003还可以包括标准的有线接口、无线接口。网络接口1004可选的可以包括标准的有线接口、无线接口(如wi

fi接口)。存储器1005可以是高速ram存储器,也可以是稳定的存储器(non

volatile memory),例如磁盘存储器。存储器1005可选的还可以是独立于前述处理器1001的存储装置。
48.可选地,终端还可以包括摄像头、rf(radio frequency,射频)电路,传感器、音频电路、wifi模块等等。其中,传感器比如光传感器、运动传感器以及其他传感器。具体地,光传感器可包括环境光传感器及接近传感器,其中,环境光传感器可根据环境光线的明暗来调节显示屏的亮度,接近传感器可在移动终端移动到耳边时,关闭显示屏和/或背光。作为运动传感器的一种,重力加速度传感器可检测各个方向上(一般为三轴)加速度的大小,静止时可检测出重力的大小及方向,可用于识别移动终端姿态的应用(比如横竖屏切换、相关游戏、磁力计姿态校准)、振动识别相关功能(比如计步器、敲击)等;当然,移动终端还可配置陀螺仪、气压计、湿度计、温度计、红外线传感器等其他传感器,在此不再赘述。
49.本领域技术人员可以理解,图1中示出的终端结构并不构成对终端的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件布置。
50.如图1所示,作为一种计算机存储介质的存储器1005中可以包括操作系统、网络通
信模块、用户接口模块以及分布特性预测程序。
51.在图1所示的终端中,网络接口1004主要用于连接后台服务器,与后台服务器进行数据通信;用户接口1003主要用于连接客户端(用户端),与客户端进行数据通信;而处理器1001可以用于调用存储器1005中存储的分布特性预测程序,并执行以下操作:
52.建立钢管混凝土模型,基于所述钢管混凝土模型得到标准承载力;
53.选取分布模型,设置随机变量以通过所述分布模型生成随机缺陷类型;
54.在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型,基于所述随机缺陷钢管混凝土模型得到缺陷承载力;
55.将所述标准承载力和所述缺陷承载力进行对比,得到对比结果。
56.进一步地,处理器1001可以调用存储器1005中存储的分布特性预测程序,还执行以下操作:
57.为所述钢管混凝土模型中的材料选择本构模型和接触属性,对所述钢管混凝土模型划分网格,施加边界条件和外荷载,以得到标准钢管混凝土模型。
58.进一步地,处理器1001可以调用存储器1005中存储的分布特性预测程序,还执行以下操作:
59.所述基于所述钢管混凝土模型得到标准承载力的步骤包括:
60.基于所述标准钢管混凝土模型通过预设算法计算得到标准承载力。
61.进一步地,处理器1001可以调用存储器1005中存储的分布特性预测程序,还执行以下操作:
62.确定所述标准钢管混凝土模型的计算区域,在所述计算区域内,确定缺陷中心坐标、缺陷形状和缺陷体积。
63.进一步地,处理器1001可以调用存储器1005中存储的分布特性预测程序,还执行以下操作:
64.所述确定缺陷中心坐标的步骤包括:
65.选取位置随机函数以确定位置随机的缺陷中心坐标;
66.所述确定缺陷形状的步骤包括:
67.若确定缺陷形状为多边形,通过第一预设方法确定多边形的边数,进而确定缺陷形状;
68.若确定缺陷形状为椭圆形,通过第二预设方法确定椭圆的扁平程度,进而确定缺陷形状;
69.所述确定缺陷体积的步骤包括:
70.通过随机生成的缺陷率与所述标准钢管混凝土模型的体积,确定缺陷体积。
71.进一步地,处理器1001可以调用存储器1005中存储的分布特性预测程序,还执行以下操作:
72.所述设置随机变量以通过所述分布模型生成随机缺陷类型的步骤包括:
73.将所述缺陷体积、所述缺陷形状和所述缺陷中心坐标设置为随机变量;
74.将所述随机变量代入所述分布模型以生成随机缺陷类型。
75.进一步地,处理器1001可以调用存储器1005中存储的分布特性预测程序,还执行以下操作:
76.所述在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型的步骤包括:
77.在所述标准钢管混凝土模型内部生成与所述随机缺陷类型相对应的内部缺陷,对所述内部缺陷选定范围内的单元自动删除或赋予预设参数,以得到随机缺陷钢管混凝土模型。
78.进一步地,处理器1001可以调用存储器1005中存储的分布特性预测程序,还执行以下操作:
79.所述将所述标准承载力和所述缺陷承载力进行对比,得到对比结果的步骤包括:
80.将所述标准承载力和所述缺陷承载力进行对比,得到缺陷分布特性数据,将所述缺陷分布特性数据作为对比结果,保存所述缺陷分布特性数据和所述缺陷分布特性数据对应的所述随机缺陷钢管混凝土模型。
81.参照图2,本发明第一实施例提供一种分布特性预测方法,所述分布特性预测方法包括:
82.步骤s10,建立钢管混凝土模型,基于所述钢管混凝土模型得到标准承载力;
83.本实施例中,步骤s10中所述建立钢管混凝土模型的步骤之后包括:
84.为所述钢管混凝土模型中的材料选择本构模型和接触属性,对所述钢管混凝土模型划分网格,施加边界条件和外荷载,以得到标准钢管混凝土模型;
85.步骤s10中所述基于所述钢管混凝土模型得到标准承载力的步骤包括:
86.基于所述标准钢管混凝土模型通过预设算法计算得到标准承载力;
87.需要说明的是,本实施例中,执行主体为pc,也可以是平板电脑、便携计算机等能够支持abaqus有限元软件的终端设备,本实施例对此不加以限制。所述abaqus是一套功能强大的工程模拟的有限元软件;所述本构模型是用于表示某一种材料的本构关系的物理模型,按性质可分为弹性模型、刚塑性模型、弹塑性模型、黏弹性模型、黏塑性模型和弹黏塑性模型,与之相应的性质常由模型参数表示,包括土的弹性模型、土的超弹性模型、土的次弹性模型、土的粘弹性模型、土的弹塑性模型、土的粘弹塑性模型、土的内蕴时间塑性模型等;所述接触属性是所述钢管混凝土模型的属性,包括刚性、柔性等;所述预设算法为abaqus有限元软件针对钢管混凝土模型预置的算法,也可以由用户自行修改后使用。
88.具体实现中,在abaqus有限元软件中建立钢管混凝土模型,给模型中的两种材料分别选择合适的本构模型和接触属性,即对钢管选择合适的本构模型和接触属性和对混凝土选择合适的本构模型和接触属性,然后对整个模型划分网格,施加边界条件和外荷载,以完善钢管混凝土模型得到标准钢管混凝土模型,然后基于所述标准钢管混凝土模型通过预设算法计算得到不带缺陷的钢管混凝土承载力,即所述标准承载力。
89.本实施例中,步骤s20之前包括:
90.确定所述标准钢管混凝土模型的计算区域,在所述计算区域内,确定缺陷中心坐标、缺陷形状和缺陷体积。
91.其中,所述确定缺陷中心坐标的步骤包括:
92.选取位置随机函数以确定位置随机的缺陷中心坐标;
93.所述确定缺陷形状的步骤包括:
94.若确定缺陷形状为多边形,通过第一预设方法确定多边形的边数,进而确定缺陷
形状;
95.若确定缺陷形状为椭圆形,通过第二预设方法确定椭圆的扁平程度,进而确定缺陷形状;
96.所述确定缺陷体积的步骤包括:
97.通过随机生成的缺陷率与所述标准钢管混凝土模型的体积,确定缺陷体积。
98.需要说明的是,所述计算区域即为所述标准钢管混凝土模型的内部区域。
99.可以理解的是,缺陷需要在所述标准钢管混凝土模型的内部模拟生成才有参考价值,所述缺陷中心坐标用于确定缺陷的位置。
100.具体实现中,在所建模型的计算区域内,通过选取合适的位置随机函数确定一系列位置随机的缺陷中心坐标,以确定缺陷的位置。
101.在确定缺陷的形状时,当确定缺陷为多变形时,以缺陷中心坐标为原点坐标建立平面直角坐标系,在原点坐标附近随机确定节点数,以随机确定的节点个数来确定多边形的边数,节点的位置以节点与原点坐标的连线与x轴的夹角来决定,多边形的大小以节点与原点坐标连线的长度来决定,即所述第一预设方法;当缺陷形状为椭圆形时,选取合适的分布模型,以长轴和短轴的尺寸为随机变量来确定椭圆的扁平程度,即所述第二预设方法,从而确定其形状。
102.在确定内部缺陷的大小时,通过随机生成的缺陷率结合钢管混凝土模型的体积大小,以确定内部缺陷的体积大小。
103.步骤s20,选取分布模型,设置随机变量以通过所述分布模型生成随机缺陷类型;
104.本实施例中,步骤s20中所述设置随机变量以通过所述分布模型生成随机缺陷类型的步骤包括:
105.将所述缺陷体积、所述缺陷形状和所述缺陷中心坐标设置为随机变量;
106.将所述随机变量代入所述分布模型以生成随机缺陷类型。
107.需要说明的是,所述分布模型为软件内预置的参数模型,将前述步骤中得到的所述缺陷体积、所述缺陷形状和所述缺陷中心坐标作为随机变量代入该模型后即可生成多种不同类型的随机缺陷。
108.步骤s30,在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型,基于所述随机缺陷钢管混凝土模型得到缺陷承载力;
109.本实施例中,步骤s30中所述在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型的步骤包括:
110.在所述标准钢管混凝土模型内部生成与所述随机缺陷类型相对应的内部缺陷,对所述内部缺陷选定范围内的单元自动删除或赋予预设参数,以得到随机缺陷钢管混凝土模型。
111.可以理解的是,缺陷需要生成在标准模型内部以生成缺陷模型,即所述随机缺陷钢管混凝土模型。
112.具体实现中,基于步骤s20中得到多种不同的随机缺陷类型,进行大量抽样试验,在已建成的钢管混凝土模型内部生成相对应的缺陷,对于缺陷选定范围内的单元自动删除或者赋予较低的材料参数,即所述预设参数。计算各种缺陷对应下的钢管混凝土承载力,即所述缺陷承载力。
113.步骤s40,将所述标准承载力和所述缺陷承载力进行对比,得到对比结果。
114.本实施例中,步骤s40包括:
115.将所述标准承载力和所述缺陷承载力进行对比,得到缺陷分布特性数据,将所述缺陷分布特性数据作为对比结果,保存所述缺陷分布特性数据和所述缺陷分布特性数据对应的所述随机缺陷钢管混凝土模型。
116.可以理解的是,通过对上述步骤中得到的承载力数据进行分析对比,即将所述标准承载力与所述缺陷承载力数据进行分析对比,可得到缺陷率、缺陷位置和缺陷形状等因素对钢管混凝土承载力的影响程度,即所述缺陷分布特性数据,将所述缺陷分布特性数据和所述缺陷分布特性数据对应的所述随机缺陷钢管混凝土模型保存,即可得到大量具有随机缺陷类型的钢管混凝土模型及其缺陷分布特性数据,在此基础上还能通过对比分析得到对钢管混凝土承载力影响最大的缺陷类型。
117.在本实施例中提出了一种分布特性预测方法,通过在abaqus有限元软件中建立钢管混凝土模型,给模型中的两种材料分别选择合适的本构模型和接触属性,然后对整个模型划分网格,施加边界条件和外荷载,以完善钢管混凝土模型得到标准钢管混凝土模型,进而基于所述钢管混凝土模型得到了标准承载力;通过选取分布模型,将缺陷体积、缺陷形状和缺陷中心坐标设置为随机变量,将所述随机变量代入所述分布模型以生成随机缺陷类型,使得有限元建模方法对内部缺陷的布置具备随机性;通过在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型,基于所述随机缺陷钢管混凝土模型得到缺陷承载力;将所述标准承载力和所述缺陷承载力进行对比,得到对比结果,由于本实施例中进行了大量抽样试验,故能够得到具备随机性的大量缺陷承载力与标准承载力的对比结果,进而使得有限元模型能够与实际情况更为吻合,减小最终计算的承载力存在的偏差。
118.此外,本发明实施例还提出一种计算机可读存储介质,所述计算机可读存储介质上存储有分布特性预测程序,所述分布特性预测程序被处理器执行时实现如下操作:
119.建立钢管混凝土模型,基于所述钢管混凝土模型得到标准承载力;
120.选取分布模型,设置随机变量以通过所述分布模型生成随机缺陷类型;
121.在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型,基于所述随机缺陷钢管混凝土模型得到缺陷承载力;
122.将所述标准承载力和所述缺陷承载力进行对比,得到对比结果。
123.进一步地,所述分布特性预测程序被处理器执行时还实现如下操作:
124.为所述钢管混凝土模型中的材料选择本构模型和接触属性,对所述钢管混凝土模型划分网格,施加边界条件和外荷载,以得到标准钢管混凝土模型。
125.进一步地,所述分布特性预测程序被处理器执行时还实现如下操作:
126.所述基于所述钢管混凝土模型得到标准承载力的步骤包括:
127.基于所述标准钢管混凝土模型通过预设算法计算得到标准承载力。
128.进一步地,所述分布特性预测程序被处理器执行时还实现如下操作:
129.确定所述标准钢管混凝土模型的计算区域,在所述计算区域内,确定缺陷中心坐标、缺陷形状和缺陷体积。
130.进一步地,所述分布特性预测程序被处理器执行时还实现如下操作:
131.所述确定缺陷中心坐标的步骤包括:
132.选取位置随机函数以确定位置随机的缺陷中心坐标;
133.所述确定缺陷形状的步骤包括:
134.若确定缺陷形状为多边形,通过第一预设方法确定多边形的边数,进而确定缺陷形状;
135.若确定缺陷形状为椭圆形,通过第二预设方法确定椭圆的扁平程度,进而确定缺陷形状;
136.所述确定缺陷体积的步骤包括:
137.通过随机生成的缺陷率与所述标准钢管混凝土模型的体积,确定缺陷体积。
138.进一步地,所述分布特性预测程序被处理器执行时还实现如下操作:
139.所述设置随机变量以通过所述分布模型生成随机缺陷类型的步骤包括:
140.将所述缺陷体积、所述缺陷形状和所述缺陷中心坐标设置为随机变量;
141.将所述随机变量代入所述分布模型以生成随机缺陷类型。
142.进一步地,所述分布特性预测程序被处理器执行时还实现如下操作:
143.所述在所述钢管混凝土模型内部生成所述随机缺陷类型,得到随机缺陷钢管混凝土模型的步骤包括:
144.在所述标准钢管混凝土模型内部生成与所述随机缺陷类型相对应的内部缺陷,对所述内部缺陷选定范围内的单元自动删除或赋予预设参数,以得到随机缺陷钢管混凝土模型。
145.进一步地,所述分布特性预测程序被处理器执行时还实现如下操作:
146.所述将所述标准承载力和所述缺陷承载力进行对比,得到对比结果的步骤包括:
147.将所述标准承载力和所述缺陷承载力进行对比,得到缺陷分布特性数据,将所述缺陷分布特性数据作为对比结果,保存所述缺陷分布特性数据和所述缺陷分布特性数据对应的所述随机缺陷钢管混凝土模型。
148.需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者系统不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情况下,由语句“包括一个
……”
限定的要素,并不排除在包括该要素的过程、方法、物品或者系统中还存在另外的相同要素。
149.上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
150.通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在如上所述的一个存储介质(如rom/ram、磁碟、光盘)中,包括若干指令用以使得一台终端设备(可以是手机,计算机,服务器,空调器,或者网络设备等)执行本发明各个实施例所述的方法。
151.以上仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1