ansys算例作业

ANSYS实例分析75道(含结果)【【ANSYS算例算例】】3.4.2(1)基于图形界面的桁架桥梁结构分析基于图形界面的桁架桥梁结构分析(stepbystep)下面以一个简单桁架桥梁为例，以展示有限元分析的全过程。

背景素材选自位于密执安的“OldNorthParkBridge“(1904-1988)，见图3-22。

该桁架桥由型钢组成，顶梁及侧梁，桥身弦杆，底梁分别采用3种不同型号的型钢，结构参数见表3-6。

桥长L=32m,桥高H=5.5m。

桥身由8段桁架组成，每段长4m。

该桥梁可以通行卡车，若这里仅考虑卡车位于桥梁中间位置，假设卡车的质量为4000kg，若取一半的模型，可以将卡车对桥梁的作用力简化为P1，P2和P3，其中P1=P3=5000N,P2=10000N，见图3-23。

图3-22位于密执安的“OldNorthParkBridge“(1904-1988)图3-23桥梁的简化平面模型(取桥梁的一半)表3-6桥梁结构中各种构件的几何性能参数构件惯性矩m4横截面积m2顶梁及侧梁桥身弦梁底梁解答解答以下为基于ANSYS 图形界面(GraphicUserInterface,GUI)的菜单操作流程。

(1)进入进入ANSYS（设定工作目录和工作文件）（设定工作目录和工作文件）程序程序→→ANSYS→→ANSYSInteractive→→Workingdirectory（设置工作目录）→Initialjobname（设置工作文件名）:TrussBridge→→Run→→OK(2)设置计算类型设置计算类型：Preferences…→→Structural→→OK(3)定义单元类型定义单元类型ANSYSMainMenu：Preprocessor→→ElementType→→Add/Edit/Delete.→→Add…→→Beam:2delastic3→→OK（返回到ElementTypes窗口）→→Close(4)定义实常数以确定梁单元的截面参数定义实常数以确定梁单元的截面参数ANSYSMainMenu:Preprocessor→→RealConstants…→→Add/Edit /Delete→→Add…→→selectType1Beam3→→OK→→RealConsta ntsSetNo.:1,AREA:2.19E-3，，Izz:3.83e-6(1号实常数用于顶梁和侧梁)→→Apply→→RealConstantsSetNo.:2,AREA:1.185E-3，，Izz:1.87E-6(2号实常数用于弦杆)→→Apply→→RealConstantsSetNo.:3,AREA:3.031E-3，，Izz:8.47E-6(3号实常数用于底梁)→→OK(backtoRealConstantswindow)→Close(theRealConstant swindow)(5)定义材料参数定义材料参数ANSYSMainMenu:Preprocessor→→MaterialProps→→MaterialMo dels→→Structural→→Linear→→Elastic→→Isotropic→→EX:2.1e11,PRXY:0.3(定义泊松比及弹性模量)→→OK→→Density(定义材料密度)→DENS:7800,→→OK→→Close（关闭材料定义窗口）(6)构造桁架桥模型构造桁架桥模型生成桥体几何模型ANSYSMainMenu：Preprocessor→→Modeling→→Create→→Keypoints→→InActive CS→→NPTKeypointnumber：：1，，X，，Y，，ZLocationinactiveCS：：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→→StraightLine→→依次分别连接特征点→→OK网格划分ANSYSMainMenu:Preprocessor→→Meshing→→MeshAttributes→→PickedLines→→选择桥顶梁及侧梁→→OK→→selectREAL:1,TYPE:1→→Apply→→选择桥体弦杆→→OK→→selectREAL:2,TYPE:1→→Apply→→选择桥底梁→→OK→→selectREAL:3,TYPE:1→→OK→→ANSYSMainMen u：Preprocessor→→Meshing→→MeshTool→→位于SizeControls 下的Lines：：Set→→ElementSizeonPicked→→Pickall→→Apply→→NDIV：：1→→OK→→Mesh→→Lines→→Pickall→→OK(划分网格)(7)模型加约束模型加约束ANSYSMainMenu:Solution→→DefineLoads→→Apply→→Struct ural→→Displacement→→OnNodes→→选取桥身左端节点→→OK→→selectLab2:AllDOF(施加全部约束)→→Apply→→选取桥身右端节点→→OK→→selectLab2:UY(施加Y方向约束)→→OK(8)施加载荷施加载荷ANSYSMainMenu:Solution→→DefineLoads→→Apply→→Struct ural→→Force/Moment→→OnKeypoints→→选取底梁上卡车两侧关键点（X坐标为12及20）→→OK→→selectLab:FY，，Value:-5000→→Apply→→选取底梁上卡车中部关键点（X坐标为16）→→OK→→selectLab:FY，，Value:-10000→→OK→→ANSYSUtilityMenu：→→Select→→Everything(9)计算分析计算分析ANSYSMainMenu：Solution→→Solve→→CurrentLS→→OK(10)结果显示结果显示ANSYSMainMenu：GeneralPostproc→→PlotResults→→Deedshape→→Defshapeonly →→OK（返回到PlotResults）→→ContourPlot→→NodalSolu→→DOFSolution,Y-Componentof Displacement→→OK（显示Y方向位移UY）(见图3-24(a))定义线性单元I节点的轴力ANSYSMainMenu→GeneralPostproc→→ElementTable→→Define Table→→Add→→Lab:[bar_I],Bysequencenum:[SMISC,1]→→OK →→Close定义线性单元J节点的轴力ANSYSMainMenu→→GeneralPostproc→→ElementTable→→Def ineTable→→Add→→Lab:[bar_J],Bysequencenum:[SMISC,1]→→OK→→Close画出线性单元的受力图(见图3-24(b))ANSYSMainMenu→→GeneralPostproc→→PlotResults→→ContourPlot→→LineElemRes→→LabI:[bar_I],LabJ:[bar_J],Fact :[1]→→OK(11)退出系统退出系统ANSYSUtilityMenu：File→→Exit→→SaveEverything→→OK(a)桥梁中部最大挠度值为0.003374m(b)桥梁中部轴力最大值为25380N图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,2 4,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 算例】3.3.7(3) 三梁平面框架结构的有限元分析针对【典型例题】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支反力计算结果一致。

一、问题描述某高压容器设计压力p=16Mpa，设计温度T=200℃，材料为16MnR。

筒体内径R1=775mm，壁厚t1=100mm；封头内径R2=800mm，厚度t2=48mm，筒体削边长度L=95mm。

试对该压力容器筒体与封头的连接区域进行应力分析。

二、问题分析由于主要讨论封头与筒体过渡区域的应力状况，忽略封头上的其他结构如开孔接管等，建立如图1的有限元力学模型，其中筒体长度远远大于边缘应力的衰减长度，此处去筒体长度Lc=1200mm。

有限元计算采用Plane82单元，并设定轴对称选项。

筒体下端各节点约束轴向位移，球壳对称面上各节点约束水平方向位移，内壁施加均匀压力面载荷。

图 1 有限元分析模图2 有限元网格模型图3 均布载荷及边界设置图 三、ANSYS 分析计算 （1）建立有限元模型选取Plane82单元，建立几何模型并划分网格（四边形），结果如图2所示：（2） 施加载荷与约束安装前述分析，在网格模型上施加已知载荷与简化约束，如图3所示。

(3) 求解及结果后处理将载荷和约束正确施加后，检查无误，即可进入/SOL模块求解，求解完成后，进入/POST21后处理器进行必要的后处理操作。

1）变形（位移）图2）应力云图。

最大Mises应力为183.667Mpa。

四、分析过程命令流FINISH/CLEAR,start/Filn,E4-1/Title,VORTEX/Units,si!*******参数设定***** R1=775T1=100R2=800T2=48L=95LC=1200P=16E=2E5NU=0.3NT=5 !厚度方向剖分数NS=30 !球壳经向剖分数NC=30 !筒体轴向剖分数NL=5 !过渡段剖分数RA=0.6 !剖分比例!*******前处理******* /PREP7et,1,82keyopt,1,3,1mp,ex,1,Emp,nuxy,1,NUcyl4,,,R2,0,R2+T2,90blc4,R1,0,T1,-LC+L WPROT,,-90 WPOFF,,,LASEL,S,LOC,Y,0,R2+T2 ASBW,ALLASEL,S,LOC,Y,0,L ADELE,ALL,,,1K1=KP(R1,0,0)K2=KP(R1+T1,0,0) KSEL,S,LOC,Y,L*GET,K3,KP,,NUM,MIN *GET,K4,KP,,NUM,MAX ALLSA,K1,K2,K3,K4!网格划分LSEL,S,LOC,Y,-LC+L LSEL,A,LOC,Y,0 LSEL,A,LOC,Y,L LSEL,A,LOC,X,0 LESIZE,ALL,,,NT,,,,,1 LSEL,S,LENGTH,,LC-L LSEL,R,LOC,X,R1 LESIZE,ALL,,,NC,1/RA,,,,1 LSEL,S,LENGTH,,LC-L LSEL,R,LOC,X,R1+T1 LESIZE,ALL,,,NC,RA,,,,1 LSEL,S,RADIUS,,R2 LSEL,A,RADIUS,,R2+T2 LESIZE,ALL,,,NS,RA,,,,1LSEL,S,LOC,Y,0,LLSEL,U,LOC,Y,0LSEL,U,LOC,Y,LLESIZE,ALL,,,NL,,,,,1ALLSAMESH,ALLFINI!求解/SOLUNseL,S,LOC,Y,-LC+LD,ALL,UYNSEL,S,LOC,X,0D,ALL,UXLSEL,S,LINE,,3,8,5LSEL,A,LINE,,10NSLL,S,1SF,ALL,PRES,PALLSSOL VEFINI!后处理/POST1PLNSOL,S,INT,0,1 !显示应力云图FINI。

Ansys有限元课程设计问题一：飞机机翼振动模态分析机翼模型沿着长度方向具有不规则形状，而且其横截面是由直线和曲线构成（如图所示）。

机翼一端固定于机身上，另一端则自由悬挂。

机翼材料的常数为：弹性模量E=0.26GPa，泊松比m=0.3，密度r=886kg/m^3一、操作步骤：1.选取5个keypoint，A（0,0,0）为坐标原点，同时为翼型截面的尖点；2.B（2，0，0）为下表面轮廓截面直线上一点，同时是样条曲线BCDE的起点;3.D(1.9,0.45,0)为样曲线上一点；4.C(2.3，0.2,0)为样条曲线曲率最大点，样条曲线的顶点；5.E(1,0.25,0)与点A构成直线，斜率为0.25；6.通过点A、B做直线和点B、C、D、E作样条曲线就构成了截面的形状。

沿Z 方向拉伸，就得到机翼的实体模型；7.创建截面如图：机翼材料的常数为：弹性模量E=0.26GPa，泊松比m=0.3，密度r=886kg/m^3 8.定义网格密度并进行网格划分：选择面单元PLANE42和体单元SOLID45进行划分网格求解。

面网格选择单元尺寸为0.00625，体网格划分时按单元数目控制网格划分，选择单元数目为109.对模型施加约束，由于机翼一端固定在机身上所以在机翼截面的一端所有节点施加位移和旋转约束二、有限元处理结果及分析：机翼的各阶模态及相应的变形：一阶振动模态图：二阶振动模态图：三阶振动模态图：四阶振动模态图：五阶振动模态图：命令流：/FILNAM,MODAL/TITLE,Modal analysis of a modal airplane wing /PMETH,OFF,0KEYW,PR_STRUC,1/UIS,MSGPOP,3/PREP7ET,1,PLANE42ET,2,SOLID45MP,EX,1,380012MP,PRXY,1,0.3MP,DENS,1,1.033E-3K,1,K,2,2K,3,2.3,0.2K,4,1.9,0.45K,5,1,0.25/TRIAD,OFF/PNUM,KP,1LSTR,1,2LSTR,5,1BSPLIN,2,3,4,5,,,-1,0,,-1,-0.25,, AL,1,2,3ESIZE,0.25MSHKEY,0MSHAPE,0,2DAMESH,1SAVEESIZE,,10TYPE,2VEXT,1,,,0,0,10/SOLUANTYPE,MODAL MODOPT,SUBSP,5,,,,OFF EQSLV,SPARMXPAND,5,,,,0.001 LUMPM,0PSTRES,0ESEL,U,TYPE,,1NSEL,S,LOC,Z,0D,ALL,ALLALLSEL,ALLSOLVE/POST1SET,LISTSET,FIRSTPLDI,,ANMODE,10,0.5,,0FINISH13/EXIT,ALL问题二：内六角扳手静力分析内六角扳手在日常生产生活当中运用广泛，先受1000N的力产生的扭矩作用，然后在加上200N力的弯曲，分析算出在这两种外载作用下扳手的应力分布。

ANSYS土工结构计算案例ANSYS－CHINA广州办事处9:00 AM目录计算题目及计算要求说明 (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.6⨯0。

6⨯0.6m3,混凝土后方为钩挂式土工格栅，边坡比较陡，边坡有一定排水特性）。

具体处理方案包括：1、提供计算输入界面2、计算模型或采用本构情况3、前处理方案及网格划分技巧4、特殊材料或模型嵌入技术5、计算技巧及解决方案6、后处理提供内容具体报价方案包括:进行简单报价，涉及以上题目的各模块的综合报价（包括前后处理及解决以上问题的结构与岩土问题的模块报价。

能用通用模块计算尽量用通用模块,必须用CivilFEM模块计算的请注明.提供解决方案时间：2005年5月18日之前。

题目一高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析一、计算说明高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析，可以用空间有限元模型或平面有限元模型数值求解。

编号显示设置对话框
选择菜单Solution >Define Loads>A
pply>Structure>
Force/Moment O n Keypoints，载 荷分别如下：8 点承受Z方向集 中载荷Fz=5000 N和Y方向集中 载荷Fy=-5000N； 10点承受X方向 集中载荷Fx=50 00N；3点承受Z 方向集中载荷Fz =-5000N；6点 承受X方向集中 载荷Fx=-5000N。 施加载荷，如图 所示。
圆柱实体示意图
（1）选择单元类型 运行Preprocessor>Element Type>Add/Edit/Delete，弹出Element Types
对话框，如图所示。
由于12.0版本后对单元类型进行了 合并，之前的很多单元类型在12.0 以后在添加页面不见了，但是可以 用命令流的形式调用。格式如下：
圆柱实体示意图
（7）求解 选择Solution>Solve>Current LS， 开始计算，计算结束会弹出计算完毕对话框，单击Close。 关闭对话框计算完毕。
（8）后处理 运行 General Postproc>Plot Results>Contour Plot>Nodal Solu。 弹出如图所示对话框，运行DOF Solution>Displacement vector sum和 Stress>von Mises stress，分别显示圆柱体的位移和应力云图。
ET,1,45
单元类型对话框
然后单击Add，弹出Library of Element Types窗口，如图所 示，选择SOLID45单元，单击OK。
单元类型库对话框
（2）设置材料属性 运行Preprocessor>Material Props>Material Models，弹出如图所示对话框。

1.考虑热应力的模态分析2.ansys中抓图的命令流3.生死单元的例子4.正在建设中5.热应变与结构应变6.表面单元的使用7.修改ANSYS的默认工具条8.表格加载函数加载9.特征值屈曲分析10.包辛格效应11.塑性理论12.使用*vwrite和*vread13.ansys数值变量转字符变量14.ansys中保存视图设置15.ansys中最值问题16.apdl命令分类17.ANSYS宏加密18.用路径得到任一点的应力值19.ansys中的函数20.塑性应变实例21.残余应力22.beam4与solid45的连接23.塑性分析算例24.模态缩减法25.离心应力刚化26.ansys中泊松比效应27.比较弹性模量算例28.获得单元积分点的结果29.等向强化与随动强化小算例30.网格密度与计算结果31.ansys工况组合32.谐响应算例33.安全系数云图34.常应变单元35.ANSYS网格装配36.ANSYS子模型技术37.ansys中螺栓联接的模拟38.ANSYS接触过盈分析39.ANSYS板金冲压算例40.ANSYS节点解单元解单元表41.ANSYS单元生死42.过盈装配算例43.ANSYS超单元子结构44.ANSYS点到面热辐射算例45.ANSYS面到面热辐射46.ANSYS蠕变算例47.ANSYS蠕变算例248.ANSYS转子坎贝尔图分析49.ANSYS P单元算例50.ansys路径算例51.循环对称结构的模态分析52.ansys复合材料算例53.ansys复合材料算例254.扭矩和转动惯量对模态的影响55.ANSYS接触反力提取56.ANSYS热分析后处理57.ANSYS热结构耦合算例58.ANSYS瞬态动力学小算例59.ANSYS谐响应动力学小算例60.ANSYS模态动力学小算例61.ANSYS旋转软化应力刚化算例62.link10-gap单元63.循环对称结构的模态分析64.Ansys疲劳算例65.ANSYS材料实验小算例66.ANSYS各向异性材料算例67.ANSYS焊接模拟68.港口起重机有限元实例69.ansys中由网格生成几何模型70.ansys相变分析算例71.ANSYS输出单元刚度矩阵72.ANSYS模拟单摆运动73.ansys点点接触算例74.ANSYS弹性地基梁分析75.ANSYS平面应力plane4276.平面应变plane4277.轴对称plane4278.ANSYS螺旋的建模79.ANSYS弹簧建模80.ANSYS两个DB合并81.Hill_VS_vonmises82.ANSYS将文本数据给数组83.ansys中定义变截面梁84.ansys裂纹模拟85.ANSYS粘塑性86.shell63和shell187.ANSYS超弹性88.ANSYS剪切锁定算例89.ANSYS体积锁定算例90.ANSYS材料多线性与材料非线性91.工程应力与真实应力92.导热系数对热分布的影响93.ANSYS双线性随动/等向强化应94.ANSYS双线性随动/等向强化应95.ANSYS双线性随动/等向强化应96.ANSYS双线性随动/等向强化应97.Chaboche应力控制非循环对98.循环强化: Chaboche+非99.ANSYS数组插值。

a n s y s有限元分析作业经典案例本页仅作为文档封面，使用时可以删除This document is for reference only-rar21year.March有限元分析作业作业名称输气管道有限元建模分析姓名陈腾飞学号 3070611062班级 07机制（2）班宁波理工学院题目描述：输气管道的有限元建模与分析计算分析模型如图1所示承受内压：1.0e8 PaR1=0.3R2=0.5管道材料参数：弹性模量E=200Gpa；泊松比v=0.26。

图1受均匀内压的输气管道计算分析模型（截面图）题目分析：由于管道沿长度方向的尺寸远远大于管道的直径，在计算过程中忽略管道的断面效应，认为在其方向上无应变产生。

然后根据结构的对称性，只要分析其中1/4即可。

此外，需注意分析过程中的单位统一。

操作步骤1.定义工作文件名和工作标题1.定义工作文件名。

执行Utility Menu-File→Chang Jobname-3070611062，单击OK按钮。

2.定义工作标题。

执行Utility Menu-File→Change Tile-chentengfei3070611062，单击OK按钮。

3.更改目录。

执行Utility Menu-File→ change the working directory –D/chen 2.定义单元类型和材料属性1. 设置计算类型ANSYS Main Menu: Preferences →select Structural → OK2.选择单元类型。

执行ANSYS Main Menu→Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 8node 82 →applyAdd/Edit/Delete →Add →select Solid Brick 8node 185 →OKOptions…→select K3: Plane strain →OK→Close如图2所示，选择OK接受单元类型并关闭对话框。

输气管道受力分析（ANSYS建模）任务和要求：按照输气管道的尺寸及载荷情况,要求在ANSYS中建模，完成整个静力学分析过程。

求出管壁的静力场分布。

要求完成问题分析、求解步骤、程序代码、结果描述和总结五部分。

所给的参数如下：材料参数：弹性模量E=200Gpa; 泊松比0.26；外径R₁=0.6m；内径R₂=0.4m；壁厚t=0.2m。

输气管体内表面的最大冲击载荷P为1Mpa。

四.问题求解（一）．问题分析由于管道沿长度方向的尺寸远大于管道的直径，在计算过程中忽略管道的端面效应，认为在其长度方向无应变产生，即可将该问题简化为平面应变问题，选取管道横截面建立几何模型进行求解。

（二）．求解步骤定义工作文件名选择Utility Menu→File→Chang Jobname 出现Change Jobname对话框，在[/FILNAM] Enter new jobname 输入栏中输入工作名LEILIN10074723，并将New log and eror file 设置为YES，单击[OK]按钮关闭对话框定义单元类型1)选择Main Meun→Preprocessor→Element Type→Add/Edit/Delte命令，出现Element Type 对话框，单击[Add]按钮，出现Library of Element types对话框。

2）在Library of Element types复选框选择Strctural、Solid、Quad 8node 82,在Element type reference number输入栏中出入1，单击[OK]按钮关闭该对话框。

3. 定义材料性能参数1）单击Main Meun→Preprocessor→Material Props→Material models出现Define Material Behavion 对话框。

选择依次选择Structural、Linear、Elastic、Isotropic选项，出现Linear Isotropic Material Properties For Material Number 1对话框。

ANSYS计算算例
• 拱形结构算例 – 铁路隧道 – 公路隧道 – 地铁区间隧道(矿山法施工) 矩形框架结构 – 地下街 – 地下停车场 – 地铁车站 圆形结构 – 地铁区间隧道(盾构法施工)
•
•
拱形结构分析步骤
某单线电化铁路隧道衬砌如下图
a. 岩体特性 23 kN/m3 岩体为 III 级围岩，隧道埋深100m；岩体重度 围岩的弹性反力系数 K 1200 MPa/m，基底围岩弹性反力系数K a 1.25K 。 b. 衬砌材料 采用 C20 混凝土；重度 h 25 kN/m3，弹性模量E c 27 GPa 混凝土衬砌轴心抗压强度标准值 f ck 13.5 Mpa 1.7 MPa。 f 混凝土轴心抗拉强度标准值
计算中考虑结构自重对结构内力的影响，还应在此输入材料的密度。 Nonlinear-Density
4、模型建立 （1）第一段圆弧的输入（圆心：0，－3.41，起点： 0，0，圆心角：45º ）
Menu menu：Preprocessor>Material Props>Modeling Create>Lines Arcs>By Cent & Radius
三、加载及约束
1、荷载计算及转化 (1) 荷载计算 对深埋隧道用《隧规》推荐的统计法公式计算得： 竖向围岩压力：q=54084.25Pa； 水平围岩压力：e＝0.15 q=8112.64Pa (2) 匀布荷载转化为等效节点荷载 将匀布荷载按如下步骤转化成等效节点力。 A、列出结构上所有节点坐标 GUI：Utility Menu>List>Nodes
在命令输入框中，先输入圆心坐标（0，-3.41）并回车，再输入起点 坐标（0，0）并回车后，弹出下图交互式对话框，并输入第一段圆弧的圆 心角45º 。

材料属性 是与几何模型无关的本
构属性，例如杨氏模量、密度等. 虽然材料属性并不与单元类型联 系在一起，但由于计算单元矩阵 时需要材料属性，ANSYS为了用 户使用方便，还是对每种单元类 型列出了相应的材料类型。 根据 不同的应用，材料属性可以是线 性或非线性的. 与单元类型及实 常数类似，一个分析中可以定多 种材料. 每种材料设定一个材料 编号. 对于本问题，只须定义一 种材料，这种材料只须定义一个 材料属性—杨氏模量 2.06e11.
创建好几何模型以后，就要准备单元 类型、实常数、材料属性，然后划分 网格.
4. 设定分析模块.
a. Main Menu: Preferences b. 选择 Structural. c. 选择 OK.
Space Structure Research Center, HIT, CHINA
4
练习 - 悬壁梁(续)
可以列出所有的反作用力.
Space Structure Research Center, HIT, CHINA
16
练习 - 连续梁
交互操作
解释
17.绘制弯矩图.
a. Main Menu: General Postproc > Element Table > Define Tabel b. 点击Add... c. 在Lab. 中填入MI， Item选择By sequence num，选 择SMISC，下面填入 SMISC,3 d. 点击OK.
交互操作 5. 设定单元类型相应选项.
a. Main Menu: Preprocessor > Element Type > Add/Edit/Delete b. 选择 Add . . . c. 左边单元库列表中选择 Beam. d. 在右边单元列表中选择 2 node 188 e. 在opition中，将K3改成 Quadratic Form. 解释

ANSYS及其工程应用大作业1．如图所示三维实体支架，材料的弹性模量为200GPa，支架由两个圆孔的内表面固定，在支架表面承受1000N/cm2的均匀压力荷载，要求绘制变形后形状，找出模型的应力－应变分布规律，试分析最有可能发生屈服的位置。

（给出命令流清单及相关结果图）2．图示的屋顶桁架的横截面积为21.5in2,由绿枞木构成，弹性模量为1.9×106 lb/in2。

用ansys计算每个结合点的位移、每个杆的应力以及支座处的反作用力，并验证得出的结果。

（给出命令流清单及相关结果图）3．图示为带方孔（边长为120mm）的悬臂梁，其上受部分均布载荷（p=10Kn/m）作用，试采用一种平面单元，对图示两种结构进行有限元分析，并就方孔的布置进行分析比较，如将方孔设计为圆孔，结果有何变化？（板厚为1mm，材料为钢）4．如图（a）所示简支吊车梁，梁上有移动荷载以1.0s/m的速度从梁的一端移动到另一端，计算在此过程中吊车梁的位移和应力响应。

其中，梁材料为钢材，弹性模量为2.0×1011Pa，波松比为0.3，密度为7800kg/m3，采用焊接“工”字型组合截面，如图（b）所示，其中W1=150mm，W2=300mm，T1=20mm，T2=10mm。

(a) (b)5．如图所示，一根直的细长悬臂梁，一端固定一端自由。

在自由端施加P=1lb 的载荷。

弹性模量=1.0e4psi，泊松比=0.0；L=100in，H=5in，B=2in；对该悬臂梁做特征值屈曲分析，并进行非线性载荷和变形研究。

研究目标为确定梁发生分支点失稳（标志为侧向的大位移）的临界载荷。

6．请用Ansys中的三维梁单元为下图所示结构的各部分设计横断面尺寸。

要求用空心管。

该结构用于支撑红绿灯，它所承受的风力为80mile/hour,红绿灯灯箱重10kg。

请写一份简要报告，叙述自己的最终设计。

计算分析报告应包括以下部分：A、问题描述及数学建模；B、有限元建模（单元选择、结点布置及规模、网格划分方案、载荷及边界条件处理、求解控制）C、计算结果及结果分析（位移分析、应力分析、正确性分析评判）D、多方案计算比较（结点规模增减对精度的影响分析、单元改变对精度的影响分析、不同网格划分方案对结果的影响分析等）ANSYS及其工程应用大作业7．如图所示三维实体支架，材料的弹性模量为200GPa，支架由两个圆孔的内表面固定，在支架表面承受1000N/cm2的均匀压力荷载，要求绘制变形后形状，找出模型的应力－应变分布规律，试分析最有可能发生屈服的位置。

ansys例题一：四杆桁架结构源于《有限元分析基础教程(曾攀)》1 基于图形界面(GUI)的交互式操作(step by step)(1) 进入ANSYS(设定工作目录和工作文件)程序→ANSYS →ANSYS Interactive →Working directory (设置工作目录) →Initial jobname(设置工作文件名): planetruss→Run →OK(2) 设置计算类型ANSYS Main Menu: Preferences…→Structural →OK(3) 选择单元类型ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete…→Add…→Link：2D spar 1 →OK (返回到Element Types窗口) →Close(4) 定义材料参数ANSYS Main Menu: Preprocessor →Material Props →Material Models→Structural →Linear →Elastic→Isotropic: EX:2.95e11 (弹性模量)，PRXY: 0 (泊松比) →OK →鼠标点击该窗口右上角来关闭该窗口(5) 定义实常数以确定单元的截面积ANSYS Main Menu: Preprocessor →Real Constants…→Add/Edit/Delete →Add →Type 1→OK→Real Constant Set No: 1 (第1号实常数), AREA: 1e-4 (单元的截面积) →OK→Close(6) 生成单元ANSYS Main Menu: Preprocessor →Modeling →Creat→Nodes→In Active CS→Node number 1 →X:0,Y:0,Z:0 →Apply→Node number 2 →X:0.4,Y:0,Z:0 →Apply→Node number 3 →X:0.4,Y:0.3,Z:0→Apply→Node number 4 →X:0,Y:0.3,Z:0→OKANSYS Main Menu: Preprocessor →Modeling →Create →Elements→Elem Attributes （接受默认值）→User numbered→Thru nodes→OK→选择node 1和node2→Apply→选择node 2和node3→Apply→选择node 1和node3→Apply→选择node 4和node3→Apply→OK(7) 模型施加约束和外载添加位移的约束，分别将1节点X和Y方向、2节点Y方向、4节点的X和Y方向位移约束。

ANSYS软件进行有限元计算实例工字钢梁结构静力分析一工字钢梁两端均为固定端，其截面尺寸为：l=1.0m，a=0.16m，b=0.2m，c=0.02m，d=0.03m。

试建立该工字钢梁的三维实体模型，并在考虑重力的情况下对其进行结构静力分析。

其他已知参数如下：弹性模量E=206GPa；泊松比μ=0.3；材料密度ρ=7800kg/m3；重力加速度g=9.8m/s2；作用力作用于梁的上表面沿长度方向的中线处，其大小为F y=-5000N。

1)单元类型、几何特性、材料特性定义a）定义单元类型：Main Menu: Preprocessor→Element Type→Add/Edit/Delete弹出对话框，单击对话框中的“Add…”按钮，又弹出一对话框，选中其中的“Solid”和“Brick 8node 45”选项，单击“OK”按钮，关闭该对话框返回至上一级对话框。

单击“Close”按钮，关闭该级对话框。

b)定义材料特性：Main Menu: Preprocessor→Material Props→Material Models弹出对话框; 逐级双击右侧框中的Structural→Linear→Elastic →Isotropic，弹出下一级对话框。

在“弹性模量”（EX）文本框中输入“2.06e11”；在“泊松比”（PRXY）文本框中输入“0.3”；单击“OK”按钮，关闭该对话框返回至上一级对话框。

双击右侧框中的Density选项，在弹出的对话框中的“DENS”一栏中输入材料密度“7800”，单击“OK”按钮，关闭该对话框返回至上一级对话框。

关闭材料特性定义对话框。

2)三维实体模型的建立生成关键点●Main Menu: Preprocessor→Modeling →Create →Keypoints→In Active cs弹出对话框; 在Keypoint number 一栏中输入关键点编号“1”，在“X，Y，Z Location inactive cs”一栏中输入关键点1的坐标（-0.08，0，0），单击“Apply”按钮。

简支箱梁约束扭转算例一、工程背景已知某预应力混凝土简支箱梁，计算跨径为40m，沿梁长等截面。

截面尺寸如图1所示。

采用C40混凝土，剪切模量G=1.445×104MPa，弹性模量E=3.40×104MPa。

荷载为跨中作用一偏心荷载P=451.0kN，偏心距e=2.35m(计算约束扭转时，可以简化为集中力矩M k=1060.0kN)。

具体分两个工况进行：（1）跨中截面腹板位置作用一对对称集中竖向荷载，荷载大小为P/2=225.5kN；（2）跨中截面腹板位置作用一对反对称集中竖向荷载，荷载大小为P/2=225.5kN。

分别计算跨中截面、1/4跨位置截面上的正应力与剪应力分布，并绘制相应的正应力和剪应力分布曲线。

图1 箱梁截面尺寸(尺寸单位：cm)二、Ansys计算分析采用壳单元计算，荷载采用一对称荷载和反对称荷载加载。

1、Ansys命令流命令流见附录。

2、计算结果图a.工况一(a) 1/2跨正应力云图和应力曲线图(b) 1/2跨剪应力云图和应力曲线图b.工况二(a) 1/4跨正应力云图和应力曲线图(b) 1/4跨剪应力云图和应力曲线图注：由于路径选择时重复了一段，因此右翼缘由应力重叠现象附录命令流finish/clear/PREP7 ANTYPE,STATIC !定义箱梁厚度（单位：米）t1=0.22t2=0.30t3=0.34ET,1,SHELL63!定义单元R,1,t1R,2,t2R,3,t3!定义实常数跟所选单元有关本单元采用壳单元只需要厚度MP,EX,1,3.4E10!弹性模量1为材料参考号MP,PRXY,1,0.3!泊松比K,1,-4.75,0.955K,2,-2.35,0.955K,3,-2.35,-1.165K,4,2.35,-1.165K,5,2.35,0.955K,6,4.75,0.955K,7,2.35,0.955,40!创建关键点L,1,2$L,2,3$L,3,4L,4,5$L,5,6$L,5,2$L,5,7!连接各点ADRAG,1,6,5,,,,7ADRAG,2,3,4,,,,7!创建面根据直线编号2,3 路径线7ASEL,S,LOC,Y,0.955$AATT,1,1,1!选择面S选择新面作为子集loc按坐标值选择ASEL,S,LOC,X,-2.35$AA TT,1,2,1!Y中心坐标y值为0.955的面ASEL,S,LOC,X,2.35$AATT,1,2,1!设置面的单元属性相应定义点线体的命令不同ASEL,S,LOC,Y,-1.165$AATT,1,3,1!材料号，实常数，单元类型ALLSEL,ALL!选择所有图素NUMMRG,ALL !合并实体单元NUMCMP,ALL!压缩定义编号/PNUM,LINE,1/PNUM,AREA,1!显示面和线的编号LSEL,S,LENGTH,,40!选择长度为40的线段LESIZE,ALL,,,100!将所有的线每根分成100份即上述长度为40的线段LSEL,S,LENGTH,,2.4,4.7!选择长度在2.4-4.7长度单位内的线段LESIZE,ALL,,,5ALLSEL,ALLAMESH,ALL!划分网格FINISH!施加边界条件/SOLU!工况1对称集中竖向荷载DL,3,,UX$DL,3,,UY$DL,3,,UZDL,16,,UY!对线进行约束!F,662,FY,-225500!F,57,FY,-225500!加上集中荷载节点力F,662,FY,-225500!工况2反对称集中竖向荷载F,57,FY,225500SOLVEFINISH/POST1PATH,STRESS1,7!由于这个箱梁的特点，使得选路径时重复了一段，在应力图中有反映PPA TH,1,,-4.75,0.955,20$PPATH,2,,-2.35,0.9 55,20PPA TH,3,,-2.35,-1.165,20$PPATH,4,,2.35,-1.1 65,20PPA TH,5,,2.35,0.955,20$PPATH,6,,4.75,0.955 ,20PPA TH,7,,-2.35,0.955,20!七个路径点PDEF,STRESS1,S,ZPLPAGM,STRESS1,5.0!跨中截面正应力云图PLPATH,STRESS1!跨中截面正应力曲线PDEF,STRESS1,S,XYPLPAGM,STRESS1,5.0!跨中截面剪应力云图PLPATH,STRESS1!跨中截面剪应力曲线PATH,STRESS2,7PPA TH,1,,-4.75,0.955,10$PPATH,2,,-2.35,0.9 55,10PPA TH,3,,-2.35,-1.165,10$PPATH,4,,2.35,-1.1 65,10PPA TH,5,,2.35,0.955,10$PPATH,6,,4.75,0.955 ,10PPA TH,7,,-2.35,0.955,10PDEF,STRESS2,S,ZPLPAGM,STRESS2,5.0!1/4跨截面正应力云图PLPATH,STRESS2!1/4跨截面正应力曲线PDEF,STRESS2,S,XYPLPAGM,STRESS2,5.0!1/4跨截面剪应力云图PLPATH,STRESS2!1/4跨截面剪应力曲线。

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 IncrementsX,Y,Z Offsets 输入2.25,1.25,.75 点击ApplyXY，YZ，ZX Angles输入0，－90点击OK。

创建圆柱体Main Menu：Preprocessor>Create>Cylinder> Solid CylinderRadius输入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 -> -V olumes-Block -> By 2 corners & Z在创建实体块的参数表中输入下列数值:WP X = 0WP Y = 1Width = 1.5Height = 1.75Depth = 0.75 OK Toolbar: SA VE_DB3. 偏移工作平面到轴瓦支架的前表面Utility Menu: WorkPlane -> Offset WP to -> Keypoints +在刚刚创建的实体块的左上角拾取关键点OK Toolbar: 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.75 2). OK Toolbar: SAVE_DB5. 在轴承孔的位置创建圆柱体为布尔操作生成轴孔做准备Main Menu: Preprocessor -> Modeling-Create -> V olume-Cylinder -> Solid Cylinder +1.) 输入下列参数:WP X = 0WP Y = 0Radius = 1Depth = -0.18752.) 拾取Apply3.) 输入下列参数:WP X = 0 WP Y = 0 Radius = 0.85 Depth = -24.)拾取OK6．从轴瓦支架“减”去圆柱体形成轴孔.Main Menu: Preprocessor -> Modeling-Operate -> Subtract -> V olumes 1). 拾取构成轴瓦支架的两个体，作为布尔“减”操作的母体。

finish/clear/FILNAM,riser FSI/TITLE,the Response Analysis of 3D FSI/PREP7!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!! some parameters for riser geometry !!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! thick=0.014ro=0.5ri=ro-thickr1=(ro+ri)/2thicks=0.015rs=ro-thicksr2=(ro+rs)/2w=100*r1 !half width of fluid range depth=1400h2=20 !splash heighth1=depth-h2 !down length*dim,hh,array,7 ! coordinate expressionhh(1)=-1*depthhh(2)=-1*h2hh(3)=0le=2n=16 !组合的圆环数目esize1=2*3.14*r1/nesize2=2*3.14*r1/16dens=1030 ! DENSITY of sea water (kg/m**3)visc=1.3E-6vmax=5.50alphad=0.54964betad=1.0993FZ=8.0e8!!!############### Caculation Parameters ################!!! solidtime=0.005fluidtime=0.005fsitime=0.005fsitimetotal=1.10!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!! Element type and real contants !!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! et,1,142keyopt,1,4,1r,1,1et,2,63R,2,ThickR,3,Thicksmp,VISC,1,viscmp,dens,1,densmp,ex,2,2.1e11mp,nuxy,2,0.3mp,dens,2,7800!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!creat the modal and mesh the two type !!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!########## fluid part #################!!! k,10000,0,0,hh(1)kwpave,10000cyl4,0,0,w,0,r1,90,depthcyl4,0,0,w,90,r1,180,depth cyl4,0,0,w,180,r1,270,depth cyl4,0,0,w,270,r1,360,depth vglue,alllsel,s,radius,,r1lsel,a,radius,,wcm,arc,linelsel,s,loc,z,-depth/2cm,midline,linelsel,s,loc,z,-depth/2lsel,r,loc,x,-4*r1,4*r1 lsel,r,loc,y,-4*r1,4*r1 cm,midlinein,linelsel,s,line,,arclsel,invelsel,r,loc,z,0cm,linetop,linelsel,s,line,,1cm,linedownx,linelsel,s,line,,3lsel,a,line,,15,27,12cm,linedowny,linewpstyl,defa WPSTYLE,,,,,,,,0allscm,linetype1,linecm,areatype1,areacm,volutype1,voluasel,s,loc,x,-1.1*r1,1.1*r1 asel,r,loc,y,-1.1*r1,1.1*r1cm,areafsifluid,areaallssavelesize,arc,,,4lesize,linetop,16*esize2,,,1/30 lesize,linedownx,16*esize2,,,1/30 lesize,linedowny,16*esize2,,,1/30 lesize,midline,depth/400allsmshk,1type,1real,1mat,1vmesh,allallssave!!######### sloid part #################!!! k,5001,100*r1,100*r1,hh(1)k,5002,100*r1,100*r1,hh(2)k,5003,100*r1,100*r1,hh(3)*do,i,1,2l,5000+i,5000+i+1*enddokwpave,10000circle,10000,r1,,,,nlsel,s,line,,linetype1lsel,invelsel,r,loc,z,hh(1)Adrag,all,,,,,,13lsel,s,line,,linetype1 lsel,invelsel,r,loc,z,hh(2) Adrag,all,,,,,,14lsel,s,line,,13,14 ldele,all,,1kdele,10000,,kdele,5001,5003,1 wpstyl,defa WPSTYLE,,,,,,,,0allsasel,s,area,,areatype1 asel,invecm,areatype2,areacm,areafsisolid,area lsla,scm,linetype2,lineasel,s,area,,areatype2 asel,r,loc,z,hh(2),hh(3) cm,areatype2up,area asel,s,area,,areatype2asel,r,loc,z,hh(1),hh(2)cm,areatype2down,areaalls*do,i,1,3,1lsel,s,line,,linetype2lsel,r,loc,z,hh(i)lesize,all,,,1*enddo*do,i,1,2lsel,s,line,,linetype2lsel,r,loc,z,hh(i)+r1,hh(i+1)-r1 lesize,all,depth/400*enddoallstype,2real,2mat,2mshk,1amesh,areatype2downallstype,2real,3mat,2mshk,1amesh,areatype2upalls/VIEW,1,1,1,1/ANG,1/REP,FASTsave!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! create boundary conditions !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!######### cm each area and nodes #################!!! asel,s,loc,z,0cm,ftop,areaasel,s,loc,z,-depthcm,fdown,areaasel,s,area,,35cm,finright,areaasel,s,area,,31cm,finleft,areaasel,s,area,,27cm,foutright,areaasel,s,area,,3cm,foutleft,area!!######### fluid side #################!!! allsasel,s,area,,finleftasel,a,area,,finrightnsla,s,1cm,nodein,noded,all,ux,0d,all,uy,0d,all,uz,0d,all,vy,0d,all,vz,0*do,i,0,399nsel,s,node,,nodeinnsel,r,loc,z,-i*(depth/400) d,all,vx,-0.0137*i+vmax*enddoasel,s,area,,foutleft asel,a,area,,foutright nsla,s,1cm,nodeout,noded,all,ux,0d,all,uy,0d,all,uz,0d,all,vy,0d,all,vz,0*do,i,0,399nsel,s,node,,nodeoutnsel,r,loc,z,-i*(depth/400) d,all,vx,-0.0137*i+vmax*enddoasel,s,area,,fdownnsla,s,1d,all,vz,0asel,s,area,,ftopnsla,sd,all,pres,0allssave!!######### solid side #################!!! lsel,s,line,,linetype2lsel,r,loc,z,-1*depthnsll,s,1D,all,rotz,0,,,,UX,UY,UZ!d,all,all,0lsel,s,line,,linetype2nsll,s,1nsel,r,loc,z,0D,all,rotz,0,,,,UX,UYalls!!######### FSI interfaces #################!!! !structural sideallsasel,s,area,,areafsisolidnsla,s,1sf,all,fsin,1!fluid sideallsasel,s,area,,areafsifluidnsla,s,1sf,all,fsin,1savefinish/soluallsalphad,alphadbetad,betadACEL,0,0,9.8ASEL,S,AREA,,AREAFSISOLIDNSLA,S,1NSEL,R,LOC,Z,-H2F,ALL,FZ,FZ/16deltim,solidtime !固体时间设置***********************************************ALLSFLDATA1,SOLU,TRAN,1FLDATA1,SOLU,FLOW,1FLDATA1,SOLU,TEMP,0FLDATA1,SOLU,TURB,1FLDATA1,SOLU,COMP,0FLDATA1,SOLU,VOF,0FLDATA1,SOLU,SFTS,0FLDATA1,SOLU,IVSH,0FLDATA1,SOLU,SWRL,0FLDATA1,SOLU,SPEC,0FLDATA1,SOLU,ALE,1FLDATA1,SOLU,RDSF,1 !Activates the surface-to-surface radiation solutionFLDATA12,PROP,DENS,0FLDATA13,VARY,DENS,0FLDATA12,PROP,VISC,0FLDATA13,VARY,VISC,0FLDATA12,PROP,COND,0FLDATA13,VARY,COND,0FLDATA13,VARY,SPHT,0 FLDATA7,PROT,DENS,CONSTANT FLDATA8,NOMI,DENS,1050, FLDATA9,COF1,DENS,0 FLDATA10,COF2,DENS,0 FLDATA11,COF3,DENS,0 FLDATA7,PROT,VISC,CONSTANT FLDATA8,NOMI,VISC,visc FLDATA9,COF1,VISC,0 FLDATA10,COF2,VISC,0 FLDATA11,COF3,VISC,0 FLDATA12,PROP,IVIS FLDATA7,PROT,COND,CONSTANT FLDATA8,NOMI,COND,-1, FLDATA9,COF1,COND,0 FLDATA10,COF2,COND,0 FLDATA11,COF3,COND,0 FLDATA7,PROT,SPHT,CONSTANT FLDATA8,NOMI,SPHT,-1, FLDATA9,COF1,SPHT,0 FLDATA10,COF2,SPHT,0FLDATA15,PRES,REFE,101350,FLDATA16,BULK,BETA,36,FLDATA17,GAMM,COMP,1.4,FLDATA14,TEMP,NOMI,293,FLDATA14,TEMP,TTOT,293,FLDATA14,TEMP,BULK,293,TOFFSET,0,FLDATA4,TIME,STEP,fluidtime !流体时间设置******************************************* FLDATA4,TIME,ISTEP,0FLDATA4,TIME,NUMB,10,FLDATA4,TIME,TEND,1.0e06FLDATA4,TIME,GLOB,15,FLDATA4,TIME,VX,0.01,FLDATA4,TIME,VY,0.01,FLDATA4,TIME,VZ,0.01,FLDATA4,TIME,PRES,1e-004,FLDATA4,TIME,TEMP,1e-006,FLDATA4,TIME,ENKE,0.01,FLDATA4,TIME,ENDS,0.01,FLDATA4A,STEP,OVER,0,FLDATA4,TIME,OVER,0FLDATA4A,STEP,APPE,0,FLDATA4,TIME,APPE,1.0e6FLDATA4A,STEP,SUMF,10,FLDATA4,TIME,SUMF,1.0e6FLDATA4,TIME,BC,0FLDATA4,TIME,TEND,1000000,!!!############### FSI Options ################!!! FSAN,1FSOR,FLUIDFSTR,FLUID,TRANFSTR,SOLID,TRANFSIT,20,FSCO,ALL,0.01,FSOU,5 !!!fsou,1 FSI output frequency FSRE,ALL,0.8,FSIN,NONC !非保守插值****************************************************** !FSIN,consFSTI,fsitimetotal,fsitime !FSI时间设置******************************************************FSDT,fsitime, save,solve,ext,,solusolve/post1file,,rstset,lastprnsol,u,sum plnsol,u,sum plnsol,u,yplnsol,u,xplnsol,u,zfile,,rflset,lastprnsol,v,sum plnsol,v,sum plnsol,v,xplnsol,v,zplnsol,v,y/ctype,1!/ctype,2!/ctype,2,,,1/post26file,,rst!固体节点分析（-201）NSOL,2,112586,U,Y,UY-solid_201-in !固体入口的节点XVAR,1 !X-timePLVAR,2NSOL,3,109814,U,Y,UY-solid_201-out !固体出口的节点XVAR,1PLVAR,3NSOL,4,111002,U,y,UY-solid_201-cross !固体左侧的节点XVAR,1PLVAR,4NSOL,5,114170,U,y,UY-solid_201-cross !固体右侧的节点XVAR,1PLVAR,5NSOL,6,112586,U,x,Ux-solid_201-in !固体入口的节点XVAR,1 !X-timePLVAR,6NSOL,7,109814,U,x,Ux-solid_201-out !固体出口的节点XVAR,1PLVAR,7NSOL,8,111002,U,x,Ux-solid_201-cross !固体左侧的节点XVAR,1PLVAR,8NSOL,9,114170,U,x,Ux-solid_201-cross !固体右侧的节点XVAR,1PLVAR,9!!!!!fluidfile,,rfl!液体节点分析（-1001）NSOL,2,63199,v,Y,vY-fluid_1001-in !液体入口的节点XVAR,1 !X-timePLVAR,2NSOL,3,2678,v,Y,vY-fluid_1001-out !液体出口的节点XVAR,1PLVAR,3NSOL,4,2279,v,y,vY-fluid_1001-cross !液体左侧的节点XVAR,1PLVAR,4NSOL,5,35931,v,y,vY-fluid_1001-cross !液体右侧的节点XVAR,1PLVAR,5NSOL,6,63199,v,x,vx-fluid_1001-in !液体入口的节点XVAR,1 !X-timePLVAR,6NSOL,7,2678,v,x,vx-fluid_1001-out !液体出口的节点XVAR,1PLVAR,7NSOL,8,2279,v,x,vx-fluid_1001-cross !液体左侧的节点XVAR,1PLVAR,8NSOL,9,35931,v,x,vx-fluid_1001-cross !液体右侧的节点XVAR,1PLVAR,9。