1、finish/clear/filname,baem02!参数 输入/prep7 pi=3.1415926 l=3.100!梁纵 向长 度 lp=1.000!荷载 间距 w=0.150!梁宽 度 h=0.300!梁深 度 p=0.025!保 护层 厚度 d_ps=0.008!压 筋直 径 a_ps=0.25*pi*d_ps*2 d_ts=0.018!拉 筋直 径 a_ts=0.25*pi*d_ts*2 d_gs=0.008!箍 筋直 径 a_gs=0.25*pi*d_gs*2!选择 单元类 型 et,1,link160 et,2,solid164 et,3,solid164 et,4,solid
2、164 et,5,solid164!定义 混凝土 材料 特性 mp,ex,1,10e9 mp,dens,1,960 mp,nuxy,1,0.3!钢筋 mp,ex,2,10e9 mp,dens,2,960 mp,nuxy,2,0.3!混 凝土 mp,ex,3,10e9 mp,dens,3,960 mp,nuxy,3,0.3!支座 mp,ex,4,10e9 mp,dens,4,960 mp,nuxy,4,0.3!空气 mp,ex,5,10e9 mp,dens,5,960 mp,nuxy,5,0.3!炸药!定义 实常数 r,1,r,2,a_ts r,3,a_gs r,4,a_ps!建立 几何模 型
3、blc4,l,h,w block,0.100,0.200,0,-h/10,0,w!左支 座 block,l-0.200,l-0.100,0,-h/10,0,w!右 支座 vglue,all!粘结!切出 箍筋位 置 wprota,90*do,i,1,11 wpoffs,0.100 vsbw,all*enddo wpoffs,0.050*do,i,1,3 wpoffs,0.200 vsbw,all*enddo wpoffs,0.250*do,i,1,11 vsbw,all wpoffs,0.100*enddo wpcsys,-1!切出 支座位 置 wprota,90 wpoffs,0.150 vs
4、bw,all wpcsys,-1 wprota,90 wpoffs,l-0.150 vsbw,all wpcsys,-1!切出 纵向钢 筋（轴线）位 置 wpoffs,(p+d_ts/2)!移动 工作 平面 vsbw,all wpcsys,-1!工作 平面 与 笛卡尔 坐标 系对 齐 wpoffs,w-(p+d_ts/2)vsbw,all wpcsys,-1 wprota,90!旋 转工 作平 面 wpoffs,-(p+d_ts/2)vsbw,all wpoffs,-h+2*(p+d_ts/2)vsbw,all wpcsys,-1 vglue,all!指定 划分精 度 elemsiz=0.05
5、0!划分 纵向钢 筋并 赋予 属性 lsel,s,loc,y,p+d_ts/2!一号 拉筋 lsel,r,loc,z,p+d_ts/2 cm,ts1,line!建 组件 latt,1,2,1 lesize,all,elemsiz lsel,s,loc,y,p+d_ts/2!二号 拉筋 lsel,r,loc,z,w-(p+d_ts/2)cm,ts2,line latt,1,2,1 lesize,all,elemsiz lsel,s,loc,y,h-(p+d_ts/2)!三 号压 筋 lsel,r,loc,z,p+d_ts/2 cm,ps1,line latt,1,4,1 lesize,all,e
6、lemsiz lsel,s,loc,y,h-(p+d_ts/2)!四 号压 筋 lsel,r,loc,z,w-(p+d_ts/2)cm,ps2,line latt,1,4,1 lesize,all,elemsiz!划分 箍筋并 赋予 属性 lsel,s,loc,z,(p+d_ts/2)lsel,r,loc,y,h/2 lsel,u,loc,x,0.150 lsel,u,loc,x,l-0.150 lsel,u,loc,x,1.050 lsel,u,loc,x,l-1.050 cm,gs1,line latt,1,3,1 lesize,all,elemsiz!一、二间 箍筋 lsel,s,loc
7、,z,w-(p+d_ts/2)lsel,r,loc,y,h/2 lsel,u,loc,x,0.150 lsel,u,loc,x,l-0.150 lsel,u,loc,x,1.050 lsel,u,loc,x,l-1.050 cm,gs2,line latt,1,3,1 lesize,all,elemsiz lsel,s,loc,y,p+d_ts/2 lsel,r,loc,z,w/2 lsel,u,loc,x,0.150 lsel,u,loc,x,l-0.150 lsel,u,loc,x,1.050 lsel,u,loc,x,l-1.050 cm,gs3,line latt,1,3,1 lesi
8、ze,all,elemsiz lsel,s,loc,y,h-(p+d_ts/2)lsel,r,loc,z,w/2 lsel,u,loc,x,0.150 lsel,u,loc,x,l-0.150 lsel,u,loc,x,1.050 lsel,u,loc,x,l-1.050 cm,gs4,line latt,1,3,1 lesize,all,elemsiz cmsel,s,ts1 cmsel,a,ts2 cmsel,a,ps1 cmsel,a,ps2 cmsel,a,gs1 cmsel,a,gs2 cmsel,a,gs3 cmsel,a,gs4 lmesh,all lsel,all!划分 混凝土
9、 及支 座并 赋予 属性 vsel,s,loc,y,-h/20 vatt,3,1,3 mshkey,1!映 射划 分 esize,elemsiz vsel,all vsel,u,loc,y,-h/20 vatt,2,1,2 mshkey,1 esize,elemsiz vsel,all vmesh,all!施加 约束 lsel,s,loc,x,l-0.150 lsel,r,loc,y,-h/10 dl,all,uy lsel,s,loc,x,0.150 lsel,r,loc,y,-h/10 dl,all,uy lsel,all dtran allsel!创建 ALE 空间 block,-0.5
10、*l,1.5*l,-h,3*h,0,w vsel,s,volu,36 wprota,-90 wpoffs,3*h-0.1 vsbw,all wprota,-90 wpoffs,-l/2-0.05 vsbw,all wpoffs,0.1 vsbw,all wpcsys,-1 vsel,s,volu,68 vatt,5,1,5 esize,elemsiz vmesh,all vsel,s,volu,41 vsel,a,volu,46 vsel,a,volu,64 vsel,a,volu,52 vsel,a,volu,56 vatt,4,1,4 esize,elemsiz vmesh,all!建立
11、PART edpart,create!定义 约束 nsel,s,loc,z,0 nsel,s,loc,z,w nsel,a,loc,x,0 nsel,a,loc,x,l nsel,a,loc,y,-h/10 nsel,a,loc,y,3*h cm,nonreflect,node ednb,add,nonreflect allsel!设置 能量 选项 edenergy,1,1,1,1!设置 人工 体积 粘性 选项 edbvis,1.5,0.06!设置 时间 步长 因子 edcts,!设置 求解 时间 time,0.05!设置 文件 输出 步数 edrst,100 edhtime,00.001!设置 结果 文件 输出 类型 edopt,add,lsdyna!输出 关键 字文 件 edwrite,lsdyna,explosion,k
