Get a Connectivity List for all the Elements on a Model - How to?

Hi,

It's increadible that you are comfortable with matlab but cannot do it by matlab-a very powerful software.

 

The tcl code is like this

*createmark elems 1 'by config' quad4

set ids [hm_getvalue elems mark=1 dataname=id]

set nodes [hm_getvalue elems mark=1 dataname=nodes]

set elems {}

set conns {}

foreach id $ids {n1 n2 n3 n4} [join $nodes] {

     set i 0

     foreach edge [list '$n1 $n2' '$n2 $n3' '$n3 $n4' '$n4 $n1' '$n2 $n1' '$n3 $n2' '$n4 $n3' '$n1 $n4']  conn {1 2 3 4 1 2 3 4} {

          if {[dict exists $elems $edge]]} {

               set e_c [dict get $elems $edge]

               eval dict lappend conns $id $conn [set ids_ [dict keys $e_c]]

               foreach id_ $ids_ {

                    dict lappend conns $id_ [dict get $e_c $id_] $id

               }

          }

          if {[incr i]<5} {dict lappend elems $edge $id $conn}

     }

}

set buf ''

foreach {id c_ids} $conns {

    foreach conn {1 2 3 4} {

         if {[dict exists $c_ids $con]} {

               append buf $id,[dict get $c_ids $conn]

          } else {

               append buf $id,0

          }

     }

     append buf \n

}

puts $buf

Hi @tinh,

 

Thank you very much for your answer!

 

The model I have is big qs (almost 400 000 elements) and the way I programmed it right now (probably not the best one) I'm running out of RAM before I finish the calculations..

 

I'll try your code and get back to you!

 

Thank you so much!

 

PS - How can I learn how to write TCL scripts?

<?xml version="1.0" encoding="UTF-8"?>

 

I get this error when I run it.

 

Is there something I should be doing differently?

 

I copied your code, pasted it to a text file that I saved with a .tcl extension and then I ran it from hypermesh

Please change:

         if {[dict exists $elems $edge]]} {

to:

         if {[dict exists $elems $edge]} {

 

Thanks @Q.Nguyen-Dai 

@goncalop please change 'puts $buf' by:

set fpt [open conns.txt w]

puts $fpt $buf

close $fpt

 

 

Because output 400k lines will hang hypermesh

Hello @Q.Nguyen-Dai, @tinh

 

I just did what you guys said and now the code reads like this:

 

*createmark elems 1 'by config' quad4

set ids [hm_getvalue elems mark=1 dataname=id]

set nodes [hm_getvalue elems mark=1 dataname=nodes]

set elems {}

set conns {}

foreach id $ids {n1 n2 n3 n4} [join $nodes] {

     set i 0

     foreach edge

1.  conn {1 2 3 4 1 2 3 4} {

          if {[dict exists $elems $edge]} {

               set e_c [dict get $elems $edge]

               eval dict lappend conns $id $conn [set ids_ [dict keys $e_c]]

               foreach id_ $ids_ {

                    dict lappend conns $id_ [dict get $e_c $id_] $id

               }

          }

          if {[incr i]<5} {dict lappend elems $edge $id $conn}

     }

}

set buf ''

foreach {id c_ids} $conns {

    foreach conn {1 2 3 4} {

         if {[dict exists $c_ids $con]} {

               append buf $id,[dict get $c_ids $conn]

          } else {

               append buf $id,0

          }

     }

     append buf \n

}

set fpt [open conns.txt w]

puts $fpt $buf

close $fpt

 

After copying this to a text file, changing the extension etc and running from hypermesh, I get this:

 

<?xml version="1.0" encoding="UTF-8"?>

 

Any idea?

 

Thank you very much!

 

 

 

 if {[dict exists $c_ids $con]} {

 

change to:

 

 if {[dict exists $c_ids $conn]} {

@Pandurang,

 

Thank you! I guess it worked now since It did not give me any error message.

 

Where can i find the output file?

@Pandurang,

 

Thank you! I guess it worked now since It did not give me any error message.

 

Where can i find the output file?

 

current working directory....

 

 

current working directory....

 

I searched the folder where I have the .tcl file and the folder where I have my model:

 

<?xml version="1.0" encoding="UTF-8"?>

 

here ^

 

But I can't find any conns.txt file..

 

Is it possible that Hypermesh is still writing the file?

 

Thank you!

 

Go to View -> Click on 'Command Window'

 

You will see one command line below panels in Hm...

Just type -> pwd 

hit enter...

you will get current working directory..

please find the file there...

 

Regards,

PD

 

@Pandurang,

 

Yeah, I would never find it... Sorry..

 

But I just found the file and it's perfect!

 

Thank you so much guys!

 

@tinh  @Q.Nguyen-Dai  

 

<?xml version="1.0" encoding="UTF-8"?>