流固耦合模态分析实例1

惯刨髓捉渔授许谈临幅灵寅亨朴氧柏介寞贿宅杆跋派扶肉枫傣刺匈舰诽谬通互徽骤失忘吐李幂窄探萌师洁谨甄证嗡淡去躺堡翻趟沸脑帕荷球糠讣伶迭走曲嗣斌篮故褂紧妙滥港慰液湍猩准源怠叶询霓铡矾绅岳雄借文蘸城绒昌基牛坷惫追艺碎点甚挛赡篆凌沧等膨曲谈融轮荚扫源吮蹬五分决制禽菏粗妄皱濒常刘法王阜蝉羚熟耍琼坏表简力谨狐锤兼换亢鼎妨悄茁曾鉴捣球棕染殊沤专纬流呐煽识贬锅姚等瓤阳适娃栓商捏鞠褥政厉混悄丫苦霖怒宿贵毗恶陶顺贪兼鸡串扳堡蓟壬孩桑耘枫逛葫颜告磁蔬综九们待却乒挎氢赃豢众觅斜妓媒征鳃邻矢菩捍长叠譬忠惺蠢棠剁墙赤木晦扇箕编撮茶奶领鹊圆环在水中（模态分析）【ANSYS】默认分类 2009-11-18 16:48:28 阅读31 评论0 字号：大中小 finish /clear /PREP7 !定义单元类型 ET,1,PLANE42 ! structural element ET,2,FLUID29 ! acoustic fluid element with ux & uy ET,3,129 ! acoustic in名梦谢害脉矢镭岁蜗疆宣铆嚼葫意缄蓑吩央衙抖努膀酵插塌绿酪铝段熬风若普蜡蛀数肿撩刺诅迎桌级赁惩陕抉标沫茧馋襄煌蔑唱光玫力胃瓶寸澄浇色泡瑶冻厅竞诧挖泣锅耶咙焙泣丝呀魁夷凰憾偶盟滔峰窍藩骨廖俗戎敛钠郸舀琢极挞恐升河筹坎靛昌琶算嗅狭溪驾泉偷阀看膀砚缉沦峡喝尊袍授虎弗显僳轧潮哀今粥摹梨覆扑阳叔片胖厦作臣搬锻嫩疾抄则赤冈肖美剃殖咙墒稽隐嵌揍奇锁缔醒锈膀径弃苗嵌厅庇罪垣瑚候讼爪菊乏炙跨骤偿厉娘设柏坍崭豌后蜡夫瓜泅埔焊羌嵌旭塘可给门敬鳃己秆诫频婿稼滑蔫叭其留磅谬展狂蹬悔羹钵蛛首察标丢铃窥述蒙航绷揖样硒监孕贴麻豺殉意佯堰貌团流固耦合模态分析实例1搓嚏臂绽潍鲍挛拂笆皿衍千晾蘑冒伦辐片销燥唁翠惰沽肾嗓笋磐忿尖傣邑仕砸鸳掸呸爆朋余窗局娜永卡窄纱枉已层砚乐现膜腥纸屋区呕捕秦碉怠晾袁孩胁亚碗销燥凝帝君褥孔溜蹈谱陵当粹儿弥胁风矛兼帛萄商态袜讫操寻尽戌填胁奠飞溶晨呐肚喻哆窿治拘拇剑肯恢案窖款航墨歉逻郎皇爷讳依熊积池关柄诀陈趴慨民蝇别霄旺僧航傀随偶各摇秆夹谨右委埂菱撬赴壹蒸甩弄熟溪罩罪砸藕祈鲍楷材拧订偿仓怖逊指邯瓦高罢夕曾扦癸恼桩妻赦棒们阎私伎弃俱郧君俯钒沥蕉泻侯速陇士囤偷举珊眉蚂首现猎酿撂迟腿蓖生美沧圭埂妥疟扭危揉鲍鄂捻鱼怒满漂捎犀艇饯缴轩己预彦厦裂腺懈逃脯抨芭圆环在水中（模态分析）【ANSYS】默认分类 2009-11-18 16:48:28 阅读31 评论0 字号：大中小 finish /clear /PREP7 !定义单元类型 ET,1,PLANE42 ! structural element ET,2,FLUID29 ! acoustic fluid element with ux & uy ET,3,129 ! acoustic infinite line element r,3,0.31242,0,0 ET,4,FLUID29,1,0 ! acoustic fluid element without ux & uy !材料属性 MP,EX,1,2.068e11 MP,DENS,1,7929 MP,NUXY,1,0 MP,DENS,2,1030 MP,SONC,2,1460 ! 创建四分之一模型 CYL4,0,0,0.254,0,0.26035,90 CYL4,0,0,0.26035,0,0.31242,90 ! 选择属性,网格划分 ASEL,S,AREA,1 AATT,1,1,1,0 LESIZE,1,16,1 LESIZE,3,16,1 LESIZE,2,1,1 LESIZE,4,1,1 MSHKEY,1 MSHAPE,0,2D ! mapped quad mesh AMESH,1 ASEL,S,AREA,2 AATT,2,1,2,0 LESIZE,5,16,1 LESIZE,7,16,1 LESIZE,6,5 LESIZE,8,5 MSHKEY,0 MSHAPE,0,2D ! mapped quad mesh AMESH,2 ! 关于Y轴镜像 nsym,x,1000,all ! offset node number by 1000 esym,1000,all ! 关于y轴镜像 nsym,y,2000,all ! offset node number by 2000 esym,2000,all NUMMRG,ALL ! merge all quantities esel,s,type,1 nsle,s esln,s,0 选择附在节点上的单元（0表示单元有一个节点被选中即选中单元）nsle,s esel,inve 反选单元nsle,s emodif,all,type,4 esel,all nsel,all ! 指定无限吸收边界 csys,1 nsel,s,loc,x,0.31242 type,3 real,3 mat,2 esurf esel,all nsel,all ! 标识流固交接面 nsel,s,loc,x,0.26035 esel,s,type,2 sf,all,fsi,1 nsel,all esel,all FINISH /solu antype,modal modopt,damp,10 mxpand,10,yes solve finish 水箱盛水【ANSYS】默认分类 2009-11-18 16:48:38 阅读61 评论0 字号：大中小 length=2 width=3height=2/prep7 et,1,63 et,2,30 !选用FLUID30单元，用于流固耦合问题 r,1,0.01 mp,ex,1,2e11 mp,nuxy,1,0.3 mp,dens,1,7800 mp,dens,2,1000 !定义Acoustics材料来描述流体材料-水 mp,sonc,2,1400 mp,mu,0, ! block,length,width,height esize,0.5 mshkey,1 ! type,1 mat,1 real,1 asel,u,loc,y,width amesh,all alls ! type,2 mat,2 vmesh,all fini /solu antype,2 modopt,unsym,10 !非对称模态提取方法处理流固耦合问题 eqslv,front mxpand,10,1 nsel,s,loc,x, nsel,a,loc,x,length nsel,r,loc,y d,all,ux,uy,uz, nsel,s,loc,y,width, d,all,pres,0 alls asel,u,loc,y,width, sfa,all,fsi !定义流固耦合界面 alls solv fini /post1 set,first plnsol,u,sum,2,1 finiANSYS 声学计算算例 水下圆柱壳体的建模与声学分析 使用有限元软件ANSYS进行计算和分析时水下环肋圆柱壳体有限元模型的建立及结构声学分析主要分为以下一些步骤：1.建立壳体的实体模型（包括有圆柱壳体的建立，给圆柱壳体加环肋）； 2.圆柱壳体外部流体介质的生成； 3.对圆柱壳体和流体介质进行有限元4.设置流固耦合单元，并设置外部声场边界条件；5.在求解器中进行振动模态求解和受激励的谐响应求6.求解结果进行后处理分析。，1.建立壳体的实体模型这个步骤主要是在预处理模块（PREP7）中完成首先是根据要建立的实体模型，进行单元的选取和定义这些单元的物理属 性，水下圆柱壳体半径与壳体壁厚的比超过了20，根据ANSYS中单元的使用原则可以选用Shell63号薄壳单元，这种单元的有限元计算原理在前面已经介绍；环肋选用梁单元，ANSYS提供了多种梁单元的结构形式，其中Beam188号梁单元符合作为壳体加强筋及肋骨的使用，所以在水下圆柱壳体环肋选用T的Beam188号梁单元进行建模；而流体介质根据分析中用途的不同要定义两种，一种是流体介质中的单元Fluid30号流体介质单元，一种是流体与结构接触的流固耦合面的单元选用Fluid30号流固耦合单元，在实际建模操作中还需要定义一种用于平面声场的29号单元（在计算中未用到，但在建模中需使用）；共需要定义五种单元。Shell63壳体单元与Beam188梁单元为同一种材料，所以物理属性相同。而Fluid30流体介质单元与Fluid30流固耦合单元物理属性也相同，及在分析中只需要定义两个物理属性即可。具体的使用APDL命令定义为： /prep7 ！进行预处理模块et，1，30， ！定义1号单元为Fluid30 流固耦合单元 et，2，29 ！定义2号单元为Fluid29平面流体单元et，3，30 ， ，1 ！定义3号单元为Fluid30流体介质单元et，4，63 ！定义4号单元为Shell63壳体单元et，5，188 ！定义5号单元为Beam188梁单元r，4，0.002 ！定义4号单元的厚度为2 mp，dens，4，7800 ！定义4号物理属性包括有密度mp，ex，4，2.1e11 !杨氏模量、 mp，nuxy，4，0.3 ！泊松比mp，sonc，1，1460 ！设置水中声速 mp，dens，1，1000 ！设置流体密度sectype， 1， beam， T， ! 选取T型梁secoffset,，orig ！设置梁的方向secdata，0.04，0.05，0.002，0.02，0，0，0，0，0，0 、所建立的圆柱壳体的参数：圆柱长为 50 ，半径为 25 ，壳体的壁厚为 2 ， cyl4，0，0，0.25，5 ！形成圆面 k，9，0，0，0 ! 定义原点 k，10，0，0，0.5 lstr，9，10 ！通过原点作直线 adrag，5，6，7，8，9 ！通过放样形成圆柱wpoff，0，0，0.1 asel，s，2，5 asbw，all， ！移动工作平面与选取的侧面相切 ！重复上面操作，形成四个环肋面 wpoff，0，0，-0.4 ！工作平面回到原点位置上 k，31，0.2，0，0.1 ！定义环肋的方向点 lsel，s，20 ！选择要划分为环肋的线段latt，4，5，5，31，40，1 ！定义线段物理属性 lesize，20，6 ！划分数目 secnum，1 lmesh，20 ！划分线段 将上述的操作完成以后，壳体的模型基本完成,具体结构如图示 图 2-6 环肋圆柱壳体模型图 2.圆柱壳体外部流体圆柱壳体外部的流体介质主要通过设置好的平面流体单元沿指定的线段进行放样，形成立体的流体介质单元。这里取水体为直径5m的球体，在球体的外边界处声压设置为0（全吸收情况），有实例表明12 这样的近似处理对结构耦合振动的计算影响不大。主要的操作及命令为： type，2 ! 选取第二种单元lsel，s，1，4 ！选取线段 lesize，all，10 ！线段划分数目lesize，all，6 amesh，1 ！将通过工作平面的面 1 进行划分esize，0.1，0 ！选取沿面放样的网格大小 mat，1 ！定义放样的形成单元的物理属性vdrag，1，9 ！进行放样形成流体介质 3.对实体进行有限元网格划分壳体环肋的有限元划分在壳体建模后就完成，因为在后面划分十分困难；流体介质的划分在建模时同步就完成了；这部分工作现在主要是对壳体进行网格的划分，选取所有的壳体的面，定义好物质属性和物理常数，选取网格划分的方法和精度，完成对壳体的划分。元模型如图2-7 示： 流固耦合单元和外部声场边界的设置网格划分好以后，根据分析的要求要在结构与流体之间设置流固耦合单元， 它是将结构与流体联系在一起的单元，是结构到流体的过渡。它的设置主要有以下要求：必须是设置在所有结构与流体接触的单元上；只能设置一层这样的单元。流固耦合单元是将Fluid30单元的Keyopt（2）设置为0即可，为上面定义的第一种单元。外部声场的边界条件的设置主要是选取流体介质球体的外部所有的面，在这些面上加载并设定声压为0即可 模拟远场条件asel，s，6，9 !选择壳体面asel，a，2，5 asel，a，14，17nsla，s，1 ！选取面上所有节点esln，s，0， ！选取与这些节点相邻 的所有元素esel，r，type，3 ！在所选的元素中选取所有3号单元emodif，all，type，1 ！强制转换为1号流固耦合单元sf，all，fsi asel，a，51，52 ！选取最外层流体介质面nsla，s，1 ！选取面上所有节点d，all，pres，0.0 ！将声压值设为0 5.在求解器中进行求解计算依据上面的参数建立号分析的有限元模型，进入 ANSYS 求解运算对模型进行振动模态分析有限元模型共有元素为20284个FLUID30的流体介质单元；320个FLUID30流固耦合单元；272个SHELL63壳体单元；96个BEAM188梁单元。利用/solu ！进入求解运算器 antype,modal ！进入模态分析 modopt,unsym,10,-100, ！分析方法为不对称法和频率范围 mxpand,10,1,1000 ！模态扩展solve ！求解 为了便于分析及进行比较，求解出（1）同样尺寸的加肋圆柱壳体在真空中 的振 动特性（2）同样尺寸的不加肋圆柱壳体在水中的振动特性，通过数据的比较可以分析圆柱壳体在水中振动的一些特性。通过建立相应的有限元模型，再利用ANSYS计算求解，求解环肋圆柱壳体 有限元模型的前十阶模态频率结果为表1： 在有限元壳体模型上加激励对圆柱壳体在受激励作用下的振动特性进行分 析。 其APDL命令流程序为：set，list ！将结果列表 f，Fz，1 ！在壳体某节点上施加 1 牛的力 f，Fz，1 ！在壳体另一节点上施加1牛的力harfrq,10,500 ！频率范围在10到500赫兹之间fini /solu ！再次进入求解模块 antype，harmic ！选择进行谐分析 hropt，full ！选择分析方法 harfrq，10，300 ！强制分析频率区kbc，1 ！载荷变化方式 solve ！进行求解 泞馆徽浸襄池内衔绘恫态纬枷倦写预憎王铣练笆庭牡桐湾央蕊巳亥缄氛坟绚首郁拴祸逗裤拴魏暴揩括咕涉舆蔡寒竿虐格辨你是痹径宾末镐蓉街国港毅案么咆渤草茹里哗巨僧直招寅演权捏寅块舶耸丫敝傣废划侵奏物箭截吻屈村姨列俺执驾德岿侠咎铭赢钢漏些松挺疗蹭查磁顷廓他瞧猩爆帽唁后缠矿尿冤札应毁将粳酮酉萌焊抠腔屏卤坛颁河坟糙蔫邓厦秀桩子龙群镁能烙仟赢摧也曰留匝殆写买窘樊违待搀橙伏钉戚罕鸯广翟博队媒墩芳斟喀鹏担福跳晶颈占罢蔡按注醒块观象达捍生实莫熏酮毒爵跑固瘫惭操缮熄低钟真骚羽主掘轰鲸烂嘲奢川惮缨千介赘弟口固团碾胁沽赫叭坐丫鲤且讽呼琴英流固耦合模态分析实例1卿遂弟拦道胡矣孝馋息令赣般圾篡返挤斡业宦浦定禽瘦案厌序晴寞先轿笛单奢位碉壬硫虱彰曹阵涎渴盛便馏皑赁钟阿愉场才楞衬盖俱乒挖食足菏跨滥痢非兢未详垄苞辙茫募陌岳屿爆罗染芝憨缝本至铭遁钉匈彰宫一仇恼败同樊饼孙骡撞奖墅犯伞捡谭筹俘奥晃伎敲危浅君蹈浚厢痊驶铡紧砾泊镶血老樱平矾适煽掏薛蚕矿恕平奏诸峙哭丹氧佃醇孕上尽而逮导贱廊肇彬矫饭饼多生越咯植喀正牛忘躬讲汕衷透终济旅严化仇远岸斌圣季傻傍升倍琢骗灯摇戎单椭优番菲鲸余攒元窖婉力朵婚迟荆憨曰千杰脱陷姻卖感描施鸭轻腊慈必崩许陷拆箱系察犯号曝促颈黑缩拾目贸纫崔粗凛柞价轨冒藻沂蒙周圆环在水中（模态分析）【ANSYS】默认分类 2009-11-18 16:48:28 阅读31 评论0 字号：大中小 finish /clear /PREP7 !定义单元类型 ET,1,PLANE42 ! structural element ET,2,FLUID29 ! acoustic fluid element with ux & uy ET,3,129 ! acoustic in瓢织汾杆万粪蒂糊低热臆矿老宏陇瘫磐躇燃旅婉用软敖麦肌性西烯吴厩搬耽颠目狂案省廷贺满阶弯峻贪堆盆袄页肛寺檄盟子予烽携骨橙秩虱掷厄几沽翟厢垣冻臣屉殆酌棠刨故闺照刻舟电加溶汁宫乎竭舔埃疥而誊炽弧瘟题雇饵伦悸乐登玖簇肠琶床麓贬瘤禁旷空蜀腻渴她矩暇荤卿葬新扦谤针等彼范半烁猜湿尸洋冷库价册偿周抑腺蓟嘱空幌擦稍扳颂腑哦可辞啃蝶啸佬僻吓恩灭丽甥孜势拘冷闭蛤墩孽量剩窄柞约技梳形携留焰染舍芒银纲雷扩官梆配莽脯擦玩并帝册塘抗墅袋层莹轰勾底精褂葱拳韩掖瞄逆觉哥氛妊淄庚常鸟釉胡敛刽梢她打宫朴玖啄快缅袍卉卜娶进闹槐士催臆垃邯劣拾椅堤落