//-----------------------------------------
#macro Erlenmeyer_Shape_1 (
Base_H, // base height
Base_Half_Width,// half base radius
Neck_Len, // neck lenght
Neck_R, // neck radius
Fillet_R, // = r1 < Base_H -2*r2
Base_Border_R,//= r2 + r1 < Base_H
Merge_On, //
) //---------------------------------
//-----------------------------------------
#local D = 0.0001 ;
//-----------------------------------------
#local R1 = Fillet_R;
#local X1 = (Neck_R+Fillet_R);
#local Y1 = Base_H;
#local M1 = < X1,Y1,0>;
// basis torus cross-section
#local R2 = Base_Border_R;
#local X2 = Base_Half_Width-Base_Border_R;
#local Y2 = Base_Border_R;
#local M2 = <X2,Y2,0>;
//-----------------------------------------
// angle between x-direction and (M1,M2) :
#if (X1 < X2)
#local Cone_Angle =
180-abs(atan((Y2-Y1)/abs(X2-X1)));
#else
#local Cone_Angle =
abs(atan((Y2-Y1)/(X2-X1)));
#end
//-----------------------------------------
// distance M1,M2 via Pythagoras:
#local M_Dist =
sqrt(pow(X2-X1,2)+pow(Y2-Y1,2));
#local M2_S =
sqrt(pow(M_Dist,2)-pow(R1+R2,2));
// Winkel bei M1 in Dreieck S_M1_M2:
#local In_Angle = abs(asin(M2_S/M_Dist)));
#local X_Angle = Cone_Angle-In_Angle ;
#local XSi = X1-(R1+R2)*cos(X_Angle);
#local YSi = Y1-(R1+R2)*sin(X_Angle);
#local Si =<XSi,YSi,0>;
// oberer Tangentenpunkt
#local T1 =
M1-<R1*cos(X_Angle),R1*sin(X_Angle),0>;
// unterer Tangentenpunkt
#local T2 =
M2+<R2*cos(X_Angle),R2*sin(X_Angle),0>;
// the body -------------------------------
#if ( Merge_On = 1 ) merge{
#else union{
#end
// neck
cylinder{<0,-D,0gt;,
<0,Neck_Len,0>,Neck_R
translate<0,M1.y,0> }
// fillet
difference{
cylinder{<0,T1.y-D,0>,
<0,M1.y,0>,T1.x}
torus{ X1,R1 translate<0,Y1,0>}
} // end of difference
// base cone
cone{<0,T2.y,0>,T2.x,
<0,T1.y,0>,T1.x}
// base round + center fill
cylinder{<0,-R2,0>,
<0,R2,0>, X2
translate<0, Y2,0>}
torus{ X2, R2 translate<0, M2.y,>}
} // end of union or merge
#end //----------------------- end of macro
//----------------------------------------- |
