!オンサイト木橋の理論値と誤差、kの値、ccxのkを出すプログラム
	   
	   
	   implicit real*8(a-h, o-z)
	   dimension a(10),g(10),v(10),x(10)
         n=10
         b=0.858!断面全体のb
         b1=0.849
		 b2=0.840
		 t1=0.12!木材の幅
		 t2=0.009!鋼板の幅
	     h=0.50
	     h1=0.26
		 h1=0.5-0.12
		 h2=0.26
	   	 ell=3.0
		 ell=11.5
	   	 ell2=ell*ell
	   	 ell3=ell*ell*ell
         e=6.921d9!ヤング率
		 e=5.6d9
		 c=e/15!せん断弾性係数
		 e2=206.0d9
		 c2=e2/2.6
         p=5.0d4!荷重
		 q=b1*t1/h1/t2
		 q2=q*q
		 q3=q*q*q
		 al=e2/e
		 al2=al*al
		 al3=al*al*al
		 r=b1/h1
		 r2=r*r
		 o=20*(al+3*q)*(al+3*q)/(e/c*(60*q2*r2+60*al*q*r2)+e*2.6/e2*(180*q3+300*al*q2+144*al2*q+24*al3)&
		  &+0.3*(30*q2+6*al*q-4*al2))!k*
		  !
		 og=20*(al+3*q)*(al+3*q)/&
		 &(e/c*(60*q2*r2+60*al*q*r2)&
		 &+e*2.6/e2*(180*q3+300*al*q2+144*al2*q+24*al3)&
		  &+0.3*(30*q2*r2-50.*al*q*r2)&
		  &+0.3*(30*q2+6*al*q-4*al2))!k*
		  print*,'k*g=',og
		 s=t1*t1*14!断面積木材
		 s2=t2*h*2!断面積鋼板
		 print*,'Aw=',s
		 print*,'As=',s2
		 sall=s+s2!断面積
		 f=b2*h*h*h/12-b2*h1*h1*h1/12!断面二次モーメント
		 f=b2*(h**3-h2**3)/12.
		 f2=2*t2*h*h*h/12!断面二次モーメント
		 f3=f2-2.*t2*0.174**3/12.
		 ez=(f2*e2+f*e)/(f+f2)!全体で見た時のE
		 o2=o*ez*(s+s2)!k
		 ea=e*s+e2*s2
		 print*,'EA合成=',ea
		 eak=og*ea
		 print*,'k*EA=',eak
		 print*,'Iw=',f
		 print*,'Is=',f3
		 print*,'EwIw=',e*f
		 print*,'EsIs=',e2*f3
         eig=e*f+e2*f3
		 print*,'EI合成=',eig
		 eis=1./(1./eig+12./(eak*ell**2))
		 print*,'EI(せん断考慮)=',eis
		 onaka5=20*(al+3*q)*(al+3*q)
         onaka2=e/c2*(180*q3+300*al*q2+144*al2*q+24*al3)
         onaka21=180*q3
         onaka22=300*al*q2
         onaka23=144*al2*q
         onaka24=24*al3
         onaka1=e/c*(60*q2*r2+60*al*q*r2)
         onaka4=0.3*(30*q2+6*al*q-4*al2)
		 onaka1k=15*(60*(0.849*0.12/(h-0.24)/0.009)*(0.849*0.12/(h-0.24)/0.009)*(0.849/(h-0.24))*(0.849/(h-0.24))&
		 &+60*(206/6.921)*(0.849*0.12/(h-0.24)/0.009)*(0.849/(h-0.24))*(0.849/(h-0.24)))
		 h3=((20*(206/6.921+3*0.849*0.12/(h-0.24)/0.009)**2)/&
		 &(15*(60*(0.849*0.12/(h-0.24)/0.009)**2*(0.849/(h-0.24)**2)&
		 &+60*(206/6.921)*(0.849*0.12/(h-0.24)/0.009)*(0.849/(h-0.24))**2)&
		 &+6.921*2.6/206*(180*(0.849*0.12/(h-0.24)/0.009)**3&
		 &+300*(206/6.921)*(0.849*0.12/(h-0.24)/0.009)**2&
		 &+144*(206/6.921)**2*(0.849*0.12/(h-0.24)/0.009)+24*(206/6.921)**3)&
		 &+0.3*(30*(0.849*0.12/(h-0.24)/0.009)**2&
		 &+6*(206/6.921)*(0.849*0.12/(h-0.24)/0.009)-4*(206/6.921)**2)))
	     a(1)=3.8675d-3
	     a(2)=7.3529d-3
	     a(3)=9.6036d-3
	     a(4)=1.0459d-2
	     a(5)=1.0723d-2
	     a(6)=1.0803d-2
	     a(7)=1.0828d-2
	     a(8)=1.0836d-2
	     a(9)=1.0839d-2
	     a(10)=1.8846d-2
	   

	     v(1)=p*ell/(a(1)-p*ell3/3/(e*f+e2*f2))
	     v(2)=p*ell/(a(2)-p*ell3/3/(e*f+e2*f2))
	     v(3)=p*ell/(a(3)-p*ell3/3/(e*f+e2*f2))
	     v(4)=p*ell/(a(4)-p*ell3/3/(e*f+e2*f2))
	     v(5)=p*ell/(a(5)-p*ell3/3/(e*f+e2*f2))
	     v(6)=p*ell/(a(6)-p*ell3/3/(e*f+e2*f2))
	     v(7)=p*ell/(a(7)-p*ell3/3/(e*f+e2*f2))
	     v(8)=p*ell/(a(8)-p*ell3/3/(e*f+e2*f2))
	     v(9)=p*ell/(a(9)-p*ell3/3/(e*f+e2*f2))
	     v(10)=p*ell/(a(10)-p*ell3/3/(e*f+e2*f2))

        omagetawami=p*ell3/3/(e*f+e2*f2)
		 rironti=ell*ell*ell*p/3/(e*f+e2*f2)+ell*p/o2
		 rironti2=ell*ell*ell*p/3/(e*f+e2*f2)+ell*p/(15*o*c*s)
		 print*,h3
		 print*,s
		 print*,s2
         print*,'dannmennseki='!k
		 print*,sall
		 print*,'Iw=',f
		 print*,'Is=',f2
		 snaka=1.8*10e-3*0.174**3/12.
		 print*,'snaka=',snaka
		 print*,'Is(soto)=',f2-snaka
         print*,'k*=',o!k*
		 !print*,o
         print*,'k*AEz=',o2!k
		 !print*,o2
         print*,'Ez='!k
		 print*,ez
		 !print*,onaka21
		 !print*,onaka22
		 !print*,onaka23
		 !print*,onaka24
		 print*,onaka1
		 print*,onaka1k
		 !print*,onaka2
		 !print*,onaka4
		 !print*,onaka5
		! print*,r
		! print*,r2
		! print*,q
		! print*,q2
		! print*,q3
         print*,'rironti='!rironti
		 print*,rironti
         print*,'rironti2='!rironti
		 print*,rironti2
         print*,'magetawami='!rironti
		 print*,omagetawami
		 
         g(1)=(v(1)/o2-1)*100
         g(2)=(v(2)/o2-1)*100
         g(3)=(v(3)/o2-1)*100
         g(4)=(v(4)/o2-1)*100
         g(5)=(v(5)/o2-1)*100
         g(6)=(v(6)/o2-1)*100
         g(7)=(v(7)/o2-1)*100
         g(8)=(v(8)/o2-1)*100
         g(9)=(v(9)/o2-1)*100
         g(10)=(v(10)/o2-1)*100


         x(1)=(a(1)/rironti-1)*100
         x(2)=(a(2)/rironti-1)*100
         x(3)=(a(3)/rironti-1)*100
         x(4)=(a(4)/rironti-1)*100
         x(5)=(a(5)/rironti-1)*100
         x(6)=(a(6)/rironti-1)*100
         x(7)=(a(7)/rironti-1)*100
         x(8)=(a(8)/rironti-1)*100
         x(9)=(a(9)/rironti-1)*100
         x(10)=(a(10)/rironti-1)*100
        
		
		
	
	
         print*,'ccxtawami='!ccxtawami
		 print*,(a(i),i=1,n)
		 print*,'ccxk*AEz='!ccxk
		 print*,(v(i),i=1,n)
		 print*,'kAEzsoutaigosa1'!k相対誤差
		 print*, (g(i),i=1,n)
		 print*,'tawamisoutaigosa1'!相対誤差１
		 print*, (x(i),i=1,n)
		 end
		 
