vibration2008 发表于 2009-6-2 12:23
有没有普通 混凝土梁桥的算例, 单跨不复杂的就行,学习一下,
最好是包括桥面板的,谢谢
这个光混凝土的- !---------支座处量高
- H=2
- !---------跨中梁高
- h_m=1.3
- !---------翼板厚
- d1=0.2
- !---------顶板厚
- d3=0.2
- !---------底板厚
- d2=0.3
- !---------顶板内部倒角
- t1=0.1
- k1=0.2
- !---------顶板外部倒角
- t2=0.2
- k2=0.6
- !---------内箱宽
- w=0.8
- !---------支座处内箱高
- u=1.5
- !---------腹板宽
- bw=0.4
- !---------桥面宽
- B=3.8
- !---------第一层线性温度变化厚度
- ht1=0.1
- !---------中跨跨径
- l1=50
- !---------边跨跨径
- l2=25
- !---------支座不变高距离
- l3=2
- !
- !----------进入前处理
- /prep7
- !----------定义单元
- et,1,solid45
- !----------定义材料泊松比
- mp,prxy,1,0.2
- !----------定义密度
- mp,dens,1,2402
- !----------定义杨氏模量
- mp,ex,1,3.25e10
- !----------定义线膨胀系数
- mp,alpx,1,0.99e-5
- !!----------建立模型
- !----------建立关键点
- k,1,0,0,-(l1+2*l2)/2
- k,2,w/2,0,-(l1+2*l2)/2
- k,3,w/2+bw,0,-(l1+2*l2)/2
- k,4,w/2+bw,H-d1-t2,-(l1+2*l2)/2
- k,5,w/2+bw+k2,H-d1,-(l1+2*l2)/2
- k,6,B/2,H-d1,-(l1+2*l2)/2
- k,7,B/2,H,-(l1+2*l2)/2
- k,8,0,H,-(l1+2*l2)/2
- k,9,-B/2,H,-(l1+2*l2)/2
- k,10,-B/2,H-d1,-(l1+2*l2)/2
- k,11,-(w/2+bw+k2),H-d1,-(l1+2*l2)/2
- k,12,-(w/2+bw),H-d1-t2,-(l1+2*l2)/2
- k,13,-(w/2+bw),0,-(l1+2*l2)/2
- k,14,-w/2,0,-(l1+2*l2)/2
- k,15,w/2,H-d3-t1,-(l1+2*l2)/2
- k,16,w/2-k1,H-d3,-(l1+2*l2)/2
- k,17,0,H-d3,-(l1+2*l2)/2
- k,18,-(w/2-k1),H-d3,-(l1+2*l2)/2
- k,19,-w/2,H-d3-t1,-(l1+2*l2)/2
- !----------建立线
- l,kp(0,0,-(l1+2*l2)/2),kp(w/2,0,-(l1+2*l2)/2)
- l,kp(w/2,0,-(l1+2*l2)/2),kp(w/2+bw,0,-(l1+2*l2)/2)
- l,kp(w/2+bw,0,-(l1+2*l2)/2),kp(w/2+bw,H-d1-t2,-(l1+2*l2)/2)
- l,kp(w/2+bw,H-d1-t2,-(l1+2*l2)/2),kp(w/2+bw+k2,H-d1,-(l1+2*l2)/2)
- l,kp(w/2+bw+k2,H-d1,-(l1+2*l2)/2),kp(B/2,H-d1,-(l1+2*l2)/2)
- l,kp(B/2,H-d1,-(l1+2*l2)/2),kp(B/2,H,-(l1+2*l2)/2)
- l,kp(B/2,H,-(l1+2*l2)/2),kp(0,H,-(l1+2*l2)/2)
- l,kp(0,H,-(l1+2*l2)/2),kp(-B/2,H,-(l1+2*l2)/2)
- l,kp(-B/2,H,-(l1+2*l2)/2),kp(-B/2,H-d1,-(l1+2*l2)/2)
- l,kp(-B/2,H-d1,-(l1+2*l2)/2),kp(-(w/2+bw+k2),H-d1,-(l1+2*l2)/2)
- l,kp(-(w/2+bw+k2),H-d1,-(l1+2*l2)/2),kp(-(w/2+bw),H-d1-t2,-(l1+2*l2)/2)
- l,kp(-(w/2+bw),H-d1-t2,-(l1+2*l2)/2),kp(-(w/2+bw),0,-(l1+2*l2)/20)
- l,kp(-(w/2+bw),0,-(l1+2*l2)/2),kp(-w/2,0,-(l1+2*l2)/2)
- l,kp(-w/2,0,-(l1+2*l2)/2),kp(0,0,-(l1+2*l2)/2)
- l,kp(w/2,0,-(l1+2*l2)/2),kp(w/2,H-d3-t1,-(l1+2*l2)/2)
- l,kp(w/2,H-d3-t1,-(l1+2*l2)/2),kp(w/2-k1,H-d3,-(l1+2*l2)/2)
- l,kp(w/2-k1,H-d3,-(l1+2*l2)/2),kp(0,H-d3,-(l1+2*l2)/2)
- l,kp(0,H-d3,-(l1+2*l2)/2),kp(-(w/2-k1),H-d3,-(l1+2*l2)/2)
- l,kp(-(w/2-k1),H-d3,-(l1+2*l2)/2),kp(-w/2,H-d3-t1,-(l1+2*l2)/2)
- l,kp(-w/2,H-d3-t1,-(l1+2*l2)/2),kp(-w/2,0,-(l1+2*l2)/2)
- l,kp(0,H,-(l1+2*l2)/2),kp(0,H-d3,-(l1+2*l2)/2)
- l,kp(0,H-d3,-(l1+2*l2)/2),kp(0,0,-(l1+2*l2)/2)
- !----------建立面
- al,1,2,3,4,5,6,7,21,22
- al,22,21,8,9,10,11,12,13,14
- al,1,15,16,17,22
- al,22,18,19,20,14
- asba,1,3,,delete,delete
- asba,2,4,,delete,delete
- kwpave,kp(0,H,-(l1+2*l2)/2)
- wprota,90,90,90
- wpave,,H-ht1
- asbw,all,delete
- wpave,,H-d1
- asbw,all,delete
- wpave,,H-d1-ht1
- asbw,all,delete
- wpave,,H-d1-2*ht1
- asbw,all,delete
- aglue,all
- wprota,90,90
- wpave,w/2-k1
- asbw,all,delete
- wpave,-(w/2-k1)
- asbw,all,delete
- wpave,w/2
- asbw,all,delete
- wpave,-w/2
- asbw,all,delete
- wpave,w/2+bw/2
- asbw,all,delete
- wpave,-(w/2+bw/2)
- asbw,all,delete
- wpave,w/2+bw
- asbw,all,delete
- wpave,-(w/2+bw)
- asbw,all,delete
- wpave,w/2+bw+k2/2
- asbw,all,delete
- wpave,-(w/2+bw+k2/2)
- asbw,all,delete
- wpave,w/2+bw+k2
- asbw,all,delete
- wpave,-(w/2+bw+k2)
- asbw,all,delete
- !----------建立拉伸母线
- num_k_1=kpinqr(0,14)
- k,num_k_1+1,0,0,-(l1+2*l2)/2
- k,num_k_1+2,0,0,(l1+2*l2)/2
- l,num_k_1+1,num_k_1+2
- !----------拉伸体
- ksel,s,,,num_k_1+1,num_k_1+2,1
- lslk,s,1
- *get,num_l_1,line,,num,max
- alls
- vdrag,all,,,,,,num_l_1
- wprota,,90,90
- wpave,,,l1/2-l3/2
- vsbw,all,delete
- wpave,,,l1/2
- vsbw,all,delete
- wpave,,,l1/2+l3/2
- vsbw,all,delete
- wpave,,,-l1/2+l3/2
- vsbw,all,delete
- wpave,,,-l1/2
- vsbw,all,delete
- wpave,,,-l1/2-l3/2
- vsbw,all,delete
- alls
- !--------中跨底板下沿
- !--------抛物线模拟精度
- l_div=0.2 !-----------需能被l1-2*l3/2及l2-l3/2整除
- *dim,z_1,array,1/l_div*(l1-2*l3/2)+1,1
- *dim,y_1,array,1/l_div*(l1-2*l3/2)+1,1
- *do,i,1,1/l_div*(l1-2*l3/2)+1,1
- z_1(i)=-(l1/2-l3/2)+l_div*(i-1)
- num_temp=(H-h_m)/(l1/2-l3/2)/(l1/2-l3/2)
- y_1(i)=-num_temp*z_1(i)*z_1(i)+(H-h_m)
- *enddo
- alls
- num_k_2=kpinqr(0,14)
- *do,i,1,1/l_div*(l1-2*l3/2)+1,1
- k,num_k_2+1+i,1,y_1(i),z_1(i)
- *enddo
- *do,i,1,1/l_div*(l1-2*l3/2),1
- l,num_k_2+1+i,num_k_2+2+i
- *enddo
- alls
- num_k_3=kpinqr(0,14)
- k,num_k_3+1,-1,,-(l1/2-l3/2)
- l,num_k_2+2,num_k_3+1
- ksel,s,,,num_k_2+2
- ksel,a,,,num_k_3+1
- lslk,s,1
- *get,num_l_2,line,,num,max
- ksel,s,,,num_k_2+1,num_k_3
- lslk,s,1
- adrag,all,,,,,,num_l_2
- asll
- cm,a_sub_1,area
- alls
- vsel,s,loc,z,-l1/2,l1/2
- vsba,all,a_sub_1,,delete,delete
- ksel,s,loc,y,0,(H-h_m)
- lslk,s,1
- asll,s,1
- vsla,s,1
- vdele,all,,,1
- alls
- !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建议将体切开
- !--------------save as
- !--------边跨底板下沿(-Z侧)
- alls
- *dim,z_2_1,array,1/l_div*(l2-l3/2)+1,1
- *dim,y_2_1,array,1/l_div*(l2-l3/2)+1,1
- *do,i,1,1/l_div*(l2-l3/2)+1,1
- z_2_1(i)=-(l1/2+l2)+l_div*(i-1)
- num_temp=(H-h_m)/(l2-l3/2)/(l2-l3/2)
- y_2_1(i)=-num_temp*(z_2_1(i)+l1/2+l3/2)*(z_2_1(i)+l1/2+l2+l2-l3/2)
- *enddo
- alls
- num_k_4=kpinqr(0,14)
- *do,i,1,1/l_div*(l2-l3/2)+1,1
- k,num_k_4+1+i,1,y_2_1(i),z_2_1(i)
- *enddo
- *do,i,1,1/l_div*(l2-l3/2),1
- l,num_k_4+1+i,num_k_4+2+i
- *enddo
- alls
- num_k_5=kpinqr(0,14)
- k,num_k_5+1,-1,,-(l1/2+l3/2)
- l,num_k_5,num_k_5+1
- ksel,s,,,num_k_5
- ksel,a,,,num_k_5+1
- lslk,s,1
- *get,num_l_3,line,,num,max
- ksel,s,,,num_k_4+1,num_k_5
- lslk,s,1
- adrag,all,,,,,,num_l_3
- asll
- cm,a_sub_2,area
- alls
- vsel,s,loc,z,-(l1/2+l2),-l1/2
- vsba,all,a_sub_2,,delete,delete
- ksel,s,loc,y,0,(H-h_m)
- lslk,s,1
- asll,s,1
- vsla,s,1
- vdele,all,,,1
- !--------边跨底板下沿(+Z侧)
- alls
- *dim,z_3_1,array,1/l_div*(l2-l3/2)+1,1
- *dim,y_3_1,array,1/l_div*(l2-l3/2)+1,1
- *do,i,1,1/l_div*(l2-l3/2)+1,1
- z_3_1(i)=(l1/2+l2)-l_div*(i-1)
- num_temp=(H-h_m)/(l2-l3/2)/(l2-l3/2)
- y_3_1(i)=-num_temp*(z_3_1(i)-(l1/2+l3/2))*(z_3_1(i)-(l1/2+l2+l2-l3/2))
- *enddo
- alls
- num_k_6=kpinqr(0,14)
- *do,i,1,1/l_div*(l2-l3/2)+1,1
- k,num_k_6+1+i,1,y_3_1(i),z_3_1(i)
- *enddo
- *do,i,1,1/l_div*(l2-l3/2),1
- l,num_k_6+1+i,num_k_6+2+i
- *enddo
- alls
- num_k_7=kpinqr(0,14)
- k,num_k_7+1,-1,,l1/2+l3/2
- l,num_k_7,num_k_7+1
- ksel,s,,,num_k_7
- ksel,a,,,num_k_7+1
- lslk,s,1
- *get,num_l_4,line,,num,max
- ksel,s,,,num_k_6+1,num_k_7
- lslk,s,1
- adrag,all,,,,,,num_l_4
- asll
- cm,a_sub_3,area
- alls
- vsel,s,loc,z,l1/2+l2,l1/2
- vsba,all,a_sub_3,,delete,delete
- ksel,s,loc,y,0,(H-h_m)
- lslk,s,1
- asll,s,1
- vsla,s,1
- vdele,all,,,1
- !--------跨中底板上沿
- alls
- *dim,z_2,array,1/l_div*(l1-2*l3/2)+1,1
- *dim,y_2,array,1/l_div*(l1-2*l3/2)+1,1
- *do,i,1,1/l_div*(l1-2*l3/2)+1,1
- z_2(i)=-(l1/2-l3/2)+l_div*(i-1)
- num_temp=(H-h_m)/(l1/2-l3/2)/(l1/2-l3/2)
- y_2(i)=-num_temp*z_2(i)*z_2(i)+(H-h_m)+d2
- *enddo
- alls
- num_k1_2=kpinqr(0,14)
- *do,i,1,1/l_div*(l1-2*l3/2)+1,1
- k,num_k1_2+1+i,1,y_2(i),z_2(i)
- *enddo
- *do,i,1,1/l_div*(l1-2*l3/2),1
- l,num_k1_2+1+i,num_k1_2+2+i
- *enddo
- alls
- num_k1_3=kpinqr(0,14)
- k,num_k1_3+1,-1,d2,-(l1/2-l3/2)
- l,num_k1_2+2,num_k1_3+1
- ksel,s,,,num_k1_2+2
- ksel,a,,,num_k1_3+1
- lslk,s,1
- *get,num_l1_2,line,,num,max
- ksel,s,,,num_k1_2+1,num_k1_3
- lslk,s,1
- adrag,all,,,,,,num_l1_2
- asll
- cm,a_sub_1,area
- alls
- vsel,s,loc,z,-l1/2,l1/2
- vsba,all,a_sub_1,,delete,delete
- !--------边跨底板上沿(-Z侧)
- alls
- *dim,z_2_2,array,1/l_div*(l2-l3/2)+1,1
- *dim,y_2_2,array,1/l_div*(l2-l3/2)+1,1
- *do,i,1,1/l_div*(l2-l3/2)+1,1
- z_2_2(i)=-(l1/2+l2)+l_div*(i-1)
- num_temp=(H-h_m)/(l2-l3/2)/(l2-l3/2)
- y_2_2(i)=-num_temp*(z_2_2(i)+l1/2+l3/2)*(z_2_2(i)+l1/2+l2+l2-l3/2)+d2
- *enddo
- alls
- num_k1_4=kpinqr(0,14)
- *do,i,1,1/l_div*(l2-l3/2)+1,1
- k,num_k1_4+1+i,1,y_2_2(i),z_2_2(i)
- *enddo
- *do,i,1,1/l_div*(l2-l3/2),1
- l,num_k1_4+1+i,num_k1_4+2+i
- *enddo
- alls
- num_k1_5=kpinqr(0,14)
- k,num_k1_5+1,-1,d2,-(l1/2+l3/2)
- l,num_k1_5,num_k1_5+1
- ksel,s,,,num_k1_5
- ksel,a,,,num_k1_5+1
- lslk,s,1
- *get,num_l1_3,line,,num,max
- ksel,s,,,num_k1_4+1,num_k1_5
- lslk,s,1
- adrag,all,,,,,,num_l1_3
- asll
- cm,a_sub_2,area
- alls
- vsel,s,loc,z,-(l1/2+l2),-l1/2
- vsba,all,a_sub_2,,delete,delete
- ksel,s,loc,y,0,(H-h_m)
- lslk,s,1
- asll,s,1
- vsla,s,1
- vdele,all,,,1
- !--------边跨底板上沿(+Z侧)
- alls
- *dim,z_3_2,array,1/l_div*(l2-l3/2)+1,1
- *dim,y_3_2,array,1/l_div*(l2-l3/2)+1,1
- *do,i,1,1/l_div*(l2-l3/2)+1,1
- z_3_2(i)=(l1/2+l2)-l_div*(i-1)
- num_temp=(H-h_m)/(l2-l3/2)/(l2-l3/2)
- y_3_2(i)=-num_temp*(z_3_2(i)-(l1/2+l3/2))*(z_3_2(i)-(l1/2+l2+l2-l3/2))+d2
- *enddo
- alls
- num_k1_6=kpinqr(0,14)
- *do,i,1,1/l_div*(l2-l3/2)+1,1
- k,num_k1_6+1+i,1,y_3_2(i),z_3_2(i)
- *enddo
- *do,i,1,1/l_div*(l2-l3/2),1
- l,num_k1_6+1+i,num_k1_6+2+i
- *enddo
- alls
- num_k1_7=kpinqr(0,14)
- k,num_k1_7+1,-1,d2,l1/2+l3/2
- l,num_k1_7,num_k1_7+1
- ksel,s,,,num_k1_7
- ksel,a,,,num_k1_7+1
- lslk,s,1
- *get,num_l1_4,line,,num,max
- ksel,s,,,num_k1_6+1,num_k1_7
- lslk,s,1
- adrag,all,,,,,,num_l1_4
- asll
- cm,a_sub_3,area
- alls
- vsel,s,loc,z,l1/2+l2,l1/2
- vsba,all,a_sub_3,,delete,delete
- ksel,s,loc,y,0,(H-h_m)
- lslk,s,1
- asll,s,1
- vsla,s,1
- vdele,all,,,1
- alls
- ldele,num_l_1,,,1
- ldele,num_l_2,,,1
- ldele,num_l_3,,,1
- ldele,num_l_4,,,1
- ldele,num_l1_2,,,1
- ldele,num_l1_3,,,1
- ldele,num_l1_4,,,1
- nummrg,all
- numcmp,all
- !!!!!!!!!!!!!!!!!!!!!!!!!!!save as
- alls
- ksel,s,loc,x,w/2+bw
- ksel,r,loc,y,0
- ksel,r,loc,z,0
- *get,n_l_1,kp,,num,max
- ksel,s,loc,x,-0.8
- ksel,r,loc,y,0
- ksel,r,loc,z,0
- *get,n_l_2,kp,,num,max
- alls
- l,n_l_1,n_l_2
- ksel,s,loc,y,0,1
- ksel,r,loc,x,0.8
- lslk,s,1
- asll,s,1
- vdrag,all,,,,,,2696
- vovlap,all
- !----------定义单元尺寸
- esize,0.11,0
- !----------单元网格划分
- alls
- ksel,s,loc,y,0,1
- lslk,s,1
- asll,s,1
- vsla,s,1
- vsweep,all
- esize,0.1,0
- vsel,s,loc,y,1,2
- vsweep,all
- alls
- vsweep,all
复制代码 |