July 31, 2023 at 5:36 pm
Fraser Shields
Subscriber

Thanks for the help, I've managed to improve my script slightly, however I am now running into an error "Element type 216 requires the definition of TB,PM with the activation of pressure DOF (KEYOPT(12)). ". My script is as follows:

/PREP7
 
! Define an element type with the desired real constant (CPT216)
ET, 9991, CPT216
 
! Select the solid elements you want to modify (SOLID186)
ESEL, S, TYPE,, 13, 57
! Modify the selected elements to use the specified element type (9991)
EMODIF, ALL, TYPE, 9991
 
! Define proelastic material for solid elements
MP,EX,1,0.001 ! YOUNG'S MODULUS
MP,NUXY,1,0.0 ! POISSON'S RATIO
FPX=0.5 ! PERMEABILITY
ONE=1.0
TB,PM,1,,,PERM
TBDATA,1,FPX,FPX,FPX
TB,PM,1,,,BIOT ! BIOT COEFFICIENT
TBDATA,1,ONE
 
! Define linear elastic material for shell elements 
MP,EX,2,0.1 ! YOUNG'S MODULUS
MP,NUXY,2,0.0 ! POISSON'S RATIO
 
 
! Select the solid elements you want to modify and assign material 1
ESEL, S, TYPE,, 9991
MPCHG, 1, ALL
 
! Select the shell elements you want to modify and assign material 2
ESEL, S, TYPE,, 1, 12
MPCHG, 2, ALL
 
! Select all elements
ESEL, ALL
 
! Continue with the solution phase
/SOLU