I型双块式ANSYS轨枕模型建立.docx

!I型双块式轨枕模型建立fini/cle /FILNAME,4GP/prep7 !定义轨排混凝土单元（SOLID45)ET,1,SOLID45MP,EX,1,2.06E11MP,PRXY,1,0.3MP,DENS,1,7800MP,ALPX,1,11.8E-6ET,2,beam4MP,EX,2,1.95E11MP,PRXY,2,0.3MP,DENS,2,7850mp,ALPX,2,1.2e-5ET,3,SOLID45MP,EX,3,3.65E10MP,PRXY,3,0.2MP,DENS,3,2500MP,ALPX,3,1.0E-5!钢筋直径D1=12*1e-3D2=10*1e-3D3=7*1e-3PI=3.1415926r,1,0.25*PI*D1*2,PI*D1*4/64,PI*D1*4/64,D1,D1 !上弦钢筋r,2,0.25*PI*D2*2,PI*D2*4/64,PI*D2*4/64,D2,D2 !下弦钢筋r,3,0.25*PI*D3*2,PI*D3*4/64,PI*D3*4/64,D3,D3 !斜拉钢筋!长宽高lgz=844*1e-3kgz=314*1e-3hgz=174*1e-3Ksx=170*1e-3 !上弦钢筋间距kxx=0.062 !下线钢筋到上弦z向距离!画钢轨!定义关键点坐标k, 1, 0.108, 0.1392, -1k, 2, 0.110, 0.1412,-1k, 3, 0.109, 0.1600,-1k, 4, 0.099, 0.1738,-1k, 5, 0.083, 0.1758,-1k, 6, 0.073, 0.1760,-1k, 7, 0.063, 0.1758,-1k, 8, 0.048, 0.1738,-1k, 9, 0.038, 0.1600,-1k, 10, 0.037, 0.1412,-1k, 11, 0.038, 0.1392,-1k, 12, 0.056, 0.1334,-1k, 13, 0.060, 0.1297,-1k, 14, 0.063, 0.1203,-1k, 15, 0.063, 0.0396,-1k, 16, 0.050, 0.0226,-1k, 17, 0.032, 0.0166,-1k, 18, 0.023, 0.0148,-1k, 19, 0.002, 0.0124,-1k, 20, -0.002, 0.0084,-1k, 21, -0.002, 0.0020,-1k, 22, 0.000, 0.000,-1k, 23, 0.073, 0.000,-1k, 24, 0.146, 0.000,-1k, 25, 0.148, 0.0020,-1k, 26, 0.148, 0.0084,-1k, 27, 0.145, 0.0124,-1k, 28, 0.123, 0.0148,-1k, 29, 0.115, 0.0166,-1k, 30, 0.097, 0.0226,-1k, 31, 0.084, 0.0396,-1k, 32, 0.084, 0.1203,-1k, 33, 0.087, 0.1297,-1k, 34, 0.091, 0.1334,-1!关键点创建曲线BSPLINE,2,3,4,BSPLINE,4,5,6,7,8,BSPLINE,8,9,10BSPLINE,10,11,12,13,BSPLINE,33,34,1,2BSPLINE,13,14,15,16,BSPLINE,30,31,32,33BSPLINE,16,17,18,19BSPLINE,27,28,29,30BSPLINE,19,20,21,22BSPLINE,24,25,26,27BSPLINE,22,23,24BSPLINE,13,33BSPLINE,16,30,!已定义的边界线生成面（封闭）AL,1,2,3,4,5,13AL,6,7,13,14AL,8,9,10,11,12,14cm,GG,AREA AGEN,2,ALL,-0.146/2,0.5,1000CMSEL,S,GGLSLAKSLLADELE,ALLLDELE,ALLKDELE,ALLALLSEL,ALLCM,GG,AREAlsel,none !建立钢筋 !上弦补点 *do,i,1,13k,100+i,-1.200+0.2*(i-1),0.034,(KGZ-KSX)/2*enddo *do,i,1,13k,200+i,-1.200+0.2*(i-1),0.034,kgz-(KGZ-KSX)/2*enddo !下弦补点 *do,i,1,12k,300+i,-1.100+0.2*(i-1),-0.055,(KGZ-KSX)/2-kxx*enddo *do,i,1,12k,400+i,-1.100+0.2*(i-1),-0.055,(KGZ-KSX)/2+kxx*enddo *do,i,1,12k,500+i,-1.100+0.2*(i-1),-0.055,kgz-(KGZ-KSX)/2-kxx*enddo *do,i,1,12k,600+i,-1.100+0.2*(i-1),-0.055,kgz-(KGZ-KSX)/2+kxx*enddo !联系 *do,j,1,2 *do,i,1,12l,100+i+100*(j-1),100+i+1+100*(j-1)*enddo*enddo cm,sxgj,line!联系lsel,none *do,j,1,4 *do,i,1,11l,300+i+100*(j-1),300+i+1+100*(j-1)*enddo*enddo cm,xxgj,line!联系 lsel,none *do,j,1,2 *do,i,1,12l,100+i,300+i+100*(j-1)l,100+i+1,300+i+100*(j-1)*enddo*enddo *do,j,1,2 *do,i,1,12l,200+i,500+i+100*(j-1)l,200+i+1,500+i+100*(j-1)*enddo*enddo cm,xlgj,line blc4,0.286,0,lgz,hgz,kgzvsymm,x,alllsel,none!切出钢筋位置!ganggui H1=1.516 H2=1.516 KJJJ=0.650 LMJJ=0.117WPCSYS,-1WPROTA,90 WPOFFS,H2/2-LMJJ vsbw,all WPOFFS,LMJJ*2 vsbw,all WPCSYS,-1WPROTA,90 WPOFFS,-H2/2+LMJJ vsbw,all WPOFFS,-LMJJ*2 vsbw,all !切出钢轨扣件位置ALLSEL,ALL WPCSYS,-1WPOFFS,kgz/2vsbw,allallsel,alllsba,all,all,keep,keep !对钢轨进行分组VSEL,S,LOC,Y,HGZ+0.001,HGZ+1CM,VGG,VOLU !mesh size=0.01cmsel,s,sxgjlatt,2,1,2LESIZE,ALL,sizelmesh,all cmsel,s,xxgjlatt,2,2,2LESIZE,ALL,sizelmesh,all cmsel,s,xlgjlatt,2,3,2LESIZE,ALL,size*1.5lmesh,all ALLSEL,ALLCMSEL,U,VGGTYPE,3MAT,3ESIZE,size*2MSHAPE,0,3D !指定划分单元形状（生成六面体单元）MSHKEY,2 !划分方式（0-自由，1-映射，2-尽可能映射）Vsweep,ALL !钢轨!复制单元allsel,allegen,4,500000,all,0,0,kjjjcmsel,s,xxgjcmsel,a,sxgjcmsel,a,xlgjlgen,4,all,0,0,kjjjcm,gj,line!建立钢轨实体hGG=0.032CMSEL,S,GGAGEN,2,ALL,-H1/2,-0.5+hgz+hgg,1+kgz/2-kjjj/2,100cm,gg1,areaCMSEL,S,GGAGEN,2,ALL,H1/2,-0.5+hgz+hgg,1+kgz/2-kjjj/2,200cm,gg2,areaCMSEL,S,GGLSLAKSLLADELE,ALLLDELE,ALLKDELE,ALLCMSEL,S,GG1CMSEL,A,GG2CM,GG,AREACMSEL,S,GGVEXT,ALL,0,0,4*KJJJ !切出钢轨扣件位置vsel,s,loc,y,hgz+0.0001,1 WPCSYS,-1WPOFFS,kgz/2vsbw,all*do,i,1,3WPOFFS,kjjjvsbw,all*enddoallsel,all!meshvsel,s,lo