*윤* 스토어

*윤*

개인인증

팔로워0 팔로우

소개

등록된 소개글이 없습니다.

전문분야 등록된 전문분야가 없습니다.

판매자 정보

학교정보

입력된 정보가 없습니다.

직장정보

입력된 정보가 없습니다.

자격증

입력된 정보가 없습니다.

판매지수

판매중 자료수

15개
전체 판매량

107개
최근 3개월 판매량

0개
자료후기 점수

평균C
자료문의 응답률

-

전체자료 15개

[구조동역학] (matlab) Direct integration Method(뉴마크베타)방법을 이용한 다자유도 지진해석 평가C아쉬워요

..FILE:ELCEN_NS.DAT0.02 -0.001430.04 -0.011020.06 -0.010310.08 -0.008980.10 -0.009690.12 -0.012240.14 -0.014490.16 -0.013060.18 -0.011220.20 -0.008670.22 -0.008670.24 -0.013370.26 -0.017960.28 -0.019800.30 -0.016530.32 -0.014690.34 -0.011020.36 -0.008370.38 -0.004290.40 -0.006730.42 -0.013370.44 -0.019390.46 -0.020000.48 -0.006730.50 0.003060.52 0.014390.54 -0.005000.56 -0.013060.58 -0.014690.60 -0.020710.62 -0.026530.64 -0.033160.66 -0.031220.68 -0.017550.70 -0.020100.72 -0.016630.74 -0.016730.76 -0.006840.78 0.002550.80 0.015310.82 0.024080.84 0.025710.86 0.034290.88 0.047240.90 0.050200.92 0.042760.94 0.036630.96 0.027650.98 0.023981.00 0.034591.02 0.042041.04 0.054081.06 0.065201.08 0.074691.10 0.066531.12 0.061121.14 0.040821.16 0.040821.18 0.006431.20 -0.052551.22 -0.080311.24 -0.061531.26 -0.049391.28 -0.025511.30 -0.006021.32 0.013671.34 0.031431.36 0.050921.38 0.072451.40 0.101531.42 0.124391.44 0.156021.46 0.147861.48 0.117861.50 0.095411.52 0.091021.54 0.094491.56 0.085611.53.56 -0.077763.58 -0.057043.60 -0.021943.62 -0.012863.64 -0.068783.66 -0.033063.68 -0.034393.70 -0.011123.72 0.001733.74 0.030513.76 0.049803.78 0.062043.80 0.022653.82 -0.003273.84 -0.025003.86 0.007863.88 0.021533.90 0.057963.92 0.084293.94 0.123063.96 0.150823.98 0.177244.00 0.042964.02 0.002964.04 0.026434.06 0.029904.08 -0.005614.10 -0.015004.12 0.014594.14 0.021024.16 0.050924.18 0.065824.20 0.097654.22 0.115104.24 0.147654.26 0.166224.28 0.198474.30 0.189394.32 0.202454.34 0.180514.36 0.127554.38 -0.123164.40 -0.055314.42 -0.039184.44 -0.031734.46 -0.114084.48 -0.169494.50 -0.251434.52 -0.206634.54 -0.187244.56 -0.134394.58 -0.097964.60 -0.033164.62 0.015714.64 0.083274.66 0.134594.68 0.185514.70 -0.005924.72 -0.017244.74 0.029084.76 0.045614.78 0.100314.80 0.145314.82 0.189084.84 0.250614.86 0.171944.88 -0.140824.90 -0.101944.92 -0.111124.94 -0.092554.96 -0.047864.98 -0.127555.00 -0.215415.02 -0.165005.04 -0.172655.06 -0.133275.08 -0.113375.10 -0.078885.12 -0.052045.14 -0.05551.14 0.032657.16 0.042767.18 0.012557.20 -0.016337.22 -0.020827.24 -0.008377.26 -0.021027.28 -0.013987.30 -0.005617.32 0.005417.34 0.013677.36 0.027147.38 0.023677.40 0.008067.42 -0.000827.44 0.020417.46 0.044397.48 0.050207.50 0.019497.52 0.009397.54 -0.002247.56 -0.002147.58 0.005317.60 0.009497.62 0.026027.64 0.037557.66 0.053577.68 0.055207.70 0.043377.72 0.040617.74 0.057047.76 0.077147.78 0.037247.80 0.041947.82 0.010007.84 -0.020827.86 -0.025417.88 -0.041337.90 -0.042147.92 -0.048067.94 -0.044187.96 -0.046737.98 -0.005828.00 0.018168.02 -0.021228.04 -0.050208.06 -0.054088.08 -0.036948.10 -0.041338.12 -0.031438.14 -0.032248.16 -0.027048.18 -0.027048.20 -0.027458.22 -0.035208.24 -0.031538.26 -0.022148.28 -0.007968.30 0.008888.32 0.028678.34 0.031638.36 0.036538.38 0.034808.40 0.036538.42 0.029298.44 0.031128.46 0.011438.48 0.021848.50 0.013888.52 0.039188.54 -0.087868.56 -0.137658.58 -0.136948.60 -0.138168.62 -0.121738.64 -0.106338.66 -0.084598.68 -0.066438.70 -0.045318.72 -0.026332)];R0 = [R0 ; zeros(NPT-length(R0),1)];R0 = -m*[1 ; 1 ; 1]*R0' ;alpha = 0.25;delta = 0.5;dt = 0.02;t=[0 elcen_ns(:,1)'];a0 = 1/(alpha*dt^2);a1 = delta/(alpha*dt);a2 = 1/(alpha*dt);a3 = 1/(2*alpha) -1 ;a4 = delta/alpha -1 ;a5 = (dt/2)*(delta/alpha - 2);a6 = dt*(1-delta);a7 = delta*dt;U0 = [0 ; 0 ; 0];U1 = [0 ; 0 ; 0];U2 = [0 ; 0 ; 0];Kbar = k+ a0*m +a1*c;invK=inv(Kbar);for itmp=1:length(R0)R = R0(:,itmp) + m*(a0*U0 + a2*U1 + a3*U2) + c*(a1*U0 + a4*U1 +a5*U2);U0new(:,itmp) = (invK*R)';U2new(:,itmp) = a0*(U0new(:,itmp)-U0)-(a2*U1)-(a3*U2);U1new(:,itmp) = U1 + a6*U2 + a7*U2new(:,itmp) ;U0=U0new(:,itmp) ;U1=U1new(:,itmp) ;U2=U2new(:,itmp) ;end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 그래프%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%subplot(311) ; plot(t,U0new)title('bf{Direct intergration Method}','fontsize',10)xlabel('Time [sec]')ylabel('Displacement')set(gcf, 'color', 'w')legend('1층','2층','3층')grid onsubplot(312) ; plot(t,U1new)xlabel('Time %%%%%%%%%%%%%%%%%%%%%% 그래프%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%subplot(311) ; plot(t,U0new)title('bf{Direct intergration Method}','fontsize',10)xlabel('Time [sec]')ylabel('Displacement')set(gcf, 'color', 'w')legend('1층','2층','3층')grid onsubplot(312) ; plot(t,U1new)xlabel('Time [sec]')ylabel('Velocity')grid onsubplot(313) ; plot(t,U2new)xlabel('Time [sec]')ylabel('Acceleration')grid on..FILE:newmark.jpg..FILE:story4.asv% close allclear allclcfigure%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 뉴마크방법을 이용한 MDOF 해석%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%m=[1 0 0 0 ; 0 1 0 0 ; 0 0 1 0 ; 0 0 0 1];k=[2 -1 0 0 ; -1 2 -1 0 ; 0 -1 2 -1 ; 0 0 -1 1];xi=0.1;[vec,value]=eig(k,m);w=sqrt(value);% c=2*xi*w;c=inv((vec)') * (vec'*m*vec) * 2 * xi * w * inv(vec);NPT = 501 ;load elcen_ns.dat ;R0 = [0 ; elcen_ns(:,2)];R0 = [R0 ; zeros(NPT-length(R0),1)];R0 = -m*[1 ; 1 ; 1 ;1]*R0' ;alpha = 0.25;delta = 0.5;dt = 0.02;t=[0 elcen_ns(:,1)'];a0db

프로그램소스| 2008.06.23| 9,000원| 조회(1,737)

MatLab, 구조동역학, Direct integration M, 뉴마크, 지진..

미리보기

닫기
[구조동역학] (matlab) 주파수영역방법 (frequency domain method)을 이용한 다자유도 지진해석

..FILE:c1.jpg..FILE:c2.jpg..FILE:ELCEN_NS.DAT0.02 -0.001430.04 -0.011020.06 -0.010310.08 -0.008980.10 -0.009690.12 -0.012240.14 -0.014490.16 -0.013060.18 -0.011220.20 -0.008670.22 -0.008670.24 -0.013370.26 -0.017960.28 -0.019800.30 -0.016530.32 -0.014690.34 -0.011020.36 -0.008370.38 -0.004290.40 -0.006730.42 -0.013370.44 -0.019390.46 -0.020000.48 -0.006730.50 0.003060.52 0.014390.54 -0.005000.56 -0.013060.58 -0.014690.60 -0.020710.62 -0.026530.64 -0.033160.66 -0.031220.68 -0.017550.70 -0.020100.72 -0.016630.74 -0.016730.76 -0.006840.78 0.002550.80 0.015310.82 0.024080.84 0.025710.86 0.034290.88 0.047240.90 0.050200.92 0.042760.94 0.036630.96 0.027650.98 0.023981.00 0.034591.02 0.042041.04 0.054081.06 0.065201.08 0.074691.10 0.066531.12 0.061121.14 0.040821.16 0.040821.18 0.006431.20 -0.052551.22 -0.080311.24 -0.061531.26 -0.049391.28 -0.025511.30 -0.006021.32 0.013671.34 0.031431.36 0.050921.38 0.072451.40 0.101531.42 0.124391.44 0.156021.46 0.147861.48 0.117861.50 0.095411.52 0.0910216535.30 -0.082045.32 -0.166735.34 -0.087655.36 -0.098065.38 -0.040415.40 -0.015005.42 0.032555.44 0.066125.46 0.089395.48 0.048165.50 0.020205.52 -0.002765.54 0.029805.56 0.045415.58 0.080105.60 0.105415.62 0.137965.64 0.163885.66 0.189905.68 0.130715.70 0.065315.72 0.020825.74 0.032045.76 0.038065.78 0.050615.80 0.023985.82 -0.008575.84 -0.017145.86 -0.011535.88 -0.023375.90 -0.025315.92 -0.016025.94 -0.007045.96 0.015005.98 0.038676.00 0.059086.02 0.026026.04 -0.004186.06 -0.043676.08 -0.013576.10 0.009696.12 0.023476.14 -0.013166.16 -0.005106.18 0.008166.20 0.021436.22 0.038786.24 0.052046.26 0.016026.28 -0.003276.30 -0.011336.32 0.000516.34 0.007766.36 0.003576.38 -0.009696.40 -0.003676.42 -0.001636.44 0.003886.46 0.008676.48 -0.005716.50 -0.031026.52 -0.042966.54 -0.024906.56 -0.024086.58 -0.018066.60 -0.013166.62 -0.001846.64 0.020716.66 -0.011026.68 -0.009296.70 -0.003476.72 -0.010826.74 -0.011336.76 -0.010106.78 -0.000206.80 0.007456.82 0.023986.84 0.036226.86 0.071946.88 0.079 -1 3 -2 ; 0 -2 5]*100;K=[5 -2 0 ; -2 3 -1 ; 0 -1 1]*100;xi=0.1;[vec,value]=eig(K,m);w=sqrt(value);% c=2*xi*w;c=inv((vec)') * (vec'*m*vec) * 2 * xi * w * inv(vec);% % dt = 0.02;% % T = 10;% % t = 0:dt:T ;% % df = 1/T ;% % fmax = 1/dt;% % f = 0.1:df:fmax;% % w = 2*pi*f;% % c = 2*xi*w%%% F(w)를 구하는 루틴dt = 0.02;T = 10;t = 0:dt:T ;df = 1/T ;fmax = 1/dt;f = 0:df:fmax;load elcen_ns.dat;% t=elcen_ns(:,1)';ga=elcen_ns(:,2);ga = ga' ;ga = [0 ga] ;ft = -m*[1 ; 1 ; 1] * ga ;Fw = (fft(ft'))' * dt;%%%%%% H(w)를 구하는 루틴 %%%%%%%w= 2*pi*f;dw = w(2)-w(1) ;Hw=[];for k=1:1:length(w) ;wk=w(k);% Hw(:,:,k)= inv(K-wk.^2*m + i*wk*c);% Xw(:,k) = Hw(:,:,k) * Fw(:,k);Hw = inv(K-wk.^2*m + i*wk*c);Xw(:,k) = Hw' * Fw(:,k);Xwdot(:,k) = [(-(i*wk.)*Xw(:,k))];Xw2dot(:,k) = [(-wk.*wk.* Xw(:,k))];end%%% x축으로 time의 함수로 표현% % % N = 짝수% % % for k=(length(t))/2+2:length(t);% % % Xw(:,k) = Xw(:,length(t)+2-k)' ;% % % Xwdot(:,k) = Xwdot(:,length(t)+2-k)' ;% % % Xw2dot(:,k) = Xw2dot(:,length(t)+2-k)' ;% % % end% % % N = 홀수for k=ro; 1 ; 1] * ga ;Fw = (fft(ft'))' * dt;%%%%%% H(w)를 구하는 루틴 %%%%%%%w= 2*pi*f;dw = w(2)-w(1) ;Hw=[];for k=1:1:length(w) ;wk=w(k);% Hw(:,:,k)= inv(K-wk.^2*m + i*wk*c);% Xw(:,k) = Hw(:,:,k) * Fw(:,k);Hw = inv(K-wk.^2*m + i*wk*c);Xw(:,k) = Hw' * Fw(:,k);Xwdot(:,k) = [(-i*wk.*Xw(:,k))];Xw2dot(:,k) = [(-wk.*wk.* Xw(:,k))];end%%% x축으로 time의 함수로 표현% % % N = 짝수% % % for k=(length(t))/2+2:length(t);% % % Xw(:,k) = Xw(:,length(t)+2-k)' ;% % % Xwdot(:,k) = Xwdot(:,length(t)+2-k)' ;% % % Xw2dot(:,k) = Xw2dot(:,length(t)+2-k)' ;% % % end% % % N = 홀수for k=round(length(t)/2) +1 : length(t);Xw(:,k) = Xw(:,length(t)+2-k)' ;Xwdot(:,k) = Xwdot(:,length(t)+2-k)' ;Xw2dot(:,k) = Xw2dot(:,length(t)+2-k)' ;endxt = (ifft(Xw'))' / dt ;vt = (ifft(Xwdot'))' / dt ;at = (ifft(Xw2dot'))' / dt ;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%figure ;subplot(311); plot(f,real(Xw),'b')title('bf{X(w)}','fontsize',10)set(gcf, 'color', 'w')xlabel('Frequncy [Hz]')ylabel('real{(w)}')grid onsubplot(312); plot(f,real(XwdoXw(:,k) = Xw(:,length(freq)+2-k)' ;Xwdot(:,k) = Xwdot(:,length(freq)+2-k)' ;Xw2dot(:,k) = Xw2dot(:,length(freq)+2-k)' ;end% % % N = 홀수% % % for k=round(length(t)/2) +1 : length(t);% % % Xw(:,k) = Xw(:,length(t)+2-k)' ;% % % Xwdot(:,k) = Xwdot(:,length(t)+2-k)' ;% % % Xw2dot(:,k) = Xw2dot(:,length(t)+2-k)' ;% % % endxt = (ifft(Xw',Nfft))' / dt ;vt = (ifft(Xwdot',Nfft))' / dt ;at = (ifft(Xw2dot',Nfft))' / dt ;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Frquency Domain Method 그래프%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%subplot(311) ; plot(t,real(xt)) ;set(gca,'Xlim',[0 10])title('bf{Frquency Domain Method를 이용한 MDOF해석}','fontsize',10)xlabel('time')ylabel('Displacement')set(gcf, 'color', 'w')legend('1층','2층','3층')grid onsubplot(312) ; plot(t,real(vt))set(gca,'Xlim',[0 10])xlabel('time')ylabel('Velocity')grid onsubplot(313) ; plot(t,real(at))set(gca,'Xlim',[0 10])xlabel('time')ylabel('Acceleration')grid on..FILE:exx1_jj.m%%%%%%%%%%%%%%%%%%%%%.db

프로그램소스| 2008.06.23| 9,000원| 조회(1,193)

구조동역학, 주파수영역방법, matrab, frequency domain me..

미리보기

닫기
[구조동역학] (matlab) 모드중첩법을 이용한 다자유도의 해석

..FILE:ELCEN_NS.DAT0.02 -0.001430.04 -0.011020.06 -0.010310.08 -0.008980.10 -0.009690.12 -0.012240.14 -0.014490.16 -0.013060.18 -0.011220.20 -0.008670.22 -0.008670.24 -0.013370.26 -0.017960.28 -0.019800.30 -0.016530.32 -0.014690.34 -0.011020.36 -0.008370.38 -0.004290.40 -0.006730.42 -0.013370.44 -0.019390.46 -0.020000.48 -0.006730.50 0.003060.52 0.014390.54 -0.005000.56 -0.013060.58 -0.014690.60 -0.020710.62 -0.026530.64 -0.033160.66 -0.031220.68 -0.017550.70 -0.020100.72 -0.016630.74 -0.016730.76 -0.006840.78 0.002550.80 0.015310.82 0.024080.84 0.025710.86 0.034290.88 0.047240.90 0.050200.92 0.042760.94 0.036630.96 0.027650.98 0.023981.00 0.034591.02 0.042041.04 0.054081.06 0.065201.08 0.074691.10 0.066531.12 0.061121.14 0.040821.16 0.040821.18 0.006431.20 -0.052551.22 -0.080311.24 -0.061531.26 -0.049391.28 -0.025511.30 -0.006021.32 0.013671.34 0.031431.36 0.050921.38 0.072451.40 0.101531.42 0.124391.44 0.156021.46 0.147861.48 0.117861.50 0.095411.52 0.091021.54 0.094491.56 0.085611.5 -0.111124.94 -0.092554.96 -0.047864.98 -0.127555.00 -0.215415.02 -0.165005.04 -0.172655.06 -0.133275.08 -0.113375.10 -0.078885.12 -0.052045.14 -0.055515.16 -0.122455.18 -0.123375.20 -0.118165.22 -0.116845.24 -0.073165.26 -0.055715.28 0.006535.30 -0.082045.32 -0.166735.34 -0.087655.36 -0.098065.38 -0.040415.40 -0.015005.42 0.032555.44 0.066125.46 0.089395.48 0.048165.50 0.020205.52 -0.002765.54 0.029805.56 0.045415.58 0.080105.60 0.105415.62 0.137965.64 0.163885.66 0.189905.68 0.130715.70 0.065315.72 0.020825.74 0.032045.76 0.038065.78 0.050615.80 0.023985.82 -0.008575.84 -0.017145.86 -0.011535.88 -0.023375.90 -0.025315.92 -0.016025.94 -0.007045.96 0.015005.98 0.038676.00 0.059086.02 0.026026.04 -0.004186.06 -0.043676.08 -0.013576.10 0.009696.12 0.023476.14 -0.013166.16 -0.005106.18 0.008166.20 0.021436.22 0.038786.24 0.052046.26 0.016026.28 -0.003276.30 -0.011336.32 0.000516.34 0.007766.36 0.003576.38 -0.009696.40 -0.003676.42 -0.001636.44 0.003886.46 0.008676.48 -0.005716.50 -0.0310257869.86 0.076849.88 0.081739.90 0.060419.92 0.031029.94 0.002359.96 0.006539.98 -0.0414310.00 -0.04602..FILE:mode.asvclcclear allclose allm=[1 0 0 ; 0 1 0 ; 0 0 1];k=[; -1 3 -2 ; 0 -2 5];k=[1 -1 0 ; -1 3 -2 ; 0 -2 5];% eig(k.*m);[vec,value]=eig(k,m)% eigenvec=[ vec(:,1)./0.8433 , vec(:,2)./0.5277 , vec(:,3)./0.1019 ];H=[0 1 2 3];A=[0 0.2149 0.4927 0.8433];B=[0 -0.5049 -0.6831 0.5277];C=[0 0.8360 -0.5392 0.1019];w=sqrt(value)f=w./(2*pi)subplot(131); plot(A,H,'-Or')xlabel('it{eigen vector}','fontsize',8)ylabel('it{height}','fontsize',8)title('bf{1 Mode}','fontsize',10)grid onsubplot(132); plot(B,H,'-Og')xlabel('it{eigen vector}','fontsize',8)ylabel('it{height}','fontsize',8)title('bf{2 Mode}','fontsize',10)grid onsubplot(133);plot(C,H,'-O')xlabel('it{eigen vector}','fontsize',8)ylabel('it{height}','fontsize',8)title('bf{3 Mode}','fontsize',10)grid onset(gcf, 'color', 'w')..FILE:MODE.jpg..FILE:mode.mclcclear allclose allm=[1 0 0 ; 0 1 0 ; 0 0 1];k=[5 -2 0 ; -2 3 -1 ; 0 -1 1]*100%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 단자유도로 해석하기 위한 m,c,k,p 변환 --- 2 ----modalM=vec'*m*vecmodalC=2*modalM(5)*0.1*wmodalK=vec'*k*vecM2=modalM(5);C2=modalC(5);K2=modalK(5);q0=0;qdot0=0;load elcen_ns.dat;t=elcen_ns(:,1);ga=elcen_ns(:,2);pt=-M*elcen_ns(:,2);P=[pt pt pt] ;Pt=(vec'*P')' ;R0 = Pt(:,2) ;%%%%%% 뉴마크를 이용한 q, qdot q2dot 구하기alpha = 0.25;delta = 0.5;dt = 0.02;t=elcen_ns(:,1);a0 = 1/(alpha*dt^2);a1 = delta/(alpha*dt);a2 = 1/(alpha*dt);a3 = 1/(2*alpha) -1 ;a4 = delta/alpha -1 ;a5 = (dt/2)*(delta/alpha - 2);a6 = dt*(1-delta);a7 = delta*dt;U0 = 0;U1 = 0;U2 = 0;Kbar = K2+ a0*M2 +a1*C2;invK=inv(Kbar);for itmp=1:length(R0)R = R0(itmp) + M2*(a0*U0 + a2*U1 + a3*U2) + C2*(a1*U0 + a4*U1 +a5*U2);U0new(itmp) = invK*R ;U2new(itmp) = a0*(U0new2(:,itmp)-U0)-(a2*U1)-(a3*U2);U1new(itmp) = U1 + a6*U2 + a7*U2new2(:,itmp) ;U0=U0new2(itmp) ;U1=U1new2(itmp) ;U2=U2new2(itmp) ;end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Pt(:,1) ;%%%%%% 뉴마크를 이용한 q, qdot q2dot 구하기alpha = 0.25;delta = 0.5;dt = 0.02;t=elcen_ns(:,1);a0 = 1/(alpha*dt^2);a1 = delta/(alpha*dt);a2 = 1/(alpha*dt);a3 = 1/(2*alpha) -1 ;a4 = delta/alpha -1 ;a5 = (dt/2)*(delta/alpha - 2);a6 = dt*(1-delta);a7 = delta*dt;U0 = 0;U1 = 0;U2 = 0;Kbar = K1+ a0*M1 +a1*C1;invK=inv(Kbar);for itmp=1:length(R0)R = R0(itmp) + M1*(a0*U0 + a2*U1 + a3*U2) + C1*(a1*U0 + a4*U1 +a5*U2);U0new1(itmp) = invK*R ;U2new1(itmp) = a0*(U0new1(:,itmp)-U0)-(a2*U1)-(a3*U2);U1new1(itmp) = U1 + a6*U2 + a7*U2new1(:,itmp) ;U0=U0new1(itmp) ;U1=U1new1(itmp) ;U2=U2new1(itmp) ;end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 단자유도로 해석하기 위한 m,c,k,p 변환 --- 2 ----modalM=vec'*m*vec;modalC=2*modalM(5)*0.1*w;modalK=vec'*k*vec;M2=modalM(5);C2=modalC(5);K2=modalK(5);q0=0;qdot0=0;load elcen_ns.dat;t=elcen_ns(:,1);ga=elcen_ns(:,2);pt=-M2*elcen_ns(:,2);P=[pt pt pt] ;Pt=(vec'*P')' ;R0 = Pt(:,2.db

프로그램소스| 2008.06.23| 9,000원| 조회(1,508)

MatLab, 다자유도, 구조동역학, 지진해석, 모드중첩법

미리보기

닫기
[구조동역학] matlab과 sap2000을 이용한 단자유도, 다자유도 지진해석 발표자료

구조 동역학 기말과제 발표(단자유도, 다자유도 지진해석)목 차1. 단자유도 ( SDOF )1) Time Domain Method2) Frequency Domain Method3) SAP2000목 차2. 다자유도 ( MDOF )1) Time Domain Method2) Frequency Domain Method3) Mode Super position4) SAP20001) Time Domain MethodDirect Integral Method ( Newmark beta method )SDOF (단자유도)2) Frequency Domain MethodMatlab 결과SDOF (단자유도)2) Frequency Domain MethodMatlab 결과SDOF (단자유도)2) Frequency Domain MethodMatlab 결과SDOF (단자유도)2) Frequency Domain MethodMatlab 결과SDOF (단자유도)- Compare with two resultMatlab 결과SDOF (단자유도)3) SAP2000SAP2000 모델링SDOF (단자유도)ㅇ3) SAP2000SAP2000 모델링SDOF (단자유도)3) SAP2000SAP2000 모델링SDOF (단자유도)3) SAP2000SAP2000 결과물SDOF (단자유도)3) SAP2000Matlab결과와 SAP2000 결과 비교SDOF (단자유도)3) SAP2000Matlab결과와 SAP2000 결과 비교SDOF (단자유도)3) SAP2000SAP2000 AnimationSDOF (단자유도)목 차2. 다자유도 ( MDOF )1) Time Domain Method2) Frequency Domain Method3) Mode Super position4) SAP20005) Tuned Mass Damper1) Time Domain MethodDirect Integral Method ( Newmark beta method )MDOF (다자유도)2) Frequency Domain MethodMatlab 결과MDOF (다자유도)2) Frequency Domain MethodMatlab 결과MDOF (다자유도)2) Frequency Domain MethodMatlab 결과MDOF (다자유도)2) Frequency Domain MethodFOR문으로 인한 H(w)의 오차MDOF (다자유도)2) Frequency Domain MethodMatlab 결과MDOF (다자유도)H(w)를 transpose 해줌으로써 정확한 결과를 얻을 수 있다.2) Frequency Domain MethodMatlab 결과MDOF (다자유도)H(w) 를 그대로 쓴 경우 transpose H(w) 를 쓴 경우3) Mode Super PositionMatlab 결과MDOF (다자유도)3) Mode Super PositionMatlab 결과MDOF (다자유도)- Compare with three result각각의 방법으로 해석하여 2층에서의 변위 비교MDOF (다자유도)4) SAP2000MDOF (다자유도)SAP2000 모델링4) SAP2000SAP2000 모델링MDOF (다자유도)4) SAP2000SAP2000 모델링MDOF (다자유도)4) SAP2000SAP2000 결과물MDOF (다자유도)Mode 1Mode 2Mode 34) SAP2000SAP2000 결과물MDOF (다자유도)4) SAP2000SAP2000 결과물MDOF (다자유도)4) SAP2000SAP2000 결과물MDOF (다자유도)4) SAP2000SAP2000 결과물MDOF (다자유도)4) SAP2000SAP2000 결과물MDOF (다자유도)4) SAP2000SAP2000 결과물MDOF (다자유도)4) SAP2000SAP2000 결과물MDOF (다자유도)5) Tuned Mass DamperTuned Mass Damper 매트랩 결과MDOF (다자유도)5) Tuned Mass DamperTuned Mass Damper 매트랩 결과MDOF (다자유도)5) Tuned Mass DamperTuned Mass Damper 매트랩 결과MDOF (다자유도)5) Tuned Mass DamperTuned Mass Damper 매트랩 결과 (TMD를 설치한 4층 빌딩)MDOF (다자유도)5) Tuned Mass DamperTuned Mass Damper 매트랩 결과 (TMD를 설치한 4층 빌딩)MDOF (다자유도)5) Tuned Mass DamperTuned Mass Damper 매트랩 결과 (4층 비교)MDOF (다자유도)감사합니다.{nameOfApplication=Show}

공학/기술| 2008.06.23| 45페이지| 9,000원| 조회(995)

SAP2000, 구조동역학, matrab, 지진해석, newmark beta m..

미리보기

닫기
[구조동역학] (matlab) 매틀랩을 이용한 단자유도 지진 해석 결과(single degree of freedom)

..FILE:conclusion.mclear allclose allclc%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Direct intergration Method중에 뉴마크법을 이용한 응답%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%m=1;k=200;xi=0.1;u0=0;udot0=0;load elcen_ns.dat;t=elcen_ns(:,1)';ga=elcen_ns(:,2);p=(-m*elcen_ns(:,2))';[u,udot,u2dot]=newmark_int(t,p,u0,udot0,m,k,xi);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Frquency Domain Method를 이용한 응답%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 물성치 입력m=1;xi=0.1;k=200;c=xi * 2*sqrt(m*k);%%% F(w)를 구하는 루틴dt = 0.02;T = 10;t = 0:dt:T ;df = 1/T ;fmax = 1/dt;f = 0.1:df:fmax;load elcen_ns.dat;t=elcen_ns(:,1)';ga=elcen_ns(:,2);ft=(-m*ga)';ft = [ft zeros(1,length(t)-length(ft))] ;Fw = fft(ft) * dt;%%% H(w)를 구하는 루틴w= 2*pi*f;Hw= 1./(k-w.^2*m + i*w*c);%%% X(w)를 구하는 루틴Xw = Hw.* Fw ;Xwdot = i*w.*Xw;Xw2dot = -w.*w.* Xw;%%% x축을 time의 함수로 표현 : Xw --> Xtfor k=(length(t))/2+2:length(t);Xw(k) = Xw(length(t)+2-k)' ;endxt = ifft(Xw) / dt;for k=(length(t))eleration')grid on..FILE:ELCEN_NS.DAT0.02 -0.001430.04 -0.011020.06 -0.010310.08 -0.008980.10 -0.009690.12 -0.012240.14 -0.014490.16 -0.013060.18 -0.011220.20 -0.008670.22 -0.008670.24 -0.013370.26 -0.017960.28 -0.019800.30 -0.016530.32 -0.014690.34 -0.011020.36 -0.008370.38 -0.004290.40 -0.006730.42 -0.013370.44 -0.019390.46 -0.020000.48 -0.006730.50 0.003060.52 0.014390.54 -0.005000.56 -0.013060.58 -0.014690.60 -0.020710.62 -0.026530.64 -0.033160.66 -0.031220.68 -0.017550.70 -0.020100.72 -0.016630.74 -0.016730.76 -0.006840.78 0.002550.80 0.015310.82 0.024080.84 0.025710.86 0.034290.88 0.047240.90 0.050200.92 0.042760.94 0.036630.96 0.027650.98 0.023981.00 0.034591.02 0.042041.04 0.054081.06 0.065201.08 0.074691.10 0.066531.12 0.061121.14 0.040821.16 0.040821.18 0.006431.20 -0.052551.22 -0.080311.24 -0.061531.26 -0.049391.28 -0.025511.30 -0.006021.32 0.013671.34 0.031431.36 0.050921.38 0.072451.40 0.101531.42 0.124391.44 0.156021.46 0.147861.48 0.117861.50 0.095411.52 0.091021.54 0.094.38 -0.123164.40 -0.055314.42 -0.039184.44 -0.031734.46 -0.114084.48 -0.169494.50 -0.251434.52 -0.206634.54 -0.187244.56 -0.134394.58 -0.097964.60 -0.033164.62 0.015714.64 0.083274.66 0.134594.68 0.185514.70 -0.005924.72 -0.017244.74 0.029084.76 0.045614.78 0.100314.80 0.145314.82 0.189084.84 0.250614.86 0.171944.88 -0.140824.90 -0.101944.92 -0.111124.94 -0.092554.96 -0.047864.98 -0.127555.00 -0.215415.02 -0.165005.04 -0.172655.06 -0.133275.08 -0.113375.10 -0.078885.12 -0.052045.14 -0.055515.16 -0.122455.18 -0.123375.20 -0.118165.22 -0.116845.24 -0.073165.26 -0.055715.28 0.006535.30 -0.082045.32 -0.166735.34 -0.087655.36 -0.098065.38 -0.040415.40 -0.015005.42 0.032555.44 0.066125.46 0.089395.48 0.048165.50 0.020205.52 -0.002765.54 0.029805.56 0.045415.58 0.080105.60 0.105415.62 0.137965.64 0.163885.66 0.189905.68 0.130715.70 0.065315.72 0.020825.74 0.032045.76 0.038065.78 0.050615.80 0.023985.82 -0.008575.84 -0.017145.86 -0.011535.88 -0.023375.90 -0.025315.92 -0.016025.94 -0.007045.968.78 -0.018578.80 -0.015008.82 0.008678.84 0.016638.86 0.005108.88 0.026948.90 0.059398.92 0.088478.94 0.122458.96 0.172968.98 0.113379.00 -0.112249.02 -0.037359.04 -0.045419.06 -0.024089.08 -0.097969.10 -0.066949.12 -0.060929.14 -0.068379.16 -0.056339.18 -0.002769.20 0.038579.22 0.109399.24 0.170319.26 0.096639.28 0.041639.30 0.068069.32 0.013479.34 -0.009699.36 -0.053069.38 -0.084399.40 -0.117559.42 -0.117359.44 -0.081949.46 -0.037659.48 0.002969.50 0.055619.52 0.120209.54 0.164299.56 -0.027559.58 0.003479.60 -0.005719.62 0.002049.64 0.014909.66 0.054809.68 0.081439.70 -0.020929.72 -0.060209.74 -0.017249.76 -0.017869.78 -0.002869.80 0.007559.82 0.038989.84 0.057869.86 0.076849.88 0.081739.90 0.060419.92 0.031029.94 0.002359.96 0.006539.98 -0.0414310.00 -0.04602..FILE:FDM.mclear ;clc ;close all%%% 물성치 입력m=1;xi=0.1;k=200;c=xi * 2*sqrt(m*k);%%% F(w)를 구하는 루틴dt = 0.02;T = 10;t = 0:dt:T ;df = 1/T ;fmax = 1/dt;f = 0.1:df:fmax;load elcen_ns.dat;t=elcen_ns(:,1)';ga=elcen_ns(:,2);ft=(-m*ga)';fmag(xt),'b')title('bf{Frequency Domain Method}','fontsize',10)xlabel('Time [second]')ylabel('imag{X(t)}')set(gcf, 'color', 'w')grid onsubplot(312) ; plot(t,imag(xtdot),'b')xlabel('Time [second]')ylabel('imag{Xdot(t)}')grid onsubplot(313) ; plot(t,imag(xt2dot),'b')xlabel('Time [second]')ylabel('imag{X2dot(t)}')grid on% figure ;% subplot(311) ; plot(t,abs(xt),'b') ;% subplot(312) ; plot(t,imag(xt),'b') ;% subplot(313) ; plot(t,real(xt),'b') ;..FILE:FDM1.jpg..FILE:FDM2.jpg..FILE:FDM3.jpg..FILE:FDM4.jpg..FILE:newmark.jpg..FILE:newmark_int.mfunction [u, udot, u2dot] = newmark_int(t,p,u0,udot0,m,k,xi)%Newmark's Direct Integration Method%--------------------------------------------------------------------------% Integrates a 1-DOF system with mass "m", spring stiffness "k" and damping% coeffiecient "xi", when subjected to an external load P(t).% Returns the displacement, velocity and acceleration of the system with% respect to an inertial frame of reference.%% SYNTAX% [u, udot, u2dot] = newmb

공학/기술| 2008.06.23| 6페이지| 12,000원| 조회(1,347)

MatLab, 구조동역학, Direct intergration, Frquenc..

미리보기

닫기