ANSYS混合单元计算报告书.docxVIP

ANSYS混合单元计算报告书2011.3.31算例说明:在钢架计算中,钢架可以用梁单元进行模拟,模具表面可以用壳体单元进行模拟,此时在一个计算中同时要用到梁单元和壳体单元,也即混合单元。为了掌握混合单元的使用方法,学习了一个相关的算例。算例图如图1所示，图2给出了实体模型的尺寸信息。图1 广告牌示意图图2 实体模型尺寸信息广告牌受风载荷作用，立柱为圆柱体，底部节点全部约束，材料为钢材，E=200GPa, Poisson Ratio=0.3。算例求解过程(APDL方法)：采用的单位体系是：m、N/title, multi-element in a model %定义算例标题/filname,ele_test,1 %定义算例文件名%% 定义单元类型、实常数、截面信息和材料属性/prep7Et,1,shell181 %定义单元类型Et,2,solid45Et,3,beam188R,1,0.05 % 定义第一个实常数R,2,0Sectype,1,beam,csoli %定义梁单元的截面D,selSecoffset,centSecdata,0.2,8,0,0,0,0,0,0,0,0 %圆截面半径为0.2，分布8个节点Mp,ex,1,2e11 %定义材料属性Mp,dens,1,7800Mp,mu,1,0.3Mp,nuxy,0.3%% 现在开始建立实体模型Rect,0,4,0,3 %建立承受风载荷的广告牌Block,0,4,-0.5,0,0,0.5 %建立实体，分两部分建立Block,0,4,-0.5,0,0,-0.5K,,0,-5,0K,,4,-5,0K,,5,-5,0Nummrg,kp % merge all the key pointsL,6,21 %建立立柱L,8,22Wpoff,,-0.1 %将当前工作平面沿y轴负方向移动0.1mWprot,,90 % Wprot,thxy,thyz,thxz,绕x轴选择90度Vsbw,all % 用工作平面(xc-yc)分割体积，分割后的多个相邻体积在相接处共用相接面Wpoff,,,0.3 %将当前工作平面沿z轴正方向移动0.3mVsel,s,,,3,5,2 %选择编号是3、5的体Vsbw,allAlls %选择所有的对象Nummrg,allLsel,s,,,5,7 %选择编号为5、6、7的线Lsel,a,,,21 %将编号为21的线条添加到之前的线条集合中Latt,1,2,3,,23,1 %对线条赋材料、实常数、单元类型、方向点、截面形状的属性Asel,s,,,1,18,17 %选择编号为1、18的面Aatt,1,1,1,0 %对面赋材料、实常数、单元类型、单元坐标系Vsel,allVat,1,1,2,0 %对体赋材料、实常数、单元类型、单元坐标系Alls%%通过指定线条上的节点分布，决定单元的尺寸Wpcsys,-1,0 %将工作平面移到总体坐标系下，但不改变当前视角Lsel,s,loc,x,0.01,4-0.01 %选择0x4的线条Lesi,all,,,12,1 %将选中的所有线条等分成12 份Lsel,s,loc,y,0.01,2-0.01 %选择y坐标满足0y2的线条Lesi,all,,,8,2.5 %将选中的所有线条分成8 份,相邻线单元的长度比为2.5(注意查看会发现线条4的比值不对，应该改为1/2.5)Lesi,4,,,8,1/2.5 %重新对线条4的单元进行划分Lsel,a,loc,y,-5+0.01,-0.5-0.01 %选择y坐标满足-5y-0.5的线条Lesi,all,,,10,3AllsEsize,0.05 %指定单元边的默认长度为0.05mAsel,s,loc,z,0 %对承受风载荷的广告面板划分单元网格Asel,r,loc,y,-0.1,2%面板的一块面积（深度为0.1m）插入体内Amesh,allAllsLsel，s,loc,x,0 %对立柱进行单元划分Lsel,a,loc,x,4Lsel,r,loc,z,0Lsel,r,loc,y,-0.4,-5%立柱一小段（0.1m）插入体内Lmesh,allAllsVsel,all %对体进行单元划分Vmesh,allAllsAsel,s,loc,y,0+0.01,2 %施加均布风载荷Sfa,all,1,pres,100AllsNesl,s,loc,y,-5 %立柱底端固定约束D,all,allAllsSave/solusolve问题求解结果图3给出了广告牌的有限元网格划分图图3 广告牌的有限元模型结果