how can i get shear force & bending moment diagram?

[ufsaa2022]

hello :)
I would like to get results on Shear Force and Bending Moment.
Below are the APDL commands and their results.

/PREP7
MP,PRXY,,0.3
ET,1,BEAM188
KEYOPT,1,9,3         
KEYOPT,1,3,3         
SECT,1,BEAM,I
W_F=1.048394965
W_W=0.6856481
SECD,15,15,28+(2*W_F),W_F,W_F,W_W
MP,EX,1,30E6

K,1                         
K,2,420,1

L,1,2

ESIZE,10
LMESH,ALL

D,1,ALL,ALL   
FK,2,FY,-10000   

ALLSEL

/SOLU
SOVLE

/POST1
SET, LAST
ETABLE,MZI,SMISC,3
ETABLE,MZJ,SMISC,16

PLLS,MZI,MZJ,1,0,0

I would like to obtain a graph like this using PyMAPDL.
Is there a way to do this?

[mikerife]

Hi @ufsaa2022
Yes you would just need to include a 'show' command to redirect the plot from the terminal to a file. Much like you would running MAPDL in batch mode. So the PyMAPDL version would be (assuming you used the standard conventions):

mapdl.show('png')
mapdl.plls('mzi', 'mzj', 1, 0, 0)
mapdl.show()

Mike

[ufsaa2022]

Thank you for your reply.
Is showing using pyvista not supported yet?

[germa89]

@ufsaa2022 not natively at the moment. You can open an issue to request it.

[ufsaa2022]

@germa89 Thank you :)