Skip to content

Elas force #279

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
victorkemp merged 2 commits into from
Jul 8, 2026
Merged

Elas force #279

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
eb33902
Changed from STRESS->FORCE to FORCE->STRESS to allow zero S
victorkemp Jul 7, 2026
a7daa61
Clarifying comments for STRESS() and FCONV()
victorkemp Jul 7, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 3 additions & 3 deletions Source/EMG/EMG3/ELAS1.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,11 +75,11 @@ SUBROUTINE ELAS1 ( OPT, WRITE_WARN )
ENDIF

! **********************************************************************************************************************************
! Calculate SE1 matrix for stress recovery.
! Calculate SE1 matrix for force recovery.

IF (OPT(3) == 'Y') THEN
SE1(1,I1,1) = K*FCONV(1)
SE1(1,I2,1) = -K*FCONV(1)
SE1(1,I1,1) = K
SE1(1,I2,1) = -K
ENDIF


Expand Down
2 changes: 2 additions & 0 deletions Source/LK9/L92/ELEM_STRE_STRN_ARRAYS.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -98,6 +98,8 @@ SUBROUTINE ELEM_STRE_STRN_ARRAYS ( STR_PT_NUM )

! **********************************************************************************************************************************
! Calc stresses for 1D elements
! Warning: For ELAS1/2/3/4, the STRESS() array contains force, not stress. This is because we must calculate stress from force, not
! the other way around.

IF ((TYPE(1:3) == 'BAR') .OR. (TYPE(1:4) == 'BUSH') .OR. (TYPE(1:4) == 'ELAS') .OR. (TYPE(1:3) == 'ROD') .OR. &
(TYPE(1:5) == 'USER1')) THEN
Expand Down
34 changes: 6 additions & 28 deletions Source/LK9/L92/OFP3_ELFE_1D.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,7 @@ SUBROUTINE OFP3_ELFE_1D ( JVEC, FEMAP_SET_ID, ITE, OT4_EROW )
USE FEMAP_ARRAYS, ONLY : FEMAP_EL_NUMS, FEMAP_EL_VECS
USE PARAMS, ONLY : OTMSKIP, PRTNEU
USE MODEL_STUF, ONLY : ANY_ELFE_OUTPUT, BUSH_CID, BUSH_VVEC, EDAT, ELAS_COMP, ELEM_LEN_12, ELEM_LEN_AB, EPNT, &
ETYPE, EID, ELMTYP, ELOUT, FCONV, METYPE, NUM_EMG_FATAL_ERRS, OFFDIS_GA_GB, OFFDIS_L, &
ETYPE, EID, ELMTYP, ELOUT, METYPE, NUM_EMG_FATAL_ERRS, OFFDIS_GA_GB, OFFDIS_L, &
PE_GA_GB, PEL, PLY_NUM, STRESS, TE, TE_GA_GB, TYPE, XEL
USE LINK9_STUFF, ONLY : EID_OUT_ARRAY, MAXREQ, OGEL
USE OUTPUT4_MATRICES, ONLY : OTM_ELFE, TXT_ELFE
Expand Down Expand Up @@ -188,11 +188,7 @@ SUBROUTINE OFP3_ELFE_1D ( JVEC, FEMAP_SET_ID, ITE, OT4_EROW )
CALL ELEM_STRE_STRN_ARRAYS ( 1 )
NDUM = 0
CALL CALC_ELEM_STRESSES ( 1, NDUM, 0, 'N', 'N' )
IF (FCONV(1) > ZERO) THEN ! ELAS engr force is stress/FCONV
OGEL(NUM_OGEL,1) = STRESS(1)/FCONV(1)
ELSE
OGEL(NUM_OGEL,1) = ZERO
ENDIF
OGEL(NUM_OGEL,1) = STRESS(1) ! ELAS engr force is stored in the stress array
! end elas
! ---------------------------------------------------------------------------------------------------------------
ELSE IF (ETYPE(J)(1:4) == 'BUSH') THEN
Expand Down Expand Up @@ -553,8 +549,6 @@ SUBROUTINE OFP3_ELFE_1D ( JVEC, FEMAP_SET_ID, ITE, OT4_EROW )
ENDIF
CALL DEALLOCATE_FEMAP_DATA

! For ELAS we need to calculate elem engr forces from the stresses since there is no "local" elem coord system

! elas1 ---------------------------------------------------------------------------------------------------------------------------
NDUM = 0
NUM_FROWS= 0
Expand All @@ -577,11 +571,7 @@ SUBROUTINE OFP3_ELFE_1D ( JVEC, FEMAP_SET_ID, ITE, OT4_EROW )
CALL ELMDIS
CALL ELEM_STRE_STRN_ARRAYS ( 1 )
CALL CALC_ELEM_STRESSES ( NCELAS1, NDUM, NUM_FROWS, 'N', 'Y' )
IF (FCONV(1) > 0.D0) THEN
FEMAP_EL_VECS(NUM_FROWS,1) = STRESS(1)/FCONV(1)
ELSE

ENDIF
FEMAP_EL_VECS(NUM_FROWS,1) = STRESS(1)
ENDIF
ENDDO
IF (NUM_FROWS > 0) THEN
Expand Down Expand Up @@ -611,11 +601,7 @@ SUBROUTINE OFP3_ELFE_1D ( JVEC, FEMAP_SET_ID, ITE, OT4_EROW )
CALL ELMDIS
CALL ELEM_STRE_STRN_ARRAYS ( 1 )
CALL CALC_ELEM_STRESSES ( NCELAS2, NDUM, NUM_FROWS, 'N', 'Y' )
IF (FCONV(1) > 0.D0) THEN
FEMAP_EL_VECS(NUM_FROWS,1) = STRESS(1)/FCONV(1)
ELSE

ENDIF
FEMAP_EL_VECS(NUM_FROWS,1) = STRESS(1)
ENDIF
ENDDO
IF (NUM_FROWS > 0) THEN
Expand Down Expand Up @@ -645,11 +631,7 @@ SUBROUTINE OFP3_ELFE_1D ( JVEC, FEMAP_SET_ID, ITE, OT4_EROW )
CALL ELMDIS
CALL ELEM_STRE_STRN_ARRAYS ( 1 )
CALL CALC_ELEM_STRESSES ( NCELAS3, NDUM, NUM_FROWS, 'N', 'Y' )
IF (FCONV(1) > 0.D0) THEN
FEMAP_EL_VECS(NUM_FROWS,1) = STRESS(1)/FCONV(1)
ELSE

ENDIF
FEMAP_EL_VECS(NUM_FROWS,1) = STRESS(1)
ENDIF
ENDDO
IF (NUM_FROWS > 0) THEN
Expand Down Expand Up @@ -679,11 +661,7 @@ SUBROUTINE OFP3_ELFE_1D ( JVEC, FEMAP_SET_ID, ITE, OT4_EROW )
CALL ELMDIS
CALL ELEM_STRE_STRN_ARRAYS ( 1 )
CALL CALC_ELEM_STRESSES ( NCELAS4, NDUM, NUM_FROWS, 'N', 'Y' )
IF (FCONV(1) > 0.D0) THEN
FEMAP_EL_VECS(NUM_FROWS,1) = STRESS(1)/FCONV(1)
ELSE

ENDIF
FEMAP_EL_VECS(NUM_FROWS,1) = STRESS(1)
ENDIF
ENDDO
IF (NUM_FROWS > 0) THEN
Expand Down
9 changes: 4 additions & 5 deletions Source/LK9/L92/ONE_D_STRESS_OUTPUTS.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -33,9 +33,8 @@ SUBROUTINE ONE_D_STRESS_OUTPUTS ( SIZE_ALLOCATED, NUM1, NUM_FEMAP_ROWS, WRITE_OG
USE PENTIUM_II_KIND, ONLY : BYTE, LONG, DOUBLE
USE IOUNT1, ONLY : ERR, F06
USE SCONTR, ONLY : BLNK_SUB_NAM, FATAL_ERR
USE TIMDAT, ONLY : TSEC
USE CONSTANTS_1, ONLY : ZERO
USE MODEL_STUF, ONLY : STRESS, TYPE, ZS
USE MODEL_STUF, ONLY : STRESS, TYPE, ZS, FCONV
USE LINK9_STUFF, ONLY : MSPRNT, OGEL
USE FEMAP_ARRAYS, ONLY : FEMAP_EL_VECS
USE PARAMS, ONLY : PRTNEU
Expand Down Expand Up @@ -184,11 +183,11 @@ SUBROUTINE ONE_D_STRESS_OUTPUTS ( SIZE_ALLOCATED, NUM1, NUM_FEMAP_ROWS, WRITE_OG
FATAL_ERR = FATAL_ERR + 1
CALL OUTA_HERE ( 'Y' )
ENDIF
OGEL(NUM1,1) = STRESS(1)
OGEL(NUM1,1) = FCONV(1)*STRESS(1) ! STRESS() contains force and FCONV(1) is S.
ENDIF
IF (WRITE_NEU .AND. (WRITE_FEMAP == 'Y')) THEN
FEMAP_EL_VECS(NUM_FEMAP_ROWS,1) = STRESS(1)
FEMAP_EL_VECS(NUM_FEMAP_ROWS,2) = STRESS(1)
FEMAP_EL_VECS(NUM_FEMAP_ROWS,1) = FCONV(1)*STRESS(1)
FEMAP_EL_VECS(NUM_FEMAP_ROWS,2) = FCONV(1)*STRESS(1)
ENDIF

ELSE IF (TYPE == 'ROD ') THEN ! ROD1 elements
Expand Down
2 changes: 2 additions & 0 deletions Source/Modules/MODEL_STUF.f03
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1421,6 +1421,7 @@ MODULE MODEL_STUF

REAL(DOUBLE) :: FCONV(3) = (/(ZERO, I=1,3)/)
! Array of constants to convert stresses to engr forces for elem
! Except CELAS where it converts force to stress

REAL(DOUBLE) :: HBAR = ZERO
! For the quad elem, the dist from the mean plane to the G.P.'s
Expand Down Expand Up @@ -1511,6 +1512,7 @@ MODULE MODEL_STUF

REAL(DOUBLE) :: STRESS(9) = (/(ZERO, I=1,9)/)
! Array of elem stresses for one S/C
! Except CELAS where it is force, not stress.

REAL(DOUBLE) :: TE(3,3) = RESHAPE ( (/(ZERO, I=1,3*3)/), (/3,3/) )
! Coord system transformation matrix (3 x 3) such that UEL = TE*UEB
Expand Down
Loading