求助：复合结构弹塑性分析
小弟对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,solid65ET,2,link8ET,3,shell41keyopt,3,1,2et,4,solid45MP,EX,1,18000MP,prxy,1,0.2 tb,kinh,1,1,7,tbpt,,0.0001,1.8tbpt,,0.0004,6.66tbpt,,0.0008,11.84tbpt,,0.0012,15.54tbpt,,0.0016,17.76tbpt,,0.002,18.5tbpt,,0.0033,18.5 tb,conc,1,1,9tbdata,,0.3,0.5,1.75,-1 MP,EX,2,2.1E5MP,prxy,2,0.3tb,bkin,2,1,2,1tbdata,,235,0 MP,EX,3,2.35E5MP,prxy,3,0 MP,EX,4,2.1E5MP,prxy,4,0.3r,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 /replotlgen,2,8,,,50,,,,0lgen,2,13,,,100,,,,0adrag,13,,,,,,9 adrag,14,,,,,,9/pnum,link,0/replotvsba,1,7 vsba,3,8lsel,s,loc,z,-lf/2lsel,r,loc,y,0lgen,2,all,,,,a,,,0lsel,s,loc,z,-50lsel,r,loc,y,aadrag,all,,,,,,9allsel,allvsba,2,15vsba,4,14vsba,1,13!划分混凝土单元lsel,s,loc,y,0lsel,a,loc,y,alsel,a,loc,y,hlesize,all,50, lsel,s,loc,z,-lf/2lsel,a,loc,z,-(l/2-ls)lsel,u,loc,y,0lsel,u,loc,y,alsel,u,loc,y,hlesize,all,30,lesize,alltype,1mat,1real,1esys,0 mshape,0,3dmshkey,1vmesh,all !拉伸已有的体单元，生成单向板主体部分asel,s,loc,z,-lf/2 type,1extopt,esize,lf/50,0,extopt,aclear,1extopt,attr,0,0,0mat,1real,1esys,0vext,all,,,0,0,lf asel,s,loc,z,-lf/2type,1extopt,esize,(l/2-lf/2-ls)/50,0,extopt,aclear,1extopt,attr,0,0,0mat,1real,1esys,0vext,all,,,0,0,l/2-lf/2-lsallsel,all !拉伸已有体单元，生成单向板支座部分 asel,s,loc,z,l/2-lstype,1extopt,esize,1,0,extopt,aclear,1extopt,attr,0,0,0mat,1real,1esys,0vext,all,,,0,0,25,asel,s,loc,z,l/2-ls+25type,1extopt,esize,2,0,extopt,aclear,1extopt,attr,0,0,0mat,1real,1esys,0vext,all,,,0,0,ls,asel,s,loc,z,-(l/2-ls)type,1extopt,esize,1,0,extopt,aclear,1extopt,attr,0,0,0mat,1real,1esys,0vext,all,,,0,0,-25,asel,s,loc,z,-(l/2-ls+25)type,1extopt,esize,2,0,extopt,aclear,1extopt,attr,0,0,0mat,1real,1esys,0vext,all,,,0,0,-ls,allsel,alleplot!划分受力钢筋单元lsel,s,loc,x,-200lsel,a,loc,x,-100lsel,r,loc,y,atype,2mat,2real,2esys,0lmesh,all lsel,s,loc,x,0lsel,r,loc,y,atpye,2mat,2real,3esys,0lmesh,allallsel,all!划分碳纤维单元ksel,s,loc,x,850ksel,r,loc,y,0ksel,r,loc,x,0kgen,2,all,,,-cb/2,,,,0 lstr,41,97,allsel,alllstr,41,4 adrag,221,,,,,,222asel,s,,,168lsla,s lesize,all,50,,,,,,,1tpye,3mat,3real,4esys,0mshape,1,2damesh,allallsel,allnummrg,node,,,,lownumcmp,node!生成单向板钢支座asel,s,loc,y,0asel,r,loc,z,l/2-25,l/2+25type,4extopt,esize,1,0,extopt,aclear,1extopt,attr,0,0,0mat,4real,4esys,0vext,all,,,0,-20asel,s,loc,y,0asel,r,loc,z,-(l/2-25),-(l/2+25)type,4extopt,esize,1,0,extopt,aclear,1extopt,attr,0,0,0mat,4real,4esys,0vext,all,,,0,-20allsel,allfini !进入求解过程!设置边界条件/solunsel,s,loc,x,0/godsym,symm,x, nsel,s,loc,z,l/2nsel,r,loc,y,-20/god,all,,,,,,uy, nsel,s,loc,z,-l/2nsel,r,loc,y,-20/god,all,,,,,,uy,uz,!第一载荷步antype,0 nlgeom,1nropt,full,eqslv,spar,,0,asel,s,loc,y,h/gosfa,all,1,pres,f1allsel,alloutres,all,all, time,1autots,1nsubst,10,20,5,1cnvtol,f,,0.05,2,, neqit,30,pred,on,,onesel,mat,3ekill,allallsel,allsolve !第二载荷步asel,s,loc,y,hasel,r,loc,z,-lf/2,lf/2/gosfa,all,1,pres,f2outres,all,all,time,2autots,0nsubst,30,,,1pred,-1arclen,1,10,0, esel,mat,3ealive,allallsel,allsolvefini/post1esel,mat,1set,1,last,1,plnsol,s,z,0,1set,2,last,1,plnsol,s,z,0,1fini /post26allsel,allnsel,s,loc,x,0nsel,r,loc,y,0nsel,r,loc,z,0*get,n1,node,0,num,maxnsol,2,n1,u,y,prod,3,2,,,,,,-1,1,1,xvar,3plvar,1fini
页:
[1]