mp ( 'dens', 2, 0 ) # Simple cylinder for i in range ( 4 ): mapdl. r ( 2 ) # Aluminum properties (or something) mapdl. # Define higher-order SOLID186 # Define surface effect elements SURF154 to apply torque # as a tangential pressure mapdl.
Here is the corresponding PyMAPDL script: 0001 asel, r, loc, x, RADIUS local, 11, 1 csys, 0 aatt, 2, 2, 2, 11 amesh, all finish 1 / 386.1 mp, dens, 2, 0 !- ! Simple cylinder !- *do, ICOUNT, 1, 4 cylind, RADIUS, HEIGHTH_TIP, HEIGHT, 90 *( ICOUNT - 1 ), 90 * ICOUNT *enddo nummrg, kp lsel, s, loc, x, 0 lsel, r, loc, y, 0 lsel, r, loc, z, 0, HEIGHT - H_TIP lesize, all, ELEMSIZE * 2 mshape, 0 mshkey, 1 esize, ELEMSIZE allsel, all VSWEEP, ALL csys, 1 asel, s, loc, z, HEIGHT - H_TIP +0.0001, HEIGHT0. !- ! Define higher-order SOLID186 ! Define surface effect elements SURF154 ! which is used to apply torque ! as a tangential pressure !- /prep7 et, 1, 186 et, 2, 154 r, 1, r, 2, !- ! Aluminum properties (or something) !- mp, ex, 1, 10e6 mp, nuxy, 1.