非线性有限元抗震分析-洞察及研究_第1页
非线性有限元抗震分析-洞察及研究_第2页
非线性有限元抗震分析-洞察及研究_第3页
非线性有限元抗震分析-洞察及研究_第4页
非线性有限元抗震分析-洞察及研究_第5页
已阅读5页,还剩32页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1/1非线性有限元抗震分析第一部分非线性有限元抗震基本原理 2第二部分抗震分析模型构建方法 6第三部分材料非线性特性处理 11第四部分接触与边界条件设置 15第五部分动力响应分析技术 20第六部分抗震性能评价指标 23第七部分算法优化与效率提升 27第八部分案例对比与分析 31

第一部分非线性有限元抗震基本原理

非线性有限元抗震分析是抗震工程领域的一项重要技术,它通过建立结构、地震波和地震作用之间的非线性关系,对结构在地震作用下的动力响应进行精确预测。本文将对非线性有限元抗震基本原理进行介绍。

一、非线性有限元抗震分析的基本概念

非线性有限元抗震分析是指在有限元分析过程中,考虑结构非线性行为的抗震性能分析。它主要针对以下几种非线性现象:

1.材料非线性:包括材料的弹塑性、损伤、断裂等特性。

2.几何非线性:指结构在加载过程中,由于变形引起的几何形状的变化。

3.接触非线性:指结构中节点或单元之间的接触问题。

二、非线性有限元抗震分析的基本原理

1.建立结构模型

首先,根据实际工程情况,建立结构的几何模型和材料模型。在建立模型时,应充分考虑结构的几何尺寸、材料性能、边界条件和支承条件等因素。

2.材料本构关系

根据材料的非线性特性,建立材料本构关系。常见的材料本构关系包括弹塑性、损伤、断裂等模型。在有限元分析中,需要根据实际情况选择合适的本构关系。

3.几何非线性处理

在有限元分析中,几何非线性主要通过以下方法进行处理:

(1)大变形理论:在有限元分析中,采用大变形理论描述结构的变形过程。

(2)几何非线性修正:对单元和节点的位移进行修正,以消除由于几何非线性引起的误差。

4.接触非线性处理

接触非线性主要包括以下处理方法:

(1)接触算法:通过搜索、判断和处理结构中节点或单元之间的接触,实现接触计算。

(2)接触单元:在有限元分析中,引入接触单元来描述节点或单元之间的接触。

5.地震波模拟

地震波模拟是非线性有限元抗震分析的核心部分。常见的地震波模拟方法包括:

(1)时程分析法:根据地震加速度记录,分析结构在不同地震作用下的动力响应。

(2)反应谱分析法:将地震波加速度谱转化为结构动力响应谱,分析结构在不同地震作用下的抗震性能。

6.抗震性能评估

根据非线性有限元抗震分析的结果,对结构的抗震性能进行评估。常用的抗震性能指标包括:

(1)最大位移:结构在地震作用下的最大位移值。

(2)最大加速度:结构在地震作用下的最大加速度值。

(3)破坏概率:结构在地震作用下的破坏概率。

三、非线性有限元抗震分析的应用

非线性有限元抗震分析在抗震工程领域具有广泛的应用,主要包括以下方面:

1.结构抗震设计:通过对结构进行非线性有限元抗震分析,优化结构设计,提高结构的抗震性能。

2.结构抗震加固:针对既有建筑物,通过非线性有限元抗震分析,确定加固方案,提高建筑物的抗震能力。

3.地震灾害评估:利用非线性有限元抗震分析,对地震灾害进行评估,为抗震救灾提供科学依据。

4.新兴结构抗震研究:针对新型结构,通过非线性有限元抗震分析,研究其抗震性能,为新型结构的推广应用提供理论支持。

总之,非线性有限元抗震分析是一种重要的抗震技术,通过对结构非线性行为的分析,为抗震工程领域提供了有力支持。随着计算机技术和数值方法的不断发展,非线性有限元抗震分析将在我国抗震工程领域发挥越来越重要的作用。第二部分抗震分析模型构建方法

非线性有限元抗震分析模型构建方法

摘要:非线性有限元抗震分析是结构抗震设计中的重要手段,其准确性和可靠性直接影响着结构的安全性。本文针对非线性有限元抗震分析模型构建方法进行了深入研究,详细介绍了模型构建的基本步骤、关键技术以及相关参数选取。

1.引言

随着现代建筑技术的不断发展,结构形式日益复杂,其抗震性能也受到广泛关注。非线性有限元抗震分析作为一种有效的计算方法,能够充分考虑结构非线性特性,为结构抗震设计提供有力支持。本文将针对非线性有限元抗震分析模型构建方法进行探讨。

2.模型构建基本步骤

2.1结构几何模型建立

首先,根据实际结构特点,利用CAD软件建立结构几何模型。在建立模型时,应充分考虑结构的几何尺寸、形状以及材料性能等因素。

2.2材料本构关系确定

根据结构材料特性,选取合适的本构模型。常用的本构模型有线性弹性、弹塑性、粘弹性等。在实际工程中,应根据材料性能和结构受力特点,合理选择本构模型。

2.3单元类型及材料属性确定

根据结构形式和受力特点,选取合适的单元类型。常用的单元类型有六节点六面体、八节点六面体、二十节点二十面体等。单元类型确定后,需要赋予单元材料属性,如弹性模量、泊松比、屈服强度等。

2.4荷载与边界条件施加

根据实际工程情况,确定结构所受的各种荷载,如重力、风荷载、地震荷载等。同时,根据结构支撑条件,设置相应的边界条件,如固定端、铰支端等。

2.5非线性求解器选择

非线性有限元抗震分析过程中,需要选择合适的非线性求解器,如Newton-Raphson法、增量法等。选择求解器时,应考虑计算效率、精度等因素。

3.关键技术

3.1网格划分

网格划分是有限元分析中的关键技术之一。合理的网格划分能够提高计算精度和效率。在实际工程中,应遵循以下原则:

(1)网格划分应满足几何和物理条件,保证单元内应力分布均匀;

(2)网格密度应适当,过密或过稀都会影响计算结果;

(3)网格划分应避免出现奇异单元,如三角形单元、退化单元等。

3.2材料非线性处理

材料非线性是影响结构抗震性能的重要因素。在实际计算中,应充分考虑材料非线性特性,如屈服、硬化、软化等。常用的处理方法有:

(1)初始应力法:在分析前先进行一次静力分析,得到材料屈服点的初始应力状态;

(2)增量法:将材料非线性分析过程划分为几个阶段,逐步进行计算;

(3)迭代法:通过迭代求解器,逐步逼近材料非线性响应。

3.3地震波输入

地震波是影响结构抗震性能的重要因素。在实际计算中,应选取合适的地震波,如天然地震波、人工地震波等。地震波的选择应遵循以下原则:

(1)地震波应具有代表性,能够反映实际地震情况;

(2)地震波应具有足够的持续时间,以保证结构在地震作用下的响应能够充分体现;

(3)地震波的峰值加速度应与实际地震情况相吻合。

4.参数选取

4.1材料参数

材料参数是影响分析结果的重要因素。在实际工程中,应根据材料的实际性能和结构受力特点,选取合适的材料参数,如弹性模量、泊松比、屈服强度等。

4.2荷载参数

荷载参数是影响结构抗震性能的关键因素。在实际计算中,应根据工程实际情况,选取合适的荷载参数,如重力、风荷载、地震荷载等。

4.3分析参数

分析参数包括网格划分、求解器选择、计算方法等。在实际计算中,应根据计算效率、精度等因素,合理选取分析参数。

5.结论

非线性有限元抗震分析模型构建方法对于结构抗震设计具有重要意义。本文详细介绍了模型构建的基本步骤、关键技术以及相关参数选取,为非线性有限元抗震分析提供了理论依据。在实际工程中,应根据具体情况进行调整和优化,以提高分析结果的准确性和可靠性。第三部分材料非线性特性处理

非线性有限元抗震分析是结构抗震设计的重要内容,其中材料非线性特性处理是保证分析精度和可靠性的关键环节。本文将针对《非线性有限元抗震分析》中介绍的“材料非线性特性处理”进行详细的阐述。

一、材料非线性特性概述

材料非线性特性是指在结构受力过程中,材料的应力与应变之间存在非线性关系。这是由于材料特性决定的,如弹塑性、粘弹性、各向异性等。对于高烈度地震区及特殊结构,材料非线性特性对结构的抗震性能影响较大,因此,在非线性有限元抗震分析中,必须对材料非线性特性进行处理。

二、材料非线性特性处理方法

1.材料本构方程的建立

建立合适的材料本构方程是处理材料非线性特性的基础。常见的材料本构方程包括:

(1)弹塑性本构方程:描述材料在弹性状态下和塑性状态下的应力与应变关系。常用的弹塑性本构模型有Mises屈服准则、Drucker-Prager屈服准则等。

(2)粘弹性本构方程:描述材料在受力过程中,应力与应变关系随时间变化的特性。常用的粘弹性本构模型有Maxwell模型、Kelvin模型等。

(3)损伤力学本构方程:描述材料在受力过程中,由于损伤而产生力学性能退化的特性。常用的损伤力学本构模型有连续损伤力学、离散损伤力学等。

2.材料非线性特性的有限元离散化

在有限元分析中,将材料非线性特性离散化是保证分析精度的重要环节。常用的离散化方法有:

(1)应力-应变关系离散化:将材料的应力-应变关系离散化为有限单元内的应力-应变关系。

(2)时间离散化:将材料的粘弹性特性离散化为有限时间步长内的应力-应变关系。

(3)空间离散化:将材料的各向异性特性离散化为有限单元内的应力-应变关系。

3.材料非线性特性的迭代求解

在非线性有限元抗震分析中,由于材料非线性特性的存在,导致方程组的非线性,因此需要采用迭代求解方法。常用的迭代求解方法有:

(1)增量法:将非线性问题离散化为一系列线性问题,逐步求解。

(2)全步骤法:直接求解非线性方程组。

(3)Newton-Raphson法:利用泰勒展开,将非线性方程组线性化,然后求解。

4.材料非线性特性的参数识别与验证

在非线性有限元抗震分析中,材料非线性特性的参数识别与验证至关重要。常用的参数识别方法有:

(1)理论分析:根据材料特性,确定材料参数的计算公式。

(2)实验测试:通过实验,确定材料参数的数值。

(3)对比验证:将有限元分析结果与实验结果进行对比,验证参数的准确性。

三、结论

非线性有限元抗震分析中,材料非线性特性处理是保证分析精度和可靠性的关键环节。本文针对《非线性有限元抗震分析》中介绍的“材料非线性特性处理”进行了详细的阐述,包括材料非线性特性概述、处理方法以及参数识别与验证等方面。在实际应用中,应根据结构特点、材料特性和分析要求,选择合适的非线性特性处理方法,以确保非线性有限元抗震分析的准确性和可靠性。第四部分接触与边界条件设置

在非线性有限元抗震分析中,接触与边界条件的设置是至关重要的环节,它直接影响到计算结果的准确性和可靠性。本文将对接触与边界条件的设置进行详细阐述。

一、接触条件的设置

接触条件是有限元分析中的一种特殊边界条件,它描述了两个或多个接触面之间的相互作用。在非线性有限元抗震分析中,接触条件的设置主要包括以下几个方面:

1.接触类型的选择

接触类型的选择是接触条件设置的第一步。常见的接触类型有:

(1)硬接触:当两个接触面之间的距离小于一定的阈值时,它们之间的相互作用为硬接触,即法向和切向力均为无穷大。

(2)软接触:当两个接触面之间的距离小于一定的阈值时,它们之间的相互作用为软接触,即法向和切向力均为有限值。

(3)滑动接触:当两个接触面之间的相对运动超过一定的阈值时,它们之间的相互作用为滑动接触。

选择合适的接触类型对分析结果的准确性具有重要影响。在实际应用中,应根据具体情况选择适合的接触类型。

2.接触表面的处理

接触表面的处理主要包括以下几个方面:

(1)网格划分:接触表面的网格划分应尽可能细密,以提高接触分析的精度。

(2)单元类型:接触表面的单元类型应选择具有良好接触性能的单元,如四面体单元。

(3)接触面的处理:接触面应进行适当的预处理,如去除倒角、倒圆等,以提高接触分析的准确性。

3.接触搜索算法

接触搜索算法是接触条件设置的核心部分,它负责查找并确定接触点的位置。常见的接触搜索算法有:

(1)距离判定法:通过计算两个接触面之间的距离,判断是否存在接触。

(2)投影法:将一个接触面投影到另一个接触面上,判断是否存在接触。

(3)迭代搜索法:逐步缩小搜索范围,直到找到接触点。

4.接触参数的设置

接触参数的设置主要包括以下几个参数:

(1)接触刚度:接触刚度反映了接触面之间的相互作用强度,可通过实验或经验公式确定。

(2)法向刚度:法向刚度反映了接触面之间的法向相互作用强度。

(3)切向刚度:切向刚度反映了接触面之间的切向相互作用强度。

(4)摩擦系数:摩擦系数反映了接触面之间的摩擦作用强度。

二、边界条件的设置

边界条件是有限元分析中的一种特殊约束条件,它描述了结构的边界约束情况。在非线性有限元抗震分析中,边界条件的设置主要包括以下几个方面:

1.边界类型的选择

边界类型的选择主要包括固定边界、自由边界和滑动边界。固定边界限制了结构的位移和转动;自由边界允许结构的位移和转动;滑动边界限制了结构的位移,但不限制其转动。

2.边界约束的施加

边界约束的施加主要包括以下几种方式:

(1)位移约束:限制结构的位移,如固定边界。

(2)转动约束:限制结构的转动,如固定边界。

(3)法向约束:限制结构的法向位移,如固定边界。

(4)切向约束:限制结构的切向位移,如固定边界。

3.边界参数的设置

边界参数的设置主要包括以下几种参数:

(1)固定边界的位移:固定边界的位移值应与实际情况一致。

(2)固定边界的转动:固定边界的转动值应与实际情况一致。

(3)法向约束的刚度:法向约束的刚度应与实际情况一致。

(4)切向约束的刚度:切向约束的刚度应与实际情况一致。

综上所述,在非线性有限元抗震分析中,接触与边界条件的设置对分析结果的准确性具有重要影响。在实际应用中,应根据具体情况选择合适的接触类型、接触表面处理、接触搜索算法和接触参数,以及合适的边界类型、边界约束和边界参数,以确保分析结果的可靠性和准确性。第五部分动力响应分析技术

《非线性有限元抗震分析》中的动力响应分析技术

动力响应分析技术在结构工程领域具有重要意义,尤其是在进行非线性有限元抗震分析时。该技术旨在模拟结构在地震作用下的动态响应,以评估结构的抗震性能和安全性。以下是对动力响应分析技术在非线性有限元抗震分析中的详细介绍。

一、动力响应分析的基本原理

动力响应分析基于牛顿第二定律和结构动力学的基本原理。在地震作用下,结构受到周期性激励,导致结构产生加速度、速度和位移等动态响应。动力响应分析的目标是预测结构在地震作用下的动态行为,包括位移、速度、加速度和内力等。

二、动力响应分析的主要步骤

1.建立结构模型:首先,根据实际结构的特点,建立相应的有限元模型。模型应考虑结构的几何特征、材料性能和边界条件等因素。

2.材料非线性分析:由于地震作用下结构的材料可能进入非线性状态,因此需要对材料进行非线性分析。常用的材料本构模型有弹性模型、弹塑性模型和粘弹性模型等。

3.地震输入:地震动是动力响应分析中的重要输入条件。根据实际地震记录或地震动参数,确定地震动的时程曲线。

4.计算动力响应:采用数值方法,如Newmark-β法、Wilson-θ法等,求解结构在地震作用下的动力响应方程。计算过程中,需考虑阻尼、质量、刚度等因素。

5.验证分析结果:通过对比实际地震记录或相关规范要求,对计算结果进行验证。若分析结果满足规范要求,则认为结构具有较好的抗震性能。

三、动力响应分析技术在非线性有限元抗震分析中的应用

1.结构抗震性能评估:动力响应分析可以评估结构在地震作用下的抗震性能,包括位移、速度、加速度和内力等。通过对不同地震动输入和参数的敏感性分析,可以确定结构的关键薄弱环节。

2.结构优化设计:根据动力响应分析结果,可以对结构进行优化设计。通过调整结构参数(如截面尺寸、材料性能等),提高结构的抗震性能。

3.结构加固设计:对于已建结构,动力响应分析可以确定加固部位和加固措施。通过加固,提高结构的抗震能力。

4.震后评估:在地震发生后,动力响应分析可以用于评估结构的震后状态,为灾后重建提供依据。

四、动力响应分析技术的展望

随着计算技术的不断发展,动力响应分析技术在非线性有限元抗震分析中的应用将更加广泛。以下是对动力响应分析技术未来发展的展望:

1.高性能计算:提高计算效率,实现大规模结构动力响应分析。

2.新型材料分析:研究新型材料的本构模型,提高动力响应分析的准确性。

3.智能化分析:结合人工智能技术,实现动力响应分析过程的智能化。

4.跨学科研究:加强结构工程、地震工程和计算力学等学科的交叉研究,推动动力响应分析技术的发展。

总之,动力响应分析技术在非线性有限元抗震分析中具有重要作用。通过不断完善分析方法和理论,提高分析精度,为结构抗震设计和安全评估提供有力支持。第六部分抗震性能评价指标

非线性有限元抗震分析中的抗震性能评价指标

摘要:非线性有限元抗震分析作为结构抗震设计的核心技术之一,在工程实践中具有重要作用。本文针对非线性有限元抗震分析中的抗震性能评价指标进行综述,旨在为相关研究提供参考。

一、引言

随着我国城市化进程的加快,建筑工程规模不断扩大,结构形式也越来越复杂。为确保建筑物在地震作用下的安全性和可靠性,非线性有限元抗震分析成为结构抗震设计的重要手段。在非线性有限元抗震分析过程中,抗震性能评价指标的选择与确定至关重要。本文对非线性有限元抗震分析中的抗震性能评价指标进行综述,包括位移、屈服、破坏以及反应谱等指标。

二、位移指标

位移指标是衡量结构在地震作用下变形程度的重要指标,主要包括以下几种:

1.节点位移:节点位移指的是结构中各个节点的位移值,通常以最大位移值作为评价指标。节点位移反映了结构在地震作用下的整体变形情况,对评估结构的抗震性能具有重要意义。

2.屈曲位移:屈曲位移是指结构在地震作用下发生屈曲时的变形量,通常以屈曲位移与结构尺寸的比值作为评价指标。屈曲位移反映了结构在地震作用下的局部变形情况,对评估结构的整体抗震性能具有重要作用。

3.极限位移:极限位移是指结构在地震作用下达到极限状态时的变形量,通常以极限位移与结构尺寸的比值作为评价指标。极限位移反映了结构的最大变形能力,对评估结构的抗震性能具有重要意义。

三、屈服指标

屈服指标是衡量结构在地震作用下屈服程度的重要指标,主要包括以下几种:

1.屈服位移:屈服位移是指结构在地震作用下发生屈服时的位移值,通常以屈服位移与结构尺寸的比值作为评价指标。屈服位移反映了结构的屈服能力,对评估结构的抗震性能具有重要意义。

2.屈服弯矩:屈服弯矩是指结构在地震作用下发生屈服时的弯矩值,通常以屈服弯矩与结构尺寸的比值作为评价指标。屈服弯矩反映了结构的抗弯能力,对评估结构的抗震性能具有重要意义。

3.屈服剪力:屈服剪力是指结构在地震作用下发生屈服时的剪力值,通常以屈服剪力与结构尺寸的比值作为评价指标。屈服剪力反映了结构的抗剪能力,对评估结构的抗震性能具有重要意义。

四、破坏指标

破坏指标是衡量结构在地震作用下破坏程度的重要指标,主要包括以下几种:

1.部分破坏:部分破坏是指结构在地震作用下部分构件发生破坏,但整体结构仍能保持一定的承载能力。部分破坏通常以破坏构件数量占比作为评价指标。

2.全部破坏:全部破坏是指结构在地震作用下所有构件均发生破坏,结构失去承载能力。全部破坏通常以结构破坏状态作为评价指标。

3.破坏极限荷载:破坏极限荷载是指结构在地震作用下达到破坏状态时的荷载值,通常以破坏极限荷载与结构设计荷载的比值作为评价指标。

五、反应谱指标

反应谱指标是衡量结构在地震作用下动力响应特性的重要指标,主要包括以下几种:

1.基本周期:基本周期是指结构在地震作用下的振动周期,通常以基本周期与结构自振周期的比值作为评价指标。

2.响应谱峰值:响应谱峰值是指结构在地震作用下的最大动力响应值,通常以响应谱峰值为评价指标。

3.频率响应:频率响应是指结构在不同频率地震作用下的动力响应,通常以频率响应与结构自振频率的比值作为评价指标。

六、结论

非线性有限元抗震分析中的抗震性能评价指标主要包括位移、屈服、破坏以及反应谱等指标。在实际工程中,应根据具体情况进行指标选择与确定,以全面评估结构的抗震性能。本文对非线性有限元抗震分析中的抗震性能评价指标进行了综述,为相关研究提供了参考。第七部分算法优化与效率提升

在《非线性有限元抗震分析》一文中,算法优化与效率提升是提高非线性有限元分析准确性和计算效率的重要方面。以下是对该部分内容的简明扼要介绍:

一、算法优化

1.前处理算法优化

(1)网格划分优化:采用自适应网格划分技术,根据应力、应变等物理量的变化,动态调整网格密度,提高计算精度。

(2)单元类型选择:针对不同结构特点,选择合适的单元类型,如三角形、四边形、六面体等,以适应复杂几何形状。

(3)材料属性处理:采用分段线性或非线性插值方法,将复杂材料属性离散化,提高计算效率。

2.后处理算法优化

(1)结果提取优化:针对不同需求,提取关键节点、线段、面的位移、应力、应变等物理量,提高数据提取效率。

(2)结果可视化优化:采用高性能渲染技术,实现结果的高清展示和动画演示,便于分析者直观理解分析结果。

3.迭代算法优化

(1)收敛性提高:通过改进迭代算法,如高斯-赛德尔法、共轭梯度法等,提高算法的收敛速度。

(2)计算精度优化:采用预处理器、后处理器等技术,提高计算精度。

二、效率提升

1.并行计算

(1)多核处理器利用:针对多核处理器,采用并行算法,如多线程、GPU计算等,实现计算任务并行化。

(2)分布式计算:利用高性能计算集群,通过分布式计算技术,实现大规模问题的快速求解。

2.优化计算方法

(1)预处理器:通过预处理器技术,如参数化设计、子结构分解等,减少计算量和计算时间。

(2)后处理器:通过后处理器技术,如结果提取、可视化等,提高数据处理效率。

3.模型简化与近似

(1)模型简化:针对实际工程问题,对模型进行简化,如去除冗余节点、单元等,降低计算量。

(2)近似计算:采用有限元近似方法,如数值积分近似、边界元法等,提高计算效率。

4.优化算法参数

(1)算法参数调整:针对不同问题,调整算法参数,如迭代次数、收敛精度等,提高计算效率。

(2)计算资源分配:根据计算任务的特点,合理分配计算资源,提高计算效率。

综上所述,非线性有限元抗震分析中的算法优化与效率提升是提高计算精度和效率的关键。通过前处理、后处理、迭代算法等方面的优化,以及并行计算、优化计算方法、模型简化与近似、优化算法参数等方面的改进,可以有效提高非线性有限元分析的精确性和计算效率。第八部分案例对比与分析

在文章《非线性有限元抗震分析》中,案例对比与分析部分主要针对不同抗震性能的建筑结构进行了深入探讨。以下是

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论