search

ansys / pymapdl

Pythonic interface to MAPDL
https://mapdl.docs.pyansys.com
MIT License
424 stars 120 forks source link

Cannot plot rst file #198

Closed Shannon-John closed 4 years ago

Shannon-John commented 4 years ago

I have a rst file,but can't show it in pyansys. I don't known where the problem is. Thank you! Follwing are the APDLs: /TITLE,Truss Bridge Static Analysis /COM,Structural /PREP7 ET,1,BEAM4 ET,2,SHELL63 SECTYPE, 1, BEAM, I, duan, 0 SECOFFSET, CENT SECDATA,0.4,0.4,0.4,0.016,0.016,0.016,0,0,0,0 SECTYPE, 2, BEAM, I, XIANHENG, 0 SECOFFSET, CENT SECDATA,0.4,0.4,0.4,0.012,0.012,0.012,0,0,0,0 SECTYPE, 3, BEAM, I, FU, 0 SECOFFSET, CENT SECDATA,0.3,0.3,0.4,0.012,0.012,0.012,0,0,0,0 R,1,0.0187,0.17e-3,0.54e-3,0.4,0.4,0, R,2,0.0141,0.128e-3,0.415e-3,0.4,0.4,0, R,3,0.0117,0.541e-4,0.324e-3,0.3,0.4,0, R,4,0.3, , , , , , MP,EX,1,2.1E11 MP,PRXY,1,0.3 MP,DENS,1,7850 MP,EX,2,3.5E10 MP,PRXY,2,0.1667 MP,DENS,2,2500 N, ,0,0,-5, , , , NGEN,4,4,ALL, , ,12, , ,1, NGEN,2,1,ALL, , , , ,10,1, NGEN,2,1,2,10,4, ,16, ,1, NGEN,2,1,3,11,4, , ,-10,1, TYPE,1 MAT,1 REAL,1 ESYS,0 SECNUM,1 TSHAP,LINE E,11,14 E,12,13 TYPE,1 MAT,1 REAL,2 ESYS,0 SECNUM,2 TSHAP,LINE E,2,6 E,6,10 E,10,14 E,1,5 E,5,9 E,9,13 E,3,7 E,7,11 E,4,8 E,8,12 E,1,2 E,3,4 E,5,6 E,7,8 E,9,10 E,11,12 E,13,14 TYPE,1 MAT,1 REAL,3 ESYS,0 SECNUM,3 TSHAP,LINE E,3,6 E,6,11 E,4,5 E,5,12 E,2,3 E,1,4 E,6,7 E,5,8 E,10,11 E,9,12 TYPE,2 MAT,2 REAL,4 ESYS,0 TSHAP,QUAD E,1,2,6,5 E,5,6,10,9 E,9,10,14,13 NSYM,X,14,ALL ESYM, ,14,ALL NUMMRG,ALL, , , ,LOW NUMCMP,ALL FINISH /SOL NSEL,S, , ,23,24 D,ALL, ,0, , , ,UX,UY,UZ, , , , NSEL,S, , ,13,14 D,ALL, ,0, , , , ,UY,UZ, , , , NSEL,S, , ,1,2 F,ALL,FY,-10000 ALLSEL,ALL ACEL,0,10,0, ANTYPE,0 SOLVE FINISH /POST1 PLDISP,2 PLNSOL, U,SUM, 0,1.0 PLVECT,U, , , ,VECT,ELEM,ON,0 ETABLE,ZHOU_I,SMISC, 1 ETABLE,ZHOU_J,SMISC, 7 ETABLE,JIAN_I,SMISC, 2 ETABLE,JIAN_J,SMISC, 8 ETABLE,WAN_I,SMISC, 6 ETABLE,WAN_J,SMISC, 12 PRETAB,ZHOU_I,ZHOU_J,JIAN_I,JIAN_J,WAN_I,WAN_J PLLS,ZHOU_I,ZHOU_J,1,0 PLLS,JIAN_I,JIAN_J,1,0 PLLS,WAN_I,WAN_J,1,0 PRNSOL,U,COMP FINISH SAVE

akaszynski commented 4 years ago

Could you also please the python code you tried to use?

Shannon-John commented 4 years ago

Sorry! I list the python code below.

""" Script generated by pyansys version 0.39.14 """
import pyansys
ansys = pyansys.Mapdl(loglevel="INFO",override=True)
ansys.run("/TITLE,Truss Bridge Static Analysis")
ansys.run("/COM,Structural")
ansys.run("/PREP7")
ansys.et(1, "BEAM4")
ansys.et(2, "SHELL63")
ansys.sectype(1, "BEAM", "I", "duan", 0)
ansys.secoffset("CENT")
ansys.secdata(0.4, 0.4, 0.4, 0.016, 0.016, 0.016, 0, 0, 0, 0)
ansys.sectype(2, "BEAM", "I", "XIANHENG", 0)
ansys.secoffset("CENT")
ansys.secdata(0.4, 0.4, 0.4, 0.012, 0.012, 0.012, 0, 0, 0, 0)
ansys.sectype(3, "BEAM", "I", "FU", 0)
ansys.secoffset("CENT")
ansys.secdata(0.3, 0.3, 0.4, 0.012, 0.012, 0.012, 0, 0, 0, 0)
ansys.r(1, 0.0187, 0.17e-3, 0.54e-3, 0.4, 0.4, 0)
ansys.r(2, 0.0141, 0.128e-3, 0.415e-3, 0.4, 0.4, 0)
ansys.r(3, 0.0117, 0.541e-4, 0.324e-3, 0.3, 0.4, 0)
ansys.r(4, 0.3)
ansys.mp("EX", 1, 2.1E11)
ansys.mp("PRXY", 1, 0.3)
ansys.mp("DENS", 1, 7850)
ansys.mp("EX", 2, 3.5E10)
ansys.mp("PRXY", 2, 0.1667)
ansys.mp("DENS", 2, 2500)
ansys.n("", 0, 0, -5, "", "", "")
ansys.ngen(4, 4, "ALL", "", "", 12, "", "", 1)
ansys.ngen(2, 1, "ALL", "", "", "", "", 10, 1)
ansys.ngen(2, 1, 2, 10, 4, "", 16, "", 1)
ansys.ngen(2, 1, 3, 11, 4, "", "", -10, 1)
ansys.type(1)
ansys.mat(1)
ansys.real(1)
ansys.esys(0)
ansys.secnum(1)
ansys.tshap("LINE")
ansys.e(11, 14)
ansys.e(12, 13)
ansys.type(1)
ansys.mat(1)
ansys.real(2)
ansys.esys(0)
ansys.secnum(2)
ansys.tshap("LINE")
ansys.e(2, 6)
ansys.e(6, 10)
ansys.e(10, 14)
ansys.e(1, 5)
ansys.e(5, 9)
ansys.e(9, 13)
ansys.e(3, 7)
ansys.e(7, 11)
ansys.e(4, 8)
ansys.e(8, 12)
ansys.e(1, 2)
ansys.e(3, 4)
ansys.e(5, 6)
ansys.e(7, 8)
ansys.e(9, 10)
ansys.e(11, 12)
ansys.e(13, 14)
ansys.type(1)
ansys.mat(1)
ansys.real(3)
ansys.esys(0)
ansys.secnum(3)
ansys.tshap("LINE")
ansys.e(3, 6)
ansys.e(6, 11)
ansys.e(4, 5)
ansys.e(5, 12)
ansys.e(2, 3)
ansys.e(1, 4)
ansys.e(6, 7)
ansys.e(5, 8)
ansys.e(10, 11)
ansys.e(9, 12)
ansys.type(2)
ansys.mat(2)
ansys.real(4)
ansys.esys(0)
ansys.tshap("QUAD")
ansys.e(1, 2, 6, 5)
ansys.e(5, 6, 10, 9)
ansys.e(9, 10, 14, 13)
ansys.nsym("X", 14, "ALL")
ansys.esym("", 14, "ALL")
ansys.nummrg("ALL", "", "", "", "LOW")
ansys.numcmp("ALL")
ansys.finish()
ansys.run("/SOL")
ansys.nsel("S", "", "", 23, 24)
ansys.d("ALL", "", 0, "", "", "", "UX", "UY", "UZ", "", "")
ansys.nsel("S", "", "", 13, 14)
ansys.d("ALL", "", 0, "", "", "", "", "UY", "UZ", "", "")
ansys.nsel("S", "", "", 1, 2)
ansys.f("ALL", "FY", -10000)
ansys.allsel("ALL")
ansys.acel(0, 10, 0)
ansys.antype(0)
ansys.solve()
ansys.finish()
ansys.run("/POST1")
ansys.pldisp(2)
ansys.plnsol("U", "SUM", 0, 1.0)
ansys.plvect("U", "", "", "", "VECT", "ELEM", "ON", 0)
ansys.etable("ZHOU_I", "SMISC", 1)
ansys.etable("ZHOU_J", "SMISC", 7)
ansys.etable("JIAN_I", "SMISC", 2)
ansys.etable("JIAN_J", "SMISC", 8)
ansys.etable("WAN_I", "SMISC", 6)
ansys.etable("WAN_J", "SMISC", 12)
ansys.pretab("ZHOU_I", "ZHOU_J", "JIAN_I", "JIAN_J", "WAN_I", "WAN_J")
ansys.plls("ZHOU_I", "ZHOU_J", 1, 0)
ansys.plls("JIAN_I", "JIAN_J", 1, 0)
ansys.plls("WAN_I", "WAN_J", 1, 0)
ansys.prnsol("U", "COMP")
ansys.finish()
ansys.eplot()()
ansys.exit()
akaszynski commented 4 years ago

Looks like there's two bugs that I've got to fix for this. Let me patch this and release a new version on pypi.

akaszynski commented 4 years ago

Fixed in pyvista==0.39.18