Skip to content

Commit f8a9e3d

Browse files
jginternationaljginternational
committed
Force creation of point elements
1 parent 8d5cabe commit f8a9e3d

4 files changed

Lines changed: 135 additions & 54 deletions

File tree

‎kratos.gid/kratos.tcl‎

Lines changed: 11 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -9,7 +9,6 @@ namespace eval ::Kratos {
99
variable must_write_calc_data
1010
variable must_exist_calc_data
1111

12-
variable tmp_init_mesh_time
1312
variable namespaces
1413

1514
variable mesh_criteria_forced
@@ -200,7 +199,7 @@ proc Kratos::LoadCommonScripts { } {
200199
uplevel #0 [list source [file join $kratos_private(Path) scripts Writing $filename.tcl]]
201200
}
202201
# Common scripts
203-
foreach filename {Utils Launch Applications spdAuxiliar Menus Deprecated Logs} {
202+
foreach filename {Utils Launch Applications spdAuxiliar Mesh Menus Deprecated Logs} {
204203
uplevel #0 [list source [file join $kratos_private(Path) scripts $filename.tcl]]
205204
}
206205
# Common controllers
@@ -427,19 +426,10 @@ proc Kratos::Event_BeforeMeshGeneration {elementsize} {
427426
set tmp_init_mesh_time $inittime
428427
Kratos::Log "Mesh BeforeMeshGeneration start"
429428

430-
GiD_MeshData mesh_criteria to_be_meshed 1 lines [GiD_Geometry list line]
431-
GiD_MeshData mesh_criteria to_be_meshed 1 surfaces [GiD_Geometry list surface]
432-
GiD_MeshData mesh_criteria to_be_meshed 1 volumes [GiD_Geometry list volume ]
433-
434-
# We need to mesh every line and surface assigned to a group that appears in the tree
435-
foreach group [spdAux::GetAppliedGroups] {
436-
GiD_MeshData mesh_criteria to_be_meshed 2 lines [GiD_EntitiesGroups get $group lines]
437-
GiD_MeshData mesh_criteria to_be_meshed 2 surfaces [GiD_EntitiesGroups get $group surfaces]
438-
GiD_MeshData mesh_criteria to_be_meshed 2 volumes [GiD_EntitiesGroups get $group volumes]
439-
}
429+
Mesh::PrepareMeshGeneration $elementsize
440430

441431
# Change the mesh settings depending on the element requirements
442-
if {[Kratos::CheckMeshCriteria $elementsize]<0} {
432+
if {[Mesh::CheckMeshCriteria $elementsize]<0} {
443433
return "-cancel-"
444434
}
445435

@@ -459,9 +449,17 @@ proc Kratos::Event_MeshProgress { total_percent partial_percents_0 partial_perce
459449
proc Kratos::Event_AfterMeshGeneration {fail} {
460450
variable tmp_init_mesh_time
461451

452+
if {$fail} {
453+
Kratos::Log "Mesh generation failed"
454+
return
455+
}
456+
462457
# Change the mesh settings depending on the element requirements. Reset previous settings
463458
# catch {Kratos::ResetMeshCriteria $fail}
464459

460+
Mesh::AddPointElementsIfNeeded
461+
462+
465463
# Maybe the current application needs to do some extra job
466464
apps::ExecuteOnCurrentApp AfterMeshGeneration $fail
467465
set endtime [clock seconds]

‎kratos.gid/scripts/Mesh.tcl‎

Lines changed: 101 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,101 @@
1+
##################################################################################
2+
3+
namespace eval ::Mesh {
4+
Kratos::AddNamespace [namespace current]
5+
6+
7+
}
8+
proc Mesh::PrepareMeshGeneration {elementsize} {
9+
10+
GiD_MeshData mesh_criteria to_be_meshed 1 lines [GiD_Geometry list line]
11+
GiD_MeshData mesh_criteria to_be_meshed 1 surfaces [GiD_Geometry list surface]
12+
GiD_MeshData mesh_criteria to_be_meshed 1 volumes [GiD_Geometry list volume ]
13+
14+
# We need to mesh every line and surface assigned to a group that appears in the tree
15+
foreach group [spdAux::GetAppliedGroups] {
16+
GiD_MeshData mesh_criteria to_be_meshed 2 lines [GiD_EntitiesGroups get $group lines]
17+
GiD_MeshData mesh_criteria to_be_meshed 2 surfaces [GiD_EntitiesGroups get $group surfaces]
18+
GiD_MeshData mesh_criteria to_be_meshed 2 volumes [GiD_EntitiesGroups get $group volumes]
19+
}
20+
}
21+
22+
23+
proc Mesh::CheckMeshCriteria { elementsize } {
24+
25+
set force_mesh_order [dict create]
26+
set elements_used [spdAux::GetUsedElements]
27+
set forced_mesh_order -1
28+
foreach element_id $elements_used {
29+
set element [Model::getElement $element_id]
30+
if {[$element hasAttribute "MeshOrder"]} {
31+
set element_forces [$element getAttribute "MeshOrder"]
32+
if {$element_forces eq "Quadratic"} {
33+
set element_forces 1
34+
} else {
35+
set element_forces 0
36+
}
37+
dict set force_mesh_order $element_id $element_forces
38+
if {$forced_mesh_order eq -1} {
39+
set forced_mesh_order $element_forces
40+
} else {
41+
if {$forced_mesh_order ne $element_forces} {
42+
# W "The element $element_id requires a different mesh order"
43+
W "Incompatible mesh orders in elements"
44+
return -1
45+
}
46+
}
47+
}
48+
}
49+
50+
if {$forced_mesh_order ne -1} {
51+
52+
set element [lindex [dict keys $force_mesh_order] 0]
53+
set previous_mesh_order [write::isquadratic]
54+
set current_mesh_type [Kratos::GetMeshOrderName $previous_mesh_order]
55+
set desired_mesh_type [Kratos::GetMeshOrderName $forced_mesh_order]
56+
if {$previous_mesh_order ne $forced_mesh_order} {
57+
W "The element $element requires a different mesh order: $desired_mesh_type"
58+
W "Currently the mesh order is $current_mesh_type. please change it in the menu Mesh > Quadratic type"
59+
return -1
60+
}
61+
}
62+
return 0
63+
}
64+
65+
proc Mesh::AddPointElementsIfNeeded {} {
66+
# Foreach groups assigned in tree
67+
set condition_groups [spdAux::GetUsedConditions]
68+
69+
# condition_groups is a dict of conditionid -> list of group nodes (tdom)
70+
foreach condid [dict keys $condition_groups] {
71+
set cond [Model::getCondition $condid]
72+
if {$cond eq ""} {
73+
continue
74+
}
75+
set topology [$cond getTopologyFeature "Point" 1]
76+
if {$topology eq ""} {
77+
continue
78+
}
79+
W " Groups assigned to condition $condid will be meshed with Point elements."
80+
81+
set group_nodes [dict get $condition_groups $condid]
82+
foreach node_tdom $group_nodes {
83+
set group_id [get_domnode_attribute $node_tdom n]
84+
# Get the goodname of the group
85+
set group_id [write::GetWriteGroupName $group_id]
86+
87+
set node_ids [GiD_EntitiesGroups get $group_id nodes]
88+
set new_nodeids [list]
89+
foreach nodeid $node_ids {
90+
set new_nodeid [GiD_Mesh create element append Point 1 [list $nodeid]]
91+
# Add to same groups as the node
92+
lappend new_nodeids $new_nodeid
93+
}
94+
GiD_EntitiesGroups assign $group_id elements $new_nodeids
95+
}
96+
}
97+
98+
# if group element is point and has topology for points
99+
100+
# foreach nodeid [GiD_Mesh] GiD_Mesh create element append Point 1 [list $nodeid]
101+
}

‎kratos.gid/scripts/Utils.tcl‎

Lines changed: 0 additions & 41 deletions
Original file line numberDiff line numberDiff line change
@@ -292,47 +292,6 @@ proc Kratos::GetMeshBasicData { } {
292292
return $result
293293
}
294294

295-
proc Kratos::CheckMeshCriteria { elementsize } {
296-
297-
set force_mesh_order [dict create]
298-
set elements_used [spdAux::GetUsedElements]
299-
set forced_mesh_order -1
300-
foreach element_id $elements_used {
301-
set element [Model::getElement $element_id]
302-
if {[$element hasAttribute "MeshOrder"]} {
303-
set element_forces [$element getAttribute "MeshOrder"]
304-
if {$element_forces eq "Quadratic"} {
305-
set element_forces 1
306-
} else {
307-
set element_forces 0
308-
}
309-
dict set force_mesh_order $element_id $element_forces
310-
if {$forced_mesh_order eq -1} {
311-
set forced_mesh_order $element_forces
312-
} else {
313-
if {$forced_mesh_order ne $element_forces} {
314-
# W "The element $element_id requires a different mesh order"
315-
W "Incompatible mesh orders in elements"
316-
return -1
317-
}
318-
}
319-
}
320-
}
321-
322-
if {$forced_mesh_order ne -1} {
323-
324-
set element [lindex [dict keys $force_mesh_order] 0]
325-
set previous_mesh_order [write::isquadratic]
326-
set current_mesh_type [Kratos::GetMeshOrderName $previous_mesh_order]
327-
set desired_mesh_type [Kratos::GetMeshOrderName $forced_mesh_order]
328-
if {$previous_mesh_order ne $forced_mesh_order} {
329-
W "The element $element requires a different mesh order: $desired_mesh_type"
330-
W "Currently the mesh order is $current_mesh_type. please change it in the menu Mesh > Quadratic type"
331-
return -1
332-
}
333-
}
334-
return 0
335-
}
336295

337296
proc Kratos::GetMeshOrderName {order} {
338297
switch $order {

‎kratos.gid/scripts/spdAuxiliar.tcl‎

Lines changed: 23 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -423,6 +423,29 @@ proc spdAux::GetUsedElements {{alt_un ""}} {
423423
return $lista
424424
}
425425

426+
# returns ["condition1" {gNode1 gNode2} "condition2" {gNode3 gNode4}]
427+
proc spdAux::GetUsedConditions {{root ""}} {
428+
set resultDict [dict create]
429+
430+
431+
set xp1 "./condition/group"
432+
if {$root eq "" } {
433+
set root [customlib::GetBaseRoot]
434+
set xp1 "//condition/group"
435+
}
436+
437+
foreach gNode [$root selectNodes $xp1] {
438+
set condition_node [$gNode parent]
439+
set condition_name [$condition_node @n]
440+
441+
set cond [Model::getCondition $condition_name]
442+
if {$cond eq ""} {continue}
443+
dict lappend resultDict $condition_name $gNode
444+
445+
}
446+
return $resultDict
447+
}
448+
426449
proc spdAux::LoadIntervalGroups { {root ""} } {
427450
customlib::UpdateDocument
428451
variable GroupsEdited

0 commit comments

Comments
 (0)