Welcome to Adams/Solver Subroutines - Kxcad.net

SUBROUTINE SFOSUB (ID, TIME, PAR, NPAR, DFLAG, IFLAG, VALUE)
INTEGER ID
DOUBLE PRECISION TIME
DOUBLE PRECISION PAR( * )
INTEGER NPAR
LOGICAL DFLAG
LOGICAL IFLAG
DOUBLE PRECISION VALUE
DOUBLE PRECISION DM, VR, V(3)
INTEGER IPAR(2)
LOGICAL ERRFLG, PARFLG
IPAR(1)=12
IPAR(2)=9
CALL ADAMS_NEEDS_PARTIALS(PARFLG)
CALL SYSFNC('DM', IPAR, 2, DM, ERRFLG )
CALL SYSFNC('VR', IPAR, 2, VR, ERRFLG )
IF(PARFLG) THEN
CALL IMPACT(400.D0-DM,-VR,100D0,1.D2,1.5D0, 1D1,1.D0,1,V,ERRFLG)
CALL SYSPAR('DM', IPAR, 2, V(1), 1, ERRFLG )
CALL SYSPAR('VR', IPAR, 2, V(2), 1, ERRFLG )
ENDIF
Welcome to Adams/Solver Subroutines
CALL IMPACT(400-DM,-VR,100D0,100.0D0,1.5D0,10D0,1.0D0,0,V,ERRFLG)
VALUE=-V(1)
CALL ERRMES(ERRFLG,'ERROR FOR SFORCE',ID,'STOP')
RETURN
END
Using FIESUB with SYSPAR
This example, involving FIESUB, shows how SYSPAR has closed the functionality gap between
GFOSUB and FIESUB. In the FORTRAN 77 version of Adams/Solver, the ability to define partial
derivatives of the FIELD element in a FIESUB user subroutine gives it a computational advantage over
the GFORCE force element and a GFOSUB user subroutine. The FIESUB user subroutine is defined as
follows:
CALL FIESUB(ID, TIME, PAR, NPAR, DISP, VELO,
DFLAG, IFLAG, FIELD, DFDDIS, DFDVEL)
where DISP and VELO are the displacement and velocity of the I and J markers of the field element that
are input to FIESUB. FIESUB computes the 6D value of the force, FIELD. Also, when DFLAG is
TRUE, you must provide the partial derivatives with respect to displacement, DFDDIS, and the partial
derivatives with respect to velocity, DFDVEL. Learn more about FIELD (C++ or FORTRAN).
115