StructureHelper/StructureHelperLogics/Models/BeamShears/Logics/GetLongitudinalForceFactorLogic.cs

using StructureHelperCommon.Infrastructures.Exceptions;
using StructureHelperCommon.Infrastructures.Interfaces;
using StructureHelperCommon.Models;
using StructureHelperCommon.Models.Loggers;

namespace StructureHelperLogics.Models.BeamShears.Logics
{
    public class GetLongitudinalForceFactorLogic : IGetLongitudinalForceFactorLogic
    {
        private const double fstRatioInCompression = 0.25;
        private const double sndRatioInCompression = 0.75;
        private double sectionArea;
        private ICheckEntityLogic<IInclinedSection> checkInclinedSectionLogic;

        public IShiftTraceLogger? TraceLogger { get; set; }
        public IInclinedSection InclinedSection { get; set; }
        public double LongitudinalForce { get; set; }

        public GetLongitudinalForceFactorLogic(IShiftTraceLogger? traceLogger)
        {
            TraceLogger = traceLogger;
        }


        public double GetFactor()
        {
            TraceLogger?.AddMessage(LoggerStrings.LogicType(this), TraceLogStatuses.Service);
            TraceLogger?.AddMessage("Logic of calculating of factor of influence of longitudinal force according to SP 63.13330.2018");
            Check();
            if (LongitudinalForce == 0)
            {
                TraceLogger?.AddMessage("Longitudinal force is zero", TraceLogStatuses.Service);
                return 1;
            }
            sectionArea = InclinedSection.WebWidth * InclinedSection.FullDepth;
            TraceLogger?.AddMessage($"Area of cross-section Ac = {InclinedSection.WebWidth} * {InclinedSection.FullDepth} = {sectionArea}(m^2)");
            if (LongitudinalForce < 0)
            {
                TraceLogger?.AddMessage($"Longitudinal force N={LongitudinalForce}(N) is negative (compression)", TraceLogStatuses.Service);
                return GetNegForceResult();
            }
            else
            {
                TraceLogger?.AddMessage($"Longitudinal force N={LongitudinalForce}(N) is positive (tension)", TraceLogStatuses.Service);
                return GetPosForceResult();
            }
        }

        private void Check()
        {
            checkInclinedSectionLogic ??= new CheckInclinedSectionLogic(TraceLogger);
            checkInclinedSectionLogic.Entity = InclinedSection;
            if (checkInclinedSectionLogic.Check() == false)
            {
                string errorString = checkInclinedSectionLogic.CheckResult;
                TraceLogger?.AddMessage(errorString, TraceLogStatuses.Error);
                throw new StructureHelperException(errorString);
            }
        }

        private double GetPosForceResult()
        {
            double stressInConcrete = LongitudinalForce / sectionArea;
            TraceLogger?.AddMessage($"Average stress in concrete (positive in tension) Sigma = {LongitudinalForce}(N) / {sectionArea}(m^2) = {stressInConcrete}(Pa)");
            double concreteStrength = InclinedSection.ConcreteTensionStrength;
            TraceLogger?.AddMessage($"Concrete strength Rbt = {concreteStrength}(Pa)");
            double stressRatio = stressInConcrete / concreteStrength;
            TraceLogger?.AddMessage($"Stress ratio rt = {stressInConcrete} / {concreteStrength} = {stressRatio}(dimensionless)");
            double factor = 1 - 0.5 * stressRatio;
            factor = Math.Max(factor, 0);
            TraceLogger?.AddMessage($"Factor value fi_n = {factor}(dimensionless)");
            return factor;
        }

        private double GetNegForceResult()
        {
            double stressInConcrete = (-1) * LongitudinalForce / sectionArea;
            TraceLogger?.AddMessage($"Average stress in concrete (positive in compression) Sigma = {LongitudinalForce}(N) / {sectionArea}(m^2) = {stressInConcrete}(Pa)");
            double concreteStrength = InclinedSection.ConcreteCompressionStrength;
            TraceLogger?.AddMessage($"Concrete strength Rb = {concreteStrength}(Pa)");
            double stressRatio = stressInConcrete / concreteStrength;
            TraceLogger?.AddMessage($"Stress ratio rc = {stressInConcrete} / {concreteStrength} = {stressRatio}(dimensionless)");
            double factor;
            if (stressRatio < fstRatioInCompression)
            {
                TraceLogger?.AddMessage($"Stress ratio rc = {stressRatio} < {fstRatioInCompression}");
                factor = 1 + stressRatio;
            }
            else if (stressRatio > sndRatioInCompression)
            {
                factor = 5 * (1 - stressRatio);
                factor = Math.Max(factor, 0);
            }
            else
            {
                factor = 1 + fstRatioInCompression;
            }
            TraceLogger?.AddMessage($"Factor value fi_n = {factor}(dimensionless)");
            return factor;
        }
    }
}