# [blockMesh] Include list of points

 Register Blogs Members List Search Today's Posts Mark Forums Read

June 20, 2011, 09:03
Include list of points
#1
New Member

Join Date: Jun 2011
Posts: 3
Rep Power: 14
Hello,

I'm trying to mesh a volume with curved edges. I know I have to use the Spline function to do this.

My problem is that I have first to calculate the coordonates of the points of the curves, and I don't know how to do this on blockMesh.

Here is the file I use for Gambit, but I don't know how to do the same on blockMesh : how should I declare \$Y for example?

Quote:
 / Journal File for GAMBIT 2.4.6, Database 2.4.4, ntx86 SP2007051421 / Identifier "default_id4412" / File opened for write Tue Mar 29 16:53:30 2011. \$c = 2.68384 \$R=357.741 \$Q=44.0621 declare \$Y [1:25] \$Y [1] = 0.772 \$Y [2] = 0.932 \$Y [3] = 1.169 \$Y [4] = 1.574 \$Y [5] = 2.177 \$Y [6] = 2.647 \$Y [7] = 3.040 \$Y [8] = 3.666 \$Y [9] = 4.143 \$Y [10] = 4.503 \$Y [11] = 4.760 \$Y [12] = 4.924 \$Y [13] = 4.996 \$Y [14] = 4.963 \$Y [15] = 4.812 \$Y [16] = 4.530 \$Y [17] = 4.146 \$Y [18] = 3.682 \$Y [19] = 3.156 \$Y [20] = 2.584 \$Y [21] = 1.987 \$Y [22] = 1.385 \$Y [23] = 0.810 \$Y [24] = 0.306 \$Y [25] = 0 declare \$X [1:25] \$X [1] = 0.5 \$X [2] = 0.75 \$X [3] = 1.25 \$X [4] = 2.5 \$X [5] = 5 \$X [6] = 7.5 \$X [7] = 10 \$X [8] = 15 \$X [9] = 20 \$X [10] = 25 \$X [11] = 30 \$X [12] = 35 \$X [13] = 40 \$X [14] = 45 \$X [15] = 50 \$X [16] = 55 \$X [17] = 60 \$X [18] = 65 \$X [19] = 70 \$X [20] = 75 \$X [21] = 80 \$X [22] = 85 \$X [23] = 90 \$X [24] = 95 \$X [25] = 100 declare \$YF [1:25] \$YF [1]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[1]/50)^2)-\$R*COS(\$Q/2) \$YF [2]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[2]/50)^2)-\$R*COS(\$Q/2) \$YF [3]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[3]/50)^2)-\$R*COS(\$Q/2) \$YF [4]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[4]/50)^2)-\$R*COS(\$Q/2) \$YF [5]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[5]/50)^2)-\$R*COS(\$Q/2) \$YF [6]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[6]/50)^2)-\$R*COS(\$Q/2) \$YF [7]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[7]/50)^2)-\$R*COS(\$Q/2) \$YF [8]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[8]/50)^2)-\$R*COS(\$Q/2) \$YF [9]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[9]/50)^2)-\$R*COS(\$Q/2) \$YF [10]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[10]/50)^2)-\$R*COS(\$Q/2) \$YF [11]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[11]/50)^2)-\$R*COS(\$Q/2) \$YF [12]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[12]/50)^2)-\$R*COS(\$Q/2) \$YF [13]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[13]/50)^2)-\$R*COS(\$Q/2) \$YF [14]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[14]/50)^2)-\$R*COS(\$Q/2) \$YF [15]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[15]/50)^2)-\$R*COS(\$Q/2) \$YF [16]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[16]/50)^2)-\$R*COS(\$Q/2) \$YF [17]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[17]/50)^2)-\$R*COS(\$Q/2) \$YF [18]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[18]/50)^2)-\$R*COS(\$Q/2) \$YF [19]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[19]/50)^2)-\$R*COS(\$Q/2) \$YF [20]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[20]/50)^2)-\$R*COS(\$Q/2) \$YF [21]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[21]/50)^2)-\$R*COS(\$Q/2) \$YF [22]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[22]/50)^2)-\$R*COS(\$Q/2) \$YF [23]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[23]/50)^2)-\$R*COS(\$Q/2) \$YF [24]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[24]/50)^2)-\$R*COS(\$Q/2) \$YF [25]=SQRT(\$R^2-(\$R*SIN(\$Q/2)-\$R*SIN(\$Q/2)*\$X[25]/50)^2)-\$R*COS(\$Q/2) vertex create "p000" coordinates (0.0) ( 0.0) vertex create "p001" coordinates (0.5*\$c) ( \$YF [1]+\$Y [1]*\$c) vertex create "p101" coordinates (0.5*\$c) (\$YF [1]-\$Y [1]*\$c) vertex create "p002" coordinates (0.75*\$c) (\$YF [2]+\$Y [2]*\$c) vertex create "p102" coordinates (0.75*\$c) (\$YF [2]-\$Y [2]*\$c) vertex create "p003" coordinates (1.25*\$c) (\$YF [3]+\$Y [3]*\$c) vertex create "p103" coordinates (1.25*\$c) (\$YF [3]-\$Y [3]*\$c) vertex create "p004" coordinates (2.5*\$c) (\$YF [4]+\$Y [4]*\$c) vertex create "p104" coordinates (2.5*\$c) (\$YF [4]-\$Y [4]*\$c) vertex create "p005" coordinates (5*\$c) (\$YF [5]+\$Y [5]*\$c) vertex create "p105" coordinates (5*\$c) (\$YF [5]-\$Y [5]*\$c) vertex create "p006" coordinates (7.5*\$c) ( \$YF [6]+\$Y [6]*\$c) vertex create "p106" coordinates (7.5*\$c) (\$YF [6]-\$Y [6]*\$c) vertex create "p007" coordinates (10*\$c) ( \$YF [7]+\$Y [7]*\$c) vertex create "p107" coordinates (10*\$c) (\$YF [7]-\$Y [7]*\$c) vertex create "p008" coordinates (15*\$c) ( \$YF [8]+\$Y [8]*\$c) vertex create "p108" coordinates (15*\$c) (\$YF [8]-\$Y [8]*\$c) vertex create "p009" coordinates (20*\$c) ( \$YF [9]+\$Y [9]*\$c) vertex create "p109" coordinates (20*\$c) (\$YF [9]-\$Y [9]*\$c) vertex create "p010" coordinates (25*\$c) ( \$YF [10]+\$Y [10]*\$c) vertex create "p110" coordinates (25*\$c) (\$YF [10]-\$Y [10]*\$c) vertex create "p011" coordinates (30*\$c) ( \$YF [11]+\$Y [11]*\$c) vertex create "p111" coordinates (30*\$c) (\$YF [11]-\$Y [11]*\$c) vertex create "p012" coordinates (35*\$c) ( \$YF [12]+\$Y [12]*\$c) vertex create "p112" coordinates (35*\$c) (\$YF [12]-\$Y [12]*\$c) vertex create "p013" coordinates (40*\$c) ( \$YF [13]+\$Y [13]*\$c) vertex create "p113" coordinates (40*\$c) (\$YF [13]-\$Y [13]*\$c) vertex create "p014" coordinates (45*\$c) ( \$YF [14]+\$Y [14]*\$c) vertex create "p114" coordinates (45*\$c) (\$YF [14]-\$Y [14]*\$c) vertex create "p015" coordinates (50*\$c) ( \$YF [15]+\$Y [15]*\$c) vertex create "p115" coordinates (50*\$c) (\$YF [15]-\$Y [15]*\$c) vertex create "p016" coordinates (55*\$c) ( \$YF [16]+\$Y [16]*\$c) vertex create "p116" coordinates (55*\$c) (\$YF [16]-\$Y [16]*\$c) vertex create "p017" coordinates (60*\$c) ( \$YF [17]+\$Y [17]*\$c) vertex create "p117" coordinates (60*\$c) (\$YF [17]-\$Y [17]*\$c) vertex create "p018" coordinates (65*\$c) ( \$YF [18]+\$Y [18]*\$c) vertex create "p118" coordinates (65*\$c) (\$YF [18]-\$Y [18]*\$c) vertex create "p019" coordinates (70*\$c) ( \$YF [19]+\$Y [19]*\$c) vertex create "p119" coordinates (70*\$c) (\$YF [19]-\$Y [19]*\$c) vertex create "p020" coordinates (75*\$c) ( \$YF [20]+\$Y [20]*\$c) vertex create "p120" coordinates (75*\$c) (\$YF [20]-\$Y [20]*\$c) vertex create "p021" coordinates (80*\$c) ( \$YF [21]+\$Y [21]*\$c) vertex create "p121" coordinates (80*\$c) (\$YF [21]-\$Y [21]*\$c) vertex create "p022" coordinates (85*\$c) ( \$YF [22]+\$Y [22]*\$c) vertex create "p122" coordinates (85*\$c) (\$YF [22]-\$Y [22]*\$c) vertex create "p023" coordinates (90*\$c) ( \$YF [23]+\$Y [23]*\$c) vertex create "p123" coordinates (90*\$c) (\$YF [23]-\$Y [23]*\$c) vertex create "p024" coordinates (95*\$c) ( \$YF [24]+\$Y [24]*\$c) vertex create "p124" coordinates (95*\$c) (\$YF [24]-\$Y [24]*\$c) vertex create "p025" coordinates (100*\$c) ( \$YF [25]+\$Y [25]*\$c) vertex create "p125" coordinates (100*\$c) (\$YF [25]-\$Y [25]*\$c) edge create nurbs "p001" "p002" "p003" "p004" "p005" "p006" "p007" "p008" \ "p009" "p010" "p011" "p012" "p013" "p014" "p015" "p016" "p017" "p018" \ "p019" "p020" "p021" "p022" "p023" "p024" "p025" interpolate edge create nurbs "p101" "p102" "p103" "p104" "p105" "p106" "p107" "p108" \ "p109" "p110" "p111" "p112" "p113" "p114" "p115" "p116" "p117" "p118" \ "p119" "p120" "p121" "p122" "p123" "p124" "p125" interpolate edge create nurbs "p001" "p000" "p101" interpolate face create wireframe "edge.1" "edge.2" "edge.3" real face summarize "face.1"
I hope somebody can help me!

Thanks

Last edited by Hikachu; June 22, 2011 at 18:19.