Program Package For Analysis Research Articles

Building a sustainable future: Exploring the nexus between climate awareness and thermal insulation choices in new home construction

Climate change is considered a major threat to the environment and society as a whole. Reducing building energy consumption through the use of thermal insulation is one way of combating climate change. This study examines the impact of climate change awareness on the pre-construction decision of Lebanese homeowners to use thermal insulation in their new homes. A model building consisting of a typical 160 m2 single-family home was used and its thermal loads were calculated in two cases: insulated and uninsulated case, in three different locations in Lebanon: Beirut, Chtaura, and Bcharre-Arz. The thermal loads were calculated using two methods: the degree-days empirical method, and the Hourly Analysis Program software numerical method. The computations showed energy reduction ranging between 23% in the coastal region to 52% in high mountains when thermal insulation of thickness 5 cm was used. The investment payback period ranged from 4.3 years in the coastal region to 1.5 years in the high mountains. A survey addressed to 100 respondents was used to understand how Lebanese homeowners’ decision to use thermal insulation was affected by knowledge of the obtained results. The survey results showed that climate change awareness has a high impact on this decision; however, the cost seems to play an important role in the owners’ final decision on whether or not to use thermal insulation. Spreading awareness about climate change and the payback periods of building thermal insulation will increase its use since the number of respondents refusing to use thermal insulation drastically decreased after being informed about its payback period.

Optimization of air-conditioning load for a residential building using Taguchi method

Economic growth has been raised in many countries in recent years, prompting people to seek out human comforts. For human comfort, residential buildings consume a significant amount of energy for air-conditioning systems. Building materials are the most significant factor in deciding the air-conditioning load. In this research work, a typical residential flat model is prepared using AutoCAD software and the building energy simulation is done for this model using hourly analysis program software. The results obtained from the building simulation are optimized using the Taguchi method. For optimizing the air-conditioning load, five influencing factors such as the overall heat transfer coefficient of wall (A), thickness of roof insulation material (B), heat transfer coefficient of window glass (C), bypass factor (D) and cooling coil supply temperature (E) of air-conditioning equipment at three levels are selected. For the experimental tests, an orthogonal L27 arrangement is chosen for experimental trials and in every trial, air-conditioning load of the typical residential flat has been calculated. After that, using Minitab 17 statistical software, the signal-to-noise ratios are obtained for different factors by selecting smaller-is-better option. After studying the response table for Taguchi tests, the optimum levels for all five factors are obtained. The analysis of variance method was also used to identify the main factors affecting the results significantly by finding out the contribution of all five factors in percentage. The results indicated that the most efficient set of all five factors is A1B3C1D2E3 and the most affecting factor is the wall material of the residential flat having a contribution of 91.91%.

Efficacy of increasing levels of exposure therapy in the treatment of maladaptive behaviors and anxiety.

The present study examined the effectiveness of increasing levels of exposure therapy, which is applied for the treatment of maladaptive behaviors and anxiety. A total of 16 sessions were applied to the study group in the experimental group three times a week for 10 weeks. Patients aged ≥18 years whom the referring clinician evaluated as meeting the criteria for the Diagnostic and Statistical Manual of Mental Disorders, fourth edition (DSM-V-TR) Generalized Anxiety Disorder were included in the study. For the control group, demographic characteristics and Spielberger's State-Trait Anxiety Inventory were applied in the first session, followed by Spielberger's State-Trait Anxiety Inventory as a post-test and follow-up. Electroencephalography (EEG) recordings of the study group were obtained at the cortical level. Electrodes for EEG measurements were recorded using the International 10/20 Electrode Placement System. EEG data were obtained using the EEG Analysis Program software. Following the data collection phase, all data were entered into cells based on items using SPSS 25 software. When the findings obtained in the study were examined, it was determined that the increasing levels of exposure and behavioral therapy applied for maladaptive anxiety decreased the anxiety levels compared to those before therapy. This finding can be interpreted as that the cortical function-oriented application method for anxiety effectively reduced the anxiety levels of the study group. However, EEG asymmetry revealed a change in the data before and after the application. These findings demonstrate that the application affects the EEG asymmetry changes at the cortical level.

Adaptive capabilities of the body and aggressive personality manifestations of modern students at various stages of education in a higher education institution

The determination of the subsequent comprehensive assessment of the peculiarities of the course of adaptation processes, which are characteristic of modern student youth who acquire a certain profession, are closely related to the establishment of the regularities of the course of numerous adaptive mechanisms in the plane of implementation of psychophysiological, mental, and socio-psychological adaptation. The purpose of the work is to carry out a comprehensive hygienic assessment of the adaptation capabilities of the organism based on the definition of its social and psychological component and aggressive personality manifestations of modern students at various stages of education in a higher education institution. Scientific studies were conducted using questionnaires by Rogers and Diamond and Bass-Darky. 307 students, including 150 young women and 157 young men, who studied in the first, third and sixth years of the medical faculty were under supervision during the observation period. The analysis of the obtained data involved the application of descriptive statistics procedures based on the application of the statistical analysis program package “Statistica 6.1” (license number ВXXR901E245722FA). It was established that the highest level of expression of integral indicators of social and psychological adaptation of students, regardless of the time of study in a higher education institution, is characteristic of the characteristics of adaptation manifestations regarding the desire to dominate and accept others; at the same time, the lowest level of expression is for the characteristics of internality and therefore, determines the presence of pronounced signs of transformations of the adaptive content, which mark the tense course of the processes of social and psychological adaptation of students to the conditions of stay in a medical institution of higher education. It was found that among the overwhelming number of students who were studied, the indicators of aggressive personality manifestations exceed the values typical for normative age-sex indicators. This situation determines the need for the development of methods for a comprehensive assessment of the state of adaptive resources of the body of young men and women studying and scientific substantiation of health-preserving technologies, based on which measures of psychohygienic correction should take the leading place.

Tooth Color Change and Erosion: Hydrogen Peroxide Versus Non-peroxide Whitening Strips.

The study evaluated the efficacy and potential erosion of non-peroxide strips compared to hydrogen peroxide (HP) whitening strips (WSs). Color evaluation samples (N=64) were distributed into four groups and treated according to manufacturer's directions. NC: Negative control treated with water; BT: Non-peroxide Brilliant Dissolving Strips; FM: Non-peroxide Fancymay Teeth WSs; WS: Crest 3D Brilliance HP White Strips. A contact-type spectrophotometer was used to measure color at baseline (T1), 1-day posttreatment (T2), and 1-week posttreatment (T3). Teeth were cut to a rectangular block for micro-CT erosion assessment. The samples (N=30) were divided into five groups. In addition to the four groups for color assessment, a positive control (PC) treated with 0.25% citric acid was added. The samples were scanned, reconstructed, and measured for erosion depth using a micro-CT analysis program software. Kruskal-Wallis test was used to determine differences in color change and erosion depth among the groups. Tests of hypotheses were two-sided with an alpha level of 0.05. The mean ΔE*ab at 1-day/1-week posttreatment were 2.4/2.5, 2.8/2.9, 2.8/3.2, and 8.6/11.0 for NC, BT, FM, and WS, respectively. There was a statistically significant difference for ΔE*ab at 1-day and 1-week posttreatment (p<0.001). Group WS had the highest color change, while the other three groups did not differ from each other (p>0.05). Mean erosion depths in microns were 0.52, 0.58, 0.42, 0.49, and 29.55 for NC, BT, FM, WS, and PC, respectively. There was a statistically significant difference among the groups (p=0.004). Group PC had the greatest erosion, while the other groups had negligible erosion that did not differ from each other (p>0.05). Peroxide WSs had superior whitening efficacy compared to non-peroxide strips. None of the tested products compromised tooth structure integrity through enamel erosion.

Protection Coordination Method Using Symmetrical Components in Loop Distribution System

Power utilities worldwide commonly use the radial distribution system because of its advantages of being simple in structure and having relatively inexpensive installation costs. It has a disadvantage in that its power supply reliability is low because the load side of the fault section will suffer from an outage in the event of a fault in the system. However, recently, with ICT (Information and Communication Technologies) development, system reliability is required to be high as the outage-susceptible loads increase. In addition, the increase in the connection of distributed resources such as renewable energy and electric vehicles is making it impossible to predict the power flow and reducing line utilization. Therefore, a loop power distribution system is proposed as a measure to solve this problem. Because all buses (nodes) in a loop distribution system have two or more power supply routes, they are more reliable than the radial system. It allows them to improve line utilization by connecting lines with different load peak times. However, in the case of a fault in the loop distribution system, the fault current is supplied from both directions, making it impossible to properly isolate the fault section with the protection method of the conventional distribution system. The permissive overreach transfer trip (POTT) method using communication to compensate for the limitations of conventional protection devices, and the other method using directional distance relay, is proposed. However, these methods operate by determining the direction of the fault current but have a disadvantage. It is difficult to detect a fault due to the effects of ground faults and distributed generation (DG) occurring in other lines. Therefore, in this paper, we propose a protection coordination algorithm that uses the negative-sequence component of voltage and current that occur when an unbalanced fault occurs, rather than the determination of the directionality and use of communication. To validate this, we configured a system using PSCAD/EMTDC (Manitoba Hydro International Ltd., Winnipeg, Manitoba, Canada), a system analysis program package and verified the results depending on the type of faults with the proposed algorithm.

Efficiency And Productivity Analysis In Turkish Banking Sector With Data Envelopment Analysis And Malmquist Index

The aim of this study is to analyze the efficiency and productivity of 3 public, 6 private and 6 foreign deposit banks operating in the Turkish banking sector with the help of data envelopment analysis and Malmquist index.For this purpose, the efficiency of 15 deposit banks operating in the Turkish banking sector between 2014 and 2018 was measured and whether the efficiency of Malmquist productivity index changed over the years. In the study, input-oriented Charnes Cooper and Rhodes (CCR) model was used under the assumption of constant return to scale and 4 input 2 output variables were selected. Inputs in efficiency and productivity measurement are defined as; personnel expenses / total assets (%), total loans / total assets (%), equity / total assets (%), total deposits / total assets (%), outputs are defined as; the earning power of assets (net profit / total assets), the earning power of equity (net profit / equity) (%). The Windows Data Envelopment Analysis Program (Win4Deap) package program was used in the analysis and brokerage approach was adopted. While 4 banks were active under the Constant Return to Scale(CRS) assumption between 2014 and 2018, 8 banks were found active under the assumption of Variable Return to Scale (VRS).Inefficient banks; target values were calculated by slacks movement and radial movement values to their original values and it was found that the lambda values calculated by Win4Deap program and which banks are peers. The changes observed in Malmquist total factor productivity, technical efficiency, technological efficiency, pure efficiency, scale efficiency and total factor productivity were analyzed as a whole and decision-making units experienced improvement in the 2016-2017 period.

Energetic and Exergetic Analysis of the Low-temperature Heating System in a model library

Objectives: Energetic and exergetic analysis of the low-temperature heating system in a model library are conducted numerically in this research. Methods: The method of numerical investigation and simulation are employed to compute the energetic and exergetic analysis of the low-temperature heating system in a model library. The use of C++ programming and Hourly Analysis Program software were done the simulation and optimization. The results of this investigation are presented in this article. Findings: Simulation presents that the most of energy which is required for heating is capable to be declined using low-temperature heating system in a model library. It is shown that gas consumption decreases when lowtemperature radiators are applied and the value which is obtained is about 61,458,480 m3/yr. This means that the performance of the lowtemperature radiation is better than high-temperature radiator and it is about 57,150,000 m3/yr. The use of low-temperature radiators and high-temperature radiators are compared in this article. The radiator quality factor of the low-temperature radiators and hightemperature radiators are obtained 0.13 and 0.22, respectively. By using the low-temperature method, radiator quality factor is able to be declined 7% of the energy efficiency. Based on these cases, using the low-temperature system can be used in some regions in Iran. The economic investigation shows that the IRR is evaluated based on discount 30%, and interest rate 22% in 20 years. Application and improvement: The declination of energy consumption in building is significant. Therefore, a method which can use to reduce energy consumption is low-temperature heating in winter. Keywords: Low Temperature Heating, Low Temperature Radiator, Heating, Economic evaluation.

Algorithmic support of automatic yield control system of a gas well cluster

Gas wells operation is associated with the problem of regulating their yield. Debit value must provide for meeting various criteria with some of them being mutually exclusive. Initial debit value is set under a deposit development project and then is adjusted periodically at best based on the results of quarterly hydrodynamic well exploration. However, in order to enhance operation yield of a gas well cluster debit regulation must be carried out in real time. This problem hasn’t been solved yet. The article offers a yield control algorithm for a gas well cluster consisting of three separate units: yield enhancement, yield decreasing and yield distribution among wells. Transfer from one unit to another is carried out based on the results of measuring current gas consumption rate in wells and the pressure on isolation valve building’s plumes. The gap between the allowed values of well consumption (maximum and minimum) is identified on the basis of the research of a well’s operation. Algorithm approbation was carried out in a program package for simulation, design and analysis of the operation of above-ground infrastructure and PipeSim pipeline system. The results confirmed the feasibility of the proposed algorithm and the possibility of managing operation regime of every well independently.

Development and automation of algorithm for determining basis of nonlinear parameter functions of kinetic constants

Mathematical modeling of catalytic processes is necessary for the complete and accurate description, as well as for controlling the quality and physicochemical studied of catalysts. In the paper, theoretical issues of industrial catalysis are discussed. The work is devoted to theoretical graph analysis of informativity of kinetic parameters of the model of a complex chemical reaction. The aim is the development and automation of algorithm for determining basis of nonlinear parameter functions in solving inverse problems of chemical kinetics in order to define the number and form of independent combinations of rate constants of elementary stages. A program package for analysis of informativity of kinetic parameters of the mathematical model of a complex catalytic reaction is developed and described. The obtained functional relations between the kinetic parameters can be useful for experimentalists in physicochemical interpretation and analysis of mechanisms of chemical reactions. In other words, the proposed method allows independent combinations of kinetic constants to be distinguished that results in shortening the number of the model parameters and, as a consequence, enhance the accuracy of the mathematical model. The mechanism of hydrogen oxidation over a platinum catalyst is given as an example of the use of the software.

Quarter milking parameters by lactation in dairy cows

The study was conducted on 229 dairy cows in Lithuania. DeLaval milking robots were used for cows’ milking. Milk yield (kg), milking duration (min), milk flow (kg/min), peak flow rates (kg/min) by cow lactation in separate udder quarters were analyzed. The statistical analysis of the research data was performed using the data collection and analysis program package Statistical Package for Social Science 22 for Windows, Kolmogorov-Smirnov test and One-Factor Dispersion Analysis (ANOVA) model. Distribution of the udder quarters is very important in the cow milking process. Milk yield from cow front udder quarters was about 4.6 kg (42.2 %) and from rear udder quarters 6.32 kg (57.8 %) (P&amp;lt;0.05). The milk content in front quarters was significantly lower than in the rear udder quarters in all dairy cows (P&amp;lt;0.05). Cows of the first lactation showed lower milk production per milking in all udder quarters when compared to other lactation cows (P&amp;lt;0.05). The duration of milking of cows in third and more lactation in front quarters was shorter than that of the first lactation cows (P&amp;lt;0.005). Along the increase in the number of cows’ lactations, the milking time in the front quarters decreased. The average milk flow and the peak flow in the first lactation cows were found to be lower than that in older cows (P&amp;lt;0.05).

Optimal Allocation of Transformer and Its Effects on Distribution System Loss Minimization

The periodic growth in megacities has incremented the power consumption in the name of modernized techniques and hustle-free lifestyles. In the state of Tamil Nadu, the voltage and power regulation has created a drastic decrements in the efficiency for the tail end consumers. Thus, resulting in shortage of voltage and vehemently opposing the performance factors of home appliances, micro-industries and minute gadgets used on a day-to-day basis. In this paper, a solution is provided to magnificently ameliorate by catering a high voltage distribution system, simulating its results for serving, good to the Tamil Nadu Electricity Board having fed by the Electrical Transient Analysis Program Software. Therefore, in reference to the contingent PSO algorithm enabling to convert Low Voltage Distribution System to High Voltage Distribution Systems, sufficiently acquiring a redemption for utmost performance usage of the transformer. The relational composition of the supply in the LVDS and HDVS magnitudes demonstrate the necessary output demands. The line losses vividly broadcast the allocated transformer and the HDVS holds good in showcasing power loss in advanced tail and voltage prescriptions.

Impact of Network Reconfiguration: Case Study of Port-Harcourt Town 132/33kV Sub-transmission Substation and Its 33/11kV Injection Substation Distribution Networks

This paper examined the power flow status of the Port Harcourt Town (Zone 4) distribution networks to improve the performance. The network consists of 18 injection substations fed from 4 different sizes of transformers with a total power rating of 165 MVA, 132/33kV at the Port Harcourt Town sub-transmission substation. Gauss-seidel power flow algorithm was used to analyse the network in Electrical Transient Analyzer Program software (ETAP 12.6) to determine the various bus operating voltages, power flow, and over or under-loaded Transformers’ units. From the base-case simulation results obtained, it shows that these injection distribution transformers (PH Town 106.3%, RSU 90.5%, Marine Base 86.5%, UTC 87.9%, Nzimiro 89.5%, and Borokiri 88.7%) were overloaded on the network and the operating voltages observed for (PH Town 95.1%, RSU 83.0%, Marine Base 83.4%, UTC 82.8%, Nzimiro 85.2%, and Borokiri 82.1%) indicates low voltage profile. However, using network reconfiguration technique as proposed in this paper; there was reduction in the percentage loading of the said Transformers as it was upgraded to affect positively on its lifespan with (PH Town 44.1%, RSU 65.3%, Marine Base 60.7%, UTC 47.3%, Nzimiro 61.3%, and Borokiri 52.0%) loading, and the bus voltage profiles was improved for (PH Town 100%, RSU 98.4%, Marine Base 98.8%, UTC 98.2%, Nzimiro 98.6%, and Borokiri 99.1%) with additional facilities. It is recommended that the power infrastructure facilities in Port Harcourt Town distribution network be immediately upgraded to reduce losses and improve the electricity supply to consumers. Also, in regard to these analyses, the sub-transmission substation requires 240 MW of power for effective power delivery.

Embracing chemical and structural diversity with UCONGA: A universal conformer generation and analysis program

Molecular properties depend on molecular structure, so the first step in any computational chemistry investigation is to generate all thermally accessible conformers. Typically it is necessary to make a trade-off between the number of conformers to be explored and the accuracy of the method used to calculate their energies. Ab initio potential energy surface scans can, in principle, be applied to any molecule, but their conformational cost scales poorly with both molecular size and dimensionality of the search space. Specialized conformer generation techniques rely on parameterized force fields and may also use knowledge-based rules for generating conformers, and are typically only available for drug-like organic molecules. Neither approach is well-suited to generating or identifying chemically sensible conformers for larger non-organic molecules. The Universal CONformer Generation and Analysis (UCONGA) program package fills this niche. It requires no parameters other than built-in atomic van der Waals radii to generate comprehensive ensembles of sterically-allowed conformers, for molecules of arbitrary composition and connectivity. Analysis scripts are provided to identify representative structures from clusters of similar conformers, which may be further refined by subsequent geometry optimization. This approach is particularly useful for molecules not described by parameterized force fields, as it minimizes the number of computationally intensive ab initio calculations required to characterize the conformer ensemble. We anticipate that UCONGA will be particularly useful for computational and structural chemists studying flexible non-drug-like molecules.

Three-dimensional reconstruction for coherent diffraction patterns obtained by XFEL.

The three-dimensional (3D) structural analysis of single particles using an X-ray free-electron laser (XFEL) is a new structural biology technique that enables observations of molecules that are difficult to crystallize, such as flexible biomolecular complexes and living tissue in the state close to physiological conditions. In order to restore the 3D structure from the diffraction patterns obtained by the XFEL, computational algorithms are necessary as the orientation of the incident beam with respect to the sample needs to be estimated. A program package for XFEL single-particle analysis based on the Xmipp software package, that is commonly used for image processing in 3D cryo-electron microscopy, has been developed. The reconstruction program has been tested using diffraction patterns of an aerosol nanoparticle obtained by tomographic coherent X-ray diffraction microscopy.

دور اسلوب الموجه فی تنمیة رأس المال الفکری: دراسة میدانیة فی قطاع التعلیم الأساسی فی العراق

الهدف الرئيسي من هذه الدراسة هو معرفة دور الاسلوب الموجه في تنمية رأس المال الفکري في قطاع التعليم الأساسي محل الدراسة (بغداد، بابل، ذي قار) في العراق. ولتحقيق أهداف الدراسة استخدم الباحث المنهج الوصفي التحليلي، وقوائم الاستقصاء کأداة رئيسية لجمع البيانات من مجتمع الدراسة البالغ (20798) کما بلغت عينة الدراسة (377) مفردة وکانت عشوائية حيث تم استخدام برنامج التحليل الإحصائي (SPSS) کحزمه لإدخال ومعالجة وتحليل البيانات. أهم النتائج التي توصلت إليها الدراسة، وجود علاقة معنوية واثر معنوي بين الاسلوب الموجه وتنمية رأس المال الفکري. أهم توصيات الدراسة هو، تفعيل دور القطاعات المبحوثة متمثلة بقياداتها وملاکها التعليمي والتدريسي مع المجتمع، من خلال حث القيادات والملاک التعليمي والتدريسي بإقامة الندوات بما يسهم وخدمة المجتمع. Abstract: The main objective of this study is to know the role that style directed at intellectual capital development in the basic education sector under study (Baghdad, Babil, Dhi Qar) in Iraq. To achieve the objectives of the study, the researcher used descriptive and analytical approach, and lists of the survey as a key tool to collect data from the study population of (20798) also reached the study sample (384) Single and were randomized where the use of statistical analysis (SPSS) program package for entry, processing and analysis of data. The most important findings of the study, there were no significant relationship between the significant effect-oriented style and the development of intellectual capital. The most important recommendations of the study is, activating the role of the surveyed sectors represented Bkiedadtha educational, teaching and owners with the community, through its leaders and urged owners to establish educational and teaching seminars so as to contribute and community service.

Analysis and visualization of large networks with program package Pajek

AbstractPajek is a program package for analysis and visualization of large networks (networks containing up to one billion of vertices, there is no limit—except the memory size—on the number of lines). It has been available for 20 years. The program, documentation and supporting material can be downloaded and used for free for noncommercial use from its web page:http://mrvar.fdv.uni-lj.si/pajek/

A COMPARATIVE STUDY ON FORCE BASED DESIGN AND DIRECT DISPLACEMENT BASED DESIGN OF REINFORCED CONCRETE FRAMES

In recent times, many research works related to seismic resistant design of structures were done and various methods of seismic resistant design were predicted and this study explains about two major seismic design methods (i.e) Force Based Design(FBD) and Direct Displacement Based Design(DDBD) in which former is a conventional method while later one is a performance approach of design. Design and Analysis were done on two dimensional bare frames of four, eight and twelve stories based on following codes IS 456, IS 1893:2000, FEMA 356 and the two design approaches were studied. Analysis and design for this study was done using Structural Analysis Program software(SAP 2000). Both design approaches is validated using non linear time history analysis for 16 different ground motion of PGA =0.32g. Structural parameters like Drift Ratio, Ductility Demand and Base shear were compared within the frames of different stories and between design approaches.

Longitudinal analysis of concrete U-girder bridge decks

U-shaped prestressed concrete bridge decks, simply supported, are now being increasingly used in railways and highways. Simplified methods of analysis are commonly used in design practice. In the longitudinal direction, the U-girder is modelled as a beam. However, under eccentric loading, simple beam analysis is found to yield unconservative estimates of longitudinal forces, mainly due to torsional warping. An understanding of U-girder bridge deck behaviour under torsional loading is facilitated by the application of Vlasov's thin-walled beam theory. This theory can be easily invoked as a substitute for simple beam analysis, resulting in vastly improved predictions of longitudinal stresses. The two U-girder decks have also been modelled in Structural Analysis Program software using shell elements, and the results of three-dimensional finite-element analysis validate the thin-walled beam theory. The scope of the study is limited to a straight and thin-walled U-girder bridge deck of uniform section. The effects of variation in load position have also been studied, and it is established that the thin-walled beam theory is well suited to estimate and design the U-girder for longitudinal stresses.

Modeling of BHA dynamic behaviors

Lateral and torsional vibrations of drill strings used in oil well operations are considered. This work presents finite-element-based bottom hole assembly model which addresses natural frequencies determination, and torsional and lateral oscillations on the near vertical well. In terms of the excitations of the developed finite-element model, forces are applied on the bit with time history generated using an Auto-Regressive Moving-Average digital filter. Due to the random pattern of the forces at the drill bit, a stochastic dynamic approach is adopted in investigating the problem. The method of Newmark-beta direct integration is used. It is hoped that this study will enhance the interest in using stochastic dynamic techniques in drilling system analysis and design, as they can capture quite appropriately the inherent uncertainty of the bit forces and, potentially, of other sources. The implementation of torsional degree of freedom during the development of MATLAB code has been considered new attribute regarding previous works based on the same research line. The results from the code developed attained a good agreement with the commercial SPA2000 Advanced 14.0.0 Structural Analysis Program software.