!=======1=========2=========3=========4=========5=========6=========7=========8=========9=========10========11
|
|
- Marion Robyn O’Neal’
- 6 years ago
- Views:
Transcription
1 C:\files\classes\fem\program code\truss\truss solve\truss.f ! truss.f90! Kurt Gramoll! 2 Feb 2018! Requires input file as part of the command string (assumes.dat)! Format! line 1: number of nodes, number of elements! line for each node! node#, constraint(0,1,2,3)(none,x,y,both), locx, locy, forcex, forcey! line for each element! element#, left node#, right node#, E, Area! Sample From Logan book, example 3 5! 4 3! ! ! ! ! e6 2.0! e6 2.0! e6 2.0!=======1=========2=========3=========4=========5=========6=========7=========8=========9=========10========11 program truss implicit none!! NOT USED, example of set dimension for array (old Fortran method)!! assuming maximum nodes=50, max elements=50!! real*8 :: gstiff(1:100, 1:101)! global stiffness matrix with RHS, size = numnodes*2, array start at 1!! real*8 :: gstiff(100, 101)! global stiffness matrix with RHS, size = numnodes*2, array start at 1 real*8, allocatable, dimension(:,:) :: gstiff, gstifforig!! dynamic array method real*8, DIMENSION(1:4,1:4) :: elemstiff!! one for each element, old syntax real*8, allocatable, dimension(:) :: nodelocx, nodelocy!! location of node real*8, allocatable, dimension(:) :: load, disp!! node force (RHS) and displacement real*8, allocatable, dimension(:) :: elemlen, stress, angrad!! element length, stress, orientation angle (in radians) real*8, allocatable, dimension(:) :: youngmod, area!! element modulus and area real*8 :: temp, temp1, sum!! real = real*4 (single), real*8 (double, 64 bits), real*16 (quad)
2 C:\files\classes\fem\program code\truss\truss solve\truss.f real*8 :: cosang, sinang!! cosine and sine of each element bar real*8 :: elapsedtime, t1, t2!! used for calculation time integer, allocatable, dimension(:) :: noderest!! type of restriction, 0=non, 1=disp x, 2=dips y, 3=both integer, allocatable, dimension(:) :: leftnode, rightnode!! node order, left to right integer :: i, j, k, n, m, irow, icol, ielem, inode!! counters integer :: numnode, numelem, dof!! node, element and total degree of freedom character (len=100) :: strbuffer100!! string that holds command line argument character (len=30) :: filename, fileinname, fileoutname! file name (no ext), file name in and out w/ ext!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Get command line input parameters, call for each item call GETARG(1, strbuffer100)! special function to read in argurments read (strbuffer100,*) filename! reading in file name (no extension) fileinname = TRIM(fileName) // ".txt"! file in name with extension fileoutname = TRIM(fileName) // "_out.txt"! file out name with extension!! write to console that program (1st *) has started, free format (2nd *) write (*,*) 'truss program start reading in data from: ', fileinname CALL CPU_TIME(t1)! set timer, used to determine total time program runs!! Read in data file ======================================================================== open (UNIT=5, FILE=fileInName)! file number 5 standard for input, could use 'open (5, fileinname)' read (5,*) numnode, numelem! Just one line, needed to figure DOF (used for matrix setup!! need to allocate memory for dynamic matrices, needed for load and disp vector size before reading dof = 2*numNode allocate (gstiff(dof, dof))!! used when reducing global stiffness matrix to solvable set of equations (b.c. applied) allocate (gstifforig(dof, dof)) allocate (load(dof))! external loads at each node, b.c. (maybe unknown) allocate (disp(dof))! unknown node displacements (some known, b.c.) allocate (nodelocx(numnode)) allocate (nodelocy(numnode)) allocate (noderest(numnode))! node restriction (i.e. B.C.), 0 = none, 1 = x, 2 = y, 3 = both allocate (leftnode(numelem))! node 1 in local numbering allocate (rightnode(numelem))! node 2 in local numbering
3 C:\files\classes\fem\program code\truss\truss solve\truss.f allocate (youngmod(numelem)) allocate (area(numelem)) allocate (angrad(numelem))! calculate from nodelocx and nodelocy allocate (elemlen(numelem))! calculate from nodelocx and nodelocy allocate (stress(numelem))! to be found leftnode ( : ) = 0 rightnode ( : ) = 0!! Finish reading in data now that matrices sizes are set read (5,*) (n, noderest(i), nodelocx(i), nodelocy(i), load(2*(i 1)+1), load(2*(i 1)+2), i=1, numnode) do i = 1, numelem!! long way to read in read (5,*) n, leftnode(i), rightnode(i), youngmod(i), area(i)!! Write out copy of inputted data ============================================================ open (UNIT=6, FILE=fileOutName)!! file number 6 standard for output to file 'Number of points: ', numnode, '; elements: ', numelem "Nodes" write (6,62) (i, noderest(i), nodelocx(i), nodelocy(i), load(2*(i 1)+1), load(2*(i 1)+2), i=1, numnode) 62 format (" Pt Fixed loc X loc Y Force X Force Y" / (2I5, 4E12.4)) "Elements" do i = 1, numelem elemlen(i) = sqrt( (nodelocx(rightnode(i)) nodelocx(leftnode(i)))**2 + & (nodelocy(rightnode(i)) nodelocy(leftnode(i)))**2 ) write (6,63) (i, leftnode(i), rightnode(i), youngmod(i), area(i), elemlen(i), i=1, numelem) 63 format (" elem i j Modulus Area Length" / (3I5, 3E12.4))! Initialize global stiffness matrix ================================================= gstiff ( :, : ) = 0.0! or you can loop on matrix (slower)! Construct element stiffness ====================================================== do ielem = 1, numelem! atan or datan (real*8) returns radians temp = nodelocy(rightnode(ielem)) nodelocy(leftnode(ielem)) temp1 = nodelocx(rightnode(ielem)) nodelocx(leftnode(ielem)) angrad(ielem) = datan( temp / temp1 )! put in array, use again for stress
4 C:\files\classes\fem\program code\truss\truss solve\truss.f cosang = cos(angrad(ielem)) sinang = sin(angrad(ielem)) temp = area(ielem) * youngmod(ielem) / elemlen(ielem) elemstiff(1,1) = temp*cosang**2 elemstiff(1,2) = temp*cosang*sinang elemstiff(1,3) = elemstiff(1,1) elemstiff(1,4) = elemstiff(1,2) elemstiff(2,1) = elemstiff(1,2) elemstiff(2,2) = temp*sinang**2 elemstiff(2,3) = elemstiff(2,1) elemstiff(2,4) = elemstiff(2,2) elemstiff(3,1) = elemstiff(1,1) elemstiff(3,2) = elemstiff(1,2) elemstiff(3,3) = elemstiff(1,3) elemstiff(3,4) = elemstiff(1,4) elemstiff(4,1) = elemstiff(2,1) elemstiff(4,2) = elemstiff(2,2) elemstiff(4,3) = elemstiff(2,3) elemstiff(4,4) = elemstiff(2,4)!write (*,71) ( (elemstiff(n, i), i=1, 4), n=1, 4) 71 format ((4F9.3))! put element stiffness into global matrix i = leftnode(ielem) j = rightnode(ielem) do n = 1, 4 do k = 1, 4!! do each pair differently if (n < 3) irow = 2*(i 1) + n if (n >= 3) irow = 2*(j 1) + n 2 if (k < 3) icol = 2*(i 1) + k if (k >= 3) icol = 2*(j 1) + k 2 gstiff(irow, icol) = gstiff(irow, icol) + elemstiff(n, k)
5 C:\files\classes\fem\program code\truss\truss solve\truss.f gstifforig = gstiff!! save orginal for later to get all forces at nodes!! testing output, small matix only!write (6, 65) ( (gstifforig(n, i), i=1, dof), n=1, dof) 65 format (/"Stiffness" / (8E9.2) )! Apply displacement boundary conditions ================================================! set fixed node direction to zero (row and column) except when i=j, then 1.0 do inode = 1, numnode! ZERO X deflection set row and column to zero except diagonal term to 1 if (noderest(inode) == 1.OR. noderest(inode) == 3) then! x dir fixed irow = 2*(iNode 1) + 1 do icol = 1, dof if (icol /= irow) then gstiff(irow, icol) = 0.0 gstiff(icol, irow) = 0.0 else gstiff(irow, icol) = 1.0! can be anything except 0 load(2*(inode 1)+1) = 0.0! load is also set to zero (when solving, gives defl = 0)! ZERO Y deflection set row and column to zero except diagonal term to 1 if (noderest(inode) == 2.OR. noderest(inode) == 3) then! y dir fixed irow = 2*(iNode 1) + 2 do icol = 1, dof if (icol /= irow) then gstiff(irow, icol) = 0.0 gstiff(icol, irow) = 0.0 else gstiff(irow, icol) = 1.0! can be anything except 0 load(2*(inode 1)+2) = 0.0!! Testing output!write (6, 65) ( (gstiff(n, i), i=1, dof), n=1, dof)
6 C:\files\classes\fem\program code\truss\truss solve\truss.f !write (6, 72) ( load(i), i=1, dof) 72 format (/"force" / (8F8.0))!================== Solve using Cholesky ======================================! Find all displacements, known displaces for zero still fine! Cholesky is a matrix factor method, decomposition is stored in lower triangle! RHS not effected, solution is stored in disp array!! MISSING SOLVER CODE! debugging console output!write (*, 74) ( disp(i), i=1, dof) 74 format (/"displacement" / (8F8.3))!! get all forces (known and unknown), most will be zero do i = 1, dof load(i) = 0.0 do j = 1, dof load(i) = load(i) + gstifforig(i, j) * disp(j)! Find stress in each member do ielem = 1, numelem cosang = cos(angrad(ielem)) sinang = sin(angrad(ielem)) i = leftnode(ielem) j = rightnode(ielem) temp = youngmod(ielem) / elemlen(ielem) stress(ielem) = temp * ( cosang*disp(2*(i 1)+1) sinang*disp(2*(i 1)+2) & + cosang*disp(2*(j 1)+1) + sinang*disp(2*(j 1)+2) ) CALL CPU_TIME(t2) write (*,"('elasped time:', E14.8)") (t2 t1)
7 C:\files\classes\fem\program code\truss\truss solve\truss.f ! Results Output to file " Results " "Global Node Displacement and Force" write (6,64) (i, load(2*(i 1)+1), disp(2*(i 1)+1), load(2*(i 1)+2), disp(2*(i 1)+2), i = 1, numnode) 64 format (" Node X Load X Disp Y Load Y Disp " / (I5, 4E12.4))! Results Output to file "Element Axial Load and Stress" write (6,75) (ielem, stress(ielem)*area(ielem), stress(ielem), ielem = 1, numelem) 75 format (" Element Load Stress " / (I7, 1E14.4, 1E14.4)) write (*,*) 'truss End' end program truss
Numerical Methods for PDEs : Video 9: 2D Finite Difference February 14, Equations / 29
22.520 Numerical Methods for PDEs Video 9 2D Finite Difference Equations February 4, 205 22.520 Numerical Methods for PDEs Video 9 2D Finite Difference February 4, Equations 205 / 29 Thought Experiment
More informationStatics of the truss with force and temperature load - test problem Nr 1
Statics of the truss with force and temperature load - test problem Nr E := GPa - Young modulus for truss material - steel α t := - - thermal expanssion coeficient - steel D := cm - Cross section (pipe)
More informationIsoparametric Constant Strain Triangle, CST
function Plane_Stress_T3_RS %... % Plane stress displacements for constant strain triangle, T3 % ISOPARAMETRIC VERSION %... % Set given constants n_g = 2 ; % number of DOF per node n_n = 3 ; % number of
More informationMatlab Plane Stress Example
Matlab Plane Stress Example (Draft 2, April 9, 2007) Introduction Here the Matlab closed form element matrices for the T3 element (3 node triangle, constant stress) is illustrated for a square plate, 2
More informationFinite Element Analysis Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology Madras. Module - 01 Lecture - 15
Finite Element Analysis Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology Madras Module - 01 Lecture - 15 In the last class we were looking at this 3-D space frames; let me summarize
More informationc1=e*a/l; c2=12*e*i/(l^3); c3=6*e*i/(l^2); c4=2*e*i/l; c=cos(theta); s=sin(theta); K=[]; Verification
UNIVERSITY OF ILLINOIS AT URBANA-CHAMPAIGN Civil and Environmental Engineering CEE 361 - Methods of Structural Analysis Fall Semester, 2 Problem 1: Writing the functions Laboratory #6 Solutions There were
More informationCONTENT. A. Plane truss (2D) B. Space truss (3D) C. Observed steps of running programming... 02
CONTENT A. Plane truss (2D)... 01 B. Space truss (3D)... 01 C. Observed steps of running programming... 02 D. Details... 06 1. Part 1: Algorithm and program for truss structures by fem... 06 1.1 : Algorithm...
More information2D-truss-removed-sections.f90
Page:1/11 program problem2_1! this program solves problem 2.1 page 37! it reads in the element connectivity stored in the file elements.dat! and the node coordinates stored in file nodes.dat integer ::
More information7. Procedures and Structured Programming
7. Procedures and Structured Programming ONE BIG PROGRAM external procedure: separated small and reusable program units to conduct individual subtasks smaller main program Each program unit can be debugged
More informationFortran Programming: Plane Framework
- 15 Fortran Programming: Plane Framework 15 1 hapter 15: FORTRAN PROGRAMMING: PLANE FRAMEWORK 15 2 15.1 INTRODUTION The GFRAME2 program analyzes general plane framework structures under static loads.
More informationThe Application of EXCEL in Teaching Finite Element Analysis to Final Year Engineering Students.
The Application of EXCEL in Teaching Finite Element Analysis to Final Year Engineering Students. Kian Teh and Laurie Morgan Curtin University of Technology Abstract. Many commercial programs exist for
More informationPractice Reading for Loops
ME 350 Lab Exercise 3 Fall 07 for loops, fprintf, if constructs Practice Reading for Loops For each of the following code snippets, fill out the table to the right with the values displayed when the code
More informationSupplementary information
Modern Structural Analysis - Introduction to Modelling Supplementary information Chapter 3 Section 5.10 Equivalent beam for parallel chord trusses The cross references in the form n.m or n.m.p are to sub-sections
More informationStep by step set of instructions to accomplish a task or solve a problem
Step by step set of instructions to accomplish a task or solve a problem Algorithm to sum a list of numbers: Start a Sum at 0 For each number in the list: Add the current sum to the next number Make the
More informationRobot Inverse Kinematics Asanga Ratnaweera Department of Mechanical Engieering
PR 5 Robot Dynamics & Control /8/7 PR 5: Robot Dynamics & Control Robot Inverse Kinematics Asanga Ratnaweera Department of Mechanical Engieering The Inverse Kinematics The determination of all possible
More informationComputational Methods of Scientific Programming. Lecturers Thomas A Herring Chris Hill
12.010 Computational Methods of Scientific Programming Lecturers Thomas A Herring Chris Hill Review of Lecture 2 Examined computer hardware Computer basics and the main features of programs Program design:
More informationE (N/mm^2 ) and A ( mm^2)
E (N/mm^2 ) and A ( mm^2) File:truss.m truss.m LINEAR STATIC ANALYSIS OF A TRUSS STRUCTURE clc Clear screen clear Clear all variables in memory Make these variables global so they can be shared by other
More informationWinFElt Manual. Problem description
WinFElt Manual Problem description The problem description section is used to define the problem title and the number of nodes and elements in the problem. The problem description section is the only section
More informationIntroduction. Co-rotational Concept
1 2D Co-rotational Truss Formulation by Louie L. Yaw Walla Walla University April 23, 29 key words: geometrically nonlinear analysis, 2d co-rotational truss, corotational truss, variationally consistent,
More informationComputational Methods of Scientific Programming. Lecturers Thomas A Herring Chris Hill
12.010 Computational Methods of Scientific Programming Lecturers Thomas A Herring Chris Hill Review of last lecture Start examining the FORTRAN language Development of the language Philosophy of language:
More informationThis is NOT a truss, this is a frame, consisting of beam elements. This changes several things
CES 44 - Stress Analysis Spring 999 Ex. #, the following -D frame is to be analyzed using Sstan (read the online stan intro first, and Ch-6 in Hoit) 5 k 9 E= 9000 ksi 8 I= 600 in*in*in*in 5 A= 0 in*in
More informationThe report problem for which the FEA code was written is shown below. The Matlab code written to solve this problem is shown on the following pages.
The report problem for which the FEA code was written is shown below. The Matlab code written to solve this problem is shown on the following pages. %TrussGEP Script %Solves the Generalized Eigenvalue
More information(Type your answer in radians. Round to the nearest hundredth as needed.)
1. Find the exact value of the following expression within the interval (Simplify your answer. Type an exact answer, using as needed. Use integers or fractions for any numbers in the expression. Type N
More informationTrussMaster: Educational Truss Analysis Software Version 3 User Manual Dublin Institute of Technology, Ireland
TrussMaster: Educational Truss Analysis Software Version 3 User Manual Dublin Institute of Technology, Ireland 1 Acknowledgements This program has been developed over 10 years from 2001 by Dr Colin Caprani.
More informationLecture #9 Matrix methods
Lecture #9 Matrix methods METHODS TO SOLVE INDETERMINATE PROBLEM Small degree of statical indeterminacy Force method Displacement methods Displacement method in matrix formulation Large degree of statical
More informationFigure 6.1: Truss topology optimization diagram.
6 Implementation 6.1 Outline This chapter shows the implementation details to optimize the truss, obtained in the ground structure approach, according to the formulation presented in previous chapters.
More informationCITY AND GUILDS 9210 UNIT 135 MECHANICS OF SOLIDS Level 6 TUTORIAL 15 - FINITE ELEMENT ANALYSIS - PART 1
Outcome 1 The learner can: CITY AND GUILDS 9210 UNIT 135 MECHANICS OF SOLIDS Level 6 TUTORIAL 15 - FINITE ELEMENT ANALYSIS - PART 1 Calculate stresses, strain and deflections in a range of components under
More informationAllocating Storage for 1-Dimensional Arrays
Allocating Storage for 1-Dimensional Arrays Recall that if we know beforehand what size we want an array to be, then we allocate storage in the declaration statement, e.g., real, dimension (100 ) :: temperatures
More information1. Carlos A. Felippa, Introduction to Finite Element Methods,
Chapter Finite Element Methods In this chapter we will consider how one can model the deformation of solid objects under the influence of external (and possibly internal) forces. As we shall see, the coupled
More informationIntroduction to Programming with Fortran 90
Introduction to Programming with Fortran 90 p. 1/?? Introduction to Programming with Fortran 90 Array Concepts Nick Maclaren Computing Service nmm1@cam.ac.uk, ext. 34761 November 2007 Introduction to Programming
More informationTorsion Mobile App for Engineering Education Using a High Performance Computer (HPC) Cluster
Abstract Torsion Mobile App for Engineering Education Using a High Performance Computer (HPC) Cluster Kurt Gramoll Hughes Professor Aerospace and Mechanical Engineering University of Oklahoma, Norman,
More informationData Structures and Algorithms(12)
Ming Zhang "Data s and Algorithms" Data s and Algorithms(12) Instructor: Ming Zhang Textbook Authors: Ming Zhang, Tengjiao Wang and Haiyan Zhao Higher Education Press, 2008.6 (the "Eleventh Five-Year"
More informationBeams. Lesson Objectives:
Beams Lesson Objectives: 1) Derive the member local stiffness values for two-dimensional beam members. 2) Assemble the local stiffness matrix into global coordinates. 3) Assemble the structural stiffness
More informationVariable Definition and Statement Suppression You can create your own variables, and assign them values using = >> a = a = 3.
MATLAB Introduction Accessing Matlab... Matlab Interface... The Basics... 2 Variable Definition and Statement Suppression... 2 Keyboard Shortcuts... More Common Functions... 4 Vectors and Matrices... 4
More informationDocumentation for LISP in BASIC
Documentation for LISP in BASIC The software and the documentation are both Copyright 2008 Arthur Nunes-Harwitt LISP in BASIC is a LISP interpreter for a Scheme-like dialect of LISP, which happens to have
More informationCOMPUTER OPTIMIZATION
COMPUTER OPTIMIZATION Storage Optimization: Since Normal Matrix is a symmetric matrix, store only half of it in a vector matrix and develop a indexing scheme to map the upper or lower half to the vector.
More informationAssignment 2 Simulation and modeling, Spring 2010
Assignment 2 Simulation and modeling, Spring 2010 1 Background The consultant company B&W (Besser & Wisser) AB has decided to develop an in-house FEcode. With this code the company intends to offer their
More information01. Function Description and Limitation 02. Fortran File Preparation 03. DLL File Preparation 04. Using the USSR material model in midas FEA
midas FEA User Supplied Subroutine User Manual 01. Function Description and Limitation 02. Fortran File Preparation 03. DLL File Preparation 04. Using the USSR material model in midas FEA MIDAS Information
More informationIntroduction to GNU-Octave
Introduction to GNU-Octave Dr. K.R. Chowdhary, Professor & Campus Director, JIETCOE JIET College of Engineering Email: kr.chowdhary@jietjodhpur.ac.in Web-Page: http://www.krchowdhary.com July 11, 2016
More informationProgramming 1. Script files. help cd Example:
Programming Until now we worked with Matlab interactively, executing simple statements line by line, often reentering the same sequences of commands. Alternatively, we can store the Matlab input commands
More informationBasic matrix math in R
1 Basic matrix math in R This chapter reviews the basic matrix math operations that you will need to understand the course material and how to do these operations in R. 1.1 Creating matrices in R Create
More information7 Control Structures, Logical Statements
7 Control Structures, Logical Statements 7.1 Logical Statements 1. Logical (true or false) statements comparing scalars or matrices can be evaluated in MATLAB. Two matrices of the same size may be compared,
More informationPractice Test - Chapter 6
1. Write each system of equations in triangular form using Gaussian elimination. Then solve the system. Align the variables on the left side of the equal sign. Eliminate the x-term from the 2nd equation.
More informationAC : DEVELOPMENT AND IMPLEMENTATION OF A HIGH PERFORMANCE COMPUTER (HPC) CLUSTER FOR ENGINEERING EDUCATION SIMULATIONS
AC 2012-3853: DEVELOPMENT AND IMPLEMENTATION OF A HIGH PERFORMANCE COMPUTER (HPC) CLUSTER FOR ENGINEERING EDUCATION SIMULATIONS Dr. Kurt C. Gramoll, University of Oklahoma Kurt Gramoll is Hughes Professor
More information(Refer Slide Time 01:41 min)
Programming and Data Structure Dr. P.P.Chakraborty Department of Computer Science and Engineering Indian Institute of Technology, Kharagpur Lecture # 03 C Programming - II We shall continue our study of
More informationMATLAB Lesson I. Chiara Lelli. October 2, Politecnico di Milano
MATLAB Lesson I Chiara Lelli Politecnico di Milano October 2, 2012 MATLAB MATLAB (MATrix LABoratory) is an interactive software system for: scientific computing statistical analysis vector and matrix computations
More information0_PreCNotes17 18.notebook May 16, Chapter 12
Chapter 12 Notes BASIC MATRIX OPERATIONS Matrix (plural: Matrices) an n x m array of elements element a ij Example 1 a 21 = a 13 = Multiply Matrix by a Scalar Distribute scalar to all elements Addition
More informationLecture 2 FORTRAN Basics. Lubna Ahmed
Lecture 2 FORTRAN Basics Lubna Ahmed 1 Fortran basics Data types Constants Variables Identifiers Arithmetic expression Intrinsic functions Input-output 2 Program layout PROGRAM program name IMPLICIT NONE
More informationLS-DYNA s Linear Solver Development Phase 2: Linear Solution Sequence
LS-DYNA s Linear Solver Development Phase 2: Linear Solution Sequence Allen T. Li 1, Zhe Cui 2, Yun Huang 2 1 Ford Motor Company 2 Livermore Software Technology Corporation Abstract This paper continues
More informationAN INTRODUCTION TO FORTRAN 90 LECTURE 2. Consider the following system of linear equations: a x + a x + a x = b
AN INTRODUCTION TO FORTRAN 90 LECTURE 2 1. Fortran 90. Arrays, functions and subroutines. 2. Scientific plotting. Gnuplot 1 Each coefficient and variable is a scalar. Lengthy and cumbersome! Program Scalar
More informationIsothermal Batch Reactor Modeling
Instructions for use of the tutorial: Download the compressed file Example1.zip and store it on a folder of your choice on your desktop, or in a location where you have rights to read and write. Open the
More informationGoals for This Lecture:
Goals for This Lecture: Learn about multi-dimensional (rank > 1) arrays Learn about multi-dimensional array storage Learn about the RESHAPE function Learn about allocatable arrays & the ALLOCATE and DEALLOCATE
More informationTrig Functions, Equations & Identities May a. [2 marks] Let. For what values of x does Markscheme (M1)
Trig Functions, Equations & Identities May 2008-2014 1a. Let. For what values of x does () 1b. [5 marks] not exist? Simplify the expression. EITHER OR [5 marks] 2a. 1 In the triangle ABC,, AB = BC + 1.
More informationRevised Sheet Metal Simulation, J.E. Akin, Rice University
Revised Sheet Metal Simulation, J.E. Akin, Rice University A SolidWorks simulation tutorial is just intended to illustrate where to find various icons that you would need in a real engineering analysis.
More informationSubprograms. FORTRAN 77 Chapter 5. Subprograms. Subprograms. Subprograms. Function Subprograms 1/5/2014. Satish Chandra.
FORTRAN 77 Chapter 5 Satish Chandra satish0402@gmail.com When a programs is more than a few hundred lines long, it gets hard to follow. Fortran codes that solve real research problems often have tens of
More informationChallenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure
Challenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure In the final year of his engineering degree course a student was introduced to finite element analysis and conducted an assessment
More informationixcube 4-10 Brief introduction for membrane and cable systems.
ixcube 4-10 Brief introduction for membrane and cable systems. ixcube is the evolution of 20 years of R&D in the field of membrane structures so it takes a while to understand the basic features. You must
More informationChapter 01 Arrays Prepared By: Dr. Murad Magableh 2013
Chapter 01 Arrays Prepared By: Dr. Murad Magableh 2013 One Dimensional Q1: Write a program that declares two arrays of integers and fills them from the user. Then exchanges their values and display the
More informationAMath 483/583 Lecture 8
AMath 483/583 Lecture 8 This lecture: Fortran subroutines and functions Arrays Dynamic memory Reading: class notes: Fortran Arrays class notes: Fortran Subroutines and Functions class notes: gfortran flags
More informationMAT 275 Laboratory 2 Matrix Computations and Programming in MATLAB
MATLAB sessions: Laboratory MAT 75 Laboratory Matrix Computations and Programming in MATLAB In this laboratory session we will learn how to. Create and manipulate matrices and vectors.. Write simple programs
More informationCS1073 Exam 3, Fall 2009 page 1. Important test instructions code fragments
CS1073 Exam 3, Fall 2009 page 1 Name (please print): Important test instructions code fragments Throughout this exam, you will be asked to write a code fragment based on certain assumptions (eg., assume
More informationSami Ilvonen Pekka Manninen. Introduction to High-Performance Computing with Fortran. September 19 20, 2016 CSC IT Center for Science Ltd, Espoo
Sami Ilvonen Pekka Manninen Introduction to High-Performance Computing with Fortran September 19 20, 2016 CSC IT Center for Science Ltd, Espoo All material (C) 2009-2016 by CSC IT Center for Science Ltd.
More informationA very brief Matlab introduction
A very brief Matlab introduction Siniša Krajnović January 24, 2006 This is a very brief introduction to Matlab and its purpose is only to introduce students of the CFD course into Matlab. After reading
More informationSelf-study session 1, Discrete mathematics
Self-study session 1, Discrete mathematics First year mathematics for the technology and science programmes Aalborg University In this self-study session we are working with time complexity. Space complexity
More informationDEPARTMENT - Mathematics. Coding: N Number. A Algebra. G&M Geometry and Measure. S Statistics. P - Probability. R&P Ratio and Proportion
DEPARTMENT - Mathematics Coding: N Number A Algebra G&M Geometry and Measure S Statistics P - Probability R&P Ratio and Proportion YEAR 7 YEAR 8 N1 Integers A 1 Simplifying G&M1 2D Shapes N2 Decimals S1
More informationASSIGNMENT 1 INTRODUCTION TO CAD
Computer Aided Design(2161903) ASSIGNMENT 1 INTRODUCTION TO CAD Theory 1. Discuss the reasons for implementing a CAD system. 2. Define computer aided design. Compare computer aided design and conventional
More informationMath12 Pre-Calc Review - Trig
Math1 Pre-Calc Review - Trig Multiple Choice Identify the choice that best completes the statement or answers the question. 1. Which of the following angles, in degrees, is coterminal with, but not equal
More informationCE371 Structural Analysis II Lecture 5:
CE371 Structural Analysis II Lecture 5: 15.1 15.4 15.1) Preliminary Remarks 15.2) Beam-Member Stiffness Matrix 15.3) Beam-Structure Stiffness Matrix 15.4) Application of the Stiffness Matrix. 15.1) Preliminary
More information3D Coordinate Transformation Calculations. Space Truss Member
3D oordinate Transformation alculations Transformation of the element stiffness equations for a space frame member from the local to the global coordinate system can be accomplished as the product of three
More informationME1107 Computing Y Yan.
ME1107 Computing 1 2008-2009 Y Yan http://www.staff.city.ac.uk/~ensyy About Fortran Fortran Formula Translation High level computer language Basic, Fortran, C, C++, Java, C#, (Matlab) What do we learn?
More informationData Types and Basic Calculation
Data Types and Basic Calculation Intrinsic Data Types Fortran supports five intrinsic data types: 1. INTEGER for exact whole numbers e.g., 1, 100, 534, -18, -654321, etc. 2. REAL for approximate, fractional
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 36
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 36 In last class, we have derived element equations for two d elasticity problems
More informationIntroduction to Programming in C Department of Computer Science and Engineering. Lecture No. #16 Loops: Matrix Using Nested for Loop
Introduction to Programming in C Department of Computer Science and Engineering Lecture No. #16 Loops: Matrix Using Nested for Loop In this section, we will use the, for loop to code of the matrix problem.
More information17 USING THE EDITOR AND CREATING PROGRAMS AND FUNCTIONS
17 USING THE EDITOR AND CREATING PROGRAMS AND FUNCTIONS % Programs are kept in an m-file which is a series of commands kept in the file that is executed from MATLAB by typing the program (file) name from
More informationOrdinary Differential Equation Solver Language (ODESL) Reference Manual
Ordinary Differential Equation Solver Language (ODESL) Reference Manual Rui Chen 11/03/2010 1. Introduction ODESL is a computer language specifically designed to solve ordinary differential equations (ODE
More informationFinal Exam: Precalculus
Final Exam: Precalculus Apr. 17, 2018 ANSWERS Without Notes or Calculators Version A 1. Consider the unit circle: a. Angle in degrees: What is the angle in radians? What are the coordinates? b. Coordinates:
More informationElasto-Plastic Deformation of a Truss Structure
WORKSHOP PROBLEM 8 Elasto-Plastic Deformation of a Truss Structure Objectives: Demonstrate the use of elastic-plastic material properties. Create an enforced displacement on the model. Run an MSC/NASTRAN
More informationPhysics 326 Matlab Primer. A Matlab Primer. See the file basics.m, which contains much of the following.
A Matlab Primer Here is how the Matlab workspace looks on my laptop, which is running Windows Vista. Note the presence of the Command Window in the center of the display. You ll want to create a folder
More informationRepetition Structures Chapter 9
Sum of the terms Repetition Structures Chapter 9 1 Value of the Alternating Harmonic Series 0.9 0.8 0.7 0.6 0.5 10 0 10 1 10 2 10 3 Number of terms Objectives After studying this chapter you should be
More informationScilab Programming. The open source platform for numerical computation. Satish Annigeri Ph.D.
Scilab Programming The open source platform for numerical computation Satish Annigeri Ph.D. Professor, Civil Engineering Department B.V.B. College of Engineering & Technology Hubli 580 031 satish@bvb.edu
More informationROSE-HULMAN INSTITUTE OF TECHNOLOGY
EXAM 2 WRITTEN PORTION NAME SECTION NUMBER CAMPUS MAILBOX NUMBER EMAIL ADDRESS @rose-hulman.edu Written Portion / 40 Computer Portion / 60 Total / 100 USE MATLAB SYNTAX FOR ALL PROGRAMS AND COMMANDS YOU
More informationMPI Version of the Stommel Code with One and Two Dimensional Decomposition
MPI Version of the Stommel Code with One and Two Dimensional Decomposition Timothy H. Kaiser, Ph.D. tkaiser@sdsc.edu 1 Overview We will choose one of the two dimensions and subdivide the domain to allow
More informationChapter 1 Introduction
Chapter 1 Introduction GTU Paper Analysis (New Syllabus) Sr. No. Questions 26/10/16 11/05/16 09/05/16 08/12/15 Theory 1. What is graphic standard? Explain different CAD standards. 2. Write Bresenham s
More informationStarting Matlab. MATLAB Laboratory 09/09/10 Lecture. Command Window. Drives/Directories. Go to.
Starting Matlab Go to MATLAB Laboratory 09/09/10 Lecture Lisa A. Oberbroeckling Loyola University Maryland loberbroeckling@loyola.edu http://ctx.loyola.edu and login with your Loyola name and password...
More informationChapter 8. Computer Implementation of 2D, Incompressible N-S Solver
Chapter 8 Computer Implementation of 2D, Incompressible N-S Solver 8.1 MATLAB Code for 2D, Incompressible N-S Solver (steadynavierstokes2d.m) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
More informationPreliminary remarks. Preliminary remarks. Preliminary remarks. CHAPTER 7 BM Analysis using Stiffness Method
CHAPTER 7 BM Analysis using Stiffness Method Objectives เข าใจว ธ ของ stiffness method ก บ โครงสร างประเภทคาน Member & node identification In general each element must be free from load & have a prismatic
More informationModeling Skills Stress Analysis J.E. Akin, Rice University, Mech 417
Introduction Modeling Skills Stress Analysis J.E. Akin, Rice University, Mech 417 Most finite element analysis tasks involve utilizing commercial software, for which you do not have the source code. Thus,
More informationLesson: An introduction to TFEM
Lesson: An introduction to TFEM January 10, 2018 Contents 1 Introduction 2 2 A first run 3 3 Mesh generation 3 4 Mesh reading 4 5 Problem definition 6 6 Building the system matrix and right hand side and
More informationINTRODUCTION TO FORTRAN PART II
INTRODUCTION TO FORTRAN PART II Prasenjit Ghosh An example: The Fibonacci Sequence The Fibonacci sequence consists of an infinite set of integer nos. which satisfy the following recurrence relation x n
More informationWhat is MATLAB and howtostart it up?
MAT rix LABoratory What is MATLAB and howtostart it up? Object-oriented high-level interactive software package for scientific and engineering numerical computations Enables easy manipulation of matrix
More informationHomework 5: Transformations in geometry
Math 21b: Linear Algebra Spring 2018 Homework 5: Transformations in geometry This homework is due on Wednesday, February 7, respectively on Thursday February 8, 2018. 1 a) Find the reflection matrix at
More informationA quick guide to Fortran
A quick guide to Fortran Sergiy Bubin Department of Physics Nazarbayev University History of Fortran One of the oldest general purpose high-level computer languages First developed in 1957 at IBM in the
More informationFORTRAN Basis. PROGRAM LAYOUT PROGRAM program name IMPLICIT NONE [declaration statements] [executable statements] END PROGRAM [program name]
PROGRAM LAYOUT PROGRAM program name IMPLICIT NONE [declaration statements] [executable statements] END PROGRAM [program name] Content in [] is optional. Example:- PROGRAM FIRST_PROGRAM IMPLICIT NONE PRINT*,
More informationCE2351-STRUCTURAL ANALYSIS II
CE2351-STRUCTURAL ANALYSIS II QUESTION BANK UNIT-I FLEXIBILITY METHOD PART-A 1. What are determinate structures? 2. What is meant by indeterminate structures? 3. What are the conditions of equilibrium?
More informationLecture 1 arithmetic and functions
Lecture 1 arithmetic and functions MATH 190 WEBSITE: www.math.hawaii.edu/190 Open MATH 190 in a web browser. Read and accept the Terms of Acceptable Lab Use. Open Lecture 1. PREREQUISITE: You must have
More informationCurriculum Map: Mathematics
Curriculum Map: Mathematics Course: Honors Advanced Precalculus and Trigonometry Grade(s): 11-12 Unit 1: Functions and Their Graphs This chapter will develop a more complete, thorough understanding of
More informationFEA of Composites Classical Lamination Theory Example 1
FEA of Composites Classical Lamination Theory Example 1 22.514 Instructor: Professor James Sherwood Author: Dimitri Soteropoulos Revised by Jacob Wardell Problem Description: A four layer [0/90] s graphite-epoxy
More informationNO CALCULATOR ALLOWED!!
CPSC 203 500 EXAM TWO Fall 2005 NO CALCULATOR ALLOWED!! Full Name (Please Print): UIN: Score Possible Points Prog Points Part One 33 pts Part Two 30 pts Part Three 20 pts Part Four 25 pts Total 108 pts
More information3+2 3*2 3/2 3^2 3**2 In matlab, use ^ or ** for exponentiation. In fortran, use only ** not ^ VARIABLES LECTURE 1: ARITHMETIC AND FUNCTIONS
LECTURE 1: ARITHMETIC AND FUNCTIONS MATH 190 WEBSITE: www.math.hawaii.edu/ gautier/190.html PREREQUISITE: You must have taken or be taking Calculus I concurrently. If not taken here, specify the college
More informationIntroduction to Finite Element Analysis using ANSYS
Introduction to Finite Element Analysis using ANSYS Sasi Kumar Tippabhotla PhD Candidate Xtreme Photovoltaics (XPV) Lab EPD, SUTD Disclaimer: The material and simulations (using Ansys student version)
More information