% Mondai51.m
% 
% 5章の問題（１）全体枠法による要素力計算
% 
% Copyright (C) 2007 清水 信行, 曽我部 潔, 藤川 猛 All Rights Reserved.

% =================================================================================
% プログラムの使用条件
% 1. 本プログラムおよびデータを商用で，頒布してはならない．
% 2. 本プログラムおよびデータを著作者の許可なく営利目的で利用してはならない．
% 3. 本プログラムおよびデータの改変は，営利目的でない限り自由である．
% 4. 本プログラムおよびデータを使用することにより生じた損害などについては，著者，
%    プロクラム著作者，(社)日本機械学会，および(株)コロナ社は一切の責任を負わない．
% 5. 本プログラムの著作権は，本書編著者ならびにプログラム著作者に帰す．
% =================================================================================

com='Problem 5.1 Total frame method'
L0=0.5*sqrt(2);EA=1.2e8;EI=4e3;                     %データ
Qi=[0.1; 0.1;pi/4]  ;Qj=[0.59;0.61;0.001+pi/4];     %座標    
Qij=Qj(1:2)-Qi(1:2) ;L=norm(Qij)                    %長さ
alpha=atan(Qij(2)/Qij(1))                           %姿勢角
epsL=L-L0
epsi=Qi(3)-alpha
epsj=Qj(3)-alpha                                    %変形量                                   
fxj=(EA/L0)*epsL*cos(alpha)+6*EI/(L0*L)*(epsi+epsj)*sin(alpha)
fyj=(EA/L0)*epsL*sin(alpha)-6*EI/(L0*L)*(epsi+epsj)*cos(alpha)
Mj=2*EI/L0*(epsi+2*epsj)                            %要素力
T=(EA/L0)*epsL                                      %軸力
V=6*EI/(L0*L)*(epsi+epsj)                           %横力
%Sol=[fxj;fyj;Mj];
