I was having issues when splitting the strip, when i would change is position, for central position, it works fine.
reset # Clear Prep
reset aprepro # Clear old variables
# Create the ground plane(model for antenna)
create surface rectangle width 100 height 250 zplane
create vertex -39.52 12.16 0 on surface 1
create vertex 39.52 12.16 0 on surface 1
create vertex 30.4 88.16 0 on surface 1
create vertex -30.4 88.16 0 on surface 1
create vertex 12.16 101.84 0 on surface 1
create vertex -12.16 101.84 0 on surface 1
create vertex 41.04 -48.64 0 on surface 1
create vertex -41.04 -48.64 0 on surface 1
create vertex 15.2 -104.88 0 on surface 1
create vertex -15.2 -104.88 0 on surface 1
create vertex -30.4 -89.68 0 on surface 1
create vertex 30.4 -89.68 0 on surface 1
create vertex -38 -71.44 0 on surface 1
create vertex 38 -71.44 0 on surface 1
create vertex 0 104.88 0 on surface 1
create vertex -41.04 -7.6 0 on surface 1
create vertex 41.04 -7.6 0 on surface 1
create curve vertex 5,8
create curve vertex 8,10
create curve vertex 10,19
create curve vertex 19,9
create curve vertex 9,7
create curve vertex 7,6
create curve vertex 6,21
create curve vertex 21,11
create curve vertex 11,18
create curve vertex 18,16
create curve vertex 16,13
create curve vertex 13,14
create curve vertex 14,15
create curve vertex 15,17
create curve vertex 17,12
create curve vertex 12,20
create curve vertex 20,5
create surface curve 5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21 on surface 1
intersect surface 1,2
view top
#{freq = 0.9}
#{lam = 0.3/freq*1000} # Units in mm
#{a = 0.005*lam}
#{w = 4*a}
#{L = 0.23*lam} # Quarter-wave monopole
# Create the monopole
create surface rectangle width {w} height {L} yplane
#{v_monopole = Id("volume")}
move volume {v_monopole} z {L/2}
rotate surface 3 angle 90 about z include_merged
move surface 3 location 0 -19.76 {L/2} include_merged
# Split the monopole surface near the ground plane to create a voltage source curve
split surface in volume {v_monopole} across location position 0 {-w} {w} location position 0 {w} {w}
#{s_monopole_neg = Id("surface")}
# Imprint and merge to create common curves among surfaces at intersections
imprint all
merge all
# Assign surfaces with PEC material and create the voltage source
nsem assign surface all material "PEC"
nsem voltage source "port" pos surface {s_monopole_neg-1} neg surface {s_monopole_neg} impedance 50
block 1 surface all
#{mesh_size = lam/10}
surface all scheme pave
surface {s_monopole_neg} size {w/2}
mesh surface {s_monopole_neg}
surface {s_monopole_neg-1} size {w}
mesh surface {s_monopole_neg-1}
surface 2 size auto factor 5
mesh surface 2
surface not is_meshed size {mesh_size}
mesh surface not is_meshed
nsem print status
#save cub5 "monopole.cub5" overwrite journal
Also when trying multiples strips, same problem arises and mesh of ground plane is not done in correct way, i think is related to splitting for voltage source function.
reset # Clear Prep
reset aprepro # Clear old variables
# Create the ground plane(model for antenna)
create surface rectangle width 100 height 250 zplane
create vertex -39.52 12.16 0 on surface 1
create vertex 39.52 12.16 0 on surface 1
create vertex 30.4 88.16 0 on surface 1
create vertex -30.4 88.16 0 on surface 1
create vertex 12.16 101.84 0 on surface 1
create vertex -12.16 101.84 0 on surface 1
create vertex 41.04 -48.64 0 on surface 1
create vertex -41.04 -48.64 0 on surface 1
create vertex 15.2 -104.88 0 on surface 1
create vertex -15.2 -104.88 0 on surface 1
create vertex -30.4 -89.68 0 on surface 1
create vertex 30.4 -89.68 0 on surface 1
create vertex -38 -71.44 0 on surface 1
create vertex 38 -71.44 0 on surface 1
create vertex 0 104.88 0 on surface 1
create vertex -41.04 -7.6 0 on surface 1
create vertex 41.04 -7.6 0 on surface 1
create curve vertex 5,8
create curve vertex 8,10
create curve vertex 10,19
create curve vertex 19,9
create curve vertex 9,7
create curve vertex 7,6
create curve vertex 6,21
create curve vertex 21,11
create curve vertex 11,18
create curve vertex 18,16
create curve vertex 16,13
create curve vertex 13,14
create curve vertex 14,15
create curve vertex 15,17
create curve vertex 17,12
create curve vertex 12,20
create curve vertex 20,5
create surface curve 5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21 on surface 1
intersect surface 1,2
view top
#variables at 900 MHz
#{freq_1 = 0.9}
#{lam_1 = 0.3/freq_1*1000} # Units in mm
#{a_1 = 0.005*lam_1}
#{w_1 = 4*a_1}
#{L_1 = 0.23*lam_1} # Quarter-wave monopole
# Create the monopole at 900 MHz
create surface rectangle width {w_1} height {L_1} yplane
#{v_monopole_1 = Id("volume")}
#{s_monopole_1 = Id("surface")}
move volume {v_monopole_1} z {L_1/2}
rotate surface 3 angle 90 about z include_merged
move surface 3 location 0 -19.76 {L_1/2} include_merged
# Split the monopole at 900MHz surface near the ground plane to create a voltage source curve
split surface in volume {v_monopole_1} across location position 0 {-w_1} {w_1} location position 0 {w_1} {w_1}
#{s_monopole_neg_1 = Id("surface")}
#variables at 1900 MHz
#{freq_2 = 1.9}
#{lam_2 = 0.3/freq_2*1000} # Units in mm
#{a_2 = 0.005*lam_2}
#{w_2 = 4*a_2}
#{L_2 = 0.23*lam_2} # Quarter-wave monopole
# Create the monopole at 1900 MHz
create surface rectangle width {w_2} height {L_2} yplane
#{v_monopole_2 = Id("volume")}
#{s_monopole_2 = Id("surface")}
move volume {v_monopole_2} z {L_2/2}
rotate surface 4 angle 90 about z include_merged
move surface 4 location 0 34.96 {L_2/2} include_merged
# Split the monopole at 1900MHz surface near the ground plane to create a voltage source curve
split surface in volume {v_monopole_2} across location position 0 {-w_2} {w_2} location position 0 {w_2} {w_2}
#{s_monopole_neg_2 = Id("surface")}
# Imprint and merge to create common curves among surfaces at intersections
imprint all
merge all
# Assign surfaces with PEC material and create the voltage source
nsem assign surface all material "PEC"
nsem voltage source "port_1" pos surface {s_monopole_neg_1-1} neg surface {s_monopole_neg_1} impedance 50
nsem voltage source "port_2" pos surface {s_monopole_neg_2-1} neg surface {s_monopole_neg_2} impedance 50
block 1 surface all
#{mesh_size = lam_2/10}
surface all scheme pave
surface {s_monopole_neg_1} size {w_1/2}
mesh surface {s_monopole_neg_1}
surface {s_monopole_neg_1-1} size {w_1}
mesh surface {s_monopole_neg_1-1}
surface {s_monopole_neg_2} size {w_2/2}
mesh surface {s_monopole_neg_2}
surface {s_monopole_neg_2-1} size {w_2}
mesh surface {s_monopole_neg_2-1}
surface 2 size auto factor 7
mesh surface 2
surface not is_meshed size {mesh_size}
mesh surface not is_meshed
nsem print status
#save cub5 "monopole.cub5" overwrite journal
Thanks for attention