noying.zip_noying

%Rounded Square Column Section Design %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % material & %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% f_ck=35; f_y=400; E_s=200000; if (f_ck>28) beta1=0.85-(f_ck-28)*0.007; %베타 값 설정 else beta1=0.85; end epsilon_y=f_y/E_s; %철근의 항복 변형률 pi_factor=0.7; %[나선철근] %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % cross-section % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% A_sd=30; %주철근 지름[mm] A_ssd=10; %스트럽 지름[mm] D_As=706.85; %A_sd*A_sd*pi/4; %철근 1개 면적 (mm2-D30] R=200; %모서리 원 반지름 horizontal = 800; %가로 vertical = 800; %세로 A_g=pi*R*R+(horizontal-2*R)*2+(vertical-2*R)*2+(horizontal-2*R)*(vertical-2*R); %모서리가 둥근 사각형 기둥의 단면적(mm2) A_st=D_As*8; %전체 철근의 면적(mm2) pitch=480; %전단보강 철근 간격 min(주근직격*16,띠철근직격*48,최소 단면길이) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % position % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Cc=130; %피복두께 A_sd_dis=245; %주철근 간격 d=compute_dist(A_sd_dis,Cc,A_sd,A_ssd); %d(1)=155; %[mm] %d(2)=400; %[mm] %d(3)=645; %[mm] %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % P=0,Mn calculation % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% start_a=compute_eq(beta1,R,d,E_s,f_ck,f_y,D_As,horizontal); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % P-M curve % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% step=1; for a=start_a:1:d(3)+25 theata=acos(1-2*a/R); t=horizontal-2*R; if (a<R) A=R^2*(theata-sin(theata)*cos(theata))+t*a; Ay=abs((A*2*R^3/3*sin(theata)^3+t*a*(a/2+(R-a)))/(R^2*(theata-sin(theata)*cos(theata))+t*a)+t/2); elseif if (y_>0) yy=y_+t/2; else (y_<0) yy=y_-t/2; end Ay=abs(A*yy); else(t+R<a<t+2*R) A=(pi+sin(theata)*cos(theata)-theata)*R^2+2*t*R+t*a; Ay=abs(((pi/2*R^2*4*R/(3*pi)+2*R*t*t/2+t*a*(a-R)/2+(pi*R^2/2-R^2*(theata-sin(theata)*cos(theata))))*(pi*R^2/2*4*R/(3*pi)-2*R^3/3*sin(theata)^3)/(pi*R^2/2-R^2*(theata-sin(theata)*cos(theata)))+t)/((pi+sin(theata)*cos(theata)-theata)*R^2+2*t*R+t*a)+t/2); end ca1; for i=1:1:3 if(d(i)>c) epsilon_s(i)=-((d(i)-c)*0.003)/c; else epsilon_s(i)=((c-d(i))*0.003)/c; end end for j=1:1:3 f_sigma(j)=epsilon_s(j)*E_s; if(f_sigma(j)>f_y) f_sigma(j)=f_y; end if(f_sigma(j)<0&abs(f_sigma(j))>f_y) %&둘다 동시에 참 f_sigma(j)=-f_y; end if(d(j)<c) f_sigma(j)=f_sigma(j)-0.85*f_ck; end end sum_f_forceP=0; sum_f_forceM=0; for k=1:1:3 if(k==1|k==3) f_force(k)=f_sigma(k)*3*D_As/1000; else f_force(k)=f_sigma(k)*2*D_As/1000; end sum_f_forceP=sum_f_forceP+f_force(k); sum_f_forceM=sum_f_forceM+f_force(k)*(vertical/2-d(k))/1000; end P_n(step)=0.85*f_ck*A*0.001+sum_f_forceP; M_n(step)=0.85*f_ck*Ay*0.000001+sum_f_forceM; if(abs(epsilon_s(3))<epsilon_y) PP_n(step)=0.7*P_n(step); PM_n(step)=0.7*M_n(step); elseif((abs(epsilon_s(3))>epsilon_y&abs(epsilon_s(3))<0.005)|(abs(epsilon_s(3))==epsilon_y)) PP_n(step)=(0.7+(abs(epsilon_s(3))-epsilon_y)*50)*P_n(step); PM_n(step)=(0.7+(abs(epsilon_s(3))-epsilon_y)*50)*M_n(step); else PP_n(step)=0.85*P_n(step); PM_n(step)=0.85*M_n(step); end if(ceil(start_a)==ceil(a)) MS=PM_n(step) end step=step+1; end P_n(step)=(0.85*f_ck*(A_g-A_st)+f_y*A_st)/1000; %[나선철근-0.85, 띠철근-0.8] M_n(step)=0; P_max(step)=0.7*0.85*P_n(step); step2=1; for i2=1:1:step-1 if (PP_n(i2)<P_max(step)) PP_n_plot(i2)=PP_n(i2); PM_n_plot(i2)=PM_n(i2); end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % plot & %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% set_line1=[0 3000]; set_line2=[P_max(step) P_max(step)]; plot(M_n,P_n,'b-','linewidth',2) hold on set(gca,'XTick',[0:100:1300]) set(gca,'YTick',[0:1000:16000]) set(gca,'fontsize',18) ylabel('Pn(kN)','Fontsize',30) xlabel('Mn(kN*m)','Fontsize',30) plot(PM_n_plot,PP_n_plot,'r-','linewidth',2) plot(set_line1,set_line2,'k-','linewidth',2) hold off