sewkokot / opsvis

OpenSeesPy postprocessing and plotting module
GNU General Public License v3.0
35 stars 21 forks source link

Question of Simple 2D Analysis #40

Closed hiroaki415 closed 1 year ago

hiroaki415 commented 1 year ago

Now I tried very simple structural analysis which has only single element. The element has uniformed load and point load which is apllied at the center of the element. Setting of loads shows below. load

And after analysis, I got these plots of force distribution by opsvis module. shear moment

I think these plots look like a little bit strange. The element has uniformed load, therefore the distribusion of shear force must have some gradient. Moreover the moment at node 2 must be nearly equal of zero because node 2 is free to rotate.

The values of element's force gotten by openseespy.opensees.eleFroce(1) shows below. [0.0, 0.9999999999999999, -5.551115123125783e-17, 0.0, 1.0, -2.7755575615628914e-17] The opsvis module seem to display values different from result of analysis.

Actually I think the analysis itself goes well because the values gotten by openseespy module are correct with theory. Is this any bug in opsvis module? Or did I get it a wrong way to use this module?

test2D.zip

sewkokot commented 1 year ago

Thank you for spotting this nasty bug. I fixed it now (figures attached below), but it needs more testing and in a few days I will release a new version of opsvis. test_2d_V test_2d_M

hiroaki415 commented 1 year ago

sewkokot

Thank you for early treatment. I hope the new stable version would come soon.

sewkokot commented 1 year ago

@hiroaki415 Please update the package and let me know if it works.

hiroaki415 commented 1 year ago

@sewkokot Finally, I could confirm that the new module work perfectly. Thank you very much!

By the way, I'm sorry but I have new question which is not related to this thread directly. The opsvis module can display value of element's cenrtral moment, but I could not find the way to get it by openseespy module. (For example, "0.375" is displayed in the plot showing above.) Maybe, is it calculated in opsvis module?

sewkokot commented 1 year ago

@hiroaki415 Thank you for your feedback! Indeed, the extreme values of section forces between element ends are calculated in opsvis.

hiroaki415 commented 1 year ago

@sewkokot Thank you for answering me. I solved a little problem cofusing me these days.