Excuse me for my very tiny experience in APDL commands. I have read the command and I think I'm understanding each row from the article, but I can't determine where should I input the value of ARG1, ARG2, ARG3 and which row is the star of the command?
! Commands inserted into this file will be executed just prior to the ANSYS SOLVE command.
! These commands may supersede command settings set by Workbench.
! Active UNIT system in Workbench when this object was created: Metric (mm, t, N, s, mV, mA)
! Generate an elastic foundation in normal and both tangent directions
! on Named Selection "Elastic_Here", a set of nodes on face(s).
! These commands have been tested on a 3D solid model only.
! ARG1 is the NORMAL stiffness (Force per Unit Length per Unit Area)
! ARG2 is TANGENTIAL stiffness (Force per Unit Length per Unit Area)
! ARG3 is Pinball Radius, in solver units
! ARG1 & ARG2 must be in the solver units!
! If ARG2 is blank or zero, it is set to ARG1
! ARG1 must be non-zero.
! If ARG3 is zero, then zero will be used, which should activate default
ARG1 for Normal Stiffness on XYZ Elastic Foundation must be positive
/COM,######## ARG2 was made equal to ARG1 ########
*get,nodemax,NODE,,NUM,MAX ! highest node number in model
cmsel,s,Elastic_Here ! nodes of the component "Elastic_Here"
esln ! select contacting elements
! undselect surface effect, contact, MPC and beam elements
*get,maxtype,ETYP,,NUM,MAX ! highest element type
*get, maxmat,MAT,,NUM,MAX ! highest material type
*get,maxreal,RCON,,NUM,MAX ! highest real constant
! Set maxtype, maxmat and maxreal to the highest value of all three
! Create required element types and real constant
ET,maxtype+1,CONTA174,,1,,0,3 ! Pure Penalty contact algorithm (stiffness)
KEYOPT,maxtype+1,9,1 ! Exclude initial geometrical penetration or gap and offset
KEYOPT,maxtype+1,12,5 ! Bonded Always
ET,maxtype+2,TARGE170,,1 ! Constraints by user
R,maxreal+1,0,0,-ARG1,,,-abs(ARG3) ! FKN Absolute Number, ARG3=Absolute Pinball Radius
RMODIF,maxreal+1,12,-ARG2 ! FKT as Absolute Number
TYPE,maxtype+1 ! CONTA174 elements
ESURF ! Mesh CONTA174 over underlying element faces
esln,r,1 ! Select only these CONTA174 elements
esel,r,ename,,174 ! Ensure no other elements
! Make a copy of the currently selected nodes
NGEN,2,(nodemax-current_nodemin)+1,ALL,,,0,0,0 ! Copy of nodes at same location
EGEN,2,(nodemax-current_nodemin)+1,ALL,,,0,1,0 ! Copy elements, increment TYPE by 1
esel,r,type,,maxtype+2 ! Select these new TARGE170 elements
ENSYM,0,,0,ALL ! Reverse TARGE170 elements to face contacts
nsle ! Select nodes on these target elements
d,all,all ! Constrain all nodes on target elements
! Continue with the analysis