收藏 分销(赏)

钢管砼分体球型节点风电平面塔架的损伤分析.pdf

上传人:自信****多点 文档编号:1467053 上传时间:2024-04-28 格式:PDF 页数:10 大小:4.67MB
下载 相关 举报
钢管砼分体球型节点风电平面塔架的损伤分析.pdf_第1页
第1页 / 共10页
钢管砼分体球型节点风电平面塔架的损伤分析.pdf_第2页
第2页 / 共10页
钢管砼分体球型节点风电平面塔架的损伤分析.pdf_第3页
第3页 / 共10页
亲,该文档总共10页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述

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、斌,程亚超,牛丽华 钢管混凝土风电塔架 型焊接管板节点参数研究 重庆理工大学学报,():,():()闻洋,张军,于洋,等 格构式钢管混凝土风电塔架螺栓球万向节点损伤机制分析 建筑结构学报,(增刊):,():()闻洋,屈琳琳,管丽佩 格构式钢管混凝土风电塔架 型节点受力性能试验研究 建筑结构学报,():,():()闻洋,于蛟,孟春才 格构式钢管砼风力发电塔架包裹球板分支节点屈服机制研究 重庆大学学报,():,():()熊文 格构式钢管混凝土风电塔架万向包裹装配式节点受力性能研究 包头:内蒙古科技大学,:,()韩林海 钢管混凝土结构 理论与实践 版 北京:科学出版社,:,()刘威 钢管混凝土局部受压时的工作机理研究 福州:福州大学,:,()张劲,王庆扬,胡守营,等 混凝土损伤塑性模型参数验证 建筑结构,():,():()第 期闻洋,等:钢管砼分体球型节点风电平面塔架的损伤分析 (),(,):,:,;,:;引用格式:闻洋,徐广茂,李兆建 钢管砼分体球型节点风电平面塔架的损伤分析 西北工业大学学报,():,(),():()(:),

展开阅读全文
部分上传会员的收益排行 01、路***(¥15400+),02、曲****(¥15300+),
03、wei****016(¥13200+),04、大***流(¥12600+),
05、Fis****915(¥4200+),06、h****i(¥4100+),
07、Q**(¥3400+),08、自******点(¥2400+),
09、h*****x(¥1400+),10、c****e(¥1100+),
11、be*****ha(¥800+),12、13********8(¥800+)。
相似文档                                   自信AI助手自信AI助手
百度文库年卡

猜你喜欢                                   自信AI导航自信AI导航
搜索标签

当前位置:首页 > 学术论文 > 论文指导/设计

移动网页_全站_页脚广告1

关于我们      便捷服务       自信AI       AI导航        获赠5币

©2010-2024 宁波自信网络信息技术有限公司  版权所有

客服电话:4008-655-100  投诉/维权电话:4009-655-100

gongan.png浙公网安备33021202000488号   

icp.png浙ICP备2021020529号-1  |  浙B2-20240490  

关注我们 :gzh.png    weibo.png    LOFTER.png 

客服