Canic
Civil/Environmental
- Aug 22, 2023
- 2
Please help,
Objective of the project: is to calculate vehicle-structure interaction. So technically, the DLOAD Force will be depended on the displacement of the node from the previous increment.
Can anybody know how to obtain the global node on the following script of DLOAD, so that I can make a common block for the URDFIL to extract the displacement and velocity of the global node?
Thank you in advance
Best Regards
Cank
Objective of the project: is to calculate vehicle-structure interaction. So technically, the DLOAD Force will be depended on the displacement of the node from the previous increment.
Code:
SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,
1 COORDS,JLTYP,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION TIME(2), COORDS (3)
CHARACTER*80 SNAME
C
X0 = 0.25! Initial position of the crane center
GB = 2.6! Initial position of the crane center
TG = 1! Initial position of the crane center
BB = 0.5! Initial position of the crane center
V = 28! m/s speed of crane
XL = 0.5 !lentgth of crane
W = F!N crane weight
C
! Position of the crane center
Xcenter = X0 + V*TIME(1)
Ycenter1 = GB/2 - TG/2
Ycenter2 = GB/2 + TG/2
! Calculating F
X = COORDS (1)
Y = COORDS (2)
F = 0.0
if (X>(Xcenter-XL/2.0) .and. X < (Xcenter+XL/2.0) .and. Y > (Ycenter1-BB/2.0) .and. Y < (Ycenter1+BB/2.0)) then
F = 170000/BB/BB
write(*,*) NOEL,F
End if
if (X>(Xcenter-XL/2.0) .and. X < (Xcenter+XL/2.0) .and. Y > (Ycenter2-BB/2.0) .and. Y < (Ycenter2+BB/2.0)) then
F = 170000/BB/BB
End if
RETURN
END
Can anybody know how to obtain the global node on the following script of DLOAD, so that I can make a common block for the URDFIL to extract the displacement and velocity of the global node?
Thank you in advance
Best Regards
Cank