ANSYS在土木工程中的应用

ANSYS在土建结构中的应用分析

ANSYS在土建结构中的应用分析摘要：本文详细介绍了ANSYS软件在预应力结构中的应用，并通过与实验结果的对比，论证了ANSYS软件可以替代具体实验，进行结构性能分析。

关键词：有限元；ANSYS；预应力；无粘结前言长期以来钢筋混凝土结构的设计方法是用线弹性理论来研究的应力或内力，显然，这是不太合理的，特别是混凝土，成分复杂，表现出明显的非线性行为。

而我们知道并非所有的结构都能通过试验来取得数据的，所以有限元分析软件便在研究中显示出了优越性。

ANSYS是当今国际主流的有限元分析工具，本文详细介绍了ANSYS在土建结构中的应用。

1 试验内容1.1 试件梁的截面为矩形，尺寸为200mm×400mm，跨长3000mm，折线布筋形式，钢管抽芯成孔。

混凝土强度为C40。

1.2 主要测试内容为了了解预制预应力连续梁的内力重分布，主要测量支座反力，支座位移和跨中挠度、混凝土的应变变化、截面曲率变化等。

采用逐级加载。

2 ANSYS有限元模型的具体建模过程2.1 单元类型的定义及实常数定义SOLID65单元模拟混凝土，LINK8单元模仿预应力筋，SOLID65单元中的自定义配筋模拟普通钢筋，COMBIN39非线性弹簧单元模拟钢筋和混凝土之间的滑移。

接触单元模拟支座端面与梁端面接缝关系，面与面接触形式，Conta174为接触面单元，Targe170为目标单元。

2.2 本模型所采用的本构关系混凝土的本构关系为单轴应力-应变关系。

预应力钢筋本构模型采用双线性随动强化（BKIN），其包括弹性斜率和塑性斜率。

2.3 混凝土建模混凝土模型采用了自上而下和自下而上两种建模方式相结合的办法来进行创建。

2.4 普通钢筋和无粘结预应力筋的建模为方便建模和计算，无粘结预应力钢筋采用单独实体建模，用SOLID65自带的自定义配筋方式为普通钢筋和箍筋进行弥散式配筋。

2.5 混凝土与预应力钢筋间关系的建立混凝土和预应力钢筋之间关系采用双弹簧单元来模拟，把预应力钢筋节点和与对应的最近的混凝土节点连接起来，分别设置成切向和法向的弹簧即双弹簧，耦合另一方向节点自由度。

ANSYS土木工程应用实例命令流(第二版)》

ANSYS土木工程应用实例命令流（第二版）》!隧道模拟开挖命令流!希望能够得到加分!!土体的弹性摸量增加/comfini/cle*set,x1,-12*set,y1,-12*set,w1,28.9*set,h1,30.15*set,x2,-25*set,y2,-12*set,w2,13*set,h2,30.15*set,x3,16.9*set,y3,-12*set,w3,13*set,h3,30.15*set,x4,-25*set,y4,-30*set,w4,54.9*set,h4,18*set,th,0.4 !厚度*set,length_z,50 !洞的进深 /prep7et,1,mesh200,2 !用于3-D的2节点线et,2,mesh200,6 !用于3-D的4节点四边形et,3,shell63et,4,solid45r,1,th !壳的厚度mp,ex,1,3.0e10 !壳的材料,C30混凝土mp,prxy,1,0.2 mp,dens,1,2500 mp,ex,2,4.5e8 !保留岩石的材料mp,prxy,2,0.32 !tb,dp,2!tbdata,1,20,30, mp,dens,2,2700mp,ex,3,4.51e8 !挖去岩石的材料 mp,prxy,3,0.32 !tb,dp,3!tbdata,1,20,30, mp,dens,3,2700k,,0,0k,,0,3.85k,,0.88,5.5k,,2.45,6.15k,,4.02,5.5k,,4.9,3.85k,,4.9,0larc,1,2,6,8.13 !定义两点之间的圆弧线，larc,p1,p2,pc,radlarc,2,3,6,3.21 larc,3,4,6,2.22 larc,4,5,2,2.22 larc,5,6,2,3.21 larc,6,7,2,8.13 larc,7,1,4,6a,1,2,3,4,5,6,7 !产生面1!blc4,x1,y1,w,h1 !block4,xcorner,ycorner,width,height,depth!blc4,x2,y2,w,h2 !block4,xcorner,ycorner,width,height,depth 产生面3 blc4,x1,y1,w1,h1 !产生面2 blc4,x2,y2,w2,h2 !产生面3blc4,x3,y3,w3,h3 !产生面4 blc4,x4,y4,w4,h4 !产生面5 /pnum,area,1 aovl,1,2,3,4,5 !布尔操作重叠，得到面3 nummrg,all,,,,low numcmp,all l,1,8 !从四个角点上连接出四条线 l,7,9l,6,10l,2,11lsel,s,line,,21,22,1 !用线分割面 lsel,a,line,,7asbl,5,alllsel,s,line,,21,24,3 lsel,a,line,,1asbl,7,alllsel,s,line,,22,23,1 lsel,a,line,,6asbl,8,allnummrg,all,,,,low numcmp,alllsel,s,line,,2,5,1 LCCAT,alllesize,all,,,3lsel,s,line,,9,11,2 lsel,a,line,,6lsel,a,line,,1lesize,all,,,8lsel,s,line,,8,10,2 lsel,a,line,,7lesize,all,,,12lsel,s,line,,21,24,1 lesize,all,,,10,2type,2asel,s,area,,5,8,1 amesh,all!mopt,split,on !指定网格划分选项，四边形单元如果违背了单元形状则将其分解为三角形单元!mopt,split,err !mopt,split,warn !mopt,qmesh,main !选用四边形网格划分器，当这种操作失败时，不调用备用的网格划分器asel,s,area,,1amesh,1lsel,s,line,,12,13,1 lesize,all,,,8lsel,s,line,,15,18,1 lesize,all,,,6,2 asel,s,area,,2,3,1 amesh,all lsel,s,line,,14 lesize,all,,,24 lsel,s,line,,19,20,1lesize,all,,,6,2 lsel,s,line,,15,17,2 lsel,a,line,,8LCCAT,allasel,s,area,,4amesh,allLSEL,s,LCCALDELE,allnummrg,all,,,,low numcmp,allallsel!以下开始拉伸成实体单元!首先拉伸成壳单元k,1000,,,-length_z l,1,1000/view,1,1,1,1/replotEXTOPT,ESIZE,10,0, LSEL,S,LINE,,1,7,1 ADRAG,all,,,,,,25 gplot type,3real,1mat,1ASEL,S,loc,z,-25 APLOTlsel,s,loc,z,-25 lesize,all,,,10 MSHAPE,0,2DMSHKEY,1amesh,all!拉伸岩石的实体ASEL,invertaplotEXTOPT,ESIZE,10,0, EXTOPT,ACLEAR,1 TYPE,4MAT,2asel,r,area,,2,8,1 VDRAG,all,,,,,,25 allsel!挖去部分岩石的实体MAT,3VDRAG,1,,,,,,25 EPLOTnummrg,all,,,,low numcmp,all!约束两侧面的X方向的约束asel,s,loc,x,x2 asel,a,loc,x,x2+w4 da,all,ux,0alls!约束地面的Y方向的约束asel,s,loc,y,y4 da,all,uy,0allsasel,s,loc,z,-length_zasel,a,loc,z,0 da,all,uz,0allselacel,,10fini/soluantype,staticdeltim,0.1,0.05,0.2 autots,on !使用自动时间步 pred,on !打开时间步长预测器 lnsrch,on !打开线性搜索nlgeom,on !打开大位移效果nropt,full !设定牛顿-拉普森选项cnvtol,f,,0.02,2,0.5esel,s,type,,3 !选择梁单元，杀死ekill,allesel,allesel,s,live !选择活的单元nsle,s !选择活单元上的节点nsel,invert !反向选择，即选择了死单元上的节点 d,all,all,0 !将死单元上的节点约束所有位移，使其不参与矩阵运算nsel,allesel,all/PBC,ALL,,1gplotsolve*do,ii,1,10,1esel,s,mat,,3 !选择挖去的岩石，杀死nsle,snsel,r,loc,z,0.1-(ii-1)*5,-(5.1+(ii-1)*5)esln,r,1ekill,allesel,s,type,,3 !激活挖去的岩石对应的壳单元,并将其节点上的约束删除nsle,snsel,r,loc,z,0.1-(ii-1)*5,-(5.1+(ii-1)*5)esln,r,1ealive,allnsle,sddele,all,allesel,allesel,s,live !选择活单元，此时应该包含两部份，一是梁单元，二是未挖去的岩石单元nsle,snsel,invert !反向选择，将死单元上的节点约束所有自由度d,all,all,0nsel,allesel,allsolve*enddo/post1/DEVICE,VECTOR,1 ESEL,S,type,,3 SET,1,LAST,1, PLNSOL,U,Y,0,1SET,2,LAST,1, PLNSOL,U,Y,0,1 SET,6,LAST,1, PLNSOL,U,Y,0,1 SET,11,LAST,1, PLNSOL,U,Y,0,1!*dim ,VUY_1!*dim ,VUY_2!*dim ,VUY_6!*dim ,VUY_11!SET,1,LAST,1, !PLNSOL,U,Y,0,1 !PLNSOL,U,SUM,0,1 !ETABLE,EUY_1,U,Y ! PLETAB,EUY_1,AVG !PRETAB,EUY_1!*vget,VUY_1,node,387,U,Y!SET,2,LAST,1, !PLNSOL,U,Y,0,1 !PLNSOL,U,SUM,0,1 !ETABLE,EUY_2,U,Y ! PLETAB,EUY_2,AVG !PRETAB,EUY_2!*vget,VUY_2,node,387,U,Y!SET,6,LAST,1, !PLNSOL,U,Y,0,1 !PLNSOL,U,SUM,0,1 !ETABLE,EUY_6,U,Y ! PLETAB,EUY_6,AVG !PRETAB,EUY_6!*vget,VUY_6,node,323,U,Y!SET,3,LAST,1,!PLNSOL,U,Y,0,1!PLESOL,SMISC,4 !MX !PLESOL,SMISC,5 !MY !PLESOL,SMISC,6 !MXY!ETABLE,p_i,SMISC, 13 !ETABLE,p_j,SMISC, 14 !ETABLE,p_k,SMISC,15 !ETABLE,p_l,SMISC, 16!ETABLE,p_i_3,SMISC, 17 !ETABLE,p_j_3,SMISC, 18!ETABLE,p_k_4,SMISC, 20 !ETABLE,p_l_4,SMISC, 19!ETABLE,p_k_5,SMISC, 21 !ETABLE,p_l_5,SMISC, 22 !fini!/post1!etable,if,smisc,1 !etable,jf,smisc,7 !etable,im,smisc,6 !etable,jm, smisc,12大家可以在这个的基础上进一步发掘。

ANSYS桥梁工程应用实例分析(详细)(图文)

本章介绍桥梁结构的模拟分析。

桥梁是一种重要的工程结构，精确分析桥梁结构在各种受力方式下的响应有较大的工程价值。

模拟不同类型的桥梁需要不同的建模方法，分析内容包括静力分析、动荷载响应分析、施工过程分析等等。

在本章中着重介绍桁架桥、刚架桥和斜拉桥三种类型桥梁。

内容提要第6章 ANSYS 桥梁工程应用实例分析本章重点结构分析具体步骤结构静力分析桁架结构建模方法结构模态分析本章典型效果图6.1 引言ANSYS通用有限元软件在土木工程应用分析中可发挥巨大的作用。

我们用它来分析桥梁工程结构,可以很好的模拟各种类型桥梁的受力、施工工况、动荷载的耦合等。

ANSYS程序有丰富的单元库和材料库,几乎可以仿真模拟出任何形式的桥梁。

静力分析中,可以较精确的反应出结构的变形、应力分布、内力情况等；动力分析中,也可精确的表达结构的自振频率、振型、荷载耦合、时程响应等特性。

利用有限元软件对桥梁结构进行全桥模拟分析,可以得出较准确的分析结果。

本章介绍桥梁结构的模拟分析。

作为一种重要的工程结构,桥梁的精确分析具有较大的工程价值。

桥梁的种类繁多,如梁桥、拱桥、钢构桥、悬索桥、斜拉桥等等,不同类型的桥梁可以采用不同的建模方法。

桥梁的分析内容又包括静力分析、施工过程模拟、动荷载响应分析等。

可以看出桥梁的整体分析过程比较复杂。

总体上来说,主要的模拟分析过程如下:(1) 根据计算数据,选择合适的单元和材料,建立准确的桥梁有限元模型。

(2) 施加静力或者动力荷载,选择适当的边界条件。

(3) 根据分析问题的不同,选择合适的求解器进行求解。

(4) 在后处理器中观察计算结果。

(5) 如有需要,调整模型或者荷载条件,重新分析计算。

桥梁的种类和分析内容众多,不同类型桥梁的的分析过程有所不同,分析侧重点也不一样。

在这里仅仅给出大致的分析过程,具体内容还要看具体实例的情况。

6.2 典型桥梁分析模拟过程6.2.1 创建物理环境建立桥梁模型之前必须对工作环境进行一系列的设置。

ANSYS土木工程毕业设计计算书

1目录1.工程概况........................................................................................................................................... 3 2. 分析依据.......................................................................................................................................... 3 3荷载与计算工况................................................................................................................................ 3 3.1荷载............................................................................................................................................. 3 3.2计算工况..................................................................................................................................... 4 4有限元模型的建立............................................................................................................................ 4 4.1 基本假定.................................................................................................................................... 4 4.2几何模型..................................................................................................................................... 5 4.3 主要物理参数取值.................................................................................................................... 6 4.4 单元选取.................................................................................................................................... 6 4.5 建模与分网................................................................................................................................ 6 4.5.1 建模思路........................................................................................................................... 6 4.5.2分网.................................................................................................................................... 7 5 空间结构模型分析........................................................................................................................... 8 5.1刚度校核..................................................................................................................................... 8 5.1.1工况1 可变荷载组合....................................................................................................... 8 5.1.2 工况2水平地震力........................................................................................................... 8 5.2强度校核................................................................................................................................... 10 5.2.1 工况1 可变荷载组合.................................................................................................... 10 5.2.2 工况2水平地震力......................................................................................................... 11 6.平面结构模型分析.......................................................................................................................... 12 6.1 刚度校核.................................................................................................................................. 12 6.1.1 工况1可变荷载组合..................................................................................................... 12 6.1.2工况2 水平地震力......................................................................................................... 13 6.2 强度校核.................................................................................................................................... 14 6.2.1工况1 可变荷载组合.. (14)26.2.2 工况2 水平地震力........................................................................................................ 15 6.2.3 工况3 重力荷载代表值组合........................................................................................ 18 6.3 抗震内力组合.......................................................................................................................... 19 6.3.1 轴力组合......................................................................................................................... 19 6.3.2弯矩组合.......................................................................................................................... 19 7.solid65单元建模 ............................................................................................................................. 20 8.番外篇基于ANSYS 、PKPM 、手算的误差分析....................................................................... 21 8.1 计算原理.................................................................................................................................. 21 8.2PKPM 与 ANSYS 的优点与不足 (21)3RC 框架结构分析1.工程概况工程名称：西安某省技术研发中心办公楼结构设计建设规模：总建筑面积约为4142.76m 2，主体为六层框架结构体系（无地下室），结构层高为21.2m ，底层结构高度为4.7m ，其它层层高3.3米。

ansys workbench2020工程实例解析

Ansys Workbench 2020是一款强大的工程仿真软件，广泛应用于工程领域的结构、流体、热传导等多个领域的仿真分析。

本文将以Ansys Workbench 2020为工具，通过几个典型的工程实例，解析其在工程实践中的应用和优势，帮助读者更好地了解和使用该软件。

1. 车身结构优化在汽车制造领域，车身结构的设计和优化是一个复杂而又关键的问题。

通过Ansys Workbench 2020的结构分析模块，可以对车身结构进行强度、刚度、振动等方面的仿真分析，进而优化结构设计，提高车身的整体性能和安全性。

通过对车身材料、连接结构、受力情况等多个方面的仿真分析，工程师可以更好地指导实际设计，提高设计效率和成功率。

2. 风力发电机叶片设计风力发电机的叶片设计是风力发电领域的核心问题之一。

Ansys Workbench 2020的流体仿真模块可以对风力发电机叶片的气动性能进行仿真分析，包括气动力、气流分布等多个方面的参数。

通过对叶片的材料、形状、尺寸等进行仿真分析和优化，可以提高风力发电机的发电效率和稳定性，降低能量损耗，对提高风力发电机的整体性能具有重要意义。

3. 燃烧室热传导分析在航天、航空发动机等领域，燃烧室的热传导分析是一个关键的问题。

Ansys Workbench 2020的热传导分析模块可以对燃烧室内部的温度场、热应力等进行仿真分析，帮助工程师优化燃烧室的结构设计、材料选择和冷却系统设计。

通过仿真分析，可以提高燃烧室的工作效率和寿命，确保燃烧室的安全可靠性。

4. 桥梁结构静动力分析在土木工程领域，桥梁结构的设计和分析是一个重要的问题。

Ansys Workbench 2020的静动力分析模块可以对桥梁结构在静载荷和动载荷作用下的响应进行仿真分析，包括应力、挠度、疲劳寿命等多个方面的参数。

通过仿真分析，工程师可以对桥梁的结构设计、材料选择和荷载标准进行优化，确保桥梁的安全可靠性和经济性。

Ansys Workbench 2020作为一款强大的工程仿真软件，在工程实践中具有广泛的应用前景和优势。

ANSYS土木工程应用实例命令流(第二版)》

ANSYS土木工程应用实例命令流（第二版）》!隧道模拟开挖命令流!希望能够得到加分!!土体的弹性摸量增加/comfini/cle*set,x1,-12*set,y1,-12*set,w1,28.9*set,h1,30.15*set,x2,-25*set,y2,-12*set,w2,13*set,h2,30.15*set,x3,16.9*set,y3,-12*set,w3,13*set,h3,30.15*set,x4,-25*set,y4,-30*set,w4,54.9*set,h4,18*set,th,0.4 !厚度*set,length_z,50 !洞的进深 /prep7et,1,mesh200,2 !用于3-D的2节点线et,2,mesh200,6 !用于3-D的4节点四边形et,3,shell63et,4,solid45r,1,th !壳的厚度mp,ex,1,3.0e10 !壳的材料,C30混凝土mp,prxy,1,0.2 mp,dens,1,2500 mp,ex,2,4.5e8 !保留岩石的材料mp,prxy,2,0.32 !tb,dp,2!tbdata,1,20,30, mp,dens,2,2700mp,ex,3,4.51e8 !挖去岩石的材料 mp,prxy,3,0.32 !tb,dp,3!tbdata,1,20,30, mp,dens,3,2700k,,0,0k,,0,3.85k,,0.88,5.5k,,2.45,6.15k,,4.02,5.5k,,4.9,3.85k,,4.9,0larc,1,2,6,8.13 !定义两点之间的圆弧线，larc,p1,p2,pc,radlarc,2,3,6,3.21 larc,3,4,6,2.22 larc,4,5,2,2.22 larc,5,6,2,3.21 larc,6,7,2,8.13 larc,7,1,4,6a,1,2,3,4,5,6,7 !产生面1!blc4,x1,y1,w,h1 !block4,xcorner,ycorner,width,height,depth!blc4,x2,y2,w,h2 !block4,xcorner,ycorner,width,height,depth 产生面3 blc4,x1,y1,w1,h1 !产生面2 blc4,x2,y2,w2,h2 !产生面3blc4,x3,y3,w3,h3 !产生面4 blc4,x4,y4,w4,h4 !产生面5 /pnum,area,1 aovl,1,2,3,4,5 !布尔操作重叠，得到面3 nummrg,all,,,,low numcmp,all l,1,8 !从四个角点上连接出四条线 l,7,9l,6,10l,2,11lsel,s,line,,21,22,1 !用线分割面 lsel,a,line,,7asbl,5,alllsel,s,line,,21,24,3 lsel,a,line,,1asbl,7,alllsel,s,line,,22,23,1 lsel,a,line,,6asbl,8,allnummrg,all,,,,low numcmp,alllsel,s,line,,2,5,1 LCCAT,alllesize,all,,,3lsel,s,line,,9,11,2 lsel,a,line,,6lsel,a,line,,1lesize,all,,,8lsel,s,line,,8,10,2 lsel,a,line,,7lesize,all,,,12lsel,s,line,,21,24,1 lesize,all,,,10,2type,2asel,s,area,,5,8,1 amesh,all!mopt,split,on !指定网格划分选项，四边形单元如果违背了单元形状则将其分解为三角形单元!mopt,split,err !mopt,split,warn !mopt,qmesh,main !选用四边形网格划分器，当这种操作失败时，不调用备用的网格划分器asel,s,area,,1amesh,1lsel,s,line,,12,13,1 lesize,all,,,8lsel,s,line,,15,18,1 lesize,all,,,6,2 asel,s,area,,2,3,1 amesh,all lsel,s,line,,14 lesize,all,,,24 lsel,s,line,,19,20,1lesize,all,,,6,2 lsel,s,line,,15,17,2 lsel,a,line,,8LCCAT,allasel,s,area,,4amesh,allLSEL,s,LCCALDELE,allnummrg,all,,,,low numcmp,allallsel!以下开始拉伸成实体单元!首先拉伸成壳单元k,1000,,,-length_z l,1,1000/view,1,1,1,1/replotEXTOPT,ESIZE,10,0, LSEL,S,LINE,,1,7,1 ADRAG,all,,,,,,25 gplot type,3real,1mat,1ASEL,S,loc,z,-25 APLOTlsel,s,loc,z,-25 lesize,all,,,10 MSHAPE,0,2DMSHKEY,1amesh,all!拉伸岩石的实体ASEL,invertaplotEXTOPT,ESIZE,10,0, EXTOPT,ACLEAR,1 TYPE,4MAT,2asel,r,area,,2,8,1 VDRAG,all,,,,,,25 allsel!挖去部分岩石的实体MAT,3VDRAG,1,,,,,,25 EPLOTnummrg,all,,,,low numcmp,all!约束两侧面的X方向的约束asel,s,loc,x,x2 asel,a,loc,x,x2+w4 da,all,ux,0alls!约束地面的Y方向的约束asel,s,loc,y,y4 da,all,uy,0allsasel,s,loc,z,-length_zasel,a,loc,z,0 da,all,uz,0allselacel,,10fini/soluantype,staticdeltim,0.1,0.05,0.2 autots,on !使用自动时间步 pred,on !打开时间步长预测器 lnsrch,on !打开线性搜索nlgeom,on !打开大位移效果nropt,full !设定牛顿-拉普森选项cnvtol,f,,0.02,2,0.5esel,s,type,,3 !选择梁单元，杀死ekill,allesel,allesel,s,live !选择活的单元nsle,s !选择活单元上的节点nsel,invert !反向选择，即选择了死单元上的节点 d,all,all,0 !将死单元上的节点约束所有位移，使其不参与矩阵运算nsel,allesel,all/PBC,ALL,,1gplotsolve*do,ii,1,10,1esel,s,mat,,3 !选择挖去的岩石，杀死nsle,snsel,r,loc,z,0.1-(ii-1)*5,-(5.1+(ii-1)*5)esln,r,1ekill,allesel,s,type,,3 !激活挖去的岩石对应的壳单元,并将其节点上的约束删除nsle,snsel,r,loc,z,0.1-(ii-1)*5,-(5.1+(ii-1)*5)esln,r,1ealive,allnsle,sddele,all,allesel,allesel,s,live !选择活单元，此时应该包含两部份，一是梁单元，二是未挖去的岩石单元nsle,snsel,invert !反向选择，将死单元上的节点约束所有自由度d,all,all,0nsel,allesel,allsolve*enddo/post1/DEVICE,VECTOR,1 ESEL,S,type,,3 SET,1,LAST,1, PLNSOL,U,Y,0,1SET,2,LAST,1, PLNSOL,U,Y,0,1 SET,6,LAST,1, PLNSOL,U,Y,0,1 SET,11,LAST,1, PLNSOL,U,Y,0,1!*dim ,VUY_1!*dim ,VUY_2!*dim ,VUY_6!*dim ,VUY_11!SET,1,LAST,1, !PLNSOL,U,Y,0,1 !PLNSOL,U,SUM,0,1 !ETABLE,EUY_1,U,Y ! PLETAB,EUY_1,AVG !PRETAB,EUY_1!*vget,VUY_1,node,387,U,Y!SET,2,LAST,1, !PLNSOL,U,Y,0,1 !PLNSOL,U,SUM,0,1 !ETABLE,EUY_2,U,Y ! PLETAB,EUY_2,AVG !PRETAB,EUY_2!*vget,VUY_2,node,387,U,Y!SET,6,LAST,1, !PLNSOL,U,Y,0,1 !PLNSOL,U,SUM,0,1 !ETABLE,EUY_6,U,Y ! PLETAB,EUY_6,AVG !PRETAB,EUY_6!*vget,VUY_6,node,323,U,Y!SET,3,LAST,1,!PLNSOL,U,Y,0,1!PLESOL,SMISC,4 !MX !PLESOL,SMISC,5 !MY !PLESOL,SMISC,6 !MXY!ETABLE,p_i,SMISC, 13 !ETABLE,p_j,SMISC, 14 !ETABLE,p_k,SMISC,15 !ETABLE,p_l,SMISC, 16!ETABLE,p_i_3,SMISC, 17 !ETABLE,p_j_3,SMISC, 18!ETABLE,p_k_4,SMISC, 20 !ETABLE,p_l_4,SMISC, 19!ETABLE,p_k_5,SMISC, 21 !ETABLE,p_l_5,SMISC, 22 !fini!/post1!etable,if,smisc,1 !etable,jf,smisc,7 !etable,im,smisc,6 !etable,jm, smisc,12大家可以在这个的基础上进一步发掘。

ANSYS二次开发技术及其在土木工程中的应用_林能辉

2
基于 USERMAT 的用户滞回特性 ANSYS 的二次开发
ANSYS 二次开发平台
UEL 和 USERMAT 都是基于 FORTRAN 语言的源程序，要
2． 1
使用 UEL 和 USERMAT 对 ANSYS 进行二次开发，用户首先需要一个编译平台，本文使用的平台组合为：（ 1 ） Microsoft Visual Studio 2005 Professional Edition （ 2 ） Intel Visual Fortran Compiler Professional v11． 1 （ 3 ） ANSYS 12． 0 调整环境参数并修改 Intel Fortran Compiler 和 ANSYS 自带的批处理文件调用语句，使得开发平台实现整合，同时修改相应的环境参数。
3
（ Beijing Key Lab of Earthquake Engineering and Structural Retrofit， Beijing University of Technology， Beijing 100124 ， China）（ Beijing CCI Engineerring Design and Consulting Co．，Ltd．，Beijing 100044 ， China）
Abstract
A numerical simulation technology based on secondarydevelopment function of ANSYS is presented herein． It can be used to
simulate userdefined nonlinear stressstrain relationship． The operating procedure of development approach for userdefined hysteretic model based on ANSYS is given． The hysteretic characteristics of EDR damper and pseudolinear hysteretic curve damper with frictionvariable characteristics are realized using this kind of technology． The negativestiffness characteristics of damper elements are also realised based on secondarydevelopment function of ANSYS． The validity of this kind of technology is verified in numerical analysis． The very well practicability would be expected． Keywords Secondarydevelopment of ANSYS Numerical simulation Negativestiffness characteristics EDR damper Pseudolinear

ANSYS在土石坝防渗加固分析中的应用

ANSYS在土石坝防渗加固分析中的应用土石坝是一种常见的水利工程结构，用于蓄水或防洪目的。

土石坝防渗加固是土石坝设计中非常重要的一环，其主要目的是防止坝体内的水从坝体内部泄漏出去，影响坝体的稳定性和使用寿命。

在土石坝防渗加固设计中，ANSYS软件可以发挥重要的作用。

1.坝体内水压分析：土石坝在蓄水或雨水侵蚀的情况下，坝体内部会产生水压力。

通过ANSYS软件可以对坝体内部的水压力进行分析，包括水压的大小、分布情况等。

这些信息对于设计合适的防渗加固措施非常重要。

2.渗流路径分析：土石坝防渗加固的设计需要准确地预测坝体内水的渗流路径。

通过ANSYS软件可以建立数值模型，模拟坝体内水通过不同部位的渗流路径，找出可能的渗漏点和路径，为加固设计提供参考。

3.加固材料力学性能分析：土石坝防渗加固通常涉及到各种加固材料，如混凝土、聚合物和土工合成材料等。

通过ANSYS软件可以对这些材料的力学性能进行分析，包括强度、刚度、变形等，为加固设计提供依据。

4.加固结构设计优化：在土石坝防渗加固设计过程中，需要确定合适的加固结构形式和参数。

通过ANSYS软件可以建立不同加固结构的有限元模型，进行受力分析和设计优化，找出最优的加固方案。

5.加固效果评估：在加固完成后，需要对加固效果进行评估。

通过ANSYS软件可以对加固后的土石坝进行受力分析和模拟，评估加固效果，判断是否满足设计要求。

综上所述，ANSYS在土石坝防渗加固设计中的应用具有重要意义。

通过ANSYS软件的分析和模拟，可以帮助工程师更准确地了解土石坝的受力和渗流情况，设计出更科学、安全、经济的防渗加固方案。

在实际工程应用中，工程师们可以借助ANSYS软件的强大功能，提高土石坝防渗加固设计的效率和可靠性，确保土石坝的安全运行和长期稳定性。

ansys土木工程应用实例

第三章回弹/units,si/prep7et,1,plane42keyopt,1,3,2mp,ex,1,220e9mp,nuxy,1,0.3TB,MISO,1,1,25, TBTEMP,0TBPT,,0.000814,175.88e6, TBPT,,0.001531,324.98e6, TBPT,,0.002819,591.26e6, TBPT,,0.003846,801.457e6, TBPT,,0.004239,873.22e6, TBPT,,0.004594,924.399e6, TBPT,,0.004955,956.88e6, TBPT,,0.00531,976.99e6, TBPT,,0.005699,990.2138e6, TBPT,,0.006261,1000e6, TBPT,,0.006957,1008.34e6, TBPT,,0.007804,1015.5e6, TBPT,,0.008779,1021.185e6, TBPT,,0.009827,1026.47e6, TBPT,,0.010929,1031.289e6, TBPT,,0.012066,1035.444e6, TBPT,,0.0132,1039.693e6, TBPT,,0.014354,1044.886e6, TBPT,,0.015494,1049.796e6, TBPT,,0.016631,1053.951e6, TBPT,,0.017785,1058.106e6, TBPT,,0.019039,1062.827e6, TBPT,,0.020276,1067.643e6, csys,1k,1,0.00512k,2,0.0064kgen,36,1,2,1,,10*do,j,1,69,2a,j,j+1,j+3,j+2*enddoa,71,72,2,1savelesize,all,,,4amesh,allnummrg,allnumcmp,allsavecp,1,uy,166,170,171,172,173finish/solunsel,s,loc,y,270d,all,allnplotallselnplotnsel,s,loc,y,90d,all,uxnplotallseloutres,all,alltime,1NSUBST,100nlgeom,ond,166,uy,-0.00128solve*get,def1,node,166,U,Ynsel,s,loc,y,90ddele,all,uyallself,166,fy,0solve*get,def2,node,166,U,Yfinish/post1 !进入后处理器set,list !显示所有计算的结果set,1,last,1 !读取第一个加载步最后子步的结果pldisp,1 !同时显示变形前与变形后结构的形状。

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

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