Python API Scripts

2»

Comments

  • harryvanlangenharryvanlangen
    Hi, I wrote a script that warps a reference plane to a surface and keeps off-planar nodes perpendicular to this surface - a pictures is worth more than a thousand words:



    This way I can create a simple mesh using the move/copy mesh tool and warp this in the desired shape.

  • harryvanlangenharryvanlangen
    PS: the reference plane should be at z=0
  • disladisla
    SUUPER!!!. Great contribution. Thank you @harryvanlangen

    It sames to manage Shells and Solids too. !! :o It is like a super "project into surface".
    I can't wait to try on cicular platforms/handrails.


  • harryvanlangenharryvanlangen
    I made this fully parametric now so that I can also model surfaces that cannot be expressed as z=f(x,y), e.g. a cylinder. I attach the scripts for a cylinder, hypar and paraboloid as examples. I also attach the base grid again, but now properly put at z=0.0 so it will work out of the box.






  • disladisla
    edited April 2023
    Thanks again harryvanlangen!!
    Yesterday I was thinking that this tool could be phenomenal to fillet corners . ¿Do you think that could be possible? I would say the final shape should be the cilinder after selecting a set of nodes in the corner. Probably the normal at that point can be very hard to solve as there are too planes involved.
  • VictorVictor
    That's fantastic @harryvanlangen. I notice it seems to preserve the equal sizes of the elements, so it it kind of "wrapping" the mesh over the surface without stretching it more than it has to?
  • harryvanlangenharryvanlangen
    @disla - yes filleting would be among the applications, but will require a dedicated script with a parametric description of the fillet surface.

    @Victor - how the reference plane gets stretched all depends on its parametric description. The off-plane nodes remain perpendicular to the surface and at their original distance to it.
  • harryvanlangenharryvanlangen
    ... and a script for a spiral



  • harryvanlangenharryvanlangen
    Two more scripts for mapping a cube to a sphere. The base cube needs to be origin-centered.



  • VictorVictor
    Converts shell elements to line elements along their edges.


  • GermanGerman
    This scripts does a semiautomatic adaptive mesh refinement for models with all quad8 elements.
    Each refinement has to be solved manually and then run the scrips that estimates the error on each element from Zienkiewicz-Zhu Von Mises stress and makes a selection to be refined with elements exceeding a tolerance.
    Then Menu Mesh Tools > Refine > x3 should be used, and the refined mesh may be smoothed and begin the process again.

    In the picture you can see 6 refinements of a Kirsch problem.
    Blue elements are the elements selected to refine.



    Unfortunately in this example convergence is slow due to bad shape of the elements on the right side.
    In the next picture there are not problem of element bad shapes.



    Preocess could be further automated if there were an api command to do x3 refinement and if there were somehow to know using command mw.solve() if solver ends OK or not.
  • harryvanlangenharryvanlangen
    Mecway-Paraview workflow to display principal stress directions and magnitude, with Python script embedded in the Paraview State File:

    https://hvlanalysis.blogspot.com/2023/07/principal-stress-plots-with-mecway-and.html

    Example:


  • harryvanlangenharryvanlangen
    edited July 2023
    Mecway-Paraview workflow to display reinforcement requirements for concrete, with Python script embedded in the Paraview State File:

    https://hvlanalysis.blogspot.com/2023/07/reinforcement-of-concrete-structures.html

    Example:


  • VictorVictor
    Measure the length of line2 elements.

  • GermanGerman
    edited December 2023
    I have a named selection of elements to be refined and I want to refine them x3 in a script, but mw.refine_x3 requires the node_ids list.
    How can I get the list of node_id of the nodes belonging to the elements in a named elements selection or to a list of elements?

    I found this way, but perhaps there is an easier way:

    nodosref=[]
    for e in mw.named_selection("Refine"):
    for n in mw.nodes(e):
    if n not in nodosref:
    nodosref.append(n)

    mw.refine_x3(nodosref)
  • VictorVictor
    @German it's OK to have duplicates for refine_x3() so you could omit the if n not in nodosref: line. Or if you need no duplicates for something else, use a dictionary with node numbers as keys like: 
    nodosref={}
for e in mw.named_selection("Refine"):
  for n in mw.nodes(e):
    nodosref[n] = None
mw.refine_x3(nodosref.keys())
  • GermanGerman
    Thank you @Victor .
    Is it possible to run a script from another script?

    Regards.
  • VictorVictor
    @German You can use the Python import functionality. I would think this sort of helper function makes more sense as a user-defined thing like that to keep the core API simpler.
  • GermanGerman
    Is there any way to use arrays (like matrices) in a script?
  • VictorVictor
    You can't use Numpy arrays but you can use a list of lists. Eg. to create a 4x3 "array" initialized to all zeros:
    a = [[0] * 4 for i in range(3)]
    and index it like
    a[0][0] = 3
  • GermanGerman
    edited December 2023
    I improved the adaptive mesh refinement script.
    Now, before refinex3 is applied, it improves the list of nodes to be refined, adding those vertex that have 2 neighbor vertex (from different elements) to be refined. This avoids the generation of high valence vertex (and highly distorted elements) during refinement.
    Compare the old refinement sequence.


    with the new sequence


    Note the nodes that are added to be refined that didn't belong to the original elements to be refined.

  • SebastianmaklarySebastianmaklary
    @German That's really cool! (Semi-) Automatic mesh refinement is a great feature. Will you upload the improved script as well? :smile:
  • GermanGerman
    Im attaching 2 scrips.
    One works with CCX solver and uses ZZ error estimation.
    The other works with internal solver ald usas VM stress difference at nodes as error estimation.
  • cwharpecwharpe
    edited March 8
    FLIPLINE2: A Line2 coordinate ranking script to adjust U-axis directions. A time-saving alternative to Invert Tool requiring less selection scrutiny.



    Any constructed mesh might result in random alignments of element axis (through construction order, mirroring, extrusion direction, etc). If (U,V,W) output vectors are of interest, it is best to align the element axis.

    Notably, Line2 element U-axis directions are not changeable thru Element Properties (U-axis locked out). Fortunately, the Invert (Mesh) Tool will reverse selected elements. However, as beam arrays grow larger so does the time to selectively pick out the random "backward" elements.

    FlipLine2 aims to save time by aligning blanket selections of elements at once, regardless of initial orientation, thereby eliminating the selective pick. For every Line2 element selected, this script geometrically ranks the end-node coordinates w.r.t. a chosen Global Axis, and re-assigns their order to reflect a positive U-axis direction.

    RECAP: Invert Tool will reverse EVERY element selected. FlipLine2 will ONLY reverse the "backward" elements (per chosen direction), leaving the other elements in the selection untouched.

    Test out the script on the included file below. Pre-select ALL the elements.
Sign In or Register to comment.

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Sign In Register