搜索
!!!!!梁柱的端板连接分析（1/2模型）Finish/clear/prep7BT1=12$BT2=8$BH=200$BW=120$BL=2*BHCT1=16$CT2=12$CH=240$CW=180EPT=20$EPW=160$EPH=280$CL=2*CH+EPHBHD=21$BD=20$BHT=13$NT=18$BHND=33DISH=80$DISV1=72$DISV2=106$DISV3=40DISPP1=10$NTX1=CH+EPT+NTNTY1=DISH/2$NTZ1=CH+DISPP1+BT1+DISV3!2.创建几何模型!柱/view,1,1,1,1$/ang,1,-120,zs,1$BLC4,,,ct1,cw/2,cl$vgen,2,1,,,ch-ct1$blc4,,,ch,ct2/2,clwpoff,,,ch+dispp1$blc4,,,ch,cw/2,bt1$vgen,2,4,,,,,bh-bt1$vovlap,all$numcmp,all!梁与端板vsel,none$wpcsys,-1$wpoff,ch,,ch+dispp1$blc4,,,bl+ept,bw/2,bt1$v1=vlinqr(0,14)vgen,2,v1,,,,,bh-bt1$blc4,,,bl+ept,bt2/2,bh$wpoff,,,-dispp1$blc4,,,ept,epw/2,ephvovlap,all$numcmp,all!螺栓孔vsel,none$wpcsys,-1$wpoff,ch-2*ct1,dish/2,ntz1$wprota,,,90$cyl4,,,bhd/2,,,,blv1=vlinqr(0,14)$vgen,2,v1,,,,,disv2$vgen,2,v1,,,,,disv2+disv1$cm,hvolu,voluvsel,s,loc,x,ch-ct1,ch+ept$cm,vls,volu$vsel,allvsbv,vls,hvolu$vsel,all$numcmp,all$vplot!切分柱（为了划分映射网格）vsel,s,loc,x,0,ch$wpcsys,-1$wpoff,,,ch-ct1$vsbw,all$wpoff,,,ntz1-ch+ct1$vsbw,allwpoff,,,disv2/2$vsbw,all$wpoff,,,disv2/2$vsbw,all$wpoff,,,disv1$vsbw,allwpcsys,-1$wpoff,,,ch+eph+ct1$vsbw,all$wpoff,,dish/2$wprota,,90$vsbw,allnumcmp,all!切分柱的有孔部分wpcsys,-1$kwpave,27$wprota,,143$vsbw,84$kwpave,174$wprota,,-23$vsbw,82wpcsys,-1$kwpave,59$wprota,,51$vsbw,81$kwpave,179$wprota,,-4$vsbw,85wpcsys,-1$kwpave,179$wprota,,133$vsbw,88$kwpave,69$wprota,,7$vsbw,83wpcsys,-1$kwpave,174$wprota,,60$vsbw,86$kwpave,31$wprota,,-32$vsbw,87wpcsys,-1$kwpave,35$wprota,,152$vsbw,92$kwpave,222$wprota,,-32$vsbw,91wpcsys,-1$kwpave,54$wprota,,40$vsbw,89$kwpave,223$wprota,,9$vsbw,90vsel,all$numcmp,all!切分梁vsel,s,loc,x,ch,2*bl$wpcsys,-1$wpoff,,,ntz1$vsbw,all$qxd=(bw-dish)/2wpoff,,,qxd$vsbw,all$wpoff,,,-2*qxd$vsbw,all$wpoff,,,disv2$vsbw,allwpoff,,,qxd$vsbw,all$wpoff,,,qxd$vsbw,all$wpoff,,,disv1-2*qxd$vsbw,all!???wpoff,,,qxd$vsbw,all$wpoff,,,qxd$vsbw,all$wpoff,,dish/2$wprota,,90$vsbw,allwpoff,,,-qxd$vsbw,all$wpoff,,,2*qxd$vsbw,all!切分端板vsel,s,loc,x,ch,ch+ept$wpcsys,-1$wpoff,,,ch+dispp1$vsbw,allwpoff,,,bt1$vsbw,all$wpoff,,,bh-2*bt1$vsbw,all$wpoff,,,bt1$vsbw,allwpoff,,bt2/2$wprota,,90$vsbw,all$numcmp,all!再切分端板上有孔的部分lsel,s,radius,,bhd/2$asll,s$vsla,s$lsel,all$asel,all$vsel,r,loc,x,ch,ch+eptwpcsys,-1$kwpave,438$wprota,,45$vsbw,all$kwpave,440$vsbw,all!????kwpave,444$vsbw,all$wpcsys,-1$kwpave,432$wprota,,-45$vsbw,allkwpave,436$vsbw,all$kwpave,416$vsbw,all$wpcsys,-1$allsel,allnumcmp,all!以上命令流没有问题!!!!!!!!定义组件——柱、梁与端板vsel,s,loc,x,0,ch$cm,vcolu,voluvsel,s,loc,x,ch,ch+ept$cm,vep,voluvsel,s,loc,x,ch+ept,2*bl$cm,vbeam,volu!螺栓vsel,none$wpoff,ntx1,nty1,ntz1$wprota,,,-90$cyl4,,,bhnd/2,,,,ntcyl4,,,bd/2,,,,nt+ept+ct1+bht$wpoff,,,nt+ept+ct1$cyl4,,,BHND/2,,,,bht$vovlap,allwprota,,90$vsbw,all$wprota,,,90$vsbw,all$cm,vbolt,volu$vgen,2,vbolt,,,,,disv2vgen,2,vbolt,,,,,disv2+disv1$cm,vbolt,volu$allsel,allnumcmp,all$wpcsys,-1!定义实体单元类型和材料常数，划分单元网格et,1,solid185$mp,ex,1,2.1e5$mp,prxy,1,0.3$tb,bkin,1$tbdata,1,235!梁和柱mp,ex,2,2.2e5$mp,prxy,2,0.25$tb,bkin,2$tbdata,1,400!螺栓组mp,ex,3,2.1e5$mp,prxy,3,0.28$tb,bkin,3$tbdata,1,325$R,1$R,2$R,3!端板cmsel,s,vbolt$vatt,2,2,1$mshkey,1$esize,8$vmesh,allcmsel,s,vep$vatt,3,3,1$esize,10$vmesh,allcmsel,s,vbeam$vatt,1,1,1$lsel,s,length,,bl$lesize,all,50$lsel,allesize,10$vmesh,allcmsel,s,vcolu$lsel,s,length,,ch-ct1$lesize,all,30$lsel,s,length,,ch-2*ct1lesize,all,40$lsel,all$esize,10$vmesh,all$allsel,all!定义接触单元和目标单元et,2,conta174$et,3,targe170$et,4,conta174$et,5,conta174keyopt,4,12,3$keyopt,5,5,3dfkn=1.0$dfto=0.01!创建接触对!端板与立柱之间（标准接触）r,4,,,dfkn,dfto$cmsel,s,vep$nslv,s,1$nsel,r,loc,x,ch$real,4$type,2$esurfcmsel,s,vcolu$nslv,s,1$nsel,r,loc,x,ch$nsel,r,loc,z,ch-ct1,ch+eph+ct1type,3$esurf!螺栓接触，采用循环定义各个螺栓的接触对ntz2=ntz1+disv2$ntz3=ntz2+disv1*do,ibolt,1,3$ntzi=ntz%ibolt%mni1=2+3*ibolt$mni2=3+3*ibolt$mni3=4+3*ibolt!螺母与端板之间（标准接触）r,mni1,,,dfkn,dfto$lsel,s,radius,,bhnd/2$lsel,r,loc,z,ntzi-bhnd/2,ntzi+bhnd/2lsel,r,loc,x,ntx1-nt$asll,s$asel,r,loc,x,ntx1-nt$lsel,allnsla,s,1$real,mni1$type,2$esurflsel,s,radius,,bhd/2$lsel,r,loc,z,ntzi-bhd/2,ntzi+bhd/2$lsel,r,loc,x,ntx1-ntasll,s$asel,r,loc,x,ntx1-nt$lsel,all$nsla,s,1$type,3$esurf!栓杆与孔壁之间（标准接触）r,mni2,,,dfkn,dfto$lsel,s,radius,,bd/2$lsel,r,loc,z,mtzi-bd/2,ntzi+bd/2asll,s$asel,r,loc,x,ch-ct1,ch+ept$asel,u,loc,x,ch-ct1$asel,u,loc,x,ch+eptlsel,all$nsla,s,1$real,mni2$type,5$esurflsel,s,radius,,bhd/2$lsel,r,loc,z,ntzi-bhd/2,ntzi+bhd/2$asll,s$lsel,allasel,u,loc,x,ch$asel,u,loc,x,ch+ept$asel,u,loc,x,ch-ct1$nsla,s,1$type,3$esurf!栓头与立柱之间（绑定接触）r,mni3,,,dfkn,dfto$lsel,s,radius,,bhnd/2$lsel,r,loc,z,ntzi-bhnd/2,ntzi+bhnd/2lsel,r,loc,x,ch-ct1$asll,s$asel,r,loc,x,ch-ct1$lsel,allnsla,s,1$real,mni3$type,4$esurflsel,s,radius,,bhd/2$lsel,r,loc,z,ntzi-bhd/2,ntzi+bhd/2$lsel,r,loc,x,ch-ct1asll,s$asel,r,loc,x,ch-ct1$lsel,all$nsla,s,1$type,3$esurf$*enddo!施加载荷和约束allsel,allasel,s,loc,z,0$asel,a,loc,z,2*ch+eph$da,all,allasel,s,loc,y,0$da,all,symmasel,s,loc,z,ch+dispp1+bh$asel,r,loc,x,ch+ept,2*bl$sfa,all,1,pres,6.5$allsel,all!定义求解参数并求解/solu$nlgeom,on$outres,all,all$time,1$nsubst,50$pred,off$solve!进入后处理查看结果/post1$/expand,2,rect,half,,1e-6$pldisp,1cmsel,s,vbeam$eslv$plnsol,s,xcmsel,s,vbolt$e