土木工程专业ANSYS初学者经典例题
上机实例练习ANSYS.
《弹性力学有限元》上机作业一、一带孔条形板形状、尺寸如图所示,l=300mm,w=200mm,r=20mm,h=30mm,t=5mm,板两端受均布拉力作用,q=100MPa。
已知材料的弹性模量和泊松比分别为E=210GPa,μ=0.3,求板有孔截面上的理论应力集中系数k。
报告要求:1.简述建模过程;2.画出应力分布规律图;3.求出应力集中系数;4.思考:变形最大的地方是不是应力也最大?提示:1.平面应力问题,取1/4结构进行建模,采用N、mm、MPa单位制。
2.创建矩形:(0,0)(150,100),Main Menu-Preprocessor-Modeling-Create-Areas-Rectangle-3.创建矩形:(0,0)(20,30)4.创建直角扇形:(0,30)(R=0-20,θ=0-90°)-Modeling-Create-Areas-Circle-Partial Annulus5.从大矩形减法小矩形和扇形:Main Menu-Preprocessor-Operate-Booleans-Subtract6. 创建单元类型:-Preprocessor-Element Type-Add/Edit/Delete-Solid-Plane42\82\182\183, Option-Plane Stress with thickness input7. 建立实常数:-Preprocessor-Real constant- Add/Edit/Delete-给单元厚度5mm8. 定义材料模型及其属性:-Preprocessor-Material- Material Models-左窗口Material Model Number 1, 右窗口Structural-Linear-Elastic-Isotropic给出弹性模量和泊松比及密度。
9. 网格剖分:Main Menu-Preprocessor-Meshing Tools-, 注意随时保存模型。
ansys工程实例(经典例子)
管道支架结构分析一问题描述该结构用于支撑管道,如图所示。
该结构需要有很好的长时间的支撑性,且在支撑时,变形不能过大,否则会由于支撑力不够,造成管道变形,严重的话会造成管道的泄露。
另外,所用的材料也要满足屈服条件,设计时不能造成结构的破坏。
如何设计该支撑的结构和所用的材料成了其中的关键。
材料参数为7E+008,泊松比为0.33,边界条件为最下端为固定端,载荷为管道所在弧面上,方向为垂直且指向弧面的均布面力。
二求解步骤定义工作文件名Utility Menu-->File-->Change Jobname 该工作名为yangxin10054554定义单元类型Main Menu --> Preprocessor--> Element Type --> Add/Edit/Delete…创建mesh200和brick 20node 95单元。
(mesh200还需设置options选择面单元,否则分网时会提示出问题)材料参数设定main menu-->preferences-->…选中结构类选项。
Main menu-->preprocessor-->material props-->material models-->在material models available 分组框中依次选取structural/linear/elastic/isotropic选项,设置弹性模量EX=0.7e9,泊松比=0.33。
4.生成几何模型、划分网格Main menu-->preprocessor-->modeling-->create-->keypoints-->in active cs 选项,输入关键点号和相应的坐标,如下:2)连线Main menu-->preprocessor-->modeling-->create-->lines-->lines-->straightline-->…3) 倒角Main menu-->preprocessor-->modeling-->create-->lines-->line fillet-->...4)对称Main menu-->preprocessor-->modeling-->reflect-->lines-->…之后将所有面add在一起。
ANSYS分析题目
一、工程概况:
某隧道工程,隧洞位于地下15米处,洞直径10米,直墙高10米,拱高5米,其具体尺寸见图1。
根据工程地质勘察报告,岩土各参数为:ρ= 2500kg/m3,Ε =1.4×109Ρa,v=0.27,c = 2.72×106Ρa, = 35°。
地面上将修建一30层的高层建筑,根据估计每层每米有2.5吨的荷载直接作用于地基上。
二、计算要求:
1. 分两种情况进行分析
(1)先开挖隧洞,在修建地面建筑,分析全过程隧洞洞顶下沉位移、塑性区、拉坏区和压坏区;
(2)地面上先修建好以后,再开挖隧道,分析隧洞洞顶下沉位移、塑性区、拉坏区和压坏区;
2. 本分析考虑整个施工过程的影响,为更好地体现施工过程对隧道的影响,在分析中加荷分为三次进行,每次加荷10层(即250kN/m);
3. 确定边界范围;
4. 用有限元分析软件ANSYS建模、求解、分析。
ansys上机考试题及参考答案
ansys上机考试题及参考答案一、选择题(每题2分,共20分)1. 在ANSYS中,以下哪个选项不是材料属性的参数?A. 弹性模量B. 泊松比C. 密度D. 摩擦系数答案:D2. ANSYS中,以下哪个命令用于定义几何体?A. /PREP7B. /POST1C. /SOLUD. /FINI答案:A3. 在ANSYS中,以下哪个选项不是施加载荷的方式?A. 压力B. 温度C. 重力D. 惯性答案:D4. ANSYS中,以下哪个选项不是网格划分的控制参数?A. 网格大小B. 网格形状C. 网格密度D. 材料属性答案:D5. 在ANSYS中,以下哪个选项不是边界条件?A. 固定支撑B. 压力载荷C. 位移约束D. 温度变化答案:B6. ANSYS中,以下哪个命令用于求解?A. /SOLUB. /PREP7C. /POST1D. /FINI答案:A7. 在ANSYS中,以下哪个选项不是后处理的功能?A. 查看应力分布B. 查看变形图C. 计算体积D. 定义材料属性答案:D8. ANSYS中,以下哪个命令用于退出求解器?A. /SOLUB. /FINIC. /PREP7D. /POST1答案:B9. 在ANSYS中,以下哪个选项不是网格划分的类型?A. 自由网格B. 映射网格C. 扫掠网格D. 材料属性答案:D10. ANSYS中,以下哪个命令用于退出前处理器?A. /PREP7B. /FINIC. /SOLUD. /POST1答案:B二、填空题(每空1分,共20分)1. 在ANSYS中,定义材料属性的命令是 _ 。
答案:MP2. 施加固定支撑的边界条件的命令是 _ 。
答案:D3. ANSYS中,查看应力分布的后处理命令是 _ 。
答案:PLNSOL, S4. 定义几何体的坐标系的命令是 _ 。
答案:CSYS5. 在ANSYS中,施加重力载荷的命令是 _ 。
答案:CE6. 退出后处理器并返回到前处理器的命令是 _ 。
ansys上机练习
第一题:如图等腰直角三角形薄板,直角边AB 、BC 长度为1m ,板厚0.01m ,板上有一半径为0.1m 的通孔。
薄板BC 边固定,AC 边施加压强为1000Pa 的载荷,方向与AC 边垂直。
薄板的弹性模量为200GPa ,泊松比为0.3。
求薄板的变形形状、应力分布。
第二题:如图所示空心平板,一端固定,另一端受两个集中力,平板为铝材,厚度1cm ,弹性模量为80GPa ,泊松比为0.33。
选用PLANE82单元,单元尺寸设为5mm ,求平板受力后的变形和应力分布。
A第三题:如图所示空心平板,一端固定,另一端受一个0~10Pa 分布力,平板为铝材,厚度2mm ,弹性模量为80GPa ,泊松比为0.33。
选用PLANE82单元,单元尺寸设为5mm ,求平板受力后的变形和应力分布。
第四题:一个带三个圆孔的平板,板厚1mm ,弹性模量200GPa ,泊松比为0.25。
模型左侧的两个小圆孔上施加所有位移约束,右侧大圆孔边界最下面中心施加载荷F=1000N ,要求采用平面8节点四边形单元,单元尺寸2mm ,求平板变形后应力分布。
(图中单位:mm )第五题:如图平面支架,左端固定,上端受10kN/m 的分布力,选用PLANE82单元,单元尺寸设为0.15m ,求平面支架受力后的变形和应力分布。
平面支架的弹性模量为210GPa ,泊松比为0.3。
第六题:悬臂梁长2m ,截面形状如图,梁末端作用有20kN 的集中力,梁的弹性模量为200GPa ,泊松比为0.3。
梁长度方向划分为10个单元,求解梁的挠度及应力分布。
10kN/m第七题:简支梁长2m ,截面形状如图,梁作用有5000Pa 的分布力,梁的弹性模量为200GPa ,泊松比为0.3。
梁长度方向划分为10个单元,求解梁的挠度及应力分布。
第八题:悬臂梁长2m ,梁左侧一半作用有1N/m 的分布力,梁的弹性模量为210GPa ,泊松比为0.3,横截面是0.04×0.04m2的正方形截面。
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。
ANSYS案例——20例ANSYS经典实例】
→ OK ( 返回到 Plot Results)
(10) 退出系统
ANSYS Utility Menu: File→ Exit … → Save Everything→ OK
工作文件名 ) : beam3→ Run → OK
(2) 设置计算类型 ANSYS Main Menu: Preferences …→ Structural → OK
(3) 选择单元类型 ANSYS Main Menu: Preprocessor → Element Type → Add/Edit/Delete
(7) 模型施加约束和外载 左边加 X 方向的受力 ANSYS Main Menu: Solution → Define Loads → Apply → Structural → Force/Moment → On
Nodes → 选择节点 1→ apply →Direction of force: FX → VALUE : 3000 → OK →
上,完成相应的力学分析。
3-20 。在 ANSYS 平台
图 3-19 框架结构受一均布力作用
( a) 节点位移及单元编号
(b) 等效在节点上的外力
图 3-20 单元划分、节点位移及节点上的外载
形界面的交互式操作 (step by step)
(1) 进入 ANSYS( 设定工作目录和工作文件 ) 程序 → ANSYS → ANSYS Interactive →Working directory (设置工作目录 ) → Initial jobname ( 设置
ANSYS岩土计算例子
ANSYS土工结构计算案例ANSYS-CHINA广州办事处2020年7月14日目录计算题目及计算要求说明 (1)题目一 (4)一、计算说明 (4)二、计算所用ANSYS邓肯-张的E-B模型说明 (5)三、计算有限元模型及计算结果 (6)题目二 (7)一、用三维有限元模型计算 (7)二、用三维有限元模型计算 (8)题目三 (10)一、计算说明 (10)二、计算有限元模型及边界条件 (10)三、强夯地基固结计算 (10)题目四 (17)一、计算说明 (17)二、计算几何模型和有限元模型 (17)三、计算结果 (18)1、计算边界条件 (18)2、计算结果 (19)3、结论 (20)计算题目及计算要求说明题目一:高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析题目二:大圆筒结构、波浪与地基的相互作用分析(大圆筒作为重力式码头结构,波浪为水平动荷载,门吊为竖向动荷载,地基为三层以上地基包括自抛碎石垫层、粘土层、粉细砂层和岩层,粉细砂层可能在波浪动荷载作用下液化造成圆筒倾覆)题目三:(冲击荷载下)强夯地基固结有限元分析(提供固结方程或固结方程处理方案,孔隙水压力消散计算方案、沉降计算方案及其他一些处理技巧)题目四:在降雨情况下土工格栅加筋土挡墙边坡上公路稳定分析(由上至下为公路面层,垫层,挡墙,挡墙面板采用预制混凝土块0.60.60.6m3,混凝土后方为钩挂式土工格栅,边坡比较陡,边坡有一定排水特性)。
具体处理方案包括:1、提供计算输入界面2、计算模型或采用本构情况3、前处理方案及网格划分技巧4、特殊材料或模型嵌入技术5、计算技巧及解决方案6、后处理提供内容具体报价方案包括:进行简单报价,涉及以上题目的各模块的综合报价(包括前后处理及解决以上问题的结构与岩土问题的模块报价. 能用通用模块计算尽量用通用模块,必须用CivilFEM模块计算的请注明.提供解决方案时间:2005年5月18日之前。
题目一高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析一、计算说明高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析,可以用空间有限元模型或平面有限元模型数值求解。
ansys土木经典算例
ansys土木经典算例ANSYS是一种常用的有限元软件,广泛应用于土木工程领域。
下面列举了一些经典的ANSYS土木算例。
1. 悬臂梁的静力分析:通过给定的悬臂梁的几何尺寸和加载条件,计算梁的应力、位移和变形情况,以评估梁的承载能力和结构安全性。
2. 桁架结构的模态分析:对给定的桁架结构进行模态分析,确定结构的固有频率和振型,以预测结构的动力响应和抗震性能。
3. 地基基础的强度分析:通过模拟地基基础的加载情况,计算地基的应力分布和变形情况,以评估地基的稳定性和承载能力。
4. 混凝土梁的弯曲分析:对给定的混凝土梁进行弯曲分析,确定梁的应力分布和变形情况,以评估梁的承载能力和结构安全性。
5. 土体的固结分析:通过模拟土体的加载和固结过程,计算土体的应力、应变和固结沉降情况,以评估土体的稳定性和变形特性。
6. 钢结构的塑性极限分析:对给定的钢结构进行塑性极限分析,确定结构的塑性极限强度和塑性变形分布,以评估结构的安全性和可靠性。
7. 水坝的渗流分析:通过模拟水坝的渗流过程,计算水坝的渗流速度、渗流压力和渗流路径,以评估水坝的渗漏风险和抗渗性能。
8. 地下管道的地震响应分析:对给定的地下管道进行地震响应分析,确定管道的动力响应和变形情况,以评估管道的抗震性能和结构安全性。
9. 地铁隧道的地表沉降分析:通过模拟地铁隧道的开挖和施工过程,计算地表沉降和沉降分布,以评估隧道施工对周围建筑物和地下管线的影响。
10. 高层建筑的风载分析:对给定的高层建筑进行风载分析,确定建筑物的风压分布和结构的响应,以评估建筑物的抗风性能和结构安全性。
以上是一些经典的ANSYS土木算例,涵盖了土木工程中的静力分析、动力分析、地基基础分析、固结分析、渗流分析、地震响应分析、地表沉降分析和风载分析等多个方面,可以帮助工程师们更好地理解和评估土木结构的性能和安全性。
ansys练习题及参考答案
ANSYS 上机试题题目1:如图,带倒角的长方形平板有三个销钉孔,尺寸如图所示(单位mm),忽略板的厚度。
材料E=2E11Pa,μ= 0.3。
在板的三个销钉孔施加全约束,板的上表面受到均匀压力q=1000Pa。
求,该板的总位移、SEQV当量应力。
请按照以下要求完成考核内容:1.按照所给尺寸正确地将几何模型建立出来;2.对几何模型进行规则的四边形网格划分;3.加边界条件及载荷;4.求解,并将模型的总位移和总应力强度的图片粘在word文档里,题目注明班级、姓名及学号;5.在45min内完成作业。
提示:单元类型:plane42,建立2-D模型;注意保存数据文件。
q=1000Pa题目1的参考答案:(1)命令流fini/clear/PREP7ET,1,42 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2e11 MPDATA,PRXY,1,,0.3 RECTNG,,0.14,,0.1, LFILLT,4,3,0.02, , LFILLT,3,2,0.02, , LFILLT,2,1,0.02, , LFILLT,1,4,0.02, ,FLST,3,4,4,ORDE,4 FITEM,3,5FITEM,3,8FITEM,3,11FITEM,3,14ASBL, 1,P51X FLST,2,4,5,ORDE,2 FITEM,2,2FITEM,2,-5ADELE,P51X, , ,1 wpoff,0.02,0.02CYL4, , ,0.01FLST,3,1,5,ORDE,1 FITEM,3,1AGEN,2,P51X, , , ,0.06, , ,0 FLST,3,1,5,ORDE,1 FITEM,3,2AGEN,2,P51X, , ,0.1, , , ,0 FLST,3,3,5,ORDE,2 FITEM,3,1FITEM,3,-3ASBA, 6,P51XESIZE,0,0.1,MSHAPE,0,2DMSHKEY,0AMESH,all/SOLFLST,2,12,4,ORDE,8FITEM,2,6FITEM,2,-7FITEM,2,9FITEM,2,-10FITEM,2,12FITEM,2,-13FITEM,2,15FITEM,2,-20DL,P51X, ,ALL,FLST,2,1,4,ORDE,1FITEM,2,3SFL,P51X,PRES,1000,SOLVE(2)考察结果几何模型:有限元模型:加载及约束:位移响应:应力响应:题目2:如下图所示,一个承受单向拉伸的大平薄板(可忽略厚度方向尺寸),在其中心位置有一个小圆孔。
Ansys例题2
第一题
建模思路
1.在空白区域生成block,注意这里的block是最终底座的一半大,我们先完成一半的建模,之后镜像即可。
2.在画好的block上面画一个cylinder圆柱体,通过平移复制出另一个圆柱,绘制完毕利用布尔相减指令得到孔。
3.在孔绘制完毕之后在上平面一次画一个block与一个四分之一cylinder,注意在此处长度尺均为最终结果的一半,因为我们需要先画出来一半即可。
4.下一步我们要打孔,思路与打底板的孔类似,先画出两个cylinder圆柱出来,然后利用布尔相减的方法打孔,在这里要注意因为实体较多,选择的时候注意不要选错,否则打孔会失败,效果如下。
5.画到这一步,我们还差肋板没有画,利用现有的几个KeyPoint可以相对容易的画出肋板的一个面(如下图),画完之后拉伸即得到肋板。
6.到这里,基本上大功告成,最后利用Reflect Volumes做出另一半的镜像体即可。
ANSYS案例——20例ANSYS经典实例】
三梁平面框架结构的有限元分析针对【典型例题】3.3.7(1)的模型,即如图3-19所示的框架结构,其顶端受均布力作用,用有限元方法分析该结构的位移。
结构中各个截面的参数都为:113.010Pa E =⨯,746.510m I -=⨯,426.810m A -=⨯,相应的有限元分析模型见图3-20。
在ANSYS 平台上,完成相应的力学分析。
图3-19 框架结构受一均布力作用(a ) 节点位移及单元编号 (b ) 等效在节点上的外力图3-20 单元划分、节点位移及节点上的外载解答 对该问题进行有限元分析的过程如下。
1.基于图形界面的交互式操作(step by step)(1) 进入ANSYS(设定工作目录和工作文件)程序 →ANSYS → ANSYS Interactive →Working directory (设置工作目录) →Initial jobname (设置工作文件名): beam3→Run → OK(2) 设置计算类型ANSYS Main Menu: Preferences… → Structural → OK(3) 选择单元类型ANSYS Main Menu: Preprocessor →Element Type →Add/Edit/Delete… →Add… →beam :2D elastic 3 →OK (返回到Element Types 窗口) →Close(4) 定义材料参数ANSYS Main Menu:Preprocessor →Material Props →Material Models→Structural →Linear →Elastic→Isotropic: EX:3e11 (弹性模量) →OK →鼠标点击该窗口右上角的“ ”来关闭该窗口(5) 定义实常数以确定平面问题的厚度ANSYS Main Menu: Preprocessor →Real Constant s… →Add/Edit/Delete →Add →Type 1 Beam3→OK→Real Constant Set No: 1 (第1号实常数), Cross-sectional area:6.8e-4 (梁的横截面积) →OK →Close(6) 生成几何模型生成节点ANSYS Main Menu: Preprocessor →Modeling →Creat→Nodes→In Active CS→Node number 1 →X:0,Y:0.96,Z:0 →Apply→Node number 2 →X:1.44,Y:0.96,Z:0 →Apply→Node number 3 →X:0,Y:0,Z:0→Apply→Node number 4 →X:1.44,Y:0,Z:0→OK生成单元ANSYS Main Menu: Preprocessor →Modeling →Create →Element →Auto Numbered →Thru Nodes →选择节点1,2(生成单元1)→apply →选择节点1,3(生成单元2)→apply →选择节点2,4(生成单元3)→OK(7)模型施加约束和外载左边加X方向的受力ANSYS Main Menu:Solution →Define Loads →Apply →Structural →Force/Moment →On Nodes →选择节点1→apply →Direction of force: FX →V ALUE:3000 →OK→上方施加Y方向的均布载荷ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Beams →选取单元1(节点1和节点2之间)→apply →V ALI:4167→V ALJ:4167→OK左、右下角节点加约束ANSYS Main Menu:Solution →Define Loads →Apply →Structural →Displacement →On Nodes →选取节点3和节点4 →Apply →Lab:ALL DOF →OK(8) 分析计算ANSYS Main Menu:Solution →Solve →Current LS →OK →Should the Solve Command be Executed? Y→Close (Solution is done! ) →关闭文字窗口(9) 结果显示ANSYS Main Menu: General Postproc →Plot Results →Deformed Shape … →Def + Undeformed →OK (返回到Plot Results)(10) 退出系统ANSYS Utility Menu: File→Exit …→Save Everything→OK(11) 计算结果的验证与MA TLAB支反力计算结果一致。
ANSYS练习题
ANSYS练习题习题1:矩形截面简支梁,跨度18米,高3米,厚度0.4米,弹性模量E=2×106吨/米2,泊松比μ=0.167。
在梁的上顶面承受均匀载荷q=10吨/米2。
分别用矩形单元、六节点三角形单元和八节点等参单元计算,求轴对称截面(x=0)上的x方向正应力(σx)分布。
1)四节点矩形单元:单元划分与施加的约束和荷载2)六节点三角形单元:单元划分与施加的约束和荷载单元划分与施加的约束和荷载4)四节点矩形单元加密2倍:单元划分与施加的约束和荷载单元划分与施加的约束和荷载6)结果分析根据矩形截面简支梁在均布荷载作用下应力分量的弹性力学解答,可以求得轴对称截面(x=0)上下边缘x方向的正应力精确解为 2.72MPa。
从矩形单元、六节点三角形单元和八节点等参单元的计算结果可见,结构采用八节点等参数单元在较少的单元数量下就能够达到较高的计算精度,而采用四节点矩形单元的计算精度相对较低。
采用矩形单元加密分析,可以得出较为精确的结果。
习题2:对于如图所示的平面圆孔问题,通过数值方法(分别采用三角形单元和四边形单元)研究不同板宽的孔边应力集中问题,与弹性力学的解析解进行比较。
分析应力集中系数与相对孔径尺度的关系,应力集中影响区的范围。
板宽与孔径2R的相对尺度自己定义。
固定B=0.6,长度L=1,2R/B=1/20、1/10、1/5、1/3,薄板弹性模量为69GPa,泊松比为0.3,两端的均布荷载为0.5MPa。
其中2R/B=1/20、1/5的采用三角形单元分析,其余采用四边形单元。
1)2R/B=1/202)2R/B=1/103)2R/B=1/54)2R/B=1/3根据数值计算的结果,圆孔边沿的应力一般是施加的平均压力的3倍左右。
2R/B越大,其应力集中越明显,应力集中的影响范围越大。
这一结果与弹性力学的理论解一致。
ANSYS基础培训练习1
ANSYS基础培训练习1轴承座轴瓦轴四个安装孔径向约束 (对称) 轴承座底部约束 (UY=0) 沉孔上的推力(1000 psi.)向下作用力 (5000 psi.) ANSYS 基础培训练习题第一日练习主题:实体建模EX1:轴承座的实体建模、网格划分、加载、求解及后处理练习目的:创建实体的方法,工作平面的平移及旋转,布尔运算(相减、粘接、搭接,模型体素的合并,基本网格划分。
基本加载、求解及后处理。
问题描述:具体步骤:首先进入前处理(/PREP7) 1. 创建基座模型生成长方体Main Menu :Preprocessor>Create>Block>By Dimensions 输入x1=0,x2=3,y1=0,y2=1,z1=0,z2=3 平移并旋转工作平面Utility Menu>WorkPlane>Offset WP by Increments X,Y ,Z Offsets 输入2.25,1.25,.75 点击Apply XY ,YZ ,ZX Angles 输入0,-90点击OK 。
创建圆柱体Main Menu :Preprocessor>Create>Cylinder> Solid Cylinder Radius输入0.75/2, Depth 输入-1.5,点击OK 。
载荷拷贝生成另一个圆柱体Main Menu:Preprocessor>Copy>V olume拾取圆柱体,点击Apply, DZ输入1.5然后点击OK 从长方体中减去两个圆柱体Main Menu:Preprocessor>Operate>Subtract V olumes首先拾取被减的长方体,点击Apply,然后拾取减去的两个圆柱体,点击OK。
使工作平面与总体笛卡尔坐标系一致Utility Menu>WorkPlane>Align WP with> Global Cartesian2. 创建支撑部分Utility Menu: WorkPlane -> Display Working Plane (toggle on)Main Menu: Preprocessor -> -Modeling-Create -> -Volumes-Block -> By 2 corners & Z在创建实体块的参数表中输入下列数值:WP X = 0WP Y = 1Width = 1.5Height = 1.75Depth = 0.75OKToolbar: SA VE_DB3. 偏移工作平面到轴瓦支架的前表面Utility Menu: WorkPlane -> Offset WP to -> Keypoints +1. 在刚刚创建的实体块的左上角拾取关键点2. OKToolbar: SAVE_DB4.创建轴瓦支架的上部Main Menu: Preprocessor -> Modeling-Create -> Volumes-Cylinder -> Partial Cylinder + 1). 在创建圆柱的参数表中输入下列参数:WP X = 0WP Y = 0Rad-1 = 0Theta-1 = 0Rad-2 = 1.5Theta-2 = 90Depth = -0.752). OKToolbar: SAVE_DB5. 在轴承孔的位置创建圆柱体为布尔操作生成轴孔做准备Main Menu: Preprocessor -> Modeling-Create -> V olume-Cylinder -> Solid Cylinder +1.) 输入下列参数:WP X = 0WP Y = 0Radius = 1Depth = -0.18752.) 拾取Apply3.) 输入下列参数:WP X = 0WP Y = 0Radius = 0.85Depth = -24.)拾取OK6.从轴瓦支架“减”去圆柱体形成轴孔.Main Menu: Preprocessor -> Modeling-Operate -> Subtract -> Volumes +1. 拾取构成轴瓦支架的两个体,作为布尔“减”操作的母体。
ANSYS应用精华30例
第五例杆系结构的静力学分析实例——平面桁架/clearL=0.1A=1e-4/PREP7ET,1,BEAM3R,1,AMP,EX,1,2E11MP,PRXY,1,0.3N,1N,2,LN,3,2*LN,4,L,LE,1,2E,2,3E,1,4E,2,4E,3,4FINISH/SOLUD,1,UXD,1,UYD,3,UYF,4,FY,-2000SOLVEFINISH/POST1ETABLE,FA,SMISC,1ETABLE,SA,LS,1PRETAB,FA,SAFINISH第六例干系结构的静力学分析——悬臂梁一工字悬臂梁,分析其在集中力p作用下自由端的变形,一直梁的材料为10号热轧工字钢,其横截面面积为A=14.345cm2,截面高度H=100mm,惯性矩I=245cm4,梁的长度L=1m,集中力p=10000N,钢的弹性模量E=2e11N/m2,泊松比为0.3/clear/prep7Et,1,beam3R,1,14.345e-4,245e-8,0.1Mp,ex,1,2e11Mp,prxy,1,0.3K,1,0,0,0Lstr,1,2Lesize,1,,,50Lmesh,1Finish/soluDk,1,uxDk,1,uyDk,1,rotzFk,2,fy,-10000SolveFinish/post1PldispFinish第七例平面问题的求解实例——厚壁圆筒问题一厚壁圆筒,其内半径r1=50mm,外半径r2=100mm,作用在内孔上的压力p=10MPa,无轴向压力,轴向长度很大可视为无穷,要求计算厚壁圆筒的径向应力和切向应力沿半径方向的分布。
/clear/filname,example7/prep7Et,1,plane183,,,2Mp,ex,1,2e11Mp,prxy,1,0.3Pcirc,0.1,0.05,0,90Lesize,4,,,6Lesize,3,,,8Mshape,0Mshkey,1Amesh,1Finish/soluDl,4,,uyDl,2,,uxSfl,3,pres,10e6Solve/post1Path,p1,2Ppath,1,30Ppath,2,1pdef,sr,s,xPdef,st,s,yPlpath,sr,stFinish第八例扳手的受力分析/clear/filname,example8/prep7Et,1,plane42Et,2,solid45Mp,ex,1,2e11Mp,prxy,1,0.3Rpr4,6,0,0,0.01K,7,0,0,0K,8,0,0,0.05K,9,0,0.1,0.05Lstr,7,8Lstr,8,9Lfillt,7,8,0.015Lstr,1,4AsbL,1,10Lesize,2,,,3Lesize,3,,,3Lesize,4,,,3Lesize,7,0.01Lesize,8,0.01Lesize,9,0.01 Mshape,0Mshkey,1Amesh,allVdrag,all,,,,,,7,9,8 Aclear,2,3,1Finish/soluDa,2,allDa,3,allKsel,s,,,24,29,1Fk,all,fx,100Ksel,allSolveFinish/post1/view,1,1,1,1Pldisp,2Plnsol,s,eqv,0,1Nwpave,159/type,1,sect/cplane,1/replotfinish第九章圆轴扭转分析设等直圆轴的圆截面直径D=50mm,长度L=120mm,作用在圆轴两端上的转矩Mn=1.5xe3N.m/clear/filname,example9/prep7Et,1,plane183Et,2,solid186Mp,ex,1,2.08e11Mp,prxy,1,0.3Rectng,0,0.025,0,0.12Lesize,1,,,5Lesize,2,,,8Mshape,0Mshkey,1Amesh,1Extopt,esize,5Extopt,aclear,1Vrotat,1,,,,,,1,4,360/view,1,1,1,1Wprot,0,-90Cswpla,11,1,1,1Nsel,s,loc,x,0.025Rpotat,allFinish/soluD,all,uxNsel,r,loc,z,0.12F,all,fy,1500Allsel,allDa,6,allDa,10,allDa,14,allSolveFinish/post1Pldisp,1Rsys,11Nsel,s,loc,z,0.045Nsel,r,loc,y,0Prnsol,u,ynsel,s,loc,z,0,0.045esln,r,1plesol,s,yzFinish第十例模态分析实例——均匀直杆的固有频率分析一根长度为L的等截面直杆,一端固定,一段自由。
ANSYS岩土计算例子
ANSYS土工结构计算案例ANSYS-CHINA广州办事处2020年4月24日目录计算题目及计算要求说明 (1)题目一 (4)一、计算说明 (4)二、计算所用ANSYS邓肯-张的E-B模型说明 (5)三、计算有限元模型及计算结果 (6)题目二 (7)一、用三维有限元模型计算 (7)二、用三维有限元模型计算 (8)题目三 (10)一、计算说明 (10)二、计算有限元模型及边界条件 (10)三、强夯地基固结计算 (10)题目四 (17)一、计算说明 (17)二、计算几何模型和有限元模型 (17)三、计算结果 (18)1、计算边界条件 (18)2、计算结果 (19)3、结论 (20)计算题目及计算要求说明题目一:高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析题目二:大圆筒结构、波浪与地基的相互作用分析(大圆筒作为重力式码头结构,波浪为水平动荷载,门吊为竖向动荷载,地基为三层以上地基包括自抛碎石垫层、粘土层、粉细砂层和岩层,粉细砂层可能在波浪动荷载作用下液化造成圆筒倾覆)题目三:(冲击荷载下)强夯地基固结有限元分析(提供固结方程或固结方程处理方案,孔隙水压力消散计算方案、沉降计算方案及其他一些处理技巧)题目四:在降雨情况下土工格栅加筋土挡墙边坡上公路稳定分析(由上至下为公路面层,垫层,挡墙,挡墙面板采用预制混凝土块,混凝土后方为钩挂式土工格栅,边坡比较陡,边坡有一定排水特性)。
具体处理方案包括:1、提供计算输入界面2、计算模型或采用本构情况3、前处理方案及网格划分技巧4、特殊材料或模型嵌入技术5、计算技巧及解决方案6、后处理提供内容具体报价方案包括:进行简单报价,涉及以上题目的各模块的综合报价(包括前后处理及解决以上问题的结构与岩土问题的模块报价. 能用通用模块计算尽量用通用模块,必须用CivilFEM模块计算的请注明.提供解决方案时间:2005年5月18日之前。
题目一高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析一、计算说明高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析, 可以用空间有限元模型或平面有限元模型数值求解。
ANSYS实例分析75道(含结果)
【ANSYS 算例】3.4.2(1) 基于图形界面的桁架桥梁结构分析(step by step)下面以一个简单桁架桥梁为例,以展示有限元分析的全过程。
背景素材选自位于密执安的"Old North Park Bridge" (1904 - 1988),见图3-22。
该桁架桥由型钢组成,顶梁及侧梁,桥身弦杆,底梁分别采用3种不同型号的型钢,结构参数见表3-6。
桥长L=32m,桥高H=5.5m 。
桥身由8段桁架组成,每段长4m 。
该桥梁可以通行卡车,若这里仅考虑卡车位于桥梁中间位置,假设卡车的质量为4000kg ,若取一半的模型,可以将卡车对桥梁的作用力简化为P 1 ,P 2和P 3 ,其中P 1= P 3=5000 N, P 2=10000N ,见图3-23。
图3-22位于密执安的"Old North Park Bridge" (1904 - 1988)图3-23 桥梁的简化平面模型(取桥梁的一半)表3-6 桥梁结构中各种构件的几何性能参数构件 惯性矩m 4 横截面积m 2顶梁及侧梁(Beam1) 643.8310m -⨯322.1910m -⨯ 桥身弦梁(Beam2) 61.8710-⨯31.18510-⨯ 底梁(Beam3)68.4710-⨯ 33.03110-⨯解答 以下为基于ANSYS 图形界面(Graphic User Interface , GUI)的菜单操作流程。
(1) 进入ANSYS (设定工作目录和工作文件)程序 → ANSYS → ANSYS Interactive → Working directory (设置工作目录)→ Initial jobname (设置工作文件名):TrussBridge → Run → OK(2) 设置计算类型ANSYS Main Menu :Preferences … → Structural → OK(3) 定义单元类型ANSYS Main Menu:Preprocessor → Element Type → Add/Edit/Delete... → Add…→ Beam: 2d elastic 3 → OK(返回到Element Types窗口)→ Close(4) 定义实常数以确定梁单元的截面参数ANSYS Main Menu: Preprocessor → Real Constants…→ Add/Edit/Delete → Add…→ select Type 1 Beam 3 → OK → input Real Constants Set No. : 1 , AREA: 2.19E-3,Izz: 3.83e-6(1号实常数用于顶梁和侧梁) → Apply → input Real Constants Set No. : 2 , AREA: 1.185E-3,Izz: 1.87E-6 (2号实常数用于弦杆) → Apply → input Real Constants Set No. : 3, AREA: 3.031E-3,Izz: 8.47E-6 (3号实常数用于底梁) → OK (back to Real Constants window) →Close(the Real Constants window) (5) 定义材料参数ANSYS Main Menu: Preprocessor → Material Props → Material Models → Structural → Linear → Elastic → Isotropic → input EX: 2.1e11, PRXY: 0.3(定义泊松比及弹性模量) → OK → Density (定义材料密度) →input DENS: 7800, → OK → Close(关闭材料定义窗口)(6) 构造桁架桥模型生成桥体几何模型ANSYS Main Menu:Preprocessor → Modeling → Create → Keypoints → In Active CS → NPT Keypoint number:1,X,Y,Z Location in active CS:0,0 → Apply →同样输入其余15个特征点坐标(最左端为起始点,坐标分别为 (4,0), (8,0), (12,0), (16,0), (20,0), (24,0), (28,0), (32,0), (4,5.5), (8,5.5), (12,5.5), (16.5.5), (20,5.5), (24,5.5), (28,5.5))→Lines →Lines →Straight Line →依次分别连接特征点→ OK网格划分ANSYS Main Menu: Preprocessor → Meshing → Mesh Attributes → Picked Lines →选择桥顶梁及侧梁→ OK → select REAL: 1, TYPE: 1 → Apply →选择桥体弦杆→ OK → select REAL: 2, TYPE: 1 → Apply →选择桥底梁→ OK → select REAL: 3, TYPE:1 → OK →ANSYS Main Menu:Preprocessor → Meshing → MeshTool →位于Size Controls下的Lines:Set → Element Size on Picked → Pick all → Apply → NDIV:1 → OK → Mesh → Lines → Pick all → OK (划分网格) (7) 模型加约束ANSYS Main Menu: Solution → Define Loads → Apply → Structural→ Displacement → On Nodes →选取桥身左端节点→ OK → select Lab2: All DOF(施加全部约束)→ Apply →选取桥身右端节点→ OK → select Lab2: UY(施加Y方向约束)→ OK(8) 施加载荷ANSYS Main Menu: Solution → Define Loads → Apply → Structural → Force/Moment → On Keypoints →选取底梁上卡车两侧关键点(X坐标为12及20)→ OK → select Lab: FY,Value: -5000 → Apply →选取底梁上卡车中部关键点(X坐标为16)→ OK → select Lab: FY,Value: -10000 → OK →ANSYS Utility Menu:→ Select → Everything(9) 计算分析ANSYS Main Menu:Solution → Solve → Current LS → OK(10) 结果显示ANSYS Main Menu:General Postproc → Plot Results → Deformed shape → Def shape only →OK(返回到Plot Results)→ Contour Plot → Nodal Solu → DOF Solution, Y-Component of Displacement → OK(显示Y方向位移UY)(见图3-24(a))定义线性单元I节点的轴力ANSYS Main Menu →General Postproc →Element Table → Define Table →Add →Lab: [bar_I], By sequence num: [SMISC,1] → OK → Close定义线性单元J节点的轴力ANSYS Main Menu →General Postproc →Element Table → Define Table →Add →Lab: [bar_J], By sequence num: [SMISC,1] → OK → Close画出线性单元的受力图(见图3-24(b))ANSYS Main Menu → General Postproc → Plot Results → Contour Plot → Line Elem Res →LabI: [ bar_I], LabJ: [ bar_J], Fact: [1] → OK(11) 退出系统ANSYS Utility Menu:File → Exit → Save Everything → OK(a)桥梁中部最大挠度值为0.003 374m (b)桥梁中部轴力最大值为25 380N图3.24 桁架桥挠度UY以及单元轴力计算结果【ANSYS算例】3.4.2(2) 基于命令流方式的桁架桥梁结构分析!%%%%% [ANSYS算例]3.4.2(2) %%%%% begin %%%%%%!------注:命令流中的符号$,可将多行命令流写成一行------/prep7 !进入前处理/PLOPTS,DATE,0 !设置不显示日期和时间!=====设置单元和材料ET,1,BEAM3 !定义单元类型R,1,2.19E-3,3.83e-6, , , , , !定义1号实常数用于顶梁侧梁R,2,1.185E-3,1.87e-6,0,0,0,0, !定义2号实常数用于弦杆R,3,3.031E-3,8.47E-6,0,0,0,0, !定义3号实常数用于底梁MP,EX,1,2.1E11 !定义材料弹性模量MP,PRXY,1,0.30 !定义材料泊松比MP,DENS,1,,7800 !定义材料密度!-----定义几何关键点K,1,0,0,, $ K,2,4,0,, $ K,3,8,0,, $K,4,12,0,, $K,5,16,0,, $K,6,20,0,, $K,7,24,0,, $K,8,28,0,, $K,9,32,0,, $K,10,4,5.5,, $K,11,8,5.5,, $K,12,12,5.5,, $K,13,16,5.5,, $K,14,20,5.5,, $K,15,24,5.5,, $K,16,28,5.5,,!-----通过几何点生成桥底梁的线L,1,2 $L,2,3 $L,3,4 $L,4,5 $L,5,6 $L,6,7 $L,7,8 $L,8,9!------生成桥顶梁和侧梁的线L,9,16 $L,15,16 $L,14,15 $L,13,14 $L,12,13 $L,11,12 $L,10,11 $L,1,10!------生成桥身弦杆的线L,2,10 $L,3,10 $L,3,11 $L,4,11 $L,4,12 $L,4,13 $L,5,13 $L,6,13 $L,6,14 $L,6,15 $L,7,15 $L,7,16 $L,8,16!------选择桥顶梁和侧梁指定单元属性LSEL,S,,,9,16,1,LATT,1,1,1,,,,!-----选择桥身弦杆指定单元属性LSEL,S,,,17,29,1,LATT,1,2,1,,,,!-----选择桥底梁指定单元属性LSEL,S,,,1,8,1,LATT,1,3,1,,,,!------划分网格AllSEL,all !再恢复选择所有对象LESIZE,all,,,1,,,,,1 !对所有对象进行单元划分前的分段设置LMESH,all !对所有几何线进行单元划分!=====在求解模块中,施加位移约束、外力,进行求解/soluNSEL,S,LOC,X,0 !根据几何位置选择节点D,all,,,,,,ALL,,,,, !对所选择的节点施加位移约束AllSEL,all !再恢复选择所有对象NSEL,S,LOC,X,32 !根据几何位置选择节点D,all,,,,,,,UY,,,, !对所选择的节点施加位移约束ALLSEL,all !再恢复选择所有对象!------基于几何关键点施加载荷FK,4,FY,-5000 $FK,6,FY,-5000 $FK,5,FY,-10000/replot !重画图形Allsel,all !选择所有信息(包括所有节点、单元和载荷等)solve !求解!=====进入一般的后处理模块/post1 !后处理PLNSOL, U,Y, 0,1.0 !显示Y方向位移PLNSOL, U,X, 0,1.0 !显示X方向位移!------显示线单元轴力------ETABLE,bar_I,SMISC, 1ETABLE,bar_J,SMISC, 1PLLS,BAR_I,BAR_J,0.5,1 !画出轴力图finish !结束!%%%%% [ANSYS算例]3.4.2(2) %%%%% end %%%%%%【ANSYS算例】3.2.5(3) 四杆桁架结构的有限元分析下面针对【典型例题】3.2.5(1)的问题,在ANSYS平台上,完成相应的力学分析。
ANSYS操作题练习
一、ANSYS操作题练习
建立如图所示的二维框架结构的有限元模型,其高度为5m,梯形下边长为6m,上边长为4m,中间层间隔2.5m,整体坐标系以图为准,支撑腿沿Z轴方向,2个支撑腿在Z=0处约束全部自由度,钢架杆件采用pipe16单元,四根桩腿的截面尺寸为外径0.2m x 厚度0.02m,其余杆件的尺寸为0.05m x 0.005m,材料为低碳钢,其弹性模量为2.06x1011Pa,泊松比为0.3,密度为7800kg/m3。
问题1:在如图2所示在静载荷作用下(Fx=100N)绘出结构的变形图以及等效应力图、求出支撑腿上端2个节点的位移以及支撑腿下端2个节点的支反力。
在E盘建立一个以自己姓名命名的文件夹并保存模型db文件和输出结果文件。
问题2:计算这个模型在没有载荷情况下的前6阶模态频率,模态分析方法选择subspace,把计算结果保存成一个文本文件存入上述文件夹。
ansys例题1
第二题
建模与仿真步骤
1.由于SolidWorks提供了十分强大的三维建模功能,并且输出的ParaSolid与
Ansys可以无缝对接,故此次水桶建模我利用SolidWorks2013来绘制,在绘制过程中要注意的一个细节是SW环境下的尺寸导入到ANSYS之后会被自动转换成米单位制,如原图中的100mm在ANSYS中会被转换成0.1(m),故建模与仿真均使用国际标准SI单位制即可。
首先选中正视面平面,添加草图,绘制剖面图。
2.绘制完毕之后旋转实体,得到水桶模型,建模效果图如下图所示。
之后输出
ParaSolid格式模型。
3.打开ANSYS软件,File-Inport导入刚才的ParaSolid模型,导入完毕后要
调整显示模式,Plot-Style-SolidModelFacets中,显示模式改为
NormalFaceting,之后可观察导入的模型。
4.下一步开始为仿真做准备,通过查表可得灰口铸铁的弹性模量及泊松比,注意本例中使用SI单位制,在单位换算后输入即可。
5.设置好element type为10node187,之后直接利用自适应网格尺寸划分得到划分结果。
6.固定桶的下表面。
加载,载荷大小为10个大气压,即1010Kpa。
7.运行后即可得到最终效果。
桶在10个大气压下的受力情况,图示为范米塞斯应力值。
下图为桶的应变值,由于桶边有卷起来的加强结构,故顶端应变值有减小,如果细心可以发现,一般生活中看到的水桶都会在桶沿与桶的中上部加固,我想看到这个应变结果,会发现在桶沿与桶的中上部加固是十分合理的。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
巷道开挖过程的有限元模型与力学分析学院:班级:学号:姓名:指导老师:XXX老师XX大学2009-2010年度上学期目录巷道开挖过程的有限元模型与力学分析 (3)有限元模型概述 (3)1 建立有限元模型 (4)1.1 定义工作文件和工作标题 (4)1.2 定义单元类型、实常数和材料 (4)1.3 建立几何模型 (5)2 网格划分 (7)3 加载与初始地应力模拟 (8)3.1 设置分析类型 (8)3.2施加边界条件 (9)3.3 施加上部面压力: (10)3.4 施加重力加速度: (10)3.4 设置加载步骤: (10)3.5 求解初始地应力: (11)3.6 保存分析结果: (11)4 浏览初始地应力的计算结果(后处理) (12)4.1显示变形形状: (12)4.2查看节点结果等值线图: (12)5 开挖巷道求解 (16)5.1杀死巷道对应的单元: (16)5.2输入加载步文件: (16)5.3查看杀死巷道单元后的受力情况: (16)5.4求解开挖后的有限元模型: (16)6 浏览查看开挖后的计算结果(后处理) (17)6.1显示变形形状: (17)6.2查看节点结果等值线图: (17)7 应力集中分析 (22)8 总结体会 (22)巷道开挖过程的有限元模型与力学分析某半圆形拱巷道断面,参数见图,其所处地质条件为IV级围岩,上覆盖层厚度为100米后,各材料的力学参数见表。
有限元模型概述本题采用ANSYS有限元分析软件模拟巷道开挖过程。
由于地下巷道属于细长结构物,即巷道的横断面相对于纵向的长度来说很小,可且假定在围着荷载作用下,在其纵向没有位移,只有横向发生位移. 所以,巷道的力学分析可以采用弹性力学理论中的平面应变模型进行,这是一个较复杂的非线性力学问题。
采用ANSYS有限元分析软件对巷道开挖进行模拟时,应首先根据地质条件建立合适的地下有限元分析模型,由于巷道对整个地下空间来说是属于“小孔口问题”,巷道周围出现孔口应力集中,并且应力集中区域影响范围约大于1.5倍的孔口尺寸,因此建几何模型时要选择合适的尺寸,其次采用PLANE42单元类型来分析平面应变问题。
接下来设定单元尺寸大小划分网格。
然后分步求解载荷,第一步求解初始地应力,第二步杀死巷道对应的单元后求解开挖后的载荷,最后分析巷道周围岩石的、位移、应力、应变的变化。
1建立有限元模型1.1 定义工作文件和工作标题(1)启动ANSYS 程序.(2)定义工作文件名:执行Utility Menu >File >Change Jobname命令,在弹出的对话框中输入“xiangdaokaiwa”并选择YES复选框,单击OK按钮。
(3)定义工作标题:执行Utility Menu>File>Change Title命令,在弹出的对话框中输入“xiangdaokaiwamodel”单击OK。
1.2 定义单元类型、实常数和材料(1)选择分析模块:执行Main Menu>Preferences,选择<Stuctural>复选框,单击OK。
巷道开挖属于结构模块。
(2)选择单元类型:执行Main Menu>Preprocessor>Element Type>Add/Edit/Delete命令,弹出<Element Type>对话框,单击Add 按钮,弹出<Library of Element>对话框。
选择”Structural solid”和“Quad 4node 42”选项,单击OK,在返回到<Element Type>对话框中单击<option>按钮,在弹出对话框中将单元力学模型的区K3设置成“Plane Strain”选项(即平面应变模式),再单击OK,最后单击单元类型中的<Close>按钮。
如图1所示。
这里的单元模型没有几何参数,不需要设置实常数。
图1 设置平面应变模式(3)设置材料属性:执行Main Menu> Preprocessor> Material Props> Material Models命令,在弹出的窗口中双击<Material Model Available> 列表框中的“Structural\Linear\Elastic\Isotropic”选项,弹出<Linear Isotropic Material Prorerties for Material Number 1>对话框。
在<EX>和<PRXY>文本框中分别输入“1.6E9”及“0.32”,单击OK,。
然后双击<Material Model Available>列表中的<Density>,在弹出的对话框中的<DENSY> 文本框中输入“2800”,单击OK,最后退出。
如图2所示图2 设置弹性模量、泊松比及密度1.3 建立几何模型模型说明:根据巷道开挖受影响的范围,模拟地层的横向边界距离巷道墙壁可定为两倍巷道宽度(6.8m)即14m,则地层宽度为34.8m,约为35m,巷道的顶拱为三心拱,拱高为1/4墙高即1.7m,整个巷道高度8.5m,巷道埋深为100m,巷道底下设两倍巷道高度即20m。
考虑建模时简便,有限元模型的巷道顶部到模型边界为50m,模型上部再加50m的均布载荷,即p=ρgh=2800*10*50=1.400000Pa。
如图3所示。
图3 模型简图(1)创建三心拱巷道断面:由于用ANSYS创建三心拱难度较大,为了方便高效,这里采用AUTOCAD来创建三心拱巷道断面,并且用DXFTOANSYS软件将AUTOCAD输出的DXF文件转化为ANSYS能够读取的LGW文件。
将xiangdao.lgw文件导入ANSYS中得到图4。
图4 三心拱巷道断面图(2)创建边界:执行Main Menu>Preprocessor>Modeling>Create>Keypoints>In Active CS命令,弹出当前坐标创建关键点的对话框,名关键点分别是KP15(21,-20,0),KP16(-14,-20,0),KP17(21,58,0),KP18(-14,58,0)。
(3)创建直线:执行Main Menu >Preprocessor >Modeling>Create>Lines >Lines >Straight Line命令,将地层边界连接起来。
如图5所示。
图5 几何模型(4)创建地层面积:执行Main Menu >Preprocessor >Modeling >Create>Areas >Arbitrary>By Lines,创建三个叠加起来的面积。
(5)面积叠加的布尔运算:执行Main Menu>Preprocessor >Modeling >Operate >Booleans>Overlap>Areas 命令,将三个面积进行布尔叠加计算,得到三个独立的面积。
(6)压缩全部编号:执行Main Menu>Preprocessor>Numbering Ctrls>Compress Numbers命令,弹出对话框中选择<All>选项,得到编号压缩的模型。
如图6所示。
图6完整的几何模型2网格划分(1)设置单元尺寸大小:执行Main Menu>Preprocessor>Meshing>Size Cntrls >ManualSize>Layers>Picked Lines选项,弹出以线来控制单元尺寸的对话框,选择要划分尺寸的线条,弹出<Element Sizes by Picked Lines>对话框,在对应线条划分数量的文本框中输入一个合适的数字,L1线划分为8段,L9线划分为15段,L10线划分为20段。
划分尺寸时应在巷道周围划分相对较小的单元,使巷道内外的网格更加密集,计算更加精确。
(2)自由划分网格:执行Main Menu >Preprocessor >Meshing >Mesh >Areas >Free命令,弹出一个拾取框。
拾取面积A1、A2、A3,单击OK,完成自由网格单元划分。
如图7所示图7 自由网格划分(3)保存网格数据:单击ANSYS TOOLBAL中的<Save_db>按钮。
3加载与初始地应力模拟3.1 设置分析类型(1)设置求解控制:执行Main Menu>Solution>.Analysis type>Solution controls 命令,弹出如图所示对话框,如图8设置求解控制。
图8 设置求解控制(2)设定求解方法:执行Main Menu>Solution>.Analysis type>Analysis option 命令,弹出<Statis or Steady analysis>对话框,在“NROPT(NEWTON---RAPHSON OPTION)”下拉列表中选择“FULL N—P”选项,单击OK,如图9所示。
Newton-Raphson是一种求解非线性方程组的数值方法,简称N-R法。
要运用单元生死,就要设定全牛顿-拉普森选项,才能有较好的效果。
图9 在静态与稳态分析对话框中设定全牛顿-拉普森选项3.2施加边界条件模型水平左右方向位移为零,底部为固定约束,垂直方向为自由滑动,上部边界受垂直地压力作用。
(1)横向X方向约束和底部全约束:执行Main Menu>Solution>Define Loads>Apply>Structural>Displacement>On Lines 命令,弹出一个拾取框。
拾取L10和L12,单击OK,弹出如图所示的<Apply U,Rot on Lines>对话框。
选择<UX>,并在“Value”中赋值0。
单击〈Apply〉再次拾取L9,单击OK,弹出如图10所示的<Apply U,Rot on Lines>对话框。
选择<All Dof>,并在“Value”中赋值0。
最后单击Ok。
图10 横向X方向约束(2)关键点约束:执行Main Menu>Preprocessor>Loads>DefineLoads>Apply >Structural>Displacement>On Keypoints命令,弹出一个拾取框。