OpenSees安装及使用.ppt

OpenSEES（Open System for Earthquake Engineering Simulation）由美国国家自然科学基金（NSF）资助、西部大学联盟“太平洋地震工程研究中心”（Pacific Earthquake Engineering Research Center, PEER）主导、以加州大学伯克利分校为主研发而成的、用于结构和岩土方面地震反应模拟的一个较为全面且不断发展的开放式模拟平台。,OpenSees程序自1999年正式推出以来，已广泛应用于太平洋地震工程研究中心和美国其它一些大学及科研机构的科研项目中，较好地模拟了包括钢筋混凝土结构、桥梁、岩土工程在内的实际工程和振动台试验项目。 2004年，OpenSees被美国“地震工程模拟网络”（Network for Earthquake Engineering Simulation，NEES）所采用，作为NEES计划的虚拟仿真平台，证明其具有较好的非线性数值模拟精度。,OpenSees的安装,/index.php,OpenSees的安装,Edit TCL language,Toolbar,Run,TCL Editor Software,外部支持软件,Installation guide for OpenSees Navigator Unzip and saved the file to any directory. Set MATLAB path to the OpenSeesNavigator folder (add with Subfolders). Type OpenSeesNavigator in the MATLAB command window. Now you are officially an OpenSees Navigator user! Remember to register to the user list to receive future updates.,外部支持软件,OpenSees的帮助文件,OpenSees的帮助文件 例题,解决研究中的实际问题 钢筋混凝土柱抗震性能的模拟,为什么要研究柱的抗震性能？,什么是伪静力试验？拟动力试验？振动台试验？,构件的伪静力试验要得到什么试验结果？,构件伪静力试验的主要影响因素？,轴压比、剪跨比、纵筋配筋率、体积配箍率等,等效刚度（effective stiffness）的概念,等效刚度（effective stiffness）的概念,What is drift angle ?,卸载刚度（Unloading stiffness）的概念,可以考察试件的刚度退化(Stiffness degradation),累积能量（Cumulative energy）的概念,可以考察试件的耗散能量，进而判断构件的抗震性能。,Test of Reinforced Concrete (RC ) columns 试验的关键问题：轴力的随动,我院四连杆机构压弯构件试验装置,我院四连杆机构压弯构件试验装置,OPENSEES中纤维模型（Fiber Model）的应用及原理,单元划分较细，结构总刚矩阵维数偏大，对硬件的要求较高； 对于分离式模型，由于引入的单元种类较多，分析的收敛性更加 不好控制，只适于局部应力状态分析； 较细的单元网格有时仍不一定能较好地解释钢筋混凝土结构构件 受力时的某些力学特性，如塑性铰的出现和扩展，梁柱截面及节 点区域处的不连续转角等。,采用常规有限元（实体有限元）的局限性,把构件纵向分为多个微段，把构件截面划分为平面网格，每一网格的中心为数值积分点，网格的纵向微段即定义为 Fiber（纤维）。,纤维模型的基本原理,首先通过“有限条带法（分层法）”来介绍纤维模型的基本原理,利用分层法进行梁M-单调曲线的分析,分层法（纤维模型的特例）,梁M-曲线的特点,曲率延性系数,一、基本假定 (Basic Assumptions),平截面假定（ The section remaining plane section ）； 钢筋与混凝土充分粘结，无相对滑移，变形协调； 构件的剪跨比远大于3，可以忽略剪切变形； 忽略构件翘曲和扭转的影响。,平衡条件：,几何关系：,物理关系：,曲率,应变,应力,弯矩,二、分析思路,三、截面的条带划分及应变分布,四、材料的本构关系,在分层模型中，截面被划分成若干有限条带，每个条带内只考虑一维应力状态。,钢筋应力-应变关系表达式,我国混凝土结构设计规范（GBJ50010-2002）的表达式,当n 大于2.0 时 取2.0,式中:,曲线考虑了混凝土在受压时峰值应变与极限压应变的变化情况，它由抛物线的上升段和水平段两部分组成，应用简便，但是没有考虑混凝土下降段。,Hongnestad 建议在理论分析时取ecu =0.0038，在设计中 ecu =0.003，并建议e0=2s0/E0，E0 为混凝土的受压初始 弹性模量，s0建议取为0.85 fc (混凝土圆柱体抗压强度)。,美国学者 Hongnestad 提出的表达式,是目前世界上应用最广泛的普通混凝土应力应变曲线之一。但峰值点处存在尖点，其导数不连续，给数值求解带来一定的困难；,日本学者 Desayi 和 Krishnan 提出的表达式,该公式开创了用统一式子表达上升段和下降段，很多学者在此基础上进行修正和完善。,1964年Saenz 等人的表达式,Elwi 和Murray 对Saenz表达式进行修正,根据下列5个控制条件确定A、B、C、D四个参数,自然满足,求解A、B、C、D四个常数,模型参数确定,这一公式在有限元分析中应用很广，采用统一式子表达上升段和下降段，在ADINA中基本采用此应力应变本构模型。,1971年Sargin 对Saenz表达式进行修正,适合用来描述各种标号混凝土及箍筋约束混凝土的本构关系，被CEB-FIP MC90欧洲规范采用。,若采用无量纲坐标x=e/e0，y=s/fc，则混凝土应力-应变全曲线的几何特征必须满足:,清华大学过镇海提出的应力-应变全曲线表达式,a为下降段参数，一般应有1.5a3；,过镇海等人建议的表达式考虑的影响因素更加全面，体现了应力应变曲线的几何特点，峰值点处导数连续，是一种精细型的本构关系，适用于普通混凝土和高强混凝土，在我国得到广泛认可。,双参数模型,过镇海箍筋约束高强混凝土双参数本构模型,S 为箍筋间距；b 为矩形截面的短边长度,从加载直到最终破坏，分析截面应力分布、弯矩与曲率的关系 具体分析步骤如下： 给定一截面曲率f （由小到大）； 假定受压边缘混凝土应变值ec； 由平截面假定，确定截面应变分布和钢筋应变es； 利用物理关系，确定C和yc、Tc和yt、Ts； 验算是否满足轴力平衡条件，如满足，进行； 如不满足，修正ec后，重新分析；