Difference between revisions of "R3DC Report 2-08"

From New IAC Wiki
Jump to navigation Jump to search
Line 383: Line 383:
 
=3-D Analysis=
 
=3-D Analysis=
  
The term "3-D Analysis" refers to an analysis based on a model that includes both endplates, a nose plate, a bottom plate, several carbon fiber rods, and a Hexcel backing plate.  The deflection calculation was performed using several permutations involving the wire tension, the presence of gravity, and the use of carbon fiber rods, the Hexcel backing plate as well as the impact of a ball joint in the drift  chamber nose plate.  Analysis begins with the simplest case in which only the impact of gravitational forces on the endplate is deflection is evaluated for an ionization chamber structure modeled without carbon rods, a Hexcel backing plate, nor a ball joint.  The analysis then consideres the impact of adding carbon fiber rods, the Hexcell sheet, and the ball joint proceeds to more complex cases (e.g. wire & gravity loading, carbon fiber rods and Hexcel backing included). However, to date, none of these analysis cases includes the gas pressure due to the inclusion of the mylar front sheet.  
+
The term "3-D Analysis" refers to an analysis based on a model that includes both endplates, a nose plate, a bottom plate, several carbon fiber rods, and a Hexcel backing plate.  The deflection calculation was performed using several permutations involving the wire tension, the presence of gravity, and the use of carbon fiber rods, the Hexcel backing plate as well as the impact of a ball joint in the drift  chamber nose plate.  Analysis begins with the simplest case in which only the impact of gravitational forces on the endplate is deflection is evaluated for an ionization chamber structure modeled without carbon rods, a Hexcel backing plate, nor a ball joint.  The analysis then considers the case in which wires are strung between the endplates and the impact of adding carbon fiber rods, the Hexcell sheet, and the ball joint on the endplate deflection. Figure 14 shows the 3-D model with labels indicating the location of the loadings included in the FEA program. However, to date, none of these analysis cases includes the gas pressure due to the inclusion of the mylar front sheet.  
  
Figure 14 shows the 3-D model as well as the loadings included in the FEA program. Below is a list describing in detail the applied loading and boundary conditions.
+
  Below is a list describing in detail the applied loading and boundary conditions.
  
 
*The Hinge 1 labeled A is fixed in all degrees of freedom.
 
*The Hinge 1 labeled A is fixed in all degrees of freedom.

Revision as of 16:30, 25 April 2008

Introduction

The Region 3 Drift Chamber (R3 DC) conceptual design for the Hall-B 12-GeV upgrade (CLAS12)is described below. One of the design objectives is to achieve nearly full acceptance, in the lab frame, for relativistic forward-boosted final-state particles while maintaining a momentum resolution of 0.5 to 1% for a 5-GeV charged particle. Six identical R3 DCs are presently located 4.7 meters from the target and together will have a [math]2\pi[/math] azimuthal coverage and presently span from 5 to 45 degrees in polar angle as measured from the direction of the electron beam. The central plane through the R3 DC will be tilted 25 degrees from a perpendicular to the beam line towards the target. The shadow of the toroidal coils will result in a minimum gap of 16 cm between each R3 DC and will hinder reconstructing the trajectory of scattered particles in that region. This space will house the readout boards, attendant electronics, cabling, etc. The active area for each R3 DC will then be [60 - ?] in azimuth,

The R3 DC design for CLAS12 is based on the successful hexagonal cell structure of the current CLAS R3 DC. The design criteria is delineated below:

  1. Smaller cell size of 2.07 cm to accommodate a luminosity of [math]1035 /(cm^2 s^1)[/math] and maintain an accidental occupancy of less the 4% for high tracking efficiency..
  2. Sense wires will be increased to 30 um from the current 20 um diameter. This means that the wire will require a higher tension to offset droop and electrostatic forces.
  3. All wires in each superlayer will lie in a plane instead of on a cylindrical surface. The staggered by one half cell or “brick wall pattern” of the hexagonal cells will be maintained for each successive layer.
  4. To enhance the azimuthal (φ) resolution, each superlayer will be +/- stereo, i.e. the pitch will be +6 degrees in the six layers of the first superlayer and the pitch will be -6 degrees for the second superlayer array of six wire planes. The current R3 DC has 0 and 6 degree pitch respectively for the two superlayers.
  5. There will be a total of 4928 wires of which 1344 (check this) will be sense wires.
  6. The endplates will be oriented at 60 degrees with respect to one another. Precisely holes will be drilled into each endplate with tolerance of better than 100 um for the placement of the drift chamber wires. Feed through assemblies are placed in these holes to precisely position the sense, field, and guard wires. The expected overall position accuracy will be 250 um per layer, with a wire location accuracy of 200 um.
  7. The deflection along the length of the endplates may not exceed 50 um.
  8. The R3 DCs will be supported from only two points at the outer backplane. This makes the R3 DCs self supporting.
  9. The nose plate for each R3 DC will be positioned onto a precision support hexagonal ring for aligning all R3 DCs with respect to the surveyed Hall B coordinate system.

The mechanical engineering aspects of the R3 DC is discussed below in light of the above design criteria. We present a two dimensional Finite Element Analysis. The three-dimensional FEA is currently under study

R3 Description

Geometry

An exploded view of a single R3 sector is shown below. The back plate, nose plate, and two endplates (see Figure 2.1) represent the key components to the main ionization chamber structure. The back plate is approximately 5 meters long and will likely have two mount points on the backplate to support and position the drift chambers. The nose plate, more than half a meter long, will hold a ball joint that mounts to a hexagonal ring and constrains the chamber location at a third point. The drift chamber wires are held in place by the end plates which have carbon fiber rods spanning the gap near the chamber entrance in order to reduce end plate deflection arising from the wire tension, the weight of the chamber, and the differential gas pressure between the chamber and the experimental hall. The front face of the chamber has a thin window to allow the passage of scattered particles with limited multiple scattering. The exit window of the chamber is constructed from a composite material known as Hexcell clad with a thin carbon fiber sheet to reduce end plate deflections near the back of the chamber. Figure 1.2 shows a side view of a single sector; six of these sectors are oriented around a beam pipe, as shown in Figure 2.3.


R3BreakOutDiagram.jpg
label Distance (mm)
[math]Bck_{Len1}[/math] 5004.92
[math]Bck_{Len2}[/math] 4840.04
[math]End_{Len1}[/math] 4830.17
[math]End_{Len2}[/math] 4665.1
[math]Nse_{Len1}[/math] 614
Figure 2.1: The general features of a R3 Drift Chamber. The chamber wires string between the end plates are not shown. Parameter Tables for the variables in the above pictures


R3 SideView.jpg
label value
[math]D_1[/math] 4765.6 mm
[math]\theta_1[/math] 5 degrees
[math]\theta_2[/math] 45 degrees
[math]\theta_3[/math] 25 degrees
Figure 2.2: A side view of a R3 drift chamber. The distance from the target center to the first sense wire is labeled [math]D_1[/math]. The min an max acceptance angles are labeled [math]\theta_1[/math] and [math]\theta_2[/math] respectively. The tilt angle with respect to a normal to the beam line is labeled [math]\theta_3[/math].


R3AllSectors.jpg R3 Gap.jpg
Figure 2.3: Beam line view showing the corresponding arrangement of the six comprising R3.


Mestayer's Geometry Excel Spreadsheet

Material Specifications

Physical properties germane to the materials used to construct the R3 dirft chambers are shown in Table 1. End plates for the R3 drift chambers in the previous 6 GeV CLAS detector were fabricated by gluing thin (0.048 inch) stainless steel sheets onto both sides of a 2 inch Polyeurethane foam. The current design expects to use the same structure for the end plates. The use of carbon fiber rods to minimize end plate deflections is another design feature borrowed from the previous 6 GeV CLAS detector design. Aluminum is included below as a cadidate material for capping the edges of the end plates.


Table 1. Material Properties of R3 Components.

Material Youngs Modulus Density
[math]Pa (N/m^2)[/math] [math]\frac{kg}{m^3}[/math]
Polyeurethane Foam FR-3700 [math]1.4 \times 10^8[/math] 240
Stainless Steel [math]2.1 \times 10^{11}[/math] 8030
Carbon Fiber Rods [math]1.5 \times 10^{11}[/math] 580
Aluminum [math] 7\times 10^{10}[/math] 2700

End Plates and Nose Plate Drawings

Preliminary drawings of the end plates (Figure 2.4) and nose plate (Figure 2.5) are shown below. The end plate and nose plate are two critical parts of the chamber. The end plate will provide the support to constrain the detector wires in space whereas the nose plate will determine the tilt angle of the detector and provide a mounting connection.

EndPlateDrawing 2-08.jpg
Figure 2.4: A preliminary End plate design for R3 Drift Chamber.


NosePlateDrawing 2-08.jpg
Figure 2.5: A preliminary Nose plate design for R3 Drift Chamber.

Endplate Deflection

Wire tension

The drift chamber wires are oriented at 6 degrees with respect to the horizontal, as shown in the Figure 3.1. Wires which slope downward from left to right make an angle of [math]\theta_2 = 90-30-6 = 54[/math] degrees with respect to a perpendicular to the endplate and have a force component along the perpendicular of [math]T_{1} = T \sin(54)[/math]. Wires which slope upward from left to right make an angle with repect to the perpendicular to the endplate of [math]\theta_3 = 90 + 30 - 6 = 114[/math] degrees. If T represents the tension in the wire, then the component of this force perpendicular to the endplate is [math]T_{2}=T \sin(180-\theta_3) = T \sin(66)[/math]. The net force perpendicular to the endplate will be [math]T( \sin(66) + \sin(54))[/math]. There are 112 sense wires in each of the 12 wire layers ( 2 super layers with 6 wire layers each)[Mestayer07]. Each sense wire is 30 microns in diameter and tensioned using a 20.7 g weight. There are 28 Field wire layers 140 microns in diameter with each wire tensioned at 63 grams and 4 guard wire layers 140 micron in diameter with each wire tensioned at 185 grams. The total tension from all wires on the endplate will be


[math]112 \left ( 12 \cdot 0.0207 + 28 \cdot 0.063 + 4 \cdot 0.185 \right ) \times 9.8 = 3024 N = 680 lbs[/math].

Therefore,

[math]T_1 = \frac{3024}{2} \sin(54) =1223 N[/math]
[math]T_2 = \frac{3024}{2} \sin(66) =1381 N[/math]

and

[math]T_{perp} = T_1+T_2 = 2604 N[/math] = perpendicular force on an endplate
TensionForceComponents.jpg
label Value (Degrees)
[math]\theta[/math] 30
[math]\theta_1[/math] 6
[math]\theta_2[/math] 90 - [math]\theta[/math] - [math]\theta_1[/math]
[math]\theta_3[/math] 90 + [math]\theta[/math] - [math]\theta_1[/math]
Figure 3.1: Sketch identifying the direction of the sense wire tension in terms of the [math]\theta[/math] and [math]\theta_1[/math].



Point Load Deflection

A simple back-of-the-envelope calculation for deflection was performed assuming aluminum endplates and a point-load force (representing the distributed force from the wires into a single load). The force on an endplate due only to the tension from the 4,928 wires is about 5700 N /2. A simple calculation based on the definition of Young's modulus can yield an order-of-magnitude-level estimate for the deflection of an aluminum endplate in which the wires exert a force perpendicular to the endplate. The deflection [math]y[/math] of the endplate in terms of Young's Modulus [math](Y)[/math] is given by the equation:

[math]Y \equiv \frac{F}{4wy} \left( \frac{l}{t}\right )^3[/math]

where

[math]Y = 7 \times 10^{10}N/m^2 \equiv[/math] Young's Modulus for Aluminum
[math]F = 2850 N \equiv[/math]point force/load
[math]w =0.526 m \equiv[/math] width of the endplate
[math]y \equiv[/math] deflection of endplate due to point force
[math]l =4.02 m \equiv[/math] length of the endplate
[math]t = 0.05 m \equiv[/math] thickness of the endplate

By substituting in the above values and solvign the equation, a deflection of 7 mm is expected for a 5 cm thick aluminum endplate. If a 5 cm thick stainless steel [math](Y=2.1 \times 10^{11} N/m^2 )[/math] endplate were used, the deflection would decrease from 7 mm to 2.5 mm (due to the linear dependence of the deflection on the Young's Modulus). The previous work performed by Kevin Folkman at the IAC developed a similar conclusion.

Distributed Load Deflection

The endplate deflection calculations were refined to a more realistic case in which the load is distributed evenly along the endplate. The tension of the wires on the endplate was divided into three regions; W1,W2, and W3. As shown in Figure 3.3.1, the load W3 corresponds to a region along the endplate in which wires from both super layers are creating a load on the endplate. Loads W1 and W2 represent regions in which only a single super layer is causing a load. The endplate is "simply" supported in order to calculate the deflection. The term "simply supported" implies that one side (left side of an endplate for this case) may rotate but not translate (resulting in two reaction forces) whereas the other side (the right side of and endplate for this case) is supported such that it may rotate and translate (resulting in a single, normal reaction force). Figure 7 labels the loads as well as distances over which the loads are applied on the end plate. A total force of 340 lbs exerted by the drift chamber wires is decomposed into its components according to the angle made by the wire on the endplate.

EndPlateDeflectionHandCalc.jpg
Label value
a 99.02 mm
b 96.22 mm
c 4025.88 mm
d 335.17 mm
e 273.88 mm
[math]W_1[/math] 0.3=1223/(C+B) N/mm
[math]W_2[/math] 0.61 N/mm
[math]W_3[/math] 0.32 =1381/(C+D)N/mm
Figure 3.3.1: A diagram of the loads on a simply supported endplate. The table to the right gives the load values and distances.


The following calculations are in accordance to a standard approach covered in all undergraduate "Mechanics of Materials" engineering courses. This method is based on integrating the elastic curve equation for a beam. The equation of the elastic curve relates deflection (y) to Moments (M) at given linear locations along the beam (x):


[math]\frac{d^2y}{dx^2} = \frac{M(x)}{E I_A}[/math]


where [math]E[/math] represents Young's Modulus and [math]I_A[/math] is the Area Momentum of Inertia. The Area Moment of Inertis [math](I_A)[/math] is defined as the second moment of an area about a given axis:


[math]I_{A,x} = \int x^2 dA =[/math] Area moment of inertia about the x-axis.


The mass moment of Inertia is defined as:


[math]I= \int r^2 dM = \rho I_A[/math]


In order to perform these calcuations, the following assumption were made:

  1. The end plate is made of homogenous uniform elastic material with a Young's Modulus of [math]E=210 x 10^23 N/mm^2[/math] (Stainless Steel).
  2. The horizontal component of the loading is completely neglected.
  3. The deflection in the x-direction is negligible.
  4. The end plate is rectangular in shape with a width of 526 mm and thickness of 50.8 mm resulting area moment of inertia ([math]I_A[/math]) =5746.4 [math]mm^4[/math].
  5. The loading provided by wire tensions are equivalent to a uniformly distributed load.
  6. The maximum deflection occurs along the length which is represented by [math]c[/math] in the diagram.
  7. The loading applied is under the elastic limit of the material.
  8. The effects of the endplate weight are ignored.

The moments experienced along the endplate may be divided into 5 separate equations, each equation valid only for its corresponding region:


[math]M_a = -\frac{W_1 \cdot b + W_2 \cdot c + W_3 \cdot d}{2} \cdot x[/math]
[math]M_b = -\frac{W_1 \cdot b + W_2 \cdot c + W_3 \cdot d}{2} \cdot (a+x) + \frac{W_1 \cdot x^2}{2}[/math]
[math]M_c = -\frac{W_1 \cdot b + W_2 \cdot c + W_3 \cdot d}{2} \cdot (a+b+x) + W_1 \cdot b \left ( \frac{b}{2} +x \right) + \frac{W_2 \cdot x^2}{2}[/math]
[math]M_d = -\frac{W_1 \cdot b + W_2 \cdot c + W_3 \cdot d}{2} \cdot (a+b+c+ x) + W_1 \cdot b \left ( \frac{b}{2} +c+x \right) + W_2\cdot c \left ( \frac{c}{2} + x \right) + \frac{W_3 \cdot x^2}{2}[/math]
[math]M_e = -\frac{W_1 \cdot b + W_2 \cdot c + W_3 \cdot d}{2} \cdot (a+b+c+d+ x) + W_1 \cdot b \left ( \frac{b}{2} +c+d+x \right) + W_2\cdot c \left ( \frac{c}{2} + d+x \right) + W_3 \cdot d \left (\frac{d}{2} +x \right)[/math]

where [math]x[/math] represents the linear distance along the endplate. If we assume that the largest magnitude of deflection occurs in region c, then we may integrate the above 2nd order elastic curve equation once and determine the critical point which represents the point of maximum deflection ( contra flexure).

[math]\int \frac{dy}{dx}dy = \frac{dy}{dx} = \int M_c dx[/math]
= [math]-\frac{W_1 \cdot b + W_2 \cdot c + W_3 \cdot d}{2} \cdot (ax+bx+x^2/2) + W_1 \cdot b \left ( \frac{bx}{2} +x^2/2 \right) + \frac{W_2 \cdot x^3}{6} +C_{integration}[/math]

where [math]C_{integration} =[/math] a constant of integration. Integrating the equation a 2nd time will result in an equation for the amount of the deflection. The constants of integration my be found by applying boundary conditions at points where the deflection is zero (at the ends or [math]x=0[/math] and [math]x=L[/math]).

Setting [math]\frac{dy}{dx}[/math] to zero and solving for x gives the linear distance corresponding to the point of maximum deflection of x = 1654.96 mm. In other words, the maximum deflection occurs 1654.96 mm from region [math]b[/math] of the endplate or [math]1654.96 + b + a = 1850.2[/math] mm from the right hand side of the endplate. Substituting the value of x above results in a maximum deflection of 2.2 mm. Note, this deflection is only slightly smaller than that found with a point load assumption and a stainless steel endplate (2.2 mm compared to 2.5 mm, respectively).

Distributed Load FEA

A Finite Element Analysis (FEA) was performed on the same stainless steel endplate design using the commercially available computer software ANSYS. For this analysis, the loading shown previously in the "Distrubuted Load" section was used. The maximum deflection observed from the FEA was 2.14 mm, as shown in Figure 3.4.1; this value differs from the hand calculation by only 3%. In addition to validating the FEA analysis with hand calculations, the above results also reveal the challenge of reducing the endplate deflection using stiffer materials due to an essentially linear dependence of the deflection on the Young's Modulus. The cubic dependence of the deflection on the thickness [math](t)[/math] of the endplate, however, can be used to reach the minimum endplate deflection criteria of 50 [math]\mu m[/math], in addition to the us of stiffer material. Increasing the endplate thickness, however, comes at the cost of increasing the end plate mass. As a result, composite materials were considered in order to minimize weight. The performance of these composite materials is shown in the 3-D model section of this report.

EndPlateDeflection 5cmx53cmx4m ThickStainlessSteel.jpg
Figure 3.4.1: ANSYS FEA results for deflection of a stainless steel endplate design for a distrubed loading.

Composite End Plate

Effective Young's Modulus [math]Y_{eff}[/math]

The R3 endplates used in the CLAS 6 GeV detector were comprised of a polyeurethane foam core 50.8 mm thick (2 inches) with 1.2 mm thick (0.048 inch) stainless steel sheets epoxied to each side, as shown in Figure 3.5.1.

EndPlateMaterialDimensions.jpg
Figure 3.5.1: Polyeurethane foam/stainless steel composite endplate construction and dimensions.


An effective Young's Modulus ([math]Y_{eff}[/math]) was used in the FEA model in order to simplify the model. The effective Young's Modulus for this endplate, weighted by the volume of each plate, is

[math]Y_{eff} = \frac{V_p Y_p + V_s Y_s }{V_{tot}}[/math]

where

[math]Y_p = 1.4 \times 10^8 N/m^2[/math]
[math]Y_s = 2.1 \times 10^{11} N/m^2[/math]
[math]V_p = L \times W \times 2[/math] inches
[math]V_s = L \times W \times 0.096[/math] inches


Therefore, the effective Young's Modulus was found to be:

[math]Y_{eff} = \frac{2(1.4 \times 10^8) +0.096(2.1 \times 10^{11})}{2.096} = 9.75 \times 10^{9} N/m^2[/math]

FEA Deflection Results

The effective Young's Modulus of [math]9.75 \times 10^9 N/m^2[/math] was used to determine the deflection of the endplate due to a distributed load (see Figure 3.5.2). Results indicate a maximum transverse deflection of 32.4 mm which is larger than that observed with the stainless steel material. However, when gravity is included, the composite material will deflect less due to its lower density.

DtotNoRods.jpg
Configuration X max Y max Z max D max
(mm) (mm) (mm) (mm)
No rods .67 32.4 .04 32.4
Figure 3.5.2: FEA predicted deflections for the polyeurethane foam/stainless steel composite endplate.

Carbon Rod Buckling

A support brace composed of carbon rods was considered in order to reduce the larger deflections of composite endplates compared to endplates made from stainless steel. These rods would be used on the front side of the R3 drift chamber and span between the two endplates. The use of 1" diameter hollow carbon fiber rods with 2 mm thick walls in the 6 GeV R3 design did not appear to appreciably impact the tracking capabilities of the drift chamber. The deflection of these carbon fiber rods was modeled as a separate part in ANSYS and compared with hand calculation.

Compression

If one makes the assumption that the rod buckles by an external force along the length of the rod, then the compression of the rod can be calculated using the definition of Young's Modulus (Y) given by the equation:

[math]Y= \frac{\frac{F}{A}}{\frac{\Delta L}{L}}[/math]

where

[math]F =[/math] Compression force exerted along the end of the rod
[math]L =[/math] length of the rod = 2.8 m for the longest rod
[math]A =[/math] cross sections area of the rod = [math]\pi (25 mm)^2[/math]
[math]\Delta L =[/math] rod compression distance = [math]50 \times 10^{-6}[/math]


The force needed to cause the carbon fiber rod to buckle such that the endplate deflects [math]50 \mu m[/math] can be estimated by assuming that the length of the rod (2.8 m) compresses by [math]50 \mu m[/math].

[math]F = \frac{YA \Delta L }{L} = \frac{1.5 \times 10^{11} N/m^2 (\pi 0.000625 m^2) 50 \times 10^{-6} m } {2.8 m}= 5259N[/math]

The total force along the whole endplate from just the drift chamber wire tension amount to less than 3000 N which is less than a factor of 2 smaller than the force needed to compress a single carbon rod used to brace the endplates. Therefore, placing at least 6 carbon rod braces in the drift chamber should provide substantial support for the endplates.

Buckling Load Threshold

The load at which a carbon fiber rod will start to buckle is given by Euler's formula (defined in any undergraduate "Mechanics of Materials" textbook). Figure 4.2.1 depicts the load configuration on the rod for this calculation.

EulersBucklingDiagramWequation.jpg
Figure 4.2.1: Sketch showing the definition of Euler's buckling load threshold.


[math]F_{cr}[/math] represent the critical force needed to cause the rod to begin to buckle. Young's Modulus for a carbon fiber rod is 1.5 x 10^11 N/[math]m^2[/math] and the longest rod's half length (L for the above equation) is 1.4 m. The Area Moment of inertia about the end of a hollow cylinder of inner radius[math] (R_i)[/math] and outer radius [math](R_o)[/math] can be found from:

[math]I_{A,z}= \int_{Ri}^{Ro} r^2 \pi r dr = \frac{\pi}{4} \left( Ro^4 -Ri^4\right) [/math]
[math]=\frac{\pi}{4} \left( (0.0254)^4 -(0.0244)^4\right) = 4.85 * 10^-8 m^4[/math]

The buckling force is found to be 36,646N, an order of magnitude larger than the predicted endplate forces.


The above buckling threshold calculation was performed using ANSYS. The results are shown in Figure 4.2.2 and predict a threshold force of 36989 N compared to the 36646 N value found from hand calculations. ANSYS also predicts that the rod will deflect almost 2 mm if the buckling load is reached, no deflection due to buckling will occur before the load level is reached.


CarbonRodBuckling AnsysCalc1.jpg
Figure 4.2.2. FEA of the maximum (threshold) force needed to cause a carbon rod to buckle.

Actual Buckling FEA

In the previous FEA results, both ends of the carbon rods were constrained from moving from x and y directions and does not account for the deflection of the endplates. If the endplates deflect the ends of the carbon rods will be compressed. The buckling load in this case is observed to be 54720.5N as shown below in Figure 4.3.1. In the event that the buckling load is reached, ANSYS predicts a deflection of at least 11 mm will occur. It should also be noted that the buckling does not have a preferred direction making it possible for the middle of the rod to bend in any direction.

CarbonRodBuckling AnsysCalc2.jpg
Figure 4.3.1. FEA of the maximum (threshold) force needed to cause a carbon rod to buckle with the real constraints.

3-D Analysis

The term "3-D Analysis" refers to an analysis based on a model that includes both endplates, a nose plate, a bottom plate, several carbon fiber rods, and a Hexcel backing plate. The deflection calculation was performed using several permutations involving the wire tension, the presence of gravity, and the use of carbon fiber rods, the Hexcel backing plate as well as the impact of a ball joint in the drift chamber nose plate. Analysis begins with the simplest case in which only the impact of gravitational forces on the endplate is deflection is evaluated for an ionization chamber structure modeled without carbon rods, a Hexcel backing plate, nor a ball joint. The analysis then considers the case in which wires are strung between the endplates and the impact of adding carbon fiber rods, the Hexcell sheet, and the ball joint on the endplate deflection. Figure 14 shows the 3-D model with labels indicating the location of the loadings included in the FEA program. However, to date, none of these analysis cases includes the gas pressure due to the inclusion of the mylar front sheet.

Below is a list describing in detail the applied loading and boundary conditions.
  • The Hinge 1 labeled A is fixed in all degrees of freedom.
  • The Hinge 2 labeled B is fixed in 2 degrees of freedom (DOF) and 1 DOF is allowed to move horizontal.
  • The wire tensions labeled C, D, E, and F are for all four wire super layers (two per endplate).
  • Ball joint labeled G below the Nose plate.
  • Acceleration labeled H is showing upwards, showing the gravity is opposite to it. This gravity vector is applied at an imaginary direction, as shown in Figure 14. The real orientations are different to what is shown here. To obtain the real orientations, the chamber has to be rotated clockwise at 30, 90, 150, 210, 270, 330 degrees, respectively, along z-axis (being the axis of rotation) to consider all the six orientations.


Loads.jpg
Figure 14: Loads and boundary conditions for the FEA 3-D model. This image was taken from the ANSYS FEA program.


Description of Analysis Sets

Each analysis set included a specific component added to the model to slowly increase its complexity. These analysis sets are:

  1. Gravity only loading.
  2. Wire loads and the effect of the ball joint, located on the nose plate.
  3. Wire loading with carbon fiber rods included.
  4. Wire loading with the Hexcel backing sheet included.
  5. Wire loading with both the carbon fiber rods and Hexcel backing sheet included.
  6. Gravity only loading with both the carbon fiber rods and Hexcel backing sheet included.

Critical points to remember about the analysis sets includes:

  • Analysis set 1 and 6 are just self-weight analysis with no external loading.
  • The nose Plate, endplates, back plate, and hinges are made of the sandwiched polyurethane foam/stainless steel composite material.
  • The downstream-side sheet (backign) is made of Hexcel honeycomb material (sandwiched between two carbon fiber sheets).
  • The Young's Modulus of both the endplates and the Hexcel backing are obtained by the effective Young's Modulus calculation method.
  • The compressive modulus of the Hexcel backing is observed to be different along the length and width. As are result, the Young's Modulus is calculated in both directions and the least of them is considered for calculations (more conservative approach).

Summary Table of Results

Table 2 gives a summary of the analysis sets along with the resulting deflections. Column "D Max" is the maximum total deflection (vector) observed in the model; in some sets, this is observed in the carbon fiber rods and not necessarily in the endplates. Since the maximum total deflection is the result of importance, Column "E max" gives this value for the endplates.


Table 2: Summary of 3-D Analysis Results from ANSYS.

Deflection Wire Tension Rods Hexcel Mylar Gravity Ball Joint Gas Pressure D max(mm) X max(mm) Y max(mm) Z max(mm) E max(mm)
Analysis 1.1 No No No No Yes No No 1.5315 0.0034 0.1791 0.0731 1.2
Analysis 1.2 No No No No Yes yes No 0.219 2.11e-5 0.1795 0.0318 0.219
Analysis 2.1 Yes No No No No No No 19.935 17.495 9.5869 0.4435 19.935
Analysis 2.2 Yes No No No No Yes No 19.952 17.391 9.8084 0.4442 19.95
Analysis 2.3 Yes No No No Yes Yes No 18.593 16.27 9.014 0.4813 18.593
Analysis 3.1 Yes Yes No No No No No 1.5104 1.033 0.2855 0.0675 1.17
Analysis 3.2 Yes Yes No No No Yes No 1.158 1.0266 0.5869 0.0534 1.15
Analysis 3.3 Yes Yes No No Yes Yes No 1.038 0.9392 0.3618 0.1018 1.03
Analysis 4.1 Yes No Yes No No No No 0.2904 0.146 0.2479 0.0044 0.29
Analysis 4.2 Yes No Yes No No Yes No 0.2891 0.1538 0.251 0.0166 0.289
Analysis 4.3 Yes No Yes No Yes Yes No 0.301 0.156 0.264 0.0153 0.3
Analysis 5.1 Yes Yes Yes No No No No 0.1231 0.0216 0.0099 0.0031 0.09
Analysis 5.2 Yes Yes Yes No No Yes No 0.0209 0.0157 0.0116 0.0027 0.016
Analysis 5.3 Yes Yes Yes No Yes Yes No 0.0647 0.0078 0.0105 0.0122 0.021
Analysis 6.1 No Yes Yes No Yes No No 0.5851 0.051 0.0052 0.0082 0.52
Analysis 6.2 No Yes Yes No Yes Yes No 0.0738 0.0008 0.0017 0.0098 0.008

Key observations from these results include:

  • The Analysis started with a basic frame made of the composite material (foam and stainless steel). Due to the major effects of gravity on the chamber, this composite material produced good results with respect to self-weight deflection due to its ery low density and considerably high Young's Modulus.
  • Analysis set 1 is just the analysis with gravity effects with no external loading. The very first analysis (1.1) is with hinge constraints and no ball joint. The max deflection is about 1.53 mm. By adding a ball joint (1.2) reduced the deflection to 0.22 mm. With this gravity orientation, the chamber is resting its weight on the ball joint. This tells us that the ball joint is very supportive in reducing the max deflection in every case.
  • Analysis set 2 adds wire tensions to the end plates with and without the ball joint, with and without gravity. The max deflections observed are approximately 19 mm.
  • Analysis set 3 is a step further in trying to reduce the deflection by adding six carbon fiber rods to the upstream side of the chamber. The endplate deflection is reduced to about 1.2 mm.
  • Analysis set 4 includes the effects of the hexcel backing sheet without carbon fiber rods. Under the same conditions, the deflection is reduced from 19 mm (without Hexcel) to 0.3 mm.
  • Analysis set 5 is a combination of sets 3 and 4 (the inclusion of both carbon fiber rods and the Hexcel backing sheet). Without the ball joint, the maximum total deflection of the endplates was 0.09 mm, whereas by adding the ball joint, the deflection is further reduced to 0.016 mm.
  • Analysis set 6 includes gravity effects only with both the carbon fiber rods and the Hexcel backing sheet.

The next step in the analysis process (currently underway) is to include the mylar front sheet and the effects of a possitive gas pressure on both the deformation of the mylar as well as the added loading to the endplates.

Analysis 5.3

The analysis 5.3 has all the components like carbon rods, hexcel backing sheet with wire loads and gravity loads. so, analysis 5.3 is made available in this report and the rest of the analysis sets are available in the appendices section (8) at the end of this report. Every section is described with the supporting images of deflection and a complete report available in pdf at their individual sections


With the addition of gravity forces, the deformation is observed to retain the symmetric shape. Even the back plate is deformed and the bending is observed along the gravity vector. The deformation at the nose plate is very minimal and is about 7 microns. When the gravity is included, the maximum deformation is observed at the very top carbon rod like in every other gravity case. The maximum deformation observed along the end plate is about 21 microns. Figure 15 shows the ANSYS FEA results for this analysis.

A Report of the analysis is available at File:Analysis 5.3.pdf

R3 3D Analysis5.3d.jpg
R3 3D Analysis5.3x.jpg R3 3D Analysis5.3y.jpg R3 3D Analysis5.3z.jpg

Figure 15: Images of Deflection for Analysis 5.3 in D,X,Y,Z Directions respectively

Summary

The deformation observed in the model is mainly a factor of the wire tensions. The wire tensions created the maximum deformation in the mid sections of the end plates. To reduce this effect the carbon rods and the hexcel sheet were included in the geometry. As a result, they have a considerable effect in reducing the deformation of the endplates (from 18.6 mm without carbon fiber rods and a Hexcel backing to 0.021 mm with).

Buckling analysis of the longest carbon fiber rod (i.e. the rod most likely to buckle due to it having the largest slenderness ratio) indicates that the threshold force (i.e. force required to buckle the rod) is 20 times greater than the applied loading. Therefore, buckling of the carbon fiber rods will not be an issue.

The next phase of this analysis is to (1) refine and verify the loading conditions, (2) include the effects of a positive gas pressure inside the chamber (mylar front sheet), and (3) determine mechanical connections between components.

References

  1. [Mestayer07] "Tracking Detectors for CLAS12", 1/9/07, Mac Mestayer, Collaboration Document.
  2. [Mestayer00] M.D. Mestayer, et. al., NIM, A449, pg 81, 2000 Media:NIM_A449.pdf
  3. [Craig96] Mechanics of Materials, Roy R. Craig Jr., John Wiley & Sons, 1996

Appendices

Endplate Deflection with Carbon Rods

Preliminary FEA efforts were performed to determine the impact of including carbon fiber rods on the endplate deflection. This model was only for a single endplate and not the full 3-D Analysis. Figure A1a shows the reference case of a composite endplate deflection undergoing loading from the wires without carbon fiber rods. Figure A1b shows the deflection of the endplate with the six fiber rods included as rigid boundary constraints. Since the rods are only mounted on the upstream side of the endplate, the maximum deflection is then observed on the downstream side.

DtotNoRods.jpg DtotWithRods.jpg
Configuration X max Y max Z max D max
(mm) (mm) (mm) (mm)
No rods .67 32.4 .04 32.4
With rods .01 0.29 0.16 .29


Figure A1: (a, left) FEA predicted deflections for the composite endplate for a reference condition of wire loading witihout carbon fiber rods. (b, right) FEA predicted deflections when carbon fiber rods are included.


Only the dark and light blue regions near the carbon fiber rods appear to have deflections below the [math]50 \mu m[/math] specification. The maximum deflection appears to be on the back (downstream) edge of the end plate. The use of a composite material, with an effective Young's Modulus that is a factor of 20 smaller than stainless steel, resulted in deflections which were only a factor of 16 larger while decreasing the mass of the endplate by a factor of at least 16. Carbon fiber rods along the front of the drift chamber will decrease the maximum deflection to about 300 [math]\mu m[/math] which is a factor of 6 too large.

3-D Analysis

Listed below are the results, including Figures from ANSYS, for the 3-D Analysis not presented in the main text.


Analysis 1

Analysis 1 is an effort to understand the effects of gravity in its orientation with 2 different kinds of constraints

Analysis 1.1

Analysis 1.1 describes the effects of gravity on the frame when hinge 1 and hinge 2 are active with no active ball joint at the nose plate. Wire tensions are not applied i.e no external loading. In this model hexcel sheet and carbon rods are not included. A Report of the analysis is available at File:Analysis 1.1.pdf


As the nose is free to move, the Maximum Total Deflection is observed at the nose and the deflection is gradually reduced towards the hinge side. Due to self weight, the deflection in X-direction is close to zero at those hinges and maximum in the negative side of x. The Y-deflection is more at the middle of the end plates in opposite directions. The Z-deflection is more at the nose plate as the chamber is observed to be tilted towards the Target center in the negative Z-side.

R3 3D Analysis1.1d.jpg
R3 3D Analysis1.1x.jpg R3 3D Analysis1.1y.jpg R3 3D Analysis1.1z.jpg

|- | Figure A.2: Images of Deflection for Analysis 1.1 in D,X,Y,Z Directions respectively

Analysis 1.2

Analysis 1.1 describes the effects of gravity on the frame when hinge 1, hinge 2 and ball joint are active. Wire tensions are not applied i.e no external loading. In this model also, hexcel sheet and carbon rods are not included. A Report of the analysis is available at File:Analysis 1.2.pdf

In this case The nose plate is supported with a ball joint under the nose plate. so, It is like the chamber is made to rest on the ball joint with supported hinges on the top. so, the total deformation at the nose plate and the hinges is observed to be minimum and the end plates gave an impression as they sag in the middle. Other than that, the base plate on the top also sagged below due to its own self weight. Due to gravity effects, the X-direction deflection is observed to be maximum almost at the center of the endplates. The Y-deflection gave an impression as it is symmetric in the way the end plates are bowed outward. Z-deformation is maximum in the direction away away from the target center at the middle of the back plate. The back plate on each piece will not be mounted parallel to the ground. so, the deflection is in the positive z direction.

R3 3D Analysis1.2d.jpg
R3 3D Analysis1.2x.jpg R3 3D Analysis1.2y.jpg R3 3D Analysis1.2z.jpg

|- | Figure A.3: Images of Deflection for Analysis 1.2 in D,X,Y,Z Directions respectively

Analysis 2

This is the step of analysis where the wire tensions are applied to the basic frame where the carbon rods and hexcel sheet is not modeled for the analysis.

Analysis 2.1

A Report of the analysis is available at File:Analysis 2.1.pdf

In this analysis, the deformation is pretty much dominated by the wire tensions. Now the end plates are bowed inward and the base plate has minimum deflection. whatsoever, The nose plate side is slightly pulled towards the target side even without the gravity effects and it can be observed in the Z-deformation image. Overall, the deformation is too huge and is approximately about 20mm. the maximum deformation is observed at the end plates.

R3 3D Analysis2.1d.jpg
R3 3D Analysis2.1x.jpg R3 3D Analysis2.1y.jpg R3 3D Analysis2.1z.jpg

|- | Figure A.4: Images of Deflection for Analysis 2.1 in D,X,Y,Z Directions respectively

Analysis 2.2

A Report of the analysis is available at File:Analysis 2.2.pdf

The deformations observed with the ball joint is pretty close to the deformations observed without a ball joint. however, the deformations are slightly lesser than previous especially at the nose plate with the extra constraint of the ball joint. The maximum deflection is observed along the end plates.

R3 3D Analysis2.2d.jpg
R3 3D Analysis2.2x.jpg R3 3D Analysis2.2y.jpg R3 3D Analysis2.2z.jpg

|- | Figure A.5: Images of Deflection for Analysis 2.2 in D,X,Y,Z Directions respectively

Analysis 2.3

A Report of the analysis is available at File:Analysis 2.3.pdf

In this case the addition of gravity reduced the magnitude of deflection. With out the wire tensions, under gravity the end plates are bowed outward about 2mm in analysis 0.2 and with the wire tensions under no gravity, they are deflected inwards for about 19.5mm in analysis 2.2. Now by putting both the cases together, the deformation is a compromise of both and it is observed as 18.6mm bowing inwards. The maximum deflection is observed along the end plates.

R3 3D Analysis2.3d.jpg
R3 3D Analysis2.3x.jpg R3 3D Analysis2.3y.jpg R3 3D Analysis2.3z.jpg

|- | Figure A.6: Images of Deflection for Analysis 2.3 in D,X,Y,Z Directions respectively

Analysis 3

Addition of carbon rods reduced the deflection to a very great extent, from 19-20mm to about 1.5-1.0mm. The maximum deflection is observed on the middle sections of the carbon rods. They are bowed outwards. the end plate towards the carbon rod side is almost zero and the non carbon rod side of the end plate has 0.8-1.15mm. With the ball joint and the influence of gravity, they are some slight deformations on the back plate towards the gravity. However the role of carbon rods made things quite clear by reducing the deflection the greatest extent.

Analysis 3.1

This analysis is carried with no ball joint and no gravity effects. Six carbon rods were modeled to the upstream side of the chamber. The maximum deformation is observed at one of the mid sections of the carbon rods. As this analysis is without the gravity effects and the nose plate is not constrained with the ball joint, we can see a major deflection at the nose plate. however, the end plates were deflected to about 1.2mm at the non carbon rod end. A Report of the analysis is available at File:Analysis 3.1.pdf

R3 3D Analysis3.1d.jpg
R3 3D Analysis3.1x.jpg R3 3D Analysis3.1y.jpg R3 3D Analysis3.1z.jpg

|- | Figure A.7: Images of Deflection for Analysis 3.1 in D,X,Y,Z Directions respectively

Analysis 3.2

Here, the ball joint reduced the ability of the chamber to move. The maximum deformation is observed at the nose carbon rod edge of the end plate. The maximum deformation is observed to be about 1.15mm. this deflection is more of a factor of the wire tensions. A Report of the analysis is available at File:Analysis 3.2.pdf


R3 3D Analysis3.2d.jpg
R3 3D Analysis3.2x.jpg R3 3D Analysis3.2y.jpg R3 3D Analysis3.2z.jpg

|- | Figure A.8: Images of Deflection for Analysis 3.2 in D,X,Y,Z Directions respectively

Analysis 3.3

This analysis with the gravity and ball joint, gives an expected shape of deformation. The maximum deformation is observed at the mid sections of the top carbon rods and the edges of the end plates. The nose plate has a little deflection of about 0.2mm. The deformation at the hinge section is almost zero. The gravity loads actually helped in reducing the deformation at this orientation. A Report of the analysis is available at File:Analysis 3.3.pdf

R3 3D Analysis3.3d.jpg
R3 3D Analysis3.3x.jpg R3 3D Analysis3.3y.jpg R3 3D Analysis3.3z.jpg

|- | Figure A.9: Images of Deflection for Analysis 3.3 in D,X,Y,Z Directions respectively

Analysis 4

Analysis set 4 is with the hexcel sheet and with no carbon rods. hexcel sheet is mounted on the downstream side of the flow. As expected hexcel sheet reduced the deflection even better than the carbon rods, From 19mm to 0.3mm. The max deflection is observed at the non hexcel sheet side on the end plates. Symmetric deflection is observed. However, the hexcel sheet is bowed inward i.e towards the wire inside the chamber by about. 0.025mm. This can be observed in the Z-deformation images. presence of the ball joint didnt make a lot of difference.

Analysis 4.1

This analysis revealed the importance of hexcel sheet and it reduced the deformation even better than carbon rods. The drawing shows the symmetric deformation with the maximum value at the non hexcel edge of the end plate. The deformation observed on the hexcel sheet is very minimal and that very minimal value is towards the wire in the chamber. A Report of the analysis is available at File:Analysis 4.1.pdf

R3 3D Analysis4.1d.jpg
R3 3D Analysis4.1x.jpg R3 3D Analysis4.1y.jpg R3 3D Analysis4.1z.jpg

|- | Figure A.10: Images of Deflection for Analysis 4.1 in D,X,Y,Z Directions respectively

Analysis 4.2

The presence of ball joint didnt effect the location of maximum deformation at all. It is still observed at the edges of the end plate. The value is about 0.28mm A Report of the analysis is available at File:Analysis 4.2.pdf

R3 3D Analysis4.2d.jpg
R3 3D Analysis4.2x.jpg R3 3D Analysis4.2y.jpg R3 3D Analysis4.2z.jpg

|- | Figure A.11: Images of Deflection for Analysis 4.2 in D,X,Y,Z Directions respectively

Analysis 4.3

The gravity effect at this orientation increased the max deformation along the end plates by about 0.02mm. The maximum deformation is still at the edges of the end plate. A Report of the analysis is available at File:Analysis 4.3.pdf

R3 3D Analysis4.3d.jpg
R3 3D Analysis4.3x.jpg R3 3D Analysis4.3y.jpg R3 3D Analysis4.3z.jpg

|- | Figure A.12: Images of Deflection for Analysis 4.3 in D,X,Y,Z Directions respectively

Analysis 5

The Individual addition of both carbon rods and hexcel sheet in the geometry greatly reduced the maximum deformation. Now, in this analysis Both of them are modeled at the same time. In the very first analysis, without the ball joint, the maximum deflection is observed at the nose plate as it is free to move. however, after adding ball joint, the deflection is observed to be 20 microns at one of the carbon rod bowing away from the wires. The max deflection at the edge of the end plate is about 16 microns. Gravity increased the deflection at the carbon rod and made a very slight change on the deflection at the back plate and along the edge of the end plate.

Analysis 5.1

This analysis has both carbon rods and the hexcel sheet to reduce the deformation. This model is not constrained at the ball joint and the gravity wasn't modeled. So, the maximum deformation is observed at the nose plate which is about 123 microns. The deformation is observed to be reducing to the back plate side along the end plate. The carbon rods were slightly bowed out in the upward direction.

A Report of the analysis is available at File:Analysis 5.1.pdf

R3 3D Analysis5.1d.jpg
R3 3D Analysis5.1x.jpg R3 3D Analysis5.1y.jpg R3 3D Analysis5.1z.jpg

|- | Figure A.13: Images of Deflection for Analysis 5.3 in D,X,Y,Z Directions respectively

Analysis 5.2

The deformation shape is completely changed with the addition of the ball joint. The deformation lost its symmetric values. The maximum total deformation is observed at the carbon rods. The longest carbon rod at the very top is observed to hold the maximum deformation of about 20 microns. However, the maximum deformation observed along the end plate is about 16 microns. A Report of the analysis is available at File:Analysis 5.2.pdf

R3 3D Analysis5.2d.jpg
R3 3D Analysis5.2x.jpg R3 3D Analysis5.2y.jpg R3 3D Analysis5.2z.jpg

|- | Figure A.14: Images of Deflection for Analysis 5.2 in D,X,Y,Z Directions respectively

Analysis 6

This analysis is just to verify the gravity effects with no external loading. The results are very straight forward that the max deflection is at the nose plate when there is no ball joint, and with the ball joint the max deflection is at the top carbon rod as it is the longest of all the carbon rods and is heavier than the other carbon rods. The max deflection is at the center of that carbon rod.

Analysis 6.1

Gravity with no ball joint and with no external loading. The maximum deformation is observed at the nose plate and the deformation is observed to be linearly reducing towards the hinges. A Report of the analysis is available at File:Analysis 6.1.pdf

R3 3D Analysis6.1d.jpg
R3 3D Analysis6.1x.jpg R3 3D Analysis6.1y.jpg R3 3D Analysis6.1z.jpg

|- | Figure A.15: Images of Deflection for Analysis 6.1 in D,X,Y,Z Directions respectively

Analysis 6.2

With the addition of ball joint, the deformation close to nose plate is almost negligible. The maximum deflection is at the top carbon rod which is not of importance. The maximum total deflection at the end plate is about 16 microns A Report of the analysis is available at File:Analysis 6.2.pdf

R3 3D Analysis6.2d.jpg
R3 3D Analysis6.2x.jpg R3 3D Analysis6.2y.jpg R3 3D Analysis6.2z.jpg

|- | Figure A.16: Images of Deflection for Analysis 6.2 in D,X,Y,Z Directions respectively

Go Back