Modeling a Composite Slot Cross-Section for Torsional Analysis

Transcription

1 Modeling a Composite Slot Cross-Section for Torsional Analysis The cross-section in question is shown below (see also p. 432 in the textbook). Due to double symmetry, only one-quarter of the cross-section will be modeled as shown below. 1 V. N. Kaliakin

2 To creat this model using PATRAN requires the following actions. Start the PATRAN application. Create a suitable new database (journal) file (e.g. filled_slot.db) in which to record all commands and actions. Click on the Geometry icon. Using Create -> Points -> XYZ, define the coordinates of the nine (9) points shown above. These are used to describe the four (4) recatngular regions that the domain is divided into. Using Create -> Curve -> 2 Point, define the following straight lines:1-2, 2-5,5-4,4-1, 2-3,3-6, 6-5, 5-8, 8-7, 7-4, 6-9, and 9-8. Using Create -> Surface -> Edge -> 4 Edge, create the following four (4) surfaces: surface 1 (curve 1, curve 2, curve 3, curve 4), surface 2 (curve 5, curve 6, curve 7, curve 2), surface 3 (curve 7, curve 11, curve 12, curve 8), surface 4 (curve 3, curve 8, curve 9, curve 10). Note that the curves are selected sequentially (in this case counterclockwise). Click on the Elements icon. Using Create -> Mesh Seed -> Uniform, proceed to define the number of elements desired along given curves (lines) in the various surfaces. In prticular, 3 elements along curve 1, 3 elements along curve 9, 4 elements along curve 10, 6 elements along curve 4, and 8 elements along curves 5 and 12. Using Create -> Mesh -> Surface, create the meshes in surface 1, surface 2, surface 3, and finally in surface 4. In this case quad4 (linear quadrilateral elements have been selected). Since four separate rectangular surfaces have been meshed, duplicate nodes have been created at the same location along common surface boundaries. To condense these common nodes (i.e., do delete unnecessary duplicate nodes) select the Equivalence -> All -> Tolerance Cube command with the default tolerance and press Apply. To next optimize the nodal numbering, select the Optimize -> Nodes -> Cuthill- McKee command with the Profile radio button selected (UD_scalar emplys a profile solver) and press Apply. A window will pop up with the old and new bandwidth, profile, etc. displayed. Close this window by presssing OK. Click on the Loads/BCs icon. Using Create -> Temperature -> Nodal, decsribe the homogeneous essential specificatiosn aloing the external boundary. In particular, choose a New Set Name (e.g., external_boundary). Next press the Input Data button and enter a value of 0.0 for Temperature and press OK. Press the Select Application Region button and select curves 6, 9, 11 and 12. Press the OK button. Once back in the main menu, press Apply. Some red markers will next appear at vertex nodes along the chosen boundary. 2 V. N. Kaliakin

3 Click on the Properties icon. Create -> 2D -> Shell, decsribe and apply the two material types to the appropriate surfaces. In particular, enter a Proprty Set Name (e.g., material_1). Then press the Input Proerties button and fill in the material name and thickness (equal to 1.0) and press OK. Then proceed to Select Members and choose the appropriate surfaces (surfaces 2, 3, 4 for material 1). Press Apply. Repeat this procedure for the second material (applicable to surface 1). Next write the Patran neutral file by going to the File menu and scrolling down to Export. Make sure that the format selected is indeed a neutral file and then name the file (e.g., filled_slot.out) and save it to an appropriate directory. The neutral file must next be translated to UD_scalar format. To do this run the pat2ud_scalar program. The input file asked for should be the Patran neutral file (e.g., filled_slot.out). The output file name supplied must be something different from the filename used for the neutral file (e.g., filled_slot.dat). The resulting data file (e.g., filled_slot.dat) requires a bit of editing before it can be analyzed using UD_scalar. In particular, Check that the proper element type has been specified in the ELEMENT commands appearing in the input file filled_slot.dat The number of elements needs to be properly specified in the DIMENSION command block. Although the number of elements is given in the translated file, the proper DIMENSION command must reflect the proper type of element used (in this case four-node quadrilateral elements). As such, the following command must be supplied: DIMENSION MAX QS4 110, where the total number of elements (110) value was provided by Patran. The number of material descriptions must also be specified in the DIMENSION command block. Since two different materials are used in the model, it follows that the following command must be supplied: DIMENSION MAX scalar_1 2 The heat source and material descriptions must next follow the FINISHED SETTINGS command. The former is equal to the product 2GQ, while the latter must describe an identity coefficient matrix. Since two different materials are used in the model, it follows that the following command must be supplied: scalar source number 1 source_ scalar source number 2 source_ V. N. Kaliakin

4 ! scalar conduc constant number 1 & desc "this is material #1" k k scalar conduc constant number 2 & desc "this is material #2" k k For brevity, the input file associated with a 20-element mesh is shown on the next page (NOTE: some lines of input have been wrapped by the word processing software). 4 V. N. Kaliakin

5 ana title "torsion of a 'slot' cross-section" ana title "with TWO different materials" ana title "20 element mesh of QS4 elements"! ANALYSIS ACTION ANALYZE ANALYSIS IDEALIZATION PLANE_STRAIN ANALYSIS TEMPORAL STATIC! ECHO WARNINGS OFF! dim max scalar_1 2 dim max nodes = 30 dim max qs4 = 20! FINISHED SETTINGS! scalar source number 1 source_ ! scalar source number 2 source_ ! scalar conduc constant number 1 & desc "this is material #1" k k22 1.0! scalar conduc constant number 2 & desc "this is material #2" k k22 1.0! NODES LINE NUMBER 1 x1 x E+00 x3 NODES LINE NUMBER 2 x E-01 x E+00 x3 NODES LINE NUMBER 3 x E+00 x E+00 x3 NODES LINE NUMBER 4 x E+00 x E+00 x3 NODES LINE NUMBER 5 x E+00 x E+00 x3 NODES LINE NUMBER 6 x1 x E+00 x3 NODES LINE NUMBER 7 x E-01 x E+00 x3 NODES LINE NUMBER 8 x E+00 x E+00 x3 NODES LINE NUMBER 9 x E+00 x E+00 x3 NODES LINE NUMBER 10 x E+00 x E+00 x3 NODES LINE NUMBER 11 x1 x E+00 x3 NODES LINE NUMBER 12 x1 x E+00 x3 NODES LINE NUMBER 13 x1 x E+00 x3 NODES LINE NUMBER 14 x1 x2 x3 NODES LINE NUMBER 15 x E-01 x E+00 x3 5 V. N. Kaliakin

6 NODES LINE NUMBER 16 x E+00 x E+00 x3 NODES LINE NUMBER 17 x E+00 x E+00 x3 NODES LINE NUMBER 18 x E+00 x E+00 x3 NODES LINE NUMBER 19 x E-01 x E+00 x3 NODES LINE NUMBER 20 x E-01 x E+00 x3 NODES LINE NUMBER 21 x E-01 x2 x3 NODES LINE NUMBER 22 x E+00 x E+00 x3 NODES LINE NUMBER 23 x E+00 x E+00 x3 NODES LINE NUMBER 24 x E+00 x E+00 x3 NODES LINE NUMBER 25 x E+00 x E+00 x3 NODES LINE NUMBER 26 x E+00 x2 x3 NODES LINE NUMBER 27 x E+00 x E+00 x3 NODES LINE NUMBER 28 x E+00 x E+00 x3 NODES LINE NUMBER 29 x E+00 x2 x3 NODES LINE NUMBER 30 x E+00 x2 x3! ELEM SCALAR TYPE QS4 SCALAR 2 NODES ELEM SCALAR TYPE QS4 SCALAR 2 NODES ELEM SCALAR TYPE QS4 SCALAR 2 NODES ELEM SCALAR TYPE QS4 SCALAR 2 NODES ELEM SCALAR TYPE QS4 SCALAR 2 NODES ELEM SCALAR TYPE QS4 SCALAR 2 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ELEM SCALAR TYPE QS4 SCALAR 1 NODES ! SPEC CONC SCALAR NODE 1 PHI VALUE SPEC CONC SCALAR NODE 2 PHI VALUE SPEC CONC SCALAR NODE 3 PHI VALUE SPEC CONC SCALAR NODE 4 PHI VALUE SPEC CONC SCALAR NODE 5 PHI VALUE 6 V. N. Kaliakin

7 SPEC CONC SCALAR NODE 10 PHI VALUE SPEC CONC SCALAR NODE 18 PHI VALUE SPEC CONC SCALAR NODE 24 PHI VALUE SPEC CONC SCALAR NODE 28 PHI VALUE SPEC CONC SCALAR NODE 30 PHI VALUE! FINISHED DATA! SOLUTION TIME FINAL 1.0 INCR 1 OUTPUT 1:10:1! FINISHED LOADING! A plot of the mesh, obtained using Tecplot, is shown below. The results obtained analyzing the above data by the ud_scalar computer program are given on the following page (NOTE: some lines of input have been wrapped by the word processing software). 7 V. N. Kaliakin

8 _/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/ _/ _/ _/ ud_scalar: version 1.00: _/ _/ _/ _/ _/ _/ DATE OF ANALYSIS :: day:16 month:03 year:10 _/ _/ _/ _/ ANALYSIS INITIATED AT TIME :: 15:11:38 _/ _/ _/ _/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/ data file translated from PATRAN neutral file using PAT2UD_SCALAR version 1.0 P3/PATRAN Neutral File from: C:\filled_slot.db neutral file creation date: 10-Mar-08 neutral file creation time: 14:13:19 PATRAN release number: 3.0 torsion of a 'slot' cross-section with TWO different materials 80 element mesh of QS4 elements ====================================================================== D Y N A M I C S T O R A G E A L L O C A T I O N ====================================================================== Largest NODE number which can used in the mesh = 99 Max. no. of CONSTANT scalar conductivity idealizations = 2 Max. no. of 4-node quad. "scalar" (QS4) elements = 80 ====================================================================== = G E N E R A L A N A L Y S I S I N F O R M A T I O N = ====================================================================== --> analysis with SCALAR primary dependent variables shall be performed --> TWO-DIMENSIONAL solution domain assumed (PLANE STRAIN idealization) --> solver type used : SKYLINE --> storage type : SYMMETRIC --> "Isoparametric" mesh generation scheme used --> LINEAR analysis ====================================================================== = S C A L A R M A T E R I A L P A R A M E T E R S = 8 V. N. Kaliakin

9 ====================================================================== --> idealization no.: 1 ~~~~~~~~~~~~~~~ type : constant scalar conductivity coefficients info. : this is material #1 "Conductivities" (material parameters) : k_11 = 1.000E+00 k_12 = 0.000E+00 k_13 = 0.000E+00 k_22 = 1.000E+00 k_23 = 0.000E+00 k_33 = 1.000E+00 source term S_1 = 9.600E+03 source term S_2 = 0.000E+00 --> idealization no.: 2 ~~~~~~~~~~~~~~~ type : constant scalar conductivity coefficients info. : this is material #2 "Conductivities" (material parameters) : k_11 = 1.000E+00 k_12 = 0.000E+00 k_13 = 0.000E+00 k_22 = 1.000E+00 k_23 = 0.000E+00 k_33 = 1.000E+00 source term S_1 = 4.800E+03 source term S_2 = 0.000E+00 ====================================================================== = N O D A L C O O R D I N A T E S = ====================================================================== 1 x1 = 3.000E+00 x2 = 4.000E+00 2 x1 = 2.500E+00 x2 = 4.000E+00 3 x1 = 2.000E+00 x2 = 4.000E+00 4 x1 = 1.500E+00 x2 = 4.000E+00 5 x1 = 1.000E+00 x2 = 4.000E+00 6 x1 = 7.500E-01 x2 = 4.000E+00 7 x1 = 5.000E-01 x2 = 4.000E+00 8 x1 = 2.500E-01 x2 = 4.000E+00 9 x1 = 0.000E+00 x2 = 4.000E x1 = 3.000E+00 x2 = 3.750E x1 = 2.500E+00 x2 = 3.750E x1 = 2.000E+00 x2 = 3.750E x1 = 1.500E+00 x2 = 3.750E x1 = 1.000E+00 x2 = 3.750E x1 = 7.500E-01 x2 = 3.750E x1 = 5.000E-01 x2 = 3.750E x1 = 2.500E-01 x2 = 3.750E x1 = 0.000E+00 x2 = 3.750E x1 = 3.000E+00 x2 = 3.500E x1 = 2.500E+00 x2 = 3.500E x1 = 2.000E+00 x2 = 3.500E x1 = 1.500E+00 x2 = 3.500E x1 = 1.000E+00 x2 = 3.500E x1 = 7.500E-01 x2 = 3.500E x1 = 5.000E-01 x2 = 3.500E+00 9 V. N. Kaliakin

10 26 x1 = 2.500E-01 x2 = 3.500E x1 = 0.000E+00 x2 = 3.500E x1 = 3.000E+00 x2 = 3.250E x1 = 2.500E+00 x2 = 3.250E x1 = 2.000E+00 x2 = 3.250E x1 = 1.500E+00 x2 = 3.250E x1 = 1.000E+00 x2 = 3.250E x1 = 7.500E-01 x2 = 3.250E x1 = 5.000E-01 x2 = 3.250E x1 = 2.500E-01 x2 = 3.250E x1 = 0.000E+00 x2 = 3.250E x1 = 3.000E+00 x2 = 3.000E x1 = 2.500E+00 x2 = 3.000E x1 = 2.000E+00 x2 = 3.000E x1 = 1.500E+00 x2 = 3.000E x1 = 1.000E+00 x2 = 3.000E x1 = 7.500E-01 x2 = 3.000E x1 = 5.000E-01 x2 = 3.000E x1 = 2.500E-01 x2 = 3.000E x1 = 0.000E+00 x2 = 3.000E x1 = 3.000E+00 x2 = 2.500E x1 = 2.500E+00 x2 = 2.500E x1 = 2.000E+00 x2 = 2.500E x1 = 1.500E+00 x2 = 2.500E x1 = 1.000E+00 x2 = 2.500E x1 = 7.500E-01 x2 = 2.500E x1 = 5.000E-01 x2 = 2.500E x1 = 2.500E-01 x2 = 2.500E x1 = 0.000E+00 x2 = 2.500E x1 = 3.000E+00 x2 = 2.000E x1 = 3.000E+00 x2 = 1.500E x1 = 3.000E+00 x2 = 1.000E x1 = 3.000E+00 x2 = 5.000E x1 = 3.000E+00 x2 = 0.000E x1 = 2.500E+00 x2 = 2.000E x1 = 2.000E+00 x2 = 2.000E x1 = 1.500E+00 x2 = 2.000E x1 = 1.000E+00 x2 = 2.000E x1 = 7.500E-01 x2 = 2.000E x1 = 5.000E-01 x2 = 2.000E x1 = 2.500E-01 x2 = 2.000E x1 = 0.000E+00 x2 = 2.000E x1 = 2.500E+00 x2 = 1.500E x1 = 2.500E+00 x2 = 1.000E x1 = 2.500E+00 x2 = 5.000E x1 = 2.500E+00 x2 = 0.000E x1 = 2.000E+00 x2 = 1.500E x1 = 1.500E+00 x2 = 1.500E x1 = 1.000E+00 x2 = 1.500E x1 = 7.500E-01 x2 = 1.500E x1 = 5.000E-01 x2 = 1.500E x1 = 2.500E-01 x2 = 1.500E x1 = 0.000E+00 x2 = 1.500E x1 = 2.000E+00 x2 = 1.000E x1 = 2.000E+00 x2 = 5.000E x1 = 2.000E+00 x2 = 0.000E x1 = 0.000E+00 x2 = 1.000E x1 = 0.000E+00 x2 = 5.000E x1 = 0.000E+00 x2 = 0.000E x1 = 1.500E+00 x2 = 1.000E x1 = 1.000E+00 x2 = 1.000E x1 = 7.500E-01 x2 = 1.000E x1 = 5.000E-01 x2 = 1.000E V. N. Kaliakin

11 89 x1 = 2.500E-01 x2 = 1.000E x1 = 1.500E+00 x2 = 5.000E x1 = 1.500E+00 x2 = 0.000E x1 = 2.500E-01 x2 = 5.000E x1 = 2.500E-01 x2 = 0.000E x1 = 1.000E+00 x2 = 5.000E x1 = 7.500E-01 x2 = 5.000E x1 = 5.000E-01 x2 = 5.000E x1 = 1.000E+00 x2 = 0.000E x1 = 5.000E-01 x2 = 0.000E x1 = 7.500E-01 x2 = 0.000E+00 ====================================================================== = E L E M E N T I N F O R M A T I O N = ====================================================================== --> number : 1 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 2 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 3 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 4 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 5 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 6 (type : QS4 ) (kind : SCALAR ) 11 V. N. Kaliakin

12 ~~~~~~ nodes : > number : 7 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 8 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 9 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 10 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 11 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 12 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 13 (type : QS4 ) (kind : SCALAR ) 12 V. N. Kaliakin

13 ~~~~~~ nodes : > number : 14 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 15 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 16 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 17 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 18 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 19 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 20 (type : QS4 ) (kind : SCALAR ) 13 V. N. Kaliakin

14 ~~~~~~ nodes : > number : 21 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 22 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 23 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 24 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 25 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 26 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 27 (type : QS4 ) (kind : SCALAR ) 14 V. N. Kaliakin

15 ~~~~~~ nodes : > number : 28 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 29 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 30 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 31 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 32 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 33 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 34 (type : QS4 ) (kind : SCALAR ) 15 V. N. Kaliakin

16 ~~~~~~ nodes : > number : 35 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 36 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 37 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 38 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 39 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 40 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 41 (type : QS4 ) (kind : SCALAR ) 16 V. N. Kaliakin

17 ~~~~~~ nodes : > number : 42 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 43 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 44 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 45 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 46 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 47 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 48 (type : QS4 ) (kind : SCALAR ) 17 V. N. Kaliakin

18 ~~~~~~ nodes : > number : 49 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 50 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 51 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 52 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 53 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 54 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 55 (type : QS4 ) (kind : SCALAR ) 18 V. N. Kaliakin

19 ~~~~~~ nodes : > number : 56 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 57 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 58 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 59 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 60 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 61 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 62 (type : QS4 ) (kind : SCALAR ) 19 V. N. Kaliakin

20 ~~~~~~ nodes : > number : 63 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 64 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 65 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 66 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 67 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 68 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 69 (type : QS4 ) (kind : SCALAR ) 20 V. N. Kaliakin

21 ~~~~~~ nodes : > number : 70 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 71 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 72 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 73 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 74 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 75 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 76 (type : QS4 ) (kind : SCALAR ) 21 V. N. Kaliakin

22 ~~~~~~ nodes : > number : 77 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 78 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 79 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : > number : 80 (type : QS4 ) (kind : SCALAR ) ~~~~~~ nodes : ====================================================================== = N O D E P O I N T S P E C I F I C A T I O N S = ====================================================================== Node ( c o o r d i n a t e s ) Number s p e c i f i c a t i o n: ~~~~~~ ~~~~~~~~~~~~~~~~~~~~~~~~~ 1 : ( x1 = 3.000E+00, x2 = 4.000E+00 ) 2 : ( x1 = 2.500E+00, x2 = 4.000E+00 ) 3 : ( x1 = 2.000E+00, x2 = 4.000E+00 ) 22 V. N. Kaliakin

23 4 : ( x1 = 1.500E+00, x2 = 4.000E+00 ) 5 : ( x1 = 1.000E+00, x2 = 4.000E+00 ) 6 : ( x1 = 7.500E-01, x2 = 4.000E+00 ) 7 : ( x1 = 5.000E-01, x2 = 4.000E+00 ) 8 : ( x1 = 2.500E-01, x2 = 4.000E+00 ) 9 : ( x1 = 0.000E+00, x2 = 4.000E+00 ) 10 : ( x1 = 3.000E+00, x2 = 3.750E+00 ) 19 : ( x1 = 3.000E+00, x2 = 3.500E+00 ) 28 : ( x1 = 3.000E+00, x2 = 3.250E+00 ) 37 : ( x1 = 3.000E+00, x2 = 3.000E+00 ) 46 : ( x1 = 3.000E+00, x2 = 2.500E+00 ) 55 : ( x1 = 3.000E+00, x2 = 2.000E+00 ) 56 : ( x1 = 3.000E+00, x2 = 1.500E+00 ) 57 : ( x1 = 3.000E+00, x2 = 1.000E+00 ) 58 : ( x1 = 3.000E+00, x2 = 5.000E-01 ) 59 : ( x1 = 3.000E+00, x2 = 0.000E+00 ) ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ end of mathematical model data ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ====================================================================== = E L E M E N T F L U X E S = ====================================================================== --> element 1 ( type = QS4 ) : [ 23 V. N. Kaliakin

24 @(x1 = 1.250E-01, x2 = 2.500E-01 ) : flux_1 = 3.369E+02 ; flux_2 = 5.229E+02 --> element 2 ( type = QS4 ) : = 1.250E-01, x2 = 7.500E-01 ) : flux_1 = 3.212E+02 ; flux_2 = 1.600E+03 --> element 3 ( type = QS4 ) : = 1.250E-01, x2 = 1.250E+00 ) : flux_1 = 2.904E+02 ; flux_2 = 2.771E+03 --> element 4 ( type = QS4 ) : = 1.250E-01, x2 = 1.750E+00 ) : flux_1 = 2.476E+02 ; flux_2 = 4.091E+03 --> element 5 ( type = QS4 ) : = 1.250E-01, x2 = 2.250E+00 ) : flux_1 = 2.008E+02 ; flux_2 = 5.595E+03 --> element 6 ( type = QS4 ) : = 1.250E-01, x2 = 2.750E+00 ) : flux_1 = 1.572E+02 ; flux_2 = 7.281E+03 --> element 7 ( type = QS4 ) : = 3.750E-01, x2 = 2.500E-01 ) : flux_1 = 1.015E+03 ; flux_2 = 5.149E+02 --> element 8 ( type = QS4 ) : = 3.750E-01, x2 = 7.500E-01 ) : flux_1 = 9.673E+02 ; flux_2 = 1.577E+03 --> element 9 ( type = QS4 ) : = 3.750E-01, x2 = 1.250E+00 ) : flux_1 = 8.744E+02 ; flux_2 = 2.733E+03 --> element 10 ( type = QS4 ) : = 3.750E-01, x2 = 1.750E+00 ) : flux_1 = 7.441E+02 ; flux_2 = 4.043E+03 --> element 11 ( type = QS4 ) : = 3.750E-01, x2 = 2.250E+00 ) : flux_1 = 6.006E+02 ; flux_2 = 5.548E V. N. Kaliakin

25 --> element 12 ( type = QS4 ) : = 3.750E-01, x2 = 2.750E+00 ) : flux_1 = 4.719E+02 ; flux_2 = 7.241E+03 --> element 13 ( type = QS4 ) : = 6.250E-01, x2 = 2.500E-01 ) : flux_1 = 1.704E+03 ; flux_2 = 4.990E+02 --> element 14 ( type = QS4 ) : = 6.250E-01, x2 = 7.500E-01 ) : flux_1 = 1.625E+03 ; flux_2 = 1.529E+03 --> element 15 ( type = QS4 ) : = 6.250E-01, x2 = 1.250E+00 ) : flux_1 = 1.469E+03 ; flux_2 = 2.656E+03 --> element 16 ( type = QS4 ) : = 6.250E-01, x2 = 1.750E+00 ) : flux_1 = 1.245E+03 ; flux_2 = 3.943E+03 --> element 17 ( type = QS4 ) : = 6.250E-01, x2 = 2.250E+00 ) : flux_1 = 9.930E+02 ; flux_2 = 5.450E+03 --> element 18 ( type = QS4 ) : = 6.250E-01, x2 = 2.750E+00 ) : flux_1 = 7.848E+02 ; flux_2 = 7.171E+03 --> element 19 ( type = QS4 ) : = 8.750E-01, x2 = 2.500E-01 ) : flux_1 = 2.414E+03 ; flux_2 = 4.751E+02 --> element 20 ( type = QS4 ) : = 8.750E-01, x2 = 7.500E-01 ) : flux_1 = 2.302E+03 ; flux_2 = 1.458E+03 --> element 21 ( type = QS4 ) : = 8.750E-01, x2 = 1.250E+00 ) : flux_1 = 2.081E+03 ; flux_2 = 2.538E+03 --> element 22 ( type = QS4 ) : [ 25 V. N. Kaliakin

26 @(x1 = 8.750E-01, x2 = 1.750E+00 ) : flux_1 = 1.758E+03 ; flux_2 = 3.788E+03 --> element 23 ( type = QS4 ) : = 8.750E-01, x2 = 2.250E+00 ) : flux_1 = 1.369E+03 ; flux_2 = 5.284E+03 --> element 24 ( type = QS4 ) : = 8.750E-01, x2 = 2.750E+00 ) : flux_1 = 1.078E+03 ; flux_2 = 7.087E+03 --> element 25 ( type = QS4 ) : = 1.250E+00, x2 = 2.500E-01 ) : flux_1 = 4.740E+03 ; flux_2 = 4.215E+02 --> element 26 ( type = QS4 ) : = 1.250E+00, x2 = 7.500E-01 ) : flux_1 = 4.581E+03 ; flux_2 = 1.296E+03 --> element 27 ( type = QS4 ) : = 1.250E+00, x2 = 1.250E+00 ) : flux_1 = 4.259E+03 ; flux_2 = 2.267E+03 --> element 28 ( type = QS4 ) : = 1.250E+00, x2 = 1.750E+00 ) : flux_1 = 3.774E+03 ; flux_2 = 3.414E+03 --> element 29 ( type = QS4 ) : = 1.250E+00, x2 = 2.250E+00 ) : flux_1 = 3.149E+03 ; flux_2 = 4.838E+03 --> element 30 ( type = QS4 ) : = 1.250E+00, x2 = 2.750E+00 ) : flux_1 = 2.479E+03 ; flux_2 = 6.717E+03 --> element 31 ( type = QS4 ) : = 1.750E+00, x2 = 2.500E-01 ) : flux_1 = 8.781E+03 ; flux_2 = 3.278E+02 --> element 32 ( type = QS4 ) : = 1.750E+00, x2 = 7.500E-01 ) : flux_1 = 8.562E+03 ; flux_2 = 1.011E V. N. Kaliakin

27 --> element 33 ( type = QS4 ) : = 1.750E+00, x2 = 1.250E+00 ) : flux_1 = 8.112E+03 ; flux_2 = 1.780E+03 --> element 34 ( type = QS4 ) : = 1.750E+00, x2 = 1.750E+00 ) : flux_1 = 7.407E+03 ; flux_2 = 2.712E+03 --> element 35 ( type = QS4 ) : = 1.750E+00, x2 = 2.250E+00 ) : flux_1 = 6.418E+03 ; flux_2 = 3.925E+03 --> element 36 ( type = QS4 ) : = 1.750E+00, x2 = 2.750E+00 ) : flux_1 = 5.088E+03 ; flux_2 = 5.630E+03 --> element 37 ( type = QS4 ) : = 2.250E+00, x2 = 2.500E-01 ) : flux_1 = 1.304E+04 ; flux_2 = 2.085E+02 --> element 38 ( type = QS4 ) : = 2.250E+00, x2 = 7.500E-01 ) : flux_1 = 1.277E+04 ; flux_2 = 6.442E+02 --> element 39 ( type = QS4 ) : = 2.250E+00, x2 = 1.250E+00 ) : flux_1 = 1.222E+04 ; flux_2 = 1.140E+03 --> element 40 ( type = QS4 ) : = 2.250E+00, x2 = 1.750E+00 ) : flux_1 = 1.132E+04 ; flux_2 = 1.754E+03 --> element 41 ( type = QS4 ) : = 2.250E+00, x2 = 2.250E+00 ) : flux_1 = 1.001E+04 ; flux_2 = 2.577E+03 --> element 42 ( type = QS4 ) : = 2.250E+00, x2 = 2.750E+00 ) : flux_1 = 8.118E+03 ; flux_2 = 3.761E+03 --> element 43 ( type = QS4 ) : [ 27 V. N. Kaliakin

28 @(x1 = 2.750E+00, x2 = 2.500E-01 ) : flux_1 = 1.755E+04 ; flux_2 = 7.157E+01 --> element 44 ( type = QS4 ) : = 2.750E+00, x2 = 7.500E-01 ) : flux_1 = 1.726E+04 ; flux_2 = 2.215E+02 --> element 45 ( type = QS4 ) : = 2.750E+00, x2 = 1.250E+00 ) : flux_1 = 1.665E+04 ; flux_2 = 3.932E+02 --> element 46 ( type = QS4 ) : = 2.750E+00, x2 = 1.750E+00 ) : flux_1 = 1.564E+04 ; flux_2 = 6.080E+02 --> element 47 ( type = QS4 ) : = 2.750E+00, x2 = 2.250E+00 ) : flux_1 = 1.414E+04 ; flux_2 = 8.991E+02 --> element 48 ( type = QS4 ) : = 2.750E+00, x2 = 2.750E+00 ) : flux_1 = 1.191E+04 ; flux_2 = 1.326E+03 --> element 49 ( type = QS4 ) : = 2.750E+00, x2 = 3.125E+00 ) : flux_1 = 9.694E+03 ; flux_2 = 1.786E+03 --> element 50 ( type = QS4 ) : = 2.750E+00, x2 = 3.375E+00 ) : flux_1 = 7.693E+03 ; flux_2 = 2.218E+03 --> element 51 ( type = QS4 ) : = 2.750E+00, x2 = 3.625E+00 ) : flux_1 = 5.172E+03 ; flux_2 = 2.824E+03 --> element 52 ( type = QS4 ) : = 2.750E+00, x2 = 3.875E+00 ) : flux_1 = 1.880E+03 ; flux_2 = 3.760E+03 --> element 53 ( type = QS4 ) : = 2.250E+00, x2 = 3.125E+00 ) : flux_1 = 6.293E+03 ; flux_2 = 5.001E V. N. Kaliakin

29 --> element 54 ( type = QS4 ) : = 2.250E+00, x2 = 3.375E+00 ) : flux_1 = 4.738E+03 ; flux_2 = 6.115E+03 --> element 55 ( type = QS4 ) : = 2.250E+00, x2 = 3.625E+00 ) : flux_1 = 2.941E+03 ; flux_2 = 7.562E+03 --> element 56 ( type = QS4 ) : = 2.250E+00, x2 = 3.875E+00 ) : flux_1 = 9.919E+02 ; flux_2 = 9.504E+03 --> element 57 ( type = QS4 ) : = 1.750E+00, x2 = 3.125E+00 ) : flux_1 = 3.858E+03 ; flux_2 = 7.370E+03 --> element 58 ( type = QS4 ) : = 1.750E+00, x2 = 3.375E+00 ) : flux_1 = 2.865E+03 ; flux_2 = 8.843E+03 --> element 59 ( type = QS4 ) : = 1.750E+00, x2 = 3.625E+00 ) : flux_1 = 1.771E+03 ; flux_2 = 1.062E+04 --> element 60 ( type = QS4 ) : = 1.750E+00, x2 = 3.875E+00 ) : flux_1 = 6.006E+02 ; flux_2 = 1.269E+04 --> element 61 ( type = QS4 ) : = 1.250E+00, x2 = 3.125E+00 ) : flux_1 = 1.970E+03 ; flux_2 = 8.673E+03 --> element 62 ( type = QS4 ) : = 1.250E+00, x2 = 3.375E+00 ) : flux_1 = 1.531E+03 ; flux_2 = 1.041E+04 --> element 63 ( type = QS4 ) : = 1.250E+00, x2 = 3.625E+00 ) : flux_1 = 9.665E+02 ; flux_2 = 1.237E+04 --> element 64 ( type = QS4 ) : [ 29 V. N. Kaliakin

30 @(x1 = 1.250E+00, x2 = 3.875E+00 ) : flux_1 = 3.298E+02 ; flux_2 = 1.455E+04 --> element 65 ( type = QS4 ) : = 8.750E-01, x2 = 3.125E+00 ) : flux_1 = 9.761E+02 ; flux_2 = 9.054E+03 --> element 66 ( type = QS4 ) : = 8.750E-01, x2 = 3.375E+00 ) : flux_1 = 8.358E+02 ; flux_2 = 1.104E+04 --> element 67 ( type = QS4 ) : = 8.750E-01, x2 = 3.625E+00 ) : flux_1 = 5.443E+02 ; flux_2 = 1.315E+04 --> element 68 ( type = QS4 ) : = 8.750E-01, x2 = 3.875E+00 ) : flux_1 = 1.882E+02 ; flux_2 = 1.540E+04 --> element 69 ( type = QS4 ) : = 6.250E-01, x2 = 3.125E+00 ) : flux_1 = 6.530E+02 ; flux_2 = 9.120E+03 --> element 70 ( type = QS4 ) : = 6.250E-01, x2 = 3.375E+00 ) : flux_1 = 5.282E+02 ; flux_2 = 1.124E+04 --> element 71 ( type = QS4 ) : = 6.250E-01, x2 = 3.625E+00 ) : flux_1 = 3.467E+02 ; flux_2 = 1.342E+04 --> element 72 ( type = QS4 ) : = 6.250E-01, x2 = 3.875E+00 ) : flux_1 = 1.204E+02 ; flux_2 = 1.571E+04 --> element 73 ( type = QS4 ) : = 3.750E-01, x2 = 3.125E+00 ) : flux_1 = 3.790E+02 ; flux_2 = 9.204E+03 --> element 74 ( type = QS4 ) : = 3.750E-01, x2 = 3.375E+00 ) : flux_1 = 2.978E+02 ; flux_2 = 1.137E V. N. Kaliakin

31 --> element 75 ( type = QS4 ) : = 3.750E-01, x2 = 3.625E+00 ) : flux_1 = 1.925E+02 ; flux_2 = 1.359E+04 --> element 76 ( type = QS4 ) : = 3.750E-01, x2 = 3.875E+00 ) : flux_1 = 6.679E+01 ; flux_2 = 1.589E+04 --> element 77 ( type = QS4 ) : = 1.250E-01, x2 = 3.125E+00 ) : flux_1 = 1.244E+02 ; flux_2 = 9.251E+03 --> element 78 ( type = QS4 ) : = 1.250E-01, x2 = 3.375E+00 ) : flux_1 = 9.659E+01 ; flux_2 = 1.143E+04 --> element 79 ( type = QS4 ) : = 1.250E-01, x2 = 3.625E+00 ) : flux_1 = 6.189E+01 ; flux_2 = 1.367E+04 --> element 80 ( type = QS4 ) : = 1.250E-01, x2 = 3.875E+00 ) : flux_1 = 2.139E+01 ; flux_2 = 1.598E+04 maximum values of element variables : max flux_1 = x1 = 2.750E+00, x2 = 2.500E-01 max flux_2 = x1 = 1.250E-01, x2 = 3.875E+00 ====================================================================== = N O D A L Q U A N T I T I E S = ====================================================================== 1 ( x1 = 3.000E+00, x2 = 4.000E+00 ), phi = 9.401E-18 2 ( x1 = 2.500E+00, x2 = 4.000E+00 ), phi = 4.068E-17 3 ( x1 = 2.000E+00, x2 = 4.000E+00 ), phi = 6.697E-17 4 ( x1 = 1.500E+00, x2 = 4.000E+00 ), phi = 8.199E-17 5 ( x1 = 1.000E+00, x2 = 4.000E+00 ), phi = 6.717E-17 6 ( x1 = 7.500E-01, x2 = 4.000E+00 ), phi = 5.828E-17 7 ( x1 = 5.000E-01, x2 = 4.000E+00 ), phi = 5.922E-17 8 ( x1 = 2.500E-01, x2 = 4.000E+00 ), phi = 5.974E-17 9 ( x1 = 0.000E+00, x2 = 4.000E+00 ), phi = 2.995E ( x1 = 3.000E+00, x2 = 3.750E+00 ), phi = 8.939E ( x1 = 2.500E+00, x2 = 3.750E+00 ), phi = 1.880E ( x1 = 2.000E+00, x2 = 3.750E+00 ), phi = 2.872E ( x1 = 1.500E+00, x2 = 3.750E+00 ), phi = 3.473E V. N. Kaliakin

32 14 ( x1 = 1.000E+00, x2 = 3.750E+00 ), phi = 3.802E ( x1 = 7.500E-01, x2 = 3.750E+00 ), phi = 3.897E ( x1 = 5.000E-01, x2 = 3.750E+00 ), phi = 3.957E ( x1 = 2.500E-01, x2 = 3.750E+00 ), phi = 3.990E ( x1 = 0.000E+00, x2 = 3.750E+00 ), phi = 4.001E ( x1 = 3.000E+00, x2 = 3.500E+00 ), phi = 1.824E ( x1 = 2.500E+00, x2 = 3.500E+00 ), phi = 3.292E ( x1 = 2.000E+00, x2 = 3.500E+00 ), phi = 5.241E ( x1 = 1.500E+00, x2 = 3.500E+00 ), phi = 6.411E ( x1 = 1.000E+00, x2 = 3.500E+00 ), phi = 7.048E ( x1 = 7.500E-01, x2 = 3.500E+00 ), phi = 7.226E ( x1 = 5.000E-01, x2 = 3.500E+00 ), phi = 7.339E ( x1 = 2.500E-01, x2 = 3.500E+00 ), phi = 7.402E ( x1 = 0.000E+00, x2 = 3.500E+00 ), phi = 7.422E ( x1 = 3.000E+00, x2 = 3.250E+00 ), phi = 2.532E ( x1 = 2.500E+00, x2 = 3.250E+00 ), phi = 4.401E ( x1 = 2.000E+00, x2 = 3.250E+00 ), phi = 7.190E ( x1 = 1.500E+00, x2 = 3.250E+00 ), phi = 8.884E ( x1 = 1.000E+00, x2 = 3.250E+00 ), phi = 9.778E ( x1 = 7.500E-01, x2 = 3.250E+00 ), phi = 1.002E ( x1 = 5.000E-01, x2 = 3.250E+00 ), phi = 1.017E ( x1 = 2.500E-01, x2 = 3.250E+00 ), phi = 1.026E ( x1 = 0.000E+00, x2 = 3.250E+00 ), phi = 1.028E ( x1 = 3.000E+00, x2 = 3.000E+00 ), phi = 6.060E ( x1 = 2.500E+00, x2 = 3.000E+00 ), phi = 5.294E ( x1 = 2.000E+00, x2 = 3.000E+00 ), phi = 8.798E ( x1 = 1.500E+00, x2 = 3.000E+00 ), phi = 1.096E ( x1 = 1.000E+00, x2 = 3.000E+00 ), phi = 1.204E ( x1 = 7.500E-01, x2 = 3.000E+00 ), phi = 1.229E ( x1 = 5.000E-01, x2 = 3.000E+00 ), phi = 1.246E ( x1 = 2.500E-01, x2 = 3.000E+00 ), phi = 1.256E ( x1 = 0.000E+00, x2 = 3.000E+00 ), phi = 1.260E ( x1 = 3.000E+00, x2 = 2.500E+00 ), phi = 9.715E ( x1 = 2.500E+00, x2 = 2.500E+00 ), phi = 6.619E ( x1 = 2.000E+00, x2 = 2.500E+00 ), phi = 1.123E ( x1 = 1.500E+00, x2 = 2.500E+00 ), phi = 1.416E ( x1 = 1.000E+00, x2 = 2.500E+00 ), phi = 1.556E ( x1 = 7.500E-01, x2 = 2.500E+00 ), phi = 1.585E ( x1 = 5.000E-01, x2 = 2.500E+00 ), phi = 1.607E ( x1 = 2.500E-01, x2 = 2.500E+00 ), phi = 1.620E ( x1 = 0.000E+00, x2 = 2.500E+00 ), phi = 1.624E ( x1 = 3.000E+00, x2 = 2.000E+00 ), phi = 1.113E ( x1 = 3.000E+00, x2 = 1.500E+00 ), phi = 1.208E ( x1 = 3.000E+00, x2 = 1.000E+00 ), phi = 1.269E ( x1 = 3.000E+00, x2 = 5.000E-01 ), phi = 1.304E ( x1 = 3.000E+00, x2 = 0.000E+00 ), phi = 6.574E ( x1 = 2.500E+00, x2 = 2.000E+00 ), phi = 7.518E ( x1 = 2.000E+00, x2 = 2.000E+00 ), phi = 1.291E ( x1 = 1.500E+00, x2 = 2.000E+00 ), phi = 1.640E ( x1 = 1.000E+00, x2 = 2.000E+00 ), phi = 1.815E ( x1 = 7.500E-01, x2 = 2.000E+00 ), phi = 1.854E ( x1 = 5.000E-01, x2 = 2.000E+00 ), phi = 1.882E ( x1 = 2.500E-01, x2 = 2.000E+00 ), phi = 1.899E ( x1 = 0.000E+00, x2 = 2.000E+00 ), phi = 1.905E ( x1 = 2.500E+00, x2 = 1.500E+00 ), phi = 8.126E ( x1 = 2.500E+00, x2 = 1.000E+00 ), phi = 8.519E ( x1 = 2.500E+00, x2 = 5.000E-01 ), phi = 8.741E ( x1 = 2.500E+00, x2 = 0.000E+00 ), phi = 8.812E ( x1 = 2.000E+00, x2 = 1.500E+00 ), phi = 1.406E ( x1 = 1.500E+00, x2 = 1.500E+00 ), phi = 1.797E ( x1 = 1.000E+00, x2 = 1.500E+00 ), phi = 2.000E ( x1 = 7.500E-01, x2 = 1.500E+00 ), phi = 2.048E ( x1 = 5.000E-01, x2 = 1.500E+00 ), phi = 2.083E V. N. Kaliakin

33 77 ( x1 = 2.500E-01, x2 = 1.500E+00 ), phi = 2.103E ( x1 = 0.000E+00, x2 = 1.500E+00 ), phi = 2.110E ( x1 = 2.000E+00, x2 = 1.000E+00 ), phi = 1.480E ( x1 = 2.000E+00, x2 = 5.000E-01 ), phi = 1.523E ( x1 = 2.000E+00, x2 = 0.000E+00 ), phi = 1.536E ( x1 = 0.000E+00, x2 = 1.000E+00 ), phi = 2.249E ( x1 = 0.000E+00, x2 = 5.000E-01 ), phi = 2.329E ( x1 = 0.000E+00, x2 = 0.000E+00 ), phi = 2.355E ( x1 = 1.500E+00, x2 = 1.000E+00 ), phi = 1.900E ( x1 = 1.000E+00, x2 = 1.000E+00 ), phi = 2.123E ( x1 = 7.500E-01, x2 = 1.000E+00 ), phi = 2.179E ( x1 = 5.000E-01, x2 = 1.000E+00 ), phi = 2.218E ( x1 = 2.500E-01, x2 = 1.000E+00 ), phi = 2.241E ( x1 = 1.500E+00, x2 = 5.000E-01 ), phi = 1.959E ( x1 = 1.500E+00, x2 = 0.000E+00 ), phi = 1.978E ( x1 = 2.500E-01, x2 = 5.000E-01 ), phi = 2.321E ( x1 = 2.500E-01, x2 = 0.000E+00 ), phi = 2.347E ( x1 = 1.000E+00, x2 = 5.000E-01 ), phi = 2.194E ( x1 = 7.500E-01, x2 = 5.000E-01 ), phi = 2.254E ( x1 = 5.000E-01, x2 = 5.000E-01 ), phi = 2.296E ( x1 = 1.000E+00, x2 = 0.000E+00 ), phi = 2.217E ( x1 = 5.000E-01, x2 = 0.000E+00 ), phi = 2.321E ( x1 = 7.500E-01, x2 = 0.000E+00 ), phi = 2.278E+04 max phi = node 84 ( 0.000E+00, 0.000E+00) ud_scalar -> end of analysis A contour plot of the primary dependent variable, drawn using Tecplot, is provided on the following page. 33 V. N. Kaliakin

34 34 V. N. Kaliakin