|
21 | 21 |
|
22 | 22 | # %% |
23 | 23 | # Define mesh sizes |
24 | | -mesh_size_list = [50, 20, 10, 5, 3, 2, 1] |
25 | | -nr_list = [4, 8, 12, 16, 20, 24, 32, 64] |
| 24 | +mesh_size_list = [200, 100, 50, 20, 10, 5] |
| 25 | +nr_list = [4, 8, 12, 16, 20, 24, 32] |
26 | 26 |
|
27 | 27 | # %% |
28 | 28 | # Initialise result lists |
|
33 | 33 |
|
34 | 34 | # %% |
35 | 35 | # Calculate reference solution |
36 | | -geometry = steel_sections.i_section(d=203, b=133, t_f=7.8, t_w=5.8, r=8.9, n_r=64) |
37 | | -geometry.create_mesh(mesh_sizes=[0.5]) # create mesh |
| 36 | +geometry = steel_sections.i_section(d=203, b=133, t_f=7.8, t_w=5.8, r=8.9, n_r=32) |
| 37 | +geometry.create_mesh(mesh_sizes=[5]) # create mesh |
38 | 38 | section = Section(geometry) # create a Section object |
39 | 39 | section.calculate_geometric_properties() |
40 | 40 | section.calculate_warping_properties() |
41 | 41 | j_reference = section.get_j() # get the torsion constant |
42 | 42 |
|
43 | 43 | # %% |
44 | | -# Run through mesh_sizes with n_r = 16 |
| 44 | +# Run through mesh_sizes with n_r = 8 |
45 | 45 | for mesh_size in mesh_size_list: |
46 | | - geometry = steel_sections.i_section(d=203, b=133, t_f=7.8, t_w=5.8, r=8.9, n_r=16) |
| 46 | + geometry = steel_sections.i_section(d=203, b=133, t_f=7.8, t_w=5.8, r=8.9, n_r=8) |
47 | 47 | geometry.create_mesh(mesh_sizes=[mesh_size]) # create mesh |
48 | 48 | section = Section(geometry) # create a Section object |
49 | 49 | section.calculate_geometric_properties() |
|
53 | 53 | mesh_results.append(section.get_j()) |
54 | 54 |
|
55 | 55 | # %% |
56 | | -# Run through n_r with mesh_size = 3 |
| 56 | +# Run through n_r with mesh_size = 10 |
57 | 57 | for n_r in nr_list: |
58 | 58 | geometry = steel_sections.i_section(d=203, b=133, t_f=7.8, t_w=5.8, r=8.9, n_r=n_r) |
59 | | - geometry.create_mesh(mesh_sizes=[3]) # create mesh |
| 59 | + geometry.create_mesh(mesh_sizes=[10]) # create mesh |
60 | 60 | section = Section(geometry) # create a Section object |
61 | 61 | section.calculate_geometric_properties() |
62 | 62 | section.calculate_warping_properties() |
|
0 commit comments