Advanced Ansys on Fracture and others

Dear alls,

 I am intersting in computational mechanics and now post some references which relate to ansys.

May be you know command KSCON, which is use for make concentration mesh at crack tip. I known how to use this command but very intersting its Ansys APDL code. So I am looking Ansys APDL code for KSCON command. If anybody know code fot that command, could you please send me?

Thank you.

Chau Nguyen

set lstbarID {100 200 300 400 500 600 700 800 900 1000}
set lstNod1 {1 3 15 6 7 10 9 13 2 16 16}
set lstNod2 {2 5 14 9 8 11 12 14 4 17 18}
set flag {0 0 0 0 0 0 0 0 0 0}

set lenbarID [llength $lstbarID]
for {set i 0} {$i <= $lenbarID} {incr i} {
 if {[lsearch $lstNod2 [lindex $lstNod1 $i]] != -1 || [lsearch $lstNod1 [lindex $lstNod2 $i]] != -1 || [lsearch [lreplace $lstNod1 $i $i] [lindex $lstNod1 $i]] != -1 || [lsearch [lreplace $lstNod2 $i $i] [lindex $lstNod2 $i]] != -1} {
      puts "3layer [lindex $lstNod1 $i]"
      puts "3layer [lindex $lstNod2 $i]"
 } else {
      puts "2layer [lindex $lstNod1 $i]"
      puts "2layer [lindex $lstNod2 $i]"
 }
}

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

set lstbarID {100 200 300 400 500 600 700 800 900 1000 1100 1200}
set lstNod1 {1 3 15 6 7 10 9 13 2 16 16 19}
set lstNod2 {2 5 14 9 8 11 12 14 4 17 18 20}
puts "================================"
set lst [list]

set lenbarID [llength $lstbarID]
for {set i 0} {$i < $lenbarID} {incr i} {
 if {[lsearch $lstNod2 [lindex $lstNod1 $i]] != -1 || [lsearch $lstNod1 [lindex $lstNod2 $i]] != -1 || [lsearch [lreplace $lstNod1 $i $i] [lindex $lstNod1 $i]] != -1 || [lsearch [lreplace $lstNod2 $i $i] [lindex $lstNod2 $i]] != -1} {
      set lst00 [lindex $lstbarID $i]
      if {[lsearch $lst $lst00] != -1} {
         continue
      } else {
         #puts "3layer $lst00"
         lappend lst $lst00
      }
     
      if {[lsearch $lstNod2 [lindex $lstNod1 $i]] != -1} {
          set lst21 [lindex $lstbarID [lsearch $lstNod2 [lindex $lstNod1 $i]]]
          if {[lsearch $lst $lst21] != -1} {
              continue
          } else {
              #puts "3layer $lst21"
              lappend lst $lst21
          }
      }
      if {[lsearch $lstNod1 [lindex $lstNod2 $i]] != -1} {
          set lst12 [lindex $lstbarID [lsearch $lstNod1 [lindex $lstNod2 $i]]]
          if {[lsearch $lst $lst12] != -1} {
              continue
          } else {
              #puts "3layer  $lst12"
              lappend lst $lst12
          }
      }
      if {[lsearch [lreplace $lstNod1 $i $i] [lindex $lstNod1 $i]] != -1} {
          set lst11 [lsearch [lreplace $lstNod1 $i $i] [lindex $lstNod1 $i]]
          if {$i <= $lst11} {
             set lst111 [lindex $lstbarID [expr $lst11 + 1]]
             if {[lsearch $lst $lst111] != -1} {
                 continue
             } else {
                 #puts "3layer $lst111"
                 lappend lst $lst111
             }
          } else {
             set lst112 [lindex $lstbarID $lst11]
             if {[lsearch $lst $lst112] != -1} {
                 continue
             } else {
                 #puts "3layer $lst112"
                 lappend lst $lst112
             }
          }
      }
      if {[lsearch [lreplace $lstNod2 $i $i] [lindex $lstNod2 $i]] != -1} {
          set lst22 [lsearch [lreplace $lstNod2 $i $i] [lindex $lstNod2 $i]]
          if {$i <= $lst22} {
             set lst221 [lindex $lstbarID [expr $lst22 + 1]]
             if {[lsearch $lst $lst221] != -1} {
                 continue
             } else {
                 #puts "3layer $lst221"
                 lappend lst $lst221
             }
          } else {
             set lst222 [lindex $lstbarID $lst22]
             if {[lsearch $lst $lst222] != -1} {
                 continue
             } else {
                 #puts "3layer [lindex $lstbarID $lst22]"
                 lappend lst $lst222
             }
          }
      }
      #puts $lst
 } else {
      #puts "2layer [lindex $lstNod1 $i]"
      #puts "2layer [lindex $lstNod2 $i]"
 }
}

set lenlst [llength $lst]
for {set i 0} {$i < $lenlst} {incr i} {
   puts "3Layers [lindex $lst $i][lindex $lst [expr $i + 1]]"
   puts [lindex $lst $i]
   puts [lindex $lst [expr $i + 1]]
   incr i
}