小弟对ansys不是太熟练,做完硕士课题后,感觉学好ansys还是必要的,再说离开学还有一段时间,小弟就较为系统地学了下ansys,所用参考教材为:《Algor、Ansys在桥梁工程中的应用方法与实例》,《ANSYS.结构有限元高级分析方法与范例应用(第一版)》,《ANSYS.结构有限元高级分析方法与范例应用(第二版)》,《ANSYS土木工程应用实例》,《ansys高级工程有限元分析范例精选——祝效华》。这些教材,基本上学习ansys最好的教材了,但遗憾的是这些书中都有错误。如能将这些错误改正,对自己和别人都大有益处。 这些书中的错误有些地方小弟可以自己看出来,有些地方就无能为力了,在此,特向各位师长,学长,同仁请教。 《ansys高级工程有限元分析范例精选——祝效华》中复合结构弹塑性分析命令流有错误,本人始终找不出来,希望大家人尽其才,帮忙看看。 FINISH /CLEAR !参数单位,力:N
长度:m,弹性模量:N/mm 容重:N/m3 *set,h,8
! *set,b,500/2
! *set,l,1800
! *set,lf,100 ! *set,ls,50
! *set,a,20 ! *set,pi,acos(-1)
! *set,sr,pi*(8/2)**2
! *set,ro,pi*(6/2)**2/200/2 ! *set,cb,300
! *set,ch,0.167
! *set,cl,l-2*ls
! !加载参数 *set,f1,0.002 ! *set,f2,30/lf ! /PREP7 ET,1,solid65 ET,2,link8 ET,3,shell41 keyopt,3,1,2 et,4,solid45 MP,EX,1,18000 MP,prxy,1,0.2 tb,kinh,1,1,7, tbpt,,0.0001,1.8 tbpt,,0.0004,6.66 tbpt,,0.0008,11.84 tbpt,,0.0012,15.54 tbpt,,0.0016,17.76 tbpt,,0.002,18.5 tbpt,,0.0033,18.5 tb,conc,1,1,9 tbdata,,0.3,0.5,1.75,-1 MP,EX,2,2.1E5 MP,prxy,2,0.3 tb,bkin,2,1,2,1 tbdata,,235,0 MP,EX,3,2.35E5 MP,prxy,3,0 MP,EX,4,2.1E5 MP,prxy,4,0.3 r,1,2,ro,0,0,,, rmore,,,,,,, r,2,sr,, r,3,sr/2,, r,4,ch,,,,,, rmore,, block,0,-b,0,h,-lf/2,-(l/2-ls), lplot /pnum,link,1 /replot lgen,2,8,,,50,,,,0 lgen,2,13,,,100,,,,0 adrag,13,,,,,,9 adrag,14,,,,,,9 /pnum,link,0 /replot vsba,1,7 vsba,3,8 lsel,s,loc,z,-lf/2 lsel,r,loc,y,0 lgen,2,all,,,,a,,,0 lsel,s,loc,z,-50 lsel,r,loc,y,a adrag,all,,,,,,9 allsel,all vsba,2,15 vsba,4,14 vsba,1,13 !划分混凝土单元 lsel,s,loc,y,0 lsel,a,loc,y,a lsel,a,loc,y,h lesize,all,50, lsel,s,loc,z,-lf/2 lsel,a,loc,z,-(l/2-ls) lsel,u,loc,y,0 lsel,u,loc,y,a lsel,u,loc,y,h lesize,all,30, lesize,all type,1 mat,1 real,1 esys,0 mshape,0,3d mshkey,1 vmesh,all !拉伸已有的体单元,生成单向板主体部分 asel,s,loc,z,-lf/2 type,1 extopt,esize,lf/50,0, extopt,aclear,1 extopt,attr,0,0,0 mat,1 real,1 esys,0 vext,all,,,0,0,lf asel,s,loc,z,-lf/2 type,1 extopt,esize,(l/2-lf/2-ls)/50,0, extopt,aclear,1 extopt,attr,0,0,0 mat,1 real,1 esys,0 vext,all,,,0,0,l/2-lf/2-ls allsel,all !拉伸已有体单元,生成单向板支座部分 asel,s,loc,z,l/2-ls type,1 extopt,esize,1,0, extopt,aclear,1 extopt,attr,0,0,0 mat,1 real,1 esys,0 vext,all,,,0,0,25, asel,s,loc,z,l/2-ls+25 type,1 extopt,esize,2,0, extopt,aclear,1 extopt,attr,0,0,0 mat,1 real,1 esys,0 vext,all,,,0,0,ls, asel,s,loc,z,-(l/2-ls) type,1 extopt,esize,1,0, extopt,aclear,1 extopt,attr,0,0,0 mat,1 real,1 esys,0 vext,all,,,0,0,-25, asel,s,loc,z,-(l/2-ls+25) type,1 extopt,esize,2,0, extopt,aclear,1 extopt,attr,0,0,0 mat,1 real,1 esys,0 vext,all,,,0,0,-ls, allsel,all eplot !划分受力钢筋单元 lsel,s,loc,x,-200 lsel,a,loc,x,-100 lsel,r,loc,y,a type,2 mat,2 real,2 esys,0 lmesh,all lsel,s,loc,x,0 lsel,r,loc,y,a tpye,2 mat,2 real,3 esys,0 lmesh,all allsel,all !划分碳纤维单元 ksel,s,loc,x,850 ksel,r,loc,y,0 ksel,r,loc,x,0 kgen,2,all,,,-cb/2,,,,0 lstr,41,97, allsel,all lstr,41,4 adrag,221,,,,,,222 asel,s,,,168 lsla,s lesize,all,50,,,,,,,1 tpye,3 mat,3 real,4 esys,0 mshape,1,2d amesh,all allsel,all nummrg,node,,,,low numcmp,node !生成单向板钢支座 asel,s,loc,y,0 asel,r,loc,z,l/2-25,l/2+25 type,4 extopt,esize,1,0, extopt,aclear,1 extopt,attr,0,0,0 mat,4 real,4 esys,0 vext,all,,,0,-20 asel,s,loc,y,0 asel,r,loc,z,-(l/2-25),-(l/2+25) type,4 extopt,esize,1,0, extopt,aclear,1 extopt,attr,0,0,0 mat,4 real,4 esys,0 vext,all,,,0,-20 allsel,all fini !进入求解过程 !设置边界条件 /solu nsel,s,loc,x,0 /go dsym,symm,x, nsel,s,loc,z,l/2 nsel,r,loc,y,-20 /go d,all,,,,,,uy, nsel,s,loc,z,-l/2 nsel,r,loc,y,-20 /go d,all,,,,,,uy,uz, !第一载荷步 antype,0 nlgeom,1 nropt,full, eqslv,spar,,0, asel,s,loc,y,h /go sfa,all,1,pres,f1 allsel,all outres,all,all, time,1 autots,1 nsubst,10,20,5,1 cnvtol,f,,0.05,2,, neqit,30, pred,on,,on esel,mat,3 ekill,all allsel,all solve !第二载荷步 asel,s,loc,y,h asel,r,loc,z,-lf/2,lf/2 /go sfa,all,1,pres,f2 outres,all,all, time,2 autots,0 nsubst,30,,,1 pred,-1 arclen,1,10,0, esel,mat,3 ealive,all allsel,all solve fini /post1 esel,mat,1 set,1,last,1, plnsol,s,z,0,1 set,2,last,1, plnsol,s,z,0,1 fini /post26 allsel,all nsel,s,loc,x,0 nsel,r,loc,y,0 nsel,r,loc,z,0 *get,n1,node,0,num,max nsol,2,n1,u,y, prod,3,2,,,,,,-1,1,1, xvar,3 plvar,1 fini |