FEM for section properties of a beam cross section

When I calculate section properties for a beam, I usually run AutoCAD's MASSPROP command on the section and then do any additional analysis in MathCAD. I'm often doing composite beams, beams with custom shapes, etc. so I can't rely on magic formulas for specific shapes. There are a few properties that I haven't been able to calculate as easily, such as the torsional constant, warping constant, plastic section modulus about plastic section centroid, and monosymmetry constant. These are found using finite element analysis. I have a python API written by Robbie van Leeuwen called SECTIONPROPERTIES, and it's great but my colleagues don't know how to use Python. So with that context, my question is "How can I use MathCAD Prime to cut a section into discrete elements and run a finite element analysis or solve differential equations o?

Here's the DE for the monosymmetry constant:

Here's the St. Venant torsion constant:

The warping constant, from Galambos, 1968: