source: trunk/examples/advanced/cosmicray_charging/src/LISAScienceModuleStructures.icc@ 1304

Last change on this file since 1304 was 807, checked in by garnier, 17 years ago

update
File size: 25.2 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26// ********************************************************************
27// * *
28// * cosmicray_charging advanced example for Geant4 *
29// * (adapted simulation of test-mass charging in the LISA mission) *
30// * *
31// * Henrique Araujo (h.araujo@imperial.ac.uk) & Peter Wass *
32// * Imperial College London *
33// * *
34// * LISADetectorConstruction class *
35// * *
36// ********************************************************************
37//
38// HISTORY
39// 22/02/2004: migrated from LISA-V04
40// 24/11/2004: migrated to cmath
41// 08/12/2005: removed compilation warnings
42//
43// ********************************************************************
44
45
46//**************************************************************************
47// Science Module Structure (SMS)
48//**************************************************************************
49// Primary Structure
50// Lower Deck
51// Upper Deck
52// Thermal Shield
53// Solar Array
54// Optical Solar Reflectors (OSR)
55// Radiator Panels
56//**************************************************************************
57
58
59 // Structure Tubes
60 G4double StructTube_dia = 100.*mm;
61 G4double StructTube_thi = 11.*mm;
62 G4double StructTube_hei = 490.*mm;
63 G4double StructTube_sep = 1480.*mm;
64 G4double StructTubeSpacer_dia = 200.*mm;
65 G4double StructTubeSpacer_thi = 50.*mm;
66 G4double StructTubeSpacer_hei = 16.*mm;
67
68 // Diagonal Panels
69 G4double DiagPanelSkin_thi = 0.5*mm;
70 G4double DiagPanelSkin_len = 1420.*mm;
71 G4double DiagPanelSkin_hei = 490.*mm;
72 G4double DiagPanelCore_thi = 19.*mm;
73 G4double DiagPanelCore_len = DiagPanelSkin_len;
74 G4double DiagPanelCore_hei = DiagPanelSkin_hei;
75 G4double DiagPanelCutout_dia = 466.*mm;
76 G4double DiagPanelCutout_off = 50.*mm;
77
78 // Radial Panels
79 G4double RadialPanelSkin_thi = 0.5*mm;
80 G4double RadialPanelSkin_hei = 490.*mm;
81 G4double RadialPanelSkin_top = 424.*mm;
82 G4double RadialPanelSkin_bot = 159.*mm;
83 G4double RadialPanelCore_thi = 19.*mm;
84 G4double RadialPanel_roff = 910.*mm;
85
86 // Upper and Lower Decks
87 G4double UpperDeck_dia = 2700.*mm;
88 G4double UpperDeckCore_thi = 29.0*mm;
89 G4double LowerDeckCore_thi = 29.0*mm;
90 G4double UpperDeck_zoff = 0.5*mm;
91 G4double UpperDeckSkin_thi = 0.5*mm;
92 G4double LowerDeck_dia = 2174.*mm;
93 G4double LowerDeckSkin_thi = 0.5*mm;
94 G4double LowerDeck_zoff = 0.5*mm;
95 G4double MLIBlanket_dia = LowerDeck_dia;
96 G4double MLIBlanket_thi = 1.0*mm;
97
98 // Thermal Shield
99 G4double ThermalShieldFacesheet_thi = 0.6*mm;
100 G4double ThermalShieldCore_thi = 20.*mm;
101 G4double ThermalShieldFoam_thi = 20.*mm;
102 G4double ThermalShield_zoff = StructTubeSpacer_hei;
103 G4double ThermalShield_thi = 3*ThermalShieldFacesheet_thi +
104 ThermalShieldCore_thi +
105 ThermalShieldFoam_thi;
106 // Solar Array and OSR Assy
107 G4double SolarArray_i = 0.5*1683*mm;
108 G4double SolarArray_o = 0.5*UpperDeck_dia;
109 G4double SolarArray_thi = 0.5*mm;
110 G4double OSRAssy_o = SolarArray_i;
111 G4double OSRAssy_thi = 0.5*mm;
112
113 // Radiator Panels
114 G4double RadPanel_top = UpperDeck_dia+12.*mm;
115 G4double RadPanel_bot = LowerDeck_dia+12.*mm;
116 G4double RadPanel_len = 510.*mm;
117 G4double RadPanel_thi = 1.25*mm;
118 // G4double RadPanel_zoff = 20*mm;
119 G4double RadPanelCutout_roff = 0.5*(RadPanel_top+RadPanel_bot)/2.;
120 G4double RadPanelCutout1_dia = 466.*mm;
121 G4double RadPanelCutout1_len = 200.*mm;
122 G4double RadPanelCutout2_hei = 500.*mm;
123 G4double RadPanelCutout2_wid = 400.*mm;
124 G4double RadPanelCutout2_len = 400.*mm;
125 G4double RadPanelCutout2_zoff = 330.*mm;
126
127
128
129
130 // Structure Tubes **********************************************************
131
132 // Struct Tube: 2.16 kg
133 // (ID adjusted to make up for mass of missing top part of tube)
134 G4Tubs* StructTube_sol = new G4Tubs("StructTube_sol",
135 0.5*StructTube_dia-StructTube_thi, 0.5*StructTube_dia,
136 0.5*StructTube_hei, 0., 360.*deg);
137 G4LogicalVolume* StructTube_log = new G4LogicalVolume
138 (StructTube_sol, Al6061, "StructTube_log");
139 G4VPhysicalVolume* StructTube_phys;
140 StructTube_phys = new G4PVPlacement(0, spacecraft_pos +
141 G4ThreeVector(+0.5*StructTube_sep/std::cos(30.*deg),0,0),
142 "StructTube_phys", StructTube_log, wld_phys, false, 0);
143 StructTube_phys = new G4PVPlacement(0, spacecraft_pos +
144 G4ThreeVector(-0.5*StructTube_sep*std::tan(30.*deg),+0.5*StructTube_sep,0),
145 "StructTube_phys", StructTube_log, wld_phys, false, 1);
146 StructTube_phys = new G4PVPlacement(0, spacecraft_pos +
147 G4ThreeVector(-0.5*StructTube_sep*std::tan(30.*deg),-0.5*StructTube_sep,0),
148 "StructTube_phys", StructTube_log, wld_phys, false, 2);
149 StructTube_log->SetVisAttributes(sol_white_vat);
150
151 // Struct Tube Spacer: 1.0 kg
152 // (ID adjusted to include mass of missing TopCap)
153 G4Tubs* StructTubeSpacer_sol = new G4Tubs("StructTubeSpacer_sol",
154 0.5*StructTubeSpacer_dia-StructTubeSpacer_thi,
155 0.5*StructTubeSpacer_dia, 0.5*StructTubeSpacer_hei, 0., 360.*deg);
156 G4LogicalVolume* StructTubeSpacer_log = new G4LogicalVolume
157 (StructTubeSpacer_sol, Al6061, "StructTubeSpacer_log");
158 G4VPhysicalVolume* StructTubeSpacer_phys;
159 StructTubeSpacer_phys = new G4PVPlacement
160 (0, spacecraft_pos + G4ThreeVector(+0.5*StructTube_sep/std::cos(30.*deg),0,
161 +0.5*StructTube_hei+UpperDeckCore_thi+2*UpperDeckSkin_thi
162 +0.5*StructTubeSpacer_hei+UpperDeck_zoff),
163 "StructTubeSpacer_phys", StructTubeSpacer_log, wld_phys, false, 0);
164 StructTubeSpacer_phys = new G4PVPlacement
165 (0, spacecraft_pos + G4ThreeVector(-0.5*StructTube_sep*std::tan(30.*deg),
166 +0.5*StructTube_sep, +0.5*StructTube_hei+UpperDeckCore_thi
167 +2*UpperDeckSkin_thi+0.5*StructTubeSpacer_hei+UpperDeck_zoff),
168 "StructTubeSpacer_phys", StructTubeSpacer_log, wld_phys, false, 1);
169 StructTubeSpacer_phys = new G4PVPlacement
170 (0, spacecraft_pos + G4ThreeVector(-0.5*StructTube_sep*std::tan(30.*deg),
171 -0.5*StructTube_sep,+0.5*StructTube_hei+UpperDeckCore_thi
172 +2*UpperDeckSkin_thi+0.5*StructTubeSpacer_hei+UpperDeck_zoff),
173 "StructTubeSpacer_phys", StructTubeSpacer_log, wld_phys, false, 2);
174 StructTubeSpacer_log->SetVisAttributes(sol_grey_vat);
175
176
177 // Diagonal Panels **********************************************************
178
179 // Diagonal Panel: core (x3)
180 // Honeycomb core is represented by Al6061 at 0.05 g/cm3
181 G4Box* DiagPanelCore_box = new G4Box("DiagPanelCore_box",
182 0.5*DiagPanelCore_len-2.*cm, 0.5*DiagPanelCore_thi, 0.5*DiagPanelCore_hei);
183 G4Tubs* DiagPanelCutout_tub = new G4Tubs("DiagPanelCutout_tub", 0.,
184 0.5*DiagPanelCutout_dia, 2*DiagPanelCore_thi, 0., 360.*deg);
185 G4UnionSolid* DiagPanelCutout_sol = new G4UnionSolid("DiagPanelCutout_sol",
186 DiagPanelCutout_tub, DiagPanelCutout_tub,
187 G4Transform3D(G4RotationMatrix(), G4ThreeVector(DiagPanelCutout_off,0,0)));
188 G4RotationMatrix SMS_rot1; SMS_rot1.rotateX(90.*deg);
189 G4SubtractionSolid* DiagPanelCore_sol =
190 new G4SubtractionSolid("DiagPanelCore_sol", DiagPanelCore_box,
191 DiagPanelCutout_sol, G4Transform3D(SMS_rot1,
192 G4ThreeVector(-0.5*DiagPanelCutout_off,0,0)));
193 G4LogicalVolume* DiagPanelCore_log = new G4LogicalVolume
194 (DiagPanelCore_sol, AlHoneycomb, "DiagPanelCore_log");
195 G4VPhysicalVolume* DiagPanelCore_phys;
196 G4RotationMatrix SMS_rot2; SMS_rot2.rotateZ(90.*deg);
197 DiagPanelCore_phys = new G4PVPlacement(G4Transform3D(SMS_rot2,
198 spacecraft_pos + G4ThreeVector(-0.5*DiagPanelCore_len*std::tan(30.*deg),0,0)),
199 "DiagPanelCore_phys", DiagPanelCore_log, wld_phys, false, 0);
200 G4RotationMatrix SMS_rot3; SMS_rot3.rotateZ(-30.*deg);
201 DiagPanelCore_phys = new G4PVPlacement(G4Transform3D(SMS_rot3,
202 spacecraft_pos + G4ThreeVector
203 (.5*DiagPanelCore_len*std::sin(30.*deg)*std::tan(30.*deg),
204 0.5*DiagPanelCore_len*std::sin(30.*deg),0)), "DiagPanelCore_phys",
205 DiagPanelCore_log, wld_phys, false, 1);
206 G4RotationMatrix SMS_rot4; SMS_rot4.rotateZ(+30.*deg);
207 DiagPanelCore_phys = new G4PVPlacement(G4Transform3D(SMS_rot4,
208 spacecraft_pos + G4ThreeVector
209 (0.5*DiagPanelCore_len*std::sin(30*deg)*std::tan(30*deg),
210 -0.5*DiagPanelCore_len*std::sin(30.*deg),0)),
211 "DiagPanelCore_phys", DiagPanelCore_log, wld_phys, false, 2);
212 DiagPanelCore_log->SetVisAttributes(sol_white_vat);
213
214 // Diagonal Panel: skins (x6)
215 G4Box* DiagPanelSkin_box = new G4Box("DiagPanelSkin_box",
216 0.5*DiagPanelSkin_len-2.*cm, 0.5*DiagPanelSkin_thi, 0.5*DiagPanelSkin_hei);
217 G4SubtractionSolid* DiagPanelSkin_sol =
218 new G4SubtractionSolid("DiagPanelSkin_sol", DiagPanelSkin_box,
219 DiagPanelCutout_sol, G4Transform3D(SMS_rot1,
220 G4ThreeVector(-0.5*DiagPanelCutout_off,0,0)));
221 G4LogicalVolume* DiagPanelSkin_log = new G4LogicalVolume
222 (DiagPanelSkin_sol, CFRP, "DiagPanelSkin_log");
223 G4VPhysicalVolume* DiagPanelSkin_phys;
224 DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot2,
225 spacecraft_pos + G4ThreeVector(-0.5*DiagPanelCore_len*std::tan(30.*deg) +
226 0.5*(DiagPanelCore_thi+DiagPanelSkin_thi),0,0)),
227 "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 0);
228 DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot2,
229 spacecraft_pos + G4ThreeVector(-0.5*DiagPanelCore_len*std::tan(30.*deg) -
230 0.5*(DiagPanelCore_thi+DiagPanelSkin_thi),0,0)),
231 "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 1);
232 DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot3,
233 spacecraft_pos + G4ThreeVector
234 (0.5*DiagPanelCore_len*std::sin(30*deg)*std::tan(30*deg),
235 0.5*DiagPanelCore_len*std::sin(30.*deg),0) +
236 G4ThreeVector(0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::sin(30.*deg),
237 0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::cos(30.*deg),0)),
238 "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 2);
239 DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot3,
240 spacecraft_pos + G4ThreeVector
241 (0.5*DiagPanelCore_len*std::sin(30*deg)*std::tan(30*deg),
242 0.5*DiagPanelCore_len*std::sin(30.*deg),0) +
243 G4ThreeVector(-0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::sin(30.*deg),
244 -0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::cos(30.*deg),0)),
245 "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 3);
246 DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot4,
247 spacecraft_pos + G4ThreeVector
248 (0.5*DiagPanelCore_len*std::sin(30*deg)*std::tan(30*deg),
249 -0.5*DiagPanelCore_len*std::sin(30.*deg),0) +
250 G4ThreeVector(0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::sin(30.*deg),
251 -0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::cos(30.*deg),0)),
252 "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 3);
253 DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot4,
254 spacecraft_pos + G4ThreeVector
255 (0.5*DiagPanelCore_len*std::sin(30*deg)*std::tan(30*deg),
256 -0.5*DiagPanelCore_len*std::sin(30.*deg),0) +
257 G4ThreeVector(-0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::sin(30.*deg),
258 +0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::cos(30.*deg),0)),
259 "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 4);
260 DiagPanelSkin_log->SetVisAttributes(sol_white_vat);
261
262
263
264 // Radial Panels ************************************************************
265
266 // Radial Panels: skins (x3)
267 G4Trap* RadialPanelSkin_sol = new G4Trap("RadialPanelSkin_sol",
268 2*RadialPanelSkin_thi+RadialPanelCore_thi, RadialPanelSkin_hei,
269 RadialPanelSkin_top, RadialPanelSkin_bot);
270 G4LogicalVolume* RadialPanelSkin_log = new G4LogicalVolume
271 (RadialPanelSkin_sol, CFRP, "RadialPanelSkin_log");
272 G4VPhysicalVolume* RadialPanelSkin_phys;
273 G4RotationMatrix SMS_rot5; SMS_rot5.rotateX(-90.*deg);
274 RadialPanelSkin_phys = new G4PVPlacement
275 (G4Transform3D(SMS_rot5, spacecraft_pos +
276 G4ThreeVector(+RadialPanel_roff/std::cos(30.*deg),0,0)),
277 "RadialPanelSkin_phys", RadialPanelSkin_log, wld_phys, false, 0);
278 G4RotationMatrix SMS_rot6;
279 SMS_rot6.rotateX(-90.*deg); SMS_rot6.rotateZ(120.*deg);
280 RadialPanelSkin_phys = new G4PVPlacement
281 (G4Transform3D(SMS_rot6, spacecraft_pos +
282 G4ThreeVector(-RadialPanel_roff*std::tan(30.*deg),+RadialPanel_roff,0)),
283 "RadialPanelSkin_phys", RadialPanelSkin_log, wld_phys, false, 1);
284 G4RotationMatrix SMS_rot7;
285 SMS_rot7.rotateX(-90.*deg); SMS_rot7.rotateZ(-120.*deg);
286 RadialPanelSkin_phys = new G4PVPlacement
287 (G4Transform3D(SMS_rot7, spacecraft_pos +
288 G4ThreeVector(-RadialPanel_roff*std::tan(30.*deg),-RadialPanel_roff,0)),
289 "RadialPanelSkin_phys", RadialPanelSkin_log, wld_phys, false, 2);
290 RadialPanelSkin_log->SetVisAttributes(sol_white_vat);
291 // Radial Panels: core
292 // Honeycomb core is represented by Al6061 at 0.05 g/cm3
293 G4Trap* RadialPanelCore_sol = new G4Trap("RadialPanelCore_sol",
294 RadialPanelCore_thi, RadialPanelSkin_hei,
295 RadialPanelSkin_top, RadialPanelSkin_bot);
296 G4LogicalVolume* RadialPanelCore_log = new G4LogicalVolume
297 (RadialPanelCore_sol, AlHoneycomb, "RadialPanelCore_log");
298 G4VPhysicalVolume* RadialPanelCore_phys;
299 RadialPanelCore_phys = new G4PVPlacement
300 (0, G4ThreeVector(0,0,0), "RadialPanelCore_phys",
301 RadialPanelCore_log, RadialPanelSkin_phys, false, 0);
302 RadialPanelCore_log->SetVisAttributes(sol_white_vat);
303
304
305
306 // Upper and Lower Decks ****************************************************
307
308
309 // Upper Deck - skin: 2*4.56 kg
310 G4Tubs* UpperDeckSkin_sol = new G4Tubs
311 ("UpperDeckSkin_sol", 0., 0.5*UpperDeck_dia,
312 0.5*UpperDeckCore_thi + UpperDeckSkin_thi, 0., 360.*deg);
313 G4LogicalVolume* UpperDeckSkin_log = new G4LogicalVolume
314 (UpperDeckSkin_sol, CFRP, "UpperDeckSkin_log");
315 G4VPhysicalVolume* UpperDeckSkin_phys;
316 UpperDeckSkin_phys = new G4PVPlacement
317 (0, spacecraft_pos + G4ThreeVector(0,0,+0.5*StructTube_hei
318 +0.5*UpperDeckCore_thi+UpperDeckSkin_thi+UpperDeck_zoff),
319 "UpperDeckSkin_phys", UpperDeckSkin_log, wld_phys, false, 0);
320 UpperDeckSkin_log->SetVisAttributes(white_vat);
321 // upper deck - core: 8.3 kg
322 // Honeycomb core is represented by Al6061 at 0.05 g/cm3
323 G4Tubs* UpperDeckCore_sol = new G4Tubs
324 ("UpperDeckCore_sol",0.,.5*UpperDeck_dia,.5*UpperDeckCore_thi,0.,360.*deg);
325 G4LogicalVolume* UpperDeckCore_log = new G4LogicalVolume
326 (UpperDeckCore_sol, AlHoneycomb, "UpperDeckCore_log");
327 G4VPhysicalVolume* UpperDeckCore_phys;
328 UpperDeckCore_phys = new G4PVPlacement(0,G4ThreeVector(),
329 "UpperDeckCore_phys", UpperDeckCore_log, UpperDeckSkin_phys, false, 0);
330 UpperDeckCore_log->SetVisAttributes(white_vat);
331
332 // Lower Deck - skin: 2*2.95 kg
333 G4Tubs* LowerDeckSkin_sol = new G4Tubs
334 ("LowerDeckSkin_sol", 0., 0.5*LowerDeck_dia,
335 0.5*LowerDeckCore_thi + LowerDeckSkin_thi, 0., 360.*deg);
336 G4LogicalVolume* LowerDeckSkin_log = new G4LogicalVolume
337 (LowerDeckSkin_sol, CFRP, "LowerDeckSkin_log");
338 G4VPhysicalVolume* LowerDeckSkin_phys;
339 LowerDeckSkin_phys = new G4PVPlacement
340 (0, spacecraft_pos + G4ThreeVector(0,0,-0.5*StructTube_hei
341 -0.5*LowerDeckCore_thi-LowerDeckSkin_thi-LowerDeck_zoff),
342 "LowerDeckSkin_phys", LowerDeckSkin_log, wld_phys, false, 0);
343 LowerDeckSkin_log->SetVisAttributes(sol_dgrey_vat);
344 // lower deck - core: 5.38 kg
345 // Honeycomb core is represented by Al6061 at 0.05 g/cm3
346 G4Tubs* LowerDeckCore_sol = new G4Tubs
347 ("LowerDeckCore_sol",0.,.5*LowerDeck_dia,.5*LowerDeckCore_thi,0.,360.*deg);
348 G4LogicalVolume* LowerDeckCore_log = new G4LogicalVolume
349 (LowerDeckCore_sol, AlHoneycomb, "LowerDeckCore_log");
350 G4VPhysicalVolume* LowerDeckCore_phys;
351 LowerDeckCore_phys = new G4PVPlacement(0,G4ThreeVector(),
352 "LowerDeckCore_phys", LowerDeckCore_log, LowerDeckSkin_phys, false, 0);
353 LowerDeckCore_log->SetVisAttributes(sol_dgrey_vat);
354
355 // Lower Deck MLI Blanket: 5.3 kg
356 G4Tubs* MLIBlanket_sol = new G4Tubs("MLIBlanket_sol",
357 0., 0.5*MLIBlanket_dia, 0.5*MLIBlanket_thi, 0., 360.*deg);
358 G4LogicalVolume* MLIBlanket_log = new G4LogicalVolume
359 (MLIBlanket_sol, MLImat, "MLIBlanket_log");
360 G4VPhysicalVolume* MLIBlanket_phys;
361 MLIBlanket_phys = new G4PVPlacement
362 (0, spacecraft_pos + G4ThreeVector(0,0,-0.5*StructTube_hei
363 -LowerDeckCore_thi-2*LowerDeckSkin_thi-LowerDeck_zoff-.5*MLIBlanket_thi),
364 "MLIBlanket_phys", MLIBlanket_log, wld_phys, false, 0);
365 MLIBlanket_log->SetVisAttributes(sol_orange_vat);
366
367
368 // Thermal Shield ***********************************************************
369
370 // Thermal Shield (Facesheet+HoneycombCore+Facesheet+Foam+Facesheet)
371 // Thermal Shield - facesheet: 3*5.47 kg
372 G4Tubs* ThermalShieldFacesheet_sol = new G4Tubs("ThermalShieldFacesheet_sol",
373 0., 0.5*UpperDeck_dia, 0.5*ThermalShield_thi, 0., 360.*deg);
374 G4LogicalVolume* ThermalShieldFacesheet_log = new G4LogicalVolume
375 (ThermalShieldFacesheet_sol, CFRP, "ThermalShieldFacesheet_log");
376 G4VPhysicalVolume* ThermalShieldFacesheet_phys;
377 ThermalShieldFacesheet_phys = new G4PVPlacement(0,
378 spacecraft_pos + G4ThreeVector
379 (0,0,+.5*StructTube_hei+UpperDeckCore_thi+2*UpperDeckSkin_thi+
380 UpperDeck_zoff+0.5*ThermalShield_thi+ThermalShield_zoff),
381 "ThermalShieldFacesheet_phys",ThermalShieldFacesheet_log,wld_phys,false,0);
382 ThermalShieldFacesheet_log->SetVisAttributes(red_vat);
383 // Thermal Shield: core
384 // Honeycomb core is represented by Al6061 at 0.05 g/cm3
385 G4Tubs* ThermalShieldCore_sol = new G4Tubs("ThermalShieldCore_sol",
386 0., 0.5*UpperDeck_dia, 0.5*ThermalShieldCore_thi, 0., 360.*deg);
387 G4LogicalVolume* ThermalShieldCore_log = new G4LogicalVolume
388 (ThermalShieldCore_sol, AlHoneycomb, "ThermalShieldCore_log");
389 G4VPhysicalVolume* ThermalShieldCore_phys;
390 ThermalShieldCore_phys = new G4PVPlacement(0,
391 G4ThreeVector(0,0,.5*ThermalShieldFacesheet_thi+.5*ThermalShieldCore_thi),
392 "ThermalShieldCore_phys", ThermalShieldCore_log,
393 ThermalShieldFacesheet_phys, false, 0);
394 ThermalShieldCore_log->SetVisAttributes(red_vat);
395 // Thermal Shield: foam
396 G4Tubs* ThermalShieldFoam_sol = new G4Tubs("ThermalShieldFoam_sol",
397 0., 0.5*UpperDeck_dia, 0.5*ThermalShieldFoam_thi, 0., 360.*deg);
398 G4LogicalVolume* ThermalShieldFoam_log = new G4LogicalVolume
399 (ThermalShieldFoam_sol, foam, "ThermalShieldFoam_log");
400 G4VPhysicalVolume* ThermalShieldFoam_phys;
401 ThermalShieldFoam_phys = new G4PVPlacement(0,
402 G4ThreeVector(0,0,-.5*ThermalShieldFacesheet_thi-.5*ThermalShieldFoam_thi),
403 "ThermalShieldFoam_phys", ThermalShieldFoam_log,
404 ThermalShieldFacesheet_phys, false, 0);
405 ThermalShieldFoam_log->SetVisAttributes(red_vat);
406
407
408
409 // Solar Array and OSR Assy *************************************************
410
411 // Solar Array: 13.5 kg
412 // Thickness for correct mass
413 G4Tubs* SolarArray_sol = new G4Tubs("SolarArray_sol",
414 SolarArray_i, SolarArray_o, 0.5*SolarArray_thi, 0., 360.*deg);
415 G4LogicalVolume* SolarArray_log = new G4LogicalVolume
416 (SolarArray_sol, Scell, "SolarArray_log");
417 G4VPhysicalVolume* SolarArray_phys;
418 SolarArray_phys = new G4PVPlacement(0, spacecraft_pos +
419 G4ThreeVector(0,0,+.5*StructTube_hei+UpperDeckCore_thi+2*UpperDeckSkin_thi+
420 UpperDeck_zoff+ThermalShield_thi+ThermalShield_zoff+0.5*SolarArray_thi),
421 "SolarArray_phys",SolarArray_log, wld_phys, false, 0);
422 SolarArray_log->SetVisAttributes(lblue_vat);
423 // Optical Solar Reflector: NO MASS
424 G4Tubs* OSRAssy_sol = new G4Tubs("OSRAssy_sol",
425 0., OSRAssy_o, 0.5*OSRAssy_thi, 0., 360.*deg);
426 G4LogicalVolume* OSRAssy_log = new G4LogicalVolume
427 (OSRAssy_sol, vacuum, "OSRAssy_log");
428 G4VPhysicalVolume* OSRAssy_phys;
429 OSRAssy_phys = new G4PVPlacement(0, spacecraft_pos +
430 G4ThreeVector(0,0,+.5*StructTube_hei+UpperDeckCore_thi+2*UpperDeckSkin_thi+
431 UpperDeck_zoff+ThermalShield_thi+ThermalShield_zoff+0.5*OSRAssy_thi),
432 "OSRAssy_phys",OSRAssy_log, wld_phys, false, 0);
433 OSRAssy_log->SetVisAttributes(lblue_vat);
434
435
436
437 // Radiator Panels **********************************************************
438
439 // Radiator Panels: 3*3.87 kg
440 // Thickess adjusted for correct mass of FS+Core+FS
441 G4Cons* RadPanel_con = new G4Cons("RadPanel_con",
442 0.5*RadPanel_bot-RadPanel_thi, 0.5*RadPanel_bot,
443 0.5*RadPanel_top-RadPanel_thi, 0.5*RadPanel_top,
444 0.5*RadPanel_len, 0.*deg, 360.*deg);
445 // window cutouts
446 G4Tubs* RadPanelCutout_tub = new G4Tubs("RadPanelCutout_tub", 0.,
447 0.5*RadPanelCutout1_dia, 2*RadPanelCutout1_len, 0., 360.*deg);
448 G4RotationMatrix SMS_rot8; SMS_rot8.rotateY(90.*deg);
449 G4SubtractionSolid* RadPanel1_sol =
450 new G4SubtractionSolid("RadPanel1_sol", RadPanel_con,
451 RadPanelCutout_tub, G4Transform3D(SMS_rot8,
452 G4ThreeVector(-RadPanelCutout_roff,0,0)));
453 G4RotationMatrix SMS_rot9;
454 SMS_rot9.rotateY(90.*deg); SMS_rot9.rotateZ(49.*deg);
455 G4SubtractionSolid* RadPanel2_sol =
456 new G4SubtractionSolid("RadPanel2_sol", RadPanel1_sol,
457 RadPanelCutout_tub, G4Transform3D(SMS_rot9, G4ThreeVector
458 (RadPanelCutout_roff*std::sin(49.*deg),RadPanelCutout_roff*std::cos(49.*deg),0)));
459 G4RotationMatrix SMS_rot10;
460 SMS_rot10.rotateY(90.*deg); SMS_rot10.rotateZ(131.*deg);
461 G4SubtractionSolid* RadPanel3_sol =
462 new G4SubtractionSolid("RadPanel3_sol", RadPanel2_sol,
463 RadPanelCutout_tub, G4Transform3D(SMS_rot10, G4ThreeVector
464 (RadPanelCutout_roff*std::sin(131.*deg),RadPanelCutout_roff*std::cos(131.*deg),0)));
465 // box cutouts
466 G4Box* RadPanelCutout_box = new G4Box("RadPanelCutout_box",
467 0.5*RadPanelCutout2_len, 0.5*RadPanelCutout2_wid, 0.5*RadPanelCutout2_hei);
468 G4SubtractionSolid* RadPanel4_sol =
469 new G4SubtractionSolid("RadPanel4_sol", RadPanel3_sol,
470 RadPanelCutout_box, G4Transform3D(G4RotationMatrix(), G4ThreeVector
471 (RadPanelCutout_roff,0,RadPanelCutout2_zoff)));
472 G4RotationMatrix SMS_rot11; SMS_rot11.rotateZ(+120.*deg);
473 G4SubtractionSolid* RadPanel5_sol =
474 new G4SubtractionSolid("RadPanel5_sol", RadPanel4_sol,
475 RadPanelCutout_box, G4Transform3D(SMS_rot11,
476 G4ThreeVector(RadPanelCutout_roff*std::cos(120.*deg),
477 RadPanelCutout_roff*std::sin(120.*deg), RadPanelCutout2_zoff)));
478 G4RotationMatrix SMS_rot12; SMS_rot12.rotateZ(-120.*deg);
479 G4SubtractionSolid* RadPanel6_sol =
480 new G4SubtractionSolid("RadPanel6_sol", RadPanel5_sol,
481 RadPanelCutout_box, G4Transform3D(SMS_rot12,
482 G4ThreeVector(RadPanelCutout_roff*std::cos(-120.*deg),
483 RadPanelCutout_roff*std::sin(-120.*deg), RadPanelCutout2_zoff)));
484 G4LogicalVolume* RadPanel_log =
485 new G4LogicalVolume(RadPanel6_sol, Al6061, "RadPanel_log");
486 G4VPhysicalVolume* RadPanel_phys;
487 RadPanel_phys =
488 new G4PVPlacement(0, spacecraft_pos + G4ThreeVector(0,0,0),
489 "RadPanel_phys", RadPanel_log, wld_phys, false, 0);
490 RadPanel_log->SetVisAttributes(sol_white_vat);
491
492//*****************************************************************************
493
494
495 // invisible spacecraft top
496 if(0) {
497 UpperDeckSkin_log->SetVisAttributes(G4VisAttributes::Invisible);
498 UpperDeckCore_log->SetVisAttributes(G4VisAttributes::Invisible);
499 ThermalShieldFacesheet_log->SetVisAttributes(G4VisAttributes::Invisible);
500 ThermalShieldCore_log->SetVisAttributes(G4VisAttributes::Invisible);
501 ThermalShieldFoam_log->SetVisAttributes(G4VisAttributes::Invisible);
502 // SolarArray_log->SetVisAttributes(G4VisAttributes::Invisible);
503 OSRAssy_log->SetVisAttributes(G4VisAttributes::Invisible);
504 }
505
506 // solid spacecraft top
507 if(0) {
508 UpperDeckSkin_log->SetVisAttributes(sol_white_vat);
509 UpperDeckCore_log->SetVisAttributes(sol_white_vat);
510 ThermalShieldFacesheet_log->SetVisAttributes(sol_red_vat);
511 ThermalShieldCore_log->SetVisAttributes(sol_red_vat);
512 ThermalShieldFoam_log->SetVisAttributes(sol_red_vat);
513 SolarArray_log->SetVisAttributes(sol_lblue_vat);
514 OSRAssy_log->SetVisAttributes(sol_lgreen_vat);
515 }
516
517//*****************************************************************************
Note: See TracBrowser for help on using the repository browser.