MODULAR PROGRAMMECOURSEWORK ASSESSMENT SPECIFICATIONModule Details Module CodeUBGMW9-15-3Module TitleComputational Civil EngineeringModule LeaderGhassan NounuModule TutorsAndre Jesus & Xiaodong XuComponent and Element NumberComponent BWeighting (% of the Module’s assessment)75 %Element DescriptionCourseworkTotal Assignment time50 h Dates Date Issued to StudentsThursday 25th March, 2021Date to be Returned to Students(to be confirmed)Submission PlaceBlackboard onlinesubmissionSubmission Date(to be confirmed)Submission Time14:00hrs (GMT) DeliverablesWritten Report (2000 words) — File format: PDF and programming codefile(s)/spreadsheet(s). Please use a single zip file to upload your work using an appropriate file name e.g. your student number.Module Leader SignatureGhassan Nounu1PreambleAll assessments on this module are individual work. The work you submit must be your ownwork. Submitting work that is copied in part or whole from another student with or withouttheir permission is an assessment offence.You must fully attribute/reference all sources of information used during the completion ofyour submission, failure to do so constitutes plagiarism, which is an assessment offence.If you are not familiar with the definitions of plagiarism and collusion, more information canbe found here: http://www1.uwe.ac.uk/students/academicadvice/assessments/assessmentoffences.aspxPlease ensure you are familiar with assessment procedures and policies, which can befound here: http://www1.uwe.ac.uk/students/academicadvice/assessments/assessmentsguide.aspxStructure of assessmentsThis module is assessed by component B and component A, weighted as 75 % and 25 % of thefinal mark, respectively. To pass the module you must achieve a final mark of at least 40 %. (Nocomponent should be less than 35 %). This assessment brief is for component B (2000 words,excluding appendices and references). You will receive formative feedback for this assessmentduring lectures and group work sessions.The following section describes the problems you are expected to solve using MATLAB or othergeneral programming language and specific details of the tasks and outputs to feed in to yourreport.The coursework portfolio described here asks you to consider two problems entitled:1. Structural analysis under variable loads2. Calculation of second moment of areaResults must also be evidenced by source programming code file(s)/spreadsheet(s). Thesecode routines developed for both elements must be in a text selectable form (no images orscreenshots will be accepted).Structural analysis under variable loadsWhen dealing with variable (live) loads the internal forces or reactions that a structure generates will vary according to a probability distribution. In order to determine the design value of aDr Andre Jesus & Dr Xiaodong Xu 2 University of the West of Englandlive load, we multiply the characteristic load by a partial safety factor, leading to a design liveload value, which represents a worst case scenario associated with that load. However, whenmultiple loads interact together it might not be evident if applied partial safety factors shouldamplify or minimise each individual load. A probabilistic analysis has then to be performed,so as to determine how the internal forces or reactions will vary, and determine the outputvalues of these distributions which have a small probability, on an absolute basis, of beingexceeded. A workflow of this process is shown in Fig. 1.+– p1 ∼ N(µp; σp)p2 M1 M2M11 – Generate samples forinput variable UDL2 – Compute outputreactions/internal forces3 – Plot outputs histogramsand estimate the 5% thresholdoutput value-40 -35 -30 -250 M1 [kNm]50100150200250300350400450500-500 0 500 100015000 M2 [kNm]100200300400500600700800M2M1M2M1M2Figure 1: Diagram of probabilistic analysis for a supported beam subjected to twovariable uniformly distributed loads (UDL). At the centre, three bending moment diagrams illustrate the variability which can be expected due to the UDL distributions.The histograms represent the distributions of the bending moments at one supportand one mid-span of the beam, and the vertical dashed lines are the correspondingthreshold output values.Consider the ISOSTATIC1 structures shown in Figs. 2, 3, 4, 5, and the output reactions/internal forces presented in Table 1. You are asked to assess the variability of one thesestructures’ outputs when subjected to the shown loads. Each of the loads is assumed to followa normal distribution, e.g. for a uniform distributed load assume p ∼ N(µp; σp) with meanµp and standard deviation σp. The structures, geometry and load parameters shown in thefigures are an example for illustration — you have been assigned an individual structure,with a set of geometric and load properties that are individual to you and can be downloadedfrom: Blackboard > Assignments > my_values.html file, or by following the URL https://blackboard.uwe.ac.uk/bbcswebdav/pid-8051463-dt-content-rid-22271651_2/courses/UBGMW9-15-3_20jan_1/my_values.html. Note that structure 1 has a rollersupport at location 2. Similarly, the supports at location 6 structure 2, and location 1 structure3 are also rollers.Using MATLAB or other programming language generate 10 000 data points for each load ac-1An isostatic structure can be solved with simple equations of equilibrium: equilibrium of vertical forces, equilibrium of horizontal forces and equilibrium of moments.Dr Andre Jesus & Dr Xiaodong Xu 3 University of the West of EnglandPa bd 13624 5c7MFigure 2: Structure 1P2P c1pP3ba a a a123 4 56Figure 3: Structure 2p Pfa b c d 26 547 83 1eMFigure 4: Structure 3cording to its distribution parameters and compute the probability distribution of the corresponding output reactions/internal forces.Your report should include• A description of the equations and histograms for each structural output reaction/internal force. (30 %)• An estimate of the 5 % threshold output value, which is defined here as the value whichis exceeded, on an absolute basis, by only 5 % of the load combination realizations.(20 %)Dr Andre Jesus & Dr Xiaodong Xu 4 University of the West of Englanda be f 2 Pp1345dcFigure 5: Structure 4 StructureOutputs1Bending moment at section 6 and shear force atsection 42Bending moment at section 2 and shear force atsection 13Vertical reaction at 3 and bending moment at 6towards 54Reaction force at support 1 and bending moment at 4 towards 2 Table 1: Output reactions and internal forces• A pseudocode or flowchart of the algorithm that underlies your analysis. (10 %)Calculation of second moment of area• Aim: In this task you will calculate the neutral axis, second moment of area and maximum stresses for a randomly reinforced concrete beam• Theory: In a reinforced concrete beam, steel rebars are used to withstand tensile loadswhile concrete takes the compressive loads. The background theory is summarised byDr Arnaud Marmeir’s lecture notes, shown in Figure 6.Problem: The concrete beam in question has been randomly reinforced with steel rebars asshown in Figure 7, and we need to calculate its second moment of area and neutral axis usingcode.Methodology: Under tension, ignore the contribution of concrete and calculate the equivalentDr Andre Jesus & Dr Xiaodong Xu 5 University of the West of EnglandFigure 6: Theory of reinforced concrete design.MyxyFigure 7: Reinforced concrete beam.area of steel rebars according to the Theory; Under compression, both the contribution ofconcrete and steel rebars would need to be considered which means the Theory needs to beslightly revised if the steel rebars are also under compression.Inputs: The Young’s modulus of steel Es = 200 GPa, and the Young’s modulus of concreteEc = 25 GPa. It is assumed that Es is the same for tension and compression, but Ec can beignored under tension. The overall cross-sectional dimensions of the reinforced concrete beamare 600 mm in width and 140 mm in height. A total of 5 steel rebars with a diameter of 16 mmare evenly spaced (150 mm) horizontally but each of their vertical position through the depthis a random integer between 0 and 140. The direction of the moment My has been given.Expected outputs:1. Plot the distribution of 5 steel rebars. (10 %)Dr Andre Jesus & Dr Xiaodong Xu 6 University of the West of England2. Calculate the neutral axis by solving the equation according to the Theory. (10 %)3. Calculate the second moment of area Iyy. (10 %)4. Calculate the maximum tensile stress in the steel rebars and maximum compressivestress in the concrete under a constant moment of My = 9 kN m (20 %)Dr Andre Jesus & Dr Xiaodong Xu 7 University of the West of England
