ansys建模实例-混凝土箱梁(命令流)
最新ANSYS土木工程经典实例命令流大全
1.ANSYS SOLID65环向布置钢筋的例子 (3)2.混凝土非线性计算实例(1)- MISO单压 (5)3.混凝土非线性计算实例(2)-MISO约束压 (8)4.混凝土非线性计算实例(3)- KINH滞回 (11)5.混凝土非线性计算实例(4)- KINH压-拉裂 (14)6.混凝土非线性计算实例(5) (17)7.混凝土非线性计算实例(6) (20)8.混凝土非线性计算实例(7)-MISO滞回 (23)9.混凝土非线性计算实例(8) (27)10.混凝土非线性计算实例(9)-梁平面应力 (29)11.四层弹簧-质点模型的地震分析 (33)12.悬臂梁地震分析 (80)13.用beam 54单元描述变截面梁的例子 (103)14.变截面梁实例 (106)15.拱桥浇筑过程分析-单元生死应用实例 (107)16.简支梁实体与预应力钢筋分析实例 (109)17. 简单的二维焊接分析-单元生死实例 (111)18.隧道开挖(三维)的命令流 (118)19.岩土接触分析实例 (136)20.钢筋混凝土管的动力响应特性分析实例 (143)21.隧道模拟开挖命令流(入门) (150)22.螺栓连接的模拟实现问题 (154)23.道路的基层、垫层模量与应力之间的关系 (163)23.滞回分析 (185)24.模拟某楼层浇注 (187)25.在面上施加移动的面力 (190)27.在任意面施加任意方向任意变化的压力 (194)28.预紧分析 (195)29.几何非线性+塑性+接触+蠕变 (197)30.埋设在地下的排水管道 (202)32.幕墙企业玻璃简化计算 (206)33.等截面杆单元生死应用实例 (222)34.梁板建模联系 (224)36.简单的例子-如何对结构的振动控制分析 (227)37.模态分析结果的输出实例 (229)38.火车过桥动态加载实例(部分) (231)39.悬索结构的找形和计算的例题 (248)40.陶瓷杆撞击铝板的例子 (253)41.求反作用力的APDL命令法 (256)42.LS-DYNA实例(部分) (257)43.路面分层填筑对路基的影响 (258)44.一个例子(含地震影响,求振兴与频率) (262)45.接触面上的压力总和 (266)46.施加位置函数荷载 (270)47.非线性分析考虑刚度退化 (271)48.一个圆形水池的静力分析 (272)49.ANSYS中混凝土模式预应力模拟的算例 (273)50.悬臂梁受重力作用发生大变形求其固有频率 (276)51.循环对称结构模态分析 (278)52.三角平台受谐波载荷作用的结构响应 (280)53.三角平台受一地震谱激励的应力分布和支反力 (282)54.三角平台受时程载荷作用的应力分布和变形过程 (284)55.经典层合板理论 (286)56.定易圆轨迹的例子 (293)57.模拟门式刚架施工-单元生死 (293)58.钢筋混凝土整体式模型例子 (297)59.在荷载步之间改变材料属性例子 (299)60.含预应力的特征值屈曲计算 (301)61.振型叠加计算及工况组合例子 (304)62.柱子稳定分析算(预应力,特征值屈曲,初始缺陷) (307)63.module MConcrete !混凝土模板 (312)64.混凝土开裂实例 (327)65.螺栓网格划分 (328)66.自由液面的土石坝平面渗流分析 (329)67.导出刚度矩阵 (333)68.某混凝土拱坝工程施工期及运行期温度场仿真分析 (334)69.移动温度荷载计算 (340)70.SHSD用于壳-实体装配实例An (343)71.ansys显示-隐式-回弹分析实例 (347)72.工况组合的经典例子 (361)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,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,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,0 NSEL,S,LOC,Y,50 CP,1,UY,ALLD,ALL,UX,0D,ALL,UZ,0 NSEL,ALL FINISH/VIEW,1,1,1,1 /REPLOT/SOLU OUTRES,ALL,1 TIME,20 NSUBST,20,0,20 D,1,UY,-.04 LSWRITE,1 TIME,30 NSUBST,20,0,20 D,1,UY,-.06 LSWRITE,2 TIME,40 NSUBST,100D,1,UY,-.1 LSWRITE,3 LSSOLVE,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,5FINISH4.混凝土非线性计算实例(3)- KINH滞回!KINH滞回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,10 TBPT,, 150E-6, 4.5 TBPT,, 600E-6,14.8 TBPT,,1000E-6,20.25 TBPT,,1300E-6,22.3 TBPT,,1480E-6,22.8 TBPT,,1620E-6,22.8 TBPT,,1800E-6,22.3 TBPT,,2000E-6,21.4 TBPT,,3500E-6,12.8 TBPT,,5000E-6, 9.0 BLOCK,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,50 CP,1,UY,ALL NSEL,ALL FINISH/VIEW,1,1,1,1 /REPLOT/SOLU OUTRES,ALL,ALL TIME,20 NSUBST,20,0,20 D,1,UY,-.04 LSWRITE,1 TIME,35 NSUBST,15,0,15 D,1,UY,-.01 LSWRITE,2 TIME,90 NSUBST,55,0,55 D,1,UY,-.12 LSWRITE,3 TIME,125 NSUBST,35,0,35 D,1,UY,-.05 LSWRITE,4TIME,175NSUBST,50,0,50D,1,UY,-.15LSWRITE,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,5FINISH5.混凝土非线性计算实例(4)- KINH压-拉裂!KINH压——拉裂压到峰值,泄载、反向加载到拉裂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,10 TBPT,, 150E-6, 4.5 TBPT,, 600E-6,14.8 TBPT,,1000E-6,20.25 TBPT,,1300E-6,22.3 TBPT,,1480E-6,22.8 TBPT,,1620E-6,22.8 TBPT,,1800E-6,22.3 TBPT,,2000E-6,21.4 TBPT,,3500E-6,12.8 TBPT,,5000E-6, 9.0 BLOCK,0,50,0,50,0,50 ESIZE,10VMESH,ALLMSHAPE,0,3D NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,UY,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 NSEL,S,LOC,Y,50 CP,1,UY,ALL NSEL,ALL FINISH/VIEW,1,1,1,1 /REPLOT/SOLU OUTRES,ALL,ALL TIME,20 NSUBST,20,0,20 D,1,UY,-.04 LSWRITE,1 TIME,45 NSUBST,25,0,25 D,1,UY,.01 LSWRITE,2 LSSOLVE,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.9 BLOCK,0,50,0,50,0,50ESIZE,10 VMESH,ALL MSHAPE,0,3D NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,UY,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 NSEL,S,LOC,Y,50 CP,1,UY,ALL NSEL,ALL FINISH/VIEW,1,1,1,1 /REPLOT/SOLU OUTRES,ALL,ALL TIME,50 NSUBST,50,0,50 D,1,UY,-.25 SOLVEFINISH/POST26 NSOL,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,3D NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,UY,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 NSEL,S,LOC,Y,50 CP,1,UY,ALLD,ALL,UX,0D,ALL,UZ,0 NSEL,ALL FINISH/VIEW,1,1,1,1 /REPLOT/SOLU OUTRES,ALL,1 TIME,20 NSUBST,20,0,20 D,1,UY,-.04 LSWRITE,1 TIME,30 NSUBST,20,0,20 D,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,1 TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1 TB,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.3TBPT,,4600E-6,9.9 BLOCK,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,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,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,1 TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1 TB,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.8TBPT,,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,14 E,47,63E,48,64EGEN,10,88,589,590 /NUMBER,1/PNUM,TYPE,1/REPLOTNSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,0 NSEL,R,LOC,Z,740 D,ALL,UY,0NSEL,ALLFINISH/SOLUANTYPE,STATIC NSEL,S,LOC,Y,200 NSEL,R,LOC,Z,200 D,ALL,UY,-2 TIME,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.02 ET,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,4LSEL,,,,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,0FINISH/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,DISPFINISH/POST1/INP,,MCOMNSORT,U,SUM*GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形!开始地震计算FINISH/CONFIG,NRES,20000/SOLUANTYPE,TRANSTRNOPT,FULLALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926 BETAD,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.0042 -0.0066 -0.0131 -0.019 -0.0196 -0.0066 0.003 0.0141 -0.0049 -0.0128 -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.02360.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.00590.0308 0.0499 0.071 0.0995 0.1219 0.1529 0.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.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.0753 -0.0173 0.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.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.113 0.1363 0.02190.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.01090.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.0055 -0.0147 0.0143 0.0206 0.0499 0.06450.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.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.0469 -0.125 -0.2111 -0.1617 -0.1692 -0.1306 -0.1111 -0.0773 -0.051 -0.0544-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.13520.1861 0.1281 0.064 0.0204 0.0314 0.0373 0.0496 0.0235 -0.0084 -0.0168 -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。
Ansys绘制梁的弯矩图和剪力图(图文教程+APDL命令流)
Y ZX
-.625E-09 3482.91
6965.81
10448.7
13931.6
17414.5
20897.4
24380.3
27863.2
31346.2
Step17 弯矩图
9
1 LINE STRESS STEP=1 SUB =1 TIME=1 SZ_I SZ_J MIN =-15673.1 ELEM=18 MAX =15076.9 ELEM=61
6
Step11 施加集中载荷 Step12 施加梯度载荷
7
Step13 力学模型
Step14 计算
Step15 计算完成
8
1 NODAL SOLUTION
STEP=1
SUB =1
TIME=1
UY
(AVG)
RSYS=0
DMX =.073299
SMN =-.073299
Y MZX X
MN
-.073299
*do,i,31,45,1 sfbeam,i,2,pres,1e3/15*(i-31),1e3/15*(i-31+1) *enddo allsel,all !!!!计算 /solu
11
solve finish /post1 ETABLE,Mz_I,SMISC, 3 ETABLE,Mz_J,SMISC, 16 PLLS, Mz_I, Mz_J ETABLE,Sz_I,SMISC, 6 ETABLE,Sz_J,SMISC, 19 PLLS, Sz_I, Sz_J
27863.2
31346.2
图 1 结果图
1
图文教程
Step 1 选择单元类型 Step 2 设定材料属性
2
FRP混凝土用ANSYS模拟的命令流文件
FRP混凝土用ANSYS模拟的命令流文件Combin39单元学习小结(本文来源于清溪长河) 很有用,转过来慢慢学了(2011-07-25 10:01:48)转载标签: ansyscombine39单元杂谈分类: 专业学习这两天在帮助一位同学解决论文中的问题,其间涉及到组合结构不同材料的构件的连接。
原本,她想做的是二者分离的动态过程。
要做到这一点,貌似要用ls-dyna,而这我知之甚匮,加上时间非常有限,这个方案被我首先排除。
至于接触分析,在参考了一些资料后,感觉即便用点-点接触来进行模拟,也会造成计算机资源的浪费。
而且此分析经过了MTS的实验,显示连接很紧密,甚至直接耦合的话也相去不远。
经过慎重考虑和探讨,也放弃了做接触分析的想法。
在陆新征的那篇引用数高达230的文章的启发下,终于确定了分析方案:用Combin39弹簧单元来模拟组合结构之间的相对黏结滑移。
其实,用combin39模拟钢筋-混凝土之间的滑移,已经是比较成熟的选择了,可我以前并不十分清楚这一点~囧……看到网上有不少和我一样的combin新手,在问一些低级问题(正是由于他们,所以我只需要Google就OK了~),我就把最近学到的一些知识总结如下,欢迎讨论或者批评指正。
bin39是一种包含两个节点的一维的具有非线性性能的单向受力单元,在1-D, 2-D, 或3-D应用中具有纵向的和扭转的能力,单元的每个节点都具有X、Y、Z三个自由度。
2.该单元具有大变形的能力,此时需要打开NLGEOM选项时,并且使KEYOPT(4) = 1 或者3。
bin39需要通过一个荷载,位移(F—D)曲线来定义弹簧的受力性质,至于这个曲线,可以通过查询相关文献获得,比如[杨勇. 型钢混凝土粘结滑移基本理论及应用研究. 西安建筑科技大学,2003.] 无需定义材料属性。
4.和Combin40不同的是,Combin39的刚度不是用户手动输入的,而是程序自动根据F-D曲线求的。
ANSYS分析钢筋混凝土结构技巧及实例详解
0 前言利用ANSYS分析钢筋混凝土结构时,其有限元模型主要有分离式和整体式两种模型。
这里结合钢筋混凝土材料的工作特性,从模型建立到非线性计算再到结果分析的全过程讲述了利用ANSYS进行钢筋混凝土结构分析的方法与技巧,并以钢筋混凝土简支梁为例,采用分离式有限元模型,说明其具体应用。
1 单元选取与材料性质1. 1 混凝土单元ANSYS中提供了上百种计算单元类型,其中Solid65单元是专门用于模拟混凝土材料的三维实体单元。
该单元是八节点六面体单元,每个节点具有三个方向的自由度( UX , UY , UZ) 。
在普通八节点线弹性单元Solid45 的基础上,该单元增加了针对于混凝土的材性参数和组合式钢筋模型,可以综合考虑包括塑性和徐变引起的材料非线性、大位移引起的几何非线性、混凝土开裂和压碎引起的非线性等多种混凝土的材料特性。
使用Solid65 单元时,一般需要为其提供如下数据:1)、实常数(Real Constants) :定义弥散在混凝土中的最多三种钢筋的材料属性,配筋率和配筋角度。
对于墙板等配筋较密集且均匀的构件,一般使用这种整体式钢筋混凝土模型。
如果采用分离式配筋,那么此处则不需要填写钢筋实常数。
2)、材料模型(Material Model) :在输入钢筋和混凝土的非线性材料属性之前,首先必须定义钢筋和混凝土材料在线弹性阶段分析所需的基本材料信息,如:弹性模量,泊松比和密度。
3)、数据表(Data Table) :利用数据表进一步定义钢筋和混凝土的本构关系。
对于钢筋材料,一般只需要给定一个应力应变关系的数据表就可以了,譬如双折线等强硬化(bilinear isotropic hardening)或随动硬化模型( kinematic hardening plasticity)等。
而对于混凝土模型,除需要定义混凝土的本构关系外,还需要定义混凝土材料的破坏准则。
在ANSYS中,常用于定义混凝土本构关系的模型有:1)多线性等效强化模型(Multilinear isotropic hardening plas2ticity ,MISO模型),MISO模型可包括20条不同温度曲线,每条曲线可以有最多100个不同的应力-应变点;2)多线性随动强化模型(Multilinear kinematic hardening plas2ticity ,MKIN 模型),MKIN 模型最多允许5个应力-应变数据点;3)Drucker2Prager plasticity(DP)模型。
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仿真模拟软件建立三维混凝土试件实体裂纹扩展的模拟
MSHKEY,0
!*
CM,_Y,VOLU
VSEL, , , , 3
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!施加下端约束
FLST,2,1,5,ORDE,1
FITEM,2,18
!通水管道施压
FLST,2,2,5,ORDE,2
FITEM,2,12
FITEM,2,-13
/GO
!*
SFA,P51X,1,PRES,-2000000
!裂纹面施压
FLST,2,4,5,ORDE,2
FITEM,2,14
FITEM,2,-17
/GO
!*
SFA,P51X,1,PRES,2000000
Fini !退出前处理器
wpro,,90.000000, !旋转工作平面
CSWPLA,100,1,1,1, !在工作平面位置建立局部坐标100,类型为柱坐标
FLST,3,1,6,ORDE,1
FITEM,3,1
VGEN, ,P51X, , , ,45, , , ,1 !旋转长方体
wpro,,-90.000000, !旋转回原工作平面
!(4)定义材料参数
MP,EX,1,1.668E10 !弹性模量
MP,PRXY,1,0.3 !泊松比
!(5)建立剖面几何模型
BLOCK,-0.015,0.015,-0.025,0.025,-0.0005,0.0005, !建立一个长方体
WPSTYLE,,,,,,,,1
钢筋混凝土梁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命令流
finish$/clear$/config,fsplit,750$/prep7!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义材料!!!!!!!!!下面定义材料 1,Q235A型钢钢材!!!!!!!!!mp,ex,1,2.1e11$mp,prxy,1,0.3$mp,dens,1,7850tb,bkin,1tbdata,,268.85e6,0!!下面定义材料 9,Q235A型钢钢材!!!!!!!!!mp,ex,3,2.1e11$mp,prxy,3,0.3$mp,dens,3,7850tb,bkin,3$tbdata,,235e6,0!!!下面定义材料 2,混凝土C40!!!!!!!!!fck=26*0.85*0.8e6$ftk=2.3e6mp,ex,2,fck*0.44/0.0005$mp,prxy,2,0.2$mp,dens,2,2600tb,concr,2$tbdata,,0.125,0.7,FTK,-1tb,mkin,2$tbtemp,,strain ! 混凝土屈服准则Tbdata,1,0.0005,0.001,0.0015,0.002,0.0038$tbtemp,0 $Tbdata,1,0.44*fck,0.75*fck,0.94*fck,fck,0.85*fck !混凝土应力应变曲线!!!下面定义材料 4,纵钢筋直径20!!!!!!!!! mp,ex,4,2.1e11$mp,prxy,4,0.3$mp,dens,4,7850 tb,bkin,4tbdata,,365e6,0!!!下面定义材料7,箍筋直径8!!!!!!!!!mp,ex,5,2.1e11$mp,prxy,5,0.3$mp,dens,5,7850 tb,bkin,5tbdata,,243e6,0!!!!定义实常数!!!!!!Pi=acos(-1)$ r,1$ r,2$r,3,0.25*pi*0.020*0.020r,4,0.25*pi*0.008*0.008R,5!!!!定义单元!!!!!!!et,1,solid45$et,2,solid65$et,3,LINK8keyopt,2,7,1$keyopt,2,1,1$keyopt,2,3,2!!!!建模型!!!!!!!!柱子wpstyl,defa$Block,-0.2,0.2,-1.5,1.5,0,0.2cm,zhuzi,volu!工字钢vsel,none$wpstyl,defa$$block,-2,2,-0.125,0.125,0,0.003$bloc k,-2,2,0.125-0.015,0.125,0,0.078$$block,-2,2,-0.125,-0.125+0.015,0,0.078$$wpstyl,defa$wpoff,-0.2$$block,0,0.010,-0.125,0.125,0,0.078 $wpstyl,defa$wpoff,0.2$$block,-0.010,0,-0.125,0.125,0,0.078 $wpstyl,defa$wpoff,,0.125-0.015$block,-0.075,0.075,0,1.5-0.125+0.015,0,0.003$$block,-0.075,-0.075+0.008,0,1.39,0,0.078$$block,0.075-0.008,0.075,0,1.39,0,0.078$$wpstyl,defa$wpoff,,-0.125+0.015$block,-0.075,0.075,-1.39,0,0,0.003$$block,-0.075,-0.075+0.008,-1.39,0,0,0.078$$block,0.075-0.008,0.075,-1.39,0,0,0.078$$wpstyl,defa$wpoff,,-0.125+0.015$block,-0.2,0.2,-0.06-0.015,0,0,0.003$$wpstyl,defa$wpoff,,0.125-0.015$block,-0.2,0.2,0,0.060+0.015,0,0.003$$wpstyl,defa$wpoff,-0.2$$block,0,0.010,0.125,0.125+0.06,0,0 .2$wpstyl,defa$wpoff,0.2$$block,-0.010,0,0.125,0.125+0.06,0,0 .2$wpstyl,defa$wpoff,-0.2$$block,0,0.010,-0.125-0.06,-0.125,0 ,0.2$wpstyl,defa$wpoff,0.2$$block,-0.010,0,-0.125-0.06,-0.125,0 ,0.2$wpstyl,defa$wpoff,,-0.125$wpoff,,,0.2$block,-0.2,0.2,-0.06,0,-0.010,0$$wpstyl,defa$wpoff,,0.125$wpoff,,,0.2$block,-0.2,0.2,0,0.06,-0.010,0$$wpstyl,defa$wpoff,-0.2$$block,0,0.010,-0.11,0.11,0,0.2$wps tyl,defa$wpoff,0.2$$block,-0.010,0,-0.11,0.11,0,0.2 $wpstyl,defa$wpoff,,,0.2$block,-0.2,0.2,-0.11,0.11,-0.010,0$Nummrg,kp$vptn,all$cm,gzg,volu$numcmp,all$allsel,all$vptn,all$$allsel,allVsel,u,,,gzgCm,zhuzi,volu!切分gzgCmsel,s,gzg,volu$Nummrg,kp$!画纵筋$wpstyl,defa$wpoff,0.2$wpoff,,,0.160$vsbw,all$ $wpstyl,defa$wpoff,0.160$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,-0.16$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,0.110$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,-0.11$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,0.2$wpoff,,,0.078$vsbw,all$$numcmp,all!画箍筋$wpstyl,defa$wprota,,90 $wpoff,,,0.04$vsbw,all $wpoff,,,0.070$vsbw,all$wpoff,,,0.015$vsbw,all$wpoff,,,0.030$vsbw,all$wpoff,,,0.045$vsbw,all$wpoff,,,0.05$vsbw,all$*do,i,1,12,1$wpoff,,,0.10$vsbw,all$*enddo$wpstyl,defa$wprota,,-90 $$wpoff,,,0.04$vsbw,all $wpoff,,,0.070$vsbw,all$wpoff,,,0.015$vsbw,all$wpoff,,,0.030$vsbw,all$wpoff,,,0.045$vsbw,all$wpoff,,,0.05$vsbw,all$*do,i,1,12,1$wpoff,,,0.10$vsbw,all$*enddo$$wpstyl,defa$wprota,,,90$wpoff,,,0.075-0.008$vsbw,all $wpoff,,,0.008$vsbw,all$wpstyl,defa$wprota,,,-90$wpoff,,,0.075-0.008$vsbw,all $wpoff,,,0.008$vsbw,all$wpstyl,defa$wpoff,,,0.003$vsbw,all$wpstyl,defa$wpoff,,,0.078$vsbw,all$wpstyl,defa$wprota,,,90$vsbw,all$wpstyl,defa$wprota,,90 $wpoff,,,0.125+0.06$vsbw,all $wpstyl,defa$wprota,,-90 $wpoff,,,0.125+0.06$vsbw,all$numcmp,allcm,gzg,volu!切分柱cmsel,s,zhuzi,volu$!画纵筋$wpstyl,defa$wpoff,0.2$wpoff,,,0.160$vsbw,all$ $wpstyl,defa$wpoff,0.160$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,-0.16$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,0.110$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,-0.11$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,,,0.078$vsbw,all$$numcmp,all!画箍筋$wpstyl,defa$wprota,,90 $wpoff,,,0.04$vsbw,all $wpoff,,,0.070$vsbw,all$wpoff,,,0.015$vsbw,all$wpoff,,,0.030$vsbw,all$wpoff,,,0.045$vsbw,all$wpoff,,,0.05$vsbw,all$*do,i,1,12,1$wpoff,,,0.10$vsbw,all$*enddo$wpstyl,defa$wprota,,-90 $$wpoff,,,0.04$vsbw,all $wpoff,,,0.070$vsbw,all$wpoff,,,0.015$vsbw,all$wpoff,,,0.030$vsbw,all$wpoff,,,0.045$vsbw,all$wpoff,,,0.05$vsbw,all$*do,i,1,12,1$wpoff,,,0.10$vsbw,all$*enddo$$wpstyl,defa$wprota,,,90$wpoff,,,0.075-0.008$vsbw,all$wpoff,,,0.008$vsbw,all$wpstyl,defa$wprota,,,-90$wpoff,,,0.075-0.008$vsbw,all $wpoff,,,0.008$vsbw,all$wpstyl,defa$wpoff,,,0.003$vsbw,all$wpstyl,defa$wprota,,,90$vsbw,all$wpstyl,defa$wpoff,0.2-0.01$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,-0.2+0.01$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,0.2$wpoff,,,0.2-0.01$vsbw,all$ $wpstyl,defa$wprota,,90 $wpoff,,,0.125+0.06$vsbw,all $wpstyl,defa$wprota,,-90 $wpoff,,,0.125+0.06$vsbw,all$numcmp,allcm,zhuzi,volu!!!!!!!!!!赋钢筋单元属性并划分网格!!!!!!!!!elemsiz=50e-3lsel,s,loc,x,-0.16$lsel,r,loc,z,0.078$cm,waizj1,line$latt,4 ,3,3,0$lesize,all,elemsizlsel,s,loc,x,-0.16$lsel,r,loc,z,0.16$cm,waizj2,line$latt,4,3,3,0$lesize,all,elemsizlsel,s,loc,x,-0.11$lsel,r,loc,z,0.16$cm,waizj3,line $latt,4,3,3,0$lesize,all,elemsizlsel,s,loc,x,0.16$lsel,r,loc,z,0.078$cm,ywaizj1,line$latt,4 ,3,3,0$lesize,all,elemsizlsel,s,loc,x,0.16$lsel,r,loc,z,0.16$cm,ywaizj2,line $latt,4,3,3,0$lesize,all,elemsizlsel,s,loc,x,0.11$lsel,r,loc,z,0.16$cm,ywaizj3,line $latt,4,3,3,0$lesize,all,elemsiz$cmsel,s,waizj1,line$$cmsel,a,waizj2,line$$cmsel,a,waizj3,l ine$$cmsel,a,ywaizj1,line$$cmsel,a,ywaizj2,line$$cmsel,a,ywaizj 3,line$$cm,zj,line!节点外域箍筋lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161 $lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.295,0.295 $lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,1.5$lsel,u,loc,y,-1.5$lsel,r,loc,z,0,0.16$lsel,r,loc,x,-0.16$cm,sgj1,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161 $lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.295,0.295$lsel,u,loc,y,1.5$lsel,u,loc,y,-1.5$lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,r,loc,z,0,0.16$lsel,r,loc,x,0.16$cm,sgj2,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161 $lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.295,0.295 $lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,1.5$lsel,u,loc,y,-1.5$lsel,r,loc,z,0.16$cm,sgj3,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161$lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015$lsel,u,loc,y,-0.295,0.295$lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,1.5$lsel,u,loc,y,-1.5$Lsel,r,loc,z,0.078$cm,sgj4,line$$cmsel,s,sgj1,line$$cmsel,a,sgj2,line$ cmsel,a,sgj3,line$$cmsel,a,sgj4,line$cm,guj1,line$!节点区域箍筋!节点区域上下四肢箍箍筋lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,-0.095,0.095 $lsel,u,loc,y,0.251,1.5$lsel,u,loc,y,-1.5,-0.251$lsel,u,loc,y, 0.125$lsel,u,loc,y,-0.125$Lsel,r,loc,z,0.16$cm,ssgj1,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161 $lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,0.251,1.5$lsel,u,loc,y,-1.5,-0.251$lsel,u,loc,y, 0.125$lsel,u,loc,y,-0.125$Lsel,r,loc,z,0.078$cm,ssgj2,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161 $lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015$lsel,u,loc,y,-0.095,0.095 $lsel,u,loc,y,0.251,1.5$lsel,u,loc,y,-1.5,-0.251$lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y, 0.125$lsel,u,loc,y,-0.125$Lsel,r,loc,x,-0.16$Lsel,r,loc,z,0,0.16$cm,ssgj3,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161$lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015$lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,-0.095,0.095$lsel,u,loc,y,0.251,1.5$lsel,u,loc,y,-1.5,-0.251$lsel,u,loc,y, 0.125$lsel,u,loc,y,-0.125$Lsel,r,loc,x,0.16$Lsel,r,loc,z,0,0.16$cm,ssgj4,line$$cmsel,s,ssgj1,line$$cmsel,a,ssgj2,line$cmsel,a,ssgj3,line$cmsel,a,ssgj4,line$$cm,guj3,line$$cmsel,s,guj1,line$$cmsel,a,guj3,line$$cm,guj,line$latt,5,4 ,3,0$lesize,all,elemsiz$cmsel,s,guj,line$cmsel,a,zj,line$cm,gj,lineCmsel,s,gzg,volu$VSEL,r,LOC,Y,-0.125+0.015,0.125-0.015 $VSEL,r,LOC,z,0,0.003$cm,fb,VOLUCmsel,s,gzg,volu$VSEL,u,LOC,Y,-0.125,0.125$cm,GZhu,VOLU Cmsel,s,gzg,volu$$cmSEL,U,fb,VOLU$VSEL,r,LOC,Y,-0.125+0.015 ,0.125-0.015$cm,mianban,VOLUCmsel,s,gzg,volu$$cmSEL,U,fb,VOLU$$cmSEL,U,gzhu,VOLU$$cmSEL ,U,mianban,VOLU$$cm,yy,VOLU!!!!!!!!!!!!!!!!!!!!!!!! Mesh钢筋!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!cmsel,s,gj,line$lmesh,all $lsel,all$/eshape,1$eplot !钢筋划分完毕!!!!!!!!!!!!!!!!!!!! mesh 钢筋混凝土柱!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Cmsel,s,zhuzi,volu $esize,50e-3$vatt,2,1,2,0$MSHAPE,0,3D$MSHKEY,1$VMESH,allallsel,all!!!!!!!!!!!!!!!!!! Mesh 翼缘!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Cmsel,s,yy,volu$$LSEL,s,length,,1.8LESIZE,all, , ,12vatt,1,1,1,0$esize,50e-3$$MSHAPE,0,3D$MSHKEY,1$VMESH,all$allsel,all!!!!!!!!!!!!!!!!!! Mesh fb!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Cmsel,s,fb,volu$vatt,1,1,1,0$esize,50e-3$$MSHAPE,0,3D$MSHKEY,1$VMESH,all$allsel,all!!!!!!!!!!!!!!!!!! Mesh mianban gzhu!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Cmsel,s,mianban,volu$Cmsel,a,gzhu,volu$vatt,3,5,1,0$esize,50e-3$$MSHAPE,0,3D$MSHKEY,1 $VMESH,all$allsel,all$nummrg,node$numcmp,all!定义远距离加载点$ET,4,TARGE170ET,5,CONTA174R,6keyopt,4,2,1$keyopt,5,2,2$keyopt,5,4,2$keyopt,5,12,5 asel,s,loc,x,2$nsla,s,1$real,6$type,5$esurfnmax=ndinqr(0,14)+1$n,nmax,2,,type,4$tshap,pilo$e,nmaxR,7keyopt,4,2,1$keyopt,5,2,2$keyopt,5,4,2$keyopt,5,12,5 asel,s,loc,x,-2$nsla,s,1$real,7$type,5$esurfnmax=ndinqr(0,14)+1$n,nmax,-2,,type,4$tshap,pilo$e,nmaxR,8keyopt,4,2,1$keyopt,5,2,2$keyopt,5,4,2$keyopt,5,12,5 asel,s,loc,y,1.5$nsla,s,1$real,8$type,5$esurfnmax=ndinqr(0,14)+1$n,nmax,0,1.5,type,4$tshap,pilo$e,nmaxAllsel,all!施加约束和荷载Allsel,allAsel,s,loc,z,0Da,all,symm$D,8665,uy,,,,,uz$D,8666,uy,,,,,uz$Asel,s,loc,y,-1.5Da,all,allAllsel,all/solu$ANTYPE,0$OUTRES,ALL,ALL$$pred,off$neqit,30$lnsrch,on$nlgenom,off$Cnvtol,f,,0.05$kbc,0$eqslv,sparse,1.0E-5,2 $p=0.01$time,1$nsubst,50,,10$Asel,s,loc,y,1.5f,8667,fy,-582000$$allsel,all$OUTRES,ALL,ALL $ALLSEL,ALL $LSWRITE,1TIME,2$D,8667,Ux,19*p$$kbc,0$ NSUBST,100,1000000,1$AUTOTS,o n$OUTRES,ALL,ALL $ALLSEL,ALL$LSWRITE,2$Lssolve,1,2$save$finish/post1SET,LAST$SET,PREVIOUS*GET,RTIME,ACTIVE,0,SET,TIME !获取最终收敛时间Pldisp$plnsol,s,x !显示变形图与X方向应力分布图Plnsol,epel,x !X弹性应变分布图Plnsol,eppl,x !x塑性应变分布图Plnsol,epto,x !X方向总应变分布图Plnsol,nl,srat !显示塑性区域/device,vector,off$plcrack,Cmsel,s,gzgEslv,s/POST26nsel,s,loc,x,0.975NUMVAR,200$NSOL,2,8667,U,x$nsel,s,loc,y,-1.55$RFORCE,3,8667,F,xcmsel,s,gj,line$Esll,s$/output,Z1,csv$prod,4,2,,,,,,1000$Prod,5,3,,,,,,1/1000PRVAR,4,5!!! 钢梁应力的显示/post1$Cmsel,s,gzg$cmsel,a,FUBAN$vsel,r,loc,z,0,0.0045$vsel,r,l oc,y,-0.125+0.006,0.125-0.006$vsel,r,loc,x,-0.175,0.175$eslv,s$nsle,splnsol,s,eqvPlnsol,epel,x !X弹性应变分布图Plnsol,eppl,x !x塑性应变分布图Plesol,epto,1 !X方向总应变分布图/post1$Cmsel,s,gzg$$eslv,s$nsle,s$plnsol,s,eqv/post1$Cmsel,s,zhuzi$eslv,s$nsle,splnsol,s,eqv$Cmsel,s,zhuzi$cmsel,a,FUBAN$vsel,r,loc,z,0,0.0045$vsel,r ,loc,y,-0.125+0.006,0.125-0.006$vsel,r,loc,x,-0.175,0.175$e slv,s$nsle,splesol,s,eqv/post1$Cmsel,s,zhuzi$eslv,s$nsle,splnsol,s,eqv/post1SET,LAST$SET,PREVIOUS*GET,RTIME,ACTIVE,0,SET,TIME !获取最终收敛时间Pldisp$plnsol,s,x !显示变形图与X方向应力分布图Plnsol,epel,x !X弹性应变分布图Plnsol,eppl,x !x塑性应变分布图Plnsol,epto,eqv !X方向总应变分布图Plnsol,nl,srat !显示塑性区域/post1Cmsel,s,gj$esll,sETABLE,rebar,ls,1PLETAB,rebar/post1$wpstyl,defa$$path,23,2,30,120$ppath,1,,0.185,0.125,0$ppath ,2,,0.975,0.125,0pdef,myseqv,s,eqv$pdef,mysy,s,1plpath,myseqv,mysyPLPAGM,myseqv,1,blankPLPAGM,myseqv,1,nodeprpath,myseqv,,xg,yg,zg/device,vector,on$plcrack,Cmsel,s,gzgEslv,snsel,s,loc,x,0.975NUMVAR,200$NSOL,2,9989,U,Y$nsel,s,loc,y,-1.55$RFORCE,3,9989,F,Y/output,Z1,csv$prod,4,2,,,,,,-1000$Prod,5,3,,,,,,-1/1000 PRVAR,4,5/POST26nsel,s,loc,x,0.975NUMVAR,200$NSOL,6,10678,U,Y$nsel,s,loc,y,-1.55$RFORCE,7,10678,F,Y/output,Z2,csv$prod,8,6,,,,,,1000$Prod,9,7,,,,,,1/1000!后处理阶段/post1SET,LAST$SET,PREVIOUS*GET,RTIME,ACTIVE,0,SET,TIME !获取最终收敛时间Pldisp$plnsol,s,x !显示变形图与X方向应力分布图Plnsol,epel,x !X弹性应变分布图Plnsol,eppl,x !x塑性应变分布图Plnsol,epto,eqv !X方向总应变分布图Plnsol,nl,srat !显示塑性区域esel,s,type,,2cmsel,s,gj,line$esll,sETABLE,rebar,ls,1PLETAB,rebar!混凝土裂缝的显示!外混凝土$Cmsel,s,zhuzi,$vsel,r,loc,z,0.0625+0.001,0.175$vsel,r,lo c,y,-0.125+0.006,0.125-0.006$eslv,s$$/dscale,,off$/device,vector,on$plcrack,1,1 !第一条裂纹!内混凝土$Cmsel,s,zhuzi$vsel,r,loc,z,0,0.0625$vsel,r,loc,y,-0.125+ 0.006,0.125-0.006$eslv,s$/dscale,,off$$/device,vector,on$pl crack/post1!混凝土应力的显示$/gformat,defa !恢复缺省格式$Cmsel,s,zhuzi$eslv,s$plnsol,s,1 !主应力$Cmsel,s,zhuzi$eslv,s$plnsol,s,3 !最大压应力$Cmsel,s,zhuzi$eslv,s$plnsol,s,eqv !等效应力$Cmsel,s,gzg$Cmsel,a,fb$eslv,s$plnsol,s,eqv/post1!内混凝土矢量图$Cmsel,s,zhuzi,$vsel,r,loc,z,0,0.0625$vsel,r,loc,y,-0.125 +0.006,0.125-0.006$eslv,s$$/device,vector,off$/dscale,,off$$plnsol,s,eqv $Plvect,u$$Plvect,s$/vscale,,1$plvect,s,,,,vect !调整箭头大小$/gline,,-1$plvect,s,,,,vect ! 取消单元边界$/dscale,,off$$plnsol,s,1 !关闭结构变形$/dscale,defas$/gline,1 !恢复变形比例缺省,显示单元边界$/Device,vector,1 !设置矢量显示模式$/clabel,,5 !设置等高线上的文字间隔$/contour,,18,-16,,500 !设置等高线数目$Plnsol,s,x!外混凝土矢量图$Cmsel,s,zhuzi,$vsel,r,loc,z,0.0625+0.001,0.175$vsel,r,lo c,y,-0.125+0.006,0.125-0.006$eslv,s$/device,vector,off$/dscale,,off$$plesol,s,3$Plvect,u$$Plvect,s$/vscale,,1$plvect,s,,,,vect !调整箭头大小$/gline,,-1$plvect,s,,,,vect ! 取消单元边界$/dscale,,off$$plnsol,s,1 !关闭结构变形$/dscale,defas$/gline,1 !恢复变形比例缺省,显示单元边界$/Device,vector,1 !设置矢量显示模式$/clabel,,5 !设置等高线上的文字间隔$/contour,,18,-16,,500 !设置等高线数目$Plnsol,s,x!混凝土主压应变与等效应变曲线$Cmsel,s,zhuzi,$vsel,r,loc,z,0,0.175$vsel,r,loc,y,-0.125+ 0.006,0.125-0.006$vsel,r,loc,x,-0.135+0.03,0.175$eslv,s$nsle,s/POST26NUMVAR,200$$Esol,3,1080,26,s,eqv$esol,4,1080,26,epel,eqv$esol,5,1080,26,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa)$PRVAR,4,5/POST26NUMVAR,200$$Esol,3,1080,45,s,eqv$esol,4,1080,45,epel,eqv$esol,5,1080,45,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa) $PRVAR,4,5/POST26NUMVAR,200$$Esol,3,2706,46,s,eqv$esol,4,2706,46,epel,eqv$esol,5,2706,46,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa) $PRVAR,4,5/POST26NUMVAR,200$$Esol,3,3157,296,s,eqv$esol,4,3157,296,epel,eqv$esol,5,3157,296,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa) $PRVAR,4,5/POST26NUMVAR,200$$Esol,3,3157,291,s,eqv$esol,4,3157,291,epel,eqv$esol,5,3157,291,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa)$PRVAR,4,5/POST26NUMVAR,200$$Esol,3,3160,1169,s,eqv$esol,4,3160,1169,epel,eqv$esol,5,3160,1169,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa)$PRVAR,4,5! 剪应变与梁端荷载曲线$Cmsel,s,gzg$vsel,r,loc,z,0,0.0045$vsel,r,loc,y,-0.125+0. 006,0.125-0.006$eslv,s$nsle,s/POST26NUMVAR,200$$RFORCE,3,7150,F,Y$esol,4,6351,2661,epel,xy$esol,5,6351,2661,eppl,xy$add,4,4,5$prod,5,3,,,,,,-1/1000$PRVAR,4,5/post1$Cmsel,s,gzg$vsel,r,loc,z,0,0.0045$vsel,r,loc,y,-0.125+0. 006,0.125-0.006$$vsel,r,loc,x,-0.175,0.175$eslv,s$nsle,sSET,LAST$SET,PREVIOUS*GET,RTIME,ACTIVE,0,SET,TIME !获取最终收敛时间Pldisp$plnsol,s,x !显示变形图与X方向应力分布图Plnsol,epel,xy !Xy弹性应变分布图Plnsol,eppl,xy !xy塑性应变分布图Plnsol,epto,xy !Xy方向总应变分布图/post1Allsel,all!腹板剪力$Cmsel,s,gzg$Nslv,s,1$Nsel,r,loc,x,-0.175,0.175$Nsel,r,loc,z,-0.0045,0.0045 $Nsel,r,loc,y,0,0.05$Esln,s,1$Nsel,r,loc,y,0$fsum!内混凝土剪力Allsel,all$Cmsel,s,zhuzi$Nsel,r,loc,x,-0.175,0.175$Nsel,r,loc,z,-0.0625,0.0625$Nsel,r,loc,y,0,0.05$Esln,s,1$Nsel,r,loc,y,0$fsum/post1Allsel,all!外混凝土剪力$Cmsel,s,zhuzi$Nslv,s,1$Nsel,r,loc,x,-0.175,0.175$Nsel,u,loc,z,-0.0625+0.001,0.0625-0.001 $Nsel,r,loc,y,0,0.05$Esln,s,1$Nsel,r,loc,y,0$fsum!混凝土总剪力Allsel,all$Cmsel,s,zhuzi$Nsel,r,loc,x,-0.175,0.175 $Nsel,r,loc,y,0,0.05$Esln,s,1$Nsel,r,loc,y,0$fsum!总剪力/post1Allsel,all$Nsel,s,loc,x,-0.175,0.175 $Nsel,r,loc,y,0,0.05$Esln,s,1$Nsel,r,loc,y,0$fsum。
混凝土梁的ANSYS命令流
! 水***** * *****定义支座垫板的几何特性及材料特性****** * *
E , O I4 T 3 S LD 5 , M, 33O 1 PE, . 0 X , E
SV AE
!***** 建立混凝土梁单元**** ****** *** *
附录 2
B O K 0 L 0 2 , . LC , 0 . 0 0 2 , , 5 , 0 BO K 0 L0 2 , 50 0 2 LC , , 5 0 3 , . , . . , 0 BO KO L 2 , 5- . 0 LC , 0 50 3 ,0 1, , . . 5 BO K 0 L0 2 , 5 02 , 5 LC , , 50 3 , 0 0 3 , . . . . V E , . S LA I L VL EA I G U, L
SV AE
左翼缘构造筋
L 部构造钢筋划分网格 [E ,, E,6 1, SlS LN, 1, 1 I 6 LE, IE, 71, S LA LN ,1, 1 , 7 LT , , AT 2 3 2 , LS Z , , 5 E IE A L 00 L .
LE IA L M SI L , l lAL 名E , L
K 2 , , 5 0 1, , 3 5 0 0 , 5 4 . . . L T , SR 2 , 1 2 2 L T , SR 2 , 3 2 4
!*****下面纵向钢筋划分网格* ** * ****** * **
LE , , E , 3 1, S L S LN , 1, 2 I 5
R1 ,,
! 定义混凝土的常数
M, , e0弹性模量要根据其他公式求 PE, 3 l! X12 4 .
M , Y 10 2 P P X , . R , 0 M , S l2E 密度 PDN , 53! E ,
钢筋混凝土结构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命令流实例
ansys命令流实例1.结构静力分析--梁分析/PREP7 !进入前处理器ET,1,BEAM3 !定义单元类型R,1,0.25,0.0052,0.5 !定义实常数MP,EX,1,210E6MP,PRXY,1,0.3 !定义材料属性N,1,0N,2,1N,3,5N,4,7N,5,7N,6,9N,7,11N,8,11N,9,13N,10,14 !定义节点E,1,2E,2,3E,3,4E,5,6E,6,7E,8,9E,9,10 !生成单元CP,1,UX,4,5CP,2,UY,4,5CP,3,UX,7,8CP,4,UY,7,8 !耦合节点FINISH/SOLU !进入求解器D,2,UXD,2,UYD,3,UYD,6,UYD,9,UY !施加位移约束F,10,FY,-4 !施加集中约束SFBEAM,1,1,PRES,4,4SFBEAM,2,1,PRES,4,4 !施加均布力SOLVE !求解/POST1 !进入后处理器PLDISP !绘制结构变形图PRDISP !列出各节点的位移ETABLE,IMOMENT,SMISC,6ETABLE,JMOMENT,SMISC,12ETABLE,ISHEAR,SMISC,2ETABLE,JSHEAR,SMISC,8 !将节点弯矩、剪力制表PRETAB !列表显示单元的弯矩、剪力/TITLE,SHEAR FORCE DISTRIBUTION !设置剪力分布图的标题PLLS,ISHEAR,JSHEAR !绘制剪力分布图/TITLE,BENDING MOMENT IDSTRIBUTION !设置弯矩分布图的标题PLLS,IMOMENT,JMOMENT !绘制弯矩分布图流固耦合实例RAD=0.8 !底面半径H=1G=9.8OMEGAR=2ROU=1000 !定义参数变量/PREP7 !进入前处理器ET,1,FLUID79 !选择单元类型KEYOPT,1,3,1 !设置单元关键字MP,EX,1,2E9 !设置杨氏模量MP,DENS,1,ROU !设置材料密度K,1K,2,RADK,3,RAD,HK,4,,H !生成关键点A,1,2,3,4 !连接关键点生成面积LESIZE,ALL,,,10 !设置网格划分精度AMESH,ALL !将面积划分网络/SOLU !进入求解器DL,2,,UXDL,1,,UYNSEL,S,LOC,XDSYM,SYMM,XD,ALL,UXNSEL,ALL !施加位移约束ACEL,,GOMEGA,,OMEGAR !施加惯性力SOLVE !求解/POST1SET,LAST !进入通用后处理器PLNSOL,U,X,0,1 !绘制应力云图UCENT=UY(22)UEDGE=UY(12)UELEV=UEDGE-UCENT !提取节点位移结构静力分析--壳结构内力分析LENGTH=100YOUNG=200000THICKNESS=2FORCE=1000DENSITY=9E-6 !将材质、载荷、板的几何尺寸等参数化/PREP7 !进入前处理器MP,EX,1,YOUNGMP,NUXY,1,0.3MP,DENS,1,DENSITY !定义材质ET,1,SHELL63 !定义单元类型R,1,THICKNESS,THICKNESS,THICKNESS,THICKNESS !定义实常数!构建结构的几何模型K,1,0,0K,2,LENGTH,0K,3,LENGTH,LENGTHK,4,0,LENGTH !定义关键点A,1,2,3,4!通过关键点生成面LSEL,ALLLESIZE,ALL,,,16AMESH,ALL !设定网格划分参数,划分网格FINISH/SOLU !进入求解器NSEL,S,LOC,X,0,0D,ALL,ALL,0 !选择X=0的节点将其固定NSEL,S,LOC,X,LENGTH,LENGTHD,ALL,ALL,0 !选择X=LENGTH的节点将其固定NSEL,S,LOC,X,0.5*LENGTH,0.5*LENGTHNSEL,R,LOC,Y,0.5*LENGTH,0.5*LENGTHF,ALL,FZ,FORCEALLSEL !捕捉板的中心点并在中心点处施加集中力荷载SOLVE !求解FINISH/POST1 !进入后处理器/DSC,,10PLNSOL,U,Z,0,1 !绘图显示板的竖向变形NSEL,ALL !提取板的最大竖向变形NSORT,U,Z,1,1 !将节点的位移绝对值以升序排序*GET,MAXDEFLECTION,SORT,0,MAX !提取位移最大值并赋给变量NSEL,S,LOC,X,0NSEL,A,LOC,X,LENGTH,LENGTH !选择固定边节点NSORT,S,EQV,1,1 !将节点等效应力的绝对值以升序排序*GET,MAXSTRESS,SORT,0,MAX !提取等效应力最大值并赋给变量结构静力分析--平面桁架分析/PREP7 !进入前处理器ET,1,LINK1 !选择单元R,1,0.1 !定义实常数MP,EX,1,30E6MP,PRXY,1,0.3 !定义材料属性N,1,0N,2,4N,3,8N,4,12N,5,0,3N,6,4,3N,7,8,3 !生成节点E,1,2E,2,3E,3,4E,4,7E,3,7E,2,7E,2,6E,2,5E,1,5E,5,6E,6,7 !生成单元FINISH/SOLU !进入求解器D,1,UXD,1,UYD,5,UX !实加位移约束F,2,FY,-15F,3,FY,-15F,4,FY,-15 !施加集中力SOLVE !求解/POST1 !进入能用后处理器PRESOL,FORC !列表显示反力。
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建模apdl命令流实例应用
E,49*25-1+i+c,49*25+i+c *enddo *do,d,1,2*25,1 E,72*25-1+i+d,72*25+i+d *enddo !4 截面 type,1 mat,1 secnum,11 *do,a,1,3*25,1 E,9*25-1+i+a,9*25+i+a *enddo *do,b,1,4*25,1 E,26*25-1+i+b,26*25+i+b *enddo *do,c,1,3*25,1 E,51*25-1+i+c,51*25+i+c *enddo *do,d,1,4*25,1 E,68*25-1+i+d,68*25+i+d *enddo !3 截面 type,1 mat,1 secnum,10 *do,a,1,1*25,1 E,7*25-1+i+a,7*25+i+a *enddo *do,b,1,5*25,1 E,12*25-1+i+b,12*25+i+b *enddo *do,c,1,4*25,1 E,22*25-1+i+c,22*25+i+c *enddo *do,d,1,4*25,1 E,54*25-1+i+d,54*25+i+d *enddo *do,e,1,5*25,1 E,63*25-1+i+e,63*25+i+e *enddo !2 截面 type,1 mat,1 secnum,9 *do,a,1,1*25,1
/nerr,0 ! 将参数读入数据库 finish /clear /PREP7 ! 单元类型编号 et,1,beam44 et,2,shell63 et,3,link8 ! 材料类型编号 Q420/Q370/ 拉索 -1 ,混凝土 C40/C50-2 ,刚臂 -3 mp,ex,1,2.1e11 mp,prxy,1,0.3 mp,dens,1,7850 mp,alpx,1,1e-5
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方向的应力。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
/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.46,5329.65, 5433.84,5538.03,5642.22,5746.41*SET,z(55),5850.6,5954.79,6058.98,6163.17,6267.36,6371.55,6475.74,6579.93,6682.47,6785.01, *SET,z(65),6887.55,6990.09,7092.63,7195.17,7297.71,7400.25,7502.79,7605.33*SET,z(73),7707.87,7810.41,7912.95,8024.28,8135.61,8246.94,8358.27,*SET,z(80),8469.6,8580.93,8680.93,8780.93,8880.93,8980.93,9080.93,9180.93,9280.93,9380.93, 9480.93*SET,z(91),9580.93,9691.93,9802.93,9913.93,10024.93,10135.93,10246.93,10346.93,10446.93, *SET,z(100),10546.93,10646.93,10746.93,10846.93,10946.93,11046.93,11146.93,11246.93,1135 8.26*SET,z(109),11469.59,11580.92,11692.25,11803.58,11914.91,12005.91,*SET,z(115),12096.91,12187.91,12278.91,12369.91,12460.91,12551.91,12642.91,12733.91,1282 4.91,12915.91,13006.91,13125.91*SET,z(127),13244.91,13363.91,13482.91,13601.91,13720.91,13839.91,13939.91,14039.91,*SET,z(135),14139.91,14239.91,14339.91,14439.91,14539.91,14639.91,14739.91,14839.91,1494 2.41,15044.91*SET,z(145),15147.41,15249.91,15352.41,15454.91,15557.41,15659.91,15770.91,15881.91,1599 2.91,16109.91!!定义横隔梁的梯度参数表*set,nd,18*dim, del, array, nd*SET,del(1), 2,5,15,31,47,63,76,82,92,95,98,108,114,125,132,142,*SET,del(17),150,153!!定义初始面座标参数表*set,ni,36*dim, initial, array,ni, 3*SET,initial(1,1),-633,-533,-399.6667,-266.333,-133,0,133.333,266.6666667,400,533.333,633,-63 3,-533,-399.667,-266.3333333,-133,*SET,initial(17,1),0,133.3333333,266.6666667,400,533.3333333,633,-633,-633,633,633,-766,-69 9.5,-633,-566.5,-500,500,566.5,633,*SET,initial(35,1),699.5,766*SET,initial(1,2),0,0,0,0,0,0,0,0,0,0,0,0*SET,initial(12,2),60,60,60,60,60,60,60,60,60,60,60,*SET,Initial(23,2),184,276,184,276,368,368,368,368,368,368,368,368,368,368,*SET,initial(1,3),115.82 ,97.53 ,73.13 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-73.19 ,-97.59 ,-115.82 , *SET,initial(12,3),115.82 ,97.53 ,73.13 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-73.19 ,-97.59 ,-115.82 , 115.82 ,*SET,initial(24,3),115.82 ,-115.82 ,-115.82 ,140.16 ,127.99 ,115.82 ,103.66 ,91.49 ,-91.49 ,-103.66 ,-115.82 ,-127.99 ,-140.16!定义横隔梁参数表*set,nv,25*dim,vertical,array,nv,3*SET,vertical(1 ,3),97.53 ,73.13 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-73.19 ,-97.59 ,97.53 ,73.13 ,48 .73 ,24.34*SET,vertical(14,3),0.00 ,-24.40 ,-48.79 ,-73.19 ,-97.59 ,73.19 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-73.19*SET,vertical(1,2),184,184,184,184,184,184,184,184,184,*SET,vertical(10,2),276,276,276,276,276,276,276,276,276,*SET,vertical(19,2),368,368,368,368,368,368,368*SET,vertical(1,1),-533,-399.667 ,-266.333 ,-133 ,0 ,133.333 ,266.667 ,400.000 ,533.333 ,-533.000 ,-399.667 ,-266.333 ,-133*SET,vertical(14,1),0,133.333 ,266.667 ,400.000 ,533.333 ,-400.000 ,-266.333 ,-133.000 ,0.000 ,13 3.333 ,266.667 ,400.000!!定义混凝土的节点参数表*set,nc,135*dim,concrete,array,nc,3*SET,concrete(1,1),-1301.625 ,-1301.625 ,-1301.625 ,-1248.225 ,-1248.225 ,-1248.225 ,-1194.825 ,-1194.825 ,-1194.825 ,-1141.425*SET,concrete(11,1),-1141.425 ,-1141.425 ,-1088.025 ,-1088.025 ,-1088.025 ,-1034.625 ,-1034.62 5 ,-1034.625 ,-981.225 ,-981.225*SET,concrete(21,1),-981.225 ,-927.825 ,-927.825 ,-927.825 ,-874.425 ,-874.425 ,-874.425 ,-821.0 25 ,-821.025 ,-821.025*SET,concrete(31,1),-766.000 ,-766.000 ,-766.000 ,-699.500 ,-699.500 ,-699.500 ,-633.000 ,-633.0 00 ,-633.000 ,-566.500*SET,concrete(41,1),-566.500 ,-566.500 ,-500.000 ,-500.000 ,-500.000 ,-450.000 ,-450.000 ,-450.0 00 ,-400.000 ,-400.000*SET,concrete(51,1),-400.000 ,-333.167 ,-350.000 ,-333.167 ,-266.333 ,-250.000 ,-266.333 ,-199.6 67 ,-199.667 ,-199.667*SET,concrete(61,1),-133.000 ,-133.000 ,-133.000 ,-66.500 ,-66.500 ,-66.500 ,0.000 ,0.000 ,0.000 , 66.500*SET,concrete(71,1),66.500 ,66.500 ,133.000 ,133.000 ,133.000 ,199.667 ,199.667 ,199.667 ,266. 333 ,250.000*SET,concrete(81,1),266.333 ,333.167 ,350.000 ,333.167 ,400.000 ,400.000 ,400.000 ,450.000 ,45 0.000 ,450.000*SET,concrete(91,1),500.000 ,500.000 ,500.000 ,566.500 ,566.500 ,566.500 ,633.000 ,633.000 ,63 3.000 ,699.500*SET,concrete(101,1),699.500 ,699.500 ,766.000 ,766.000 ,766.000 ,821.025 ,821.025 ,821.025 ,8 74.425 ,874.425*SET,concrete(111,1),874.425 ,927.825 ,927.825 ,927.825 ,981.225 ,981.225 ,981.225 ,1034.625 , 1034.625 ,1034.625*SET,concrete(121,1),1088.025 ,1088.025 ,1088.025 ,1141.425 ,1141.425 ,1141.425 ,1194.825 ,1 194.825 ,1194.825 ,1248.225*SET,concrete(131,1),1248.225 ,1248.225 ,1301.625 ,1301.625 ,1301.625*SET,concrete(1, 2),468.25,443.219 ,418.188 ,468.250 ,440.716 ,413.181 ,468.250 ,443.219 ,408.175 ,468.250*SET,concrete(11, 2), 443.219 ,403.169 ,468.250 ,430.703 ,398.163 ,468.250 ,430.703 ,393.156 ,468.250 ,418.188*SET,concrete(21, 2), 388.150 ,468.250 ,418.188 ,383.144 ,468.250 ,418.188 ,378.138 ,468.250 ,418.188 ,373.131*SET,concrete(31, 2), 468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250*SET,concrete(41, 2), 418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188*SET,concrete(51, 2), 368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000*SET,concrete(61, 2), 468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250*SET,concrete(71, 2), 418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188*SET,concrete(81, 2), 368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000*SET,concrete(91, 2), 468.25,418.188 ,368.000 ,468.25,418.188,368.000 ,468.25,418.188,368,468.25*SET,concrete(101, 2), 418.188,368,468.25,418.188,368,468.250 ,418.188,373.131,468.25,418.188*SET,concrete(111, 2), 378.138,468.25,418.188,383.144,468.25,418.188 ,388.15,468.25,430.703,393.156*SET,concrete(121, 2), 468.25,430.703,398.163,468.25,443.219,403.169,468.25,443.219,408.175,468.25*SET,concrete(131, 2), 440.716,413.181,468.25,443.219,418.188,*SET,concrete(1, 3), 238.1649941,238.165 ,238.165 ,228.394 ,228.394 ,228.394 ,218.623 ,218.623 ,218.623 ,208.852 *SET,concrete(11, 3), 208.852 ,208.852 ,199.082 ,199.082 ,199.082 ,189.311 ,189.311 ,189.311 ,179.540 ,179.540*SET,concrete(21, 3), 179.540 ,169.769 ,169.769 ,169.769 ,159.998 ,159.998 ,159.998 ,150.227 ,150.227 ,150.227*SET,concrete(31, 3), 140.159 ,140.159 ,140.159 ,127.991 ,127.991 ,127.991 ,115.823 ,115.823 ,115.823 ,103.655*SET,concrete(41, 3), 103.655 ,103.655 ,91.488 ,91.488 ,91.488 ,82.339 ,82.339 ,82.339 ,73.190 ,73.190*SET,concrete(51, 3), 73.190 ,60.961 ,64.041 ,60.961 ,48.732 ,45.744 ,48.732 ,36.534 ,36.534 ,36.534*SET,concrete(61, 3), 24.336 ,24.336 ,24.336 ,12.168 ,12.168 ,12.168 ,0.000 ,0.000 ,0.000 ,-12.168*SET,concrete(71, 3), -12.168 ,-12.168 ,-24.336 ,-24.336 ,-24.336 ,-36.534 ,-36.534 ,-36.534 ,-48.732 ,-45.744*SET,concrete(81, 3), -45.744 ,-60.961 ,-64.041 ,-60.961 ,-73.190 ,-73.190 ,-73.190 ,-82.339 ,-82.339 ,-82.339*SET,concrete(91, 3), -91.48756137,-91.488 ,-91.488 ,-103.65,-103.65,-103.655 ,-115.8232527,-115.8232527,-115.82,-127.99*SET,concrete(101, 3), -127.991,-127.9910,-140.159,-140.159,-140.159,-150.227 ,-150.23,-150.23,-159.9980217,-159.99 80217*SET,concrete(111, 3), -159.998,-169.768,-169.768,-169.768,-179.539,-179.540 ,-179.5397,-189.3106,-189.3106364,-18 9.3106364*SET,concrete(121, 3), -199.0815,-199.0815,-199.0815,-208.852,-208.8523,-208.8523,-218.62,-218.623,-218.623251,-22 8.3941226*SET,concrete(131, 3), -228.3941226,-228.3941226,-238.1649941,-238.1649941,-238.1649941,!节点编号,生成箱梁的节点*do, i, 1,1,1*do, j, 1, nin, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)*enddo*enddo*do, i, nz,nz,1*do, j, 1, nin, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)*enddo*enddo*do, i, 2,nz-1,1*do, j, 1, nin, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)*enddo*enddo!生成横隔板节点*set, infl, ni*(nz+1)*set, nfl, 39*do, i, 1, 18, 1*SET,k,del(i)*do, j, 1, 25, 1n, infl+i*nfl+j, vertical(j,1), vertical(j,2), vertical(j,3)+z(k)*enddo*do, j, 26, 32n, infl+i*nfl+j, vertical(j-7,1), vertical(j-7,2), vertical(j-7,3)+z(k-1)*enddo*do, j, 33, nfln, infl+i*nfl+j, vertical(j-14,1), vertical(j-14,2), vertical(j-14,3)+z(k+1)*enddo*enddo!节点编号,生成混凝土的节点*set,incr,infl+18*nfl+1*do, i, 1, 1,*do, j, 1, ncn, incr+j+nc*i, concrete(j,1),concrete(j,2),concrete(j,3)+z(i)*enddo*enddo*do, i, 2, nz*do, j, 1, ncn, incr+j+nc*i, concrete(j,1),concrete(j,2),concrete(j,3)+z(i)*enddo*enddo!生成钢绞线的节点*set,insty1,incr+nc*(nz+1)!Y1*do, i, 37, 127, 1n, insty1+(i-36), concrete(53,1),concrete(53,2),concrete(53,3)+z(i)n, insty1+(i-36)+127-37+1, concrete(68,1),concrete(68,2),concrete(68,3)+z(i) n, insty1+(i-36)+2*(127-37)+2, concrete(83,1),concrete(83,2),concrete(83,3)+z(i) *enddo!Y2*set,insty2,insty1+3*(127-37+1)*do, i, 66,114,1n, insty2+(i-65), concrete(56,1),concrete(56,2),concrete(56,3)+z(i)n, insty2+(i-65)+141-66+1, concrete(80,1),concrete(80,2),concrete(80,3)+z(i)*enddo!生成钢绞线锚固区节点*set,instb1,insty2+2*(141-66)*set,instb2,instb1+40!Y1*do, i, 1, 4*SET,k,27+in, instb1+i, -350.000 , 298.188+12*i, 64.0413+z(k)n, instb2+i, -350.000 , 298.188+12*i, 64.0413+z(k)*enddo*do, i, 1, 4*SET,k,27+in, instb1+i+4, 0 , 298.188+12*i, z(k)n, instb2+i+4, 0 , 298.188+12*i, z(k)*enddo*do, i, 1, 4*SET,k,27+in, instb1+i+8, 350.000 , 298.188+12*i, -64.0413+z(k)n, instb2+i+8, 350.000 , 298.188+12*i, -64.0413+z(k)*enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+12, -350.000 , 418.188-12*(i+6), 64.0413+z(k) n, instb2+i+12, -350.000 , 418.188-12*(i+6), 64.0413+z(k) *enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+16, 0 , 418.188-12*(i+6), z(k)n, instb2+i+16, 0 , 418.188-12*(i+6), z(k)*enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+20, 350.000 , 418.188-12*(i+6), -64.0413+z(k) n, instb2+i+20, 350.000 , 418.188-12*(i+6), -64.0413+z(k) *enddo!Y2*do, i, 1, 4*SET,k,56+in, instb1+i+24, -250.000 , 298.188+12*i, 45.743+z(k)n, instb2+i+24, -250.000 , 298.188+12*i, 45.743+z(k)*enddo*do, i, 1, 4*SET,k,56+in, instb1+i+28, 250.000 ,298.188+12*i, -45.743+z(k)n, instb2+i+28, 250.000 ,298.188+12*i, -45.743+z(k)*enddo*do, i, 1, 4*SET,k,115+i+6n, instb1+i+32, -250.000 ,418.188-12*(i+6), 45.743+z(k) n, instb2+i+32, -250.000 ,418.188-12*(i+6), 45.743+z(k) *enddo*do, i, 1, 4*SET,k,115+i+6n, instb1+i+ni, 250.000 ,418.188-12*(i+6), -45.743+z(k) n, instb2+i+ni, 250.000 ,418.188-12*(i+6), -45.743+z(k) *enddo!生成箱梁的单元*do, i, 1,nz-1,1!底板,shell63,钢,厚度8type,1mat,1real,2*do, j, 1, 10, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j,*enddo!顶板,shell63,钢,厚度7type,1mat,1real,1*do, j, 27, 30, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j,*enddo*do, j, 32, 35, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j,*enddo!腹板、底板加劲肋,shell63,钢,厚度4type,1mat,1real,3*do, j, 1, 11, 1e, i*ni+j, (i+1)*ni+j, (i+1)*ni+j+11, i*ni+j+11*enddoe, i*ni+12, (i+1)*ni+12, (i+1)*ni+23, i*ni+23e, i*ni+23, (i+1)*ni+23, (i+1)*ni+24, i*ni+24e, i*ni+24, (i+1)*ni+24, (i+1)*ni+29, i*ni+29e, i*ni+22, (i+1)*ni+22, (i+1)*ni+25, i*ni+25e, i*ni+25, (i+1)*ni+25, (i+1)*ni+26, i*ni+26e, i*ni+26, (i+1)*ni+26, (i+1)*ni+34, i*ni+34*enddo!生成横隔板单元type,1mat,1real,3*do, i, 1, 18, 1!生成横隔板单元*do, l, 0, 9, 1e, del(i)*ni+1+l, del(i)*ni+2+l, del(i)*ni+13+l, del(i)*ni+12+l, *enddo*do, l, 0, 7,1e, infl+i*nfl+l+1, infl+i*nfl+l+2, infl+i*nfl+l+11, infl+i*nfl+l+10,e, del(i)*ni+l+13,del(i)*ni+l+14, infl+i*nfl+l+2, infl+i*nfl+l+1,*enddo*do, l, 0, 5, 1e, infl+i*nfl+l+11, infl+i*nfl+l+12, infl+i*nfl+l+20, infl+i*nfl+l+19,*enddo*do, l, 0, 5, 1e, infl+i*nfl+l+19, infl+i*nfl+l+20, infl+i*nfl+l+27, infl+i*nfl+l+26e, infl+i*nfl+l+33, infl+i*nfl+l+34, infl+i*nfl+l+20, infl+i*nfl+l+19*enddoe, del(i)*ni+31, infl+i*nfl+19, infl+i*nfl+26, (del(i)-1)*ni+31,e, (del(i)+1)*ni+31, infl+i*nfl+33, infl+i*nfl+19, del(i)*ni+31,e, infl+i*nfl+25, del(i)*ni+32, (del(i)-1)*ni+32, infl+i*nfl+32,e, infl+i*nfl+nfl, (del(i)+1)*ni+32, del(i)*ni+32, infl+i*nfl+25,e, del(i)*ni+12, del(i)*ni+13, infl+i*nfl+1,del(i)*ni+23e, del(i)*ni+23, infl+i*nfl+1, infl+i*nfl+10, del(i)*ni+24,e, del(i)*ni+24, infl+i*nfl+10, del(i)*ni+31, del(i)*ni+29,e, infl+i*nfl+10, infl+i*nfl+11, infl+i*nfl+19, del(i)*ni+31,e, del(i)*ni+21, del(i)*ni+22, del(i)*ni+25, infl+i*nfl+9,e, infl+i*nfl+9, del(i)*ni+25, del(i)*ni+26, infl+i*nfl+18,e, infl+i*nfl+18, del(i)*ni+26, del(i)*ni+34, del(i)*ni+32,e, infl+i*nfl+17, infl+i*nfl+18, del(i)*ni+32, infl+i*nfl+25,*enddo!生成混凝土单元type,2mat,2real, 4*do, i, 1, 153, 1*do, j, 0, 130, 3*do, k, 2, 3, 1e,incr+nc*i+j+k,incr+nc*(i+1)+j+k,incr+nc*(i+1)+j+k+3,incr+nc*i+j+k+3,incr+nc*i+j+k-1,incr+nc*(i+ 1)+j+k-1,incr+nc*(i+1)+j+k+2, incr+nc*i+j+k+2,*enddo*enddo*enddo!生成Y1钢绞线单元type,3mat,1real, 5*do, i, 37, 126, 1e, insty1+i-36,insty1+i-ni+1e, insty1+i-36+127-37+1, insty1+i-ni+127-37+2e, insty1+(i-36)+2*(127-37)+2,insty1+(i-ni)+2*(127-37)+3*enddoe, instb1+4, insty1+1,e, instb1+8, insty1+1+127-37+1e, instb1+12, insty1+1+2*(127-37)+2e, instb1+13, insty1+(127-36)e, instb1+17, insty1+(127-36)+127-37+1,e, instb1+21, insty1+(127-36)+2*(127-37)+2*do, i, 1, 3*do, j ,0, 5e, instb1+i+j*4, instb1+i+j*4+1*enddo*enddofinish/SOL!进入求解器施加荷载ESEL,S,TYPE,,3 !选择现有单元类型为LINK8的单元EPLOTBFE,all,TEMP,1,tem1 !对目前显示的所有单元施加温度荷载,大小为-tem1 ESEL,none !全不选择/PREP7!进入前处理器,生成Y2钢绞线的单元*set,insty2,insty1+3*(127-37+1)*do, i, 66,114,1n, insty2+(i-65), concrete(56,1),concrete(56,2),concrete(56,3)+z(i)n, insty2+(i-65)+141-66+1, concrete(80,1),concrete(80,2),concrete(80,3)+z(i) *enddo*do, i, 66, 113, 1e, insty2+i-65,insty2+i-65+1e, insty2+(i-65)+141-66+1, insty2+(i-65)+141-66+2*enddoe, instb1+28, insty2+1e, instb1+32, insty2+1+141-66+1e, instb1+33, insty2+(113-65)e, instb1+37,insty2+(113-65)+141-66+1*do, i, 1, 3e, instb1+i+j*4+24, instb1+i+j*4+1+24*enddo*enddoFINISH/SOL!进入求解器施加荷载BFE,all,TEMP,1,tem2 !对目前选择的单元即Y2施加温度荷载,大小为-tem2ESEL,all !选择所有的单元FINISH/PREP7!生成钢筋锚固区type,2mat,2real, 4*do, i, 1, 3e, instb2+i, instb2+4+i, instb2+5+i, instb2+1+i, incr+54+nc*(i+27),incr+69+nc*(i+27),incr+69+nc*(i+28), incr+54+nc*(i+28),e, instb2+4+i,instb2+8+i, instb2+9+i,instb2+5+i, incr+nc*(i+27)+69,incr+nc*(i+27)+84, incr+nc*(i+28)+84, incr+nc*(i+28)+69e, instb2+12+i, instb2+16+i, instb2+17+i, instb2+13+i, incr+54+nc*(i+134),incr+69+nc*(i+134),incr+69+nc*(i+135),incr+54+nc*(i+135)e, instb2+16+i, instb2+20+i, instb2+21+i, instb2+17+i, incr+69+nc*(i+134),incr+84+nc*(i+134),incr+84+nc*(i+135),incr+69+nc*(i+135)*enddo*do, i,1, 3e, instb2+24+i, instb2+28+i, instb2+29+i, instb2+25+i,incr+57+nc*(i+56),incr+81+nc*(i+56), incr+81+nc*(i+57),incr+57+nc*(i+57)e, instb2+32+i, instb2+36+i, instb2+37+i, instb2+33+i,incr+57+nc*(i+121),incr+81+nc*(i+121), incr+81+nc*(i+122),incr+57+nc*(i+122)*enddo!加载FINISH/SOL*do, i, 1, 2*do, j, 4,8d, ni*i+j, all,*enddo*do, i, 94, 96 *do, j, 4,8d, ni*i+j, all, *enddo*enddo*do, i, 153, 154 *do, j, 4,8d, ni*i+j, ux, *enddo*enddo*do, i, 153, 154 *do, j, 4,8d, ni*i+j, uz, *enddo*enddo FINISH/PREP7 cpintf, ux, 0.01 cpintf, uy, 0.01 cpintf, uz, 0.01nsel, alleplot/solusolve。