弹性与塑性力学基础 第六章塑性力学解题方法及应用举例ppt课件.ppt

上传人:小飞机 文档编号:1920774 上传时间:2022-12-26 格式:PPT 页数:52 大小:2.55MB
返回 下载 相关 举报
弹性与塑性力学基础 第六章塑性力学解题方法及应用举例ppt课件.ppt_第1页
第1页 / 共52页
弹性与塑性力学基础 第六章塑性力学解题方法及应用举例ppt课件.ppt_第2页
第2页 / 共52页
弹性与塑性力学基础 第六章塑性力学解题方法及应用举例ppt课件.ppt_第3页
第3页 / 共52页
弹性与塑性力学基础 第六章塑性力学解题方法及应用举例ppt课件.ppt_第4页
第4页 / 共52页
弹性与塑性力学基础 第六章塑性力学解题方法及应用举例ppt课件.ppt_第5页
第5页 / 共52页
点击查看更多>>
资源描述

《弹性与塑性力学基础 第六章塑性力学解题方法及应用举例ppt课件.ppt》由会员分享,可在线阅读,更多相关《弹性与塑性力学基础 第六章塑性力学解题方法及应用举例ppt课件.ppt(52页珍藏版)》请在三一办公上搜索。

1、,弹性与塑性力学基础,第 六 章塑性力学解题方法及应用举例,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,1、塑性力学问题求解现状 (1) 在塑性状态物体内应力的大小与分布求解比较弹性状态困难; (2) 非线性塑性应力应变关系方程; (3) 联解平衡方程和屈服准则,补充必要的物理方程和几何方程,在 一定的边界条件下可以求得变形体内的应力大小及分布; (4) 某些特殊情况下能够数学解析,一般空间问题,数学上极其困难, 甚至不可能解。2、塑性力学问题求解方法 (1) 塑性理论基础上,引进各种简化假设,提出求解的近似解析方法; (2) 主应力法(切块法)、滑移线法和上限法,数值模拟方法。,

2、6-1 平衡微分方程和屈服准则联立求解及其应用 6.1.1 求解方法 6.1.2 受内压塑性圆筒及受内拉的塑性圆环应力计算 6-2 圆柱体镦粗变形力计算的主应力法 6.2.1 主应力法的特点 6.2.2 主应力法的要点平行模板间圆柱体镦粗 6.2.3 平行模板间圆柱体镦粗,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法 6.3.2 亨盖(Hencky)应力方程 6.3.3 常见的滑移线场类型 6.3.4 用滑移线场理论求解的二个典型问题 6-4 上限法及其在平面问题中的应用 6.4.1 上限法的概

3、念 6.4.2 上限法原理 6.4.3 上限法在平面变形问题中的应用,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-1 平衡微分方程和屈服准则联立求解及其应用 6.1.1 求解方法平衡微分方程和屈服准则联立求解,求出物体塑性变形时的应力分布;联解过程积分常数根据自由表面和接触面上的边界条件确定。一般只能求解平面轴对称等简单塑性问题。,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-1 平衡微分方程和屈服准则联立求解及其应用 6.1.2 受内压塑性圆筒及受内拉的塑性圆环应力计算,弹性与塑性力学基础,受力状态:圆筒的内壁作用有均匀压力p几何尺寸:筒的尺寸如图所示变形类型:

4、平面应变(圆筒很长,相当于 压力容器、管道、挤压凹模等) 轴对称平面问题,、,、,应力分析: rz、r为零 、 r为主应力,仅随 r 变化;平衡微分方程: (6-1),第六章 塑性力学解题方法及应用举例,弹性与塑性力学基础,6-1 平衡微分方程和屈服准则联立求解及其应用 6.1.2 受内压塑性圆筒及受内拉的塑性圆环应力计算解题方法 根据Mises屈服准则,有 - r = s (6-2) 式中 为中间主应力影响系数,对于平面应变问题 , 代入式(6-1),得 积分得 利用边界条件确定积分常数C,当r=b, r =0,则,第六章 塑性力学解题方法及应用举例,弹性与塑性力学基础,6-1 平衡微分方程

5、和屈服准则联立求解及其应用 6.1.2 受内压塑性圆筒及受内拉的塑性圆环应力计算联立式(6-2)得到塑性圆筒的应力解 (6-3) 分析上式知,截面的 总为拉应力, r总为压应力。 当r=a时,有最大的压力p,所以,第六章 塑性力学解题方法及应用举例,弹性与塑性力学基础,6-1 平衡微分方程和屈服准则联立求解及其应用 6.1.2 受内压塑性圆筒及受内拉的塑性圆环应力计算对于圆环受拉问题,平衡微分方程法依旧,由于是平面应力问题,屈服准则为 r - = s ,可取=1.1,将边界条件代入后可得 (6-4),第六章 塑性力学解题方法及应用举例,弹性与塑性力学基础,弹性与塑性力学基础,6-2 圆柱体镦粗

6、变形力计算的主应力法 6.2.1 主应力法的要点 将问题近似地按轴对称问题或平面问题来处理 假设非接触面上仅有均布的正应力即主应力 假设接触面上的正应力即为主应力(即忽略摩擦切应力的影响) max - min = s 将上述的平衡方程与近似屈服准则联解,以求接触面上的应力分布,这就是主应力法。由于该方法需要截取基元块,又形象地称为切块法。,第六章 塑性力学解题方法及应用举例,弹性与塑性力学基础,弹性与塑性力学基础,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 适合主应力法求解物体几何上轴对称,受载荷也是轴对称的,属于轴对称问题,,第六章 塑性力学解题方法及应用举例,

7、弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 解题步骤,切取基元块 列平衡方程(沿r向) 整理并略去高次项得平衡微分方程 (6-7),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 解题步骤,找 r与的关系 可以从 r与的关系和应力应变关系式判别 实心圆柱镦粗的径向应变为 切向应变为 两者相等,根据应力应变关系理论必然有(6-8) 将式(6-8)代入式(6-7),可得 (6-9),弹性与塑性

8、力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 解题步骤,代入边界摩擦条件 边界上可能存在的摩擦条件为 设边界上选最大值,即 (6-10) 超过此数值工件与模板间的摩擦由剪切所代替 将式(6-10)代入式(6-9),可得 (6-11),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 解题步骤,写出屈服准则的表达式 由应变状态可见 r = 0 z (- z ) 此时的屈服准则 max - min = s

9、即视 r 、 z为主应力,有 (- r ) - (- z ) = s 即: z - r = s (6-12),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 解题步骤,从式(6-12)微分可得 (6-13) 由此可见只要=常数,式(6-13)总是成立的联立求解 将屈服准则式(6-13)与微分方程式(6-11)联解(6-11) 得(6-14),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 解题

10、步骤,积分上式,得(6-15) 定积分常数,当r=d/2 时, r =0 代入式(6-12)得 z = s,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 解题步骤,求接触面上压力分布公式 将C代入式(6-15)得圆柱体镦粗压力分布公式(6-16) 若边界摩擦=z (6-17) 若边界摩擦=s(6-18),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 解题步骤,镦粗单位压力分布,=s/2,=

11、s,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 解题步骤,求总变形力 沿接触平面积分即可得总变形力(6-19)(6-20) (6-21),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗 解题步骤,求平均单位力 =s/2时(6-22) 如热锻等 =z时(6-23) 摩擦系数较小的冷变形情况 =s时(6-24) 摩擦系数较小的冷变形情况,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学

12、解题方法及应用举例,6-2 圆柱体镦粗变形力计算的主应力法 6.2.3 平行模板间圆柱体镦粗解题步骤,随着及d/h增加,平均单位变形力将迅速增加有相对厚度越小平均单位变形力越大的概念 =s/2时(6-22) =z时(6-23) =s时(6-21),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法滑移线的定义 塑性变形体内各点最大剪应力的轨迹称为滑移线 由于最大剪应力成对正交,因此滑移线在变形体内成两族互相正交的线网,组成所谓滑移线场。,弹性与塑性力学基础,弹性与塑性力学基础,第六

13、章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法滑移线法 移线法就是利用滑移线场特性来求解变形问题严格地说仅适用于处理理想刚塑性体的平面应变问题,但在一定条件下,也可推广到平面应力和轴对称问题以及硬化材料;与其它方法相比,在数学上比较严谨,理论上比较完整;可以计算变形力、确定塑性变形区内的应力分布和速度分布、接触面上的应力分布及等静压迹线等; 构成滑移线的网格为在计算机上采用数值方法求解提供了自然单元,为该方法在工程计算上应用提供了方便。,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3

14、滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法滑移线的基本概念 塑性变形体(或变形区)内任一点的应力状态如图所示,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法滑移线的基本概念 塑性变形体(或变形区)内任一点的应力状态如图所示,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法滑移线的基本概念 两个彼此正交的最大剪应力面(即点和点所

15、代表的物理平面)与塑性流动平面相垂直,最大剪应力为,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法滑移线的基本概念 作用于最大剪应力面上的正应力13恰等于平均应力m或中间主应力2 ,即,任一点应力状态可用静水压(平均 应力)与最大剪切力K相叠加来表 示,即有 (6-25),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法滑移线的基本概念 由应力莫尔圆可知,应力分量

16、x、y、xy与静水压m和最大切应力K的关系(6-26),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法滑移线的基本概念 塑性变形体内每一点都能找到一对正交的最大剪应力方向,将无限接近的最大剪应力方向连接起来,即得两族正交的曲线,线上任一点的切线方向即为该点最大剪应力方向。,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法滑移线的基本概念,与线形成右手坐标系轴,则最

17、大主应力1位于第一与第三象限,线两旁最大剪应力顺时针方向,而线两旁最大剪应力逆时针方向。线的切线方向与Ox轴的夹角以表示,并规定以Ox轴为角的度量起始线,逆时针旋转形成的角为正值,顺时针旋转形成的角为负值。,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.1 滑移线的定义与滑移线法滑移线的基本概念 按应力边界条件判断边界上某点1与3方向(2必垂直于物理平面),确定最大切应力K的方向,根据最大切应力K确定该线的族别。,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线

18、场概念及其在平冲头镦粗半无限体中的应用 6.3.2 亨盖(Hencky)应力方程对于理想刚塑性材料中处于塑性平面应变状态下的各点,其应力状态不同实质只是各点的平均正应力m不同,而各点处的最大切应力K为材料常数都是一样的(无强化的塑性变形)各点的应力状态不同只是归结为各点的平均正应力m不同,沿滑移线上各点的平均正应力的变化规律即是亨盖应力方程(沿线) (6-27)(沿线) (6-28)亨盖方程描述了平均正应力m沿滑移线的变化规律,沿某线(或线)上任一点平均正应力m减去(或加上)该点的倾角同2K的乘积为一常数。从一条滑移线转到同族另一条滑移线上时,这一常数才改变。,弹性与塑性力学基础,弹性与塑性力

19、学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.2 亨盖(Hencky)应力方程各点不同的应力状态反映在应力平面-上的几何图形是一系列直径均为2K的应力莫尔圆。只是由于各点所对应的单元体上的平均正应力m不同(即圆心不同),而使诸应力圆沿轴分布,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.2 亨盖(Hencky)应力方程滑移线场的若干性质直接或间接由亨盖方程导出,沿线: =常数,沿线: =常数,同一条滑移线上任意两点间平均正应力 m的变化与该滑移

20、线在这两点间切线的 转角成正比。滑移线方向变化越剧烈(转 角越大),则平均正应力变化亦越大。直 线滑移线沿线各点的平均正应力不变。,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.3 常见的滑移线场类型两族正交的直线构成的滑移线场(塑性区中的均匀应力场)一族为直线,另一族为曲线所构成的中心扇形场与无中心扇形场(亦称简单滑移线场),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.3 常见的滑移线场类型其滑移线场为两族正

21、交的对数螺线所构成一族为直线(当圆弧边界面为自由表面或作用有均匀载荷时)由两等半径或不等半径圆弧所决定的有心扇形场。特别是等半径的有心扇形场是求解金属塑性加工问题中一种最常见、最重要的滑移线场类型,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 受内压无限长厚壁圆筒问题 无限长厚壁圆筒内表面沿全长受均匀径向内压p作用,计算使整个壁厚都进入塑性状态时p值的大小。设圆筒外径D2b,内径d=2a,材料屈服应力为m,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题

22、方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 受内压无限长厚壁圆筒问题(轴对称塑性平面应变问题),r=r=0,筒壁上任一点的径向正应力r和切向正应力都是主应力,主应力迹线为一系列的同心圆与径向线,又因为滑移线为最大切应力的迹线,故每一条滑移线与径向线相交成,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 受内压无限长厚壁圆筒问题(轴对称塑性平面应变问题),取单元体 积分得因为所以即 两族正交的对数螺

23、线方程,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 受内压无限长厚壁圆筒问题(轴对称塑性平面应变问题),根据边界上B点的应力状态,判断哪一族是线族 AB线为族滑移线,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 受内压无限长厚壁圆筒问题(轴对称塑性平面应变问题),求边界上点A的m和值 在内表面A点: 取r=a,r =-p 屈服准则: 1

24、 - 3= - r=2k 所以: =2K p 内表面各点的平均正应力: mA =(1+3)/2=(+r)/2=k-p,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 受内压无限长厚壁圆筒问题(轴对称塑性平面应变问题),求边界上点B的m和值 在外表面B点: 取r=b,r =0 屈服准则: 1 - 3= - r=2k 所以: =2K 外表面各点的平均正应力: mB =(1+3)/2=(+r)/2=k,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举

25、例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 受内压无限长厚壁圆筒问题(轴对称塑性平面应变问题),由Hencky应力方程,沿线AB有 将求得的mA、mB、A、B值代入 所以:,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 受内压无限长厚壁圆筒问题(轴对称塑性平面应变问题),由Trecas屈服准则,K= s/2由Mises屈服准则,K= 以上结果与前述平衡方程与屈服准则联解结果相同,弹性与塑性力学基础,弹性与塑

26、性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 平冲头压入半无限高坯料问题 这一问题有两种滑移线场解即普朗特场与希尔场(两种场p值相同) 普朗特场由三个等腰直角三角形的均匀场与二个的中心扇形场组成 希尔场由四个小等腰直角三角形与二个的中心扇形场组成,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 平冲头压入半无限高坯料问题 按照定族规则,从A点或B点分析,判断AHB线是

27、族滑移线 求边界上点B的m和值 :B=/ 4, 1=0 屈服准则: 1 - 3= 2k 所以: 3 = -2K B点的平均正应力: mB =(1+3)/2= -k,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 平冲头压入半无限高坯料问题 按照定族规则,从A点或B点分析,判断AHB线是族滑移线 求边界上点A的m和值 :A= -4, 3= -p 屈服准则: 1 - 3= 2k 所以: 1 =2K - p A点的平均正应力: mA =(1+3)/2=k - p,弹性与塑性

28、力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-3 滑移线场概念及其在平冲头镦粗半无限体中的应用 6.3.4 用滑移线场理论求解的二个典型问题 平冲头压入半无限高坯料问题 因AHB是族滑移线,由亨盖方程可得 mA 2kA= mB 2kB ; mA = mB 2k(B A) 所以: K p = -K-2K(/ 4 +/ 4 ) p = K( + 2 ) 5.14K 平冲头单位长度上的极限压力: P = 2b1p=2bK( + 2 ),弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-4 上限法及其在平面问题中的应用 6.4.1 上限法的概念上限

29、法是确定金属塑性变形时近似载荷的一种有效方法,确定的载荷(上限解)大于或等于真实载荷,也称为“高估法”,好的上限解是略高于真实值。金属塑性成形上限解能确保塑性变形发生,在选择成形设备与设计模具时,为了安全可靠,一般也取比真实变形力较大的载荷,上限法正适合于这一要求。上限法的力学基础是虚功原理。上限法难以克服的缺点,高估量是未知的,不能给出变形区的应力场。,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-4 上限法及其在平面问题中的应用 6.4.1 上限法原理塑性变形区虚设一个动可容的速度场 ,该速度场应满足以下三个条件:符合位移边界条件;在变形区内保持连续,不产生重叠和拉开;保持体积不变;速度场 对应的应力场 则不一定要求满足平衡条件和力的边界条件。,弹性与塑性力学基础,弹性与塑性力学基础,第六章 塑性力学解题方法及应用举例,6-4 上限法及其在平面问题中的应用 6.4.1 上限法原理动可容速度场 条件下的虚功方程,动可容速度场的变形体,根据最大散逸功原理有凡真实的变形功率总小于在假想速度场情况下所作的功率。这就是上限原理。,

展开阅读全文
相关资源
猜你喜欢
相关搜索
资源标签

当前位置:首页 > 生活休闲 > 在线阅读


备案号:宁ICP备20000045号-2

经营许可证:宁B2-20210002

宁公网安备 64010402000987号