ABAQUS钢框架结构抗震仿真分析之欧阳家百创编
- 格式:doc
- 大小:105.51 KB
- 文档页数:14
ABAQUS钢框架结构抗震仿真分析D才建好的ac加速度,OK。
7、MeshSeed Part Instance;Assign Mesh Controls:Quard,Free,Advancing front;Assign Element Type:Standard,Linear,Reduced Integration,壳单元S4R;Mesh Part;最后Verify Mesh检查网格是否错误。
8、JobCreate Job;Write Input(输出inp文件),Submit;打开Monitor观察是否有错误、警报;点击Result进入可视化窗口Visualization。
9、VisualizationResult,Step/Frame分析步/帧,查看各时程变形情况;Tools,XY Data,Create,ODB field output,plot顶层楼板位移时程曲线,各层楼板加速度时程曲线,各层楼板位移曲线。
10、数据后处理在Visualization,XY Data Manager中点击Edit,输出各层楼板位移数据到Excel中,利用Excel:相邻两层的位移作差后得到层间相对位移,画出以Time为横轴、层间相对位移为纵轴画出各层层间位移时程曲线;然后层间相对位移除以层高得到层间位移角,画出以Time为横轴、层间位移角为纵轴画出各层层间位移角时程曲线。
四、结果展示钢框架内力云图1钢框架内力云图2钢框架内力云图3顶层楼板位移时程曲线底端加速度时程曲线一层楼板加速度时程曲线二层楼板加速度时程曲线三层楼板加速度时程曲线四层楼板加速度时程曲线五层楼板加速度时程曲线顶层楼板加速度时程曲线底端与一层楼板层间位移一层与二层楼板层间位移二层与三层楼板层间位移三层与四层楼板层间位移四层与五层楼板层间位移五层与顶层楼板层间位移底端与一层楼板层间位移角一层与二层楼板层间位移角二层与三层楼板层间位移角三层与四层楼板层间位移角四层与五层楼板层间位移角五层与六层楼板层间位移角五、总结与分析通过分析数据:①该钢框架模型顶层楼板最大位移发生在9.904615s时,值为0.521279;根据《钢结构设计规范》:多层框架结构柱顶位移H/500,对于该模型为18m/500=0.036m,顶层楼板最大位移超出规范要求。
ABAQUS时程分析法计算地震反应的简单实例地震反应分析是工程结构设计中非常重要的一项内容,ABAQUS是一款常用的有限元分析软件,可以用于进行地震反应分析。
下面将通过一个简单的实例来介绍ABAQUS中的时程分析法计算地震反应的方法。
假设我们要计算一个两层楼的简化结构在地震作用下的反应。
该结构由两根弹性细杆组成,分别代表两层楼板的水平支撑,两层楼板分别由质点模拟。
结构的初始状态为静力平衡,无任何位移和速度,没有外力作用。
我们将利用ABAQUS中的时程分析法来计算结构在地震作用下的运动。
1.组建模型首先,我们需要在ABAQUS中组建结构的有限元模型。
在ABAQUS中,我们可以使用CONNECTION命令来定义模型的节点和单元。
在本例中,我们定义了两个节点表示两层楼板的连接处,再使用ELEMENT命令定义两根弹性细杆,将两个节点连接起来。
这样就创建了一个简化的两层楼模型。
2.定义材料和截面属性在ABAQUS中,我们还需要定义模型的材料和截面属性。
在本例中,楼板部分可以使用线性弹性材料进行模拟,我们可以使用ELASTIC命令定义材料的弹性模量和泊松比。
而弹性细杆可以使用线性弹性的截面,我们可以使用SECTION命令定义截面的几何特性和材料特性。
3.定义边界条件为了模拟地震作用下的结构反应,我们需要定义结构的边界条件。
在本例中,我们希望模型底部节点固定,即节点的位移和旋转均为零。
我们可以使用BOUNDARY命令来定义节点的边界条件。
4.定义地震荷载在ABAQUS中,我们可以使用ACCELERATION命令来定义地震的加速度时程。
地震时程是描述地震波动过程的函数,通常由科学家通过实测数据进行获取。
在本例中,假设我们已经获取到地震波动的加速度时程,并将其转换为ABAQUS可用的格式。
5.设置时程分析参数在进行时程分析之前,我们还需要设置模型的时程分析参数。
在ABAQUS中,我们可以使用STEP命令设置分析的步数、步长和加载参数。
基于ABAQUS的框架结构抗震分析朱秀亭;于广杰;杨晨辉【摘要】介绍了应用ABAQUS对框架结构进行时程分析法和反应谱分析法的具体的建模及分析要点,通过对框架结构的应力云图、位移云图和加速度云图的分析,以及后处理的数据分析结果可以看出,运用ABAQUS分析抗震分析的正确性及可靠性都是可以保证的,为研究人员运用ABAQUS进行更为复杂的分线性分析工作奠定了基础.【期刊名称】《山西建筑》【年(卷),期】2013(039)011【总页数】3页(P36-38)【关键词】时程分析法;地震波;反应谱【作者】朱秀亭;于广杰;杨晨辉【作者单位】扬中市市政园林工程处,江苏扬中212200【正文语种】中文【中图分类】TU375.40 引言目前ABAQUS由于其强大的计算器功能已成为极为流行的有限元分析工具,现已广泛的应用于机械、化工、土木、航空及船舶等各个工程和科研领域[1-3]。
地震由于对人们的生命财产造成了极其严重的破坏,对其的研究分析已成为重要的课题[4-7]。
本文将介绍用ABAQUS进行建筑抗震分析的方法,以作为从事建筑工程设计及学习人员的学习参考。
1 抗震时程分析时程分析法是对结构动力方程直接进行逐步积分求解的一种动力分析方法。
时程分析法将地震波按时段进行数值化后,输入结构体系的振动微分方程,采用直接积分法计算出结构在整个强震时域中的振动状态全过程,给出各个时刻各个杆件的内力和变形。
1.1 数值分析模型本文的模型为一榀钢结构框架结构,框架柱子为300 mm×300 mm的矩形钢柱,框架梁为工字形钢梁,模型的具体尺寸如图1所示。
钢柱采用二级钢HRB335,弹性模量为Es=2.1×105MPa,钢梁采用一级钢HPB300,弹性模量为Es=2.0×105MPa。
模型采用Beam单元,并进行B31单元离散。
ABAQUS进行时程分析时选择DYNAMIC IMLICITY进行分析。
具体的数值分析模型如图2所示。
ABAQUS钢框架结构抗震仿真分析首先,我们需要建立结构的有限元模型。
钢框架结构主要由柱、梁、节点和连接件组成,我们需要根据实际情况进行建模。
在ABAQUS中,我们可以使用节点(节点)和单元(单元)建立结构模型。
其次,我们需要定义结构的材料特性。
在钢框架结构中,材料的弹性模量(E)和泊松比(ν)是两个重要参数。
根据实际材料的特性,我们可以在ABAQUS中定义这些参数。
接下来,我们需要定义结构的边界条件。
抗震仿真分析通常需要在地震力作用下进行,我们需要定义结构的固定支撑条件,以模拟垂直方向上的地震力。
在ABAQUS中,我们可以将结构的底部或其他特定地方固定支撑。
然后,我们需要定义地震载荷。
地震力通常由地震加速度谱表示,在ABAQUS中,我们可以通过载荷定义来输入这些数据。
根据地震保护设计准则,我们可以计算出地震力对结构的作用。
在进行抗震仿真分析之前,我们还需要进行网格划分和网格优化。
钢框架结构通常具有较高的刚度和复杂的形状,我们需要根据结构的实际情况进行网格划分,并使用ABAQUS的网格优化工具来确保网格质量。
最后,我们可以进行抗震仿真分析。
在此过程中,我们可以将地震载荷应用于结构,并模拟结构在地震力作用下的响应。
ABAQUS可以计算出结构的位移、应力和变形等参数,并可生成相应的结果报告。
总结起来,ABAQUS是一种强大的有限元分析工具,可以用于钢框架结构的抗震仿真分析。
通过建立模型、定义材料特性、边界条件和地震载荷,进行网格划分和网格优化,并进行仿真分析,我们可以获取结构在地震力作用下的响应情况,评估结构的抗震性能,并指导实际工程设计。
前言有限元法是工程中广泛使用的一种数值计算方法。
它是力学、计算方法和计算机技术相结合的产物。
在工程应用中,有限元法比其它数值分析方法更流行的一个重要原因在于:相对与其它数值分析方法,有限元法对边界的模拟更灵活,近似程度更高。
所以,伴随着有限元理论以及计算机技术的发展,大有限元软件的应用证变得越来越普及。
ABAQUS软件一直以非线性有限元分析软件而闻名,这也是它和ANSYS,Nastran等软件的区别所在。
非线性有限元分析的用处越来越大,因为在所用材料非常复杂很多情况下,用线性分析来近似已不再有效。
比方说,一个复合材料就不能用传统的线性分析软件包进行分析。
任何与时间有关联,有较大位移量的情况都不能用线性分析法来处理。
多年前,虽然非线性分析能更适合、更准确的处理问题,但是由于当时计算设备的能力不够强大、非线性分析软件包线性分析功能不够健全,所以通常采用线性处理的方法。
这种情况已经得到了极大的改善,计算设备的能力变得更加强大、类似ABAQUS这样的产品功能日臻完善,应用日益广泛。
非线性有限元分析在各个制造行业得到了广泛应用,有不少大型用户。
航空航天业一直是非线性有限元分析的大客户,一个重要原因是大量使用复合材料。
新一代波音 787客机将全部采用复合材料。
只有像 ABAQUS这样的软件,才能分析包括多个子系统的产品耐久性能。
在汽车业,用线性有限元分析来做四轮耐久性分析不可能得到足够准确的结果。
分析汽车的整体和各个子系统的性能要求(如悬挂系统等)需要进行非线性分析。
在土木工程业, ABAQUS能处理包括混凝土静动力开裂分析以及沥青混凝土方面的静动力分析,还能处理高度复杂非线性材料的损伤和断裂问题,这对于大型桥梁结构,高层建筑的结构分析非常有效。
瞬态、大变形、高级材料的碰撞问题必须用非线性有限元分析来计算。
线性分析在这种情况下是不适用的。
以往有一些专门的软件来分析碰撞问题,但现在ABAQUS在通用有限元软件包就能解决这些问题。
基于BQUS下钢筋水泥房屋结构抗震能力分析徐绍娟李良江谢圣富张敏谢顺添摘要:本文通过对XX地区钢筋混泥土房屋结构的调查,展开了钢筋混泥土房屋结构抗震特点的讨论。
基于bqus有限元分析软件建立模型,然后接受模态分析计算出该模型的前四阶振型及固有频率,并与公式的计算频率进行对比分析,预报该模型的还原程度及可信度。
选取XX省地震波为地震动的输入时程,分别计算在不同峰值加速度时程下的钢筋混泥土结构的地震反应,得到相应的拉伸损伤XX图,并分析机构的破坏特征,给出该模型抗震能力薄弱之处。
结果说明:当钢筋混泥土房屋结构输入东西和南北地震动时,首先由窗户上横墙和纵墙相交汇的地方开裂,且两相邻交汇处裂纹渐渐靠拢,又由于整体重力过大导致该房屋结构的下半部分领先损坏,直至房屋倒塌。
关键词:钢筋混泥土房屋结构;抗震特点 bqus有限元分析;模态分析;频率;拉伸损伤XX图引言本文从XX地区现有的农居结构出发,选出具有代表性的建筑作为讨论对象,利用有限元软件bqus建模,并输入不同的地震波,探究不同的地震等级对房屋结构的影响。
1 模型建立本文针对XX地区进行了抽样调查。
最终接受bqus建立调查后的钢筋混泥土房屋结构模型。
1.1 调查分析依据调查状况,本文选用了XX常见钢筋混泥土房屋为主要的参考对象。
通过对其结构尺寸进行提取,可得到本文需建立的模型基本参数如表所示:在模型的建立中依据上述基本尺寸进行建立,且房屋的样式与主要参考对象略有不同。
1.2 有限元模型的建立为完成钢筋混泥土结构的抗震性能分析。
本文首先使用bqus 有限元分析软件中的部件〔prt〕以及互相作用关系〔interction〕模塊完成了基础模型的建立。
其次,在模型进行有限元分析之前还还需要对模型的材料以及XX格进行划分,这里则需要使用的是bqus中的属性〔property〕和XX格〔Mesh〕模块,且在设置属性的过程中需要对钢筋混泥土材料的本构关系进行确定,而XX格〔Mesh〕的作用则是将模型划分为小的单元,以便分析。
第18卷第6期2020年12月水利与建筑工程学报JournalofWaterResourcesandArchitecturalEngineeringVol.18No.6Dec.,2020DOI:10.3969/j.issn.1672-1144.2020.06.025收稿日期:2020 08 09 修稿日期:2020 09 01基金项目:国家自然科学基金青年项目(51408223);国家自然科学基金面上项目(51679091;51979109)作者简介:程学斌(1995—),男,江西上饶人,硕士研究生,研究方向为工程结构抗震。
E mail:lbjcheng@163.com通讯作者:马 颖(1982—),女,河南郑州人,博士,硕士生导师,主要从事水工、桥梁等工程结构抗震研究工作。
E mail:maying198208@163.com基于ABAQUS的钢筋混凝土柱抗震数值模拟分析程学斌,马 颖,袁子淇(华北水利水电大学水利学院,河南郑州450045)摘 要:为了研究ABAQUS软件中实体单元和纤维梁单元在不同破坏模式下(RC)柱滞回性能数值模拟的适用性,从美国PEER数据库中收集了9根钢筋混凝土矩形截面柱的拟静力试验数据,柱试件分别发生了弯曲、弯剪或剪切不同模式的破坏。
在ABAQUS中分别建立柱试件的实体单元模型和纤维梁单元模型并进行往复荷载作用下RC柱滞回性能的数值模拟,将模拟结果与试验数据进行了对比分析。
结果表明:对于弯曲破坏RC柱,适合采用纤维梁单元模拟,而对于弯剪破坏和剪切破坏RC柱,基于实体单元的模拟结果与试验结果更为接近;纤维梁单元能够更准确地模拟RC柱滞回曲线的捏拢效应。
关键词:钢筋混凝土柱;ABAQUS;实体单元;纤维梁单元;滞回性能中图分类号:TU375.3 文献标识码:A 文章编号:1672—1144(2020)06—0146—07SeismicNumericalSimulationAnalysisofReinforcedConcreteColumnsBasedonABAQUSCHENGXuebin,MAYing,YUANZiqi(SchoolofWaterConservancy,NorthChinaUniversityofWaterResourcesandElectricPower,Zhengzhou,He’nan450045,China)Abstract:InordertoassesstheapplicabilityofnumericalsimulationofhystereticbehaviorofRCcolumnswithsolidelementandfiberbeamelementinABAQUSsoftwareunderdifferentfailuremodes.Thepseudo statictestdataof9re inforcedconcreterectangularcross sectioncolumnswerecollectedfromthePEERdatabaseintheUnitedStates.Thecolumnspecimenswerefailedindifferentmodesofflexure,flexure shearorshear.Basedontheforcetestdata,thesolidelementmodelandthefiberbeamelementmodelofthespecimenwereestablishedinABAQUStosimulatethehystereticperformanceoftheRCcolumnunderthereciprocatingload.Thesimulationresultswerecomparedwiththetestdata.TheresultsshowthatforflexurefailureRCcolumns,fiberbeamelementsimulationissuitable,whileforflexure shearfailureandshearfailureRCcolumns,thesimulationresultsbasedonsolidelementsareclosertothetestresults,andfiberbeamelementcanmoreaccuratelysimulatethepincheffectofRCcolumnhystereticcurve.Keywords:reinforcedconcretecolumns;ABAQUS;solidelement;fiberbeamelement;hystereticperformance 在地震作用下,钢筋混凝土柱作为水工、桥梁、房屋等结构的主要竖向承重与水平抗力构件,承载着整个结构的竖向荷载和由地震引起的水平荷载。
ABAQUS钢框架结构抗震仿真分析报告【摘要】本文使用ABAQUS软件对一幢钢框架结构进行了抗震仿真分析。
首先,建立了结构的有限元模型,并对其进行了网格划分。
然后,加载了地震波荷载,进行了静力分析和动力分析。
最后,通过比较结构的位移响应和内力分布,评估了结构的抗震性能。
结果表明,结构具有较好的抗震能力。
【关键词】ABAQUS;钢框架结构;抗震仿真;有限元模型;地震波荷载1.引言随着城市化进程的推进,建设规模越来越大的钢框架结构变得越来越普遍。
然而,地震是一个常见的自然灾害,在一些地区频繁发生。
因此,钢框架结构的抗震性能成为了设计的重要考虑因素。
为了评估钢框架结构的抗震性能,可以通过抗震仿真分析来模拟地震情况,得到结构的位移、应力等响应。
2.方法2.1建立有限元模型首先,根据结构的几何形状和材料性质,建立了合适的有限元模型。
钢框架结构主要由梁柱组成,因此可以使用梁单元和柱单元来建模。
在建立模型时,需要考虑结构的几何非线性和材料非线性。
2.2网格划分在建立有限元模型后,需要对结构进行网格划分。
合理的网格划分能够提高计算精度和计算效率。
一般来说,细小的单元可以更好地模拟结构的性能,但也会增加计算量,因此需要权衡。
3.分析3.1静力分析首先,按照建筑物受到的地震荷载大小进行静力分析。
静力分析是为了确定结构在地震荷载下的受力状态。
通过静力分析,可以获得结构的位移响应和内力分布。
3.2动力分析在静力分析的基础上,进行动力分析。
动力分析是为了模拟地震时结构的动态响应。
在动力分析中,需要加载地震波荷载,并设置一定的计算时间。
通过动力分析,可以获得结构在地震中的动态位移响应和内力分布。
4.结果与讨论通过比较静力分析和动力分析的结果,可以评估钢框架结构的抗震性能。
如果位移响应较小,内力分布均匀,说明结构具有较好的抗震能力。
反之,则说明结构抗震能力较差。
5.结论本文使用ABAQUS软件对一幢钢框架结构进行了抗震仿真分析。
一、引言欧阳引擎(2021.01.01)时程分析法是对结构动力方程直接进行逐步积分求解的一种动力分析方法。
时程分析法将地震波按时段进行数值化后,输入结构体系的振动微分方程,采用直接积分法计算出结构在整个强震时域中的振动状态全过程,给出各个时刻各个杆件的内力和变形。
现已成为多数国家抗震设计规范或规程的分析方法之一。
二、有限元软件ABAQUS简介ABAQUS是美国ABAQUS公司(原名HKS公司.即Hibbitt,Karlsson&Sorensen,INC.2005年被法国达索公司收购,2007年公司更名为SIMULIA)。
ABAQUS已成为国际上最先进的大型通用有限元力学分析软件之一。
ABAQUS是一套功能强大的进行工程模拟的有限元软件。
其解决问题的范围从相对简单的线性分析到许多复杂的非线性问题。
ABAQUS拥有CAE工业领域最为广泛的材料模型,它可以模拟绝大部分工程材料的线形和非线形行为,可以进行结构的静态和动态分析,如应力、变形、振动、热传导以及对流等。
也可以模拟广泛的材料性能,如金属、橡胶、塑料、弹性泡沫等,而且任何一种材料都可以和任何一种单元或复合材料的层一起用于任何合适的分析类型。
三、模型建立与求解1、PartCreate Part:Name:Ban,3D,Deformable,Shell ,Planar,输入坐标创建一个18X9m的壳部件,作为混凝土楼板部件;Create Part:Name:Zhu,3D,Deformable,Wire ,Planar,输入坐标创建一个长3m线部件,作为柱部件;Create Part:Name:Liang,3D,Deformable,Wire ,Planar,输入坐标创建一个长6X3m,宽4.5X2m的线网部件,作为梁网部件;2、 SectionCreate Material:Name:steel,General,Density 7800;Mechanical,Elasticity,Young’s Modulus 2.1e11,Poisson’Ratio 0.3;Create Material:Name:concrete,General,Density 2500;Mechanical,Elasticity,Young’s Modulus 3e10,Poisson’Ratio 0.3。
基于ABAQUS的楼板对钢筋混凝土框架的抗震影响分析钢筋混凝土框架是一种常用于建筑结构中的结构形式,其在抗震设计中起到了至关重要的作用。
然而,由于楼板的存在,框架结构的抗震性能可能会受到一定的影响。
因此,通过使用ABAQUS软件对楼板的抗震影响进行分析,能够更好地评估钢筋混凝土框架结构的抗震性能。
首先,为了进行抗震分析,需要使用ABAQUS软件进行有限元建模。
在建模过程中,需要考虑框架结构的几何形状、材料性质和边界条件等因素。
通过细分建模,可以更精确地模拟结构的受力状态。
同时,应根据实际情况合理设置初始条件,如初始位移和初始速度等。
在建模完成后,需要进行地震荷载的施加。
ABAQUS软件可以通过施加动力荷载或地震时间历程来模拟地震作用。
在分析过程中,可以通过定义响应谱等指标来评估结构的抗震性能。
此外,还可以通过初始刚度矩阵的求解,确定结构的初始刚度分布,为后续的动力分析提供准确的初始条件。
在模拟地震作用后,需要进行动力分析以评估结构的受力性能。
通过ABAQUS软件的求解器,可以得到结构的响应结果,如位移、加速度和内力等。
通过对这些结果的分析,可以了解结构在地震作用下的受力情况,对结构的稳定性和安全性进行评估。
此外,需要注意的是,在进行抗震分析时,还需要考虑楼板对结构的影响。
由于楼板的存在,可能会产生一定的刚度和阻尼效应,对结构的抗震性能产生一定的影响。
因此,在建模过程中,需要合理设置楼板的材料属性和几何参数,以准确模拟楼板对结构的影响。
综上所述,基于ABAQUS的楼板对钢筋混凝土框架的抗震影响分析是一种重要的抗震设计方法。
通过细致的建模和准确的分析,可以更好地评估钢筋混凝土框架结构的抗震性能,为工程设计提供依据。
在实际应用中,还可以进一步优化模型和参数,提高分析结果的准确性和可靠性,以确保工程的安全性和稳定性。
ABAQUS在抗震性能分析中的应用作者:郭开发来源:《大东方》2018年第09期摘要:空间框架结构整体性、刚度较好,仍广泛应用于我国大部分地区。
但即使严格执行抗震设计规范的空间框架结构也可能发生不同程度的倒塌破坏,因此其抗震性能是保障人民生命财产安全的重要指标。
本文主要简述ABAQUS有限元分析软件对地震荷载下的空间框架结构抗震性能分析的优势及意义。
关键词:抗震理论;abaqus分析优势;意义一、抗震理论的发展早期抗震工程研究主要侧重于如何减少巨大地震下的建筑物倒塌和人员伤亡,基于承载力和构造保证延性的传统抗震设计方法并不完善,已不能适应现代社会对结构抗震性能的要求。
由此引起了各国工程界对现有抗震设计思想和方法进行深刻的反思,迫使工程人员寻求更加完善的设计思想,设工程结构在各种可能遇到的地震作用下的反应和损伤状态控制在设计预期要求的范围内。
不仅确保生命安全,而且确保经济损失最少。
我国抗震规范(GB50011-2010)采用两阶段的抗震设计思想:第一阶段,对于小震和中震作用下的一般结构采用弹性分析方法计算构件截面及配筋;第二阶段,对于特殊结构(如超限高层结构)需使用弹塑性分析方法进行罕遇地震作用下结构的损伤和抗倒塌能力验算。
目前,结构弹塑性分析的方法有两种,分别为静力弹塑性分析和动力弹塑性分析。
静力弹塑性分析又称静力推覆分析或Pushover分析,它是一种简便易行,能够在一定程度上反映结构弹塑性性能的方法。
动力弹塑性分析是利用结构动力学方程和构件弹塑性恢复力模型,采用逐步积分技术获得各时刻下结构各质点和构件地震响应的方法。
其中静力弹塑性分析理论得到了较多的完善和发展,虽然较动力弹塑性理论有就较多的缺陷,但它已成为了一种完善的方法。
而使用动力弹塑性理论分析具有以下优势:(1)可以获得整个地震过程中各时刻结构的地震响应;(2)对结构的简化假设较少,可以完整、真实、准确地反应结构地震响应。
无论是静力弹塑性理论还是动力弹塑性理论,应用有限元分析更为简便,得到的结果也更能接近结构受荷时的真实内力分布状况。
(一)创建部件1:模块:部件2:从菜单栏中选择部件→创建,弹出创建部件对话框,将名称修改为Basis,模型空间选择三维,类型选择可变形,形状选择实体,类型为拉伸。
采用SI (mm)量纲,考虑到所创建部件的最大尺寸为100mm,大约尺寸改为200,为最大尺寸的2倍。
3:点击继续,进入草绘模式,为实体拉伸绘制界面草图。
4:点击创建构造工具,创建2条构造线,分别与X 轴、Y轴重合5:点击添加约束工具,弹出添加约束对话框,选择固定约束,按住shift点选2条构造线,按鼠标中键。
6:点击创建矩形工具,在绘图区内随意绘制一个矩形7:点击添加尺寸工具,修改矩形的长为100,宽为608:点击添加约束工具,使矩形分别关于2条构造线对称9:点击创建圆工具,先绘制1个圆,半径为5,圆心到长边的距离为15,圆心到宽边的距离为1010:点击镜像工具,选择镜像操作的类型为复制,选择水平方向的构造线为镜像轴,选择上一步绘制的第1个圆为待镜像的对象,按下鼠标中键完成镜像11:使用同样的方法绘制另外2个圆12:按下鼠标中键完成截面草图的绘制,弹出编辑基本拉伸对话框,将深度修改为20,点击确定13:拉伸出的实体如下图14:点击创建实体:拉伸工具,为实体拉伸选择一个合适的平面,点选一条合适的边作为草绘的参照,进入草绘模式15:绘制如下界面草图16:按下鼠标中键退出草绘模式,弹出编辑拉伸对话框,类型为指定深度,深度设为20,由于绘图区中的拉伸方向不是想要的方向,点击翻转方向按钮17:点击确定,完成拉伸18:点击创建切削:拉伸工具,为拉伸切削选择一个合适的平面,选择一条合适的边作为草绘参照19:绘制如下界面草图:20:按下鼠标中键完成草图绘制,弹出编辑切削拉伸对话框,类型选择指定深度,深度修改为5,选择默认的拉伸切削方向,点击确定,结束切削拉伸21:点击创建实体:拉伸工具,为实体拉伸选择一个合适的平面,选择一条边作为草绘参照,进入草绘模式,绘制如下界面草图22:按下鼠标中键退出草绘模式,弹出编辑拉伸对话框,类型设为指定深度,深度设为523:点击确定,完成拉伸24:创建一个基本平面,作为三角形肋板的镜像平面点击创建基准平面工具,点选偏移所参照的平面,通过输入大小来设定偏移,确定偏移方向,输入偏移大小为30,点击确定完成基准平面的创建25:点击创建镜像工具,选择刚刚创建的基准平面为镜像平面,并选择保留原几何,完成镜像至此已经完成部件模型的创建(二)定义材料和界面属性1:模块:属性2:点击创建材料工具,弹出编辑材料对话框,名称改为Steel,通用→密度,输入密度为7.85e-93:力学→弹性→弹性,输入弹性模量为2.1e5,输入泊松比为0.3,点击确定4:点击创建截面工具,弹出创建截面对话框,将名称修改为Basis_Section,类别为实体,类型为均质,点击继续,在弹出的编辑截面对话框中选择确定5:点击指派截面工具,选择整个部件为要指派截面的区域,点击完成,弹出编辑截面指派对话框,点击确定(三)生成装配件1:模块:装配2:点击创建实例工具,弹出创建实例对话框,点击确定(四)定义分析步和指定输出要求1:模块:分析步2:点击创建分析步工具,弹出创建分析步对话框,修改名称为Basic_Load,程序类型选择通用,静力、通用。
一、引言时程分析法是对结构动力方程直接进行逐步积分求解的一种动力分析方法。
时程分析法将地震波按时段进行数值化后,输入结构体系的振动微分方程,采用直接积分法计算出结构在整个强震时域中的振动状态全过程,给出各个时刻各个杆件的内力和变形。
现已成为多数国家抗震设计规范或规程的分析方法之一。
二、有限元软件ABAQUS简介ABAQUS是美国ABAQUS公司(原名HKS公司.即Hibbitt,Karlsson&Sorensen,INC.2005年被法国达索公司收购,2007年公司更名为SIMULIA)。
ABAQUS已成为国际上最先进的大型通用有限元力学分析软件之一。
ABAQUS是一套功能强大的进行工程模拟的有限元软件。
其解决问题的范围从相对简单的线性分析到许多复杂的非线性问题.ABAQUS拥有CAE工业领域最为广泛的材料模型,它可以模拟绝大部分工程材料的线形和非线形行为,可以进行结构的静态和动态分析,如应力、变形、振动、热传导以及对流等.也可以模拟广泛的材料性能,如金属、橡胶、塑料、弹性泡沫等,而且任何一种材料都可以和任何一种单元或复合材料的层一起用于任何合适的分析类型.三、模型建立与求解1、PartCreate Part:Name:Ban,3D,Deformable, Shell ,Planar,输入坐标创建一个18X9m的壳部件,作为混凝土楼板部件;Create Part:Name:Zhu,3D,Deformable, Wire ,Planar,输入坐标创建一个长3m线部件,作为柱部件;Create Part:Name:Liang,3D,Deformable, Wire ,Planar,输入坐标创建一个长6X3m,宽4.5X2m的线网部件,作为梁网部件;2、 SectionCreate Material:Name:steel,General,Density 7800;Mechanical,Elasticity,Young’s Modulus 2.1e11,Poisson’ Ratio 0。
一、引言欧阳家百(2021.03.07)时程分析法是对结构动力方程直接进行逐步积分求解的一种动力分析方法。
时程分析法将地震波按时段进行数值化后,输入结构体系的振动微分方程,采用直接积分法计算出结构在整个强震时域中的振动状态全过程,给出各个时刻各个杆件的内力和变形。
现已成为多数国家抗震设计规范或规程的分析方法之一。
二、有限元软件ABAQUS简介ABAQUS是美国ABAQUS公司(原名HKS公司.即Hibbitt,Karlsson&Sorensen,INC.2005年被法国达索公司收购,2007年公司更名为SIMULIA)。
ABAQUS已成为国际上最先进的大型通用有限元力学分析软件之一。
ABAQUS是一套功能强大的进行工程模拟的有限元软件。
其解决问题的范围从相对简单的线性分析到许多复杂的非线性问题。
ABAQUS拥有CAE工业领域最为广泛的材料模型,它可以模拟绝大部分工程材料的线形和非线形行为,可以进行结构的静态和动态分析,如应力、变形、振动、热传导以及对流等。
也可以模拟广泛的材料性能,如金属、橡胶、塑料、弹性泡沫等,而且任何一种材料都可以和任何一种单元或复合材料的层一起用于任何合适的分析类型。
三、模型建立与求解1、PartCreate Part:Name:Ban,3D,Deformable, Shell ,Planar,输入坐标创建一个18X9m的壳部件,作为混凝土楼板部件;Create Part:Name:Zhu,3D,Deformable,Wire ,Planar,输入坐标创建一个长3m线部件,作为柱部件;Create Part:Name:Liang,3D,Deformable,Wire ,Planar,输入坐标创建一个长6X3m,宽4.5X2m的线网部件,作为梁网部件;2、 SectionCreate Material:Name:steel,General,Density 7800;Mechanical,Elasticity,Young’s Modulus 2.1e11,Poisson’ Ratio 0.3;Create Material:Name:concrete,General,Density 2500;Mechanical,Elasticity,Young’s Modulus 3e10,Poisson’Ratio 0.3。
Create Profile:分别创建梁剖面Liang Profile和柱剖面Zhu Profile为0.3X0.3m,厚0.1m和0.6X0.6m,厚0.1m。
Create Section:板截面:Name:Ban Section,Shell,Homogeneous,Shell Thickness Value:0.1,Material:concrete;梁截面:Name:Liang Section,Profile name:Liang Profile,Material name:steel;柱截面:Name:Zhu Section,Profile name:Zhu Profile,Material name:steel。
Assign Section:分别把创建好的板、梁和柱截面指派到板、梁、柱部件上,梁和柱截面指派时同时指派部件方向。
3、AssemblyCreate Instance:选择板、梁、柱部件,用Dependent方式生成实体;将柱子旋转,使其与板和梁平面垂直;阵列柱子4X3=12根,将12根柱子分别移到梁网交点处,完成一层框架的构造;将构造好的一层框架阵列成3X2=6个框架层(阵列时,注意各个框架层间的距离应该稍大,这样利于平移框架),平移各层框架层,完成6层框架模型的构造。
4、StepCreate Step:Step-1,General,Static General,重力分析,结果作为时程分析数据;Step-2,General,Dynamic Implicit,Time period:10。
5、InteractionCreate constraint:Tie,选择板、梁、柱接触点、面,绑定。
6、LoadCreate Boundary Condition:Name:BC-1,Step:Step-2,Mechanical,Displacement/Rotation,选择柱低端,约束住柱低端除X方向外的5个自由度;Create Boundary Condition:Name:BC-2,Step:Step-2,Mechanical,Acceleration/Angular acceleration,选择柱低端,确定;勾选A1:1(加速度单位m/s²),在Amplitude中创建一个Name:ac,Type:Tabular,Time span:Toal time,复制EICentro数据到Amplitude Data,OK,Amplitude右框下拉选择刚才建好的ac加速度,OK。
7、MeshSeed Part Instance;Assign Mesh Controls:Quard,Free,Advancing front;Assign Element Type:Standard,Linear,Reduced Integration,壳单元S4R;Mesh Part;最后Verify Mesh检查网格是否错误。
8、JobCreate Job;Write Input(输出inp文件),Submit;打开Monitor 观察是否有错误、警报;点击Result进入可视化窗口Visualization。
9、VisualizationResult,Step/Frame分析步/帧,查看各时程变形情况;Tools,XY Data,Create,ODB field output,plot顶层楼板位移时程曲线,各层楼板加速度时程曲线,各层楼板位移曲线。
10、数据后处理在Visualization,XY Data Manager中点击Edit,输出各层楼板位移数据到Excel中,利用Excel:相邻两层的位移作差后得到层间相对位移,画出以Time为横轴、层间相对位移为纵轴画出各层层间位移时程曲线;然后层间相对位移除以层高得到层间位移角,画出以Time为横轴、层间位移角为纵轴画出各层层间位移角时程曲线。
四、结果展示钢框架内力云图1钢框架内力云图2钢框架内力云图3顶层楼板位移时程曲线底端加速度时程曲线一层楼板加速度时程曲线二层楼板加速度时程曲线三层楼板加速度时程曲线四层楼板加速度时程曲线五层楼板加速度时程曲线顶层楼板加速度时程曲线底端与一层楼板层间位移一层与二层楼板层间位移二层与三层楼板层间位移三层与四层楼板层间位移四层与五层楼板层间位移五层与顶层楼板层间位移底端与一层楼板层间位移角一层与二层楼板层间位移角二层与三层楼板层间位移角三层与四层楼板层间位移角四层与五层楼板层间位移角五层与六层楼板层间位移角五、总结与分析通过分析数据:①该钢框架模型顶层楼板最大位移发生在9.904615s时,值为0.521279;根据《》:多层柱顶位移H/500,对于该模型为18m/500=0.036m,顶层楼板最大位移超出规范要求。
②该钢框架模型最大加速度:底端发生在2.13928s,值为3.40922,一层楼板发生在 4.53904s,值为7.02656,二层楼板发生在2.64224s,值为7.58957,三层楼板发生在 3.22835s,值为8.61059,四层楼板发生在4.78149s,值为11.1651,五层楼板发生在 4.78149s,值为12.4091,顶层楼板发生在 4.53904,值为13.9828。
③该钢框架模型层间位移:底端与一层楼板最大值为0.006104,一层与二层楼板最大值为0.025158,二层与三层楼板最大值为0.016052,三层与四层楼板最大值为0.021461,四层与五层楼板最大值为0.012741,五层与六层楼板最大值为0.016035;根据《》:多层层间位移限值h/400,对于该模型为3m/400=0.0.0075m,最大层间位移超出规范要求。
④该钢框架模型层间位移角:底端与一层楼板最大值为0.002035,一层与二层楼板最大值为0.008386,二层与三层楼板最大值为0.005356,三层与四层楼板最大值为0.007154,四层与五层楼板最大值为0.004247,五层与六层楼板最大值为0.005345;根据《》GB50011-2010第5.5.1条规定:“地震作用下,多、弹性限值为1/250”,最大层间位移角超出规范要求。
⑤由以上分析,顶层楼板位移、层间位移、层间位移角均超出规范要求,作者对此结果进行了简单分析,造成此结果可能原因有:没有对时程地震波EICentro波进行系数调整,以满足抗震规范;超过加速度时程的加速度没有按规范将其最大值调整为0.35m/s²;对于输入的地震加速度时程曲线没有满足抗震规范,要满足地震动三要素:频谱特性、有效峰值和持续时间要符合规定的要求。
⑥对于此钢框架结构,其抗侧刚度不足,应通过增加柱截面、梁高和梁宽等措施提高结构抗侧刚度。
⑦以上分析可能不够准确,望读者批评指出。
参考文献[1]王勖成,有限单元法.北京:清华大学出版社,2003[2] 马晓峰,ABAQUS6.11有限元分析从入门到精通. 北京:清华出版社,2013[3] 曹金凤等,ABAQUS有限元分析常见问题与解答. 北京:机械工业出版社,2010[4]徐珂.ABAQUS建筑结构分析应用.北京:中国建筑出版社出版社,2013命令流*Heading** Job name: gkjdzfzfx Model name: Model-1** Generated by: Abaqus/CAE 6.13-1*Preprint, echo=NO, model=NO, history=NO, contact=NO**** PARTS***Part, name=Ban*Node1, 9., 4.5, 0.2, 7.19999981, 4.5, 0.3, 5.4000001, 4.5, 0.4, 3.5999999, 4.5, 0.5, 1.79999995, 4.5, 0.6, 0., 4.5, 0.7, -1.79999995, 4.5, 0.8, -3.5999999, 4.5, 0.9, -5.4000001, 4.5, 0.10, -7.19999981, 4.5, 0.11, -9., 4.5, 0.12, 9., 2.70000005, 0.13, 7.19999981, 2.70000005, 0.14, 5.4000001, 2.70000005, 0.15, 3.5999999, 2.70000005, 0. 16, 1.79999995, 2.70000005, 0. 17, 0., 2.70000005, 0.18, -1.79999995, 2.70000005, 0. 19, -3.5999999, 2.70000005, 0. 20, -5.4000001, 2.70000005, 0. 21, -7.19999981, 2.70000005, 0. 22, -9., 2.70000005, 0.23, 9., 0.899999976, 0.24, 7.19999981, 0.899999976, 0. 25, 5.4000001, 0.899999976, 0. 26, 3.5999999, 0.899999976, 0. 27, 1.79999995, 0.899999976, 0. 28, 0., 0.899999976, 0.29, -1.79999995, 0.899999976, 0. 30, -3.5999999, 0.899999976, 0. 31, -5.4000001, 0.899999976, 0. 32, -7.19999981, 0.899999976, 0. 33, -9., 0.899999976, 0.34, 9., -0.899999976, 0.35, 7.19999981, -0.899999976, 0. 36, 5.4000001, -0.899999976, 0. 37, 3.5999999, -0.899999976, 0. 38, 1.79999995, -0.899999976, 0. 39, 0., -0.899999976, 0.40, -1.79999995, -0.899999976, 0. 41, -3.5999999, -0.899999976, 0. 42, -5.4000001, -0.899999976, 0. 43, -7.19999981, -0.899999976, 0. 44, -9., -0.899999976, 0.45, 9., -2.70000005, 0.46, 7.19999981, -2.70000005, 0. 47, 5.4000001, -2.70000005, 0. 48, 3.5999999, -2.70000005, 0. 49, 1.79999995, -2.70000005, 0. 50, 0., -2.70000005, 0.51, -1.79999995, -2.70000005, 0. 52, -3.5999999, -2.70000005, 0. 53, -5.4000001, -2.70000005, 0. 54, -7.19999981, -2.70000005, 0. 55, -9., -2.70000005, 0.56, 9., -4.5, 0.57, 7.19999981, -4.5, 0.58, 5.4000001, -4.5, 0. 59, 3.5999999, -4.5, 0. 60, 1.79999995, -4.5, 0. 61, 0., -4.5, 0.62, -1.79999995, -4.5, 0. 63, -3.5999999, -4.5, 0. 64, -5.4000001, -4.5, 0. 65, -7.19999981, -4.5, 0. 66, -9., -4.5, 0.*Element, type=S4R1, 1, 2, 13, 122, 2, 3, 14, 133, 3, 4, 15, 144, 4, 5, 16, 155, 5, 6, 17, 166, 6, 7, 18, 177, 7, 8, 19, 188, 8, 9, 20, 199, 9, 10, 21, 2010, 10, 11, 22, 2111, 12, 13, 24, 2312, 13, 14, 25, 2413, 14, 15, 26, 2514, 15, 16, 27, 2615, 16, 17, 28, 2716, 17, 18, 29, 2817, 18, 19, 30, 2918, 19, 20, 31, 3019, 20, 21, 32, 3120, 21, 22, 33, 3221, 23, 24, 35, 3422, 24, 25, 36, 3523, 25, 26, 37, 3624, 26, 27, 38, 3725, 27, 28, 39, 3826, 28, 29, 40, 3927, 29, 30, 41, 4028, 30, 31, 42, 4129, 31, 32, 43, 4230, 32, 33, 44, 4331, 34, 35, 46, 4532, 35, 36, 47, 4633, 36, 37, 48, 4734, 37, 38, 49, 4835, 38, 39, 50, 4936, 39, 40, 51, 5037, 40, 41, 52, 5138, 41, 42, 53, 5239, 42, 43, 54, 5340, 43, 44, 55, 5441, 45, 46, 57, 5642, 46, 47, 58, 5743, 47, 48, 59, 5844, 48, 49, 60, 5945, 49, 50, 61, 6046, 50, 51, 62, 6147, 51, 52, 63, 6248, 52, 53, 64, 6349, 53, 54, 65, 6450, 54, 55, 66, 65*Nset, nset=Set-1, generate1, 66, 1*Elset, elset=Set-1, generate1, 50, 1*Nset, nset=Set-4, generate1, 66, 1*Elset, elset=Set-4, generate1, 50, 1*Nset, nset=Set-5, generate1, 66, 1*Elset, elset=Set-5, generate1, 50, 1** Section: Ban Section*Shell Section, elset=Set-5, material=concrete 0.1, 5*End Part***Part, name=Liang*Node1, -3., 4.5, 0.2, 3., 4.5, 0.3, 9., 4.5, 0.4, 9., 0., 0.5, 9., -4.5, 0.6, 3., -4.5, 0.7, -3., -4.5, 0.8, -9., -4.5, 0.9, -9., 0., 0.10, -3., 0., 0.11, 3., 0., 0.12, -9., 4.5, 0.13, -2.4000001, 4.5, 0. 14, -1.79999995, 4.5, 0. 15, -1.20000005, 4.5, 0. 16, -0.600000024, 4.5, 0. 17, 0., 4.5, 0.18, 0.600000024, 4.5, 0. 19, 1.20000005, 4.5, 0. 20, 1.79999995, 4.5, 0. 21, 2.4000001, 4.5, 0. 22, 3.5999999, 4.5, 0. 23, 4.19999981, 4.5, 0. 24, 4.80000019, 4.5, 0. 25, 5.4000001, 4.5, 0. 26, 6., 4.5, 0.27, 6.5999999, 4.5, 0. 28, 7.19999981, 4.5, 0. 29, 7.80000019, 4.5, 0. 30, 8.39999962, 4.5, 0. 31, 9., 3.9375, 0.32, 9., 3.375, 0.33, 9., 2.8125, 0.34, 9., 2.25, 0.35, 9., 1.6875, 0.36, 9., 1.125, 0.37, 9., 0.5625, 0.38, 9., -0.5625, 0.39, 9., -1.125, 0.40, 9., -1.6875, 0.41, 9., -2.25, 0.42, 9., -2.8125, 0.43, 9., -3.375, 0.44, 9., -3.9375, 0.45, 8.39999962, -4.5, 0. 46, 7.80000019, -4.5, 0. 47, 7.19999981, -4.5, 0. 48, 6.5999999, -4.5, 0. 49, 6., -4.5, 0.50, 5.4000001, -4.5, 0.51, 4.80000019, -4.5, 0. 52, 4.19999981, -4.5, 0. 53, 3.5999999, -4.5, 0. 54, 2.4000001, -4.5, 0. 55, 1.79999995, -4.5, 0. 56, 1.20000005, -4.5, 0. 57, 0.600000024, -4.5, 0. 58, 0., -4.5, 0.59, -0.600000024, -4.5, 0. 60, -1.20000005, -4.5, 0. 61, -1.79999995, -4.5, 0. 62, -2.4000001, -4.5, 0. 63, -3.5999999, -4.5, 0. 64, -4.19999981, -4.5, 0. 65, -4.80000019, -4.5, 0. 66, -5.4000001, -4.5, 0. 67, -6., -4.5, 0.68, -6.5999999, -4.5, 0. 69, -7.19999981, -4.5, 0. 70, -7.80000019, -4.5, 0. 71, -8.39999962, -4.5, 0. 72, -9., -0.5625, 0.73, -9., -1.125, 0.74, -9., -1.6875, 0.75, -9., -2.25, 0.76, -9., -2.8125, 0.77, -9., -3.375, 0.78, -9., -3.9375, 0.79, -3.5999999, 0., 0. 80, -4.19999981, 0., 0. 81, -4.80000019, 0., 0. 82, -5.4000001, 0., 0. 83, -6., 0., 0.84, -6.5999999, 0., 0. 85, -7.19999981, 0., 0. 86, -7.80000019, 0., 0. 87, -8.39999962, 0., 0. 88, -3., -3.9375, 0.89, -3., -3.375, 0.90, -3., -2.8125, 0.91, -3., -2.25, 0.92, -3., -1.6875, 0.93, -3., -1.125, 0.96, -3., 1.125, 0.97, -3., 1.6875, 0.98, -3., 2.25, 0.99, -3., 2.8125, 0.100, -3., 3.375, 0.101, -3., 3.9375, 0.102, 2.4000001, 0., 0. 103, 1.79999995, 0., 0. 104, 1.20000005, 0., 0. 105, 0.600000024, 0., 0. 106, 0., 0., 0.107, -0.600000024, 0., 0. 108, -1.20000005, 0., 0. 109, -1.79999995, 0., 0. 110, -2.4000001, 0., 0. 111, 3., -3.9375, 0.112, 3., -3.375, 0.113, 3., -2.8125, 0.114, 3., -2.25, 0.115, 3., -1.6875, 0.116, 3., -1.125, 0.117, 3., -0.5625, 0.118, 3., 0.5625, 0.119, 3., 1.125, 0.120, 3., 1.6875, 0.121, 3., 2.25, 0.122, 3., 2.8125, 0.123, 3., 3.375, 0.124, 3., 3.9375, 0.125, 8.39999962, 0., 0. 126, 7.80000019, 0., 0. 127, 7.19999981, 0., 0. 128, 6.5999999, 0., 0. 129, 6., 0., 0.130, 5.4000001, 0., 0. 131, 4.80000019, 0., 0. 132, 4.19999981, 0., 0. 133, 3.5999999, 0., 0. 134, -9., 3.9375, 0.135, -9., 3.375, 0.139, -9., 1.125, 0.140, -9., 0.5625, 0.141, -8.39999962, 4.5, 0. 142, -7.80000019, 4.5, 0. 143, -7.19999981, 4.5, 0. 144, -6.5999999, 4.5, 0. 145, -6., 4.5, 0.146, -5.4000001, 4.5, 0. 147, -4.80000019, 4.5, 0. 148, -4.19999981, 4.5, 0. 149, -3.5999999, 4.5, 0. *Element, type=B311, 1, 132, 13, 143, 14, 154, 15, 165, 16, 176, 17, 187, 18, 198, 19, 209, 20, 2110, 21, 211, 2, 2212, 22, 2313, 23, 2414, 24, 2515, 25, 2616, 26, 2717, 27, 2818, 28, 2919, 29, 3020, 30, 321, 3, 3122, 31, 3223, 32, 3324, 33, 3425, 34, 3526, 35, 3627, 36, 3728, 37, 431, 39, 40 32, 40, 41 33, 41, 42 34, 42, 43 35, 43, 44 36, 44, 5 37, 5, 45 38, 45, 46 39, 46, 47 40, 47, 48 41, 48, 49 42, 49, 50 43, 50, 51 44, 51, 52 45, 52, 53 46, 53, 6 47, 6, 54 48, 54, 55 49, 55, 56 50, 56, 57 51, 57, 58 52, 58, 59 53, 59, 60 54, 60, 61 55, 61, 62 56, 62, 7 57, 7, 63 58, 63, 64 59, 64, 65 60, 65, 66 61, 66, 67 62, 67, 68 63, 68, 69 64, 69, 70 65, 70, 71 66, 71, 8 67, 9, 72 68, 72, 73 69, 73, 74 70, 74, 75 71, 75, 7674, 78, 8 75, 10, 79 76, 79, 80 77, 80, 81 78, 81, 82 79, 82, 83 80, 83, 84 81, 84, 85 82, 85, 86 83, 86, 87 84, 87, 9 85, 7, 88 86, 88, 89 87, 89, 90 88, 90, 91 89, 91, 92 90, 92, 93 91, 93, 94 92, 94, 1093, 10, 95 94, 95, 96 95, 96, 97 96, 97, 98 97, 98, 99 98, 99, 100 99, 100, 101 100, 101, 1 101, 11, 102 102, 102, 103 103, 103, 104 104, 104, 105 105, 105, 106 106, 106, 107 107, 107, 108 108, 108, 109。