SDOF_SFS_MODEL.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	sdof_sfs_model (sid, gr_acc, direction)
	Response computation for 4DOF model.

Function/Subroutine Documentation

sdof_sfs_model()

subroutine sdof_sfs_model	(	integer*4	sid,
		real*8, dimension(3)	gr_acc,
		integer*4	direction
	)

Response computation for 4DOF model.

Author

Aline Herlin

Date

February, 2021

Version

1.0

Parameters

[in]	sID	oscillator ID
[in]	gr_acc	(negative) ground acceleration
[in]	direction	motion direction

Definition at line 27 of file SDOF_SFS_MODEL.f90.

 
  use speed_sci
 
  implicit none
 
  integer*4 :: sID, direction
  real*8, dimension(3) :: gr_acc
  real*8, dimension(4) :: tn, tn1, tn2, tn3, tn4, temp, du, df, vec, a_old, v_old, u_old, f_old
 
  tn = 0; tn1 = 0; tn2 = 0; tn3 = 0; tn4 = 0; temp = 0; du = 0; df = 0; vec = 0
  a_old = sys(sid)%a(:,direction); v_old = sys(sid)%v(:,direction); u_old = sys(sid)%u(:,direction); f_old = sys(sid)%f(:,direction)
 
  tn1(1) = sys(sid)%Ms(1,1)*gr_acc(direction)
  tn1(2) = sys(sid)%Mf*gr_acc(direction)
  tn1(4) = (sys(sid)%Ms(1,1) + sys(sid)%Mf)*gr_acc(3)
 
  temp = (1 - 2*sys(sid)%beta_newmark)/(2*sys(sid)%beta_newmark)*a_old + 1/(sys(sid)%beta_newmark*sys(sid)%dt)*v_old
  tn2 = matmul(sys(sid)%MAT_M,temp)
 
  temp = (sys(sid)%gamma_newmark - 2*sys(sid)%beta_newmark)/(2*sys(sid)%beta_newmark)*sys(sid)%dt*a_old + &
          (sys(sid)%gamma_newmark - sys(sid)%beta_newmark)/(sys(sid)%beta_newmark)*v_old
  tn3 = matmul(sys(sid)%MAT_C,temp)
 
  tn4 = matmul(sys(sid)%MAT_KS,u_old)
 
  tn = tn1 + tn2 + tn3 - tn4 - f_old
 
  du = matmul(sys(sid)%MAT_MCinv,tn)
 
  df = matmul(sys(sid)%MAT_KI,du)
 
  sys(sid)%u(:,direction) = u_old + du
 
  sys(sid)%f(:,direction) = f_old + df
 
  vec = matmul(sys(sid)%MAT_KS,sys(sid)%u(:,direction))
  sys(sid)%fs(direction) = vec(1)
  sys(sid)%fb(direction) = sys(sid)%f(2,direction) + sys(sid)%Mf*(gr_acc(direction) - a_old(2))
 
  sys(sid)%v(:,direction) = v_old + (2*sys(sid)%beta_newmark - sys(sid)%gamma_newmark)/(2*sys(sid)%beta_newmark)*sys(sid)%dt*a_old + &
    sys(sid)%gamma_newmark/(sys(sid)%beta_newmark*sys(sid)%dt)*(du - v_old*sys(sid)%dt)
 
  sys(sid)%a(:,direction) = (2*sys(sid)%beta_newmark - 1)/(2*sys(sid)%beta_newmark)*a_old + &
    1/(sys(sid)%beta_newmark*sys(sid)%dt2)*(du - v_old*sys(sid)%dt)
 
  sys(sid)%IntForce(1,direction) = sys(sid)%f(2,direction)
 
  ! UnComment this after verifying
  !sys(sID)%SDOFItF(3) = sys(sID)%f(4,direction)
 
  return

References speed_sci::sys.