第N章随机振动案例下面介绍对一个任意模型加载中国军用标准中振动试验标准所规定的功率密度谱来演示ANSYSWOKBENCH14.0机械设模块动力学分析中随机振动分析模块的基本操作过程。1.5.1案例介绍本案例主要参考了GJB150.16-1986:《中华人民共和国国家军用标准--军用设备环境试验方法--振动试验》。其部分内容如下:“本标准规定了军用设备振动试验方法,是制订军用设备技术条件或产品标准等技术文件的相应部分的基础和选则依据”。根据标准第2.3.1条规定,“作为固紧货物的设备通过陆地、海上或空中运输时都将遇到这种环境。陆上运输环境比海上或空中更为严重,而且所有海上或空中运输的前后都将包括陆上运输,因此以陆上运输来作为基本运输环境。陆上运输环境包括公路运输和铁路运输,而公路运输比铁路运输更为严重因此以公路运输来作为运输环境。公路运输的环境是一种宽带振动,它是由于车体的支撑、结构与路面平度的综合作用产生的。设备的运输一般是指从制造厂到用户以及用户之间所经受的典型环境。这些运输科分为两个阶段,公路运输和野战任务运输野战任务运输通常是由双轮拖车,2.5~10T的卡车,半拖车和(或)履带车来完成,典型举例是500KM。路面条件差,在战斗环境下将经历恶劣的路面和原始地形”。由于野战运输环境下的功率密度谱的振动更强,故笔者选用了标准中规定了第一类设备在“基本运输环境”中第98,99页“双轮拖车环境”的功率密度谱。本次通过使用军用标准中激励相对较强的功率密度谱进行分析,可以体现较为严格的环境从而更完整的展示结构在随机振动激励下的各种响应情况。1.5.2启动Workbench并建立分析项目(1)打开ANSYSWORKBENCH14.0,并双击Toolbox(工具箱)→AnalysisSystem(分析系统)→Model(模态分析),如图-1所示。(2)单击RandomVibration(随机震动)模块,按住鼠标并将其拖动到项目管理区分析项目A6Solution(分析)中。如图-2所示。这样模态分析结果会作为随机振动分析的条件,两个分析数据也会互相连接共用。操作完毕后会如图-3所示。图-1启动模态分析图-2启动随机振动分析图-3打开分析项目图-4导入模型1.5.3导入模型(1)在项目管理区双击Model(模态分析)中的A3Geometry(几何模型)项目进入DM模块。如图-3所示。(2)单击DM模块菜单栏上的File(文件)→ImportExternalGeometryFile(输入外部几何模型文件),如图-4所示。然后选择合适的模型文件→单击“打开”按钮,如图-5所示。(3)单击菜单栏上的(刷新)按钮来刷新此次操作完成对外部模型的导入。刷新模型过程可能需要等待几分钟时间。模型导入后如图-6所示。图-5打开模型文件图-6导入后的模型(4)保存项目文件。单击菜单栏中的(保存)按钮先将此次分析项目保存,出现如图-7所示的对话框。在合适的文件夹处,我们暂时命名为“1”然后单击“保存”按钮。注意:对于ANSYSWORKBENCH14.0,随着分析规模的不同可能需要几百兆到几万兆的硬盘空间,保存项目文件时请尽量选择磁盘空间较大的分区。而且有限元分析过程中将会出现大量的数据交换工作。硬盘作为计算机中数据交换速度最慢的部件,有时会成为整机的瓶颈,需要尽可能高速的硬盘。由于大部分用户的硬盘都是传统的机械硬盘,其转速一定。磁盘外圈的线速度较大从而其读写速度相对较快。常规情况下磁盘外圈的分区一般为C盘。为了减少数据交换瓶颈,笔者建议如有可能尽量将项目文件保存在C盘中。为了节约宝贵的C盘空间,完成分析后请将项目文件转移到其他分区。如有更多投资,为了提高数据硬盘整体读写效率及可靠性,还建议使用更高转速(10000转或15000转)的硬盘以及使用基于RAID技术的一些冗余方案。如果用户有足够的资金,笔者强烈建议使用更高速的固态硬盘。对于SATA3.0接口最新的固态硬盘其最高读写速度普遍能比机械硬盘快2倍或以上。而部分PCI-E4X接口的固态甚至能快近10倍,而4K小文件的随机读写速度甚至可以快近100倍。从而带来比普通机械硬盘快几倍甚至几十倍的效率提升。对于时间紧迫的有限元分析工作来说,同样投资额下增加固态硬盘投资带来的性能提升比例会大大优于提高其他部件的性能带来的效果。对于固态硬盘来说其不存在内外圈速度不同的问题,故可以保存在任意固态硬盘分区下。但是由于固态硬盘的本质原因,其在整盘磁盘空间接近饱和时读写速度会下降,故请尽量使其预留一定的空闲空间。如笔者2012年春节前后进行的某核电站空调结构抗震分析,由于电脑硬盘速度较慢,在整个连续近80小时运行的分析过程中,有70多个小时都是CPU占用率5%左右而硬盘占用率接近100%。硬盘读写效率瓶颈明显,极大的影响了电脑其他部件性能的发挥。分析完成后,最终结果文件共23G。故笔者购买了某品牌容量为120G的固态硬盘,更换后笔者的Windows7系统硬盘部分评分从5.4分提高到了7.8分,直逼系统默认硬件最高评分7.9分。取得了良好的加速效果。由于使用固态硬盘尽可能的消除了硬盘部分的性能瓶颈,使得计算机其他硬件性能更集中发挥,这会带来发热量增多问题,故请适当考虑加强散热。图-7保存项目文件图-8网格划分1.5.3网格划分(1)下面将实体模型进行离散化从而划分出有限数量的节点与单元,方便程序分析计算。回到项目管理区双击A4Model(模型),如图-8所示。(2)单击Outline(分析树)→Mesh(网格)。我们先使用自动划分网格查看质量后再考虑是否进行网格控制。单击Mesh(网格)菜单栏上的(刷新网格)。稍等几分钟后网格划分完毕,如图-9所示。图-8划分网格图-9划分后的网格(3)查看发现网格较为粗大,我们指定网格尺寸后再此划分网格。单击Detailsof“Mesh”(网格详细信息)中的ElementSize(网格尺寸)输入5,单位是毫米。如图-10所示。再次单击(刷新)按钮查看网格。同时在图-10中的Statistice(统计)中也可以看见第一次网格划分时的节点数为31129个,单元数为16631个,也相对较少。图-11显示了设定较小网格尺寸后的划分状况。再次查看Detailsof“Mesh”(网格详细信息)中的Statistice(统计)。节点数变为312169个,单元数为164609个,数量较第一次划分的网格增多了近10倍。图-10设定网格尺寸图-11第二次划分后的网格图-12第二次划分后的网格1.5.4施加约束(1)我们在四个模型底脚处设定4个固定约束来模拟设备实际固定状况。首先查看模型空间中的坐标系,鼠标指针移动到Y轴的反向,系统会自动生成一个灰色的与Y轴方向相反的指向箭头。单击其,将模型旋转到Y轴反向来方便我们选择约束面。如图-13所示。图-13旋转方向图-14施加固定约束(2)按住键盘上的Ctrl(控制)键同时分别单击模型4个底面。选择完毕后的目标面会变成绿色。→单击菜单栏中的Supports(支撑)→单击第一项FixedSupport(固定约束)如图-14所示。然后回到Detailsof“FixedSupport”(固定约束的详细信息)单击Apply(确定)按钮。如图-15所示。图-15确定约束面1.5.5求解模态(1)单击Outline(分析树)→Model(A5)→AnalysisSettings(分析设置)→在Detailsof“AnalysisSettings”(分析设置详细信息)中的MaxModestoFind(最大模态数)中设定10。来计算本模型的前10阶固有频率以及阵型。如图-16所示。注意:一般来说外界激励对模型的较低阶的固有频率影响较大。读者们如果有精力可以参看相关的建筑抗震设计规范(GB50011)以及《化工设备设计全书-塔设备》等资料。其提供各种简化后的手工计算公式。由于完整的结构计算的计算量过于巨大,手工计算基本无法完成,标准中基本都是仅仅计算结构在第一阶固有频率下的响应情况。而对于有限元分析来说一般而言计算前6阶的固有频率基本可以达到一般要求,计算前10阶固有频率基本都能满足工程计算精度。故除特殊情况,不需要计算过多的固有频率。(2)单击(求解)按钮。由于网格数量达到了16万多个,相对计算规模较大,故此求解过程需要较长时间。求解过程中会出现ANSYSWorkbenchSolutionStatus对话框,我们也可以打开系统的任务管理器监视整个求解过程。如图-17所示。图-16设定模态数图-17正在求解(3)求解完毕后我们可以查看分析信息。单击Outline(分析树)→SolutionA6(分析)→SolutionInformation(求解信息),如图-18所示。查看右边的Worksheet(数据表)。数据表最后面的ElapsedTime(共用时间)显示了模态分析所耗费的时间,共1320秒。如图-19所示。注意:在求解过程中查看Worksheet(数据表)实时渐刷新并显示求解过程的信息以及提示各种Warming(警告)和Error(错误)。虽然大部分时候分析都是成功的,但是依然会不时出现各种问题尤其是各种的Error(错误)的报告。当出现这种严重错误导致分析无法进行下去的时候就需要仔细查看Worksheet(数据表)中的各种提示查找出错原因并适当修改各种分析设置或是模型、网格尺寸来消除错误。对于有限元分析来说如果不是特别针对于某些倒角,钻孔等特征进行分析,我们可以在建立分析用的模型的时候适当删除掉原始模型中尺寸较小的倒角,棱边,槽,孔等特征。这样一方面可以减少计算规模,更重要的是可以划分出较高质量的网格从而提高分析的精度。图-18查看分析信息图-19查看所费时间1.5.6模态分析结果后处理(1)单击Outline(分析树)→SolutionA6(分析)项→单击后处理菜单栏里面的Deformation(变形)→Total(全部的)。如图-20所示。(2)由于我们在AnalysisSettings(分析设置)里面设置了MaxModestoFind(最大模态数)为10,那么我们也需要分别查看这10阶阵型的结果,所以需要连续单击10次Total(全部的)并在每一个Total(全部的)结果里面依次设定显示1到10阶阵型结果。比如,要显示第二阶阵型则现需要在Detailsof“TotalDeformation2”(第二阶阵型)→Model(模态数)里设定为2。如图-21所示。设定完毕后单击菜单栏上的(求解)按钮对结果进行求解。图-20提取结果图-21显示模态结果(3)查看模型的前10阶阵型。分别单击Outline(分析树)→SolutionA6(分析)→从“TotalDeformation”一直到“TotalDeformation10”。如图-22到图-31所示。图-22一阶阵型图-23二阶阵型图-24三阶阵型图-25四阶阵型图-26五阶阵型图-27六阶阵型图-28七阶阵型图-29八阶阵型图-30九阶阵型图-31十阶阵型(4)我们也可以查看TabularData(数据列表)来看看模型的前10阶固有频率的情况。如图-32所示。以及查看Graph(图表)来宏观上看模型固有频率的分布。如图-33所示。图-32频率列表图-33图表显示1.5.7随机振动分析-加载功率密度谱(1)我们需要分别加载X、Y、Z三个方向的功率密度谱。单击RandomVibration(B5)(随机振动)→单击PSDBaseExcitation(基础的功率密度谱激励)→PSDGAcceleration(功率密度谱G加速度激励)按钮。由于我们需要施加三个方向的功率密度谱,那么需要连续单击三次PSDGAcceleration(功率密度谱G加速度激励)按钮。如图-34所示。图-34施加G加速度激励图-35选择边界条件(2)单击Outline(分析树)→PSDGAcceleration。在Detailsof“PSDGAcceleration”(功率密度谱G加速度激励的详细信息)中单击BoundaryCondition(边界条