Axisymmetric 2D pipe x-section with multiple loads

Is this possible to solve using pycalculix? Axisymmetric 2D with multiple loads at the same time. I’m looking for some sort of python code example for similar task as shown in the image. I need to make a simple programmable solution where user input OD, ID, Elastic Modulus, Loads, and etc. And the code does solve and visualise von-mises and strain preferably in bokeh library or other. Running localhost with Flask.
Is anyone can share a similar setup if there is some :slight_smile: