钢筋混凝土梁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命令流

ＳＶＡＥ
！水＊＊＊＊＊＊＊＊＊＊＊定义支座垫板的几何特性及材料特性＊＊＊＊＊＊＊＊
Ｅ，ＯＩ４Ｔ３ＳＬＤ５，Ｍ，３３Ｏ１ＰＥ，．０Ｘ，Ｅ
ＳＶＡＥ
！＊＊＊＊＊建立混凝土梁单元＊＊＊＊＊＊＊＊＊＊＊＊＊＊
附录２
ＢＯＫ０Ｌ０２，．ＬＣ，０．００２，，５，０ＢＯＫ０Ｌ０２，５００２ＬＣ，，５０３，．，．．，０ＢＯＫＯＬ２，５－．０ＬＣ，０５０３，０１，，．．５ＢＯＫ０Ｌ０２，５０２，５ＬＣ，，５０３，００３，．．．．ＶＥ，．ＳＬＡＩＬＶＬＥＡＩＧＵ，Ｌ
ＳＶＡＥ
左翼缘构造筋
Ｌ部构造钢筋划分网格［Ｅ，，Ｅ，６１，ＳｌＳＬＮ，１，１Ｉ６ＬＥ，ＩＥ，７１，ＳＬＡＬＮ，１，１，７ＬＴ，，ＡＴ２３２，ＬＳＺ，，５ＥＩＥＡＬ００Ｌ．
ＬＥＩＡＬＭＳＩＬ，ｌｌＡＬ名Ｅ，Ｌ
Ｋ２，，５０１，，３５００，５４．．．ＬＴ，ＳＲ２，１２２ＬＴ，ＳＲ２，３２４
！＊＊＊＊＊下面纵向钢筋划分网格＊＊＊＊＊＊＊＊＊＊＊＊＊
ＬＥ，，Ｅ，３１，ＳＬＳＬＮ，１，２Ｉ５
Ｒ１，，
！定义混凝土的常数
Ｍ，，ｅ０弹性模量要根据其他公式求ＰＥ，３ｌ！Ｘ１２４．
Ｍ，Ｙ１０２ＰＰＸ，．Ｒ，０Ｍ，Ｓｌ２Ｅ密度ＰＤＮ，５３！Ｅ，

钢筋混凝土结构Ansys命令流

钢筋混凝土结构Ansys命令流整体式finish/clear/config,nres,2000/prep7!定义单元与材料性质et,1,solid65,,,,,,,1mp,ex,1,13585mp,prxy,1,0.2mp,ex,2,2.0e5mp,prxy,2,0.3fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,fctb,miso,1,,11tbpt,,0.0002,fc*0.19tbpt,,0.0004,fc*0.36tbpt,,0.0006,fc*0.51tbpt,,0.0008,fc*0.64tbpt,,0.0010,fc*0.75tbpt,,0.0012,fc*0.84tbpt,,0.0014,fc*0.91tbpt,,0.0016,fc*0.96tbpt,,0.0018,fc*0.99tbpt,,0.002,fctbpt,,0.0033,fc*0.85tb,biso,2tbdata,,300,0pi=acos(-1)v1=0.25*pi*22*22*1.5/(75*60) v2=0.25*pi*10*10/(75*60)v3=0.25*pi*10*10*9.5/(1000*75) v4=0.25*pi*10*10*9.5/(1000*60) r,1,2,v1,0,90,2,v4,0,0,2,v3,90,0 r,2,2,v2,0,90,2,v4,0,0,2,v3,90,0 r,3,2,v3,90,0！几何建模blc4,,,150/2,300,2000/2 wprota,,-90wpoff,,,60vsbw,allwpoff,,,180vsbw,allwpcsys,-1wpoff,,,50vsbw,allwpoff,,,850vsbw,allwpcsys,-1！定义属性划分单元网格vsel,s,loc,y,0,60vatt,1,1,1vsel,s,loc,y,240,300vatt,1,2,1vsel,s,loc,y,60,240vatt,1,3,1esize,50mshkey,1vmesh,all!device,vector,on/eshape,1 ！此两个命令查看钢筋设置！施加荷载和约束lsel,s,loc,y,0lsel,r,loc,z,900dl,all,,uyasel,s,loc,z,0da,all,symmasel,s,loc,x,75da,all,symmp0=180000q0=p0/150/100asel,s,loc,z,0,50asel,r,loc,y,300sfa,all,1,pres,q0allsel,all!求解控制/soluantype,0nsubst,50outres,all,allautots,onneqit,50cnvtol,u,,0.015solve！进入post1查看结果/post1pldisp,1etable,reb1,smisc,2pletab,reb1/device,vector,onplcrack!时程后处理/post26nsol,2,69,u,yprod,3,2,,,,,,-1prod,4,1,,,,,,p0/1000xvar,3plvar,4分离式! 关闭压碎,keyopt(1)=0,keyopt(7)=1 finish/clear/config,nres,2000/prep7! 定义单元与材料性质et,1,solid65,,,,,,,1et,2,link8mp,ex,1,13585mp,prxy,1,0.2fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,fc tb,miso,1,,11 tbpt,,0.0002,fc*0.19 tbpt,,0.0004,fc*0.36 tbpt,,0.0006,fc*0.51 tbpt,,0.0008,fc*0.64 tbpt,,0.0010,fc*0.75 tbpt,,0.0012,fc*0.84 tbpt,,0.0014,fc*0.91 tbpt,,0.0016,fc*0.96 tbpt,,0.0018,fc*0.99 tbpt,,0.002,fc tbpt,,0.0033,fc*0.85 mp,ex,2,2.0e5mp,prxy,2,0.3tb,biso,2tbdata,,300,0pi=acos(-1)r,1,0.25*pi*22*22r,2,0.25*pi*22*22/2r,3,0.25*pi*10*10r,4,0.25*pi*10*10/2!创建几何模型blc4,,,150/2,300,2000/2 *do,i,1,9 wpoff,,,100vsbw,all*enddowpcsys,-1 wpoff,,,50vsbw,all wpcsys,-1 wprota,,-90 wpoff,,,30 vsbw,all wpoff,,,240 vsbw,all wpcsys,-1 wpoff,30 wprota,,,90 vsbw,all wpcsys,-1!划分钢筋网格elemsiz=50 lsel,s,loc,x,30 lsel,r,loc,y,30 cm,zj,line latt,2,1,2 lesize,all,elemsiz lsel,s,loc,x,75 lsel,r,loc,y,30 cm,zjb,line latt,2,2,2 lesize,all,elemsiz lsel,s,loc,x,30 lsel,r,loc,y,270 cm,jlj,line latt,2,3,2 lesize,all,elemsiz lsel,s,tan1,zlsel,r,loc,y,30,270 lsel,r,loc,x,30,70 lsel,u,loc,z,50 cm,gj,line latt,2,3,2 lesize,all,elemsiz lsel,s,loc,z,0 lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 cm,gjb,line latt,2,4,2 lesize,all,elemsiz lsel,all cmsel,s,zj cmsel,a,zjb cmsel,a,jlj cmsel,a,gj cmsel,a,gjb cm,gj,line lmesh,alllsel,all/eshape,1 eplot!划分混凝土网格vatt,1,,1 mshkey,1 esize,elemsiz vmesh,all allsel,all!施加荷载和约束lsel,s,loc,y,0lsel,r,loc,z,900dl,all,,uyasel,s,loc,z,0da,all,symmasel,s,loc,x,75da,all,symmp0=180000q0=p0/150/100 asel,s,loc,z,0,50 asel,r,loc,y,300 sfa,all,1,pres,q0 allsel,all!求解控制/soluantype,0nsubst,100 outres,all,all autos,onneqit,50cnvtol,u,,0.015 solve！进入post1查看结果/post1set,lastpldisp,1esel,s,type,,2 etable,saxl,ls,1 plls,saxl,saxlesel,s,type,,1/device,vector,on plcrack!时程后处理/post26nsol,2,205,u,y prod,3,2,,,,,,-1 prod,4,1,,,,,,p0/1000 xvar,3plvar,4预应力混凝土结构!实体切分法finish/clear/prep7!定义变量emst=1.95e5as=139tf=180e3denss=7921e-12 emcon=3.45e4 densc=2300e-12r0=9000b=150h=200d0=40dd=200-2*d0 spanl=3000d1=(39-3*sqrt(29))/35*dd-d0 !定义单元与材料性质et,1,solid95et,2,link8mp,ex,1,emconmp,prxy,1,0.2mp,dens,1,denscmp,ex,2,emstmp,prxy,2,0.3mp,dens,2,denssr,1r,2,as,tf/(emst*as)*1.0271!切分体形成力筋线blc4,,,b,h,spanllsel,nonek,9,-10,h-d0k,10,-10,-d1,spanl/2k,11,-10,h-d0,spanll,9,10l,10,11lfillt,13,14,r0lcomb,allk,14,b+10,h-d0l,9,14adrag,13,,,,,,14lsel,allvsba,1,7wpoff,b/2wprota,,,90vsbw,allwpcsys,-1 wpoff,,,spanl/2 vsbw,all!划分单元网格lsel,s,radius,,r0 lsel,r,loc,x,b/2latt,2,2,2lesize,all,40 lmesh,allesize,40vatt,1,1,1 mshape,0 mshkey,1 vmesh,allallsel,all!定义约束和荷载并求解lsel,s,loc,z,0lsel,a,loc,z,spanl lsel,r,loc,y,0dl,all,,uydk,kp(0,0,0),ux,,,,uz dk,kp(0,0,spanl),ux allsel,all/soluacel,,9800solvefinish!查看结果/post1pldisp,1etable,sigs,ls,1 plls,sigs,sigs,1 path,kzhx,2 ppath,1,,0,0,spanl/2 ppath,2,,b,0,spanl/2 pdef,sigc,s,z plpath,sigcpath,kzsx,2ppath,1,,b/2,0,spanl/2 ppath,2,,b/2,h,spanl/2 pdef,sigc,s,z plpath,sigc!节点耦合法finish/clear/prep7!定义变量emst=1.95e5as=139tf=180e3denss=7921e-12 emcon=3.45e4 densc=2300e-12r0=9000b=150h=200d0=40dd=200-2*d0spanl=3000d1=(39-3*sqrt(29))/35*dd-d0 !定义单元与材料性质et,1,solid95et,2,link8mp,ex,1,emconmp,prxy,1,0.2mp,dens,1,denscmp,ex,2,emstmp,prxy,2,0.3mp,dens,2,denssr,1r,2,as,tf/(emst*as)*1.0271!定义单元与材料性质et,1,solid95et,2,link8mp,ex,1,emconmp,prxy,1,0.2mp,dens,1,denscmp,ex,2,emstmp,prxy,2,0.3mp,dens,2,denssr,1r,2,as,tf/(emst*as)*1.01225!分别创建体和力筋线blc4,,,b,h,spanllsel,nonek,9,b/2,h-d0k,10,b/2,-d1,spanl/2k,11,b/2,h-d0,spanll,9,10l,10,11lfillt,13,14,r0lcomb,all wpoff,,,spanl/2lsbw,allcm,sline,linelsel,allvsbw,allwpcsys,-1!分别划分单元网格cmsel,s,slinelatt,2,2,2lesize,all,40lmesh,allesize,40vatt,1,1,1lsel,s,length,,h lesize,all,,,10lsel,allmshape,0mshkey,1vmesh,all!定义节点耦合自由度cmsel,s,slinensll,s,1cm,gsnode,node*get,senum,node,,count *dim,snode,,senum*dim,cnode,,senum*get,nodei,node,,num,min snode(1)=nodei*do,i,2,senumsnode(i)=ndnext(snode(i-1)) *enddoallsel,allcmsel,u,gsnode*do,i,1,senumnodei=snode(i)nodej=nnear(nodei) cnode(i)=nodejnsel,u,,,nodej*enddoallsel,all*do,i,1,senumcp,next,all,snode(i),cnode(i) *enddo!定义约束和荷载并求解lsel,s,loc,z,0lsel,a,loc,z,spanllsel,r,loc,y,0dl,all,,uydk,kp(0,0,0),ux,,,,uzdk,kp(0,0,spanl),ux allsel,all/soluacel,,9800solvefinish!查看结果/post1pldisp,1etable,sigs,ls,1 plls,sigs,sigs,1 path,kzhx,2 ppath,1,,0,0,spanl/2 ppath,2,,b,0,spanl/2 pdef,sigc,s,z plpath,sigcpath,kzsx,2ppath,1,,b/2,0,spanl/2 ppath,2,,b/2,h,spanl/2 pdef,sigc,s,z plpath,sigc!约束方程法finish/clear/prep7!定义变量emst=1.95e5as=139tf=180e3denss=7921e-12 emcon=3.45e4 densc=2300e-12r0=9000b=150h=200d0=40dd=200-2*d0spanl=3000d1=(39-3*sqrt(29))/35*dd-d0 !定义单元与材料性质et,1,solid95et,2,link8mp,ex,1,emconmp,prxy,1,0.2mp,dens,1,denscmp,ex,2,emstmp,prxy,2,0.3mp,dens,2,denssr,1r,2,as,tf/(emst*as)*1.0271!定义单元与材料性质et,1,solid95et,2,link8mp,ex,1,emconmp,prxy,1,0.2mp,dens,1,denscmp,ex,2,emstmp,prxy,2,0.3mp,dens,2,denssr,1r,2,as,tf/(emst*as)*1.01225!分别创建体和力筋线blc4,,,b,h,spanllsel,nonek,9,b/2,h-d0k,10,b/2,-d1,spanl/2k,11,b/2,h-d0,spanl l,9,10l,10,11lfillt,13,14,r0 lcomb,all wpoff,,,spanl/2 lsbw,allcm,sline,linelsel,allvsbw,allwpcsys,-1!分别划分单元网格cmsel,s,slinelatt,2,2,2lesize,all,40 lmesh,allesize,40vatt,1,1,1lsel,s,length,,h lesize,all,,,10lsel,allmshape,0 mshkey,1 vmesh,all!定义约束方程cmsel,s,slinensll,s,1ceintf,,ux,uy,uz allsel,all!定义约束和荷载并求解lsel,s,loc,z,0lsel,a,loc,z,spanl lsel,r,loc,y,0dl,all,,uydk,kp(0,0,0),ux,,,,uz dk,kp(0,0,spanl),ux allsel,all/soluacel,,9800solvefinish!查看结果/post1pldisp,1etable,sigs,ls,1 plls,sigs,sigs,1 path,kzhx,2 ppath,1,,0,0,spanl/2 ppath,2,,b,0,spanl/2 pdef,sigc,s,z plpath,sigcpath,kzsx,2ppath,1,,b/2,0,spanl/2 ppath,2,,b/2,h,spanl/2 pdef,sigc,s,z plpath,sigc。

ANSYS工程结构数值分析---命令流实例

!------------------------------------------------------!EX8.26 钢筋混凝土简支梁数值分析!分离式模型,关闭压碎,keyopt(1)=0,keyopt(7)=1 !力加载,位移收敛准则，误差1.5%,1/4模型分析!--------------------------------------------finish/clear/config,nres,2000/prep7!1.定义单元与材料性质--------------------et,1,solid65,,,,,,,1 !K1=0,k7=1et,2,link8mp,ex,1,13585mp,prxy,1,0.2fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,-1tb,miso,1,,11tbpt,,0.0002,fc*0.19tbpt,,0.0004,fc*0.36tbpt,,0.0006,fc*0.51tbpt,,0.0008,fc*0.64tbpt,,0.001,fc*0.75tbpt,,0.0012,fc*0.84tbpt,,0.0014,fc*0.91tbpt,,0.0016,fc*0.96tbpt,,0.0018,fc*0.99tbpt,,0.002,fctbpt,,0.0033,fc*0.85mp,ex,2,2.0e5mp,prxy,2,0.3tb,biso,2tbdata,,300,0pi=acos(-1)r,1,0.25*pi*22*22r,2,0.25*pi*22*22/2r,3,0.25*pi*10*10r,4,0.25*pi*10*10/2!2.创建几何模型blc4,,,150/2,300,2000/2*do,i,1,9wpoff,,,100vsbw,all*enddo wpcsys,-1 wpoff,,,50 vsbw,all wpcsys,-1 wprota,,-90 wpoff,,,30 vsbw,all wpoff,,,240 vsbw,all wpcsys,-1 wpoff,30 wprota,,,90 vsbw,all wpcsys,-1!3.划分钢筋网格elemsiz=50 lsel,s,loc,x,30 lsel,r,loc,y,30 cm,zj,linelatt,2,1,2 lesize,all,elemsiz lsel,s,loc,x,75 lsel,r,loc,y,30 cm,zjb,line latt,2,2,2 lesize,all,elemsiz lsel,s,loc,x,30 lsel,r,loc,y,270 cm,jlj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,tan1,z lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 lsel,u,loc,z,50 cm,gj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,loc,z,0 lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 cm,gjb,line latt,2,4,2lesize,all,elemsizlsel,allcmsel,s,zjcmsel,a,zjbcmsel,a,jljcmsel,a,gjcmsel,a,gjbcm,gj,linelmesh,alllsel,all!4.划分混凝土网格vatt,1,,1mshkey,1esize,elemsizvmesh,allallsel,all!5.施加荷载和约束lsel,s,loc,y,0lsel,r,loc,z,900dl,all,,uyasel,s,loc,z,0da,all,symmasel,s,loc,x,75da,all,symm!-----------------p0=180000q0=p0/150/100asel,s,loc,z,0,50asel,r,loc,y,300sfa,all,1,pres,q0allsel,all!6.求解控制设置/soluantype,0nsubst,80outres,all,allautos,onneqit,50cnvtol,u,,0.015solve!7.进入POST1查看结果/post1set,lastesel,s,type,,2etable,saxl,ls,1plls,saxl,saxlesel,s,type,,1/device,vector,onplcrackeall!8.进入时程后处理/post26nsol,2,205,u,yprod,3,1,,,,,,p0/1000prod,4,2,,,,,,-1xvar,4plvar,3!EX2.6D 60kg/m钢轨实体建模FINISH/CLEAR/PREP7!1.定义参数H1=30.5H2=48.5H3=48.5B1=45.75B2=29.25B3=16.5B4=36.5B5=10XD1=1/3XD2=1/9XD3=1/3XD4=1/20R0=400R1=2R2=4R3=40R4=20R5=25R7=13R8=80R9=300R10=2FAI=43H=H1+2*H2+H3!2.创建轨底主要组成部分的线K,1K,2,B1+B2K,3,B1+B2,H1-B1*XD1-B2*XD2K,4,B1,H1-B1*XD1K,5,0,H1*DO,I,1,4L,I,I+1*ENDDO!3.计算确定轨头下三个圆弧的位置CTA=A TAN(XD3)BTA=(90-FAI)*ACOS(-1)/180-CTAA=1/TAN(CTA)C1=B3/2+R0-R6*SIN(CTA)-(R5-R6)*COS(BTA)C2=R5*SIN(BTA)+R6*COS(CTA)-R6*SIN(BTA)-H2C=(C1-C2/TAN(CTA))/(R0-R5)C3=A*C-SQRT(A*A-C*C+1)C3=C3/(A*A+1)REFA=ASIN(C3)C1=R6*SIN(CTA)+(R5-R6)*COS(BTA)+(R0-R5)*COS(REFA) !4.创建轨腰和轨头结构线ROUX=B3/2+R0-C1ROUY=ROUX*TAN(CTA)ANGI=ASIN((H1+H2)/R0)K,6,B3/2+R0-R0*COS(ANGI)K,7,B3/2+R0-R0*COS(REFA),H1+H2+R0*SIN(REFA) LARC,6,7,2,R0DX1=R6*(COS(BTA)-SIN(CTA))DY1=R6*(COS(CTA)-SIN(BTA))K,8,ROUX-DX1,H1+2*H2+ROUY-DY1K,9,ROUX,H1+2*H2+ROUYLARC,7,8,2,R5LARC,8,9,2,R6K,10,B4,H1+2*H2+B4*XD3K,11,B4-(H3-B4*XD3)*XD4,HL,9,10L,10,11ANGI=ASIN(B5/R9)K,12,0,HK,13,B5,H-(R9-R9*COS(ANGI)) LARC,12,13,1,R9DX1=B4-B5+R8*SIN(ANGI)ANGI1=ASIN(DX1/R8)Y1=R8*COS(ANGI)-R8*COS(ANGI1) K,14,B4,KY(13)-Y1LARC,13,14,1,R8!5.倒角或弧线连接LFILLT,1,2,R1LFILLT,2,3,R2LFILLT,3,4,R3LPTN,4,5LFILLT,17,18,R4LFILLT,8,9,R10LPTN,9,11LFILLT,20,22,R7LDELE,15,16,1,1LDELE,19,21,2,1LSYMM,X,ALLNUMMRG,ALLNUMCMP,ALLAL,ALLASUMVOFFST,1,300!EX2.9B 斜向圆台建模finish/clear/prep7r0=3r1=2h=4cyl4,,,r0wpoff,r1,,hcyl4,,,r1*do,i,1,4askin,i,i+4*enddova,all!斜向延伸创建圆台vgen,1,1,,,,3*r0,,,,1wpcsyscyl4,,,r0vext,7,,,r1,,h,r1/r0,r1/r0!EX2.9B 已知函数方程时的建模!玫瑰线1==================================== !极坐标方程为ρ=asin(4Φ/3)finish/clear/prep7a=20csys,1*do,i,0,360*3fei=i*acos(-1)/180r0=a*sin(fei*4/3)k,,r0,i*enddo*do,i,1,360*3l,i,i+1*enddo!玫瑰线2=================================== !极坐标方程为ρ=asin(4Φ)finish/clear/prep7a=20csys,1*do,i,0,360fei=i*acos(-1)/180r0=a*sin(fei*4)k,,r0,i*enddo*do,i,1,360l,i,i+1*enddo!长辐圆内旋轮线============================ !参数方程x=(a-b)cost+λcos(a-b)/bt!参数方程y=(a-b)sint-λsin(a-b)/btfinish/clear/prep7a=50b=5lmda=15*afun,deg*do,i,0,360x1=(a-b)*cos(i)+lmda*cos((a-b)/b*i)y1=(a-b)*sin(i)-lmda*sin((a-b)/b*i)k,,x1,y1*enddo*do,i,1,360l,i,i+1*enddo!椭圆抛物面================================ !方程为z=x*x/a/a+y*y/b/b!使用蒙皮方法finish/clear/prep7a=20b=10n=20*do,i,1,niz=i/10lsel,nonewpoff,,,izcyl4,,,1adele,alllsscale,all,,,a*sqrt(iz),b*sqrt(iz),,,,1*enddoallsel*do,i,1,4*(n-1),4askin,i,i+4askin,i+1,i+5askin,i+2,i+6askin,i+3,i+7*enddo!EX2.9C 旋转图素建模!类花瓶建模--------------------finish/clear/prep7K,1K,2,100K,3,300,400K,4,150,550k,5,130,880K,6,300,1150k,7,400,1100k,8,500,1150bsplin,allk,1001,,1000arotat,all,,,,,,1,1001,360,5!类轮建模----------------------finish/clear/prep7!先创建类轮结构的一个断面，然后旋转之blc4,,,4,6cyl4,2,6,1.5asba,1,2k,100,-3k,101,-3,10vrotat,all,,,,,,100,101!EX3.1 自定义箱形截面finish/clear/prep7k,1k,2,2k,3,2,2.2k,4,3,2.3k,5,3,2.5k,6,0,2.5a,1,2,3,4,5,6blc4,,0.2,1.7,2asba,1,2wpoff,1.7wprota,,,90asbw,allwpoff,,,0.3asbw,allwpoff,,0.2wprota,,90asbw,allwpoff,,,-2asbw,allarsym,x,allwpcsysaglue,allet,1,plane82lesize,all,,,1amesh,allsecwrite,myboxfinish/clear/prep7et,1,beam189mp,ex,1,3.0e10mp,prxy,1,0.167sectype,1,beam,mesh secread,mybox,,,meshk,1k,2,10k,100,5,5l,1,2lesize,all,,,20latt,1,,1,,100,,1lmesh,all/eshape,1Eplot!EX3.2 自定义多种材料截面finish/clear/prep7Ro=1.5Ri=1.0csys,1cyl4,,,ricyl4,,,roaptn,allwprota,,90asbwa,allwprota,,,90asbw,allwpcsyset,1,plane82mymat1=4mymat2=7mp,ex,mymat1,1.0mp,ex,mymat2,2.0asel,s,loc,x,0,riaatt,mymat1,,1asel,s,loc,x,ri,roaatt,mymat2,,1allselesize,0.25mshape,0,2dmshkey,1amesh,allsecwrite,mycsolid,sectfinish/clear/prep7et,1,beam189mym1=4mym2=7mp,ex,mym1,3.0e10mp,prxy,mym1,0.167 mp,ex,mym2,2.1e11mp,prxy,mym2,0.3 sectype,1,beam,mesh secread,mycsolid,sect,,mesh k,1k,2,,,10l,1,2lesize,all,,,20latt,,,1,,,,1lmesh,all/eshape,1/pnum,mat,1eplot!EX3.3A 单个方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8K,1K,2,,,1000L,1,2K,100,,500,500LATT,1,,1,,100,,1LGEN,4,1,,,500LESIZE,ALL,,,10LMESH,ALL/ESHAPE,1EPLOT!EX3.3B 同时使用KB和KE方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0L,1,2K,100,,DLK,200,DXC,-DLK,300,2*DXC,DLK,301,2*DXC+DLK,400K,500,8*DXCLGEN,5,1,,,DXCLATT,1,,1,,100,,1LSEL,S,,,2LATT,1,,1,,200,,1LSEL,S,,,3LATT,1,,1,,300,301,1LSEL,S,,,4LATT,1,,1,,400,,1LSEL,S,,,5LATT,1,,1,,500,,1LSEL,ALLLESIZE,ALL,,,50LMESH,ALL/ESHAPE,1EPLOT!EX3.3C 双方位关键点翻转示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,I SECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0K,3,DXCK,4,DXC,,L0/2K,5,DXC,,L0L,1,2L,3,4L,4,5K,100,,DLk,101,,-DLK,200,DXC,DLk,201,2*DXCK,202,DXC,-DLLATT,1,,1,,100,101,1lesize,all,,,100LSEL,S,,,2LATT,1,,1,,200,201,1lesize,all,,,50LSEL,S,,,3LATT,1,,1,,201,202,1lesize,all,,,50LSEL,ALLLMESH,ALL/ESHAPE,1EPLOT!EX3.4A 两种单元形状和两种网格划分比较finish/clear/prep7et,1,plane82k,1k,2,8k,3,7,6k,4,1,6a,1,2,3,4esize,1mshape,0mshkey,0!mshape,0!mshkey,1!mshape,1!mshkey,0!mshape,1!mshkey,1amesh,all!EX3.4B 中间节点位置控制网格划分比较finish/clear/prep7cyl4,,,4,,8,60lesize,all,,,2mshape,0mshkey,1mshmid,0!mshmid,1!mshmid,2amesh,all!EX3.5 线上单元尺寸设置示例!下边密上边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10lsel,s,loc,x,0lesize,all,,,9,1/8lsel,s,loc,x,10lesize,all,,,9,8lsel,allmshape,0mshkey,1amesh,all!中间密外边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10,-1/5lsel,s,tan1,xlesize,all,,,9,-1/8lsel,allmshape,0mshkey,1amesh,all!EX3.6 合并线和连接线以进行映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,3mshape,0mshkey,1lccat,1,2lcomb,4,5amesh,all!EX3.7 合并线或连接线的网格划分设置finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,10lesize,6,,,8lesize,1,,,4lesize,2,,,3lesize,4,,,5lesize,5,,,2mshape,0mshkey,1lcomb,1,2!lccat,1,2!如采用lcatt,1,2则采用7个划分数lcomb,4,5amesh,all!EX3.7B 简化映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,6mshape,0amap,1,2,5,3,4!EX3.8A 过渡四边形映射网格finish/clear/prep7k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,8lesize,3,,,3lesize,4,,,7lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.8B 过渡四边形映射网格finish/clear/prep7et,1,plane42k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,11lesize,3,,,3lesize,4,,,2lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.9A 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,3lesize,5,,,15allselmshape,0,3dvmesh,all!EX3.9B 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,4lesize,7,,,12mshape,0,3dmshkey,1vmesh,all!EX3.9C 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lsel,s,,,1,8,7lesize,all,,,6lsel,s,,,2,7,5lesize,all,,,3lsel,s,,,4,5lesize,all,,,7lsel,s,,,3,6,3lesize,all,,,2lsel,s,,,9,12,1lesize,all,,,5allselmshape,0,3dmshkey,1vmesh,all!EX3.9D 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,2,,,2lesize,4,,,2lesize,5,,,2lesize,1,,,3lesize,3,,,3lesize,7,,,6lesize,6,,,7lesize,11,,,7lesize,9,,,3lesize,10,,,3lesize,12,,,3allselmshape,0,3dmshkey,1vmesh,all!EX3.10 带两孔的长方体的扫掠网格划分finish/clear/prep7a=10r=2et,1,mesh200,6et,2,solid45blc4,,,2*a,a,acyl4,a/2,a/2,r,,,,awprota,,90cyl4,1.5*a,a/2,r,,,,-avsbv,1,2vsbv,4,3wprota,,,90wpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpcsyswpoff,,a/2,a/2vsbw,allwprota,,90vsbw,allwpcsysesize,1amap,105,15,16,26,63amap,107,16,13,60,26amap,108,13,14,28,60amap,103,15,14,28,63lesize,94,,,4lesize,79,,,5vsel,s,loc,x,0,avsweep,allasel,s,loc,y,aasel,r,loc,x,a,2*alccat,2,45lccat,57,71lccat,68,78lccat,65,72mshape,0,2dmshkey,1amesh,allvsel,s,loc,x,a,2*avsweep,allallsel/view,1,1,2,3Eplot!EX3.11 拉伸类命令生成体单元网格finish/clear/prep7et,1,82et,2,95blc4,,,4,4blc4,6,,4,4blc4,12,,4,4blc4,18,,4,4esize,1amesh,allesize,,8vrotat,1,,,,,,1,4,90vext,2,,,,,10,0.5,0.5voffst,3,10vdrag,4,,,,,,35!EX3.12 圆的网格划分finish/prep7et,1,plane82r0=10cyl4,,,r0cyl4,3*r0,,,,r0 wprota,,90asbw,allwprota,,,90asbw,allwpoff,,,3*r0asbw,allwpcsys,-1asel,s,loc,x,-r0,r0lsla,slesize,all,,,8mshape,0,2dmshkey,1amesh,allallselasel,s,loc,x,2*r0,4*r0 lesize,all,,,8lsel,r,length,,r0lesize,all,,,8,0.1,1 amesh,allallsel!EX3.13 圆环的网格划分finish/clear/prep7et,1,plane82r0=10cyl4,,,r0/3,,r0,90cyl4,2*r0,,r0/10,,r0,90 asel,s,loc,x,-r0,r0lsla,slesize,all,,,8lsel,r,length,,r0*2/3 lesize,all,,,3,,1 mshape,0,2dmshkey,1amesh,allasel,s,loc,x,2*r0,4*r0lesize,5,,,12lesize,7,,,6lsel,s,length,,r0*9/10lesize,all,,,7amesh,all!EX3.14 圆柱面和圆柱体的网格划分!圆柱面finish/clear/prep7r0=10h0=50et,1,shell63cyl4,,,r0adele,1cm,l1cm,linek,50k,51,,,h0l,50,51adrag,l1cm,,,,,,5lsel,s,loc,z,0lesize,all,,,6lsel,s,length,,h0lesize,all,,,8mshape,0,2dmshkey,1amesh,all!圆柱体finish/clear/prep7r0=10h0=50et,1,solid95cyl4,,,r0,,,,h0wprota,,90vsbw,allwprota,,,90vsbw,allmshape,0,3dmshkey,1lsel,s,loc,z,0lsel,s,length,,h0lesize,all,,,8vmesh,all!EX3.15 圆锥的网格划分!扫掠网格划分finish/clear/prep7cone,10,,,15,,90et,1,200,7et,2,95lsel,s,,,5,6lesize,all,,,12,0.5lsel,alllesize,3,,,6mshape,0,2dmshkey,1amesh,3vsweep,1,3,4vsymm,x,allvsymm,y,allvglue,all!六面体映射网格划分finish/clear/prep7cone,10,,,15,,90et,2,95lesize,all,,,8mshape,0,3dmshkey,1vmesh,allvsymm,x,allvsymm,y,allvglue,all!EX3.16 正多边形面的通用网格划分finish/clear/prep7ns=11rpr4,ns,,,10kp0=100+nsk,kp0*if,mod(ns,2),eq,0,then*do,i,1,ns/2l,kp0,2*i-1*enddo*else*do,i,1,nsl,kp0,i*enddo*endiflsel,s,,,ns+1,2*nscm,l1cm,lineallselasbl,1,l1cm*if,mod(ns,2),eq,0,then cmsel,s,l1cmlesize,all,,,8lsel,invelesize,all,,,4*elselesize,all,,,8*endifallselmshape,0mshkey,1amesh,all!EX3.17 球体及球面网格划分finish/clear/prep7r0=10sphere,,r0,,90vsbw,allvdele,2,,,1esize,2et,1,solid95mshape,0,3dmshkey,1vsymm,x,all vsymm,y,all vsymm,z,allvglue,allfinish/clear/prep7r0=10sphere,,r0,,90 vsbw,allvdele,2,,,1vdele,allasel,s,loc,x,0asel,a,loc,y,0asel,a,loc,z,0adele,all,,,1allselesize,2et,1,shell63 mshape,0,2d mshkey,1amesh,allarsym,x,allarsym,y,allarsym,z,all nummrg,all!EX3.18 孔板网格划分finish/clear/prep7a0=300b0=800r0=15blc4,,,a0,b0cyl4,a0/4,b0/8,r0 agen,2,2,,,a0/2 agen,2,2,3,1,,b0/8 agen,2,2,5,1,,b0*5/8 asel,s,,,2,9,1cm,a2cm,areaallselwprota,,-90*do,i,1,5wpoff,,,b0/16asbw,all*enddowpoff,,,b0*5/16*do,i,1,5wpoff,,,b0/16asbw,all*enddowprota,,,90*do,i,1,3wpoff,,,a0/4asbw,all*enddowpcsys,-1numcmp,alllsel,s,radius,,r0lesize,all,,,8lsel,invelesize,all,,,4lsel,allet,1,82mshape,0,2dmshkey,1asel,u,loc,y,b0/16,b0*5/16 asel,u,loc,y,b0*11/16,b0*15/16 lsla,slsel,r,tan1,xlesize,all,50,,,,1amesh,allallselamap,21,30,31,54,62 amap,22,31,32,54,64 amap,32,29,30,62,76 amap,33,29,32,64,76amap,42,34,35,76,82 amap,43,35,36,76,84 amap,41,33,34,53,82 amap,44,33,36,53,84amap,19,22,23,49,64 amap,20,23,24,49,60 amap,9,21,22,64,74 amap,31,21,24,60,74amap,38,26,27,74,84 amap,39,27,28,74,80 amap,37,25,26,50,84 amap,40,25,28,50,80amap,17,14,15,43,59 amap,18,15,16,43,63 amap,7,13,14,59,73 amap,8,13,16,63,73amap,34,18,19,73,79 amap,30,17,18,79,44 amap,35,19,20,73,83 amap,36,17,20,44,83amap,15,6,7,39,63amap,16,7,8,39,58amap,4,6,5,63,71amap,5,5,8,58,71amap,27,10,11,71,83 amap,28,11,12,71,78 amap,26,9,10,83,40 amap,29,9,12,40,78!EX3.19 角支架的网格划分finish/clear/prep7blc4,,,150,50blc4,100,,50,-50cyl4,,25,25cyl4,125,-50,25aadd,allnumcmp,allcyl4,,25,10cyl4,125,-50,10 asel,s,,,2,3cm,a1cm,area asel,allasba,1,a1cm lcomb,1,6lfillt,1,2,20 asbl,4,6adele,1,,,1lfillt,3,4,20al,18,19,20 aadd,all numcmp,all!---------------------- wprota,,90 wpoff,,,-25 asbw,all wpoff,,,75 asbw,all wprota,,,90 asbw,all wpoff,,,125 asbw,all wpcsys,-1 wpoff,25 wprota,,,90 asbw,all kwpave,18 asbw,all kwpave,3 wprota,,90 asbw,all kwpave,21 wprota,,-45 asbw,8 wpcsys,-1!-------------------- et,1,plane82 mshape,0,2d mshkey,1 esize,6lesize,33,,,6 lesize,37,,,6 lesize,42,,,6amesh,3,5,2amap,6,9,10,4,23amap,7,9,12,1,23amesh,11,13,2amesh,1,9,8amesh,2,4,2amap,16,3,18,26,28amap,12,14,15,5,28amap,15,13,14,6,28lcomb,19,22lcomb,27,46amesh,10,14,4!========================= !EX3.20 大板小孔的网格划分A finish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1wprota,,90,45asbw,allet,1,plane82mshape,0,2dmshkey,1lsel,s,loc,x,0,a0/20lesize,all,,,8lsel,alllesize,7,,,12,0.1lesize,8,,,12,0.1lesize,9,,,12,10amesh,allcsys,0arsym,x,allarsym,y,allnummrg,all!============================= !EX3.20 大板小孔的网格划分Bfinish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1et,1,plane82mshape,0,2dmshkey,1lesize,5,,,8lesize,1,,,10lesize,4,,,8,5lesize,6,,,8,5amesh,1lesize,7,,,20,0.1lesize,8,,,20,0.1amap,2,50,51,2,4csys,0arsym,x,allarsym,y,allnummrg,all!EX3.21 Z形面网格划分及局部细分finish/clear/prep7blc4,,,15,10blc4,10,6,14,12aadd,allwprota,,-90wpoff,,,6asbw,allwpoff,,,4wprota,,,90wpoff,,,10asbw,allwpoff,,,5asbw,allwpcsys,-1et,1,82esize,2mshape,0,2dmshkey,1amesh,allkrefine,9,10,1,1,,off !不光滑处理（节点位置不变）aclear,allamesh,allkrefine,9,10,1,1,,smooth !光滑处理，节点位置可改变，但既有单元不删除aclear,allamesh,allkrefine,9,10,1,1 !光滑处理，可删除既有单元重细分，节点位置改变!EX3.22 分布小孔环形面的网格划分finish/clear/prep7r1=15r2=25r3=3n=8cyl4,,,r1,,r2,180/ncyl4,0.5*(r1+r2),,r3asba,1,2kl,2larc,6,8,4,0.5*(r1+r2)asbl,all,4lesize,4,,,4lesize,10,,,4lesize,6,,,8lesize,3,,,6lesize,8,,,4lesize,5,,,8lesize,9,,,4lesize,7,,,4et,1,plane82mshape,0,2damap,2,6,7,4,8amap,1,5,6,8,1!lrefine,5,6,1,1,smootharsym,y,allcsys,1agen,n,all,,,,360/nnummrg,all!EX3.23 同心圆轴体的网格划分finish/clear/prep7cylind,50,0,100,150cylind,40,0,200,250cylind,30,0,100,250cylind,20,0,50,250cylind,10,0,0,250vptn,allwprota,,,90vsbw,allwprota,,90vsbw,allet,1,95csys,1lsel,s,length,,50lesize,all,,,5lsel,invelsel,u,radius,,10,50lesize,all,,,4n1=6n2=6n3=8n4=10n5=12lsel,r,loc,x,0,10lesize,all,,,n1,,1lsel,s,radius,,10lesize,all,,,n1lsel,s,radius,,20lsel,s,radius,,30lesize,all,,,n3lsel,s,radius,,40lesize,all,,,n4lsel,s,radius,,50lesize,all,,,n5allselmshape,0,3dmshkey,1vmesh,all!EX3.24A 圆柱与长方体组合模型的网格划分!采用六面体映射网格划分finish/clear/prep7a0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,allaccat,4,6et,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.24B 圆柱与长方体组合模型的网格划分!采用过渡六面体映射网格划分finish/cleara0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,alllsel,s,radius,,r0lsel,a,length,,r0lesize,all,,,6lsel,s,loc,z,0lsel,a,loc,z,h1lesize,all,,,7asel,s,loc,x,a0/2asel,a,loc,y,a0/2accat,allallselet,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.25-任意圈数圆柱形螺旋弹簧的扫略网格划分FINISH/CLEAR/PREP7!1.定义弹簧参数（同上）D=4C=8DZ=C*DT=DZ/6N=4.7*IF,T,LT,D,THEN*ENDIFTKPD=90!2.创建全部螺旋线CSYS,1TDEG=N*360TDEG1=MOD(TDEG,TKPD) N0=(TDEG-TDEG1)/TKPD+1 *AFUN,DEG*DO,I,1,N0CTA=(I-1)*TKPDZ=T/360*CTAK,I,DZ/2,CTA,Z*ENDDO*IF,TDEG1,LT,1.0E-2,THEN *ELSEN0=N0+1CTA=CTA+TDEG1Z=T/360*CTAK,N0,DZ/2,CTA,Z*ENDIF*DO,I,1,N0-1L,I,I+1*ENDDOCM,L1,LINE!3.在螺旋线端部创建簧丝截面CSYS,0WPOFF,DZ/2 WPROTA,,90CYL4,,,D/2wprota,,90asbw,allwprota,,,90asbw,allcm,a1cm,areaVDRAG,a1cm,,,,,,L1et,1,mesh200,7et,2,solid95cmsel,s,a1cmlsla,slesize,all,,,6amesh,allallselesize,d/2vsweep,all!EX3.26 具坑缺陷圆柱的网格划分finish/clear/prep7l0=50r0=10r1=5r2=1h0=5et,1,solid45cyl4,,,r0,,,,l0wpave,0,r0,l0/2wprota,,90con4,,,r1,r2,h0vsbv,1,2wprota,,90vsbw,allwpoff,,,r1+3vsbw,allwpoff,,,-2*(r1+3)vsbw,allwprota,,,90vsbw,allwpoff,,,r1+1vsbw,allwpoff,,,-2*(r1+1)vsbw,allwpcsys,-1wprota,,90vsbw,allwpoff,,,-r1vsbw,allwpoff,,,2*r1vsbw,allwpcsys,-1numcmp,allksel,s,loc,z,l0/2ksel,r,loc,x,0ksel,u,loc,y,0*get,kp1,kp,0,num,min kp2=kpnext(kp1) ksel,s,loc,z,l0/2+r2*get,kp3,kp,0,num,min allselk,1000,,-r0,l0/2+5l,kp3,1000asel,none*get,l1,line,0,num,max arotat,l1,,,,,,kp1,kp2 cm,a1cm,areaallselvsba,all,a1cm!vsel,s,loc,z,17,33!vsel,r,loc,x,-6,6!aslv,s!aplotaccat,128,61accat,160,138 accat,224,202 accat,183,35accat,31,134accat,151,164 accat,215,231 accat,82,195accat,130,5accat,162,145 accat,229,209 accat,192,89accat,12,123accat,117,155 accat,109,219 accat,14,178mshape,0,3d mshkey,1esize,1vmesh,all!EX3.27长方体开圆柱槽的网格划分finish/clear/prep7blc5,,,20,10,30wpoff,,5,5cyl4,,,8,,,,10wpoff,,,10cyl4,,,4,,,,12vsel,s,,,2,3cm,v1cm,voluallselwpcsys,-1vsbv,1,v1cmvoffst,15,5voffst,16,3voffst,17,15vptn,allwprota,,,90vsbw,allwpcsys,-1wpoff,,,5vsbw,allwpoff,,,10vsbw,allwpoff,,,12vsbw,allet,1,mesh200,7esize,2mshape,0,2dmshkey,1asel,s,loc,z,30lccat,6,74lccat,8,68amesh,allasel,s,loc,z,0lccat,3,73lccat,1,67amesh,allallselet,2,solid95mshape,0,3dmshkey,1vsweep,all!ex4.1累加节点自由度约束值finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,1lmesh,alld,1,ux,1e-2,,,,uyd,1,rotzd,2,alldlistnsel,s,d,u,0,0.1dofsel,s,uydscale,2.5dlistdofsel,alldcum,addd,1,ux,2e-2dlistdcum,ignod,2,uy,1.0dlist!ex4.2 对线施加约束并转换finish/clear/prep7et,1,95blc4,,,10,10,10dl,7,,ux,0.1dl,5,,alldl,11,6,symmdl,10,6,asymdl,6,,symmdllistesize,2vmesh,alldtrandlist!ex4.3 施加集中荷载与节点坐标系finish/clear/prep7et,1,beam4k,1k,2,5k,3,10l,1,2l,2,3local,12,0,,,,90nrotat,alldk,1,allfk,2,fy,-1000esize,1lmesh,allnrotat,alllplotfk,3,fy,1000f,6,fx,-1000sbctraneplot!ex4.4a 3d单元sf加载示例finish/clear/prep7et,1,95blc4,,,10,10,20esize,,4vmesh,allasel,s,loc,y,10sf,all,pres,1000asel,s,loc,z,20nsla,s,1sf,all,pres,1000!ex4.4b 2d单元sf加载示例finish/clear/prep7et,1,82blc4,,,100,200blc4,30,60,40,80asba,1,2wprota,,-90wpoff,,,60asbw,allwpoff,,,80asbw,allwprota,,,90wpoff,,,30asbw,allwpoff,,,40asbw,allwpcsys,-1esize,5amesh,all/psf,pres,norm,2sf,all,pres,100sfdele,all,presnsel,s,loc,x,0sf,all,pres,100nsel,s,loc,x,15,20esln,s,1nsel,r,loc,x,15sf,all,pres,110nsel,s,loc,x,40,60nsel,r,loc,y,10,30esln,s,1sf,all,pres,100lsel,s,loc,x,100nsll,s,1esln,snsel,s,loc,x,95sf,all,pres,-100eplot!ex4.5 节点号及其荷载函数finish/clear/prep7et,1,45blc4,,,10,10,20esize,5vmesh,all*dim,mypres,,100*do,i,1,100mypres(i)=i*10.0*enddosffun,pres,mypres(1)nsel,s,loc,y,10sf,all,pres,10sflist*do,i,1,100mypres(i)=i*50.0*enddonsel,s,loc,z,20sf,all,pres,0allselsflist!ex4.6a 2d平面单元plane82 finish/clear/prep7et,1,82blc4,,,10,50esize,2amesh,allsfe,1,4,pres,,100,50sfe,6,4,pres,,200,40nsel,s,loc,x,10esln,ssfe,all,2,pres,,-100/psf,pres,norm,2,0,1eplot!ex4.6b 3d体单元solid95 finish/clear/prep7et,1,95blc4,,,10,10,30esize,5vmesh,all/psf,pres,norm,2,0,1sfe,22,3,pres,,100sfe,22,6,pres,,100,110,120,130 sfelist!ex4.6c 3d壳单元shell63 finish/clear/prep7et,1,63wprota,,90blc4,,,10,10esize,5amesh,all/psf,pres,norm,2,0,1sfe,1,1,pres,,100sfe,1,3,pres,,100sfe,1,6,pres,,100sfe,4,1,pres,,100,110,120,130!ex4.7 在梁单元上施加荷载finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,,10lmesh,all/pnum,elem,1sfbeam,3,1,pres,50,100 sfbeam,5,1,pres,100 sfbeam,7,1,pres,50,100,,,0.2,0.1sfbeam,9,1,pres,100,,,,0.4,-1sfbeam,3,2,pres,50,100,,,0.2,0.1!ex4.8a 在线上施加面荷载finish/clear/prep7et,1,82blc4,,,10,30esize,5amesh,all/psf,pres,norm,2sfl,4,pres,10,60sfl,2,pres,60sftraneplot!ex4.8b 在线上施加面荷载finish/clear/prep7et,1,shell63wprota,,90blc4,,,10,30esize,5amesh,allsfl,3,pres,100sftraneplot!ex4.9 利用荷载梯度在直角坐标系下的施加方法finish/clear/prep7et,1,82blc4,,,10,60esize,2amesh,all/psf,pres,norm,2sfgrad,pres,,y,0,-5nsel,s,loc,x,0nsel,r,loc,y,0,40sf,all,pres,600sfgrad,pres,,y,30,-20nsel,s,loc,x,10sf,all,pres,0allseleplot!ex4.10利用荷载梯度在柱坐标系下的施加方法finish/clear/prep7csys,1k,1,10,-90k,2,10,90k,3,10,90,30l,1,2l,2,3adrag,1,,,,,,2ldele,2,,,1numcmp,allet,1,63esize,2amesh,allsfgrad,pres,1,y,-90,1sfa,all,2,pres,400sftran!可考察下述命令及其施加方法!①!sfgrad,pres,1,y,270,1sfa,all,2,pres,400sftran!②!local,12,1cscir,12,1sfgrad,pres,12,y,270,1sfa,all,2,pres,400sftran!ex4.11 ldread命令的示例finish/clear/filname,ldtest1/prep7et,1,plane82mp,ex,1,2.1e5mp,prxy,1,0.3blc4,,,10,40esize,2amesh,alld,all,allsfl,3,pres,100/solusolve/post1prrsolfinish/filename,ldtest2/solulsclear,allnsel,s,loc,y,0d,all,allallselldread,reac,,,,,ldtest1,rst solve/post1plnsol,s,y!ex4.12 初应力荷载finish/clear/filname,colu1/prep7et,1,plane82mp,ex,1,2e5mp,nuxy,1,0.3blc4,,,1,10esize,2amesh,allnsel,s,loc,y,0d,all,uyd,1,uxnsel,s,loc,y,10sf,all,pres,-10。

用Ansys或Abaqus分析钢管混凝土结构或构件

用Ansys或Abaqus分析钢管混凝土结构或构件用Ansys或Abaqus分析钢管混凝土结构或构件以上两个软件国外都有人用来分析钢管混凝土结构，但建模的方法不尽相同。

关键在于钢管和混凝土本构关系的选取以及两者之间的界面处理方法，各位有没有这方面的经验能向我们大家介绍一下。

==========程序中大概只有Drucker-Prager比较适合描述受约束混凝土的本构关系，因为这个模型可以考虑 hydrostatic stress （流体静应力）的影响。

在程序中，需要输入cohesion, angle of internal friction，(one more for ANSYS is the angle of dilatancy)。

值得注意的是，两个软件确定这几个参数的公式各不相同，很是令人头疼。

其实user manuals不可能给出明确的表达式，因为到目前为止，好像没有研究把钢管的强度，混凝土的强度，含钢率等等因素（i.e. the confinement）全部在Drucker-Prager 中考虑进去。

至于两种材料的界面，日本的 Hanbin Ge曾用link element来模拟，但在他的文章中，没有详细的描述。

轴压状况下，好像可以忽略滑移。

偏压可能情况有所不同。

==========韩教授书上的混凝土应力-应变关系，可以简单理解为单向受力的混凝土本构关系（考虑了钢管的约束），因此不能用于多向应力状态下混凝土的有限元分析。

材料非线性有限元分析，需要定义材料的屈服面，流动准则，强化准则，等等。

对受约束的混凝土，还要考虑体积膨胀，钢管对它的约束等因素。

显然，不是一个简单的应力-应变曲线所能概括的。

==========三向有限元分析，需要定义屈服面、流动准则和强化准则等等，而考虑钢管约束的混凝土本构关系，只是应力-应变关系。

对钢管混凝土的有限元分析，主要困难是如何定义屈服面，和模拟两个材料之间的滑移，我曾经用过接触分析（contact analysis）来求轴压构件的承载力，发现最大承载力能够比较精确地求得，但是精确的荷载-位移曲线很难获得，因为商用软件（Ansys\Marc）里面的D-P模型是塑性模型。

ANSYS结构非线性分析相应步骤及命令流

ANSYS结构非线性分析相应步骤及命令流屈服准则概念：1.理想弹性材料物体发生弹性变形时，应力与应变完全成线性关系，并可假定它从弹性变形过渡到塑性变形是突然的。

2.理想塑性材料（又称全塑性材料）材料发生塑性变形时不产生硬化的材料，这种材料在进入塑性状态之后，应力不再增加，也即在中性载荷时即可连续产生塑性变形。

3.弹塑性材料在研究材料塑性变形时，需要考虑塑性变形之前的弹性变形的材料这里可分两种情况：Ⅰ.理想弹塑性材料在塑性变形时，需要考虑塑性变形之前的弹性变形，而不考虑硬化的材料，也即材料进入塑性状态后，应力不再增加可连续产生塑性变形。

Ⅱ.弹塑性硬化材料在塑性变形时，既要考虑塑性变形之前的弹性变形，又要考虑加工硬化的材料，这种材料在进入塑性状态后，如应力保持不变，则不能进一步变形。

只有在应力不断增加，也即在加载条件下才能连续产生塑性变形。

4.刚塑性材料在研究塑性变形时不考虑塑性变形之前的弹性变形。

这又可分两种情况：Ⅰ.理想刚塑性材料在研究塑性变形时，既不考虑弹性变形，又不考虑变形过程中的加工硬化的材料。

Ⅱ.刚塑性硬化材料在研究塑性变形时，不考虑塑性变形之前的弹性变形，但需要考虑变形过程中的加工硬化材料。

屈服准则的条件：1.受力物体内质点处于单向应力状态时，只要单向应力大到材料的屈服点时，则该质点开始由弹性状态进入塑性状态，即处于屈服。

2.受力物体内质点处于多向应力状态时，必须同时考虑所有的应力分量。

在一定的变形条件（变形温度、变形速度等）下，只有当各应力分量之间符合一定关系时，质点才开始进入塑性状态，这种关系称为屈服准则，也称塑性条件。

它是描述受力物体中不同应力状态下的质点进入塑性状态并使塑性变形继续进行所必须遵守的力学条件，这种力学条件一般可表示为）＝Cf（σij又称为屈服函数，式中C是与材料性质有关而与应力状态无关的常数，可通过试验求得。

屈服准则是求解塑性成形问题必要的补充方程。

1.1 什么是结构非线性在日常生活中，经常会遇到结构非线性。

用ANSYS建立钢筋混凝土梁模型 (2)

用ANSYS建立钢筋混凝土梁模型问题描述：钢筋混凝土梁在受到中间位移荷载的条件下的变形以及个组成部分的应力情况。

一、用合并节点的方法模拟钢筋混凝土梁1．用solid65号单元以及beam188单元时材料特性钢材的应力应变关系混凝土的弹性模量采用线弹性B=150mmH=300mmL=2000mmP=5mm位移图1钢筋混凝土结构尺寸图建立钢筋线对钢筋线划分网格后形成钢筋单元建立混凝土单元合并单元节点后施加约束以及位移载荷进入求解器进行求解钢筋单元的受力云图混凝土的应力云图混凝土开裂2 使用单元solid45号单元与beam188钢筋的应力应变关系不变，而混凝土应力应变关系为：混凝土单元钢筋单元力与位移曲线合并节点时的命令流:fini/clear,nostart/config,nres,5000/prep7/title,rc-beamb=150h=300a=30l=2000fcu=40ec=2.85e4displacement=10!定义单元类型et,1,solid45et,2,beam188et,3,plane42!定义截面类型sectype,1,beam,csolid,,0 secoffset,centsecdata,8,0,0,0,0,0,0,0,0,0 sectype,2,beam,csolid,,0 secoffset,centsecdata,4,0,0,0,0,0,0,0,0,0!定义材料属性，混凝土材料属性mp,ex,1,ecmp,prxy,1,0.2tb,kinh,1,,16tbpt,,0.000179067,5.10tbpt,,0.000358133,9.67tbpt,,0.0005372,1.37e1tbpt,,0.000716267,1.72e1 tbpt,,0.000895333,2.01e1 tbpt,,0.0010744,2.26e1tbpt,,0.001253467,2.44e1 tbpt,,0.001432533,2.58e1 tbpt,,0.0016116,2.66e1tbpt,,0.001790667,2.69e1 tbpt,,0.0019916,2.65e1tbpt,,0.002393467,2.57e1tbpt,,0.002795333,2.48e1 tbpt,,0.0031972,2.40e1 tbpt,,0.003599067,2.32e1 tbpt,,0.0038,2.28e1tb,conc,1,1,9tbdata,,0.4,1,3,-1!纵向受拉钢筋mp,ex,2,2e5mp,prxy,2,0.3tb,bkin,2,1,2,1tbdata,,350!横向箍筋，受压钢筋材料属性mp,ex,3,2e5mp,prxy,3,0.25tb,bkin,3,1,2,1tbdata,,200!生成钢筋线k,,k,,bkgen,2,1,2,,,hk,,a,ak,,b-a,akgen,2,5,6,,,h-2*akgen,21,5,8,,,,-100*do,i,5,84,1l,i,i+4*enddo*do,i,5,85,4l,i,i+1l,i,i+2*enddo*do,i,8,88,4l,i,i-1l,i,i-2*enddo!受拉钢筋lsel,s,loc,y,alsel,r,loc,x,alsel,a,loc,x,b-alsel,r,loc,y,acm,longitudinal,linetype,2mat,2.secnum,1lesize,all,50lmesh,allallscmsel,u,longitudinalcm,hooping reinforcement,line!箍筋，受压钢筋type,2mat,2secnum,2lesize,all,50lmesh,all/eshape,1!将钢筋节点建为一个集合cm,steel,node!生成面单元，以便拉伸成体单元a,1,2,4,3lsel,s,loc,y,0lsel,a,loc,y,hlesize,all,,,8lsel,alllsel,s,loc,x,0lsel,a,loc,x,blesize,all,,,10type,3amesh,all!拉伸成混凝土单元type,1real,3mat,1extopt,esize,20extopt,aclear,1vext,all,,,,,-lalls!合并节点nummrg,allnumcmp,all!边界条件约束nsel,s,loc,y,0nsel,r,loc,z,0d,all,uyd,all,uxnsel,s,loc,y,0nsel,r,loc,z,-l d,all,uyd,all,ux!施加外部荷载/solunsel,allnsel,s,loc,y,hnsel,r,loc,z,-1000d,all,uy,-displacementalls!求解nlgeom,onnsubst,50outres,all,allneqit,50pred,oncnvtol,f,,0.05,2,0.5allselsolvefinish/post1allsel/device,vector,1!时间历程后处理/post26nsel,s,loc,z,-l/2*get,Nmin,node,0,num,min nsol,2,nmin,u,yprod,3,2,,,,,,-1nsel,s,loc,y,0nsel,r,loc,z,0*get,Nnum,node,0,count *get,Nmin,node,0,num,min n0=Nminrforce,5,Nmin,f,y*do,i,2,ndinqr(1,13)ni=ndnext(n0)rforce,6,ni,f,yadd,5,5,6n0=ni*enddoprod,7,5,,,,,,1/1000/axlab,x,uy/axlab,y,p(kn)xvar,3plvar,7二、用约束方程法模拟钢筋混凝土梁1．用solid65号单元以及beam188单元时混凝土以及钢筋采用线弹性关系：建立钢筋线对钢筋线划分网格后形成钢筋单元建立混凝土单元对钢筋线节点以及混凝土节点之间建立约束方程后施加约束以及位移载荷进入求解器进行求解；钢筋单元的受力云图混凝土的应力云图混凝土开裂2 使用单元solid45号单元与beam188使用混凝土的本构关系曲线钢材的本构关系曲线钢筋的von mises 应力混凝土的应力用在solid45号单元下，用合并节点法、约束方程法建立模中钢筋与混凝土之间的关系的时候的一个力与位移全程曲线的比较。

ANSYS应用实例：钢筋混凝土简支梁数值模拟

（ii ）纵向钢筋：PIPE20 （iii ）横向箍筋：PIPE202.2 材料性质（i ）、混凝土材料表5-4 混凝土材料的输入参数一览表[16~19]·单轴受压应力-应变曲线（εσ-曲线）在ANSYS ○R程序分析中，需要给出混凝土单轴受压下的应力应变曲线。

在本算例中，混凝土单轴受压下的应力应变采用Sargin 和Saenz 模型[17,18]：221⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛-+=c c s c c E E E εεεεεσ （5-30）式中取4'4')108.0028.1(c c c f f -=ε；断面图配筋图断面图配筋图断面图配筋图RCBEAM-01 RCBEAM-02 RCBEAM-03图5-12 各梁FEM模型断面图（a）单元网格图（b）钢筋单元划分图图5-13 算例（一）的FEM模型图2.4 模型求解在ANSYS○R程序中，对于非线性分析，求解步的设置很关键，对计算是否收敛关系很大，对于混凝土非线性有限元分析，在计算时间容许的情况下，较多的求解子步（Substeps）或较小的荷载步和一个非常大的最大子步数更容易导致收敛[2]。

在本算例中，设置了100个子步。

最终本算例收敛成功，在CPU为P41.6G、内存为256MB的微机上计算，耗时约为8小时。

2.5 计算结果及分析2.5.1 荷载—位移曲线图5-14为ANSYS○R程序所得到的各梁的荷载-跨中挠度曲线，从图中可以看出：（i）、梁RCBEAM-01：曲线形状能基本反映钢筋混凝土适筋梁剪切破坏的受力特点，而且荷载-跨中挠度曲线与钢筋混凝土梁的弯剪破坏形态非常类似，即当跨中弯矩最大截面的纵筋屈服后，由于裂缝的开展，压区混凝土的面积逐渐减小，在荷载几乎不增加的情况下，压区混凝土所受的正应力和剪应力还在不断增加，当应力达到混凝土强度极限时，剪切破坏发生，荷载突然降低。

（ii）、梁RCBEAM-02：荷载-跨中挠度曲线与超筋梁的试验荷载-跨中挠度曲线很相似，在荷载达到极限情况下，没有出现屈服平台，而是突然跌落。

钢筋混凝土梁的ansys分析(最新整理)

基于ANSYS 的钢筋混凝土力学分析摘要本文介绍ANSYS 模拟钢筋混凝土梁的过程，讨论了有限元模型的建立以及在ANSYS 中的实现，给出了用分离式配筋方法对混凝土梁的分析的一般过程。

并给出了详细的命令流过程。

并在此基础上对混凝土梁进行了分析，讨论了在力的作用下混凝土梁的塑形变形和裂缝的发展过程。

关键词 Ansys 混凝土梁分离式配筋The analysis of mechanics of a reinforced concrete based on ANSYSAbstract This paper introduces ANSYS simulation of the reinforced concrete beam process, discusses the establishment of the finite element model and the realization, and gives the ANSYS reinforcement method with separate the analysis of concrete beams of the general process. And gives the detailed command flow process. Based on the analysis of concrete beams, and discussed the concrete beam under the action of forces of the body deformation and fracture process.Keywords Ansys concrete beams reinforced separated1 引言由于钢筋混凝上材料性质复杂，使其表现出明显的非线性行为[1]。

长期以来采用线弹性理论的设计方法来研究钢筋混凝上结构的应力或内力，显然不太合理，尽管有此理论是基于人量试验数据上的经验公式，还是不能准确反映混凝上的力学性能，特别是受力复杂的重要结构，必须采用三维钢筋混凝上非线性有限元方法才能很好地掌握其力学性能。

ansys建模实例-混凝土箱梁(命令流)

/TITLE, BRIDGE XX/PREP7! 定义单元类型ET,1,SHELL63ET，2，SOLID45ET,3,LINK8!定义1,钢,2，混凝土，3,无限刚度MP，EX,1,2.1E11MP,DENS,1,7。

8E3MP,PRXY,1，0.3mp,alpx,1,1.0e—5MP，EX,2，3。

5E10MP，DENS，2,2.6E3MP,PRXY，2，0。

1667MP，EX，3，3。

5E15MP,DENS,3，7.8E3MP，PRXY，3,0。

3*set,nz，154!定义箱梁N1R,1,7，7，7,7!定义箱梁N2R,2,8,8,8，8!定义箱梁N3R,3，4，4,4，4!定义混凝土solid65的实常数R, 4!定义钢绞线的实常数R， 5， 140!降温=预加力/线膨胀系数*钢筋截面积*弹性模量tem1=—1190.4e6/（1。

0e-5*2。

11e11）tem2=-1339。

2e6/（1。

0e—5＊2.11e11）！定义纵向梯度参数表＊dim， z， array， nz*SET,z(1)，130,260，371,482，593，702,811,920,1029，1138，1247,1356，1465，1574，1683,1787.19,1891。

38，1995.57＊SET,z(19），2099。

76，2203。

95，2308。

14,2412.33,2516.52，2620.71,2724。

9,2829。

09，2933.28，3037。

47,3141.66,3245.85，＊SET，z(31），3350.04,3454。

23,3558.42,3662.61，3766。

8，3870。

99*SET，z（37),3975。

18,4079.37,4183.56,4287.75，*SET,z（41),4391.94,4496.13,4600。

32,4704。

51，4808.7，4912。

89,5017.08,5121.27，5225。

钢混凝土组合柱子ANSYS分析建模命令流

!注参数命名：h（高度），w（宽度），t（厚度），sec（代表截面信息）!单位：mm!柱截面尺寸HW300(截面高度)*300(翼缘宽度)*10(腹板厚度)*15(翼缘厚度) hw_sec_h = 300hw_sec_w = 300hw_sec_t1 = 10hw_sec_t2 = 15!梁截面尺寸HN250*125*6*9hn_sec_h = 250hn_sec_w = 125hn_sec_t1 = 6hn_sec_t2 = 9!T型肋加强板尺寸HT75*125*6*9ht_sec_h = 75ht_sec_w = 125ht_sec_t1 = 6ht_sec_t2 = 9ht_l = 1000 !长度!混凝土板厚度尺寸con_sec_t = 120 !(修改）!模型整体尺寸mod_h = 3000 !模型的高mod_w1 = 640 !模型横向宽度原3000mod_w2 = 5100 !模型纵向宽度!其他参数rebar_vr = 0.014 !体积率（配筋率）为1.4％（修改）!横向钢筋？？？6@200??nail_dis = 100 !栓钉间距为100mm （修改）gridding = 100 !网格steelbar = 200 !钢筋间隔concrete.mac!-------------------------------------------------------------------!model1,钢筋全部由实常数定义的配筋率!先运行参数定义命令define.mac/FILNAME,model,0/prep7!定义属性et,1,solid65,,,,,,3 !定义混凝土板单元KEYOPT(6) 混凝土非线性解输出控制:3 -- 同时还给出积分点的解et,2,SHELL181 !定义钢梁单元et,3,combin39,,,0 !定义弹簧单元combin39et,4,solid45 !定义垫板单元(增加垫块以消除混凝土单元的应力集中)et,5,link8 !定义link8单元mp,ex,1,3E4 !定义混凝土板的弹性模量N/mm2mp,prxy,1,0.2 !定义混凝土板的泊松比mp,dens,1,2.6e-9 !定义混凝土板的密度mp,ex,2,2.06E5 !定义钢梁的弹性模量mp,prxy,2,0.3 !定义钢梁的泊松比mp,dens,2,7.8e-9mp,ex,3,2.0E5 !定义钢筋的弹性模量mp,prxy,3,0.3 !定义钢筋的泊松比mp,dens,3,7.85e-9mp,ex,4,2.02E5 !定义垫板材料特性mp,prxy,4,0.3tb,miso,1 !混凝土屈服准则*do,i,1,17,1*if,i,LE,7,THENx=0.0001*2*iy=(2.07*(x/1.58e-3)+(3-2*2.07)*(x/1.58e-3)*(x/1.58e-3)+0.07*(x/1.58e-3)*(x/1.58e-3)*(x/1.58e-3))*26.1*ELSExtemp=0.0001*2*i/1.58e-3x=0.0001*2*iy=xtemp*26.1/(1.13*(xtemp-1)*(xtemp-1)+xtemp)*ENDIFtbpt,,x,y*enddotbpt,,0.0002,6!该处将0.0033改为了0.0035tbpl,miso,1tb,concr,1 !混凝土破坏准则tbdata,,0.6,1.0,2.61,-1!混凝土材料的前四个实常数的含义!1 裂缝张开剪力传递系数.!2 裂缝闭合剪力传递系数!3 单轴抗拉强度!4 单轴抗压强度! 当变量3（4）被设为-1时表示混凝土无开裂（无压碎）tb,bkin,2,1,2 !激活钢梁的温度特性？？？mkintbdata,,235,2000tbpl,bkin,2tb,bkin,3,1,2 !激活钢筋的温度特性tbdata,,335,2000 !320屈服应力，10代表屈服后的切线模量tbpl,bkin,3r,1,3,rebar_vr,0,0,0,0 !体积率，方向角3是什么意思？!r,1,3,0.0089,90,90,3,0.007819rmore,0,0 !定义混凝土板实常数（三个方向的钢筋）!定义钢梁腹板厚度r,2,hw_sec_t1,hw_sec_t1,hw_sec_t1,hw_sec_t1 !定义钢梁腹板厚度(柱子)r,3,hn_sec_t1,hn_sec_t1,hn_sec_t1,hn_sec_t1 !定义钢梁腹板厚度(水平梁)!定义钢梁上下翼缘的厚度r,4,hw_sec_t2,hw_sec_t2,hw_sec_t2,hw_sec_t2 !定义钢梁上下翼缘的厚度(柱子) r,5,hn_sec_t2,hn_sec_t2,hn_sec_t2,hn_sec_t2 !定义钢梁上下翼缘的厚度(水平梁)r,6,0,0,0.02,4607.0,0.04,6765.5 !定义弹簧的荷载变形曲线rmore,0.05,7651.1,0.06,8457.1,0.07,9201.7rmore,0.08,9896.9,0.09,10551.2,0.1,11170.8rmore,0.2,16154.6,0.3,19890.4,0.4,22932.2rmore,0.5,25506.4,0.6,27734.7,0.7,29692.3rmore,0.8,31430.3,0.9,32985.4,1,34385.1rmore,2,43095.0,4,100000,!3,46937.7,4,50000,r,7,12 !定义link8单元的面积!*KEYOPT,3,1,0KEYOPT,3,2,0KEYOPT,3,3,2KEYOPT,3,4,0KEYOPT,3,6,0!*!********************************************************************** !开始建立几何模型!********************************************************************** !先建柱子!********************************************************************** wpcsys,-1,0csys,4 !在局部坐标系中建立柱子模型wpoffs,-(mod_w1+hw_sec_w*0.5),-mod_w2*0.5,0wprota,,90 !(positive Y toward Z)wpoffs,,,-hw_sec_h*0.5blc4,0,0,hw_sec_w,mod_h+con_sec_twpoffs,,,hw_sec_hblc4,0,0,hw_sec_w,mod_h+con_sec_twprota,,,90 !(positive Z toward X)wpoffs,,,hw_sec_w*0.5blc4,0,0,hw_sec_h,mod_h+con_sec_taptn,allasel,allwpcsys,-1,0wpoffs,,,mod_h-hn_sec_t2*0.5asbw,allwpoffs,,,-hn_sec_hasbw,all!!!!!纵梁wpcsys,-1,0csys,4 !在局部坐标系中建立纵梁模型asel,invert !反选wpoffs,-mod_w1,-mod_w2*0.5,mod_h-hn_sec_t2*0.5blc4,-hn_sec_w*0.5,-hw_sec_h*0.5,hn_sec_w,(mod_w2+hw_sec_h)*0.5wpoffs,,,-hn_sec_hblc4,-hn_sec_w*0.5,-hw_sec_h*0.5,hn_sec_w,(mod_w2+hw_sec_h)*0.5 wprota,,,-90 !(positive Z toward X)blc4,0,-hw_sec_h*0.5,hn_sec_h,(mod_w2+hw_sec_h)*0.5cm,hn_temp,areawpcsys,-1,0!agen,3,hn_temp,,,mod_w1aptn,allasel,all!!切割梁wpcsys,-1,0wpoffs,-mod_w1,-(mod_w2+hn_sec_w)*0.5,mod_h-hn_sec_t2*0.5 wpro,,90wpro,,,90wpoffs,,,-hn_sec_w*0.5asbw,allwpoffs,,,hn_sec_wasbw,allaptn,all!混凝土板wpcsys,-1,0csys,4 !在局部坐标系中建立柱子模型wpoffs,-mod_w1,-(mod_w2+hw_sec_h)*0.5,mod_hblc4,-mod_w1*0.5,0,mod_w1,(mod_w2+hw_sec_h)*0.5,con_sec_t!切体wpcsys,-1,0wpoffs,-mod_w1,-mod_w2*0.5,mod_h-hn_sec_t2*0.5wpro,,90!vsbw,all!wpoffs,,,-hw_sec_h*0.5!vsbw,allwpro,,,90vsbw,allwpoffs,,,-hn_sec_w*0.5vsbw,allwpoffs,,,hn_sec_wvsbw,all!!!付属性ALLSEL,BELOW,VOLUasel,invertcm,gangjiegou,areawpcsys,-1,0wpoffs,-mod_w1,-(mod_w2+hn_sec_h)*0.5,mod_h-hn_sec_t2*0.5 asel,r,loc,z,wpoffs,,,-hn_sec_hasel,a,loc,z,cm,hnsect2,areaAA TT, 2, 5, 2, 0, !!!!梁翼缘cmsel,s,gangjiegou,areacmsel,u,hnsect2,areawpcsys,-1,0wpoffs,-mod_w1,-mod_w2*0.5,mod_h-hn_sec_t2*0.5asel,r,loc,y,hw_sec_w*0.5+0.001,mod_w2cm,hnsect1,areaAA TT, 2, 3, 2, 0, !!!!!梁腹板cmsel,s,gangjiegou,areacmsel,u,hnsect2,areacmsel,u,hnsect1,areacm,hwall,areaasel,r,loc,x,cm,hwsect1,areaAA TT, 2, 2, 2, 0, !!!!!柱子腹板cmsel,s,hwall,areacmsel,u,hwsect1,areaAA TT, 2, 4, 2, 0, !!!!!柱子翼缘!!!!钢结构网格cmsel,s,gangjiegou,areaESIZE,gridding,0,AMESH,allnummrg,all!!!!!!!选钢筋切体wpcsys,-1,0allsel,allKWPA VE, 33wpoffs,,,-20vsbw,allwpcsys,-1,0KWPA VE, 33wpoffs,20,wpro,,,90vsbw,all*do,i,1,mod_w1/steelbar,1 wpoffs,,,steelbar !!!!!!改了vsbw,all*enddo!!!!!!!!!!!!!!!!!!!!!体模型save.amesh!!!!!!!!!!!!!!!!!!!!!选钢筋的线wpcsys,-1,0allsel,allKWPA VE, 33wpoffs,20,,-20lsel,r,loc,zcm,ltemp,line !选出那一层的线lsel,r,loc,x,0cm,ltemp1,linecmsel,s,ltemp,line*do,i,1,mod_w1/steelbar,1lsel,r,loc,x,i*steelbarcmsel,a,ltemp1,linecm,ltemp1,linecmsel,s,ltemp,line*enddocmsel,s,ltemp1,lineLATT,3,7,5, , , ,ESIZE,gridding,0,LMESH,all!!!!!!!!!!!!!!!!!体单元allsel,allESIZE,gridding,0,V ATT, 1, 1, 1, 0 !体属性vmesh,all!**************************************************** !加弹簧、耦合!上混泥土板!***************************type,3 !设定弹簧单元real,6wpcsys,-1,0wpoffs,-mod_w1,nsel,s,loc,xnsel,r,loc,y,0,-(mod_w2-hw_sec_w)*0.5 !!1400nsel,r,loc,z,mod_h-0.001,mod_h+0.001cm,enod,node*get,max0,node,0,count !max1=15*dim,ojd0,,max0*dim,jd0,,max0*get,nod0,node,0,num,minojd0(1)=nod0*do,i,2,max0ojd0(i)=ndnext(ojd0(i-1))*enddoallsel,allnsel,allcmsel,u,enod*do,i,1,max0nod0=ojd0(i)j0=nnear(nod0)jd0(i)=j0*enddonsel,all*do,i,1,max0e,ojd0(i),jd0(i)*enddoallsel,all!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!node(x,y,z) !返回最近节点的编号!nnear(n) !返回最接近n的节点编号!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!耦合纵梁节点wpcsys,-1,0wpoffs,-mod_w1,!nsel,r,loc,y,0,-2400 !!1400nsel,r,loc,z,mod_h-0.001,mod_h+0.001nsel,r,loc,x,-hn_sec_w*0.5,hn_sec_w*0.5cm,cmljnod,node*get,max1,node,0,count !max1=15*dim,ojd,,max1*dim,jd,,max1*get,nod1,node,0,num,minojd(1)=nod1*do,i,2,max1ojd(i)=ndnext(ojd(i-1))*enddoallsel,allnsel,allcmsel,u,cmljnod*do,i,1,max1nod1=ojd(i)j=nnear(nod1)jd(i)=j*enddonsel,allji = 1*do,i,1,max1cp,ji,ux,ojd(i),jd(i)cp,ji+1,uz,ojd(i),jd(i)ji=ji+2*enddoallsel,all!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!node(x,y,z) !返回最近节点的编号!nnear(n) !返回最接近n的节点编号!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!耦合端部节点wpcsys,-1,0wpoffs,-mod_w1,nsel,s,loc,x,-hn_sec_w*0.5,hn_sec_w*0.5nsel,r,loc,y,-(mod_w2+hw_sec_w)*0.5 ,-(mod_w2-hw_sec_w)*0.5 !!1400nsel,r,loc,z,mod_h-0.001,mod_h+0.001cm,egnod,node*get,max2,node,0,count !max1=15*dim,ojd2,,max2*dim,jd2,,max2*get,nod2,node,0,num,minojd2(1)=nod2*do,i,2,max2ojd2(i)=ndnext(ojd2(i-1))*enddoallsel,allnsel,allcmsel,u,egnod*do,i,1,max2nod2=ojd2(i)j2=nnear(nod2)jd2(i)=j2*enddonsel,allji2 = ji + 1*do,i,1,max2cp,ji2,uy,ojd2(i),jd2(i)ji2=ji2+1*enddoallsel,all!!!!!!!!!!!载荷wpcsys,-1,0nsel,s,loc,z,D,all, , , , , ,ALL, , , , ,!ACEL,0,0,0,nsel,s,loc,y,D,all, , , , , ,UY, , , , ,allsel,allwpcsys,-1,0wpoffs,-mod_w1,-(mod_w2-hw_sec_h)*0.5,mod_h-hn_sec_t2*0.5 nsel,r,loc,z,nsel,r,loc,x,nsel,r,loc,y,0,mod_h-hn_sec_t2*0.5cm,N_load,nodeallsel,allwpcsys,-1,0!cmsel,s,N_load,node!*get,Nnod,node,0,count!*Dim,Nodes,array,Nnod!*get,Nd,node,0,num,min!*do,I,1,Nnod,1! Nodes(I)=Nd! F,Nd,FZ,-11000 !!!!!2000-20000 ! Nd=NDNEXT(Nd)!*ENDDO!!SFTRAN!allsel,all!wpcsys,-1,0!!/sol!cnvtol,f,,0.05,2 !定义收敛条件,使用缺省的V ALUE!!nsubst,50 !定义子步数!outres,all,all !输出每一子步的结果!autots,1 !打开自动时间步控制!lnsrch,1 !打开线性搜索!ncnv,2 !如果不收敛时结束而不退出!!neqit,50 !每一子步中方程的迭代次数限值!pred,on !打开预测器!!ANTYPE,0!NLGEOM,1!NSUBST,100,0,0!OUTRES,ERASE!OUTRES,ALL,LAST!AUTOTS,-1.0!PSTRES,1!TIME,20。

ANSYS简支梁计算命令流

ANSYS简支梁计算命令流finish/clear/PREP7et,1,link8 !定义link8单元et,2,solid65 !定义solid65单元keyopt,2,7,1r,1,314 !定义link8单元的面积r,2 !定义solid65的实常数号mp,ex,1,2e5 !定义link8单元的弹性模量mp,prxy,1,0.3 !定义link8单元的泊松系数mp,ex,2,4e4 !定义solid65单元的弹性模量mp,prxy,2,0.3 !定义solid65单元的泊松系数tb,concr,2 !定义2号为混凝土tbdata,,0.9,1,1.8,50 !定义混凝土的c1,c2,Rl,Ra!----------建立几何模型-------------------------blc4, , ,50,200,1500 !定义梁体/view,1,1,1,1 !定义ISO查看/ang,1vplot !绘制梁体kwpave,5 !工作平面移动到关键点5wpoff,25 !工作平面移动25mmwprot,0,0,90 !工作平面旋转vsbw,1 !分割梁体wpoff,0,40 !工作平面移动40mmwprot,0,90 !工作平面旋转vsbw,all !分割梁体wpoff,300 !再将梁体分割为三个区域wprot,0,0,90 !(原为控制单元密度而设)vsbw,allwpoff,0,0,900vsbw,allwpstyl !关闭工作平面显示nummrg,all !整理编号numcmp,all!------------划分单元网格---------------------------lsel,s,,,41,57,16 !定义line41,line57为新的选择集lsel,a,,,9 !定义line9也在选择集中latt,1,1,1 !定义选择集的属性lesize,9,,,20,,,,,1 !将line9划分为20个单元lesize,41,,,20,,,,,1 !将line41划分为20个单元lesize,57,,,60,,,,,1 !将line57划分为60个单元lmesh,all !对线划分单元!-------------allsel,all !新的选择集为所有的实体vsel,s,,,all !定义所有体为选择集vatt,2,2,2 !定义选择集的属性lesize,22,,,2,,,,,0 !定义体线的单元数以控制体单元lesize,23,,,2,,,,,0lesize,33,,,2,,,,,0lesize,28,,,4,,,,,0mshape,0,3d !将体划分单元的形状定位HEX mshkey,1 !采用MAPPED划分器vmesh,all !对体进行划分单元finish!--------------定义约束及荷载-------------------------/soluasel,s,,,2,22,20asel,a,,,14,17,3da,all,uy !对支撑端面施加Uy方向约束asel,s,,,3,23,20asel,a,,,15,19,4da,all,uz !对正对称面施加Uz方向约束asel,s,,,32,36,4asel,a,,,50,52,2asel,a,,,41,43,2da,all,ux !对侧对称面施加Ux方向约束asel,s,,,21,27,6asel,a,,,42,46,4asel,a,,,35,51,16sfa,all,1,pres,0.04 !对顶面施加面荷载(N/mm*mm) allsel,all!--------------定义非线性分析选项-------------------- cnvtol,f,,0.05,2 !定义收敛条件,使用缺省的VALUE nsubst,25 !定义子步数outres,all,all !输出每一子步的结果autots,1 !打开自动时间步控制lnsrch,1 !打开线性搜索ncnv,2 !如果不收敛时结束而不退出neqit,50 !每一子步中方程的迭代次数限值pred,on !打开预测器time,50 !定义时间solve !求解!--------------进入后处理过程----------------------- finish/post1pldisp,0 !绘制变形图(缺省为最后一步的) etable,sig,ls,1 !定义应力表plls,sigi,sigi,1 !绘制上述应力plnsol,s,z,0,1 !绘制Z方向的应力。

ANSYS土木工程经典实例命令流大全解析

1．ANSYS SOLID65环向布置钢筋的例子 (3)2.混凝土非线性计算实例（1）- MISO单压 (5)3.混凝土非线性计算实例（2）-MISO约束压 (6)4.混凝土非线性计算实例（3）- KINH滞回 (9)5.混凝土非线性计算实例（4）- KINH压-拉裂 (11)6.混凝土非线性计算实例（5） (12)7.混凝土非线性计算实例（6） (14)8.混凝土非线性计算实例（7）-MISO滞回 (16)9.混凝土非线性计算实例（8） (18)10.混凝土非线性计算实例（9）-梁平面应力 (20)11.四层弹簧－质点模型的地震分析 (22)12.悬臂梁地震分析 (48)13.用beam 54单元描述变截面梁的例子 (72)14.变截面梁实例 (73)15.拱桥浇筑过程分析-单元生死应用实例 (74)16.简支梁实体与预应力钢筋分析实例 (75)17. 简单的二维焊接分析－单元生死实例 (77)18.隧道开挖（三维）的命令流 (84)19.岩土接触分析实例 (101)20.钢筋混凝土管的动力响应特性分析实例 (109)21.隧道模拟开挖命令流（入门） (116)22.螺栓连接的模拟实现问题 (119)23.道路的基层、垫层模量与应力之间的关系 (129)23.滞回分析 (151)24.模拟某楼层浇注 (153)25.在面上施加移动的面力 (155)27.在任意面施加任意方向任意变化的压力 (159)28.预紧分析 (160)29.几何非线性+塑性+接触+蠕变 (162)30.埋设在地下的排水管道 (167)32.幕墙企业玻璃简化计算 (172)33.等截面杆单元生死应用实例 (188)34.梁板建模联系 (189)36.简单的例子－如何对结构的振动控制分析 (192)37.模态分析结果的输出实例 (194)38.火车过桥动态加载实例（部分） (196)39.悬索结构的找形和计算的例题 (213)40.陶瓷杆撞击铝板的例子 (218)41.求反作用力的APDL命令法 (221)42.LS-DYNA实例（部分） (222)43.路面分层填筑对路基的影响 (223)44.一个例子（含地震影响，求振兴与频率） (227)45.接触面上的压力总和 (231)46.施加位置函数荷载 (235)47.非线性分析考虑刚度退化 (236)48.一个圆形水池的静力分析 (237)49.ANSYS中混凝土模式预应力模拟的算例 (238)50.悬臂梁受重力作用发生大变形求其固有频率 (240)51.循环对称结构模态分析 (242)52.三角平台受谐波载荷作用的结构响应 (244)53.三角平台受一地震谱激励的应力分布和支反力 (246)54.三角平台受时程载荷作用的应力分布和变形过程 (248)55.经典层合板理论 (250)56．定易圆轨迹的例子 (257)57．模拟门式刚架施工－单元生死 (257)58．钢筋混凝土整体式模型例子 (260)59．在荷载步之间改变材料属性例子 (262)60．含预应力的特征值屈曲计算 (263)61．振型叠加计算及工况组合例子 (265)62．柱子稳定分析算（预应力，特征值屈曲，初始缺陷） (268)63．module MConcrete !混凝土模板 (271)64.混凝土开裂实例 (279)65．螺栓网格划分 (280)66．自由液面的土石坝平面渗流分析 (281)67．导出刚度矩阵 (285)68．某混凝土拱坝工程施工期及运行期温度场仿真分析 (286)69．移动温度荷载计算 (293)70．SHSD用于壳-实体装配实例An (295)71．ansys显示－隐式－回弹分析实例 (299)72．工况组合的经典例子 (314)1．ANSYS SOLID65环向布置钢筋的例子! 一个管道，环向配筋率为1%，纵向配筋率为0.5%，径向配筋率为0.1%!FINISH/CLEAR/PREP7!*! 单元属性ET,1,SOLID65!*KEYOPT,1,1,0KEYOPT,1,5,0KEYOPT,1,6,0KEYOPT,1,7,1!*!实参数1：不同方向配筋R,1,2,.001, , ,2, .01,RMORE, 90, ,2,.005 ,90 ,90 ,!Adds real constants to a set.!材料属性!混凝土基本材料属性MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,30E2MPDATA,PRXY,1,,.2!屈服准则TB,MISO,1,1,5,TBTEMP,0TBPT,,0.0005,15TBPT,,0.001,21TBPT,,0.0015,24TBPT,,0.002,27TBPT,,0.003,24!破坏准则TB,CONC,1,1,9,TBTEMP,0TBDATA,,.5,.9,3,30,,!抗拉和抗压是10 倍的关系TBDATA,,,,1,,,MPTEMP,,,,,,,,MPTEMP,1,0!钢材基本属性MPDATA,EX,2,,200E3MPDATA,PRXY,2,,.27!屈服准则TB,BISO,2,1,2,TBTEMP,0TBDATA,,310,2E3,,,,! 管道内径和外径CYL4,0,0,3000, , , ,10000CYL4,0,0,2000, , , ,10000VSBV, 1, 2! 定义局部柱坐标CSWPLA,11,1,1,1,! Defines a local coordinate system at the origin of the working plane. KWPAVE, 11!move to keypoint11WPRO,,-90.000000,VSBW, 3WPCSYS,-1,0!Defines the working plane location based on a coordinate system. KWPAVE, 1WPRO,,,-90.000000VSBW,ALLESIZE,500,0,! 注意：设定单元局部坐标VATT, 1, 1, 1, 11!*VSWEEP,ALL/DEVICE,VECTOR,1/ESHAPE,1.0!Displays elements with shapes determined from the real constants or section/REPLO! 注意：红色代表最大配筋方向，绿色代表其次，蓝色表示最小配筋方向2.混凝土非线性计算实例（1）- MISO单压!MISO单压FINISH/CLEAR/PREP7ET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8TBPT,,3200E-6,14.7TBPT,,3800E-6,12.3TBPT,,4600E-6,9.9!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50MSHAPE,0,3DMSHKEY,1ESIZE,10VMESH,ALLNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUANTYPE,STATICAUTOTS,OFFOUTRES,ALL,1TIME,50NSUBST,50D,1,UY,-.25NSEL,ALLSOLVEFINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH3.混凝土非线性计算实例（2）-MISO约束压!MISO约束压FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8TBPT,,3200E-6,14.7TBPT,,3800E-6,12.3TBPT,,4600E-6,9.9!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLD,ALL,UX,0D,ALL,UZ,0NSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,1TIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,30NSUBST,20,0,20D,1,UY,-.06LSWRITE,2TIME,40NSUBST,100D,1,UY,-.1LSWRITE,3LSSOLVE,1,3FINISH/POST26NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50 ADD,5,3,,,STRESS,,,-1/2500 /AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH4.混凝土非线性计算实例（3）- KINH滞回!KINH滞回FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,KINH,1,1,10TBPT,, 150E-6, 4.5TBPT,, 600E-6,14.8TBPT,,1000E-6,20.25TBPT,,1300E-6,22.3TBPT,,1480E-6,22.8TBPT,,1620E-6,22.8TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,3500E-6,12.8TBPT,,5000E-6, 9.0BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALLTIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,35NSUBST,15,0,15D,1,UY,-.01LSWRITE,2TIME,90NSUBST,55,0,55D,1,UY,-.12LSWRITE,3TIME,125NSUBST,35,0,35D,1,UY,-.05LSWRITE,4TIME,175NSUBST,50,0,50D,1,UY,-.15LSWRITE,5LSSOLVE,1,5FINISH/POST26NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50 ADD,5,3,,,STRESS,,,-1/2500 /AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH5.混凝土非线性计算实例（4）- KINH压-拉裂!KINH压——拉裂压到峰值，泄载、反向加载到拉裂FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,KINH,1,1,10TBPT,, 150E-6, 4.5TBPT,, 600E-6,14.8TBPT,,1000E-6,20.25TBPT,,1300E-6,22.3TBPT,,1480E-6,22.8TBPT,,1620E-6,22.8TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,3500E-6,12.8TBPT,,5000E-6, 9.0BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALLTIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,45NSUBST,25,0,25D,1,UY,.01LSWRITE,2LSSOLVE,1,2FINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH6.混凝土非线性计算实例（5）FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1 TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1 TB,KINH,1,1,18 TBPT,, 100E-6, 3.0 TBPT,, 300E-6, 8.3 TBPT,, 600E-6,14.6 TBPT,, 900E-6,19.1 TBPT,,1100E-6,21.0 TBPT,,1250E-6,22.0 TBPT,,1400E-6,22.6 TBPT,,1550E-6,22.8 TBPT,,1650E-6,22.7 TBPT,,1800E-6,22.3 TBPT,,2000E-6,21.4 TBPT,,2400E-6,19.1 TBPT,,2800E-6,16.8 TBPT,,3200E-6,14.7 TBPT,,3600E-6,13.0 TBPT,,4100E-6,11.3 TBPT,,4600E-6,9.9BLOCK,0,50,0,50,0,50 ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALLTIME,50NSUBST,50,0,50D,1,UY,-.25SOLVEFINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH7.混凝土非线性计算实例（6）!MISO约束压FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,KINH,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8TBPT,,3200E-6,14.7TBPT,,3800E-6,12.3TBPT,,4600E-6,9.9!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLD,ALL,UX,0D,ALL,UZ,0NSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,1TIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,30NSUBST,20,0,20D,1,UY,-.06LSWRITE,2TIME,40NSUBST,100D,1,UY,-.1LSWRITE,3LSSOLVE,1,3FINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH8.混凝土非线性计算实例（7）-MISO滞回!MISO滞回FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15 TBPT,, 100E-6, 3.0 TBPT,, 300E-6, 8.3 TBPT,, 600E-6,14.6 TBPT,, 900E-6,19.1 TBPT,,1100E-6,21.0 TBPT,,1250E-6,22.0 TBPT,,1400E-6,22.6 TBPT,,1550E-6,22.8 TBPT,,1650E-6,22.7 TBPT,,1800E-6,22.3 TBPT,,2000E-6,21.4 TBPT,,2800E-6,16.8 TBPT,,3200E-6,14.7 TBPT,,3800E-6,12.3 TBPT,,4600E-6,9.9BLOCK,0,50,0,50,0,50 ESIZE,10VMESH,ALL MSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50 CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALL TIME,20NSUBST,20,0,20D,1,UY,-.04 LSWRITE,1TIME,35NSUBST,15,0,15D,1,UY,-.01LSWRITE,2TIME,60NSUBST,25,0,25D,1,UY,-.06LSWRITE,3TIME,80NSUBST,20,0,20D,1,UY,-.02LSWRITE,4TIME,110NSUBST,30,0,30D,1,UY,-.08LSWRITE,5LSSOLVE,1,5FINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH9.混凝土非线性计算实例（8）!MISO单压FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,10TBPT,, 150E-6, 4.5TBPT,, 600E-6,14.8TBPT,,1000E-6,20.25TBPT,,1300E-6,22.3TBPT,,1480E-6,22.8TBPT,,1620E-6,22.8TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,3500E-6,12.8TBPT,,5000E-6, 9.0!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,1TIME,50NSUBST,100,0,100D,1,UY,-.25SOLVEFINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH10.混凝土非线性计算实例（9）-梁平面应力!梁平面应力FINISH/CLEAR/PREP7ET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8 TBPT,,3200E-6,14.7 TBPT,,3800E-6,12.3 TBPT,,4600E-6, 9.9ET,2,LINK8R,2,30R,3,6MP,EX,2,200E3MP,NUXY,2,0.3TB,BKIN,2,1TBTEMP,0TBDATA,,300,0/ESHAPE,1/TRIAD,LBOT/VIEW,1,1,1,1N,1,0,0,0N,2,20,0,0NGEN,11,2,1,2,1,0,20 NGEN,43,22,1,22,1,0,0,20 E,1,2,4,3,23,24,26,25 EGEN,10,2,1TYPE,2REAL,2MAT,2E,3,25E,4,26REAL,3E,19,41E,20,42EGEN,42,22,1,14E,47,63E,48,64EGEN,10,88,589,590/NUMBER,1/PNUM,TYPE,1/REPLOTNSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,0NSEL,R,LOC,Z,740D,ALL,UY,0NSEL,ALLFINISH/SOLUANTYPE,STATICNSEL,S,LOC,Y,200NSEL,R,LOC,Z,200D,ALL,UY,-2TIME,200NSUBST,10LSWRITE,1D,ALL,UY,-0.21TIME,21NSUBST,100LSWRITE,2D,ALL,UY,-1TIME,100NSUBST,100LSWRITE,3NSEL,ALLAUTOTS,1OUTRES,ALL,1PRED,ONLSSOLVE,1,3FINISH11.四层弹簧－质点模型的地震分析包括频率分析，谱分析和时程分析作者：陆新征，清华大学土木工程系FINI/CLEAR/UNITS,SI/PREP7DAMPRATIO=0.02 !振型阻尼系数都为0.02ET,1,BEAM188KEYOPT,1,7,1KEYOPT,1,8,1ET,2,MASS21,,,4MP,EX,1,210E9MP,NUXY,1,0.3MP,DENS,1,7800MPDATA,DAMP,1,,DAMPRATIO SECTYPE,3,BEAM,ISECDATA,0.2,0.2,0.5,0.01,0.01,0.008 R,1,1.6E2,1.6E2R,2,1.2E2,1.2E2K,1,0,0,0K,2,0,3,0K,3,0,6,0K,4,0,9,0K,5,0,12,0K,6,0,0,10000L,1,2L,2,3L,3,4L,4,5LSEL,,,,1,4LATT,1,0,1,,6,,3ALLSLESIZE,ALL, , ,3 , , , , ,1 LMESH,ALLNUMMRG,ALLNUMCMP,ALLTYPE,2REAL,1E,2TYPE,2REAL,2E,8E,14E,20*SET,NT,1001*SET,DT,0.02*DIM,AC,,NT*VREAD,AC(1),RECORD,TXT(F8.3)/SOLUD,1,ALL,ANTYPE,0SOLVEFINISH/SOLU !模态分析ANTYPE,2MODOPT,SUBSP,8MXPAND,8, , ,1SOLVEFINI!得到自振频率1*GET,FREQ1,MODE,1,FREQ/ESHAPE,1.0/SOLUANTYPE,SPECTR !谱分析SPOPT,SPRSMDAMP,1,DAMPRATIO, , , , , ,SED,1SVTYPE,2,1FREQ,0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667SV,,0.154350625,0.191590625,0.210210625,0.224175625,0.250716714,0.329 514922,0.484352764,0.598723486,0.786897371FREQ,2.222,10,100000SV,,1.034212766,1.034212766,0.3528!按7度多遇地震，地震影响系数0.08 第一组III场地TG＝0.45SRSS,0.01,DISPSOLVEFINISH/POST1/INP,,MCOMNSORT,U,SUM*GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形!开始地震计算FINISH/CONFIG,NRES,20000/SOLUANTYPE,TRANSTRNOPT,FULLALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)*DO,I,1,1001ACEL,AC(I),0,0TIME,I*0.02OUTRES,ALL,ALLSOLVE*ENDDOFINISH/POST26FILENSOL,2,20,U,X,PLVAR,2, , , , , , , , , ,STORE,MERGEXVAR,1PLVAR,2,*GET,DMAX_TRANS,VARI,2,EXTREM,VMAX, , ! 时程分析最大位移数据文件record.txt-0.0014-0.0108-0.0101-0.0088-0.0095-0.012-0.0142-0.0128-0.011-0.0085-0.0085-0.0131-0.0176-0.0194-0.0162-0.0144-0.0108-0.0082-0.0042-0.0066-0.0131-0.019-0.0196-0.00660.0030.0141-0.0049-0.0144 -0.0203 -0.026 -0.0325 -0.0306 -0.0172 -0.0197 -0.0163 -0.0164 -0.0067 0.0025 0.015 0.0236 0.0252 0.0336 0.0463 0.0492 0.0419 0.0359 0.0271 0.0235 0.0339 0.0412 0.053 0.0639 0.0732 0.0652 0.0599 0.04 0.04 0.0063 -0.0515 -0.0787 -0.0603 -0.0484 -0.025 -0.0059 0.0134 0.0308 0.0499 0.071 0.0995 0.12190.1449 0.1155 0.0935 0.0892 0.0926 0.0839 0.0901 0.0993 0.1209 0.0328 -0.1475 -0.2066 -0.1989 -0.2034 -0.1816 -0.1725 -0.1752 -0.1753 -0.1805 -0.163 -0.1347 -0.1087 -0.0782 -0.0429 -0.0017 0.036 0.0785 0.1164 0.1598 0.196 0.2412 0.2729 0.3036 0.32 0.3417 0.2821 0.2324 -0.1198 -0.2373 -0.164 -0.1865 -0.1095 -0.07530.0113 0.0533 0.0895 0.1186 0.1757 0.0576 -0.2631 -0.1547 -0.1729 -0.1012 -0.0579 0.0237 -0.067 -0.198 -0.1641 -0.1685 -0.1481 -0.1231 -0.1001 -0.0751 -0.0523 -0.0271 -0.0044 0.0188 -0.0095 -0.0433 -0.0838 -0.0951 -0.0716 -0.0599 -0.0334 -0.0108 0.0185 0.042 0.0673 -0.0097 -0.0372 -0.004 0.0011 0.0344 0.0565 0.0883 0.1130.0219 0.0241 0.0683 0.0689 0.1318 0.1353 0.204 -0.0931 -0.1308 -0.0692 -0.0546 0.0072 0.0675 -0.1067 -0.1488 -0.1071 -0.1162 -0.0762 -0.0559 -0.0215 -0.0126 -0.0674 -0.0324 -0.0337 -0.0109 0.0017 0.0299 0.0488 0.0608 0.0222 -0.0032 -0.0245 0.0077 0.0211 0.0568 0.0826 0.1206 0.1478 0.1737 0.0421 0.0029 0.0259 0.0293-0.0147 0.0143 0.0206 0.0499 0.0645 0.0957 0.1128 0.1447 0.1629 0.1945 0.1856 0.1984 0.1769 0.125 -0.1207 -0.0542 -0.0384 -0.0311 -0.1118 -0.1661 -0.2464 -0.2025 -0.1835 -0.1317 -0.096 -0.0325 0.0154 0.0816 0.1319 0.1818 -0.0058 -0.0169 0.0285 0.0447 0.0983 0.1424 0.1853 0.2456 0.1685 -0.138 -0.0999 -0.1089 -0.0907-0.125 -0.2111 -0.1617 -0.1692 -0.1306 -0.1111 -0.0773 -0.051 -0.0544 -0.12 -0.1209 -0.1158 -0.1145 -0.0717 -0.0546 0.0064 -0.0804 -0.1634 -0.0859 -0.0961 -0.0396 -0.0147 0.0319 0.0648 0.0876 0.0472 0.0198 -0.0027 0.0292 0.0445 0.0785 0.1033 0.1352 0.1606 0.1861 0.1281 0.064 0.0204 0.0314 0.0373 0.0496 0.0235 -0.0084-0.0113 -0.0229 -0.0248 -0.0157 -0.0069 0.0147 0.0379 0.0579 0.0255 -0.0041 -0.0428 -0.0133 0.0095 0.023 -0.0129 -0.005 0.008 0.021 0.038 0.051 0.0157 -0.0032 -0.0111 0.0005 0.0076 0.0035 -0.0095 -0.0036 -0.0016 0.0038 0.0085 -0.0056 -0.0304 -0.0421 -0.0244 -0.0236 -0.0177 -0.0129 -0.0018 0.0203 -0.0108 -0.0091 -0.0034-0.0111 -0.0099 -0.0002 0.0073 0.0235 0.0355 0.0705 0.0779 0.0184 -0.0263 -0.0124 -0.0042 0.0159 0.0048 -0.0219 -0.0467 -0.0428 -0.0216 -0.0043 0.0159 0.032 0.0419 0.0123 -0.016 -0.0204 -0.0082 -0.0206 -0.0137 -0.0055 0.0053 0.0134 0.0266 0.0232 0.0079 -0.0008 0.02 0.0435 0.0492 0.0191 0.0092 -0.0022 -0.0021 0.00520.0255 0.0368 0.0525 0.0541 0.0425 0.0398 0.0559 0.0756 0.0365 0.0411 0.0098 -0.0204 -0.0249 -0.0405 -0.0413 -0.0471 -0.0433 -0.0458 -0.0057 0.0178 -0.0208 -0.0492 -0.053 -0.0362 -0.0405 -0.0308 -0.0316 -0.0265 -0.0265 -0.0269 -0.0345 -0.0309 -0.0217 -0.0078 0.0087 0.0281 0.031 0.0358 0.0341 0.0358 0.0287 0.0305 0.01120.0136 0.0384 -0.0861 -0.1349 -0.1342 -0.1354 -0.1193 -0.1042 -0.0829 -0.0651 -0.0444 -0.0258 -0.006 -0.0091 -0.0182 -0.0147 0.0085 0.0163 0.005 0.0264 0.0582 0.0867 0.12 0.1695 0.1111 -0.11 -0.0366 -0.0445 -0.0236 -0.096 -0.0656 -0.0597 -0.067 -0.0552 -0.0027 0.0378 0.1072 0.1669 0.0947 0.0408 0.0667 0.0132 -0.0095-0.0827 -0.1152 -0.115 -0.0803 -0.0369 0.0029 0.0545 0.1178 0.161 -0.027 0.0034 -0.0056 0.002 0.0146 0.0537 0.0798 -0.0205 -0.059 -0.0169 -0.0175 -0.0028 0.0074 0.0382 0.0567 0.0753 0.0801 0.0592 0.0304 0.0023 0.0064 -0.0406 -0.0451 -0.0079 0.0168 0.0567 0.0093 -0.0055 0.0044 -0.0123 -0.0282 -0.0437 -0.0352 -0.02550.0205 0.0519 0.0854 0.1144 0.0733 0.0237 -0.0368 -0.0271 -0.0217 -0.0873 -0.0973 -0.0589 -0.0336 0.0077 0.0259 0.0508 0.0361 0.0081 -0.0056 -0.0209 -0.0317 -0.0238 -0.0376 -0.055 -0.0722 -0.0803 -0.0523 -0.034 -0.0011 0.0065 -0.0037 -0.0005 -0.0168 -0.041 -0.008 0.0079 0.0374 0.0615 0.0665 0.0254 -0.0057 -0.0474 -0.0356-0.0048 0.0126 0.0379 0.0241 -0.0227 -0.0428 -0.0679 -0.0661 -0.059 -0.0513 -0.0408 -0.0309 -0.0266 -0.0541 -0.0628 -0.0908 -0.1107 -0.0881 -0.077 -0.0582 -0.0473 -0.0333 -0.0199 0.002 0.0211 0.0432 0.0613 0.0767 0.0933 0.1066 0.113 0.1187 0.1247 0.1334 0.1594 0.1797 0.2037 0.1236 0.0442 -0.014 -0.0666 -0.0555 -0.0693-0.1246 -0.1179 -0.105 -0.092 -0.0743 -0.0809 -0.085 -0.086 -0.0863 -0.0873 -0.0868 -0.0885 -0.0537 0.0052 0.0215 0.0245 0.058 0.0314 0.0236 0.0485 0.0589 0.0525 0.0355 0.0197 0.0199 0.0492 0.0343 0.0288 0.0432 0.0239 0.0088 0.0077 -0.0148 -0.0077 -0.0019 0.0075 0.0044 -0.0145 -0.0316 -0.0241 -0.0028 0.0182 0.04260.0512 0.0466 0.0479 0.0193 0.0222 0.0274 0.0393 0.0504 0.0577 0.0588 0.0822 0.0797 0.0949 0.0345 0.0045 -0.0123 -0.0347 -0.0426 -0.0416 -0.0275 -0.027 0.0074 0.0428 -0.0231 -0.0387 -0.0083 0.0139 0.0445 0.0027 -0.0697 -0.0796 -0.0251 -0.0135 0.0079 -0.0115 -0.0251 -0.0333 -0.0269 -0.0301 -0.02 -0.0067 -0.0038 0.01050.0344 0.0957 0.0898 0.0179 -0.0362 -0.0994 -0.0807 -0.0744 -0.0539 -0.033 -0.0128 0.0031 0.0148 0.0508 -0.0022 -0.0489 -0.0358 -0.0691 -0.0516 -0.0371 0.0088 0.0632 0.0841 0.1276 0.1388 0.1193 0.0751 0.0225 -0.0088 -0.0227 0.0074 0.0181 0.0544 0.0399 0.0045 -0.0082 -0.0185 -0.002 0.0006 -0.0117 -0.021 -0.0303 -0.0512-0.0579 -0.0266 -0.0178 0.004 0.0098 0.0137 0.0221 0.0437 0.0091 -0.0548 -0.0555 -0.0243 -0.0081 0.025 0.041 0.0182 -0.0027 -0.0243 -0.0015 0.0247 0.0482 0.0783 0.0622 0.0331 -0.0014 -0.0195 -0.0247 -0.0212 -0.011 0.005 0.0241 -0.0034 -0.0216 -0.0471 -0.0363 -0.0195 -0.0018 0.017 -0.008 0.0005 0.023 0.0374 0.06010.0432 0.0344 0.0505 0.0653 0.0683 0.0172 -0.017 -0.0527 -0.0664 -0.0387 -0.0222 -0.0033 0.0119 -0.0128 -0.0351 -0.0514 -0.0335 -0.0218 -0.0012 0.0142 0.007 -0.0063 -0.012 -0.0322 -0.0346 -0.0091 0.0073 0.0309 0.0472 0.0603 0.0576 0.033 -0.0073 -0.0777 -0.0608 -0.0438 -0.0209 0.0031 0.035 0.0293 0.0121 0.0338 0.03170.0206 0.0198 0.0174 0.0021 -0.0144 -0.0343 -0.0339 -0.0145 -0.0028 0.017 -0.0096 -0.0255 -0.0279 -0.0388 -0.0242 -0.0215 -0.0182 -0.0174 -0.0038 -0.0027 -0.0185 -0.0123 0.0087 0.0343 0.0695 0.091 0.0853 0.076 0.0513 0.0186 0.0015 -0.019 -0.0151 -0.0073 0.0021 0.0129 0.0215 0.0024 -0.0124 -0.0329 -0.0519 -0.0708 -0.0579-0.0307 -0.0145 -0.0009 -0.018 -0.0318 -0.0465 -0.0391 -0.0345 -0.0316 -0.0435 -0.0491 -0.0475 -0.042 -0.0361 -0.0277 -0.0258 -0.0139 -0.0068 0.0507 0.0722 0.0878 0.0782 0.0765 0.0439 0.008 0.0013 -0.0126 -0.0015 0.003 0.0104 0.0104 0.0193 0.0205 0.0074 -0.0056 -0.0072 0.007 0.0106 0.0147 -0.0009 -0.0159 -0.0187 -0.00070.0105 -0.0115 -0.0302 -0.0309 -0.0095 -0.0058 0.0004 0.002 0.005 0.0057 0.0097 0.0134 0.0177 0.0218 0.0261 0.0302 0.0346 0.0386 0.0474 0.0393 0.0238 0.0115 -0.0079 -0.0124 0.0054 0.0027 -0.025 -0.0566 -0.063 -0.0591 -0.0413 -0.0068 0.0272 0.0277 -0.0021 -0.006 -0.011 -0.0221 -0.0416 -0.0519 -0.0222 0.003 0.00790.0171 0.0253 0.0323 0.0391 0.0164 -0.0136 -0.0323 -0.0291 -0.0287 -0.0304 -0.0339 -0.0245 -0.0076 0.0125 0.0376 0.0402 0.0245 0.0156 -0.004 -0.0153 -0.0289 -0.0316 -0.0111 0.0094 0.0335 0.0576 0.0424 0.0143 -0.0007 -0.0135 -0.027 -0.0341 -0.0357 -0.0396 -0.0402 -0.0488 -0.048 -0.0406 -0.0407 -0.0351 -0.0187 -0.0057 0.0044-0.0072-0.0169-0.01150.0126finish/clear12.悬臂梁地震分析!定义单元类型/prep7ET,1,BEAM3R,1,0.09,0.000675,0.3, !0.3*0.3m截面MP,EX,1,2.06e11MP,NUXY,1,0.3MP,DAMP,1,0.02 !阻尼比0.02MP,dens,1,7.85e3!--------------------------------------------------!定义节点单元和位移约束*do,i,1,11n,i,(i-1)*0.5,0,0 !5m长,分10段*enddo*do,i,1,10e,i,i+1*enddod,1,ALL,0!--------------------------------------------------!定义和读入时程曲线NT=1000 !时程曲线有NT个点DT=0.01 !时间间隔*dim,ac,,NT !读入数据,这个数据文件可以用excel等软件来写/input,tianjin,txt !天津波东西向共10秒!--------------------------------------------------/SOLUNSUBST,1, , ,1 !1个子步OUTRES,ALL,1 !输出每个子步的结果ANTYPE,TRANS !时程分析*do,i,1,NTACEL,0,ac,0solve*enddo!--------------------------------------------------/POST26NSOL,2,11,U,Y,PLVAR,2,3, , ,下面是文件tianjin.txt，这个文件可以用excel等软件制作或编写很短的程序就能写出来先贴前500个点ac( 1)= -0.06334598ac( 2)= -0.04417088ac( 3)= -0.02188456ac( 4)= -0.00622243ac( 5)= 0.01599961ac( 6)= 0.03431334ac( 7)= 0.05332774ac( 8)= 0.07603510ac( 9)= 0.09322070ac( 10)= 0.11568701ac( 11)= 0.13706695ac( 12)= 0.15446468ac( 13)= 0.17966937ac( 14)= 0.16881546ac( 15)= 0.10252541ac( 16)= 0.03081323ac( 17)= -0.02835768ac( 18)= -0.04577148ac( 19)= 0.00371867ac( 20)= 0.05160179ac( 21)= 0.08542017ac( 22)= 0.14456214ac( 23)= 0.15453540ac( 24)= 0.07974091ac( 25)= -0.00708701ac( 26)= -0.08872422ac( 27)= -0.15415934ac( 28)= -0.16546641ac( 29)= -0.15415934ac( 30)= -0.15538390ac( 31)= -0.15533569ac( 32)= -0.15030567ac( 33)= -0.14931896ac( 34)= -0.14984284ac( 35)= -0.14375542 ac( 36)= -0.14322831 ac( 37)= -0.14130631 ac( 38)= -0.07721145 ac( 39)= 0.02666387 ac( 40)= 0.06092257 ac( 41)= 0.03247169 ac( 42)= 0.01391368 ac( 43)= -0.00049497 ac( 44)= -0.02224453 ac( 45)= -0.04070290 ac( 46)= -0.05431446 ac( 47)= -0.05096541 ac( 48)= -0.03585931 ac( 49)= -0.02533004 ac( 50)= -0.01427366 ac( 51)= -0.00180630 ac( 52)= 0.00995074 ac( 53)= 0.02199384 ac( 54)= 0.03315949 ac( 55)= 0.04807596 ac( 56)= 0.05553902 ac( 57)= 0.03931121 ac( 58)= -0.00056568 ac( 59)= -0.05451052 ac( 60)= -0.07338993 ac( 61)= -0.02214490 ac( 62)= 0.05202605 ac( 63)= 0.11738081 ac( 64)= 0.17217417 ac( 65)= 0.18113498 ac( 66)= 0.15603314 ac( 67)= 0.14236373 ac( 68)= 0.12819293 ac( 69)= 0.10743007 ac( 70)= 0.09416241 ac( 71)= 0.07652363 ac( 72)= 0.05719426 ac( 73)= 0.04319058 ac( 74)= 0.02332124 ac( 75)= 0.00492073 ac( 76)= -0.00940756 ac( 77)= -0.04927159 ac( 78)= -0.13574918。

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、画出解题的流程图
选定相关的参数
选择合适的单元来模拟混凝土、钢筋、垫板
定义混凝土单元和钢筋单元及垫板的材料参数
建立几何模型
进行网格尺寸指定并赋予材料属性
转换几何模型为有限元模型
施加荷载和边界条件
选择求解方法和屈服准则进行求解
进行后处理得到图形和数据。