ansys 储罐静力学分析
ansys 静力学 显式动力学
ANSYS静力学显式动力学1. 引言ANSYS是一款多功能的工程仿真软件，广泛应用于不同行业的产品设计、分析和优化中。
其中，静力学和显式动力学是ANSYS的两个重要模块，本文将对这两个模块进行全面、详细、完整且深入的探讨。
2. 静力学2.1 概述静力学是研究物体在静止状态下受力平衡的学科。
通过静力学分析，可以确定物体的受力情况、结构的稳定性以及构件的强度等信息。
2.2 ANSYS中的静力学分析ANSYS中的静力学分析模块可以通过建立几何模型、定义材料和边界条件来进行分析。
在分析过程中，可以考虑不同的加载情况，如静力加载和重力加载。
2.3 静力学分析的步骤静力学分析通常包括以下步骤： 1. 建立或导入模型：使用ANSYS的建模工具创建几何模型或导入现有模型。
2. 定义材料和属性：为模型定义材料特性和材料属性。
3. 定义边界条件：为模型的边界定义约束和加载条件。
4. 网格划分：将模型划分为离散的网格单元。
5. 求解分析：通过求解静力学方程，得到模型的受力状态。
6. 后处理：分析结果的可视化和数据输出。
3. 显式动力学3.1 概述显式动力学是一种研究物体在动力加载作用下的运动和响应的学科。
与静力学不同，显式动力学考虑了时间因素，可以模拟和预测物体在瞬态加载情况下的动态响应。
3.2 ANSYS中的显式动力学分析ANSYS中的显式动力学分析模块可以模拟各种动力加载条件下的物体运动和响应。
该模块可以用于模拟撞击、爆炸、碰撞、结构破坏等情况，并可以为工程师提供重要的设计参考信息。
3.3 显式动力学分析的步骤显式动力学分析通常包括以下步骤： 1. 建立或导入模型：与静力学分析相同，需要建立或导入模型。
2. 定义材料和属性：为模型定义材料特性和材料属性，以便模拟加载情况下的材料响应。
3. 定义边界条件：为模型的边界定义约束和加载条件，包括初始速度和力。
4. 网格划分：将模型划分为离散的网格单元。
5. 求解分析：通过求解显式动力学方程，得到模型在不同时间步长下的运动和响应。
大型固定式储油罐壁板静载荷下的有限元分析
大型固定式储油罐壁板静载荷下的有限元分析摘要：本文采用三维有限元法，使用ANSYS软件对某20000m3的固定式储油罐壁板在静载荷下的合应力及合位移进行了分析和模拟，对壁板应力与位移成因进行了简要分析，指出了储罐罐壁可能存在问题的重点部位，为设计、施工提供了参考。
关键词：大型储油罐、有限元方法、静载荷、应力分析引言目前大型储油罐罐壁大多采用多层不等厚壁板，在施工过程中容易发生应力集中和尺寸变形。
采用常规计算方式对其进行应力分析与强度校核，很难快速直观的得出结论。
特别是在施工及检验监理现场，对于油罐存在的应力和变形，大多数施工人员都只是简单的采用在经验公式或者简单的手工计算：不仅缺乏详细的有效数据，更谈不上进行快速的模拟与分析。
针对此类大型薄壳结构，国际上已经广泛的采用有限元法进行计算机辅助分析，本文所采用的ANSYS有限元分析软件，是一款已在工程领域大量应用的成熟软件。
1.1 油罐尺寸：罐底板外径Φ40000mm；罐内径39700mm；罐壁高度17452mm 罐顶高度23235.5mm罐壁板采用不等厚焊接：底圈6300×2000×182圈6300×2000×163圈6300×2000×144圈6300×2000×125圈6300×2000×106.7圈6300×2000×88.9圈6300×1800×8 材料16MnR焊条型号J5072.1油罐有限元模型将该油罐视为不等厚薄壳结构进行分析，考虑到油罐整体结构和载荷对称性，本文选取半个罐体进行分析，确保在不影响分析结果的前提下，减少运算量，提高分析效率，便于观察罐内及罐壁的分析结果。
罐体静载荷分析的有限元模型为：罐体及罐顶按照图纸施工，采用实体单元，总单元数为101296个，总节点数为178450个；设定分析条件时取无风、雪载荷载荷的理想条件；未考虑法兰、螺栓及焊缝等影响，未对爬梯、平台等辅助设备建模，罐壁与罐底T 形接头角焊缝视为固定约束（该处应力与变形本文未做分析）。
基于ANSYS储罐的建模研究
基于ANSYS储罐的建模研究摘要: 本文介绍了液固耦合的基本概念和储罐的有限单元类型，并利用大型有限元软件ANSYS建立了储罐液固耦合的有限元模型，为进一步开展储罐结构的静力和动力性能分析奠定基础。
关键词: 液固耦合储罐ANSYS 有限元模型Abstract：This paper introduces the basic concept of solid-liquid coupling and finite element type of tank, solid-liquid coupling finite element model of tank is established by finite element software ANSYS, and it lay the foundation for further developing the static and dynamic performance analysis of tank structure.Key Words: solid-liquid couplingtankANSYSfinite element model1.引言石油是工业的血液，在国民生产生活中有着举足轻重的作用[1]。
作为石油生产加工运输的重要设备储罐的抗震性能的好坏，就决定了石油工业的安全。
储罐是由管壁、底板和储液三部分组成，受力性能较复杂，大型有限元软件ANSYS 用于分析这种复杂结构的静力、动力、线性、非线性等响应特征时具有强大优势，可以很好地反映这种结构在各种复杂因素作用下的力学特征。
本文利用大型有限元软件ANSYS建立了储罐液固耦合的有限元模型，详细介绍了整个建模过程，为进一步开展储罐的静力和动力性能分析奠定基础。
2.液固耦合的基本概念流固耦合力学的研究对象是固体在流场作用下的各种行为以及固体变形或运动对流场的影响。
其重要特征是在于两相介质之间的相互作用，固体在流体动载荷作用下产生变形或运动，而固体的变形或运动反过来影响流场，从而改变流体载荷的分布和大小。
ANSYS静力分析的简单步骤
ANSYS静力分析的简单步骤第一步，启动工作台软件，然后选择与启动DS模块弹出得界面。
第二步，导入三维模型。
根据操作步骤进行。
首先，单击“几何体”，选择“文件”，然后选择弹出窗口中的3D模型文件，如果当时catia文件格式不符，可以把三维图先转换为“.stp”的格式，即可导入。
第三步，选择零件材料：文件导入软件后，在这个时候，依次选择“几何”下的“零件”，并且在左下角的“Details of ‘Part’”中以调整零件材料属性，本次钟形壳的材料是刚。
第四步，划分网格：选择“Project”树中的“Mesh”，右键选择“Generate Mesh”即可在这一点上，你可以在左下角的“网格”对话框的细节调整网格的大小（体积元）。
第五步，添加类型分析：第一选择顶部工具栏上的“分析”按钮，添加需要的类型分析，因为我们需要做的是在这种情况下的静态分析。
所以选择结构静力。
第六步，添加固定约束：首先选择“Project”树中的“Static Structural”按钮，右键点击支持插入固定树。
这时候在左下角的“Details of ‘Fixed Support’”对话框中“Geometry”会被选中，会要求输入固定的支撑面。
在这种情况下，固定支架的类型是表面支持，确定六凹面（此时也可点击“Edge”来确定“边”）。
然后一直的按住“CTRL”键，连续选择其它几个弧面为支撑面，在点击“Apply”进行确认，第七步，添加载荷：选择“Project”树中的“结构静力”，右键选择“Insert”中的“Force”，然后在选择载荷的作用面，再次点击“Apply”按钮进行确定。
第八步，添加变形：右键点击选择“Project”树中的“Solution”，随后依次选择插入，变形，Total”，添加变形。
第九步，添加等效应变：右键单击“项目”的树，“>插入应变->解决方案->添加等效，等效应变。
第十步，添加等效应力：首先右键点击“Project”树中的“Solution—>Insert—> Stress—>Equivalent”，添加等效应力。
ansys 储罐静力学分析
一、水箱1水箱采用钢结构（Q235），内径为18.9米，高14.533米，其中箱壁高12.288米，拱顶高2.245米。
2箱壁厚度分为4段，分别为12mm，10mm，8mm，6mm。
3边缘板厚度为12mm，中幅板厚度为10mm。
4.拱顶厚度为6mm。
5.加强圈角钢100×63×8。
命令流:finish/clear/filname,cg_1/title,cg/prep7!设置单元类型et,1,shell63et,2,beam188!荷载分项系数!结构自重分项系数rg1=1.2!除结构自重外的各项永久荷载的分项系数rg2=1.27!风荷载分项系数rq1=1.4!温度作用分项系数rq2=1.4!可变荷载作用的组合系数zuhe=0.9!罐壁和底板实常数r,1,0.012r,2,0.010r,3,0.008r,4,0.006!加强圈截面参数sectype,1,beam,l,,0secoffset,user,0,0.100secdata,0.063,0.1,0.008,0.008,0,0,0,0,0,0k,1000,0,0,8!钢罐材料属性mp,ex,1,206e9mp,prxy,1,0.3mp,dens,1,7800*rg1!罐顶材料属性,将雪荷载与自重合并xueya=0.4*100/0.006zizhong=7800*rg1+xueya*rq2*zuhe mp,ex,2,206e9mp,prxy,2,0.3mp,dens,2,zizhong!罐壁材料属性mp,ex,3,206e9mp,prxy,3,0.3mp,dens,3,7800*rg1mp,alpx,3,1e-5*set,rd,8.4 !内径*set,bian,0.065 !边缘板k,1,rd,0,0k,2,rd,0,3k,3,rd,0,6k,4,rd,0,8k,5,rd,0,10.5k,6,rd,0,12.24k,7,rd+bian,0,0k,8,rd-bian,0,0k,9,0,0,0k,10,0,0,12.24*do,i,1,5,1l,i,i+1*enddo!生成罐壁几何模型ldiv,all,,,5,0arotat,all,,,,,,9,10,360,100l,7,1l,1,8l,8,9!生成底板几何模型lsel,s,loc,z,0csys,1lsel,r,loc,y,0lsel,r,loc,x,0,rd-biancsys,0lplotldiv,all,,,20,0lsel,s,loc,z,0csys,1lsel,r,loc,y,0lsel,r,loc,x,rd-bian,rdlsel,a,loc,x,rd+0.00001,rd+bian csys,0lplotldiv,all,,,1,0lsel,s,loc,z,0csys,1lsel,r,loc,y,0lsel,r,loc,x,0,rd+biancsys,0lplotarotat,all,,,,,,9,10,360,100!生成罐顶几何模型k,10000,,,-8.379k,10001,rd,,12.24k,10002,7.6423,,12.975k,10003,5.7816,,13.5517k,10004,3.8783,,13.9669k,10005,1.9463,,14.214k,10006,,,14.301l,10001,10002l,10002,10003l,10003,10004l,10004,10005l,10005,10006lsel,s,loc,z,12.24+0.0001,14.301 lplotldiv,all,,,5,0arotat,all,,,,,,9,10,360,100!边缘板网格划分asel,s,loc,z,0csys,1asel,r,loc,x,rd-bian-0.0001,rd+bian csys,0aplottype,1real,1mat,1esize,,1amesh,all!中幅板网格划分asel,s,loc,z,0csys,1asel,r,loc,x,0,rd-bian+0.00001 csys,0aplottype,1real,2mat,1esize,,1amesh,all!罐壁网格划分asel,s,loc,z,0.0001,3csys,1asel,r,loc,x,rd,rdcsys,0aplottype,1real,1mat,3esize,,1amesh,allasel,s,loc,z,3,6aplottype,1real,2mat,3amesh,allasel,s,loc,z,6,10.5aplottype,1real,3mat,3esize,,1amesh,allasel,s,loc,z,10.5,12.24aplottype,1real,4mat,3esize,,1amesh,all!加强圈网格划分lsel,s,loc,z,8.492lplotlatt,1,1,2,,1000,,1lesize,all,1lmesh,all!罐顶网格划分asel,s,loc,z,12.24+0.0001,14.301aplottype,1real,4mat,2esize,,1amesh,all!罐顶与罐壁的铰接处理csys,1*do,i,1,100,1jiaodu=(i-1)*3.6nsel,s,loc,z,12.24nsel,r,loc,y,jiaodu-0.0001,jiaodu+0.0001 cp,i,ux,allcp,i,uy,allcp,i,uz,allnsel,s,loc,z,0nummrg,nodesavefinish!求解xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx /solu!设置底板约束ksel,s,loc,z,0dk,all,uzdk,9,uxdk,9,uydk,9,uz!施加重力荷载acel,0,0,9.8!施加风荷载csys,1*afun,degasel,s,loc,z,0.0001,12.24aplotesla,replotnsle,s,nplot!提取节点编号*get,nodemax,node,,num,max*get,nodemin,node,,num,min*dim,wind,array,nodemax,1,1*do,i,nodemin,nodemax,1*if,nsel(i),eq,1,thenzn=nz(i)yn=ny(i)!风荷载体形系数*if,zn,ne,0,then*if,yn,ge,0,and,yn,lt,15,thentxxs=1.0+(0.8-1.0)/15*(yn-0)*elseif,yn,ge,15,and,yn,lt,30txxs=0.8+(0.1-0.8)/15*(yn-15)*elseif,yn,ge,30,and,yn,lt,45txxs=0.1+(-0.7-0.1)/15*(yn-30)*elseif,yn,ge,45,and,yn,lt,60txxs=-0.7+(-1.2+0.7)/15*(yn-45) *elseif,yn,ge,60,and,yn,lt,75txxs=-1.2+(-1.5+1.2)/15*(yn-60) *elseif,yn,ge,75,and,yn,lt,90txxs=-1.5+(-1.7+1.5)/15*(yn-75) *elseif,yn,ge,90,and,yn,lt,105 txxs=-1.7+(-1.2+1.7)/15*(yn-90) *elseif,yn,ge,105,and,yn,lt,120 txxs=-1.2+(-0.7+1.2)/15*(yn-105) *elseif,yn,ge,120,and,yn,lt,135 txxs=-0.7+(-0.5+0.7)/15*(yn-120) *elseif,yn,ge,135,and,yn,lt,150 txxs=-0.5+(-0.4+0.5)/15*(yn-135) *elseif,yn,ge,150,and,yn,le,180 txxs=-0.4*elseif,yn,ge,-180,and,yn,lt,-150 txxs=-0.4*elseif,yn,ge,-150,and,yn,lt,-135 txxs=-0.4+(-0.5+0.4)/15*(yn+150) *elseif,yn,ge,-135,and,yn,lt,-120 txxs=-0.5+(-0.7+0.5)/15*(yn+135) *elseif,yn,ge,-120,and,yn,lt,-105 txxs=-0.7+(-1.2+0.7)/15*(yn+120) *elseif,yn,ge,-105,and,yn,lt,-90 txxs=-1.2+(-1.7+1.2)/15*(yn+105) *elseif,yn,ge,-90,and,yn,lt,-75 txxs=-1.7+(-1.5+1.7)/15*(yn+90) *elseif,yn,ge,-75,and,yn,lt,-60 txxs=-1.5+(-1.2+1.5)/15*(yn+75) *elseif,yn,ge,-60,and,yn,lt,-45 txxs=-1.2+(-0.7+1.2)/15*(yn+60) *elseif,yn,ge,-45,and,yn,lt,-30 txxs=-0.7+(0.1+0.7)/15*(yn+45) *elseif,yn,ge,-30,and,yn,lt,-15 txxs=0.1+(0.8-0.1)/15*(yn+30)*elseif,yn,ge,-15,and,yn,lt,0txxs=0.8+(1.0-0.8)/15*(yn+15)*elsetxxs=1.0*endif!风压高度变化系数fygd=3.12*(zn/350)**0.32!风振系数fzxs=1+zn/12.24*1.47*0.47/fygd w0=0.45e3*rq1wind(i)=fzxs*txxs*fygd*w0*elsewind(i)=0*endif*elsewind(i)=0*endif*enddocsys,0asel,s,loc,z,0.0001,12.24aplotesla,replot!提取单元编号*get,ecount,elem,,count*dim,enum,,ecount*get,enum(1),elem,,num,min*do,i,2,ecountenum(i)=elnext(enum(i-1))*enddo*do,i,1,ecountpn1=wind(nelem(enum(i),1))pn2=wind(nelem(enum(i),2))pn3=wind(nelem(enum(i),3))pn4=wind(nelem(enum(i),4)) sfe,enum(i),2,pres,,pn1,pn2,pn3,pn4 *enddo/psf,pres,norm,2,0eplotcsys,0!施加油压力!将力的施加方式设置为“累计”fcum,addcsys,0asel,s,loc,z,0.0001,12.24aplotesla,replotnsle,s,nplot*get,nodemax,node,,num,max*get,nodemin,node,,num,min*dim,water,array,nodemax,1,1*do,i,nodemin,nodemax,1*if,nsel(i),eq,1,thenzn=nz(i)water(i)=1e4*(12.24-zn)*rg2*elsewater(i)=0*endif*enddoasel,s,loc,z,0.0001,12.24aplotesla,replot!提取单元编号*get,ecount,elem,,count*dim,enum,,ecount*get,enum(1),elem,,num,min*do,i,2,ecountenum(i)=elnext(enum(i-1))*enddo*do,i,1,ecountpn1=water(nelem(enum(i),1))pn2=water(nelem(enum(i),2))pn3=water(nelem(enum(i),3))pn4=water(nelem(enum(i),4)) sfe,enum(i),1,pres,,pn1,pn2,pn3,pn4 *enddo/psf,pres,norm,2,0eplotcsys,0!将力的施加方式还原为默认的“替代”fcum,repl!施加温度作用csys,1asel,s,loc,z,0.0001,12.24aplotesla,replotnsle,s,nplotcsys,0!提取节点编号*get,nodemax,node,,num,max*get,nodemin,node,,num,min*dim,wind,array,nodemax,1,1*dim,temp_out,,nodemax*dim,temp_in,,nodemax*do,i,nodemin,nodemax,1*if,nsel(i),eq,1,thentemp_out(i)=-15*rq2*zuhetemp_in(i)=15*rq2*zuhe*elsetemp_out(i)=0temp_in(i)=0*endif*enddoasel,s,loc,z,0.0001,12.24aplotesla,replot!提取单元编号*get,ecount,elem,,count*dim,enum,,ecount*get,enum(1),elem,,num,min*do,i,2,ecountenum(i)=elnext(enum(i-1))*enddo*do,i,1,ecountpn1=temp_out(nelem(enum(i),1)) pn2=temp_out(nelem(enum(i),2)) pn3=temp_out(nelem(enum(i),3)) pn4=temp_out(nelem(enum(i),4)) bfe,enum(i),temp,5,pn1,pn2,pn3,pn4 pn1=temp_in(nelem(enum(i),1))pn2=temp_in(nelem(enum(i),2))pn3=temp_in(nelem(enum(i),3))pn4=temp_in(nelem(enum(i),4)) bfe,enum(i),temp,1,pn1,pn2,pn3,pn4 *enddoallselsolvefinish!后处理/post1rsys,1plnsol,u,x,0,1.0plnsol,u,y,0,1.0plnsol,u,z,0,1.0plnsol,s,x,0,1.0plnsol,s,y,0,1.0plnsol,s,z,0,1.0etable,shellmx,smisc,4etable,shellmy,smisc,5etable,shellmxy,smisc,6etable,shelltx,smisc,1etable,shellty,smisc,2etable,shelltxy,smisc,3pletab,shellmx,noavpletab,shellmy,noavpletab,shellmxy,noavpletab,shelltx,noavpletab,shellty,noavpletab,shelltxy,noa。
压力容器ansys分析.
高压空气储气罐ANSYS 应力分析
压力容器是在冶金、化工、炼油、气体等工业生产中频繁使用，常常用来存储各类不同压力、温度、介质的气体，或被使用为干燥罐，蒸压釜、反应釜、缓冲罐、医用氧气瓶等等。
同时大部分罐都属于特种设备—压力容器，其制造和使用国家都有严格规范标准，特别是压力容器的疲劳强度和形体薄弱环节的研究对于特种设备的安全使用很重要，这里借助于ansys软件很直观精确地将其中一种压力容器—高压空气储气罐进行了疲劳分析之一—压力应力分析。
一、高压储气罐的设计条件：
①
建立几何模型
② 由于该容器形体的对称性，选择1/4 来分析：
三、加载求解
四、结果分析。
基于ANSYS WORKBENCH的刚体动力学-静力学分析
基于ANSYS Workbench的刚体动力学-静力学分析在机械系统中，大量构件处于运动状态。
在构件的运动过程中，在某些时刻，它处于最危险的工况。
那么，如何对于一个运动的机构中某个别构件进行强度分析呢？按照以往的方法，是先使用多体动力学软件例如ADAMS进行刚体动力学分析，得到铰链处的约束力，然后再在有限元软件例如ANSYS中对感兴趣的构件划分网格，并导入从ADAMS中得到的载荷，对之进行强度分析。
ANSYS提供了一套完善的解决方案，使得直接在WORKBENCH中就可以完成全过程。
其方法如下：1. 从工具箱中，拖拽一个刚体动力学模板到项目示意图中，然后按照正常步骤创建一个刚体动力学分析，施加力，力偶等，然后插入所需要的求解结果物体。
2. 在图形窗口中确定感兴趣的时间点。
3. 选择某个求解结果物体，然后在右键菜单中选择Export Motion Loads，并指定一个载荷文件名。
4. 在项目示意图中，拷贝一个rigid dynamics分析系统。
并把它用static structural 分析系统进行取代。
5.编辑static structural分析系统，压制不需要的构件，而只留下需要分析其强度刚度的构件。
6. 把该构件的刚度行为从rigid改变成flexible.7. 把网格求解器设置从ANSYS Rigid Dynamics改成ANSYS Mechanical8. 删除或者压制所有在Rigid Dynamics分析中所使用的载荷。
9.选择static structural分支，然后在其右键菜单汇总选择Insert> Motion Loads....，从而导入前面文件中的载荷。
10.删除原有的结果物体，添加新的应力，变形等物体。
11. 求解得到此时刻(t=0.49495s)构件的应力和变形。
12.返回workbench工作平面。
ansys中的静力分析
第1章 静力分析1.1 力的概念力在我们的生产和生活中随处可见，例如物体的重力、摩擦力、水的压力等，人们对力的认识从感性认识到理性认识形成力的抽象概念。
力是物体间的机械作用，这种作用可以使物体的机械运动状态或者使物体的形状和大小发生改变。
从力的定义中可以看出力是在物体间相互作用中产生的，这种作用至少是两个物体，如果没有了这种作用，力也就不存在，所以力具有物质性。
物体间相互作用的形式很多，大体分两类，一类是直接接触，例如物体间的拉力和压力；另一类是“场”的作用，例如地球引力场中重力，太阳引力场中万有引力等。
同时力有两种效应：一是力的运动效应，即力使物体的机械运动状态变化，例如静止在地面物体当用力推它时，便开始运动；二是力的变形效应，即力使物体大小和形状发生变化，例如钢筋受到横向力过大时将产生弯曲，粉笔受力过大时将变碎等。
描述力对物体的作用效应由力的三要素来决定，即力的大小、力的方向和力的作用点。
力的大小表示物体间机械作用的强弱程度，采用国际单位制，力的单位是牛顿（N ）（简称牛）或者千牛顿（kN ）（简称千牛），1kN =103N 。
力的方向是表示物体间的机械作用具有方向性，它包括方位和指向。
力的作用点表示物体间机械作用的位置。
一般说来，力的作用位置不是一个几何点而是有一定大小的一个范围，例如重力是分布在物体的整个体积上的，称体积分布力，水对池壁的压力是分布在池壁表面上的，称面分布力，同理若分布在一条直线上的力，称线分布力，当力的作用范围很小时，可以将它抽象为一个点，此点便是力的作用点，此力称为集中力。
由力的三要素知，力是矢量，记作F ，本教材中的黑体均表示矢量，可以用一有向线段表示，如图1-1所示，有向线段AB 的大小表示力的大小；有向线段AB 的指向表示力的方向；有向线段的起点或终点表示力的作用点。
图1-11.2 静力学基本原理所谓静力学基本原理是指人们在生产和生活实践中长期积累和总结出来并通过实践反复验证的具有一般规律的定理和定律。
基于ANSYS的罐式集装箱的强度分析
82随着石油化工行业的迅速发展，世界各国对液化石油气和其它液化气的需求量越来越大。
罐式集装箱是罐式集装箱运输的运输单元，主要用于运输液化气体、液态货物和粉状危险货物，具有装载量大、运输成本低等特点，应用日趋广泛[1]。
由于罐式集装箱的运输介质大多为易燃、易爆或有毒，如果发生泄漏会造成严重的后果。
所以，罐式集装箱在贮运过程中必须安全、可靠。
本文以20英尺罐式集装箱为例，建立有限元分析模型。
利用有限元软件ANSYS [2-3]计算LPG罐式集装箱在常温环境下，考虑动态惯性力影响后结构的应力水平，根据压力容器的相关规定，分析了罐式集装箱整体和结构部件在受载情况下的应力强度分布。
基于ANSYS的分析结果，得出相应的等效应力云图，验证了罐式集装箱的结构在设计上的安全性与可靠性。
通过有限元分析的计算方法，可以为罐式集装箱的设计及强度校核提供技术参考。
1 罐式集装箱主要技术参数LPG罐式集装箱的主要技术参数如表1所示。
LPG 罐式集装箱的材料参数如表2所示。
表1 LPG罐式集装箱的设计参数标准规范ASME锅炉及压力容器规范第Ⅷ卷第二册《压力容器建造另一规则》2013版框架外形尺寸（长/宽/高）/mm 6058/2438/2591容器内径/mm φ2386介质LPG（50%C 3，50%C 4）罐体容积/m324.7工作压力/MPa 1.80 工作温度/℃-40～55设计压力/MPa 1.85设计温度/℃55介质最大充装量/Kg 11040腐蚀裕量/mm 0壳体最小成形厚度/mm 10.5（筒体）/10.0（封头）主体材质SA-612（罐体）/ SPA-H（框架）表2 LPG罐式集装箱的材料参数材料牌号公称厚度/mm 弹性模量E/MPa 泊松比v抗拉强度S T /MPa 屈服强度S Y /MPa 许用应力S /MPa SA-612t≤132007720.3570345223.6SPA-H——2007720.3490355——2 计算原理及工况2.1 有限元模型的建立采用三维设计软件Pro/E建立该罐式集装箱的实体模型，在应力分析过程中对模型进行适当简化，省略了对结构的强度和相应影响很小的附属结构，如管路系统、箱体、铭牌固定架等，这样既不影响整体计算结果的准确性，又简化了计算的工作量。
ansys静力学结果导入静力学模块
ansys静力学结果导入静力学模块通过将ANSYS静力学结果导入静力学模块，可以有效地分析结构的静力学行为。
本文将介绍ANSYS静力学结果导入静力学模块的过程，并探讨这种方法在结构分析中的应用。
我们将从简单的定义和概念入手，逐步深入了解该主题。
1. ANSYS简介ANSYS是一种广泛应用于工程仿真领域的软件，提供了多种模块来模拟和分析复杂的工程问题。
其中，静力学模块专注于结构的静力学行为和应力分析。
通过将ANSYS静力学结果导入静力学模块，我们可以进一步分析和评估结构的受力情况。
2. 静力学模块概述静力学模块是ANSYS中的一个重要模块，用于静力学仿真和结构分析。
它使用有限元方法来模拟结构的行为，并计算结构的应力和变形等参数。
该模块提供了丰富的功能和工具，帮助我们深入了解结构的受力情况。
3. ANSYS静力学结果导入将ANSYS静力学结果导入静力学模块是一种方便且有效的方法，可以直接在静力学模块中进行进一步的分析。
以下是导入过程的简要步骤：a. 在ANSYS静力学模块中打开待分析的结构模型。
b. 导入ANSYS静力学结果文件，通常是一个*.rst或*.rstp文件。
c. 对导入的结果进行后处理，包括应力和变形的计算和分析。
d. 根据需要，可以进行静力学分析的其他操作，如载荷和边界条件的修改，进一步改进模型。
4. 静力学结果分析一旦完成ANSYS静力学结果导入静力学模块的过程，我们可以根据需要执行各种静力学分析。
以下是一些常见的分析方法：a. 应力分析：通过导入的结果数据，我们可以计算和分析结构的应力分布情况，确定应力集中区域和潜在的破坏点。
b. 变形分析：通过导入的位移和变形数据，我们可以了解结构在受力时的变形情况，包括形变量和结构整体的变形情况。
c. 疲劳分析：基于导入的应力数据，可以进行疲劳分析，评估结构在长期使用和循环加载下的疲劳性能。
d. 效应分析：通过改变载荷和边界条件，我们可以模拟并分析不同情况下的结构行为，进一步优化结构设计。
Ansys静力分析详细步骤
三.分析步骤
1.打开桌面Workbench14.5，进入窗口，点击Anslysis System 下的Static Structure选项。结果如图1
2
1 图1
2.右键Geometry选项，在快捷菜单中选择Import Geometry，然后选择Browse， 如图2所示：
3
4
5
图2
3.选择路径，导入事先准备好的Stp，Igs，x_t等模型档案。选中文件，点击打 开，如图3
图11
12.点击Force，选择需要添加力的面，选择Apply，来完成力的添加，在 Magnitude里输入300N，来添加力的大小。
29
31 32
30 图12
13.点击Solution，选择Deformation，选择Total，来定义总得变形量。 34 35
33
图13
14.点击Solution，选择Stress，选择Equivalent(von-Mises)，来定义等效应力。 37 38
36
图14
15.点击Solution，点击Solve，来计算结果，点击Total Deformation ，显示模型 总变形。
39
40
图15
16. 点击Equivalent(von-Mises) ，显示等效应力。
41 图16
11
12 13
图5
6.点击Mesh，点击Mesh Control，选择sizing，来插入网格。 15
16 14
图6
7.点击Body Sizing，选择Definition下的Element Size，输入5mm，来划分网格尺 寸。
17
18 图7
8.点击Mesh，选择Generate Mesh，系统开始划分网格，结果如下。 19 20
AnsysWorkbench静力分析详细实例
生成网格后的图形如下图所示：
6 添加分析类型：选择上方工具条中的“New Analysis”，添加所需 做 的 分 析 类 型 ， 此 例 中 要 做 的 是 静 力 分 析 ， 因 此 选 择 “ Static Structural”，如下图所示。
7 添加固定约束：如下图所示，选择“Project”树中的“Static Structural”，右键选择“Insert”中的“Fixed Support”。
12 查看等效应变图：点击“Project”树中的“Solution”下的 “Equivalent Strain”，可以查看零件等效应变。
13 查看等效应力图：点击“Project”树中的“Solution”下的 “Equivalent Stress”，可以查看零件等效应力。
从弹出窗口中选择三维模型文件，如果文件格式不符，可以把三维图 转换为“.stp”格式文件，即可导入，如下图所示。
4 选择零件材料：文件导入后界面如下图所示，这时，选择 “Geometry”下的“Part”，在左下角的“Details of ‘Part’”中 可以调整零件材oject”树中的“Mesh”，右键选择 “Generate Mesh”即可。【此时也可以在左下角的“Details of ‘Mesh’”对话框中调整划分网格的大小（“Element size”项）】。
Ansys 静力分析实例：
1 问题描述： 如图所示支架简图，支架材料为结构钢，厚度 10mm，支架左侧的两 个通孔为固定孔，顶面的开槽处受均布载荷，载荷大小为 500N/mm。
2 启动 Ansys Workbench，在界面中选择 Simulation 启动 DS 模块。
3 导入三维模型，操作步骤按下图进行，单击“Geometry”，选择“From File”。
ansys储液罐实验总结
ansys储液罐实验总结一、实验介绍本次实验使用了ANSYS软件对储液罐进行了模拟分析。
储液罐是一种常见的储存液体的设备，其安全性和稳定性对于生产和生活都有着重要的意义。
通过使用ANSYS软件，可以对储液罐进行各种力学性能和流体特性的分析，以便更好地理解其工作原理和优化设计。
二、实验目的本次实验的主要目的是通过ANSYS软件模拟分析储液罐，掌握如何使用该软件进行力学分析和流体特性分析，并了解储液罐在不同工作条件下的应力、变形、温度等参数变化情况。
三、实验步骤1.建立模型：首先需要根据储液罐的实际尺寸建立一个三维模型，并设置材料属性、边界条件等参数。
2.网格划分：将模型划分成小网格，以便进行数值计算。
3.加载荷载：根据实际工作条件设置荷载类型和大小，并将其施加在模型上。
4.求解方程组：利用ANSYS软件求解方程组，得到各种参数变化情况。
5.结果后处理：对于得到的结果进行后处理，如绘制应力云图、变形云图等，以便更直观地了解储液罐的工作情况。
四、实验结果通过使用ANSYS软件进行模拟分析，得到了储液罐在不同工作条件下的各种参数变化情况。
具体包括以下几个方面：1.应力分布：根据模拟结果可以看出，在荷载作用下，储液罐内壁和底部会受到较大的压力和拉力，而顶部则受到较小的压力。
此外，在不同荷载大小和方向下，应力分布也有所不同。
2.变形情况：储液罐在荷载作用下会发生一定程度的变形，主要表现为内壁和底部向外膨胀、顶部向内凹陷等。
同时，在不同荷载条件下，变形程度也有所差异。
3.温度分布：由于液体在储存过程中会产生热量，因此储液罐内部温度会随着时间推移而逐渐升高。
根据模拟结果可以看出，在不同时间点和位置处的温度分布情况。
4.流体特性：储液罐中的液体在受到荷载作用时会产生流动，因此需要对其流体特性进行分析。
通过ANSYS软件可以得到液体的速度、压力等参数变化情况。
五、实验结论通过使用ANSYS软件对储液罐进行模拟分析，可以得到储液罐在不同工作条件下的各种参数变化情况。
Ansys_Workbench_静力分析讲义详解
Training Manual
Surface Body Edge
Solid Body Face (Scope = Target) Solid Body Edge (Scope = Target) Surface Body Face (Scope = Target)
(Scope = Contact) Bonded, No Separation All formulations Asymmetric only 1 Not supported for solving Bonded, No Separation All formulations Asymmetric only Bonded, No Separation Bonded, No Separation All formulations Symmetry respected All formulations Asymmetric only Bonded only MPC formulation Asymmetric only
Training Manual
• Solver Controls（求解控制）：
– 两种求解方式（默认是Program Controlled）：
• 直接求解 （ANSYS中是稀疏矩阵法） • 迭代求解 （ ANSYS中是PGC（预共轭梯度法））.
– Weak springs：
• 尝试模拟得到无约束的模型
对称接触
非对称接触
4-8
Static Structural Analysis
…组件 – 实体接触
• 可以使用的五种接触类型：
Contact Type Bonded No Separation Frictionless Rough Frictional Iterations 1 1 Multiple Multiple Multiple Normal Behavior (Separation) Tangential Behavior (Sliding) No Gaps No Sliding No Gaps Sliding Allowed Gaps Allowed Sliding Allowed Gaps Allowed No Sliding Gaps Allowed Sliding Allowed
基于ANSYS理论的高压储气罐三维静态分析
基于ANSYS 的高压储气罐三维静力学分析罗旭 201220679 动力工程摘要：本文利用ANSYS 软件对储气罐进行了静力学应力、应变分析，得到了储气罐的应力、应变分布图，所得结果与理论计算值能够较好的吻合，为高压储气罐的设计提供依据。
关键词：储气罐；有限元；应力；应变1工程背景高压储气罐用以储存压缩气体，保证在正常或应急情况下，提供足量的气体。
高压储气罐具有建造结构简单、使用寿命长、对环境污染少、造价低等优点，随着工业生产的发展，储气罐在我国的应用越来越广泛。
为提高储气罐的储存效率，人们总是尽可能提高储气罐内气体的压力和容积，因此对高压储气罐的结构进行全面的计算分析，具有重要的理论意义和工程实用价值。
本文以如图1.1所示储气罐作为对象进行分析，储气罐材料选用16MnR 低合金钢，并采用正火处理。
储气罐筒体直径为2米，高4米，壁厚0.04米；封头为椭圆型封头，封头高度0.4米，壁厚0.04米，具体尺寸如图1所示。
储罐内介质压力5.7MPa ，工作温度为20℃。
图1.1 储气罐结构简图2 理论计算0.060 2 40.42.080.082.1 筒体应力计算传统的高压储气罐设计是利用薄壳理论进行计算的。
根据该理论，储气罐器壁中的应力计算公式为：σ1t=p c（D i+t e）2t e≤υ[σ]t其中：σ1t——校核设计温度下圆筒器壁中的计算应力，MPa；t e——筒壁有效厚度，m；p c——储气罐设计压力，MPa；D i——筒体内径，m；υ——焊缝系数；[σ]t——材料的许用应力，MPa。
表2.1 焊接接头系数υ值其中焊缝系数可根据表1选取。
假设储液罐的所用焊缝接为双面焊的全焊透对接焊缝，焊缝进行100%探伤，则取焊缝系数υ=1。
由于储液罐材料选用16MnR低合金钢，钢材采用正火处理，查机械设计手册表32.1-12可得，材料的σb=470MPa，σs=285MPa，钢材安全系数Ƞ=1.25。
则钢材的许用应力σs t=228MPa。
ansys中的静力分析
第1章 静力分析1.1 力的概念力在我们的生产和生活中随处可见，例如物体的重力、摩擦力、水的压力等，人们对力的认识从感性认识到理性认识形成力的抽象概念。
力是物体间的机械作用，这种作用可以使物体的机械运动状态或者使物体的形状和大小发生改变。
从力的定义中可以看出力是在物体间相互作用中产生的，这种作用至少是两个物体，如果没有了这种作用，力也就不存在，所以力具有物质性。
物体间相互作用的形式很多，大体分两类，一类是直接接触，例如物体间的拉力和压力；另一类是“场”的作用，例如地球引力场中重力，太阳引力场中万有引力等。
同时力有两种效应：一是力的运动效应，即力使物体的机械运动状态变化，例如静止在地面物体当用力推它时，便开始运动；二是力的变形效应，即力使物体大小和形状发生变化，例如钢筋受到横向力过大时将产生弯曲，粉笔受力过大时将变碎等。
描述力对物体的作用效应由力的三要素来决定，即力的大小、力的方向和力的作用点。
力的大小表示物体间机械作用的强弱程度，采用国际单位制，力的单位是牛顿（N ）（简称牛）或者千牛顿（kN ）（简称千牛），1kN =103N 。
力的方向是表示物体间的机械作用具有方向性，它包括方位和指向。
力的作用点表示物体间机械作用的位置。
一般说来，力的作用位置不是一个几何点而是有一定大小的一个范围，例如重力是分布在物体的整个体积上的，称体积分布力，水对池壁的压力是分布在池壁表面上的，称面分布力，同理若分布在一条直线上的力，称线分布力，当力的作用范围很小时，可以将它抽象为一个点，此点便是力的作用点，此力称为集中力。
由力的三要素知，力是矢量，记作F ，本教材中的黑体均表示矢量，可以用一有向线段表示，如图1-1所示，有向线段AB 的大小表示力的大小；有向线段AB 的指向表示力的方向；有向线段的起点或终点表示力的作用点。
1.2 静力学基本原理所谓静力学基本原理是指人们在生产和生活实践中长期积累和总结出来并通过实践反复验证的具有一般规律的定理和定律。
AnsysWorkbench静力分析详细实例
这时左下角的“Details of ‘Fixed Support’”对话框中“Geometry” 被选中，提示输入固定支撑面。本例中固定支撑类型是面支撑，因此 要确定图示 6 位置为“Face”，【此处也可选择“Edge”来选择“边”】
然后按住“CTRL”键，连续选择两个孔面为支撑面，按“Apply”确 认，如下图所示。
8 添加载荷：选择“Project”树中的“Static Structural”，右键 选择“Insert”中的“Force”，如下图所示。
同上一步骤，选择载荷的作用面，按“Apply”确认。
然后在下面的“Magnitude”中输入“500”。
选择下面的“Direction”按钮可以选择载荷的作用方向，如下图所 示。
生成网格后的图形如下图所示：
6 添加分析类型：选择上方工具条中的“New Analysis”，添加所需 做 的 分 析 类 型 ， 此 例 中 要 做 的 是 静 力 分 析 ， 因 此 选 择 “ Static Structural”，如下图所示。
7 添加固定约束：如下图所示，选择“Project”树中的“Static Structural”，右键选择“Insert”中的“F现方式做保护处理对用户上传分享的文档内容本身不做任何修改或编辑并不能对任何下载内容负责
Ansys 静力分析实例：
1 问题描述： 如图所示支架简图，支架材料为结构钢，厚度 10mm，支架左侧的两 个通孔为固定孔，顶面的开槽处受均布载荷，载荷大小为 500N/mm。
10 求解：右键选择“Project”树中的“Solution”选择“Solve” 求解。【或者点击上方工具条中的“Solve”】，如下图所示。
11 查看变形：点击“Project”树中的“Solution”下的“Total Deformation”，可以查看零件的变形，右下方可以查看其动态变形。
AnsysWorkbench静力分析详细实例
选择图中的任意一条，或一个面，即可改变载荷的作用方向，如图中 标记 10 位置所示。然后点击“Apply”确认。
9 添加要查看的结果： 9.1 添加变形：如下图所示，右键点击“Project”树中的“Solution”， 选择“—> Insert—>Deformation—>Total”，添加变形分析。
2 启动 Ansys Workbench，在界面中选择 Simulat步骤按下图进行，单击“Geometry”，选择“From File”。
从弹出窗口中选择三维模型文件，如果文件格式不符，可以把三维图 转换为“.stp”格式文件，即可导入，如下图所示。
8 添加载荷：选择“Project”树中的“Static Structural”，右键 选择“Insert”中的“Force”，如下图所示。
同上一步骤，选择载荷的作用面，按“Apply”确认。
然后在下面的“Magnitude”中输入“500”。
选择下面的“Direction”按钮可以选择载荷的作用方向，如下图所 示。
这时左下角的“Details of ‘Fixed Support’”对话框中“Geometry” 被选中，提示输入固定支撑面。本例中固定支撑类型是面支撑，因此 要确定图示 6 位置为“Face”，【此处也可选择“Edge”来选择“边”】
然后按住“CTRL”键，连续选择两个孔面为支撑面，按“Apply”确 认，如下图所示。
12 查看等效应变图：点击“Project”树中的“Solution”下的 “Equivalent Strain”，可以查看零件等效应变。
13 查看等效应力图：点击“Project”树中的“Solution”下的 “Equivalent Stress”，可以查看零件等效应力。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、水箱1水箱采用钢结构（Q235），内径为18.9米，高14.533米，其中箱壁高12.288米，拱顶高2.245米。
2箱壁厚度分为4段，分别为12mm，10mm，8mm，6mm。
3边缘板厚度为12mm，中幅板厚度为10mm。
4.拱顶厚度为6mm。
5.加强圈角钢100×63×8。
命令流:finish/clear/filname,cg_1/title,cg/prep7!设置单元类型et,1,shell63et,2,beam188!荷载分项系数!结构自重分项系数rg1=1.2!除结构自重外的各项永久荷载的分项系数rg2=1.27!风荷载分项系数rq1=1.4!温度作用分项系数rq2=1.4!可变荷载作用的组合系数zuhe=0.9!罐壁和底板实常数r,1,0.012r,2,0.010r,3,0.008r,4,0.006!加强圈截面参数sectype,1,beam,l,,0secoffset,user,0,0.100secdata,0.063,0.1,0.008,0.008,0,0,0,0,0,0k,1000,0,0,8!钢罐材料属性mp,ex,1,206e9mp,prxy,1,0.3mp,dens,1,7800*rg1!罐顶材料属性,将雪荷载与自重合并xueya=0.4*100/0.006zizhong=7800*rg1+xueya*rq2*zuhe mp,ex,2,206e9mp,prxy,2,0.3mp,dens,2,zizhong!罐壁材料属性mp,ex,3,206e9mp,prxy,3,0.3mp,dens,3,7800*rg1mp,alpx,3,1e-5*set,rd,8.4 !内径*set,bian,0.065 !边缘板k,1,rd,0,0k,2,rd,0,3k,3,rd,0,6k,4,rd,0,8k,5,rd,0,10.5k,6,rd,0,12.24k,7,rd+bian,0,0k,8,rd-bian,0,0k,9,0,0,0k,10,0,0,12.24*do,i,1,5,1l,i,i+1*enddo!生成罐壁几何模型ldiv,all,,,5,0arotat,all,,,,,,9,10,360,100l,7,1l,1,8l,8,9!生成底板几何模型lsel,s,loc,z,0csys,1lsel,r,loc,y,0lsel,r,loc,x,0,rd-biancsys,0lplotldiv,all,,,20,0lsel,s,loc,z,0csys,1lsel,r,loc,y,0lsel,r,loc,x,rd-bian,rdlsel,a,loc,x,rd+0.00001,rd+bian csys,0lplotldiv,all,,,1,0lsel,s,loc,z,0csys,1lsel,r,loc,y,0lsel,r,loc,x,0,rd+biancsys,0lplotarotat,all,,,,,,9,10,360,100!生成罐顶几何模型k,10000,,,-8.379k,10001,rd,,12.24k,10002,7.6423,,12.975k,10003,5.7816,,13.5517k,10004,3.8783,,13.9669k,10005,1.9463,,14.214k,10006,,,14.301l,10001,10002l,10002,10003l,10003,10004l,10004,10005l,10005,10006lsel,s,loc,z,12.24+0.0001,14.301 lplotldiv,all,,,5,0arotat,all,,,,,,9,10,360,100!边缘板网格划分asel,s,loc,z,0csys,1asel,r,loc,x,rd-bian-0.0001,rd+bian csys,0aplottype,1real,1mat,1esize,,1amesh,all!中幅板网格划分asel,s,loc,z,0csys,1asel,r,loc,x,0,rd-bian+0.00001 csys,0aplottype,1real,2mat,1esize,,1amesh,all!罐壁网格划分asel,s,loc,z,0.0001,3csys,1asel,r,loc,x,rd,rdcsys,0aplottype,1real,1mat,3esize,,1amesh,allasel,s,loc,z,3,6aplottype,1real,2mat,3amesh,allasel,s,loc,z,6,10.5aplottype,1real,3mat,3esize,,1amesh,allasel,s,loc,z,10.5,12.24aplottype,1real,4mat,3esize,,1amesh,all!加强圈网格划分lsel,s,loc,z,8.492lplotlatt,1,1,2,,1000,,1lesize,all,1lmesh,all!罐顶网格划分asel,s,loc,z,12.24+0.0001,14.301aplottype,1real,4mat,2esize,,1amesh,all!罐顶与罐壁的铰接处理csys,1*do,i,1,100,1jiaodu=(i-1)*3.6nsel,s,loc,z,12.24nsel,r,loc,y,jiaodu-0.0001,jiaodu+0.0001 cp,i,ux,allcp,i,uy,allcp,i,uz,allnsel,s,loc,z,0nummrg,nodesavefinish!求解xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx /solu!设置底板约束ksel,s,loc,z,0dk,all,uzdk,9,uxdk,9,uydk,9,uz!施加重力荷载acel,0,0,9.8!施加风荷载csys,1*afun,degasel,s,loc,z,0.0001,12.24aplotesla,replotnsle,s,nplot!提取节点编号*get,nodemax,node,,num,max*get,nodemin,node,,num,min*dim,wind,array,nodemax,1,1*do,i,nodemin,nodemax,1*if,nsel(i),eq,1,thenzn=nz(i)yn=ny(i)!风荷载体形系数*if,zn,ne,0,then*if,yn,ge,0,and,yn,lt,15,thentxxs=1.0+(0.8-1.0)/15*(yn-0)*elseif,yn,ge,15,and,yn,lt,30txxs=0.8+(0.1-0.8)/15*(yn-15)*elseif,yn,ge,30,and,yn,lt,45txxs=0.1+(-0.7-0.1)/15*(yn-30)*elseif,yn,ge,45,and,yn,lt,60txxs=-0.7+(-1.2+0.7)/15*(yn-45) *elseif,yn,ge,60,and,yn,lt,75txxs=-1.2+(-1.5+1.2)/15*(yn-60) *elseif,yn,ge,75,and,yn,lt,90txxs=-1.5+(-1.7+1.5)/15*(yn-75) *elseif,yn,ge,90,and,yn,lt,105 txxs=-1.7+(-1.2+1.7)/15*(yn-90) *elseif,yn,ge,105,and,yn,lt,120 txxs=-1.2+(-0.7+1.2)/15*(yn-105) *elseif,yn,ge,120,and,yn,lt,135 txxs=-0.7+(-0.5+0.7)/15*(yn-120) *elseif,yn,ge,135,and,yn,lt,150 txxs=-0.5+(-0.4+0.5)/15*(yn-135) *elseif,yn,ge,150,and,yn,le,180 txxs=-0.4*elseif,yn,ge,-180,and,yn,lt,-150 txxs=-0.4*elseif,yn,ge,-150,and,yn,lt,-135 txxs=-0.4+(-0.5+0.4)/15*(yn+150) *elseif,yn,ge,-135,and,yn,lt,-120 txxs=-0.5+(-0.7+0.5)/15*(yn+135) *elseif,yn,ge,-120,and,yn,lt,-105 txxs=-0.7+(-1.2+0.7)/15*(yn+120) *elseif,yn,ge,-105,and,yn,lt,-90 txxs=-1.2+(-1.7+1.2)/15*(yn+105) *elseif,yn,ge,-90,and,yn,lt,-75 txxs=-1.7+(-1.5+1.7)/15*(yn+90) *elseif,yn,ge,-75,and,yn,lt,-60 txxs=-1.5+(-1.2+1.5)/15*(yn+75) *elseif,yn,ge,-60,and,yn,lt,-45 txxs=-1.2+(-0.7+1.2)/15*(yn+60) *elseif,yn,ge,-45,and,yn,lt,-30 txxs=-0.7+(0.1+0.7)/15*(yn+45) *elseif,yn,ge,-30,and,yn,lt,-15 txxs=0.1+(0.8-0.1)/15*(yn+30)*elseif,yn,ge,-15,and,yn,lt,0txxs=0.8+(1.0-0.8)/15*(yn+15)*elsetxxs=1.0*endif!风压高度变化系数fygd=3.12*(zn/350)**0.32!风振系数fzxs=1+zn/12.24*1.47*0.47/fygd w0=0.45e3*rq1wind(i)=fzxs*txxs*fygd*w0*elsewind(i)=0*endif*elsewind(i)=0*endif*enddocsys,0asel,s,loc,z,0.0001,12.24aplotesla,replot!提取单元编号*get,ecount,elem,,count*dim,enum,,ecount*get,enum(1),elem,,num,min*do,i,2,ecountenum(i)=elnext(enum(i-1))*enddo*do,i,1,ecountpn1=wind(nelem(enum(i),1))pn2=wind(nelem(enum(i),2))pn3=wind(nelem(enum(i),3))pn4=wind(nelem(enum(i),4)) sfe,enum(i),2,pres,,pn1,pn2,pn3,pn4 *enddo/psf,pres,norm,2,0eplotcsys,0!施加油压力!将力的施加方式设置为“累计”fcum,addcsys,0asel,s,loc,z,0.0001,12.24aplotesla,replotnsle,s,nplot*get,nodemax,node,,num,max*get,nodemin,node,,num,min*dim,water,array,nodemax,1,1*do,i,nodemin,nodemax,1*if,nsel(i),eq,1,thenzn=nz(i)water(i)=1e4*(12.24-zn)*rg2*elsewater(i)=0*endif*enddoasel,s,loc,z,0.0001,12.24aplotesla,replot!提取单元编号*get,ecount,elem,,count*dim,enum,,ecount*get,enum(1),elem,,num,min*do,i,2,ecountenum(i)=elnext(enum(i-1))*enddo*do,i,1,ecountpn1=water(nelem(enum(i),1))pn2=water(nelem(enum(i),2))pn3=water(nelem(enum(i),3))pn4=water(nelem(enum(i),4)) sfe,enum(i),1,pres,,pn1,pn2,pn3,pn4 *enddo/psf,pres,norm,2,0eplotcsys,0!将力的施加方式还原为默认的“替代”fcum,repl!施加温度作用csys,1asel,s,loc,z,0.0001,12.24aplotesla,replotnsle,s,nplotcsys,0!提取节点编号*get,nodemax,node,,num,max*get,nodemin,node,,num,min*dim,wind,array,nodemax,1,1*dim,temp_out,,nodemax*dim,temp_in,,nodemax*do,i,nodemin,nodemax,1*if,nsel(i),eq,1,thentemp_out(i)=-15*rq2*zuhetemp_in(i)=15*rq2*zuhe*elsetemp_out(i)=0temp_in(i)=0*endif*enddoasel,s,loc,z,0.0001,12.24aplotesla,replot!提取单元编号*get,ecount,elem,,count*dim,enum,,ecount*get,enum(1),elem,,num,min*do,i,2,ecountenum(i)=elnext(enum(i-1))*enddo*do,i,1,ecountpn1=temp_out(nelem(enum(i),1)) pn2=temp_out(nelem(enum(i),2)) pn3=temp_out(nelem(enum(i),3)) pn4=temp_out(nelem(enum(i),4)) bfe,enum(i),temp,5,pn1,pn2,pn3,pn4 pn1=temp_in(nelem(enum(i),1))pn2=temp_in(nelem(enum(i),2))pn3=temp_in(nelem(enum(i),3))pn4=temp_in(nelem(enum(i),4)) bfe,enum(i),temp,1,pn1,pn2,pn3,pn4 *enddoallselsolvefinish!后处理/post1rsys,1plnsol,u,x,0,1.0plnsol,u,y,0,1.0plnsol,u,z,0,1.0plnsol,s,x,0,1.0plnsol,s,y,0,1.0plnsol,s,z,0,1.0etable,shellmx,smisc,4etable,shellmy,smisc,5etable,shellmxy,smisc,6etable,shelltx,smisc,1etable,shellty,smisc,2etable,shelltxy,smisc,3pletab,shellmx,noavpletab,shellmy,noavpletab,shellmxy,noavpletab,shelltx,noavpletab,shellty,noavpletab,shelltxy,noa。