How to create an I-beam for the FEMAP property using EXCEL VBA

Sub CreatIBeam

  'Creating a Femap model object
  Dim f As Object
  Set f = GetObject(,"femap.model")
  
  'Creating Property Objects
  Dim pr As Object
  Set pr = f.feProp
  
  'Where return codes are stored
  Dim rc As Long

  'Input of geometry（For I-shaped cross section)
  Dim shape(5) As Double
  shape(0) = 400 'Input the height. In this case, the height is 400.
  shape(1) = 400 'Top flange width. In this case, the width is 400.
  shape(2) = 400 'Bottom flange width. In this case, the width is 400.
  shape(3) = 16  'Top flange thickness. In this case, it is 16.
  shape(4) = 16  'Bottom flange thickness. In this case, it is 16.
  shape(5) = 12  'Web thickness. In this case, it is 12.
  
  'Set the orientation direction. 0: right, 1: up, 2: left, 3: down
  Dim orient As Long
  orient = 0
  
  pr.title = "400X400X12/16 I"  'Input title. In this case, “400X400X12/16 I”.
  pr.matlID = 1  'Enter the material ID. In this case, it is 1.
  pr.Type = 5    'Enter the type. For beams, use 5.
  
  'Creating the property
  rc = pr.ComputeStdShape(9, shape, orient, 0, True, True, False) 'Arguments: shape ID, cross-section definition parameters, orientation, cross-section evaluation method, shear center offset calculation, torsional constant calculation, stress extraction location
  pr.Put(1)  'Register the created property with property ID 1.
  
  'Releasing objects
  Set pr = Nothing
  Set f = Nothing
  
  if rc = 0 then
    MsgBox "Failed to create property"
  else
    MsgBox "Property created successfully"
  End If
  
End Sub

Sub CreatIBeam

  'Create Femap model object
  Dim f As Object
  Set f = GetObject(,"femap.model")
  
  'Create property object
  Dim pr As Object
  Set pr = f.feProp
  
  'Storage for the return code
  Dim rc As Long

  'Storage for the return code
  Dim shape(5) As Double
  shape(0) = Cells(2,4) 'Input height
  shape(1) = Cells(2,5) 'Top flange width
  shape(2) = Cells(2,6) 'Bottom flange width
  shape(3) = Cells(2,7) 'Top flange thickness
  shape(4) = Cells(2,8) 'Bottom flange thickness
  shape(5) = Cells(2,9) 'Web thickness
  
  'Set the orientation direction. 0: right, 1: up, 2: left, 3: down
  Dim orient As Long
  orient = Cells(2,10)
  
  pr.title = Cells(2,3)  'Input title
  pr.matlID = Cells(2,2) 'Input material ID
  pr.Type = 5            'Input material ID
  
  'Create the property
  rc = pr.ComputeStdShape(9, shape, orient, 0, True, True, False) 'Arguments: shape ID, cross-section definition parameters, orientation, cross-section evaluation method, shear center offset calculation, torsional constant calculation, stress extraction location
  pr.Put(Cells(2,1))  'Register the created property with the specified ID number
  
  'Release objects
  Set pr = Nothing
  Set f = Nothing
  
  if rc = 0 then
    MsgBox "Failed to create property"
  else
    MsgBox "Property created successfully"
  End If
  
  End Sub

Sub GetBeamProperty()

  On Error GoTo OUT1
  'Generate FEMAP model object
  Dim f As Object
  Set f = GetObject(,"femap.model")
  
  On Error GoTo OUT2
  'Generate property object
  Dim pr As Object
  Set pr = f.feProp
  
  Dim rc As Long 'Return code
  
  pr.Get(Int(Cells(1,2))) 'Retrieve the property ID entered in cell B1
  
  Cells(2,2) = pr.matlID  'Retrieve material ID
  Cells(3,2) = pr.title   'Retrieve title
  Cells(4,2) = pr.Type    'Retrieve type (should be 5 for beam elements)
 
  'If the retrieved property is a beam, retrieve the following items
  If pr.Type = 5 Then
    Cells(5,2) = pr.pval(40)   'General section height
    Cells(6,2) = pr.pval(41)   'General section top flange width
    Cells(7,2) = pr.pval(42)   'General section bottom flange width 
    Cells(8,2) = pr.pval(43)   'General section top flange thickness  
    Cells(9,2) = pr.pval(44)   'General section bottom flange thickness  
    Cells(10,2) = pr.pval(45)  'General section web thickness  
    Cells(11,2) = pr.pval(0)   'Cross-sectional area at end point A
    Cells(12,2) = pr.pval(5)   'Effective cross-sectional area in Y direction at end point A 
    Cells(13,2) = pr.pval(6)   'Effective cross-sectional area in Z direction at end point A
    Cells(14,2) = pr.pval(1)   'Moment of inertia I1 at end point A
    Cells(15,2) = pr.pval(2)   'Moment of inertia I2 at end point A   
    Cells(16,2) = pr.pval(16)  'Neutral axis offset Y at end point A
    Cells(17,2) = pr.pval(17)  'Neutral axis offset Z at end point A
  Else
    MsgBox "This property  " & Cells(1,2) & " is not a beam."
  End If  
    
  'Release objects
  Set pr = Nothing
  Set f = Nothing
  
  Exit Sub
OUT1:
  MsgBox "Please start FEMAP"
  Exit Sub
OUT2:
  MsgBox "Unknown error"
End Sub