(11.1.1)--11.1集中荷载作用下简支梁的有限元分析.pdf

《(11.1.1)--11.1集中荷载作用下简支梁的有限元分析.pdf》由会员分享，可在线阅读，更多相关《(11.1.1)--11.1集中荷载作用下简支梁的有限元分析.pdf（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、finish !结束当前操作/clear !清空当前数据库/filname,baem01 !取文件名!参数输入/prep7 pi=3.1415926 l=3100 !梁纵向长度 lp=1000 !荷载间距 w=150 !梁宽度 h=300 !梁深度 p=25 !保护层厚度 d_ps=8 !压筋直径 a_ps=0.25*pi*d_ps*2 d_ts=18 !拉筋直径 a_ts=0.25*pi*d_ts*2 d_gs=8 !箍筋直径 a_gs=0.25*pi*d_gs*2 fc=14.3e6 fy=350e6 !钢筋屈服强度 q=1e6 !均布载荷 e0=2e9 !钢筋强化系数 ex_c=1358

2、5e6 !混凝土杨氏模量 prxy_c=0.2 !混凝土泊松比 ex_ts=200e9 !钢筋杨氏模量 prxy_ts=0.27 !拉筋泊松比 ex_gs=200e9 !箍筋杨氏模量 prxy_sp=0.27 !箍筋泊松比 ex_sp=210e15 !垫板杨氏模量 prxy_sp=0.48 !垫板泊松比!建立几何模型 blc4,l,h,w block,100,200,0,-h/10,0,w !左支座 block,l-200,l-100,0,-h/10,0,w !右支座 vglue,all !粘结 !选择单元类型 et,1,solid65,1 et,2,link180!定义混凝土材料特性 mp,

3、ex,1,ex_c mp,prxy,1,prxy_c tb,miso,1,11 tbpt,0.0002,fc*0.19 tbpt,0.0004,fc*0.36 tbpt,0.0006,fc*0.51 tbpt,0.0008,fc*0.64 tbpt,0.001,fc*0.75 tbpt,0.0012,fc*0.84 tbpt,0.0014,fc*0.91 tbpt,0.0016,fc*0.96 tbpt,0.0018,fc*0.99 tbpt,0.002,fc tbpt,0.0033,fc tb,concr,1 tbdata,0.5,0.95,1.43e6,-1!定义钢筋材料特性 mp,ex,

4、2,ex_ts mp,prxy,2,prxy_ts tb,kinh,2,1,3 tbpt,fy/ex_ts,fy tbpt,fy/ex_ts+0.1,fy+e0*0.1 tbpt,fy/ex_ts+0.1001,1000 mp,ex,3,ex_gs mp,prxy,3,prxy_gs tb,bkin,3,1,2,1 tbdata,fy,e0!定义支座材料特性 mp,ex,4,ex_sp mp,prxy,4,prxy_sp!定义实常数 r,1,r,2,a_ts r,3,a_gs r,4,a_ps !切出箍筋位置 wprota,90*do,i,1,11 wpoffs,100 vsbw,all*en

5、ddo wpoffs,50*do,i,1,3 wpoffs,200 vsbw,all*enddo wpoffs,250*do,i,1,11 vsbw,all wpoffs,100*enddo wpcsys,-1!切出支座位置 wprota,90 wpoffs,150 vsbw,all wpcsys,-1 wprota,90 wpoffs,l-150 vsbw,all wpcsys,-1!切出纵向钢筋（轴线）位置 wpoffs,(p+d_ts/2)!移动工作平面 vsbw,all wpcsys,-1 !工作平面与笛卡尔坐标系对齐 wpoffs,w-(p+d_ts/2)vsbw,all wpcsy

6、s,-1 wprota,90 !旋转工作平面 wpoffs,-(p+d_ts/2)vsbw,all wpoffs,-h+2*(p+d_ts/2)vsbw,all wpcsys,-1 vglue,all!指定划分精度 elemsiz=50!划分纵向钢筋并赋予属性 lsel,s,loc,y,p+d_ts/2 !一号拉筋 lsel,r,loc,z,p+d_ts/2 cm,ts1,line !建组件 latt,2,2,2 lesize,all,elemsiz lsel,s,loc,y,p+d_ts/2 !二号拉筋 lsel,r,loc,z,w-(p+d_ts/2)cm,ts2,line latt,2,

7、2,2 lesize,all,elemsiz lsel,s,loc,y,h-(p+d_ts/2)!三号压筋 lsel,r,loc,z,p+d_ts/2 cm,ps1,line latt,2,4,2 lesize,all,elemsiz lsel,s,loc,y,h-(p+d_ts/2)!四号压筋 lsel,r,loc,z,w-(p+d_ts/2)cm,ps2,line latt,2,4,2 lesize,all,elemsiz!划分箍筋并赋予属性 lsel,s,loc,z,(p+d_ts/2)lsel,r,loc,y,h/2 lsel,u,loc,x,150 lsel,u,loc,x,l-15

8、0 lsel,u,loc,x,1050 lsel,u,loc,x,l-1050 cm,gs1,line latt,3,3,2 lesize,all,elemsiz !一、二间箍筋 lsel,s,loc,z,w-(p+d_ts/2)lsel,r,loc,y,h/2 lsel,u,loc,x,150 lsel,u,loc,x,l-150 lsel,u,loc,x,1050 lsel,u,loc,x,l-1050 cm,gs2,line latt,3,3,2 lesize,all,elemsiz lsel,s,loc,y,p+d_ts/2 lsel,r,loc,z,w/2 lsel,u,loc,x,

9、150 lsel,u,loc,x,l-150 lsel,u,loc,x,1050 lsel,u,loc,x,l-1050 cm,gs3,line latt,3,3,2 lesize,all,elemsiz lsel,s,loc,y,h-(p+d_ts/2)lsel,r,loc,z,w/2 lsel,u,loc,x,150 lsel,u,loc,x,l-150 lsel,u,loc,x,1050 lsel,u,loc,x,l-1050 cm,gs4,line latt,3,3,2 lesize,all,elemsiz cmsel,s,ts1 cmsel,a,ts2 cmsel,a,ps1 cms

10、el,a,ps2 cmsel,a,gs1 cmsel,a,gs2 cmsel,a,gs3 cmsel,a,gs4 lmesh,all lsel,all!划分混凝土及支座并赋予属性 vsel,s,loc,y,-h/20 vatt,4,1 mshkey,1 !映射划分 esize,elemsiz vsel,all vsel,u,loc,y,-h/20 vatt,1,1 mshkey,1 esize,elemsiz vsel,all vmesh,all !施加约束与荷载 lsel,s,loc,x,l-150 lsel,r,loc,y,-h/10 dl,all,uy lsel,s,loc,x,150

11、lsel,r,loc,y,-h/10 dl,all,uy dl,all,ux dl,all,uz lsel,all dtran asel,all lsel,all asel,all vsel,all asel,s,loc,y,h asel,r,loc,x,1000,1100 sfa,all,1,pres,q asel,s,loc,y,h asel,r,loc,x,2100,2200 sfa,all,1,pres,q asel,all sftran !求解!/solu antype,0 !静力分析 nlgeom,1 !打开几何大变形 nropt,modi !用修正法进行求解 eqslv,spar

12、,0 !方程求解器类型 nsubst,200,200,50,1 kbc,0 !使用线性荷载 cnvtol,u,0.015 !指定收敛值 neqit,50 !最大平衡迭代次数 outres,all,all autos,on solve finish !后处理!/post1!set,last!pldisp,1 !查看变形!esel,s,type,2 !选择钢筋!etable,saxl,ls,1!plls,saxl,saxl !定义单元表，绘制钢筋应力图!esel,s,type,1 !选择混凝土!/DEVICE,VECTOR,on!plcrack !设置矢量模型，绘制裂缝 !/post1 !pldisp,0!绘制变形图(缺省为最后一步的)!etable,sig,ls,1!定义应力表 !plls,sigi,sigi,1!绘制上述应力 !plnsol,s,y,0,1!绘制 Z 方向的应力