钢筋混凝土梁ansys分析附命令流
基于ANSYS的钢筋混凝土结构试验有限元分析共3篇
基于ANSYS的钢筋混凝土结构试验有限元分析共3篇基于ANSYS的钢筋混凝土结构试验有限元分析1混凝土结构是我们生活和工作环境中不可或缺的部分。
为了保证结构的安全性和耐久性，需要进行大量的试验和分析。
钢筋混凝土结构试验有限元分析是其中一种方法，本文将介绍如何基于ANSYS进行试验有限元分析。
1、前期准备工作进行钢筋混凝土结构试验有限元分析前，需要进行一些前期准备工作。
首先要确定模型的尺寸和几何形状，包括梁的长度、宽度和高度，钢筋的数量和材料等信息。
其次是建立材料模型。
钢筋和混凝土的本构关系可以参考各种规范和文献，例如ACI318和EHE等。
最后是进行荷载和边界条件的设置。
这些参数可以根据试验的要求进行设定。
2、建立有限元模型通过ANSYS软件建立钢筋混凝土结构的有限元模型。
其中，混凝土部分采用可压缩性线性弹性模型；钢筋采用弹塑性模型，可以考虑材料的塑性性质。
首先，选择适当的元素类型，包括梁单元和实体单元。
对于梁单元，要选择适当的截面类型和断面参数。
对于实体单元，要确定网格的大小和形状。
然后，按照模型的几何形状和材料参数设置单元类型和属性。
最后，进行单元的划分和网格生成，调整边界条件，使其与试验条件保持一致。
3、分析和结果在模型准备就绪之后，进行分析和结果的处理。
首先，定义荷载和边界条件，可以模拟多种加载模式，例如单点荷载、均布荷载、自重等。
然后，进行静态分析或动态分析。
静态分析可以计算结构的变形、应力和应变等参数；动态分析可以模拟结构在地震、风等自然灾害下的响应。
最后，进行结果的处理和分析。
包括可视化、动画演示、应力云图、位移云图等，能够对计算结果进行全方位的检查和分析。
综上所述，基于ANSYS的钢筋混凝土结构试验有限元分析是一种非常有用的手段，可以帮助工程师更准确地评估结构的安全性和耐久性。
它具有良好的可靠性和可操作性，可在较短的时间内快速建立模型和分析结果。
基于ANSYS的钢筋混凝土结构试验有限元分析2钢筋混凝土结构是目前建筑工程最常用的一种结构形式，其优点在于承载能力强、耐久性好、施工方便等。
钢筋混凝土梁的弹塑性分析ansys命令流
!(1)工作环境设置/FILENAME,RC-BEAM !指定工作文件名/TITLE,ALAL YSIS OF A RC-BEAM !指定图形标题!(2)进行前处理器/prep7!(3)定义单元类型ET,1,LINK8 !定义钢筋单元ET,2,SOLID65 !定义混凝土单元ET,3,MESH200 !用于拉伸成体单元KEYOPT,3,1,6!(4)定义钢筋截面积r,1,28.3r,2,50.3r,3,314.1!(5)为solid65单元定义一个实参数组r,4,!(6)定义混凝土材料MP,EX,2,2.55E10MP,PRXY,2,0.3TB,CONC,2,1,9, !定义混凝土的破坏参数TBDA TA,,0.3,0.55,1.55E6,-1,,TBDA TA,,,,0.6!(7)定义钢筋材料模型及参数mp,ex,2,2e5 !纵向受拉钢筋材料mp,prxy,2,0.3tb,bkin,2,1,2,1tbdata,,350mp,ex,3,2e5 !横向箍筋，架立钢筋材料mp,prxy,3,0.25tb,bkin,3,1,2,1tbdata,,200!(8) 创建以及复制节点/pnum,node,1/pnum,elem,1n,1n,9,200fill,1,9ngen,11,9,1,9,1,,30ngen,11,99,1,99,1,,,-150/view,1,1,1,1!(9)建立箍筋单元type,1real,1mat,3!水平箍筋*do,i,11,16,1e,i,i+1e,i+(83-11),i+(83-11)+1*enddo!坚直箍筋*do,i,11,74,9e,i,i+9e,i+6,i+6+9*enddo!产生整个模型的箍筋egen,11,99,all!(10)建立架立筋以及纵筋单元!创建上部的架立钢筋单元*do,i,83,node(25,270,-1500+150),99 e,i,i+99e,i+6,i+6+99*enddo!纵向受拉钢筋单元的属性type,1real,3mat,2!创建纵筋单元*do,i,11,node(25,30,-1500+150),99 e,i,i+99e,i+3,i+3+99e,i+6,i+6+99*enddo/view,1,1,1,1/pnum,elem,0/pnum,node,0/eshape,1eplot!(11)建立混凝土剖面并划分网格k,1K,2,200,k,3,200,300k,4,,300a,1,2,3,4lsel,s,loc,y,0lsel,a,loc,y,300lesize,all,,,8lsel,alllsel,s,loc,x,0lsel,a,loc,x,200lesize,all,,,10type,3amesh,all!(12)拉伸形成混凝土单元type,2real,3mat,1extopt,esize,20extopt,aclear,1vext,all,,,,,-1500/pnum,mat,1/pnum,node,0/pnum,elem,0eplotallsel!(13)合并压缩节点编号nummrg,allnumcmp,alleplot!(14)施加支座约束nsel,s,loc,y,0nsel,r,loc,z,-1500+75d,all,uyd,all,ux!(15)施加对称面约束asel,s,loc,z,0da,all,symm!(16)施加载荷nsel,allnsel,s,loc,y,300nsel,r,loc,z,-450d,all,uy,-30allsel!(17)退出前处理器fini!(18) 进入求解器/solu!(19)求解器选项设置nlgeom,onnsubst,200outres,all,allneqit,50pred,oncnvtol,f,,0.05,2,0.5allsel!(20)求解并退出求解器solvefini!(21)进入通用处理器并读入最后一个子步的结果/post1set,last!(22)后处理操作plnsol,u,y,0,1/device,vector,1plcrack,0,0!(23)退出通用后处理器finish。
钢管混凝土ansys命令流
碳纤维加固桥梁结构的健康监测摘要：提出了一种针对老化桥梁的健康监测的可靠性方法。
首先, 根据桥梁的类型提出了临界失效标准。
然后基于这些标准定义安全界限。
假定大量的失效标准是正态分布的随机变量。
基于这样的假设,基本可靠性指标和基本失效概率为每个关键失效模式估计。
已经发现基本失效概率值,桥梁系统失效概率值按这段时间考虑计算。
最后,用系统失效概率来获得桥梁系统的可靠性指标。
系统的可靠性指数用来表达的桥梁在使用期间的状况较为片面。
选择斯里兰卡国家铁路桥网络的一座铁路桥梁作为一个案例研究来说明该可靠度程序。
对于这座桥,失效的最重要的形式是疲劳和腐蚀。
应用程序的建议方法表明,桥梁在现有负载的情况下的现状是令人满意的。
对可靠性的预测,不是为了抵御破坏,而是一种对基础设施管理进行有效的积极的健康监测方式。
关键字：碳纤维结构健康监测状态评估维护可靠性指标失效概率1、介绍在当今迅速变化的世界中,人类的生活质量和国家经济发展依赖于基础设施的数量、质量、和有效性]1[。
在世界各国，民用基础设施通常是需要维护的,康复,或更换的]2[。
特别是对基础设施构成一个相当大的投资的高速公路和相关的地面运输系统，直接影响到一个国家的工业产值。
目前大多数桥梁养护策略基于由桥梁检查员在不同时间间隔进行的视觉检查]3[。
由于人类检测取决于个人检查员,其结果存在一定程度的不确定性。
因此,绝对的依赖检验报告是不可靠的，而且可能导致错误的决策和桥梁拥有者的高维护成本。
基于桥梁的剩余使用寿命维护成本的优化。
另一方面,桥梁的检查员发现结构和非结构化缺陷,破坏和恶化过程的能力不应该被忽视。
事实上,人类经验应该被纳入维护决策,但要有足够的主观的知识。
只在那时,任何可接受的战略才能成为可能。
当前桥梁结构健康监测方法存在许多缺点和普遍的不确定性。
对结构健康监测和桥梁维护的知识不足造成由于频繁的维护活动,不可预见的桥梁损害,和偶然的桥梁失效引起的重大的经济损失。
Ansys经典命令流钢管混凝土拱桥算例
Ansys经典命令流钢管混凝土拱桥算例/prep7/title,the analyse of steel arch bridge!************************************************************** *****************et,1,82 !建立辅助单元作出截面特性CYL4,0.6,0.6,0.6, ,0.586 !以下命令是用于建立一个空心钢管截面CYL4,0.6,2.4,0.6, ,0.586RECTNG,0.35,0.364,0.8,2.2,RECTNG,0.85,0.836,0.8,2.2,asel,allaptn,allADELE,5,8,1 !删除多余的面积allselaadd,allsmrtsize,5amesh,allsecwrite,gg,sect,,1SECTYPE,1,BEAM,MESH, !截面代号1SECOFFSET,CENT,,,SECREAD,'gg','SECT',,MESH !读入空心钢管截面asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replotCYL4,0.6,0.6,0.586 !以下命令是用于建立一个空心钢管内的混凝土截面CYL4,0.6,2.4,0.586RECTNG,0.364,0.836,0.8,2.2,allselaadd,allsmrtsize,5amesh,allsecwrite,hnt,sect,,1SECTYPE,2,BEAM,MESH, !截面代号2SECOFFSET,CENT,,,SECREAD,'hnt','SECT',,MESH !读入混凝土（空心钢管内的）截面asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replot !重新生成图形k,1,-0.5,0 !以下命令是用于建立一个横梁标准截面k,2,-0.75,0.25k,3,-0.75,0.5k,4,-0.55,0.5k,5,-0.55,0.85k,6,0.55,0.85k,7,0.55,0.5k,8,0.75,0.5k,9,0.75,0.25k,10,0.5,0a,1,2,3,4,5,6,7,8,9,10 !生成横梁标准截面1 A-Asmrtsize,5amesh,allsecwrite,hl1,sect,,1SECTYPE,4,BEAM,MESH, !截面代号4SECOFFSET,CENT,,,SECREAD,'hl1','SECT',,MESH !读入横梁标准截面1 A-A asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replot !重新生成图形k,1,-0.5,0 !以下命令是用于建立一个横梁标准截面b-bk,2,-0.5,1k,3,-0.75,1.25k,4,-0.75,1.5k,5,0.75,1.5k,6,0.75,1.25k,7,0.5,1k,8,0.5a,1,2,3,4,5,6,7,8 !生成横梁标准截面2 b-bsmrtsize,5amesh,allsecwrite,hl2,sect,,1SECTYPE,5,BEAM,MESH, !截面代号5SECOFFSET,CENT,,,SECREAD,'hl2','SECT',,MESH !读入横梁标准截面2asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replot !重新生成图形k,1,-0.5,0 !以下命令是用于建立一个横梁标准截面c-ck,2,-0.5,1.04k,3,-0.75,1.04+0.25k,4,-0.75,1.04+0.5k,5,-0.55,1.04+0.5k,6,-0.55,1.04+0.5+0.35k,7,0.55,1.04+0.5+0.35k,8,0.55,1.04+0.5k,9,0.75,1.04+0.5k,10,0.75,1.04+0.25k,11,0.5,1.04k,12,0.5,0k,13,-0.2,0.58k,14,-0.2,1.04+0.5-0.3k,15,0.2,1.04+0.5-0.3k,16,0.2,0.58a,1,2,3,4,5,6,7,8,9,10,11,12a,13,14,15,16asba,1,2 !横梁标准截面3 C-Csmrtsize,5amesh,allsecwrite,hl3,sect,,1SECTYPE,6,BEAM,MESH, ! 截面代号6 SECOFFSET,CENT,,,SECREAD,'hl3','SECT',,MESH !横梁标准截面3asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replot !重新生成图形ldele,allkdele,allk,1,-0.5,0 !以下命令是用于建立一个横梁标准截面d-d k,2,-0.5,1.23k,3,-0.75,1.23+0.25k,4,-0.75,1.23+0.5k,5,-0.55,1.23+0.5k,6,-0.55,1.23+0.5+0.35k,7,0.55,1.23+0.5+0.35k,8,0.55,1.23+0.5k,9,0.75,1.23+0.5k,10,0.75,1.23+0.25k,11,0.5,1.23k,12,0.5,0k,13,-0.2,0.38k,14,-0.3,0.38+0.2k,15,-0.3,1.23+0.5-0.1-0.2k,16,-0.2,1.23+0.5-0.1k,17,0.2,1.23+0.5-0.1k,18,0.3,1.23+0.5-0.1-0.2k,19,0.3,0.38+0.2k,20,0.2,0.38a,1,2,3,4,5,6,7,8,9,10,11,12a,13,14,15,16,17,18,19,20asba,1,2 !横梁标准截面4 D-Dsmrtsize,5amesh,allsecwrite,hl4,sect,,1SECTYPE,7,BEAM,MESH, !截面代号7 SECOFFSET,CENT,,,SECREAD,'hl4','SECT',,MESH !读入横梁标准截面4 asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replot ! 重新生成图形k,1,-0.5,0 !以下命令是用于建立一个横梁标准截面e-ek,2,-0.5,1.23k,3,-0.75,1.23+0.25k,4,-0.75,1.23+0.5k,5,-0.55,1.23+0.5k,6,-0.55,1.23+0.5+0.35k,7,0.55,1.23+0.5+0.35k,8,0.55,1.23+0.5k,9,0.75,1.23+0.5k,10,0.75,1.23+0.25k,11,0.5,1.23k,12,0.5,0a,1,2,3,4,5,6,7,8,9,10,11,12 !横梁标准截面5 E-E smrtsize,5amesh,allsecwrite,hl5,sect,,1SECTYPE,8,BEAM,MESH, ! 截面代号8SECOFFSET,CENT,,,SECREAD,'hl5','SECT',,MESH !横梁标准截面5asel,allaclear,alladele,all,,,1/replotetdele,1 !将辅助单元类型删除/replot!*********************************************************************!定义所有材料特性et,1,beam44 !钢管特性mp,ex,1,2.1e11mp,dens,1,7800mp,prxy,1,0.3n,90000,0,0,30 !参考点et,2,beam44 !拱肋内混凝土特性 50#混凝土mp,ex,2,3.5e10mp,dens,2,2600mp,prxy,2,0.1667et,3,beam44 ! 纵梁特性 30#混凝土mp,ex,3,3.0e10mp,dens,3,2600mp,prxy,3,0.1667et,4,beam44 !横梁特性 30#混凝土mp,ex,4,3.0e10mp,dens,4,2600mp,prxy,4,0.1667et,5,beam44 !风撑特性mp,ex,5,2.1e11mp,dens,5,7800mp,prxy,5,0.3et,6,link10 !吊杆特性（钢绞线）mp,dens,6,7800mp,prxy,6,0.3keyopt,6,3,0 !只拉吊杆et,7,beam44 !盖梁特性 30# 混凝土mp,ex,7,3.0e10mp,dens,7,2600mp,prxy,7,0.1667et,8,beam44 !墩柱特性 40#混凝土mp,ex,8,3.3e10mp,dens,8,2600mp,prxy,8,0.1667et,9,beam44 !横系梁特性 30#混凝土mp,ex,9,3.0e10mp,dens,9,2600mp,prxy,9,0.1667et,10,link10 !系杆特性（钢绞线）mp,ex,10,2.0e11mp,dens,10,7800mp,prxy,10,0.3keyopt,10,3,0 !只拉系杆et,11,beam44 !承台桩基础特性 25#混凝土mp,ex,11,2.85e10mp,dens,11,2500n,90000,0,0,30!************************************************************** ********!开始建立模型!0# 横梁(端横梁)n,1,0,2.025n,2,6.25,2.025n,3,6.25+5.1,2.025n,4,6.25+5.1+0.2+1.9/2,2.025 !12.5米位置n,5,6.25+5.1+2.3,2.025n,6,6.25+5.1+2.3+2.23,2.025nsym,x,10,2,6,1type,4 !横梁r,1(定义一组单元实常数)real,1 (划分单元时指定的实常书)mat,4secnum,7! 单元1到10为0# 横梁e,1,2 !截面7 D-Dsecnum,6e,2,3 !截面6 D-Dsecnum,5e,3,4 !截面5 D-De,4,5 !截面5 D-Dsecnum,4e,5,6 !截面4 a-asecnum,7e,1,12 !截面7 D-Dsecnum,6e,12,13 !截面6 D-Dsecnum,5e,13,14 !截面5 D-De,14,15 !截面5 D-Dsecnum,4e,15,16 !截面4 a-a!1# 端横梁egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2!其余横梁2#~17#egen,17,100,11,20,1,,,,,,0,5.3+1.5! 横梁18#(端横梁)egen,2,100,171,180,1,,,,,,0,4.7+1.5/2+1.25/2cm,hl,elemesel,none!*********************************************************** ! 纵梁单元type,3r,2real,2mat,3SECTYPE,10,beam,rect,zl !纵梁截面代号10secdata,0.4,1.2secnum,10!0# 纵梁(端纵梁）e,1,101e,2,102e,3,103e,5,105e,12,112e,13,113e,15,115!中间的纵梁1#~16#*do,i,0,15,1e,101+i*100,201+i*100e,102+i*100,202+i*100e,104+i*100,204+i*100e,112+i*100,212+i*100e,114+i*100,214+i*100*enddo!17#纵梁(端纵梁）e,1701,1801e,1702,1802e,1703,1803e,1705,1805e,1712,1812e,1713,1813e,1715,1815cm,zl,elemesel,none!************************************************************** *********!定义拱轴线 !请参照《桥梁工程》p299页相关公式*dim,y,array,16 !将纵桥向距离定义为数组*dim,cs,array,16 !**值*dim,z,array,16 !将拱轴高度定义为数组!吊杆上点!首先定义吊杆的y向（顺桥向）坐标y(1)=8.1y(2)=14.9y(3)=21.7y(4)=28.5y(5)=35.3y(6)=42.1y(7)=48.9y(8)=55.7y(9)=62.5!风撑上点!定义风撑的y向（顺桥向）坐标y(11)=21.7y(12)=28.5y(13)=42.1y(14)=48.9y(15)=59.1!0#横梁的y向（顺桥向）坐标y(16)=2.025f1=25 !矢高m1=1.1 !悬链线拱轴系数k=log(m1+sqrt(m1*m1-1))*do,i,1,16,1cs(i)=1-y(i)*2/125z(i)=25-f1/(m1-1)*(cosh(k*cs(i))-1) !拱轴高度*enddo*dim,ygd,array,23 !定义预拱度（是在后续工况中反复调算的结果）初始运行时可将预拱度全部设置为0ygd(1) = 0.40004E-02+1.29E-02ygd(2) = 0.10685E-01+2.22E-02ygd(3) = 0.19483E-01+3.03E-02ygd(4) = 0.29219E-01+3.69E-02ygd(5) = 0.39022E-01+4.23E-02ygd(6) = 0.48243E-01+4.65E-02ygd(7) = 0.56390E-01+4.95E-02ygd(8) = 0.63070E-01+5.15E-02ygd(9) = 0.68014E-01+5.21E-02ygd(10)= 0.71003E-01+5.15E-02ygd(11)= 0.71739E-01+4.95E-02ygd(12)= 0.69962E-01+4.65E-02ygd(13)= 0.65514E-01+4.23E-02ygd(14)= 0.58289E-01+3.69E-02ygd(15)= 0.48248E-01+3.04E-02ygd(16)= 0.35507E-01+2.24E-02ygd(17)= 0.20415E-01+1.23E-02ygd(18)=0.65770E-01 +0.40234E-01ygd(19)=0.69769E-01 +0.41215E-01!生成各个吊杆的Z向节点（竖直方向）ngen,2,10000,104,114,10,0,0,z(1) +ygd(1) !1#吊杆节点ngen,2,10000,204,214,10,0,0,z(2) +ygd(2) !2#吊杆节点ngen,2,10000,304,314,10,0,0,z(3) +ygd(3) !3#吊杆节点ngen,2,10000,404,414,10,0,0,z(4) +ygd(4) !4#吊杆节点ngen,2,10000,504,514,10,0,0,z(5) +ygd(5) !5#吊杆节点ngen,2,10000,604,614,10,0,0,z(6) +ygd(6) !6#吊杆节点ngen,2,10000,704,714,10,0,0,z(7) +ygd(7) !7#吊杆节点ngen,2,10000,804,814,10,0,0,z(8) +ygd(8) !8#吊杆节点ngen,2,10000,904,914,10,0,0,z(9) +ygd(9) !9#吊杆节点ngen,2,10000,1004,1014,10,0,0,z(8) +ygd(10) !10#吊杆节点ngen,2,10000,1104,1114,10,0,0,z(7) +ygd(11) !11#吊杆节点ngen,2,10000,1204,1214,10,0,0,z(6) +ygd(12) !12#吊杆节点ngen,2,10000,1304,1314,10,0,0,z(5) +ygd(13) !13#吊杆节点ngen,2,10000,1404,1414,10,0,0,z(4) +ygd(14) !14#吊杆节点ngen,2,10000,1504,1514,10,0,0,z(3) +ygd(15) !15#吊杆节点ngen,2,10000,1604,1614,10,0,0,z(2) +ygd(16) !16#吊杆节点ngen,2,10000,1704,1714,10,0,0,z(1) +ygd(17) !17#吊杆节点!说明：之所以重复定义吊杆节点是因为以后吊杆转动要进行放松ngen,2,3000,104,114,10,0,0,0 !1#吊杆下节点ngen,2,3000,204,214,10,0,0,0 !2#吊杆下节点ngen,2,3000,304,314,10,0,0,0 !3#吊杆下节点ngen,2,3000,404,414,10,0,0,0 !4#吊杆下节点ngen,2,3000,504,514,10,0,0,0 !5#吊杆下节点ngen,2,3000,604,614,10,0,0,0 !6#吊杆下节点ngen,2,3000,704,714,10,0,0,0 !7#吊杆下节点ngen,2,3000,804,814,10,0,0,0 !8#吊杆下节点ngen,2,3000,904,914,10,0,0,0 !9#吊杆下节点ngen,2,3000,1004,1014,10,0,0,0 !10#吊杆下节点ngen,2,3000,1104,1114,10,0,0,0 !11#吊杆下节点ngen,2,3000,1204,1214,10,0,0,0 !12#吊杆下节点ngen,2,3000,1304,1314,10,0,0,0 !13#吊杆下节点ngen,2,3000,1404,1414,10,0,0,0 !14#吊杆下节点ngen,2,3000,1504,1514,10,0,0,0 !15#吊杆下节点ngen,2,3000,1604,1614,10,0,0,0 !16#吊杆下节点ngen,2,3000,1704,1714,10,0,0,0 !17#吊杆下节点ngen,2,6000,104,114,10,0,0,z(1) +ygd(1) !1#吊杆节点ngen,2,6000,204,214,10,0,0,z(2) +ygd(2) !2#吊杆节点ngen,2,6000,304,314,10,0,0,z(3) +ygd(3) !3#吊杆节点ngen,2,6000,404,414,10,0,0,z(4) +ygd(4) !4#吊杆节点ngen,2,6000,504,514,10,0,0,z(5) +ygd(5) !5#吊杆节点ngen,2,6000,604,614,10,0,0,z(6) +ygd(6) !6#吊杆节点ngen,2,6000,704,714,10,0,0,z(7) +ygd(7) !7#吊杆节点ngen,2,6000,804,814,10,0,0,z(8) +ygd(8) !8#吊杆节点ngen,2,6000,904,914,10,0,0,z(9) +ygd(9) !9#吊杆节点ngen,2,6000,1004,1014,10,0,0,z(8)+ygd(10) !10#吊杆节点ngen,2,6000,1104,1114,10,0,0,z(7)+ygd(11) !11#吊杆节点ngen,2,6000,1204,1214,10,0,0,z(6)+ygd(12) !12#吊杆节点ngen,2,6000,1304,1314,10,0,0,z(5)+ygd(13) !13#吊杆节点ngen,2,6000,1404,1414,10,0,0,z(4)+ygd(14) !14#吊杆节点ngen,2,6000,1504,1514,10,0,0,z(3)+ygd(15) !15#吊杆节点ngen,2,6000,1604,1614,10,0,0,z(2)+ygd(16) !16#吊杆节点ngen,2,6000,1704,1714,10,0,0,z(1)+ygd(17) !17#吊杆节点!定义吊杆的参数（每根吊杆都采用109根直径为7毫米的高强钢丝）dgarea=109*3.1415926*0.25*0.007*0.007 !吊杆面积(m2)dgyl=0.3*1670*1e6 !吊杆初应力取用0.3倍的抗拉强度(N/m2) dgl=dgyl*dgarea !吊杆力(N)edg=1.95e11 !吊杆的弹性模量(N/m2)*dim,L0,array,17 !定义吊杆的初始0应力长度数组*dim,L,array,17 !定义吊杆的受力后长度数组*dim,deltL,array,17 !定义吊杆的受力后长度增长量数组*dim,istrn,array,17 !定义吊杆的应变数组*dim,fdg,array,17 !定义在吊杆上施加的力的数组*do,i,1,9,1L0(i)=z(i) !0应力时单元长度L0(18-i)=z(i) !吊杆是沿中轴对称的*enddofdg(1)=dgl !1#吊杆施加的预应力fdg(2)=dgl !2#吊杆施加的预应力fdg(3)=dgl !3#吊杆施加的预应力fdg(4)=dgl !4#吊杆施加的预应力fdg(5)=dgl !5#吊杆施加的预应力fdg(6)=dgl !6#吊杆施加的预应力fdg(7)=dgl !7#吊杆施加的预应力fdg(8)=dgl !8#吊杆施加的预应力fdg(9)=dgl !9#吊杆施加的预应力fdg(10)=dgl !10#吊杆施加的预应力fdg(11)=dgl !11#吊杆施加的预应力fdg(12)=dgl !12#吊杆施加的预应力fdg(13)=dgl !13#吊杆施加的预应力fdg(14)=dgl !14#吊杆施加的预应力fdg(15)=dgl !15#吊杆施加的预应力fdg(16)=dgl !16#吊杆施加的预应力fdg(17)=dgl !17#吊杆施加的预应力type,6mat,6*do,i,1,17,1deltL(i)=Fdg(i)*L0(i)/(edg*dgarea) !吊杆受力后的增长量L(i)=L0(i)+deltL(i)istrn(i)=deltl(i)/L(i) !吊杆受力后的应变r,i+1000,dgarea,istrn(i) !定义吊杆的实常数特性*enddo*do,i,1,17,1real,i+1000 !赋予吊杆的实常数特性e,i*100+4+3000,i*100+4+6000 !生成吊杆单元e,i*100+10+4+3000,i*100+10+4+6000*enddocm,dg,elemesel,noneallsel!************************************************************** ********* r,1018real,1018 !定义一个全为0的实常数便于利用截面特性时，不产生冲突!************************************************************** *********!生成风撑上的节点n,20001,0,21.7,z(11) !风撑上节点坐标n,20002,0,28.5,z(12)n,20003,0,42.1,z(13)n,20004,0,48.9,z(14)n,20005,0,59.1,z(15)n,20006,0,65.9,z(15)n,20007,0,76.1,z(14)n,20008,0,82.9,z(13)n,20009,0,96.5,z(12)n,20010,0,103.3,z(11)n,20011,12.5,21.7,z(11)n,20012,12.5,28.5,z(12)n,20013,12.5,42.1,z(13)n,20014,12.5,48.9,z(14)n,20015,12.5,59.1,z(15)+ygd(18)n,20016,12.5,65.9,z(15)+ygd(19)n,20017,12.5,76.1,z(14)n,20018,12.5,82.9,z(13)n,20019,12.5,96.5,z(12)n,20020,12.5,103.3,z(11)n,20021,-12.5,21.7,z(11)n,20022,-12.5,28.5,z(12)n,20023,-12.5,42.1,z(13)n,20024,-12.5,48.9,z(14)n,20025,-12.5,59.1,z(15)+ygd(18)n,20026,-12.5,65.9,z(15)+ygd(19)n,20027,-12.5,76.1,z(14)n,20028,-12.5,82.9,z(13)n,20029,-12.5,96.5,z(12)n,20030,-12.5,103.3,z(11)n,20000,0,62.5,z(15) !(中撑中点)!风撑1（倒K撑）type,5mat,5SECTYPE,500, beam,ctube,windc !定义风撑的形状为空心圆secdata,0.486,0.50 !赋予截面形状参数,!直径按1.0米计算厚度为1.4厘米secnum,500 !调用截面!生成风撑单元!风撑1（倒K撑）e,20002,10404,90000e,20002,10414,90000e,20002,10304,90000e,20002,10314,90000!风撑2（倒K撑）e,20004,10704,90000e,20004,10714,90000e,20004,10604,90000e,20004,10614,90000!风撑3（中撑X撑）e,20000,20015,90000e,20000,20016,90000e,20000,20025,90000e,20000,20026,90000!风撑4（k撑）e,20007,11104,90000e,20007,11114,90000e,20007,11204,90000e,20007,11214,90000!风撑5（k撑）e,20009,11404,90000e,20009,11414,90000e,20009,11504,90000e,20009,11514,90000cm,fc,elemesel,none!************************************************************** *********!连接拱轴线!定义拱轴上的坐标点n,30000,12.5,0,0 !左拱轴起点n,30001,12.5,125,0 !左拱轴端点n,30002,12.5,2.025,z(16) !左横梁起点n,30003,12.5,125-2.025,z(16) !左横梁端点n,40000,-12.5,0,0 !右拱轴起点n,40001,-12.5,125,0 !右拱轴端点n,40002,-12.5,2.025,z(16) !右横梁起点n,40003,-12.5,125-2.025,z(16) !右横梁端点*dim,zgzx,array,23 !定义拱轴线,左轴线*dim,ygzx,array,23 !定义拱轴线,右轴线zgzx(1)= 30000 !以下是拱轴线上的节点号zgzx(2)= 30002zgzx(3)= 10104zgzx(4)= 10204zgzx(5)= 10304zgzx(6)= 10404zgzx(7)= 10504zgzx(8)= 10604zgzx(9)= 10704zgzx(10)=10804zgzx(11)=20015zgzx(12)=10904zgzx(13)=20016zgzx(14)=11004zgzx(15)=11104zgzx(16)=11204zgzx(17)=11304zgzx(18)=11404zgzx(19)=11504zgzx(20)=11604zgzx(21)=11704zgzx(22)=30003zgzx(23)=30001ygzx(1)= 40000ygzx(2)= 40002ygzx(3)= 10114ygzx(4)= 10214ygzx(5)= 10314ygzx(6)= 10414ygzx(7)= 10514ygzx(8)= 10614ygzx(9)= 10714ygzx(10)=10814ygzx(11)=20025ygzx(12)=10914ygzx(13)=20026ygzx(14)=11014ygzx(15)=11114ygzx(16)=11214ygzx(17)=11314ygzx(18)=11414ygzx(19)=11514ygzx(20)=11614ygzx(21)=11714ygzx(22)=40003ygzx(23)=40001n,90001,12.5,0,30 !钢管左轴线参考点n,90002,-12.5,0,30 !钢管左轴线参考点type,1mat,1r,2000secnum,1*do,i,1,22 !生成拱轴单元（钢管）e,zgzx(i),zgzx(i+1),90001e,ygzx(i),ygzx(i+1),90002*enddocm,gzgg,elem !生成组esel,nonetype,2r,2001secnum,2*do,i,1,22 !生成拱轴单元（混凝土）e,zgzx(i),zgzx(i+1),90001e,ygzx(i),ygzx(i+1),90002*enddocm,gzhnt,elem !生成组esel,none!************************************************************** ***************!建立拱轴底座type,9 !30#混凝土r,6real,6mat,9SECTYPE,70, beam, rect,dizuo !定义截面形状secdata,1.4,4 !赋予截面形状参数secnum,70 !调用截面e,30001,30003,90001 !拱轴底座e,40001,40003,90002e,30000,30002,90001e,40000,40002,90002cm,gzdz,elem !生成组esel,none!************************************************************** ***************!建立盖梁ngen,2,2000,1,16,1,0,-2.025,0 !复制节点到拱轴线起点位置ngen,2,100,1801,1816,1,0,2.025,0 !复制节点到拱轴线终点位置mat,7r,7real,7SECTYPE,60, beam, rect,gail !定义截面形状secdata,4,1.6 !赋予截面形状参数secnum,60 !调用截面n,90005,0,0,20 !定义参考点n,90006,0,125,20e,1901,1902,90006 !盖梁单元125米处e,1902,1903,90006e,1903,1904,90006e,1904,1905,90006e,1905,1906,90006e,1901,1912,90006e,1912,1913,90006e,1913,1914,90006e,1914,1915,90006e,1915,1916,90006e,2001,2002,90005 !盖梁单元0米处e,2002,2003,90005e,2003,2004,90005e,2004,2005,90005e,2005,2006,90005e,2001,2012,90005e,2012,2013,90005e,2013,2014,90005e,2014,2015,90005e,2015,2016,90005cm,gailiang,elem !生成组esel,none!************************************************************** *********!建立墩柱n,90003,14.5,0,0 !参考点n,90004,-14.5,0,0type,8mat,8r,8real,8SECTYPE,50, beam,hrec,dun !定义截面形状secdata,3,4,0.8,0.8,0.8,0.8 !赋予截面形状参数ngen,2,10000,2004,2014,10,0,0,-10 !墩节点距离桥面-10米ngen,2,10000,1904,1914,10,0,0,-10secnum,50e,2004,12004,90003e,2014,12014,90004e,1904,11904,90003e,1914,11914,90004cm,dunzhu,elem !生成组esel,none!************************************************************** ************!承台单元type,11mat,11r,1SECTYPE,40, beam,rect,chengtaisecdata,12.7,9r,11real, 11secnum,40ngen,2,20000,2004,2014,10,0,0,-11.5ngen,2,20000,1904,1914,10,0,0,-11.5ngen,2,30000,2004,2014,10,0,0,-13ngen,2,30000,1904,1914,10,0,0,-13e,12004,22004,90003 !生成承台单元e,11904,21904,90003e,12014,22014,90004e,11914,21914,90004e,22004,32004,90003e,21904,31904,90003e,22014,32014,90004e,21914,31914,90004cm,chtai,elem !生成组esel,none!************************************************************** ******!横系梁type,9mat,9r,9real,9n,90007,12.5,0,10 !参考点n,90008,-12.5,0,10SECTYPE,90,beam, rect,hxl !定义截面形状secdata,3,3 !赋予截面形状参数secnum,90 !调用截面e,22004,22014,90007e,21904,21914,90008cm,hxl,elem !生成组esel,none!!************************************************************** *************!桩n,50000,6.15+1.85,2.7,-13 !定义桩节点n,50001,6.15+1.85,-2.7,-13n,50002,6.15+1.85+4.5,2.7,-13n,50003,6.15+1.85+4.5,-2.7,-13n,50004,6.15+12.7-1.85,2.7,-13n,50005,6.15+12.7-1.85,-2.7,-13nsel,s,,,50000,50005,1nsym,x,10,50000,50005,1allselngen,2,100,50000,50015,1,0,125,0 !首先生成平面内的24个节点! 每3米生成一层节点(共8层)ngen,10,1000,50000,51000,1,0,0,-2.5!生成桩单元esel,nonetype,11mat,11r,11real,11SECTYPE,20, beam,csolid,piersecdata,0.9secnum,20*do,i,0,8,1*do,k,0,5,1e,50000+k+i*1000,50000+k+(i+1)*1000 ,90003e,50000+k+10+i*1000,50000+k+10+(i+1)*1000 ,90003e,50000+k+100+i*1000,50000+k+100+(i+1)*1000,90003 e,50000+k+110+i*1000,50000+k+110+(i+1)*1000,90003 *enddo*enddocm,zhuang,elem !生成组esel,none!************************************************************** ***!在桩上施加弹簧单元!第一层距离承台-3米 317-340!计算弹簧的弹性系数 m的具体计算方法参照本书斜拉桥建模部分d=1.8kf=0.9k0=1+1/db1=0.6 !和每排桩数有关的系数,当为2排时,为0.6L1=3.6 !桩间净距h1=3*(d+1) !桩在地面或冲刷线下的计算深k=b1+(1-b1)/0.6*L1/h1B=Kf*k0*k*d !换算宽度!第1层土粉质量粘土 m=8 厚度5!第2层土粘土（细砂） m=15 厚度5 !第3层土中砂 m=20 厚度5!第4层土粘土 m=25 厚度5!第5层基岩 c0=300 厚度10米*dim,m,array,5m(1)=8m(2)=15m(3)=20m(4)=25m(5)=4000 !基岩系数*dim,h,array,5h(1)=5 !第一层土的厚度h(2)=5h(3)=5h(4)=5h(5)=10*dim,c,array,4c(1)=m(1)*h(1)c(2)=m(2)*(h(1)+h(2))c(3)=m(3)*(h(1)+h(2)+h(3))c(4)=m(4)*(h(1)+h(2)+h(3)+h(4)) *dim,e,array,4 !弹簧弹性系数e(1)=0.5*c(1)*h(1)*B*1e6*do,i,2,4,1e(i)=0.5*(c(i-1)+c(i))*h(i)*B*1e6*enddoet,20,combin14 !弹簧单元r,20,e(1) !第1层土的弹性系数r,30,e(2) !第2层土的弹性系数r,40,e(3) !第3层土的弹性系数r,50,e(4) !第4层土的弹性系数r,60,m(5)*3.1415926*B*B/4*1e6 !基础岩石弹性系数!第1层底部距离承台-5米ngen,2,10000,51000,51115,1,1,0,0 !x 方向ngen,2,20000,51000,51115,1,0,1,0 !y 方向!ngen,2,30000,51000,51115,1,0,0,1 !z 方向type,20real,20*do,i,0,5,1e,51000+i,51000+i+10000e,51000+i+10,51000+i+10+10000e,51100+i,51100+i+10000 ! xe,51100+i+10,51100+i+10+10000e,51000+i,51000+i+20000e,51000+i+10,51000+i+10+20000e,51100+i,51100+i+20000 !ye,51100+i+10,51100+i+10+20000*enddocm,spring1,elemesel,none!第2层底部距离承台-10米ngen,2,10000,53000,53115,1,1,0,0!x 方向ngen,2,20000,53000,53115,1,0,1,0!y 方向!ngen,2,30000,53000,53115,1,0,0,1!z 方向type,20real,30*do,i,0,5,1e,53000+i+10,53000+i+10+10000e,53000+i,53000+i+10000e,53100+i,53100+i+10000 ! xe,53100+i+10,53100+i+10+10000e,53000+i,53000+i+20000e,53000+i+10,53000+i+10+20000e,53100+i,53100+i+20000e,53100+i+10,53100+i+10+20000 !y *enddocm,spring2,elemesel,none!第3层距离承台-15米ngen,2,10000,55000,55115,1,1,0,0 !x 方向ngen,2,20000,55000,55115,1,0,1,0 !y 方向type,20real,40*do,i,0,5,1e,55000+i,55000+i+10000e,55000+i+10,55000+i+10+10000e,55100+i,55100+i+10000 ! xe,55100+i+10,55100+i+10+10000e,55000+i,55000+i+20000e,55000+i+10,55000+i+10+20000e,55100+i,55100+i+20000 !ye,55100+i+10,55100+i+10+20000*enddocm,spring3,elemesel,none!!第4层距离承台-20米ngen,2,10000,57000,57115,1,1,0,0 !x 方向ngen,2,20000,57000,57115,1,0,1,0 !y 方向type,20real,50*do,i,0,5,1e,57000+i+10,57000+i+10+10000e,57100+i,57100+i+10000 ! xe,57100+i+10,57100+i+10+10000e,57000+i,57000+i+10000e,57000+i,57000+i+20000e,57000+i+10,57000+i+10+20000e,57100+i,57100+i+20000 !ye,57100+i+10,57100+i+10+20000*enddocm,spring4,elemesel,none!!第5层距离承台-22.5米ngen,2,10000,58000,58115,1,1,0,0 !x 方向ngen,2,20000,58000,58115,1,0,1,0 !y 方向type,20real,60*do,i,0,5,1e,58000+i+10,58000+i+10+10000e,58000+i,58000+i+10000e,58100+i,58100+i+10000 ! xe,58100+i+10,58100+i+10+10000e,58000+i,58000+i+20000e,58000+i+10,58000+i+10+20000e,58100+i,58100+i+20000 !ye,58100+i+10,58100+i+10+20000*enddocm,spring5,elemesel,none!!第6层距离承台-25米ngen,2,10000,59000,59115,1,1,0,0 ngen,2,10000,59000,59115,1,1,0,0 !x 方向ngen,2,20000,59000,59115,1,0,1,0 ngen,2,20000,59000,59115,1,0,1,0 !y 方向ngen,2,30000,59000,59115,1,0,0,-1 ngen,2,30000,59000,59115,1,0,0,-1 !z 方向type,20real,60*do,i,0,5,1e,59000+i,59000+i+10000e,59000+i+10,59000+i+10+10000e,59100+i,59100+i+10000 ! xe,59100+i+10,59100+i+10+10000e,59000+i,59000+i+20000e,59000+i+10,59000+i+10+20000e,59100+i,59100+i+20000 !ye,59100+i+10,59100+i+10+20000e,59000+i,59000+i+30000e,59000+i+10,59000+i+10+30000e,59100+i,59100+i+30000 !ze,59100+i+10,59100+i+10+30000*enddocm,spring6,elemallsel!************************************************************** *********!剔除不需要进行对称复制的单元allselnsel,s,loc,z,-50,-0.01nsel,r,loc,y,-10,10cm,zdnode,nodeesln,s,,allcm,zdelem,elemallselnsel,u,,,zdnodensym,y,100000,all !节点对称allselESEL,U,,,zdelemESYM,,100000,ALL !见图11-5ALLSEL!图11-5 全桥单元模型图!************************************************************** ********* !定义系杆单元（全桥通长）!系杆单元!定义系杆的参数（每根系杆都采用8束27*7*d5高强钢丝）xgarea=8*27*7*3.1415926*0.25*0.005*0.005 !系杆集束面积exg=1.95e11*dim,fxg,array,2*dim,Llast,array,2*dim,Listrn,array,2*dim,detL,array,2xgyl=0.3*1860*1e6xgl=xgyl*xgareafxg(1)=xgl !施加的预应力fxg(2)=xgl !施加的预应力Lzero=250 !0应力时单元长度*do,i,1,2,1detL(i)=Lzero*Fxg(i)/(exg*xgarea)Llast(i)=Lzero+detL(i) !单元最终长度Listrn(i)=detl(i)/Llast(i)r,i+100,xgarea,Listrn(i)*enddotype,10mat,10real,101e,30001,130001 !系杆1 250米real,102e,40001,140001 !系杆2 250米.cm,xigan,elemesel,noneallsel!************************************************************** *******!将盖梁和拱轴交界点进行节点耦合,放松uy,和rotxcp,1,ux,2004,30000cp,1,uy,2004,30000cp,next,uz,2004,30000cp,next,roty,2004,30000 cp,next,rotz,2004,30000 cp,next,ux,1904,30001cp,next,uy,1904,30001cp,next,uz,1904,30001cp,next,roty,1904,30001 cp,next,rotz,1904,30001 cp,next,ux,2014,40000cp,next,uy,2014,40000cp,next,uz,2014,40000cp,next,roty,2014,40000 cp,next,rotz,2014,40000 cp,next,ux,1914,40001cp,next,uy,1914,40001cp,next,uz,1914,40001cp,next,roty,1914,40001 cp,next,rotz,1914,40001cp,next,ux,101914,140001 cp,next,uy,101914,140001 cp,next,uz,101914,140001 cp,next,roty,101914,140001 cp,next,rotz,101914,140001 cp,next,ux,101904,130001 cp,next,uy,101904,130001 cp,next,uz,101904,130001 cp,next,roty,101904,130001 cp,next,rotz,101904,130001 allselnsel,s,,,30000,130000,100000nsel,a,,,40000,140000,100000nummrg,nodeallsel!************************************************************* ********do,i,1,3,1cp,next,all,1800+i,1900+i !盖梁和0#横梁耦合*enddoallsel*do,i,5,6,1cp,next,all,1800+i,1900+i*enddoallsel*do,i,1,2,1cp,next,all,1800+i+11,1900+i+11cp,next,all,i+11,2000+i+11,i+100000+11*enddoallsel*do,i,4,5,1cp,next,all,1800+i+11,1900+i+11cp,next,all,i+11,2000+i+11,i+100000+11*enddoallsel*do,i,1,3,1cp,next,all,i,2000+i,100000+i*enddoallsel*do,i,5,6,1cp,next,all,i,2000+i,100000+i*enddoallsel*do,i,1,3,1cp,next,all,101800+i,101900+i !盖梁和0#横梁耦合*enddoallsel*do,i,5,6,1cp,next,all,101800+i,101900+iallsel*enddo*do,i,1,2cp,next,all,101800+i+11,101900+i+11*enddoallsel*do,i,4,5,1cp,next,all,101800+i+11,101900+i+11*enddoallsel!************************************************************** ********!!桩节点和承台节点耦合nsel,s ,loc,z,-13nsel,r ,loc,y,-3,3nsel,r ,loc,x,8,20cp,next,all,allallselnsel,s ,loc,z,-13nsel,r ,loc,y,-3,3nsel,r ,loc,x,-8,-20cp,next,all,allallselnsel,s ,loc,z,-13nsel,r ,loc,y,122,128nsel,r ,loc,x,-8,-20cp,next,all,allallselnsel,s ,loc,z,-13nsel,r ,loc,y,122,128nsel,r ,loc,x,8,20cp,next,all,allallselnsel,s ,loc,z,-13nsel,r ,loc,y,-122,-128nsel,r ,loc,x,8,20cp,next,all,allallselnsel,s ,loc,z,-13nsel,r ,loc,y,-122,-128nsel,r ,loc,x,-8,-20cp,next,all,allallsel!************************************************************** ********!将吊杆转动进行放松*do,i,1,17,1cp,next,ux,i*100+4,i*100+4+3000cp,next,ux,i*100+4+6000,i*100+4+10000cp,next,uy,i*100+4,i*100+4+3000cp,next,uy,i*100+4+6000,i*100+4+10000cp,next,uz,i*100+4,i*100+4+3000cp,next,uz,i*100+4+6000,i*100+4+10000cp,next,ux,i*100+4+10,i*100+4+3000+10cp,next,ux,i*100+4+6000+10,i*100+4+10000+10cp,next,uy,i*100+4+10,i*100+4+3000+10cp,next,uy,i*100+4+6000+10,i*100+4+10000+10cp,next,uz,i*100+4+10,i*100+4+3000+10cp,next,uz,i*100+4+6000+10,i*100+4+10000+10cp,next,ux,i*100+4+100000,i*100+4+3000+100000cp,next,ux,i*100+4+6000+100000,i*100+4+10000+100000 cp,next,uy,i*100+4+100000,i*100+4+3000+100000cp,next,uy,i*100+4+6000+100000,i*100+4+10000+100000 cp,next,uz,i*100+4+100000,i*100+4+3000+100000cp,next,uz,i*100+4+6000+100000,i*100+4+10000+100000 cp,next,ux,i*100+4+100000+10,i*100+4+3000+100000+10 cp,next,ux,i*100+4+6000+100000+10,i*100+4+10000+100 000+10cp,next,uy,i*100+4+100000+10,i*100+4+3000+100000+10cp,next,uy,i*100+4+6000+100000+10,i*100+4+10000+100 000+10cp,next,uz,i*100+4+100000+10,i*100+4+3000+100000+10 cp,next,uz,i*100+4+6000+100000+10,i*100+4+10000+100 000+10*enddoallsel!************************************************************** ********!将土弹簧的一端约束allsel!nsel,u,,,erasensel,s,loc,x,13.5nsel,a,loc,x,9nsel,a,loc,x,18nsel,a,loc,x,-11.5nsel,a,loc,x,-7nsel,a,loc,x,-16nsel,a,loc,y,3.7nsel,a,loc,y,-1.7nsel,a,loc,y,128.7nsel,a,loc,y,123.3nsel,a,loc,y,-128.7nsel,a,loc,y,-123.3nsel,r,loc,z,-40,-13d,all,allallselnsel,s,loc,z,-36.5d,all,all/color,pbak,off/replotallselnsel,s,loc,y,0nsel,r,loc,z,0nummrg,node,allesel,noneallsel!************************************************************** ********!工况1 不包括桥面板时，求横梁和纵梁的总重量/soluallselmpdele,dens,allmp,dens,3,2600mp,dens,4,2600acel,,,9.8pstres,onsolvesave! 总重量：0.65064E+08 N!************************************************************** *********!工况2 恒载下调节预拱度首先只有拱肋自重作用下的拱的变位/soluallselmpdele,dens,3,11,1mp,dens,1,7800mp,dens,2,2600*do,i,1,17,1r,1000+i,dgarea,0*enddor,101,xgarea,0r,102,xgarea,0acel,,,9.8solve!************************************************************** *********!工况3一期恒载（包括桥面板）作用下的无预应力状态/soluallsel*do,i,1,17,1。
ansys混凝土命令流
vsel,s,loc,x,-0.002,0 !选择左边建筑结构胶
vatt,2,2,1
esize,0.01
mshape,0,3d$mshkey,1
vmesh,all
vsel,s,loc,x,0.102,0.112 !选择右边钢板
vatt,3,3,1
wpoff,0.002,, !工作平面平移
blc4,,,0.01,0.05 ,0.16 !钢板
wpcsys,-1
wpoff,-0.002,0.025,0.05 !工作平面平移
blc4,,,0.002,0.05 ,0.15 !建筑结构胶
blc4,,,-0.01,0.05 ,0.16 !钢板
fini
/cle
/FILNAME, hntd, 1
/TITLE, tiegang
/prep7
!----------定义单元及材料等---------------------
Et,1,solid45 !定义混凝土单元类型
Et,2,solid45 !定义建筑结构胶单元类型
mp,prxy,1,0.2 !定义混凝土泊松比
mp,prxy,2,0.317 !定义建筑结构胶泊松比
mp,prxy,3,0.3 !定义钢板泊松比
mp,dens,1,2550 !定义混凝土密度
mp,dens,2,1700 !定义建筑结构胶密度
sfa,all,,pres,-1200e3/(2*6)
save
!-----------求解---------------------------
/solu
/nerr,10000,100000
solve
钢筋混凝土梁ansys分析附命令流
钢筋混凝⼟梁ansys分析附命令流钢筋混凝⼟⾮线性分析2015⼤作业上海交通⼤学陈明1、参数选择梁的截⾯宽度为200mm，上部配置2Φ8受压筋，混凝⼟的净保护层厚度为25 mm（从纵向钢筋外边缘算起），箍筋两端区采⽤8@100的双肢箍，中间区取8@200 双肢箍1）梁的截⾯⾼度选300mm；2）两加载间的距离选1000mm；3）混凝⼟选C30；4）纵向受拉钢筋配筋选218；2、描述选⽤的有限元模型及单元的特点采⽤ansys软件进⾏模拟计算，钢筋混凝⼟模型采⽤分离式模型，不考虑钢筋与混凝⼟之间的相对滑移。
混凝⼟采⽤solid65单元模拟，solid65⽤于模拟三维有钢筋或⽆钢筋的混凝⼟模型。
该单元能够计算拉裂和压碎。
在混凝⼟应⽤中，该单元的实体功能可以⽤于建⽴混凝⼟模型，同时，还可⽤加筋功能建⽴钢筋混凝⼟模型。
另外，该单元还可以应⽤于加强复合物和地质材料。
该单元由⼋个节点定义，每个节点有三个⾃由度:节点坐标系的x，y，z⽅向的平动。
⾄多可以定义三种不同规格的钢筋。
钢筋单元采⽤link180单元模拟，link180是⼀个适⽤于各类⼯程应⽤的三维杆单元。
根据具体情况，该单元可以被看作桁架单元、索单元、链杆单元或弹簧单元等等。
本单元是⼀个轴向拉伸⼀压缩单元，每个节点有三个⾃由度:节点坐标系的x，y，z⽅向的平动。
本单元是⼀种顶端铰接结构，不考虑单元弯曲。
本单元具有塑性、蠕变、旋转、⼤变形和⼤应变功能。
缺省时，当考虑⼤变形时任何分析中LINK180单元都包括应⼒刚化选项。
3、描述选⽤的混凝⼟与钢筋粘结滑移本构关系的具体形式、参数等。
钢筋的应⼒应变关系曲线考虑到极限塑性应变最⼤值为0.01,钢筋本构模型采⽤多线性模型kinh，初始弹性模量为Es=200000Mpa，强化系数为0.001。
混凝⼟的应⼒应变关系曲线混凝⼟选⽤各向同性的miso模型，当计⼊下降端时，程序报错，所以只取了前⾯的上升段，⽤5段折线模拟混凝⼟应⼒应变曲线。
ansys混凝土浇筑命令流
ansys混凝土浇筑命令流/CLEAR/FILNAME,CT,1/TITLE,CT01*SET,T1,1 !第一层开始浇筑时间*SET,T2,8 !第二层开始浇筑时间*SET,TZZ,30 !终止时间*SET,cdiam,0.028 !水管径*SET,carea,3.14*cdiam*cdiam/4 !水管面积*SET,massflow,28800 !管内流体质量流率（kg/d）*SET,KX1,234 !混凝土导热系数*SET,C1,0.963 !混凝土比热*SET,DENSE1,2400 !砼密度*SET,KX2,50.1 !水导热系数*SET,C2,4.2 !水比热*SET,DENSE2,1000 !水密度*SET,DLD,588 !上层表面放热系数$*SET,DLB,994 !底层放热系数$*SET,DLC,1983 !侧面放热系数*SET,Q_SUM,326*377 !单位水化热总量*SET,TEMP_q,19 !大气温度*SET,TEMP_c,22 !砼温度*SET,TEMP_w,17 !水流温度/PREP7ET,1,70 !定义单元solid70et,2,116,1,2 !单元fluid116MP,KXX,1,KX1 !单元70 材料MP,C,1,C1MP,DENS,1,DENSE1MP,KXX,2,KX2 !单元116材料MP,C,2,C2MP,DENS,2,DENSE2MP,HF,2,3574r,1,cdiam,carea,1 !单元116实常数!建立实体模型K,1,0,0,0$K,2,12.19,0,0$K,3,12.19,0,8.25$K,4,0,0,8.25K,5,12.65,0,8.25K,6,15.15,0,0$A,1,2,3,4A,2,3,5,6$ASEL,ALL$VEXT,ALL,,,,5 !通过偏移量由面生成体ALLSEL,ALL !选择所有的实体VGLUE,ALL !体粘接NUMMRG,ALL !对实体进行合并$NUMCMP,ALL !压缩所定义项的编号TYPE, 1 !划分模型 !指定单元类型属性指示器FLST,5,4,4,ORDE,4FITEM,5,1FITEM,5,3FITEM,5,8FITEM,5,10CM,_Y,LINE !由所选择的实体生成一个元件LSEL, , , ,P51X !p51x说包含你用鼠标选取的所有线CM,_Y1,LINECMSEL,,_Y !选择一组元件或部件子集LESIZE,_Y1, , ,23, , , , ,1 !对选择的线设置网格单元大小23为分段数1修改所有已选择的线FLST,5,6,4,ORDE,6 !为某一命令选择操作实体时产生的FITEM,5,2FITEM,5,4FITEM,5,6FITEM,5,11FITEM,5,17CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINE CMSEL,,_YLESIZE,_Y1, , ,33, , , , ,1 FLST,5,6,4,ORDE,4 FITEM,5,12FITEM,5,-15FITEM,5,19FITEM,5,-20CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINE CMSEL,,_YLESIZE,_Y1, , ,10, , , , ,1 FLST,5,4,4,ORDE,4 FITEM,5,5FITEM,5,7FITEM,5,16FITEM,5,18CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINE CMSEL,,_YLESIZE,_Y1, , ,5, , , , ,1 FLST,5,2,6,ORDE,2 !6为实体FITEM,5,-2CM,_Y,VOLUVSEL, , , ,P51XCM,_Y1,VOLUCHKMSH,'VOLU' !用来检查网格质量的，可查看你画的网格是否有警告单元，或err单元CMSEL,S,_YVSWEEP,_Y1 !利用与体邻近的面单元采用扫琼方式对体进行网格划分CMDELE,_Y !删除一个现存的元件或部件CMDELE,_Y1CMDELE,_Y2NSEL,S,LOC,Y,1 !选择第一层水管节点FLST,5,242,1,ORDE,212FITEM,5,1916FITEM,5,2174FITEM,5,2184FITEM,5,2200FITEM,5,2209FITEM,5,2218FITEM,5,2227FITEM,5,2236FITEM,5,2245FITEM,5,2254FITEM,5,2263FITEM,5,2272FITEM,5,2281FITEM,5,2299 FITEM,5,2308 FITEM,5,2317 FITEM,5,2326 FITEM,5,2335 FITEM,5,2344 FITEM,5,2353 FITEM,5,2362 FITEM,5,2371 FITEM,5,2380 FITEM,5,2389 FITEM,5,2414 FITEM,5,2429 FITEM,5,2438 FITEM,5,2447 FITEM,5,2456 FITEM,5,2465 FITEM,5,2474 FITEM,5,2483 FITEM,5,2492 FITEM,5,2501 FITEM,5,2510 FITEM,5,2519 FITEM,5,2528 FITEM,5,2537 FITEM,5,2546 FITEM,5,2555 FITEM,5,2564 FITEM,5,2573FITEM,5,2591 FITEM,5,2600 FITEM,5,2609 FITEM,5,2618 FITEM,5,2643 FITEM,5,2658 FITEM,5,2946 FITEM,5,2955 FITEM,5,2964 FITEM,5,2973 FITEM,5,2981 FITEM,5,3021 FITEM,5,-3052 FITEM,5,3311 FITEM,5,3326 FITEM,5,3335 FITEM,5,3344 FITEM,5,3353 FITEM,5,3360 FITEM,5,3369 FITEM,5,3378 FITEM,5,3387 FITEM,5,3396 FITEM,5,3405 FITEM,5,3414 FITEM,5,3423 FITEM,5,3432 FITEM,5,3441 FITEM,5,3450FITEM,5,3468 FITEM,5,3477 FITEM,5,3486 FITEM,5,3495 FITEM,5,3504 FITEM,5,3513 FITEM,5,3522 FITEM,5,3531 FITEM,5,3540 FITEM,5,3549 FITEM,5,3558 FITEM,5,3567 FITEM,5,3576 FITEM,5,3585 FITEM,5,3594 FITEM,5,3603 FITEM,5,3612 FITEM,5,3621 FITEM,5,3630 FITEM,5,3639 FITEM,5,3648 FITEM,5,3657 FITEM,5,3666 FITEM,5,3675 FITEM,5,3684 FITEM,5,3693 FITEM,5,3702 FITEM,5,3711 FITEM,5,3720FITEM,5,3738 FITEM,5,3747 FITEM,5,3756 FITEM,5,3765 FITEM,5,3774 FITEM,5,3783 FITEM,5,3792 FITEM,5,3801 FITEM,5,3810 FITEM,5,3819 FITEM,5,3828 FITEM,5,3837 FITEM,5,3846 FITEM,5,3855 FITEM,5,3864 FITEM,5,3873 FITEM,5,3882 FITEM,5,3891 FITEM,5,3900 FITEM,5,3909 FITEM,5,3918 FITEM,5,3927 FITEM,5,3936 FITEM,5,3945 FITEM,5,3954 FITEM,5,3963 FITEM,5,3972 FITEM,5,3981 FITEM,5,3990FITEM,5,4008 FITEM,5,4017 FITEM,5,4026 FITEM,5,4035 FITEM,5,4044 FITEM,5,4053 FITEM,5,4062 FITEM,5,4071 FITEM,5,4080 FITEM,5,4089 FITEM,5,4098 FITEM,5,4107 FITEM,5,4116 FITEM,5,4125 FITEM,5,4134 FITEM,5,4143 FITEM,5,4152 FITEM,5,4161 FITEM,5,4170 FITEM,5,4179 FITEM,5,4188 FITEM,5,4197 FITEM,5,4206 FITEM,5,4215 FITEM,5,4224 FITEM,5,4233 FITEM,5,4242 FITEM,5,4251 FITEM,5,4260FITEM,5,4278 FITEM,5,4287 FITEM,5,4296 FITEM,5,4305 FITEM,5,4314 FITEM,5,4323 FITEM,5,4332 FITEM,5,4341 FITEM,5,4350 FITEM,5,4359 FITEM,5,4368 FITEM,5,4377 FITEM,5,4386 FITEM,5,4395 FITEM,5,4404 FITEM,5,4413 FITEM,5,4422 FITEM,5,4431 FITEM,5,4440 FITEM,5,4449 FITEM,5,4458 FITEM,5,4467 FITEM,5,4476 FITEM,5,4485 FITEM,5,4494 FITEM,5,4503 FITEM,5,4518 FITEM,5,4806 FITEM,5,5094FITEM,5,5670 FITEM,5,5958 FITEM,5,6246 FITEM,5,6534 FITEM,5,6822 FITEM,5,7110 FITEM,5,7398 FITEM,5,7686 FITEM,5,7974 FITEM,5,8262FITEM,5,8550 FITEM,5,8838 FITEM,5,9126 FITEM,5,9414 FITEM,5,9702 FITEM,5,9990 FITEM,5,10278 FITEM,5,10566 NSEL,U, , ,P51X NPLOT !显示节点编号FLST,5,362,1,ORDE,333 FITEM,5,1885 FITEM,5,-1915 FITEM,5,4527 FITEM,5,4536 FITEM,5,4545 FITEM,5,4554 FITEM,5,4563FITEM,5,4581 FITEM,5,4590 FITEM,5,4599 FITEM,5,4608 FITEM,5,4617 FITEM,5,4626 FITEM,5,4635 FITEM,5,4644 FITEM,5,4653 FITEM,5,4662 FITEM,5,4671 FITEM,5,4680 FITEM,5,4689 FITEM,5,4698 FITEM,5,4707 FITEM,5,4716 FITEM,5,4725 FITEM,5,4734 FITEM,5,4743 FITEM,5,4752 FITEM,5,4761 FITEM,5,4770 FITEM,5,4779 FITEM,5,4788 FITEM,5,5112 FITEM,5,5121 FITEM,5,5130 FITEM,5,5139 FITEM,5,5148FITEM,5,5166 FITEM,5,5175 FITEM,5,5184 FITEM,5,5193 FITEM,5,5202 FITEM,5,5211 FITEM,5,5220 FITEM,5,5229 FITEM,5,5238 FITEM,5,5247 FITEM,5,5256 FITEM,5,5265 FITEM,5,5274 FITEM,5,5283 FITEM,5,5292 FITEM,5,5301 FITEM,5,5310 FITEM,5,5319 FITEM,5,5328 FITEM,5,5337 FITEM,5,5346 FITEM,5,5355 FITEM,5,5364 FITEM,5,5373 FITEM,5,5679 FITEM,5,5688 FITEM,5,5697 FITEM,5,5706 FITEM,5,5715FITEM,5,5733 FITEM,5,5742 FITEM,5,5751 FITEM,5,5760 FITEM,5,5769 FITEM,5,5778 FITEM,5,5787 FITEM,5,5796 FITEM,5,5805 FITEM,5,5814 FITEM,5,5823 FITEM,5,5832 FITEM,5,5841 FITEM,5,5850 FITEM,5,5859 FITEM,5,5868 FITEM,5,5877 FITEM,5,5886 FITEM,5,5895 FITEM,5,5904 FITEM,5,5913 FITEM,5,5922 FITEM,5,5931 FITEM,5,5940 FITEM,5,6264 FITEM,5,6273 FITEM,5,6282 FITEM,5,6291 FITEM,5,6300FITEM,5,6318 FITEM,5,6327 FITEM,5,6336 FITEM,5,6345 FITEM,5,6354 FITEM,5,6363 FITEM,5,6372 FITEM,5,6381 FITEM,5,6390 FITEM,5,6399 FITEM,5,6408 FITEM,5,6417 FITEM,5,6426 FITEM,5,6435 FITEM,5,6444 FITEM,5,6453 FITEM,5,6462 FITEM,5,6471 FITEM,5,6480 FITEM,5,6489 FITEM,5,6498 FITEM,5,6507 FITEM,5,6516 FITEM,5,6525 FITEM,5,6831 FITEM,5,6840 FITEM,5,6849 FITEM,5,6858 FITEM,5,6867FITEM,5,6885 FITEM,5,6894 FITEM,5,6903 FITEM,5,6912 FITEM,5,6921 FITEM,5,6930 FITEM,5,6939 FITEM,5,6948 FITEM,5,6957 FITEM,5,6966 FITEM,5,6975 FITEM,5,6984 FITEM,5,6993 FITEM,5,7002 FITEM,5,7011 FITEM,5,7020 FITEM,5,7029 FITEM,5,7038 FITEM,5,7047 FITEM,5,7056 FITEM,5,7065 FITEM,5,7074 FITEM,5,7083 FITEM,5,7092 FITEM,5,7407 FITEM,5,7416 FITEM,5,7425 FITEM,5,7434 FITEM,5,7443FITEM,5,7461 FITEM,5,7470 FITEM,5,7479 FITEM,5,7488 FITEM,5,7497 FITEM,5,7506 FITEM,5,7515 FITEM,5,7524 FITEM,5,7533 FITEM,5,7542 FITEM,5,7551 FITEM,5,7560 FITEM,5,7569 FITEM,5,7578 FITEM,5,7587 FITEM,5,7596 FITEM,5,7605 FITEM,5,7614 FITEM,5,7623 FITEM,5,7632 FITEM,5,7641 FITEM,5,7650 FITEM,5,7659 FITEM,5,7668 FITEM,5,7677 FITEM,5,7983 FITEM,5,7992 FITEM,5,8001 FITEM,5,8010FITEM,5,8028 FITEM,5,8037 FITEM,5,8046 FITEM,5,8055 FITEM,5,8064 FITEM,5,8073 FITEM,5,8082 FITEM,5,8091 FITEM,5,8100 FITEM,5,8109 FITEM,5, 8118 FITEM,5,8127 FITEM,5,8136 FITEM,5,8145 FITEM,5,8154 FITEM,5,8163 FITEM,5,8172 FITEM,5,8181 FITEM,5,8190 FITEM,5,8199 FITEM,5,8208 FITEM,5,8217 FITEM,5,8226 FITEM,5,8235 FITEM,5,8244 FITEM,5,8568 FITEM,5,8577 FITEM,5,8586FITEM,5,8604 FITEM,5,8613 FITEM,5,8622 FITEM,5,8631 FITEM,5,8640 FITEM,5,8649 FITEM,5,8658 FITEM,5,8667 FITEM,5,8676 FITEM,5,8685 FITEM,5,8694 FITEM,5,8703 FITEM,5,8712 FITEM,5,8721 FITEM,5,8730 FITEM,5,8739 FITEM,5,8748 FITEM,5,8757 FITEM,5,8766 FITEM,5,8775 FITEM,5,8784 FITEM,5,8793 FITEM,5,8802 FITEM,5,8811 FITEM,5,8820 FITEM,5,8829 FITEM,5,9135 FITEM,5,9144 FITEM,5,9153FITEM,5,9171 FITEM,5,9180 FITEM,5,9189 FITEM,5,9198 FITEM,5,9207 FITEM,5,9216 FITEM,5,9225 FITEM,5,9234 FITEM,5,9243 FITEM,5,9252 FITEM,5,9261 FITEM,5,9270 FITEM,5,9279 FITEM,5,9288 FITEM,5,9297 FITEM,5,9306 FITEM,5,9315 FITEM,5,9324 FITEM,5,9333 FITEM,5,9342 FITEM,5,9351 FITEM,5,9360 FITEM,5,9369 FITEM,5,9378 FITEM,5,9387 FITEM,5,9396 FITEM,5,9720 FITEM,5,9729 FITEM,5,9738FITEM,5,9756 FITEM,5,9765 FITEM,5,9774 FITEM,5,9783 FITEM,5,9792 FITEM,5,9801 FITEM,5,9810 FITEM,5,9819 FITEM,5,9828 FITEM,5,9837 FITEM,5,9846 FITEM,5,9855 FITEM,5,9864 FITEM,5,9873 FITEM,5,9882 FITEM,5,9891 FITEM,5,9900 FITEM,5,9909 FITEM,5,9918 FITEM,5,9927 FITEM,5,9936 FITEM,5,9945 FITEM,5,9954 FITEM,5,9963 FITEM,5,9972 FITEM,5,9981 FITEM,5,10287 FITEM,5,10296 FITEM,5,10305FITEM,5,10323 FITEM,5,10332 FITEM,5,10341 FITEM,5,10350 FITEM,5,10359 FITEM,5,10368 FITEM,5,10377 FITEM,5,10386 FITEM,5,10395 FITEM,5,10404 FITEM,5,10413 FITEM,5,10422 FITEM,5,10431 FITEM,5,10440 FITEM,5,10449 FITEM,5,10458 FITEM,5,10467 FITEM,5,10476 FITEM,5,10485 FITEM,5,10494 FITEM,5,10503 FITEM,5,10512 FITEM,5,10521 FITEM,5,10530 FITEM,5,10539 FITEM,5,10548 NSEL,U, , ,P51X NPLOTFLST,5,20,1,ORDE,20FITEM,5,4815FITEM,5,5085FITEM,5,5391FITEM,5,5661FITEM,5,5967FITEM,5,6237FITEM,5,6543FITEM,5,6813FITEM,5,7389FITEM,5,7965FITEM,5,8271FITEM,5,8541FITEM,5,8847FITEM,5,9117FITEM,5,9423FITEM,5,9693FITEM,5,9999FITEM,5,10269FITEM,5,10845NSEL,U, , ,P51XNSEL,A,NODE,,10573,10574,1NSEL,A,NODE,,7693,7694,1NSEL,A,NODE,,7117,7118,1NSEL,A,NODE,,2665,2666,1CM,NODE1,NODENGEN,2,100000,ALL !创建第一层水管单元ngen复制节点TYPE,2MAT,2REAL,1 !单元实常数属性指示器E,110573,110574,10573,10574 !通过节点相连生成一个单元E,110574,110575,10574,10575*do,i,1,29,1e,110575+(i-1)*9,110584+(i-1)*9,10575+(i-1)*9,10584+(i-1)*9*enddoe,110836,110557,10836,10557e,110557,110260,10557,10260*do,i,1,28,1e,110260-(i-1)*9,110251-(i-1)*9,10260-(i-1)*9,10251-(i-1)*9 *enddoe,110008,109711,10008,9711e,109711,109432,9711,9432*do,i,1,28,1e,109432+(i-1)*9,109441+(i-1)*9,9432+(i-1)*9,9441+(i-1)*9 *enddoe,109684,109405,9684,9405e,109405,109108,9405,9108*do,i,1,28,1e,109108-(i-1)*9,109099-(i-1)*9,9108-(i-1)*9,9099-(i-1)*9*enddoe,108856,108559,8856,8559e,108559,108280,8559,8280*do,i,1,28,1e,108280+(i-1)*9,108289+(i-1)*9,8280+(i-1)*9,8289+(i-1)*9 *enddoe,108532,108253,8532,8253e,108253,107956,8253,7956*do,i,1,29,1e,107956-(i-1)*9,107947-(i-1)*9,7956-(i-1)*9,7947-(i-1)*9 *enddoE,107695,107694,7695,7694E,107694,107693,7694,7693E,102665,102666,2665,2666E,102666,102667,2666,2667*do,i,1,29,1e,102667+(i-1)*9,102676+(i-1)*9,2667+(i-1)*9,2676+(i-1)*9 *enddoe,102928,104797,2928,4797e,104797,105076,4797,5076*do,i,1,28,1e,105076-(i-1)*9,105067-(i-1)*9,5076-(i-1)*9,5067-(i-1)*9 *enddoe,104824,105103,4824,5103e,105103,105400,5103,5400*do,i,1,28,1e,105400+(i-1)*9,105409+(i-1)*9,5400+(i-1)*9,5409+(i-1)*9 *enddoe,105652,105949,5652,5949e,105949,106228,5949,6228*do,i,1,28,1e,106228-(i-1)*9,106219-(i-1)*9,6228-(i-1)*9,6219-(i-1)*9 *enddoe,105976,106255,5976,6255e,106255,106552,6255,6552*do,i,1,28,1e,106552+(i-1)*9,106561+(i-1)*9,6552+(i-1)*9,6561+(i-1)*9 *enddoe,106804,107101,6804,7101e,107101,107380,7101,7380*do,i,1,29,1e,107380-(i-1)*9,107371-(i-1)*9,7380-(i-1)*9,7371-(i-1)*9 *enddoE,107119,107118,7119,7118E,107118,107117,7118,7117 ! 第一层水管单元创建完毕! 选择第二层水管节点allsel,allNSEL,S,LOC,Y,3NSEL,U,LOC,X,0NSEL,U,LOC,z,0NSEL,U,LOC,Z,7.9,8.25NSEL,R,LOC,X,0.2,12.20FLST,5,351,1,ORDE,351FITEM,5,2933FITEM,5,4523FITEM,5,4532FITEM,5,4541FITEM,5,4550FITEM,5,4559FITEM,5,4568FITEM,5,4577FITEM,5,4586FITEM,5,4595FITEM,5,4604FITEM,5,4613FITEM,5,4622FITEM,5,4640 FITEM,5,4649 FITEM,5,4658 FITEM,5,4667 FITEM,5,4676 FITEM,5,4685 FITEM,5,4694 FITEM,5,4703 FITEM,5,4712 FITEM,5,4721 FITEM,5,4730 FITEM,5,4739 FITEM,5,4748 FITEM,5,4757 FITEM,5,4766 FITEM,5,4775 FITEM,5,4784 FITEM,5,4811 FITEM,5,5081 FITEM,5,5108 FITEM,5,5117 FITEM,5,5126 FITEM,5,5135 FITEM,5,5144 FITEM,5,5153 FITEM,5,5162 FITEM,5,5171 FITEM,5,5180 FITEM,5,5189FITEM,5,5207 FITEM,5,5216 FITEM,5,5225 FITEM,5,5234 FITEM,5,5243 FITEM,5,5252 FITEM,5,5261 FITEM,5,5270 FITEM,5,5279 FITEM,5,5288 FITEM,5,5297 FITEM,5,5306 FITEM,5,5315 FITEM,5,5324 FITEM,5,5333 FITEM,5,5342 FITEM,5,5351 FITEM,5,5360 FITEM,5,5369 FITEM,5,5387 FITEM,5,5657 FITEM,5,5675 FITEM,5,5684 FITEM,5,5693 FITEM,5,5702 FITEM,5,5711 FITEM,5,5720 FITEM,5,5729 FITEM,5,5738FITEM,5,5756 FITEM,5,5765 FITEM,5,5774 FITEM,5,5783 FITEM,5,5792 FITEM,5,5801 FITEM,5,5810 FITEM,5,5819 FITEM,5,5828 FITEM,5,5837 FITEM,5,5846 FITEM,5,5855 FITEM,5,5864 FITEM,5,5873 FITEM,5,5882 FITEM,5,5891 FITEM,5,5900 FITEM,5,5909 FITEM,5,5918 FITEM,5,5927 FITEM,5,5936 FITEM,5,5963 FITEM,5,6233 FITEM,5,6260 FITEM,5,6269 FITEM,5,6278 FITEM,5,6287 FITEM,5,6296FITEM,5,6314 FITEM,5,6323 FITEM,5,6332 FITEM,5,6341 FITEM,5,6350 FITEM,5,6359 FITEM,5,6368 FITEM,5,6377 FITEM,5,6386 FITEM,5,6395 FITEM,5,6404 FITEM,5,6413 FITEM,5,6422 FITEM,5,6431 FITEM,5,6440 FITEM,5,6449 FITEM,5,6458 FITEM,5,6467 FITEM,5,6476 FITEM,5,6485 FITEM,5,6494 FITEM,5,6503 FITEM,5,6512 FITEM,5,6521 FITEM,5,6539 FITEM,5,6809 FITEM,5,6827 FITEM,5,6836 FITEM,5,6845FITEM,5,6863 FITEM,5,6872 FITEM,5,6881 FITEM,5,6890 FITEM,5,6899 FITEM,5,6908 FITEM,5,6917 FITEM,5,6926 FITEM,5,6935 FITEM,5,6944 FITEM,5,6953 FITEM,5,6962 FITEM,5,6971 FITEM,5,6980 FITEM,5,6989 FITEM,5,6998 FITEM,5,7007 FITEM,5,7016 FITEM,5,7025 FITEM,5,7034 FITEM,5,7043 FITEM,5,7052 FITEM,5,7061 FITEM,5,7070 FITEM,5,7079 FITEM,5,7088 FITEM,5,7385 FITEM,5,7403 FITEM,5,7412FITEM,5,7430 FITEM,5,7439 FITEM,5,7448 FITEM,5,7457 FITEM,5,7466 FITEM,5,7475 FITEM,5,7484 FITEM,5,7493 FITEM,5,7502 FITEM,5,7511 FITEM,5,7520 FITEM,5,7529 FITEM,5,7538 FITEM,5,7547 FITEM,5,7556 FITEM,5,7565 FITEM,5,7574 FITEM,5,7583 FITEM,5,7592 FITEM,5,7601 FITEM,5,7610 FITEM,5,7619 FITEM,5,7628 FITEM,5,7637 FITEM,5,7646 FITEM,5,7655 FITEM,5,7664 FITEM,5,7673 FITEM,5,7961FITEM,5,7988 FITEM,5,7997 FITEM,5,8006 FITEM,5,8015 FITEM,5,8024 FITEM,5,8033 FITEM,5,8042 FITEM,5,8051 FITEM,5,8060 FITEM,5,8069 FITEM,5,8078 FITEM,5,8087 FITEM,5,8096 FITEM,5,8105 FITEM,5,8114 FITEM,5,8123 FITEM,5,8132 FITEM,5,8141 FITEM,5,8150 FITEM,5,8159 FITEM,5,8168 FITEM,5,8177 FITEM,5,8186 FITEM,5,8195 FITEM,5,8204 FITEM,5,8213 FITEM,5,8222 FITEM,5,8231 FITEM,5,8240FITEM,5,8537 FITEM,5,8564 FITEM,5,8573 FITEM,5,8582 FITEM,5,8591 FITEM,5,8600 FITEM,5,8609 FITEM,5,8618 FITEM,5,8627 FITEM,5,8636 FITEM,5,8645 FITEM,5,8654 FITEM,5,8663 FITEM,5,8672 FITEM,5,8681 FITEM,5,8690 FITEM,5,8699 FITEM,5,8708 FITEM,5,8717 FITEM,5,8726 FITEM,5,8735 FITEM,5,8744 FITEM,5,8753 FITEM,5,8762 FITEM,5,8771 FITEM,5,8780 FITEM,5,8789 FITEM,5,8798 FITEM,5,8807FITEM,5,8825 FITEM,5,8843 FITEM,5,9113 FITEM,5,9131 FITEM,5,9140 FITEM,5,9149 FITEM,5,9158 FITEM,5,9167 FITEM,5,9176 FITEM,5,9185 FITEM,5,9194 FITEM,5,9203 FITEM,5,9212 FITEM,5,9221 FITEM,5,9230 FITEM,5,9239 FITEM,5,9248 FITEM,5,9257 FITEM,5,9266 FITEM,5,9275 FITEM,5,9284 FITEM,5,9293 FITEM,5,9302 FITEM,5,9311 FITEM,5,9320 FITEM,5,9329 FITEM,5,9338 FITEM,5,9347 FITEM,5,9356FITEM,5,9374 FITEM,5,9383 FITEM,5,9392 FITEM,5,9419 FITEM,5,9689 FITEM,5,9716 FITEM,5,9725 FITEM,5,9734 FITEM,5,9743 FITEM,5,9752 FITEM,5,9761 FITEM,5,9770 FITEM,5,9779 FITEM,5,9788 FITEM,5,9797 FITEM,5,9806 FITEM,5,9815 FITEM,5,9824 FITEM,5,9833 FITEM,5,9842 FITEM,5,9851 FITEM,5,9860 FITEM,5,9869 FITEM,5,9878 FITEM,5,9887 FITEM,5,9896 FITEM,5,9905 FITEM,5,9914FITEM,5,9932 FITEM,5,9941 FITEM,5,9950 FITEM,5,9959 FITEM,5,9968 FITEM,5,9977 FITEM,5,9995 FITEM,5,10265 FITEM,5,10283 FITEM,5,10292 FITEM,5,10301 FITEM,5,10310 FITEM,5,10319 FITEM,5,10328 FITEM,5,10337 FITEM,5,10346 FITEM,5,10355 FITEM,5,10364 FITEM,5,10373 FITEM,5,10382 FITEM,5,10391 FITEM,5,10400 FITEM,5,10409 FITEM,5,10418 FITEM,5,10427 FITEM,5,10436 FITEM,5,10445 FITEM,5,10454 FITEM,5,10463FITEM,5,10481 FITEM,5,10490 FITEM,5,10499 FITEM,5,10508 FITEM,5,10517 FITEM,5,10526 FITEM,5,10535 FITEM,5,10544 FITEM,5,10841 NSEL,U, , ,P51X NSEL,A,NODE,,1506 NSEL,A,NODE,,10568 NSEL,A,NODE,,10569 NSEL,A,NODE,,10570NSEL,A,NODE,,7690 NSEL,A,NODE,,7689 NSEL,A,NODE,,7688 NSEL,A,NODE,,1186 NSEL,A,NODE,,7114 NSEL,A,NODE,,7113 NSEL,A,NODE,,7112 NSEL,A,NODE,,1122 NSEL,A,NODE,,2662 NSEL,A,NODE,,2661 NSEL,A,NODE,,2660 NSEL,A,NODE,,1567 cm,node2,nodengen,2,200000,allcm,node2_,node!创建第二层水管单元TYPE,2MAT,2REAL,1E,201506,210568,1506,10568E,210568,210569,10568,10569E,210569,210570,210569,210570E,210570,210571,10570,10571*do,i,1,29,1$e,200000+10571+(i-1)*9,200000+10580+(i-1)*9,10571+(i-1)*9,10580+(i-1)*9$*enddo$e,210832,210553,10832,10553$e,210553,210256,10553,10256*do,i,1,28,1$e,200000+10256-(i-1)*9,200000+10247-(i-1)*9,10256-(i-1)*9,10247-(i-1)*9$*enddo$e,200000+10004,200000+9707,10004,9707$e,200000+9707,200000+9428,9707,9428*do,i,1,28,1$e,200000+9428+(i-1)*9,200000+9437+(i-1)*9,9428+(i-1)*9,9437+(i-1)*9$*enddo$e,200000+9680,200000+9401,9680,9401$e,200000+9401,200000+9104,9401,9104*do,i,1,28,1$e,200000+9104-(i-1)*9,200000+9095-(i-1)*9,9104-(i-1)*9,9095-(i-1)*9$*enddo$e,200000+8852,200000+8555,8852,8555$e,200000+8555,200000+8276,8555,8276*do,i,1,28,1$e,200000+8276+(i-1)*9,200000+8285+(i-1)*9,8276+(i-1)*9,8285+(i-1)*9$*enddo$e,200000+8528,200000+8249,8528,8249$e,200000+8249,200000+7952,8249,7952*do,i,1,29,1$e,200000+7952-(i-1)*9,200000+7943-(i-1)*9,7952-(i-1)*9,7943-(i-1)*9$*enddoE,207691,207690,7691,7690E,207690,207689,7690,7689E,207689,207688,7689,7688E,207688,201186,7688,1186E,201567,202660,1567,2660E,202660,202661,2660,2661E,202661,202662,2661,2662E,202662,202663,2662,2663*do,i,1,29,1$e,200000+2663+(i-1)*9,200000+2672+(i-1)*9,2663+(i-1)*9,2672+(i-1)*9$*enddo$e,200000+2924,200000+4793,2924,4793$e,200000+4793,200000+5072,4793,5072$*do,i,1,28,1$e,200000+5072-(i-1)*9,200000+5063-(i-1)*9,5072-(i-1)*9,5063-(i-1)*9$*enddo$e,200000+4820,200000+5099,4820,5099$e,200000+5099,200000+5396,5099,5396$*do,i,1,28,1$e,200000+5396+(i-1)*9,200000+5405+(i-1)*9,5396+(i-1)*9,5405+(i-1)*9$*enddo$e,200000+5648,200000+5945,5648,5945$e,200000+5945,200000+6224,5945,6224$*do,i,1,28,1$e,200000+6224-(i-1)*9,200000+6215-(i-1)*9,6224-(i-1)*9,6215-(i-1)*9$*enddo$e,200000+5972,200000+6251,5972,6251$e,200000+6251,200000+6548,6251,6548$*do,i,1,28,1$e,200000+6548+(i-1)*9,200000+6557+(i-1)*9,6548+(i-1)*9,6557+(i-1)*9$*enddo$e,200000+6800,200000+7097,6800,7097$e,200000+7097,200000+7376,7097,7376$*do,i,1,29,1$e,200000+7376-(i-1)*9,200000+7367-(i-1)*9,7376-(i-1)*9,7367-(i-1)*9$*enddoE,207115,207114,7115,7114E,207114,207113,7114,7113E,207113,207112,7113,7112E,207112,201122,7112,1122!第二层水管创建完毕!选择第三层水管节点ALLSEL,ALLNSEL,S,LOC,Y,4NSEL,U,LOC,Z,0NSEL,U,LOC,Z,7.9,8.25NSEL,R,LOC,X,0.2,12.2FLST,5,351,1,ORDE,351FITEM,5,2931FITEM,5,4521FITEM,5,4530FITEM,5,4539FITEM,5,4548FITEM,5,4557FITEM,5,4566FITEM,5,4575FITEM,5,4584FITEM,5,4593FITEM,5,4602FITEM,5,4611FITEM,5,4620FITEM,5,4629FITEM,5,4647 FITEM,5,4656 FITEM,5,4665 FITEM,5,4674 FITEM,5,4683 FITEM,5,4692 FITEM,5,4701 FITEM,5,4710 FITEM,5,4719 FITEM,5,4728 FITEM,5,4737 FITEM,5,4746 FITEM,5,4755 FITEM,5,4764 FITEM,5,4773 FITEM,5,4782 FITEM,5,4809 FITEM,5,5079 FITEM,5,5106 FITEM,5,5115 FITEM,5,5124 FITEM,5,5133 FITEM,5,5142 FITEM,5,5151 FITEM,5,5160 FITEM,5,5169 FITEM,5,5178 FITEM,5,5187 FITEM,5,5196FITEM,5,5214 FITEM,5,5223 FITEM,5,5232 FITEM,5,5241 FITEM,5,5250 FITEM,5,5259 FITEM,5,5268 FITEM,5,5277 FITEM,5,5286 FITEM,5,5295 FITEM,5,5304 FITEM,5,5313 FITEM,5,5322 FITEM,5,5331 FITEM,5,5340 FITEM,5,5349 FITEM,5,5358 FITEM,5,5367 FITEM,5,5385 FITEM,5,5655 FITEM,5,5673 FITEM,5,5682 FITEM,5,5691 FITEM,5,5700 FITEM,5,5709 FITEM,5,5718 FITEM,5,5727 FITEM,5,5736 FITEM,5,5745FITEM,5,5763 FITEM,5,5772 FITEM,5,5781 FITEM,5,5790 FITEM,5,5799 FITEM,5,5808 FITEM,5,5817 FITEM,5,5826 FITEM,5,5835 FITEM,5,5844 FITEM,5,5853 FITEM,5,5862 FITEM,5,5871 FITEM,5,5880 FITEM,5,5889 FITEM,5,5898 FITEM,5,5907 FITEM,5,5916 FITEM,5,5925 FITEM,5,5934 FITEM,5,5961 FITEM,5,6231 FITEM,5,6258 FITEM,5,6267 FITEM,5,6276 FITEM,5,6285 FITEM,5,6294 FITEM,5,6303 FITEM,5,6312FITEM,5,6330 FITEM,5,6339 FITEM,5,6348 FITEM,5,6357 FITEM,5,6366 FITEM,5,6375 FITEM,5,6384 FITEM,5,6393 FITEM,5,6402 FITEM,5,6411 FITEM,5,6420 FITEM,5,6429 FITEM,5,6438 FITEM,5,6447 FITEM,5,6456 FITEM,5,6465 FITEM,5,6474 FITEM,5,6483 FITEM,5,6492 FITEM,5,6501 FITEM,5,6510 FITEM,5,6519 FITEM,5,6537 FITEM,5,6807 FITEM,5,6825 FITEM,5,6834 FITEM,5,6843 FITEM,5,6852 FITEM,5,6861FITEM,5,6879 FITEM,5,6888 FITEM,5,6897 FITEM,5,6906 FITEM,5,6915 FITEM,5,6924 FITEM,5,6933 FITEM,5,6942 FITEM,5,6951 FITEM,5,6960 FITEM,5,6969 FITEM,5,6978 FITEM,5,6987 FITEM,5,6996 FITEM,5,7005 FITEM,5,7014 FITEM,5,7023 FITEM,5,7032 FITEM,5,7041 FITEM,5,7050 FITEM,5,7059 FITEM,5,7068 FITEM,5,7077 FITEM,5,7086 FITEM,5,7383 FITEM,5,7401 FITEM,5,7410 FITEM,5,7419FITEM,5,7437 FITEM,5,7446 FITEM,5,7455 FITEM,5,7464 FITEM,5,7473 FITEM,5,7482 FITEM,5,7491 FITEM,5,7500 FITEM,5,7509 FITEM,5,7518 FITEM,5,7527 FITEM,5,7536 FITEM,5,7545 FITEM,5,7554 FITEM,5,7563 FITEM,5,7572 FITEM,5,7581 FITEM,5,7590 FITEM,5,7599 FITEM,5,7608 FITEM,5,7617 FITEM,5,7626 FITEM,5,7635 FITEM,5,7644 FITEM,5,7653 FITEM,5,7662 FITEM,5,7671 FITEM,5,7959 FITEM,5,7977FITEM,5,7995 FITEM,5,8004 FITEM,5,8013 FITEM,5,8022 FITEM,5,8031 FITEM,5,8040 FITEM,5,8049 FITEM,5,8058 FITEM,5,8067 FITEM,5,8076 FITEM,5,8085 FITEM,5,8094 FITEM,5,8103 FITEM,5,8112 FITEM,5,8121 FITEM,5,8130 FITEM,5,8139 FITEM,5,8148 FITEM,5,8157 FITEM,5,8166 FITEM,5,8175 FITEM,5,8184 FITEM,5,8193 FITEM,5,8202 FITEM,5,8211 FITEM,5,8220 FITEM,5,8229 FITEM,5,8238 FITEM,5,8265。
混凝土梁的ANSYS命令流
！ 水＊＊＊＊＊ ＊ ＊＊＊＊＊定义支座垫板的几何特性及材料特性＊＊＊＊＊＊ ＊ ＊
Ｅ ， Ｏ Ｉ４ Ｔ ３ Ｓ ＬＤ ５ ， Ｍ， ３３Ｏ １ ＰＥ， ． ０ Ｘ ， Ｅ
ＳＶ ＡＥ
！＊＊＊＊＊ 建立混凝土梁单元＊＊＊＊ ＊＊＊＊＊＊ ＊＊＊ ＊
附录 ２
Ｂ Ｏ Ｋ ０ Ｌ ０ ２ ， ． ＬＣ ， ０ ． ０ ０ ２ ， ， ５ ， ０ ＢＯ Ｋ ０ Ｌ０ ２ ， ５０ ０ ２ ＬＣ ， ， ５ ０ ３ ， ． ， ． ． ， ０ ＢＯ ＫＯ Ｌ ２ ， ５－ ． ０ ＬＣ ， ０ ５０ ３ ，０ １， ， ． ． ５ ＢＯ Ｋ ０ Ｌ０ ２ ， ５ ０２ ， ５ ＬＣ ， ， ５０ ３ ， ０ ０ ３ ， ． ． ． ． Ｖ Ｅ ， ． Ｓ ＬＡ Ｉ Ｌ ＶＬ ＥＡ Ｉ Ｇ Ｕ， Ｌ
ＳＶ ＡＥ
左翼缘构造筋
Ｌ 部构造钢筋划分网格 ［Ｅ ，， Ｅ，６ １， ＳｌＳ ＬＮ， １， １ Ｉ ６ ＬＥ， ＩＥ， ７１， Ｓ ＬＡ ＬＮ ，１， １ ， ７ ＬＴ ， ， ＡＴ ２ ３ ２ ， ＬＳ Ｚ ， ， ５ Ｅ ＩＥ Ａ Ｌ ００ Ｌ ．
ＬＥ ＩＡ Ｌ Ｍ ＳＩ Ｌ ， ｌ ｌＡＬ 名Ｅ ， Ｌ
Ｋ ２ ， ， ５ ０ １， ， ３ ５ ０ ０ ， ５ ４ ． ． ． Ｌ Ｔ ， ＳＲ ２ ， １ ２ ２ Ｌ Ｔ ， ＳＲ ２ ， ３ ２ ４
！＊＊＊＊＊下面纵向钢筋划分网格＊ ＊＊ ＊ ＊＊＊＊＊＊ ＊ ＊＊
ＬＥ ， ， Ｅ ， ３ １， Ｓ Ｌ Ｓ ＬＮ ， １， ２ Ｉ ５
Ｒ１ ，，
！ 定义混凝土的常数
Ｍ， ， ｅ０弹性模量要根据其他公式求 ＰＥ， ３ ｌ！ Ｘ１２ ４ ．
Ｍ ， Ｙ １０ ２ Ｐ Ｐ Ｘ ， ． Ｒ ， ０ Ｍ ， Ｓ ｌ２Ｅ 密度 ＰＤＮ ， ５３！ Ｅ ，
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
钢筋混凝土非线性分析2015大作业上海交通大学陈明1、参数选择梁的截面宽度为200mm，上部配置2Φ8受压筋，混凝土的净保护层厚度为25 mm（从纵向钢筋外边缘算起），箍筋两端区采用8@100的双肢箍，中间区取8@200 双肢箍1）梁的截面高度选300mm；2）两加载间的距离选1000mm；3）混凝土选C30；4）纵向受拉钢筋配筋选218；2、描述选用的有限元模型及单元的特点采用ansys软件进行模拟计算，钢筋混凝土模型采用分离式模型，不考虑钢筋与混凝土之间的相对滑移。
混凝土采用solid65单元模拟，solid65用于模拟三维有钢筋或无钢筋的混凝土模型。
该单元能够计算拉裂和压碎。
在混凝土应用中，该单元的实体功能可以用于建立混凝土模型，同时，还可用加筋功能建立钢筋混凝土模型。
另外，该单元还可以应用于加强复合物和地质材料。
该单元由八个节点定义，每个节点有三个自由度:节点坐标系的x，y，z方向的平动。
至多可以定义三种不同规格的钢筋。
钢筋单元采用link180单元模拟，link180是一个适用于各类工程应用的三维杆单元。
根据具体情况，该单元可以被看作桁架单元、索单元、链杆单元或弹簧单元等等。
本单元是一个轴向拉伸一压缩单元，每个节点有三个自由度:节点坐标系的x，y，z方向的平动。
本单元是一种顶端铰接结构，不考虑单元弯曲。
本单元具有塑性、蠕变、旋转、大变形和大应变功能。
缺省时，当考虑大变形时任何分析中LINK180单元都包括应力刚化选项。
3、描述选用的混凝土与钢筋粘结滑移本构关系的具体形式、参数等。
钢筋的应力应变关系曲线考虑到极限塑性应变最大值为0.01,钢筋本构模型采用多线性模型kinh，初始弹性模量为Es=200000Mpa，强化系数为0.001。
混凝土的应力应变关系曲线混凝土选用各向同性的miso模型，当计入下降端时，程序报错，所以只取了前面的上升段，用5段折线模拟混凝土应力应变曲线。
不考虑混凝土与钢筋之间的相对滑移4、迭代方法和收敛标准。
使用修正的Newton-Raphson迭代方法进行求解。
收敛标准采用位移来控制（cnvtol,u,,0.05,0,），将toler选项改成0.05来加速收敛。
计算时，程序将连续进行平衡迭代直到满足收敛准则，或者直到允许平衡迭代最大次数50次时停止迭代（neqit,50）。
5、画出解题的流程图选定相关的参数选择合适的单元来模拟混凝土、钢筋、垫板定义混凝土单元和钢筋单元及垫板的材料参数建立几何模型进行网格尺寸指定并赋予材料属性转换几何模型为有限元模型施加荷载和边界条件选择求解方法和屈服准则进行求解进行后处理得到图形和数据6、给出梁的应力云图、梁逐级加载的开裂图梁的应力云图Sub=15Sub=16Sub=17 Sub=18钢筋应力云图Sub=15Sub=19梁逐级加载的开裂图（第16子步开始开裂）Sub=16Sub=17Sub=18Sub=197、梁的荷载－跨中挠度曲线8、ansys命令流finish/clear/filname,beam参数输入/prep7pi=3.1415926l=3100 !梁纵向长度lp=1000 !荷载间距w=200 !梁宽度h=300 !梁深度p=25 !保护层厚度r_ps=4 !压筋半径a_ps=0.25*pi*r_ps**2r_ts=9 !拉筋半径a_ts=0.25*pi*r_ts**2r_gs=4 !箍筋半径a_gs=0.25*pi*r_gs**2fy=350e6 !钢筋屈服强度/mpa q=1e6 !均布载荷/mpae0=2e9 !钢筋强化系数ex_c=30e9 !混凝土杨氏模量prxy_c=0.2 !混凝土泊松比ex_ts=200e9 !拉筋杨氏模量prxy_ts=0.27 !拉筋泊松比ex_gs=210e9 !箍筋杨氏模量prxy_sp=0.27 !箍筋泊松比ex_sp=210e15 !垫板杨氏模量prxy_sp=0.48 !垫板泊松比选择单元类型et,1,solid65,,,,,,,1keyopt,1,7,1et,2,link180定义混凝土材料特性mp,ex,1,ex_cmp,prxy,1,prxy_ctb,miso,1,,8tbpt,,0.0001,3e6tbpt,,0.0005,10.1907e6tbpt,,0.001,15.8935e6tbpt,,0.0015,19.1027e6tbpt,,0.002,20.1e6!tbpt,,0.0025,19.3067e6!tbpt,,0.003,17.4327e6!tbpt,,0.0033,16.08e6tb,concr,1tbdata,,0.5,0.95,1.43e6,-1定义钢筋材料特性mp,ex,2,ex_tsmp,prxy,2,prxy_tstb,kinh,2,1,3tbpt,,fy/ex_ts,fytbpt,,fy/ex_ts+0.1,350e6+e0*0.1tbpt,,fy/ex_ts+0.1001,1000mp,ex,3,ex_gsmp,prxy,3,prxy_gstb,bkin,3,1,2,1tbdata,,fy,E0定义垫板材料特性mp,ex,4,ex_spmp,prxy,4,prxy_sb定义实常数r,1,r,2,a_tsr,3,a_gsr,4,a_ps建立几何模型blc4,,,l,h,wblock,100,200,0,-h/10,0,wblock,l-200,l-100,0,-h/10,0,wvglue,allwpoffs,,,(p+r_ts/2)vsbw,allwpcsys,-1wpoffs,,,w-(p+r_ts/2)vsbw,allwpcsys,-1切出箍筋位置wprota,,,90*do,i,1,11wpoffs,,,100vsbw,all*enddowpoffs,,,50*do,i,1,3wpoffs,,,200vsbw,all*enddowpoffs,,,250*do,i,1,11vsbw,allwpoffs,,,100*enddowpcsys,-1切出支座位置wprota,,,90wpoffs,,,150vsbw,allwpcsys,-1wprota,,,90wpoffs,,,l-150vsbw,allwpcsys,-1切出纵向钢筋位置wprota,,90wpoffs,,,-(p+r_ts/2)vsbw,allwpoffs,,,-h+2*(p+r_ts/2)vsbw,allwpcsys,-1vglue,all指定划分精度elemsiz=w/4划分纵向钢筋并赋予属性lsel,s,loc,y,p+r_ts/2lsel,r,loc,z,p+r_ts/2cm,ts1,linelatt,2,2,2lesize,all,elemsizlsel,s,loc,y,p+r_ts/2lsel,r,loc,z,w-(p+r_ts/2)cm,ts2,linelatt,2,2,2lesize,all,elemsizlsel,s,loc,y,h-(p+r_ts/2)lsel,r,loc,z,p+r_ts/2cm,ps1,linelatt,2,4,2lesize,all,elemsizlsel,s,loc,y,h-(p+r_ts/2)lsel,r,loc,z,w-(p+r_ts/2)cm,ps2,linelatt,2,4,2lesize,all,elemsiz划分箍筋并赋予属性lsel,s,loc,z,(p+r_ts/2)lsel,r,loc,y,h/2lsel,u,loc,x,150lsel,u,loc,x,l-150lsel,u,loc,x,1050lsel,u,loc,x,l-1050cm,gs1,linelatt,3,3,2lesize,all,elemsizlsel,s,loc,z,w-(p+r_ts/2)lsel,r,loc,y,h/2lsel,u,loc,x,150lsel,u,loc,x,l-150lsel,u,loc,x,1050lsel,u,loc,x,l-1050cm,gs2,linelatt,3,3,2lesize,all,elemsizlsel,s,loc,y,p+r_ts/2lsel,r,loc,z,w/2lsel,u,loc,x,150lsel,u,loc,x,l-150lsel,u,loc,x,1050lsel,u,loc,x,l-1050cm,gs3,linelatt,3,3,2lesize,all,elemsizlsel,s,loc,y,h-(p+r_ts/2)lsel,r,loc,z,w/2lsel,u,loc,x,150lsel,u,loc,x,l-150lsel,u,loc,x,1050lsel,u,loc,x,l-1050cm,gs4,linelatt,3,3,2lesize,all,elemsizcmsel,s,ts1cmsel,a,ts2cmsel,a,ps1cmsel,a,ps2cmsel,a,gs1cmsel,a,gs2cmsel,a,gs3cmsel,a,gs4lmesh,alllsel,all划分混凝土及垫板并赋予属性vsel,s,loc,y,-h/20vatt,4,,1mshkey,1esize,elesizvsel,allvsel,u,loc,y,-h/20vatt,1,,1mshkey,1esize,elesizvsel,allvmesh,all施加荷载与约束asel,s,loc,y,hasel,r,loc,x,1000,1100sfa,all,1,pres,qasel,s,loc,y,hasel,r,loc,x,2100,2200sfa,all,1,pres,qasel,allsftranlsel,s,loc,x,l-150lsel,r,loc,y,-h/10dl,all,,uylsel,s,loc,x,150lsel,r,loc,y,-h/10dl,all,,uydl,all,,uxdl,all,,uzlsel,alldtranasel,alllsel,allasel,allvsel,all求解/soluantype,0 !静力分析nlgeom,1 !打开几何大变形nropt,modi !用修正法进行求解eqslv,spar,,0nsubst,200,200,50,1kbc,0 !使用线性荷载cnvtol,u,,0.05,0, !指定收敛值neqit,50 !最大平衡迭代次数outres,all,allsolvefinish。