Scalloping and Stress Concentration in DRIE-Manufactured Comb-Drives

: In the last decades, microelectromechanical systems have been increasing their number of degrees of freedom and their structural complexity. Hence, most recently designed MEMSs have required higher mobility than in the past and higher structural strength and stability. In some applications, device thickness increased up to the order of tens (or hundred) of microns, which nowadays can be easily obtained by means of DRIE Bosch process. Unfortunately, scalloping introduces stress concentration regions in some parts of the structure. Stress concentration is a dangerous source of strength loss for the whole structure and for comb-drives actuators which may suffer from side pull-in. This paper presents an analytical approach to characterize stress concentrations in DRIE micro-machined MEMS. The method is based on the linear elasticity equations, the de Saint-Venant Principle, and the boundary value problem for the case of a torsional state of the beam. The results obtained by means of this theoretical method are then compared with those obtained by using two other methods: one based on ﬁnite difference discretization of the equations, and one based on ﬁnite element analysis (FEA). Finally, the new theoretical approach yields results which are in accordance with the known value of the stress concentration factor for asymptotically null radius notches.

Since 2013 [16,17], a new class of MEMS, equipped with Conjugate Surface Flexure Hinges (CSFH) has been developed [18,19] and fabricated [20]. A sample of CSFH is depicted in Figure 1. These systems consist essentially of micro compliant mechanisms where the flexure hinges are manufactured as CSFHs. In such kind of flexure, the thickest rigid elements are linked by a certain number of CSFHs, where each CSFH consists of a curved beam together with a portion of a conjugate-profile pair, as illustrated in Figure 1. According to the CSFH design, the center of each arc is coincident with the center of the elastic weights of the curved beam [21]. The availability of CSFH hinges gave rise to the design and fabrication of new devices, such as microrobots [16,17,[22][23][24] micromechanisms [17,25], microgrippers [18,[26][27][28], and even a microtribometer [29].
In-plane MEMSs perform properly only if their flexible elements offer high compliance selectivity around the perpendicular axis. Hence, high aspect ratio structures are desirable. While surface micro-machining offers, typically, thickness from 2 to 5 µm, the Deep-Reactive Ion Etching (DRIE) process is used to obtain thickness of 50 µm and much more.
Since 1993 [30], Bosch DRIE process has introduced an effective method to provide, simultaneously, high mask selectively and very high anisotropy of the etched structure. In 1999, the scalloping observed on vertical walls during time multiplexed deep etching (TMDE), the roughness of horizontal surfaces exposed to the glow discharge, and the radius at the bottom of etched features were highlighted [31]. Few years later, it was demonstrated that important features of the samples, such as achievable anisotropy, etching uniformity, fillet radii, and, mainly, surface roughness, are strongly dependent on chamber pressure, applied coil and electrode power, and reactant gases flow rate [32]. This experimental observation lays the groundwork for the present investigation. In fact, scalloped sidewall roughness affects the performance of DRIE-manufactured comb-drives from both an electrical and a mechanical point of view. An example of DRIE-manufactured rotary comb-drive is depicted in Figure 1.
Considering the electrical aspect, scalloping determines, along the thickness of the device layer, gap variations between the facing surfaces of the fixed and movable fingers. Therefore, it affects the force or the torque exerted by the linear or rotary comb-drive, respectively. Furthermore, charge concentrations occur in sharp edges [33], determining variations of the electric field along the thickness direction.
From a mechanical perspective, it was observed that [31,32] post etch behavior for specimens with high surface roughness always indicated low fracture strength, whereas for specimens with better surface quality there was a wider distribution in sample strength. In addition, recent experimental tests [34] showed that some teeth are subject to lateral pull-in, with possible large deflection and instability.
The present investigation focuses on the mechanical aspect. As discussed in the next paragraph, one source of strength decay is the stress concentration factor, which holds in those zones where surface discontinuities appear. In high aspect ratio DRIE Bosch process, problems and challenges regarding scalloping attenuation have been reported in various investigations [35][36][37][38][39][40], although the problem presented several difficulties.
More recently, it has been shown that the sidewall roughness varies with the depth and depends on the trench width. In addition, the surface upper region exhibits a scalloping morphology, while the rougher lower region shows a curtaining morphology [41]. Although cryogenic-DRIE (see for example Refs. [42][43][44][45]) has recently been employed to fabricate trenches with highly vertical sidewalls and obtained relatively smooth surfaces, Bosch DRIE remains of great interest for the construction of many devices.
In 2012, a study [46] dedicated to Through-Silicon Via (TSV) technology showed that an etch method not based on Bosch process can get quite smooth sidewalls with no scalloping, together with a good control of the etched profile.
Since scalloping is not acceptable in specific applications, a dedicated effort has been provided to control and reduce it, mainly based on the optimization of DRIE process parameters [47][48][49][50], or, alternatively, removing it after fabrication [38].
Unfortunately, indentations and notches on the surface of a DRIE manufactured device give rise to high stress concentration that leads the samples to break at an unexpectedly low load.
The present work intends to contribute in filling this gap and to investigate about the shear stress distribution in the neighborhood of each single scalloped discontinuity on the walls device.
Hence, to develop a theoretical model, a uniform cantilever beam with circular cross section has been considered. To model the scalloping effect, a circular region, or notch, has been subtracted from the beam cross section. The center of the notch lies on the cross section profile.
By applying the linear elasticity theory to the beam with notched cross section, the shear stress distribution can be obtained analytically for the case of pure torsion. The idea of studying the beam torsional state can be justified by considering that stress concentration is actually due to tangential tensions. In addition, the shearing force is quite likely eccentric with respect to the center of shear of the cross-section. This implies a twisting moment acting on the beam, which is dangerous because of the very low torsion rigidity of thin-walled open profiles. Furthermore, the equations are integrated also by using a numerical approach based on finite discretization of the section and finite differences. Finally, the results are compared to those obtained with finite element analysis (FEA).

Motivation of This Work
The Bosh DRIE process is performed by means of pulsed regimes in two steps: in the first one, plasma etches the wafer along mainly the vertical direction, while in the second phase a passivation layer (Polytetrafluoroethylene) is chemically deposited on the silicon surface thanks to the action of a source gas (usually Octafluorocyclobutane). The two steps are iterated in such a way that, after each surface passivation phase, an etching phase follows, during which ions attack selectively the passivation layer at the bottom of the trench but not along the side walls. This iterative process permits very deep vertical etching [30]. However, the etching phase affects also the lateral walls and so the side is not perfectly planar, but it is rather undulating, as pictorially represented in Figure 2. The first image, reported in Figure 2a, represents a simplified scheme of a portion of a SOI wafer, after four cycles (note that the stop layer, represented in gray color, is reached after one more cycle only). The central cantilever beam appears after the wet etching phase, as illustrated in Figure 2b,c, during which the oxide is partially removed.
As known, any geometric irregularity on the cross section of a beam is responsible for stress concentration, which decreases its actual strength. The capability of micro-electromechanical systems (MEMS) to move is a direct consequence of their embedded flexures. According to the adopted process, a SOI (Silicon on Insulator) is used to obtain a suspended so-called device layer that can move along in-plane directions. The geometry of the device is designed and then transferred to the device layer by means of a mask. In the device layer, there are some pseudo-rigid parts that achieve motion capability because of the presence of flexible micro-beams. It is therefore clear that the strength of the micro-beams affects performance and reliability of the whole microsystem. Given the large variety of applications, micro-beams in MEMS are used for many different purposes. For example, they can be used either to create the fingers (linear or curved) of a comb-drive actuators or as a structural element (once again linear of curved) of a flexure. Therefore, these micro-beams are generally subject to local or distributed external forces that induce shear tensions on the beam cross-section. These forces may act along either an out-of-plane or an in-plane direction. Out-of-plane tensions are represented in Figure 2b, while in-plane tensions are depicted in Figure 2c. In the case of comb-drive actuators, fingers are subject to the lateral pull-in effect, and therefore they are inflected around an axis that is orthogonal to the pane of motion (in-plane inflection). On the other hand, surface forces, such as electrostatic forces, may act perpendicularly to the plane of motion, giving rise to out-of-plane inflection.
Furthermore, since the gaps among fingers can be reduced to as small as one micron, the possible inflection is responsible for stitching, which makes the MEMS practically impossible to use any further. The mechanics of adhesion this context has been extensively discussed in a recent review [51], where some fundamental parameters for microscale adhesion in MEMS are discussed in detail. Figure 3 shows scalloping over nine fingers of a circular comb-drive that the authors use for MEMS actuation. Although comb-drives fingers do not present direct contacts, mechanical impacts, vibrations or electrostatic field forces can still play a certain action on the fingers. However, a much greater stress state is generated on the curved beam which sustains the moving parts. For example, in the case of CSFH hinges [34], the curved beam is bent around an axis which is orthogonal to the plane of motion.

Theoretical Background
We refer to a prismatic beam as a region bounded by an exterior cylindrical surface and two terminal plane sections perpendicular to the generators [52]. The beam axis has a mild curvature compared to its thickness, so the theory for straight beams can be extended to this case of study. According to Ref. [53], this assumption is true when the ratio of the curvature radius to thickness is greater than 5.
The generic cross-section is depicted in Figure 5. Axes x and y, not necessarily principal of inertia, intersect each other in the centroid S and form an orthogonal Cartesian frame with z lying along the direction of the beam axis. The present formulation applies with no restriction to cross-sections of arbitrary shape and also to hollow beams. However, the case of simply-connected domain is analyzed here.
In the case of two axes of symmetry, the section centroid coincides with the center of shear, and extension, bending , and torsion are uncoupled if forces are applied to S [54][55][56][57][58].
The cross-section is a domain D of area Ω, bounded by the exterior directrix s whose direction is uniquely defined asking the associated outward normal n and tangent t unit vectors to form a right-handed system with the unit vectorê z .
In the following sections, shear stresses is determined from the theory of linear elasticity. Once the stress and strain tensors are defined, the semi-inverse method of Saint-Venant [59][60][61] is applied by setting [52,54,55,57,58] σ x = σ y = τ xy ≡ 0 , the equilibrium equations, neglecting body forces, read ∂σ z ∂z Assuming the exterior surface to be free of stress corresponds to setting on ∂D, with dy = n x ds and dx = −n y ds.
The hypotheses of isotropy and linear elastic behavior of the material lead to the stress-strain relations: where u = [u, v, w] T is the small displacement field, and E and G are Young's and shear modulus, respectively. The stress resultants on the cross-section are the normal and shear forces T z , T x , T y and the twisting and bending couples M z , M x and M y :

Shear Stresses in Torsion
As mentioned in Section 1, we focus on the shear stress concentration due to a scallop-like geometry in a beam subject to pure torsion. To compute the shear stresses, the solution of Poisson's problems, with Neumann or Dirichlet boundary conditions, is required. Such problems imply the use of Laplace operator in two dimensions [52,[54][55][56][57][58]62].
Pure torsion is defined as the partial solution of the problem of shear stresses where σ z ≡ 0 (implying the absence of shear forces and bending moments) [52]. The shear stresses τ = [τ xz , τ yz ] T satisfy the simplified equilibrium equation: When a twisting couple M z is applied to the free end of a clamped beam, the flatness of the section is preserved only if the beam has circular cross-sections. In all other cases, the longitudinal displacement w is not linear in the variables x and y, causing a warping of the section.
Assuming the torsion angle of the cross-section ω o = ω o (z) to be small, the displacement field u in pure torsion is devoid of rigid contributions, where the constant parameter θ = dω o dz is the twist of the centroidal fiber and g(x, y) is the so-called warping function.
Generally, the function g(x, y) is not linear in x and y and causes the loss of flatness of the section. It is defined only up to an additive constant and can be made unique by the constraint D g dxdy = 0.
Deriving shear strains as shown in Equations (6) and (7) allows us to write shear stresses as: By substituting Equation (11) into the equilibrium condition, Equation (9), leads to Together with the boundary condition in Equation (4) it determines the following Neumann's problem: with ∇ 2 Laplace operator in the plane of the domain D. The solution g(x, y) of this boundary value problem is unique up to an additive constant, and can be found in a closed form only for some simple cases (e.g., circular cross-sections, where g = 0 is an exact solution and fulfills the boundary condition [52]). Once the problem is solved, shear stresses are determined by means of Equation (11). Another formulation of the problem can be derived defining Prandtl's stress function Φ(x, y) such that The equilibrium condition, Equation (9), is satisfied when the stress function Φ is a solution of the problem: Since the addition of a constant to Φ does not affect stresses, the constant value Φ 0 on the outer boundary ∂D 0 can be set equal to 0. In hollow beams (multiple connected domains), the values of Φ 0i on the inner boundaries ∂D i cannot be chosen arbitrarily, but must satisfy the circuital conditions where Ω i denotes the area of the ith cavity. This constraint, which is necessary to ensure the displacement w has one value, cannot in general be solved explicitly within numerical analysis [56,63].

Circular Notched Cross-Section
As recalled in previous sections, the shear stress field is the solution of elliptic boundary value problems defined on the cross-section of the beam. Such problems are generally solved by means of finite differences or FEA.
However, the scheme presented here takes advantage of the simple form that the stress functions present for some particular geometries.
The warping function g(x, y) and its harmonic conjugate Z(x, y) = Φ(x, y) + 1 2 x 2 + y 2 are two functions which are coupled by the Cauchy-Riemann relations, i.e., ∂g ∂x Simple analytical solutions of Laplace equation are given by polynomials in the complex variable x + iy, such that the generic m-th term is In the case of negative values of m in Equation (15), the function g + iZ will present poles, i.e., singularity points where the function approaches an infinite value. Thus, negative exponents m are possible in Equation (15), provided that the generated poles fall out of the elastic region D. For example, for m = −1, there is one simple pole in the origin: In this case, the stress function is regular inside the domain D, vanishes on the boundaries x 2 + y 2 = r 2 and (x − h) 2 + y 2 = h 2 and satisfies the problem in Equation (14) (it is Prandtl's stress function for the torsion of a circular cross-section beam, of radius h, with a circular notch of radius r centred on its surface, see Figure 6). Considering the details shown in Figure 4, the circular shape of notches is assumed to be a first approximation fitting the real profile. Figure 6. Reference frame, cross-section with radius h, and notch with radius r.
The associated warping function is, up to a constant, Thus, for a beam with such a cross-section undergoing pure torsion, shear stresses can be easily evaluated with Equation (11) or (13) by partial differentiation of either g or Φ.

Theoretical Aspects in the Finite Element Formulation
The second order partial differential equation problem (with Dirichlet or Neumann boundary conditions) can be written, in its strong formulation, as: in the previous equations, in the case ϕ(x, y) = g, while The weak form of the problem, which can be implemented by most finite element solvers, can be obtained by a weighted residuals method, introducing a set of weighting functions, η. The integral sum over the domain D in the product of the weights η by the residual R = −∇ 2 ϕ + f 0 (x, y) is required to vanish Integration by parts, and Neumann's boundary condition, provide the weak form of the boundary value problem: with ∇ the gradient in the plane of the domain D.
The condition for η is η = 0 on ∂D. The integral over the boundary will then vanish in the case of Dirichlet boundary conditions. Equation (20) can be solved approximately using the finite element method; the domain D is first decomposed into a number of sub-domains, each one having n nodes.
The nodal approximations for the coordinates vector x = [x, y] T , the unknown function ϕ and the weighting function η are built using the same shape function N I . Once their values in the nodes are known, the interpolations which derive are: where the~(tilde) sign denotes the value of the indicated function.
Inserting the nodal approximations in Equation (21) into the weak form of Equation (20), and summing over the subdomains leads to the finite element problem: with the global matrix K IK and the vector P I formed by the elementary matrices and vectors: To compute shear stresses, we must solve the linear system in Equation (22). To this end, the cross-section D properties, the second order moments of inertia, I xx , I xy , I yy , and the centroid coordinates x S and y S , as well as the derivatives of the shape functions N with respect to x and y, must be known.

Stress Analysis
The presented methods can be specialized for the study of the above-mentioned class of CSFH-based microdevices. A member of this class is characterized by the presence of an embedded CSFH hinge, which is composed of a flexible curved element connecting a couple of conjugate profiles. One possible configuration is depicted in Figure 7. Relative motion among the pseudo-rigid parts i and j is possible because the thin cross-section of the curved element k behaves as a flexible rod. The technological limitations or the application complexity do not guarantee that the flexible curved beam k is subject only to pure bending. On the contrary, the tension field is rather complex and also includes shear component on the cross section. In fact, if a curved beam is considered, in-plane shear stress (as in Figure 2b) is always present, while out-of-plane, or transverse, shear stress ( Figure 2c) is likely to be present for several reasons. One of them, already pointed out in Section 1, is that the shearing force on the beam is quite likely eccentric with respect to the shear center, thus implying a twisting couple. Considering that beams are thin-walled, we know this implies remarkable shear stress.
Since the aim of this paper is restricted to the analysis of stress concentration, it is reasonable to simplify the general case of curved beam by a straight axis beam with in-plane shear at one end. For this reason, the stress distribution over a cross section of the beam is investigated, and increase in the overall tension maximum value is monitored to avoid unpleasant damages or fractures at unexpectedly low loads.
In the following paragraphs, the elementary beam theory is applied to the case of a circular cross-section shaft with a circular shaped notch undergoing pure torsion. Stresses are computed by using the analytical solution of the equations shown in Section 5. Then, two more methods, both having numerical nature, are used to compare their results with the analytical solutions: finite differences and finite elements.

Analytical Solution
For the geometry in Figure 6, it has been shown that Prandtl's stress function and the warping function have to satisfy Equations (11) and (13): The stress function and the partial derivatives can be rewritten in terms of polar coordinates, by introducing ρ 2 = x 2 + y 2 : Some physical parameters have to be imposed to evaluate the vector τ = [τ xz , τ yz ] T : the radii of the beam and of the notch h and r, as well as the twist of the centroidal fiber θ and the shear modulus G. Some plausible values are considered, although they do not resemble the exact geometry of the real device. This does not affect the accuracy of the results, since this calculations are not intended to give an exact solution, for which more accurate experimental studies are needed, in the first place. However, it is thought that, by assigning proper values to the parameters, the following computations give a close enough picture of the real stress field at a section of the hinge.
From the literature [32], we set 1 µm as fixed value for r, radius of the notch. Stresses were then studied for h linearly varying from 1 µm to 20 µm with steps of 1 µm.
We then imposed: The results were implemented and plotted using Matlab © . The stresses were initially computed deriving partially Φ in a 10,000 nodes sample grid in x and y over the domain D; the grid was then refined with 40,000 and, finally, 160,000 nodes. Figure 8a,b represents the shear vector fields in the cross section and around the notch for the two cases h = 10r and h = 20r, respectively, in different scaled grids. Figure 9a-c shows the shear stress field obtained by composing τ xz and τ yz , for different values of h r ratio, 20, 10 and 1, respectively. However, for the sake of the present investigation, the area around the notch is more interesting than the bulk of the cross-section, and so a detailed computation was performed within the same coordinate system and number of nodes over a sample grid narrowed to a 3r × 3r square. Results, for h = 20r, are shown in Figure 10a. Similar conclusions can be obtained by considering the case for which h = r, as depicted in Figure 10b. Stress is concentrated around the circular arc midpoint, being the point where the maximum stress is attained.

Numerical Solution via Finite Difference Discretization
The computation by means of finite differences has been implemented using the free version of the software FlexPDE (Version 6.36).
The resulting stress field, obtained over a 587 nodes grid, is shown in Figure 11a-c, together with the plots of its two components τ xz and τ yz , and the stress function Φ.

Numerical Solution via Finite Element Analysis
Finally, straight axis beams with circular section and circular notch were modeled by means of standard 3D modeling and FEA package ANSYS © .
Three different geometries were considered, thus three different solid models were built. Taking r = 1 µm, three cross sections with h = 20r, h = 10r and h = r, respectively, were used to extrude the solid model of the beam. Using the same physical parameters in the previous methods, and imposing again 5 • relative rotation among the end sections, the values of tension modules were obtained at the nodes. The adopted meshes, for the three h to r ratios, are represented in Figures 12a, 13a and 14a. The maximum shear stress fields are also reported in Figure 12b,c, for the case h = r, while the two groups in Figures 13b,c and 14b,c concern the two cases h = 10r and h = 20r, respectively.

Discussion about the Achieved Results
The obtained results have shown that the maximum shear stress holds in the center of the notch arc. Such point is called core point in Table 1. On the other hand, two other points have been considered, namely, the the upper and lower edge points of the notch, which correspond to the scallops spikes. In such points, called edge points in Table 1, the stress (denoted with τ edge ) is virtually null. This is consistent with the fact that such point belongs to a stagnation zone of the tension field. Of course, symmetry implies identical solution for the top and the bottom edge points. Shear stress has been computed in these points using three methods. The analytical solution yields null values (i.e., zero, to the best of the numerical approximation of the floating point accuracy), while the other two methods give values much smaller (practically null values) than the those calculated in the core. This difference with the analytical method is certainly due to the discretization approximation of the numerical approaches.
The proposed analytical method gave us the opportunity to calculate stress in many different points and sections. For example, τ edge and τ max could be calculated as a function of the ratio ρ of the notch radius r to the circular section radius h. The results obtained have shown that the maximum stress decreases as ρ = r h increases. This is consistent with the physics of the phenomenon since, as h reduces in size, the discontinuity in the field becomes less severe. Considering the edge points, τ edge also shows a similar decreasing trend as h decreases. As mentioned above, its magnitude remains significantly lower than the τ max stress achieved at the notch core. Figure 15 shows the dependency of τ max and τ edge upon the r h ratio. The diagram reported in Figure 15 concerns the tension absolute values that hold in the notch core and edge as a function of ρ. However, this result would convey rather little insight if it was not reduced to a physical model equivalent to the scalloped notch. For this purpose, there are still two particular problems to be solved: • the definition of an equivalence criterion for relating the adopted simplified geometry (see Figure 6) to the original geometry (see Figures 2 and 4); and • the introduction of a reliable index of the severity of stress concentration.

Adopted Equivalence Criterion
One advantage of the adopted cross section is that it is very simply described by means of analytical and geometrical tools. For example, referring back again to Figure 6, it is clear that the scalloped notch is represented by the circular segment having r, d and a as radius, depth and chord, respectively. One characteristic of scalloping is the dimensionless ratio δ of the notch depth d to the chord a. The two dimensionless parameters ρ and δ can be easily related using two simple, well-known geometrical equations for the circular segment, namely,

Stress Concentration Factor
The stress concentration factor K t , for the considered beam, can be calculated as the ratio of the maximum stress on the notch core τ max by the stress calculated on a point positioned on the other side of the beam section (with reference to Figure 6, this points is on the x axis, opposite to the origin O).
The adopted equivalence criterion and the introduction of K t give rise to an interesting interpretation of the iterative application of the analytical method. Figure 16 reports the values of the stress concentration factor K t as a function of the parameter η = ρ −1 = h r , which spans, in the real domain, from 1 2 (which correspond to a beam with null cross sectional area) to ∞. There are no reasons for η to go beyond the order of magnitude of tens and so 40 has been regarded as the maximum. The trend shows a sudden increase of K t for small increments of η. For values of η greater than about 5, K t tends asymptotically to 2. However, according to the equivalence criterion, the scalloping geometry is better characterized by the ratio δ = d a . By using Equations (25), it is possible to evaluate the ratio δ of the scalloping dept d to chord a as a function of η (see Figure 17). Finally, the relation between η and δ can be used to obtain the explicit dependency of K t upon δ, as reported in Figure 18. In fact, the dependency of K t on η and of the latter on δ, can be combined together to get, via function composition, K t (δ). When δ = 0, which corresponds to a smooth wall, no stress concentration occurs (K t = 1). However, for positive values of δ, the trend depicted in Figure 18 shows that the notch depth d has a great impact on stress concentration. The maximum value of K t is 2 and it corresponds to a notch which has its center on a vertical wall (with h >> r).
Thanks to the above-presented analysis, it is possible to state that the stress concentration factor tends asymptotically to the theoretical value 2 as the scalloping depth increases with respect to the scalloped chord, where the limit value is already known in the literature (see, for example, Ref. [65]).
Consequently, it can be also stated, as a general principle, that, for a DRIE micromachined beam with apparent and not treated scalloping, the strength can decrease to the half value of the strength of the same beam with no notches.

Conclusions
The wide range of applications has made MEMS more complicated. Therefore, their mechanical structure and internal stress state also become more complex. In many of those applications, MEMS are called to interact with the environment and so they need high aspect ratios, which can be achieved by manufacturing them via DRIE process. Unfortunately, the etched profiles cannot always get rid of the scalloping problem, which induces a dangerous increase of concentrated stress in the structure. Since this stress concentration gives rise to evident loss of structural strength for silicon, it is important to take into account the effects of scalloping on stress concentration in this kind of structures.
In this paper, three different methods are proposed for stress concentration characterization of silicon DRIE micromachined structures. All of them use a straight axis beam with a circular section and a circular notch. The first approach is based on the theory of linear elasticity and evaluates the shear stresses for pure torsion (stress concentration being due to shear state only). The other two methods use numerical approach: finite difference and finite element analysis. Although the three methods are very different in nature, the obtained values are fairly comparable. Furthermore, using the analytical method makes it easier to obtain diagrams that show the stress concentration factor with respect to the notch radius. The asymptotic value of the stress concentration factor, as a function of the notch relative radius, is in accordance with the literature. This leads to the conclusion that, in DRIE manufactured notched beams, strength is virtually half of the one with no notches.

Conflicts of Interest:
The authors declare no conflict of interest.