第三日 练习主题:梁、壳单元;多种后处理方法<BR> <BR>1.梁构架的受力分析<BR>/prep7<BR>k,1,-90,0,60 ! 首先创建关键点1,后面三个数值为坐标值,下同<BR>k,2,90,0,60 <BR>k,3,90,0,-60<BR>k,4,-90,0,-60<BR><BR>kgen,2,all,,,,120 !将所有关键点沿y方向拷贝120单位生成另外4个点<BR><BR>k,9,0,180,0<BR>k,100,0,200,0 ! 定义方向关键点<BR>k,101,90<BR>k,102,-90<BR>l,1,5 !第一个符号为L的小写,后面为关键点号,即通过两关键点连线<BR>l,2,6<BR>l,3,7<BR>l,4,8<BR>l,5,6<BR>l,6,7<BR>l,7,8<BR>l,8,5<BR>l,9,5<BR>l,9,6<BR>l,9,7<BR>l,9,8<BR><BR>lsel,,loc,y,0,119 ! 按位置(Y)选择立柱(线)<BR>cm,lvert,line ! 定义部件Vertical lines<BR>lsel,,loc,y,120 !选择水平横梁<BR>cm,lhoriz,line ! Horizontal lines<BR>lsel,,loc,y,121,180 !选择顶梁<BR>cm,lslope,line ! Sloping lines<BR>lsel,all !选择所有线<BR><BR>et,1,188 !定义188梁单元<BR>mp,ex,1,2e6 !定义弹性模量<BR>mp,nuxy,1,0.3<BR>mp,dens,1,7800<BR>sectype,1,beam,i,beam !定义第一种截面为工字梁<BR>secdata,6.535,6.535,8.06,.465,.465,.285 !定义梁截面尺寸<BR>sectype,2,beam,i,column !定义第二种截面为工字梁<BR>secdata,12,12,12.12,.605,.605,.39 !定义梁截面尺寸<BR>sectype,3,beam,hrec,peak !定义第三种种截面为口字梁<BR>secdata,6,6,.25,.25,.25,.25 !定义梁截面尺寸<BR>save,frame,db<BR><BR>cmsel,,lslope ! Sloping lines<BR>latt,1,,1,,100,,3 ! 顶梁的截面号为3,方向关键点为100<BR>cmsel,,lvert<BR>lsel,r,loc,x,-90 ! Left vertical lines<BR>latt,1,,1,,102,,2 ! 左侧的垂直梁(线1、4)截面号为2,方向关键点为102<BR>cmsel,,lvert<BR>lsel,r,loc,x,90 ! Right vertical lines<BR>latt,1,,1,,101,,2 ! 右侧的垂直梁(线2、3)截面号为2,方向关键点为101<BR>cmsel,,lhoriz<BR>lsel,u,loc,z,-60<BR>lsel,u,loc,x,90 ! Left and front horizontal lines<BR>latt,1,,1,,1,,1 ! 左前水平梁(线5、8)截面号为1,方向关键点为1<BR>cmsel,,lhoriz<BR>lsel,u,loc,z,60<BR>lsel,u,loc,x,-90 ! Right and back horizontal lines<BR>latt,1,,1,,3,,1 ! 右后水平梁(线6、7)截面号为1,方向关键点为3<BR>lsel,all<BR>lesize,all,20 ! 所有的线分割数都为20<BR>lmesh,all !对所有的线划分网格<BR> <BR><BR>nsel,,loc,y,0 !选择y=0处的所有节点<BR>d,all,all ! 约束全部自由度Fix bottom nodes<BR>ksel,all !选择所有关键点 !<BR>fk,9,fy,-10000 !给关键点9加-10000的Y方向的集中力<BR>alls !选择所有实体及有限元节点和单元<BR>/solu !进入SOLUTION<BR>solve !求解<BR><BR>后处理<BR>打开单元形状显示,查看变形及应力结果。<BR><BR>2.利用前面的轴承座模型进行扩展后处理的讲解<BR>具体体会路径操作,单元表的使用,等值图的画法,标注、切片图等技术<BR><BR>4.梁壳实体单元混合使用分析(命令流)具体讲解<BR>问题描述:广告牌承受风载荷的模拟<BR> <BR>采用单元:<BR>SOLID45<BR>SHELL181<BR>BEAM188<BR><BR>约束为两立柱的底部节点全部约束<BR>风载全部加在壳单元的面上<BR>材料均为钢材<BR><BR>以下为该例题的命令流<BR><BR>fini<BR>/cle<BR>/title, mult-element in a model<BR>/filname,ele_test,1<BR><BR>/prep7<BR>ET,1,SHELL181 <BR>ET,2,SOLID45<BR>ET,3,BEAM188<BR>R,1,0.05 <BR>r,2,0<BR>SECTYPE,1,BEAM,CSOLID,se1 <BR>SECOFFSET, CENT <BR>SECDATA,0.2,8,0,0,0,0,0,0,0,0 <BR>mp,ex,1,2e11<BR>mp,dens,1,7800<BR>mp,mu,1,0.3<BR>mp,nuxy,1,0.3<BR><BR>rect,0,4,,3<BR>block,,4,-.5,0,,0.5<BR>block,,4,-.5,0,,-0.5<BR>k,,0,-5,0<BR>k,,4,-5,0<BR>k,,5,-5,0<BR><BR>nummrg,kp<BR><BR>l,6,21<BR>l,8,22<BR><BR>wpoff,,-0.1<BR>wprot,,90<BR>vsbw,all<BR><BR>wpoff,,,0.3<BR>vsel,s,,,3,5,2<BR>vsbw,all<BR>alls<BR><BR>nummrg,all<BR><BR>lsel,s,,,5,7<BR>lsel,a,,,21<BR>LATT,1,2,3,,23,,1<BR><BR>asel,s,,,1,18,17<BR>AATT,1,1,1,0<BR><BR>vsel,all<BR>VATT,1,1,2,0<BR>alls<BR><BR>WPCSYS,-1,0<BR>lsel,s,loc,x,0.01,4-0.01<BR>lesi,all,,,12,1<BR><BR>lsel,s,loc,y,0.01,2-0.01<BR>lesi,all,,,8,2.5<BR><BR>lsel,a,loc,y,-.5-0.01,-5+0.01<BR>lesi,all,,,10,3<BR>alls<BR><BR>lesi,4,,,8,1/2.5<BR>esize,0.05<BR><BR>asel,s,loc,z,0<BR>asel,r,loc,y,-0.1,2<BR>amesh,all<BR>alls<BR><BR>lsel,s,loc,x,0<BR>lsel,a,loc,x,4<BR>lsel,r,loc,z,0<BR>lsel,r,loc,y,-0.4,-5<BR>lmesh,all<BR>alls<BR><BR>vsel,all<BR>vmesh,all<BR>alls<BR><BR>asel,s,loc,y,0+0.01,2<BR>sfa,all,1,pres,100<BR>alls<BR><BR>nsel,s,loc,y,-5<BR>d,all,all<BR>alls<BR><BR>save<BR>/solu<BR>solve<BR>fini<BR>/post1<BR>PLNSOL,S,EQV,0,1<BR>ANCNTR,10,0.5<BR> |