structuraid_core 0.1.3 → 0.2.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: e6a3c25bdfad6f0af2c266ad01cf0321e74780e4366a9a0b941e37c108f6d152
-  data.tar.gz: 60afe6f91b70f5acd7ac3d12b1e0c33282afed9620f2c55ee2c29ada6fe3295f
+  metadata.gz: 39b0943343ece08b91c2884a16e133ad77e663c71c8416db4d65248a20871b67
+  data.tar.gz: e9ead6c46af568fd4f7d7fda9cfdd73996f93cff85afdf602161b01caef091f7
 SHA512:
-  metadata.gz: fc6983c21cb7278bb90e9b2b3d4a68e6c873e1fd7de57205724344385200e4674c05fceb61c02342a4329c5ce1515636c614ec8c55d5dba6952f08cd595a309b
-  data.tar.gz: da91f2acb4102987f40c8d05b314727349782e0ac170859423460e131aa6b5793c4761066c583b0c6836d1db0392f47b574911faf75368fef61ebf4b4485497e
+  metadata.gz: b48678b5b1131cf438cfaf5db1d8c2dac8ba1dc89c2d8ee77f580b7c1386f64c7081fbdd1980c7ea6d8373e93fccb7b414bc4555a1ce4c64461e83e1fd40f780
+  data.tar.gz: d986996199861f01fa6d495bf13b87119d0039d9b23146e545da8e063aab47c8b6ef199a3c232e80a569f0350aec37dcfad856557217830b13a0abdac162ade2

data/.rubocop.yml

@@ -25,12 +25,6 @@ RSpec/MultipleMemoizedHelpers:
 RSpec/NestedGroups:
   Enabled: false
 
-RSpec/ContainExactly:
-  Enabled: false
-
-RSpec/MatchArray:
-  Enabled: false
-
 Lint/MissingSuper:
   Enabled: false
 

data/lib/structuraid_core/design_codes/aci_318_19/rc/footings/bending_reinforcement_ratio.rb

@@ -0,0 +1,65 @@
+module StructuraidCore
+  module DesignCodes
+    module ACI31819
+      module RC
+        module Footings
+          class BendingReinforcementRatio
+            MINIMUM_RATIO = 0.0025
+            CODE_REFERENCE = 'ACI 318-19'.freeze
+
+            include DesignCodes::Utils::CodeRequirement
+            use_schema DesignCodes::Schemas::RC::Footings::BendingReinforcementRatioSchema
+
+            def call
+              [
+                solve_cuadratic_equation_for_steel_reinforcement_ratio,
+                MINIMUM_RATIO
+              ].max
+            end
+
+            private
+
+            ## the target is to solve [(1 - sqrt(1 - 4ac))/2a, (1 + sqrt(1 - 4ac))/2a]
+
+            def solve_cuadratic_equation_for_steel_reinforcement_ratio
+              [
+                equation_solver_option_1.negative? ? Float::INFINITY : equation_solver_option_1,
+                equation_solver_option_2
+              ].min
+            end
+
+            def equation_solver_option_1
+              (1 - equation_component_root) / (2 * equation_component_a)
+            end
+
+            def equation_solver_option_2
+              (1 + equation_component_root) / (2 * equation_component_a)
+            end
+
+            def equation_component_a
+              0.59 * design_steel_yield_strength / design_compression_strength
+            end
+
+            def equation_component_c
+              section_area = width * effective_height # mm**2
+              reduced_steel_strength = capacity_reduction_factor * design_steel_yield_strength # N/(mm**2)
+              flexural_moment / (reduced_steel_strength * section_area * effective_height) # N/(mm**2)
+            end
+
+            def equation_component_root
+              unless 4 * equation_component_a * equation_component_c < 1
+                raise RequirementNotFulfilledError.new(
+                  :flexural_moment,
+                  "Moment #{flexural_moment} is too hight for this element",
+                  CODE_REFERENCE
+                )
+              end
+
+              Math.sqrt(1 - 4 * equation_component_a * equation_component_c)
+            end
+          end
+        end
+      end
+    end
+  end
+end

data/lib/structuraid_core/design_codes/aci_318_19/rc/footings/one_way_shear_capacity.rb

@@ -0,0 +1,29 @@
+module StructuraidCore
+  module DesignCodes
+    module ACI31819
+      module RC
+        module Footings
+          class OneWayShearCapacity
+            CODE_REFERENCE = 'ACI 318-19'.freeze
+
+            include DesignCodes::Utils::CodeRequirement
+            use_schema DesignCodes::Schemas::RC::Footings::OneWayShearCapacitySchema
+
+            def call
+              compute_concrete_shear_capacity
+            end
+
+            private
+
+            def compute_concrete_shear_capacity
+              effective_area = effective_height * width
+              concrete_strength = light_concrete_modification_factor * Math.sqrt(design_compression_strength)
+
+              0.17 * concrete_strength * effective_area
+            end
+          end
+        end
+      end
+    end
+  end
+end

data/lib/structuraid_core/design_codes/aci_318_19/rc/footings/punching_critical_section_perimeter.rb

@@ -0,0 +1,105 @@
+require 'byebug'
+
+module StructuraidCore
+  module DesignCodes
+    module ACI31819
+      module RC
+        module Footings
+          class PunchingCriticalSectionPerimeter
+            include DesignCodes::Utils::CodeRequirement
+            use_schema DesignCodes::Schemas::RC::Footings::PunchingCriticalSectionPerimeterSchema
+
+            EDGES_INDEXES = {
+              top: { from: :top_right, to: :top_left },
+              left: { from: :top_left, to: :bottom_left },
+              bottom: { from: :bottom_left, to: :bottom_right },
+              right: { from: :bottom_right, to: :top_right }
+            }.freeze
+
+            def call
+              column_location = footing.find_or_add_column_location(column_absolute_location, column_label)
+              perimeter_vertices = build_initial_perimeter(column_location)
+
+              edges_indexes = EDGES_INDEXES
+              edges_indexes = select_edges_indexes_into_the_footing(edges_indexes, perimeter_vertices)
+              edges_indexes = update_edges_indexes_vertices(edges_indexes, perimeter_vertices)
+              compute_perimeter(edges_indexes, perimeter_vertices)
+            end
+
+            private
+
+            def update_edges_indexes_vertices(edges_indexes, perimeter_vertices)
+              edges_indexes.each do |_edge_name, edge|
+                update_location_to_limit(perimeter_vertices[edge[:from]])
+                update_location_to_limit(perimeter_vertices[edge[:to]])
+              end
+            end
+
+            def select_edges_indexes_into_the_footing(edges_indexes, perimeter_vertices)
+              edges_indexes.select do |_edge_name, edge|
+                start_vertex_inside_footing = footing.inside_me?(perimeter_vertices[edge[:from]])
+                end_vertex_inside_footing = footing.inside_me?(perimeter_vertices[edge[:to]])
+
+                start_vertex_inside_footing || end_vertex_inside_footing
+              end
+            end
+
+            def build_initial_perimeter(column_location)
+              perimeter_vertices_relative_to_column_location.map do |label, perimeter_vertex|
+                location = local_coordinates_system.find_or_add_location_from_vector(
+                  perimeter_vertex + column_location.to_vector,
+                  label: "column_#{column_label}_punching_#{label}"
+                )
+                [label, location]
+              end.to_h
+            end
+
+            def update_location_to_limit(location)
+              location.value_1 = update_to_constraint(location.value_1, 0.5 * footing.length_1)
+              location.value_2 = update_to_constraint(location.value_2, 0.5 * footing.length_2)
+            end
+
+            def update_to_constraint(initial_value, constraint_value)
+              if initial_value > constraint_value
+                constraint_value
+              elsif initial_value < -constraint_value
+                -constraint_value
+              else
+                initial_value
+              end
+            end
+
+            def perimeter_vertices_relative_to_column_location
+              {
+                top_right: perimeter_vertices_relative_to_column_location_with(way_1: 1, way_2: 1),
+                top_left: perimeter_vertices_relative_to_column_location_with(way_1: -1, way_2: 1),
+                bottom_left: perimeter_vertices_relative_to_column_location_with(way_1: -1, way_2: -1),
+                bottom_right: perimeter_vertices_relative_to_column_location_with(way_1: 1, way_2: -1)
+              }
+            end
+
+            def perimeter_vertices_relative_to_column_location_with(way_1:, way_2:)
+              Vector[
+                way_1 * 0.5 * (column_section_length_1 + footing.effective_height),
+                way_2 * 0.5 * (column_section_length_2 + footing.effective_height),
+                0.0
+              ]
+            end
+
+            def local_coordinates_system
+              footing.coordinates_system
+            end
+
+            def compute_perimeter(edges_indexes, perimeter_vertices)
+              edges_indexes.map do |_edge_name, edge|
+                vector_from = perimeter_vertices[edge[:from]].to_vector
+                vector_to = perimeter_vertices[edge[:to]].to_vector
+                (vector_to - vector_from).magnitude
+              end.sum
+            end
+          end
+        end
+      end
+    end
+  end
+end

data/lib/structuraid_core/design_codes/aci_318_19/rc/minimum_steel_cover.rb

@@ -0,0 +1,55 @@
+module StructuraidCore
+  module DesignCodes
+    module ACI31819
+      module RC
+        class MinimumSteelCover
+          include DesignCodes::Utils::CodeRequirement
+          use_schema DesignCodes::Schemas::RC::MinimumSteelCoverSchema
+
+          CODE_REFERENCE = 'ACI 318-19 21.5.1.3.1'.freeze
+
+          def call
+            compute_cover
+          end
+
+          private
+
+          # Table 20.5.1.3.1
+          def compute_cover
+            return 75 if concrete_casting_against_soil
+            return expossed_to_environment_or_soil if !concrete_casting_against_soil && environment_exposure
+            return 50 unless structural_element
+
+            non_expossed_to_environment_or_soil
+          end
+
+          def expossed_to_environment_or_soil
+            return 40 if maximum_rebar_diameter&.<= 16
+
+            50
+          end
+
+          def non_expossed_to_environment_or_soil
+            if %i[slab wall joist].include?(structural_element)
+              slab_wall_joist_minimum_cover
+            elsif %i[beam column tensor_joint pedestal].include?(structural_element)
+              40
+            end
+          end
+
+          def slab_wall_joist_minimum_cover
+            return 20 if maximum_rebar_diameter&.<= 36
+
+            40
+          end
+
+          def shell_minimum_cover
+            return 13 if maximum_rebar_diameter&.<= 16
+
+            20
+          end
+        end
+      end
+    end
+  end
+end

data/lib/structuraid_core/design_codes/aci_318_19/rc/reduction_factor.rb

@@ -0,0 +1,88 @@
+module StructuraidCore
+  module DesignCodes
+    module ACI31819
+      module RC
+        class ReductionFactor
+          include DesignCodes::Utils::CodeRequirement
+          use_schema DesignCodes::Schemas::RC::ReductionFactorSchema
+
+          CONTROL_STRENGTH_CASES = %i[
+            compression_controlled
+            tension_controlled
+            transition_controlled
+            crushing_controlled
+            shear_nonseismic_controlled
+            torsion_controlled
+            corbel_bracket_controlled
+            strud_and_tie_controlled
+          ].freeze
+
+          MAX_STRAIN_BEFORE_TRANSITION = 0.002 # 21.2.2.1 - 21.2.2.2
+          MIN_STRAIN_AFTER_TRANSITION = 0.005 # 21.2.2.1 - 21.2.2.2
+
+          CODE_REFERENCE = 'ACI 318-19 21.2'.freeze
+
+          # ACI 318-19 21.2
+          def call
+            unless CONTROL_STRENGTH_CASES.include?(strength_controlling_behaviour)
+              raise UnrecognizedValueError.new(strength_controlling_behaviour, :strength_controlling_behaviour)
+            end
+
+            # Table 21.2.1 (a) (h)
+            return tension_controlled_factor if strength_controlling_behaviour == :tension_controlled
+            # Table 21.2.1 (a)
+            return compression_controlled_factor if strength_controlling_behaviour == :compression_controlled
+            # Table 21.2.1 (d)
+            return crushing_controlled_factor if strength_controlling_behaviour == :crushing_controlled
+            # Table 21.2.1 (b) (c) (f) (g)
+            if %i[
+              shear_nonseismic_controlled
+              torsion_controlled
+              corbel_bracket_controlled
+              strud_and_tie_controlled
+            ].include?(strength_controlling_behaviour)
+              return shear_nonseismic_controlled_factor
+            end
+
+            # Table 21.2.1 (a)
+            transition_controlled_factor
+          end
+
+          private
+
+          def shear_nonseismic_controlled_factor
+            0.75
+          end
+
+          def crushing_controlled_factor
+            0.65
+          end
+
+          def tension_controlled_factor
+            0.90
+          end
+
+          def compression_controlled_factor
+            return 0.75 if is_coil_rebar
+
+            0.65
+          end
+
+          def transition_controlled_factor
+            raise MissingParamError, :strain unless strain
+
+            transition_controlled_by_strain
+          end
+
+          def transition_controlled_by_strain
+            return compression_controlled_factor if strain < MAX_STRAIN_BEFORE_TRANSITION
+            return tension_controlled_factor if strain > MIN_STRAIN_AFTER_TRANSITION
+
+            transition_rate = is_coil_rebar ? 50 : 250 / 3
+            compression_controlled_factor + (strain - MAX_STRAIN_BEFORE_TRANSITION) * transition_rate
+          end
+        end
+      end
+    end
+  end
+end

data/lib/structuraid_core/design_codes/nsr_10/rc/footings/bending_reinforcement_ratio.rb

@@ -0,0 +1,66 @@
+module StructuraidCore
+  module DesignCodes
+    module NSR10
+      module RC
+        module Footings
+          class BendingReinforcementRatio
+            MINIMUM_RATIO = 0.0025
+            CODE_REFERENCE = 'NSR-10 C.15'.freeze
+
+            include DesignCodes::Utils::CodeRequirement
+            use_schema DesignCodes::Schemas::RC::Footings::BendingReinforcementRatioSchema
+
+            # NSR-10 C.15
+            def call
+              [
+                solve_cuadratic_equation_for_steel_reinforcement_ratio,
+                MINIMUM_RATIO
+              ].max
+            end
+
+            private
+
+            ## the target is to solve [(1 - sqrt(1 - 4ac))/2a, (1 + sqrt(1 - 4ac))/2a]
+
+            def solve_cuadratic_equation_for_steel_reinforcement_ratio
+              [
+                equation_solver_option_1.negative? ? Float::INFINITY : equation_solver_option_1,
+                equation_solver_option_2
+              ].min
+            end
+
+            def equation_solver_option_1
+              (1 - equation_component_root) / (2 * equation_component_a)
+            end
+
+            def equation_solver_option_2
+              (1 + equation_component_root) / (2 * equation_component_a)
+            end
+
+            def equation_component_a
+              0.59 * design_steel_yield_strength / design_compression_strength
+            end
+
+            def equation_component_c
+              section_area = width * effective_height # mm**2
+              reduced_steel_strength = capacity_reduction_factor * design_steel_yield_strength # N/(mm**2)
+              flexural_moment / (reduced_steel_strength * section_area * effective_height) # N/(mm**2)
+            end
+
+            def equation_component_root
+              unless 4 * equation_component_a * equation_component_c < 1
+                raise RequirementNotFulfilledError.new(
+                  :flexural_moment,
+                  "Moment #{flexural_moment} is too hight for this element",
+                  CODE_REFERENCE
+                )
+              end
+
+              Math.sqrt(1 - 4 * equation_component_a * equation_component_c)
+            end
+          end
+        end
+      end
+    end
+  end
+end

data/lib/structuraid_core/design_codes/nsr_10/rc/footings/one_way_shear_capacity.rb

@@ -0,0 +1,31 @@
+module StructuraidCore
+  module DesignCodes
+    module NSR10
+      module RC
+        module Footings
+          class OneWayShearCapacity
+            CODE_REFERENCE = 'NSR-10 C.11.11'.freeze
+
+            include DesignCodes::Utils::CodeRequirement
+            use_schema DesignCodes::Schemas::RC::Footings::OneWayShearCapacitySchema
+
+            # NSR-10 C.11.11
+
+            def call
+              compute_concrete_shear_capacity
+            end
+
+            private
+
+            def compute_concrete_shear_capacity
+              effective_area = effective_height * width
+              concrete_strength = light_concrete_modification_factor * Math.sqrt(design_compression_strength)
+
+              0.17 * concrete_strength * effective_area
+            end
+          end
+        end
+      end
+    end
+  end
+end

data/lib/structuraid_core/design_codes/nsr_10/rc/footings/punching_critical_section_perimeter.rb

@@ -0,0 +1,103 @@
+module StructuraidCore
+  module DesignCodes
+    module NSR10
+      module RC
+        module Footings
+          class PunchingCriticalSectionPerimeter
+            include DesignCodes::Utils::CodeRequirement
+            use_schema DesignCodes::Schemas::RC::Footings::PunchingCriticalSectionPerimeterSchema
+
+            EDGES_INDEXES = {
+              top: { from: :top_right, to: :top_left },
+              left: { from: :top_left, to: :bottom_left },
+              bottom: { from: :bottom_left, to: :bottom_right },
+              right: { from: :bottom_right, to: :top_right }
+            }.freeze
+
+            def call
+              column_location = footing.find_or_add_column_location(column_absolute_location, column_label)
+              perimeter_vertices = build_initial_perimeter(column_location)
+
+              edges_indexes = EDGES_INDEXES
+              edges_indexes = select_edges_indexes_into_the_footing(edges_indexes, perimeter_vertices)
+              edges_indexes = update_edges_indexes_vertices(edges_indexes, perimeter_vertices)
+              compute_perimeter(edges_indexes, perimeter_vertices)
+            end
+
+            private
+
+            def update_edges_indexes_vertices(edges_indexes, perimeter_vertices)
+              edges_indexes.each do |_edge_name, edge|
+                update_location_to_limit(perimeter_vertices[edge[:from]])
+                update_location_to_limit(perimeter_vertices[edge[:to]])
+              end
+            end
+
+            def select_edges_indexes_into_the_footing(edges_indexes, perimeter_vertices)
+              edges_indexes.select do |_edge_name, edge|
+                start_vertex_inside_footing = footing.inside_me?(perimeter_vertices[edge[:from]])
+                end_vertex_inside_footing = footing.inside_me?(perimeter_vertices[edge[:to]])
+
+                start_vertex_inside_footing || end_vertex_inside_footing
+              end
+            end
+
+            def build_initial_perimeter(column_location)
+              perimeter_vertices_relative_to_column_location.map do |label, perimeter_vertex|
+                location = local_coordinates_system.find_or_add_location_from_vector(
+                  perimeter_vertex + column_location.to_vector,
+                  label: "column_#{column_label}_punching_#{label}"
+                )
+                [label, location]
+              end.to_h
+            end
+
+            def update_location_to_limit(location)
+              location.value_1 = update_to_constraint(location.value_1, 0.5 * footing.length_1)
+              location.value_2 = update_to_constraint(location.value_2, 0.5 * footing.length_2)
+            end
+
+            def update_to_constraint(initial_value, constraint_value)
+              if initial_value > constraint_value
+                constraint_value
+              elsif initial_value < -constraint_value
+                -constraint_value
+              else
+                initial_value
+              end
+            end
+
+            def perimeter_vertices_relative_to_column_location
+              {
+                top_right: perimeter_vertices_relative_to_column_location_with(way_1: 1, way_2: 1),
+                top_left: perimeter_vertices_relative_to_column_location_with(way_1: -1, way_2: 1),
+                bottom_left: perimeter_vertices_relative_to_column_location_with(way_1: -1, way_2: -1),
+                bottom_right: perimeter_vertices_relative_to_column_location_with(way_1: 1, way_2: -1)
+              }
+            end
+
+            def perimeter_vertices_relative_to_column_location_with(way_1:, way_2:)
+              Vector[
+                way_1 * 0.5 * (column_section_length_1 + footing.effective_height),
+                way_2 * 0.5 * (column_section_length_2 + footing.effective_height),
+                0.0
+              ]
+            end
+
+            def local_coordinates_system
+              footing.coordinates_system
+            end
+
+            def compute_perimeter(edges_indexes, perimeter_vertices)
+              edges_indexes.map do |_edge_name, edge|
+                vector_from = perimeter_vertices[edge[:from]].to_vector
+                vector_to = perimeter_vertices[edge[:to]].to_vector
+                (vector_to - vector_from).magnitude
+              end.sum
+            end
+          end
+        end
+      end
+    end
+  end
+end

data/lib/structuraid_core/design_codes/nsr_10/rc/minimum_steel_cover.rb

@@ -0,0 +1,57 @@
+module StructuraidCore
+  module DesignCodes
+    module NSR10
+      module RC
+        class MinimumSteelCover
+          include DesignCodes::Utils::CodeRequirement
+          use_schema DesignCodes::Schemas::RC::MinimumSteelCoverSchema
+
+          CODE_REFERENCE = 'NSR-10 C.7.7.1'.freeze
+
+          def call
+            compute_cover
+          end
+
+          private
+
+          def compute_cover
+            return 75 if concrete_casting_against_soil # c.7.7.1.a
+            return expossed_to_environment_or_soil if !concrete_casting_against_soil && environment_exposure # c.7.7.1.b
+            return 50 unless structural_element
+
+            non_expossed_to_environment_or_soil
+          end
+
+          def expossed_to_environment_or_soil
+            return 40 if maximum_rebar_diameter&.<= 16
+
+            50
+          end
+
+          # c.7.7.1.c
+          def non_expossed_to_environment_or_soil
+            if %i[slab wall joist].include?(structural_element)
+              slab_wall_joist_minimum_cover
+            elsif %i[beam column].include?(structural_element)
+              40
+            elsif %i[shell thin_shell].include?(structural_element)
+              shell_minimum_cover
+            end
+          end
+
+          def slab_wall_joist_minimum_cover
+            return 20 if maximum_rebar_diameter&.<= 36
+
+            40
+          end
+
+          def shell_minimum_cover
+            return 13 if maximum_rebar_diameter&.<= 16
+
+            20
+          end
+        end
+      end
+    end
+  end
+end