PARAMETER gapwidth = .5 PARAMETER leadwidth = 2 PARAMETER depth = 2 PARAMETER top = 1 PARAMETER theta=45 #define gamma_effective ((cos(theta*pi/180))) // virtual tension of face on plane // implicit surface tension is one boundary 1 parameters 1 //right edge of non-wetting gap x1: gapwidth/2 x2: p1 x3: 0 constraint 11 nonnegative //edge of right lead formula: x - gapwidth/2 boundary 2 parameters 1 //left edge of non-wetting gap x1: -gapwidth/2 x2: p1 x3: 0 constraint 12 nonnegative //edge of left lead formula: -x - gapwidth/2 constraint 10 //wetting leads formula: z = 0 energy: e1: gamma_effective*y e2: 0 e3:0 vertices: /* coordinates */ 1 -(gapwidth/2+leadwidth/2) -depth/2 0 constraint 10 12 2 -depth/2 boundary 2 3 -depth/2 boundary 1 4 (gapwidth/2+leadwidth/2) -depth/2 0 constraint 10 11 11 -(gapwidth/2+leadwidth/2) depth/2 0 constraint 10 12 12 depth/2 boundary 2 13 depth/2 boundary 1 14 (gapwidth/2+leadwidth/2) depth/2 0 constraint 10 11 21 -(gapwidth/2+leadwidth/2) -depth/2 top 22 (gapwidth/2+leadwidth/2) -depth/2 top 23 (gapwidth/2+leadwidth/2) depth/2 top 24 -(gapwidth/2+leadwidth/2) depth/2 top edges: 1 1 2 constraint 10 12 2 2 3 3 3 4 constraint 10 11 11 12 11 constraint 10 12 12 13 12 13 14 13 constraint 10 11 21 21 22 22 22 23 23 23 24 24 24 21 31 11 1 constraint 10 12 33 2 12 boundary 2 34 13 3 boundary 1 36 4 14 constraint 10 11 41 1 21 42 4 22 43 14 23 44 11 24 faces: 1 21 22 23 24 color blue //top 2 1 2 3 42 -21 -41 color red //front top 3 13 12 11 44 -23 -43 color cyan //back top 10 36 43 -22 -42 //rhs 11 31 41 -24 -44 //lhs 12 33 -12 34 -2 //bottom bodies: 1 1 2 3 10 11 12 volume top*(gapwidth + leadwidth)*depth read VV:= {inc:=10;while inc> 0 do { V;inc:=inc-1}}