Harry Millwater, Ph.D. David Wagner, MS Jose Garza, MS Andrew Baines, BS Kayla Lovelady, BS Carolina Dubinsky, MS Thomas Ross, MS
|
|
- Claude Lucas
- 5 years ago
- Views:
Transcription
1 Complex Variable Finite Element Methods for Fracture Mechanics Analysis Harry Millwater, Ph.D. David Wagner, MS Jose Garza, MS Andrew Baines, BS Kayla Lovelady, BS Carolina Dubinsky, MS Thomas Ross, MS Sivirt Student Seminary, April 2013
2 Overview Development of the complex variable finite element method (ZFEM) Derivatives wrt shape, material prop, load and others Implemented into the Abaqus finite element code Verification against known solutions Applications to fracture mechanics Future efforts
3 Finite Difference Method Im Perturb along the imaginary axis Finite Differencing Re Forward Differencing h df (x o ) dx» f (x o + h) - f (x o ) h Determining h is problematic
4 Complex Taylor Series Expansion Perturb along the imaginary axis Im Subset of Fourier Differentiation h Re F(x + ih) = F(x) + ih df dx - h 2 2! - ih 3 3! d 3 F dx 3 + d 2 F dx 2 df (x o ) dx» Im( f (x o + ih)) h h can be very small ~ 10-30
5 Example: dexp[x] dx = Exp[x] Exp[(x + h)] - Exp[x] h Finite Difference Im[Exp[x +ih]]/h = Exp[x]sin(h)/h CTSE No subtraction of nearly equal numbers h can be as small as desired Exact derivative can be recovered as h becomes small
6 Finite element Implementation ì Pü í ý î 0 þ = EI é 12 6L ù ì L 3 ê ë 6L 4L 2 ú d ü í ý û î fþ ì Pü í ý î 0 þ = EI é 12 6(L + ih) ù ì (L + ih) 3 ê ë 6(L + ih) 4(L + ih) 2 ú d ü í ý û î fþ d = ì PL3 3EI - PLh 2 í î EI ü ý þ + i ì PL 2 h í î EI - Ph 3 3EI ü ý þ d L = Im[d]/h
7 Finite Element Analysis Shape sensitivities difficult to compute using finite difference method Mesh must be perturbed enough to see the perturbation but not too large to distort elements h must be kept small
8 Complex Variable Advantage Perturb nodal coordinates in imaginary axis only! Minimal mesh distortion issues!
9 Computational Issues Standard finite element EXCEPT complex nodal coordinates Stiffness matrix is now complex->complex solver required Apply an imaginary displacement to the nodal coordinates to represent a shape (domain) change Shape change is arbitrary Sensitivities for entire displacement/strain/stress available u a Imaginary displacement
10 Example: 1 st Der. of σ θ w.r.t. R i Analytic Solution CTSE UTSA Matlab code University (2011) , of doi: /j.finel Texas at San Antonio A. Voorhees, H.R. Millwater, R.L. Bagley, Complex Variable Methods for Shape Sensitivity of Finite Element Models, Finite Elem. Anal. Des., 47
11 L2 Norm Error controlled by Mesh L2 norm reduces as mesh is refined
12 Current Effort Implementation into Abaqus Abaqus user element implementation (uel) 6 dof/node (3 real, 3 imag) Abaqus cannot solve complex stiffness matrix, represent as real { } é Re P ê ê ë Im P { } ù = é Re K ú ê ú û ë ê Im K [ ] -Im[ K] [ ] Re[ K] ù é Re U ú ê û ú ê ë Im U n x n complex matrix solved as 2n x 2n real matrix { } { } ù ú ú û
13 Abaqus Example Plane Strain ½ P ½ P Real nodes (physical coordinates) Imaginary nodes (perturbation in physical coordinates Abaqus UEL (user element) Plane stress/strain Element has 8 nodes (4 real, 4 imag)
14 Abaqus UEL Example Code COMPLEX(KIND=r8),DIMENSION(MCRD,NNODE/2 ) :: zcoords COMPLEX(KIND=r8),DIMENSION(NNODE,NNODE)::zAMATRX!Complex Stiffness Matrix!Construct complex nodal coordinate array zcoords = CMPLX(COORDS(1:2,1:3), COORDS(1:2,NNODE/2+1:NNODE)*h) f1 = zcoords(1,2)*zcoords(2,3) - zcoords(1,3)*zcoords(2,2) f2 = zcoords(1,3)*zcoords(2,1) - zcoords(1,1)*zcoords(2,3) f3 = zcoords(1,1)*zcoords(2,2) - zcoords(1,2)*zcoords(2,1) b1 = zcoords(2,2) - zcoords(2,3) b2 = zcoords(2,3) - zcoords(2,1) b3 = zcoords(2,1) - zcoords(2,2) c1 = zcoords(1,3) - zcoords(1,2) c2 = zcoords(1,1) - zcoords(1,3) c3 = zcoords(1,2) - zcoords(1,1) Area = half * (f1 + f2 + f3) zamatrx = MATMUL(Bt, MATMUL(Cthick,zB_ip) ) * Area
15 Obtaining Derivatives Shape sensitivity input imaginary coordinates to represent shape change. (All Imag nodes not perturbed have a value of zero) ½ P ½ P Nodal Coordinates 1 (0,0) 8 2 (1,0) 3 (0,1) L 4 (1,1) 5 (0,0) 6 (0,0) 7 (0,1)h h = ( 0,1)h
16 Obtaining Derivatives Load sensitivity Apply perturbation in loading to IM nodes. ½ P ½ P Applied Loads 3 (0,1/2 P) 4 (0,1/2 P) 7 (0,1/2)h 8( 0,1/2)h Imaginary coordinates all zero h = 10-10
17 Obtaining Derivatives Material sensitivity Apply perturbation to constitutive matrix. ½ P ½ P E (1+n)(1-2n) é 1-n n 0 ù ê ú ê n 1-n 0 ú ë ê 0 0 1/2(1-2n) û ú E = E+ ih n =n + ih h = Imaginary coordinates all zero
18 Abaqus Verification Results Exact Solutions ABAQUS UEL Solutions Exact Solutions ABAQUS UEL Solutions u L = P EA = U L = 1 P 2 2 EA = u P = L EA = U P = PL EA = u L = U L = u P = U P = s L = 0 e L = 0 s P = 1 A = e P = 1 EA = s L = 0 e L = 0 s P = e P = u E = - PL E 2 A = U E = -1 P 2 L 2 E 2 A = u E = U E = s E = 0 e E = -P E 2 A = s E = 0 e E = u A = - PL EA = U A = -1 P 2 L 2 EA = u A = U A = s A = -P A = e A = -P EA = s A = e A =
19 Abaqus Verification Results SDVs Stress Values SDV 1 SDV 2 SDV 3 SDV 4 SDV 5 SDV 6 s xx s yy s xy s xx q s yy q s xy q SDVs SDV7 SDV8 SDV9 SDV1 0 SDV11 SDV1 2 Strain Values e xx e yy e xy e xx q e yy q e xy q SDVs SDV13 SDV14 Strain Energy Values U U q SDVs SDV1 5 SDV1 6 SDV1 7 SDV1 8 SDV1 9 SDV2 0 Displacement Values u x1 q u y1 q u x2 q u y2 q u x3 q u y3 q
20 Abaqus 2D Verification Analytical Abaqus s q r i
21 Applications to Fracture Mechanics
22 Weight Function Development Calculation of the partial derivative of the crack opening displacement with respect to crack length required m(a,x) = E u 2K A a Standard research approaches: m(x,a) = Assume 3-4 term approximations to du/da or approximate the weight fn. directly. 2 2p(a - x) 1+ M 1 (1- x /a)1/ 2 + M 2 (1- x /a) M n (1- x /a) n / 2 [ ] Use multiple reference solutions to solve for Mi
23 Perturbation of Crack Length Crack tip element can be perturbed in the imaginary domain no perturbation of real mesh Perturb a no. of elements around the crack tip Perturbation of Crack Length in Imaginary domain
24 Accuracy: Infinite Array of Cracks Crack Opening Displacement nodes on crack line u A a Accurate calculation of the crack opening displacement with respect to crack length all along the crack line
25 Infinite Array Weight Function Accuracy controlled by the user: high order expansions possible 7-term [WCTSE] 3-term [reference] 3-term [WCTSE]
26 ComparisontoDomainIntegral (K/K16) ComparisontoDomainIntegral (K/K16) Accuracy: Crack from a Hole HolewithaSingleCrack StresInensityFactor (R/W=0.5, H/W=2) HolewithaSingleCrack StresIntensityFactor (R/W=0.2, H/W=2) WCTSE ±0.15% WCTSE ±0.06% TensionRef.[Wu] TensionWCTSE Tension[Wu] TensionBEM BendingWCTSE Bending[Wu] BendingBEM CrackLength(a/B) 0.98 TensionRef.[Wu] TensionWCTSE Tension[Wu] TensionBEM BendingWCTSE Bending[Wu] BendingBEM CrackLength(a/B) WCTSE weight function consistently better than published weight fns. University Series Expansion of Method, Texas Engng at Fract San Mech 86 (2012), Antonio D. Wagner, and H. Millwater, 2D Weight Function Development using a Complex Taylor
27 Strain Energy Release Rate Double Cantilever Beam (DBC) Energy release rate (Exact) = Contours 8 4 2
28 Strain Energy Release Rate # of Contours Perturbation Method Double Cantilever Beam (DBC) Exact = 0.8 Exact = Abaqus J Integral Abaqus CTSE 20 Tip only Tip only Tip only , , Tip only , Contours Contours Contour + Quarterpoints 2 Crack tip+quarterpoints , ,
29 3D Applications u L
30 3D Applications s R 2
31 Implementation into Abaqus [ ] ===================== Real Mesh ======================= *NODE, NSET=rnodes 1, 0.0d0, 0.0d0 2, 2.0d0, 0.0d0 3, 1.2d0, 1.0d0 4, 0.2d0, 1.0d0 5, 1.0d0, 0.0d0 6, 1.6d0, 0.5d0 7, 0.6d0, 1.0d0 8, 0.1d0, 0.5d0 *ELEMENT, TYPE=CPE8, ELSET=relems 1, 1, 2, 3, 4, 5, 6, 7, 8 ===================== Material ======================== *SOLID SECTION, ELSET=relems, MATERIAL=aluminum 1d0 *MATERIAL, NAME=aluminum *ELASTIC 100d9, 0.3d0 ======== Node Sets, Element Sets, and Surfaces ======== *NSET, NSET=rbotedge 1, 5, 2 ================ Boundary Conditions ================== *BOUNDARY rbotedge, 2 1, 1 ==================== Load Step ======================== *STEP, NAME=Step-1 *STATIC, DIRECT 1d0, 1d0, 1d-05, 1d0 [ ] Abaqus Input File Changes Needed [ ] *USER ELEMENT, TYPE=U28, NODES=16, COORDINATES=2, UNSYMM, PROPERTIES=5, IPROPERTIES=2, VARIABLES=108 1, 2 ===================== Real Mesh ======================= *NODE, NSET=rnodes 1, 0.0d0, 0.0d0 2, 2.0d0, 0.0d0 3, 1.2d0, 1.0d0 4, 0.2d0, 1.0d0 5, 1.0d0, 0.0d0 6, 1.6d0, 0.5d0 7, 0.6d0, 1.0d0 8, 0.1d0, 0.5d0 ==================== Complex Mesh ===================== *NODE, NSET=inodes 101, 0.0d0, 0.0d0 102, 0.0d0, 0.0d0 103, 1.0d0, 0.0d0 104, 0.0d0, 0.0d0 105, 0.0d0, 0.0d0 106, 0.5d0, 0.0d0 107, 0.5d0, 0.0d0 108, 0.0d0, 0.0d0 *ELEMENT, TYPE=U28, ELSET=zelems 101, 1, 2, 3, 4, 5, 6, 7, 8, 101,102,103,104,105,106,107,108 ===================== Material ======================== *UEL PROPERTY, ELSET=zelems 100d9, 0.3d0, 1d-9, 1d0, 1d0, 3, 0 ======== Node Sets, Element Sets, and Surfaces ======== *NSET, NSET=rbotedge 1, 5, 2 *NSET, NSET=ibotedge 101,105,102 ================ Boundary Conditions ================== *BOUNDARY rbotedge, 2 1, 1 ibotedge, 2 101, 1 ==================== Load Step ======================== *STEP, NAME=Step-1, UNSYMM=YES *STATIC, DIRECT 1d0, 1d0, 1d-5, 1d0 [ ] Imaginary nodes UEL definition UEL properties
32 Multicomplex Mathematics Bi-complex numbers representation : Matrix representation of bi-complex numbers: Tri-complex numbers representation : 2) X «x 1 -x 2 -x 3 x 4 x 2 x 1 -x 4 -x 3 x 3 -x 4 x 1 -x 2 x 4 x 3 x 2 x 1 é ë ê ê ê ê ê ù û ú ú ú ú ú 1) Z «z 1 -z 2 z 2 z 1 é ë ê ê ù û ú ú 2)X = x 1 +i 1 x 2 +i 2 x 3 +i 1 i 2 x 4 ;x 1,, x 4 Î R, i 1 2 = i 2 2 = -1, i 1 2 i 2 2 = i 2 2 i 1 2 1) Z = z 1 +i 2 z 2 ; z 1, z 2 Î C 1, i 2 2 = -1 1,,, )]; ( ) [( )] ( ) [( 1) i i i R x x x i x i x i x i x i x i x i x X ( 2 ) C ( 3 ) C 1,, ); ( ) ( 2) i C z z z i z i z i z Z
33 Multicomplex Mathematics Example f ( x +i 1 h +i 2 h) = 2( x+i 1 h+i 2 h) 3 = 2x 3 +(6i 1 h+6i 2 h)x 2 +(12i 1 i 2 h 2-12h 2 )x +2i 1 3 h 3-6i 1 h 3-6i 2 h 3 +2i 2 3 h 3. f ( x +i 1 h +i 2 h) = 2x 3-12h 2 x+i 1 h( 6x 2-8h 2 )+i 2 h( 6x 2-8h 2 )+12i 12 h 2 x Re[ f ] = 2x 3 +O( h 2 ) f x» 1 h Im 1 f [ ] = 6x 2 +O h 2 ( ) 2 f x 2» 1 h 2 Im 12 f [ ] =12x +O h 2 ( )
34 Bicomplex Analysis Higher Order Sensitivities Thick walled cylinder s r r o s q r o 2 s r r o 2 2 s q r o 2 2 s r p r o 2 s q p r o
35 Progressive Fracture (2) Compute the curved crack path to the desired distance (Dd). (1) Perturb the crack tip in x and y. (3) Add curved crack faces, remesh, and repeat. Construct a 3 rd order Taylor series of the strain energy using tricomplex elements in front of the crack tip. Predict the crack path along the max energy release Progress the crack, remesh, repeat
36 Progressive Fracture (2) Compute the curved crack path to the desired distance (Dd). (1) Perturb the crack tip in x and y. (3) Add curved crack faces, remesh, and repeat. Construct a 3 rd order Taylor series of the strain energy using tricomplex elements. Predict the crack path along the max energy release Progress the crack, remesh, repeat
37 y/w Progressive Fracture 1 Crack Path through a Cruciform Specimen, P 2 /P 1 =1/2, C Dd/W=0.08, 12 Steps [Franc2D] Dd/W=0.04, 25 Steps [Franc2D] Dd/W=0.02, 50 Steps [Franc2D] Dd/W=0.01,100 Steps [Franc2D] Dd/W=0.24, 3 Steps [ZFEM] Dd/W=0.16, 6 Steps [ZFEM] Dd/W=0.08, 12 Steps [ZFEM] Dd/W=0.04, 24 Steps [ZFEM] x/w Similar results to Franc2D with far fewer steps
38 y/w y/w Progressive Fracture Crack Path toward a Hole, Three-point Bending, C 3 Dd/W=0.04, 8 Steps [Franc2D] Dd/W=0.02, 16 Steps [Franc2D] Dd/W=0.005, 66 Steps [Franc2D] Dd/W= 0.1, 3 Steps [ZFEM] Dd/W=0.08, 4 Steps [ZFEM] Dd/W=0.04, 8 Steps [ZFEM] Dd/W=0.02, 16 Steps [ZFEM] x/w Crack Path toward a Hole, Three-point Bending, C 3 Dd/W=0.04, 8 Steps [Franc2D] Dd/W=0.02, 16 Steps [Franc2D] Dd/W=0.005, 66 Steps [Franc2D] Dd/W= 0.1, 3 Steps [ZFEM] Dd/W=0.08, 4 Steps [ZFEM] Dd/W=0.04, 8 Steps [ZFEM] Dd/W=0.02, 16 Steps [ZFEM] x/w
39 ZFEM User Element Features Easy to Use Only a Single Input File Needed Single User Subroutine File No Abaqus Configuration Needed Analytic Extensions of Built-in Elements Identical Response of Real values 2D (6,8 noded), 3D (8, 15, 20) Modular Fortran 95 Implementation Extensible Very Good Performance abaqus job=mycomplexmodel user=zfem
40 High order sensitivities: Future Interests Multicomplex mathematics element library (complete for linear elastic statics) Extension to non-linear materials Plasticity enhancement (in progress) Composites, anisotropic materials 3D progressive fracture (in progress) Structural dynamics (in progress) Thermoelastic analysis Expanded element library Plates and shells
41 Acknowledgements Efficient Sensitivity Methods for Probabilistic Lifing and Engine Prognostics, Pat Golden, AFRL/RXLMN, Aug Sep Efficient Finite Element-based 3D Fracture Mechanics Crack Growth Analysis using Complex Variable Sensitivity Methods, DoD PETTT, Sep Aug Implementation of Complex Variable Finite Element Methods in Abaqus, DOD PETTT, Sep Aug Enhanced Fracture Mechanics Crack Growth Analysis using Complex Variable Sensitivity Methods, AFOSR (David Stargel), May
Preprocessor Geometry Properties )Nodes, Elements(, Material Properties Boundary Conditions(displacements, Forces )
در برنامه يك تدوين براي بعدي دو يك سازه محيط MATLAB Preprocessor Geometry Properties )Nodes, Elements(, Material Properties Boundary Conditions(displacements, Forces ) Definition of Stiffness Matrices
More informationBilinear Quadrilateral (Q4): CQUAD4 in GENESIS
Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms
More informationDue Monday, September 14 th, 12:00 midnight
Due Monday, September 14 th, 1: midnight This homework is considering the analysis of plane and space (3D) trusses as discussed in class. A list of MatLab programs that were discussed in class is provided
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationMethods of Analysis. Force or Flexibility Method
INTRODUCTION: The structural analysis is a mathematical process by which the response of a structure to specified loads is determined. This response is measured by determining the internal forces or stresses
More informationFinite Elements for Large Strains - A double mixed (M 2 ) Formulation
Finite Elements for Large Strains - A double mixed (M 2 ) Formulation Motivation Development of user friendly elements robustness simple treatment of incompressible materials complex geometries geometrical
More informationENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material
ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a
More information7. Hierarchical modeling examples
7. Hierarchical modeling examples The objective of this chapter is to apply the hierarchical modeling approach discussed in Chapter 1 to three selected problems using the mathematical models studied in
More informationA modified quarter point element for fracture analysis of cracks
ndian Journal of Engineering & Materials Sciences Vol. 14, February 007, pp. 31-38 A modified quarter point element for fracture analysis of cracks Sayantan Paul & B N Rao* Structural Engineering Division,
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More information46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference April 2005 Austin, Texas
th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference - April, Austin, Texas AIAA - AIAA - Bi-stable Cylindrical Space Frames H Ye and S Pellegrino University of Cambridge, Cambridge,
More informationAn exact reanalysis algorithm using incremental Cholesky factorization and its application to crack growth modeling
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 01; 91:158 14 Published online 5 June 01 in Wiley Online Library (wileyonlinelibrary.com). DOI: 10.100/nme.4 SHORT
More informationLecture 7: The Beam Element Equations.
4.1 Beam Stiffness. A Beam: A long slender structural component generally subjected to transverse loading that produces significant bending effects as opposed to twisting or axial effects. MECH 40: Finite
More informationInvestigation of Utilizing a Secant Stiffness Matrix for 2D Nonlinear Shape Optimization and Sensitivity Analysis
Print ISSN: 2322-2093; Online ISSN: 2423-6691 DOI: 10.7508/ceij.2016.02.011 Technical Note Investigation of Utilizing a Secant Stiffness Matrix for 2D Nonlinear Shape Optimization and Sensitivity Analysis
More informationExample 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.
162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
More informationFig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double cantilever beam (DCB) specimens.
a). Cohesive Failure b). Interfacial Failure c). Oscillatory Failure d). Alternating Failure Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double
More informationDevelopment of X-FEM methodology and study on mixed-mode crack propagation
Acta Mech. Sin. (2011) 27(3):406 415 DOI 10.1007/s10409-011-0436-x RESEARCH PAPER Development of X-FEM methodology and study on mixed-mode crack propagation Zhuo Zhuang Bin-Bin Cheng Received: 2 February
More informationBECAS - an Open-Source Cross Section Analysis Tool
BECAS - an Open-Source Cross Section Analysis Tool José P. Blasques and Robert D. Bitsche Presented at DTU Wind Energy stand at the EWEA 2012 conference, Copenhagen, 16.4.2012 BECAS-DTUWind@dtu.dk Motivation
More informationFEA CODE WITH MATLAB. Finite Element Analysis of an Arch ME 5657 FINITE ELEMENT METHOD. Submitted by: ALPAY BURAK DEMIRYUREK
FEA CODE WITH MATAB Finite Element Analysis of an Arch ME 5657 FINITE EEMENT METHOD Submitted by: APAY BURAK DEMIRYUREK This report summarizes the finite element analysis of an arch-beam with using matlab.
More informationM5 Simple Beam Theory (continued)
M5 Simple Beam Theory (continued) Reading: Crandall, Dahl and Lardner 7.-7.6 In the previous lecture we had reached the point of obtaining 5 equations, 5 unknowns by application of equations of elasticity
More informationWeld Fracture. How Residual Stresses Affect Prediction of Brittle Fracture. Outline. Residual Stress in Thick Welds
How Residual Stresses ffect Prediction of Brittle Fracture Michael R. Hill University of California, Davis Tina L. Panontin NS-mes Research Center Weld Fracture Defects provide location for fracture initiation
More informationFracture mechanics fundamentals. Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design
Fracture mechanics fundamentals Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design Failure modes Failure can occur in a number of modes: - plastic deformation
More informationUsing the finite element method of structural analysis, determine displacements at nodes 1 and 2.
Question 1 A pin-jointed plane frame, shown in Figure Q1, is fixed to rigid supports at nodes and 4 to prevent their nodal displacements. The frame is loaded at nodes 1 and by a horizontal and a vertical
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationStructural Analysis III Compatibility of Displacements & Principle of Superposition
Structural Analysis III Compatibility of Displacements & Principle of Superposition 2007/8 Dr. Colin Caprani, Chartered Engineer 1 1. Introduction 1.1 Background In the case of 2-dimensional structures
More informationFINITE ELEMENT MODELS WITH NODE-DEPENDENT KINEMATICS ADOPTING LEGENDRE POLYNOMIAL EXPANSIONS FOR THE ANALYSIS OF LAMINATED PLATES
FINITE ELEMENT MODELS WITH NODE-DEPENDENT KINEMATICS ADOPTING LEGENDRE POLYNOMIAL EXPANSIONS FOR THE ANALYSIS OF LAMINATED PLATES G. Li, A.G. de Miguel, E. Zappino, A. Pagani and E. Carrera Department
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationINFLUENCE OF THE LOCATION AND CRACK ANGLE ON THE MAGNITUDE OF STRESS INTENSITY FACTORS MODE I AND II UNDER UNIAXIAL TENSION STRESSES
INFLUENCE OF THE LOCATION AND CRACK ANGLE ON THE MAGNITUDE OF STRESS INTENSITY FACTORS MODE I AND II UNDER UNIAXIAL TENSION STRESSES Najah Rustum Mohsin Southern Technical University, Technical Institute-Nasiriya,
More informationComputational Analysis for Composites
Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics
More informationInfluence of impact velocity on transition time for V-notched Charpy specimen*
[ 溶接学会論文集第 35 巻第 2 号 p. 80s-84s (2017)] Influence of impact velocity on transition time for V-notched Charpy specimen* by Yasuhito Takashima** and Fumiyoshi Minami** This study investigated the influence
More informationNonlinear Thermo- Mechanics of Plates and Shallow Shells
Nonlinear Thermo- Mechanics of Plates and Shallow Shells Payam Khazaeinejad 1, Asif S. Usmani 1, Omar Laghrouche 1 IIE, School of Engineering, The University of Edinburgh IIE, School of the Built Environment,
More informationHardened Concrete. Lecture No. 16
Hardened Concrete Lecture No. 16 Fatigue strength of concrete Modulus of elasticity, Creep Shrinkage of concrete Stress-Strain Plot of Concrete At stress below 30% of ultimate strength, the transition
More informationCIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass
CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional
More informationThe Finite Element Method for Mechonics of Solids with ANSYS Applicotions
The Finite Element Method for Mechonics of Solids with ANSYS Applicotions ELLIS H. DILL 0~~F~~~~"P Boca Raton London New Vork CRC Press is an imprint 01 the Taylor & Francis Group, an Informa business
More informationThe minus sign indicates that the centroid is located below point E. We will relocate the axis as shown in Figure (1) and take discard the sign:
AOE 304: Thin Walled Structures Solutions to Consider a cantilever beam as shown in the attached figure. At the tip of the beam, a bending moment M = 1000 N-m is applied at an angle θ with respect to the
More informationThick-Walled Open Section Beam and Finite Nodal-Line Method
Solids and Structures (SAS) Volume 5, 206 doi: 0.4355/sas.206.05.00 www.seipub.org/sas/ hick-walled Open Section Beam and Finite Nodal-Line Method Yaoqing Gong *, Sai ao 2 School of Civil Engineering,
More informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationGeometry-dependent MITC method for a 2-node iso-beam element
Structural Engineering and Mechanics, Vol. 9, No. (8) 3-3 Geometry-dependent MITC method for a -node iso-beam element Phill-Seung Lee Samsung Heavy Industries, Seocho, Seoul 37-857, Korea Hyu-Chun Noh
More informationProbabilistic Fatigue and Damage Tolerance Analysis for General Aviation
Probabilistic Fatigue and Damage Tolerance Analysis for General Aviation Probabilistic fatigue and damage tolerance tool for the Federal Aviation Administration to perform risk analysis Juan D. Ocampo
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationHomework No. 1 MAE/CE 459/559 John A. Gilbert, Ph.D. Fall 2004
Homework No. 1 MAE/CE 459/559 John A. Gilbert, Ph.D. 1. A beam is loaded as shown. The dimensions of the cross section appear in the insert. the figure. Draw a complete free body diagram showing an equivalent
More informationCE 715: Advanced Strength of Materials
CE 715: Advanced Strength of Materials Lecture 1 CE 715 Course Information Instructor: Tasnim Hassan Office: Mann Hall 419 Office Hours: TTh 2:00-4:00 pm Phone: 515-8123 Email: thassan@eos.ncsu.edu 1 Course
More informationStress intensity factors determination for an inclined central crack on a plate subjected to uniform tensile loading using FE analysis
Stress intensity factors determination for an inclined central crack on a plate subjected to uniform tensile loading using FE analysis Paulo C. M. Azevedo IDMEC and Faculdade de Engenharia da Universidade
More informationFracture Mechanics of Composites with Residual Thermal Stresses
J. A. Nairn Material Science & Engineering, University of Utah, Salt Lake City, Utah 84 Fracture Mechanics of Composites with Residual Thermal Stresses The problem of calculating the energy release rate
More informationa x Questions on Classical Solutions 1. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress
Questions on Classical Solutions. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress σ xy = T is applied at infinity. Determine the value of the stress σ θθ on the edge
More informationFinite Element Method
Finite Element Method Finite Element Method (ENGC 6321) Syllabus Objectives Understand the basic theory of the FEM Know the behaviour and usage of each type of elements covered in this course one dimensional
More informationSIZE EFFECTS IN THE COMPRESSIVE CRUSHING OF HONEYCOMBS
43rd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Con 22-25 April 2002, Denver, Colorado SIZE EFFECTS IN THE COMPRESSIVE CRUSHING OF HONEYCOMBS Erik C. Mellquistand Anthony M.
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationMODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
More informationDue Monday, November 16 th, 12:00 midnight
Due Monday, November 16 th, 12:00 midnight This homework is considering the finite element analysis of transient and dynamic FEM analysis. You are asked to include transient and/dynamic effects to MatLab
More informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationLecture M1 Slender (one dimensional) Structures Reading: Crandall, Dahl and Lardner 3.1, 7.2
Lecture M1 Slender (one dimensional) Structures Reading: Crandall, Dahl and Lardner 3.1, 7.2 This semester we are going to utilize the principles we learnt last semester (i.e the 3 great principles and
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationTentamen/Examination TMHL61
Avd Hållfasthetslära, IKP, Linköpings Universitet Tentamen/Examination TMHL61 Tentamen i Skademekanik och livslängdsanalys TMHL61 lördagen den 14/10 2000, kl 8-12 Solid Mechanics, IKP, Linköping University
More informationPh.D. Preliminary Examination Analysis
UNIVERSITY OF CALIFORNIA, BERKELEY Spring Semester 2014 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials Name:......................................... Ph.D.
More informationUniversity of Sheffield The development of finite elements for 3D structural analysis in fire
The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures
More informationModeling Fracture and Failure with Abaqus
Modeling Fracture and Failure with Abaqus Day 1 Lecture 1 Lecture 2 Lecture 3 Workshop 1 Workshop 2 Basic Concepts of Fracture Mechanics Modeling Cracks Fracture Analysis Crack in a Three-point Bend Specimen
More informationRecommended teaching texts Theory:
Recommended teaching texts Theory: M. Sochor: Strength of materials II Czech Technical University Prague, 2006. In Czech: Ondráček, Vrbka, Janíček, Burša: Mechanika těles, Pružnost a pevnost II. CERM,
More informationELASTICITY AND FRACTURE MECHANICS. Vijay G. Ukadgaonker
THEORY OF ELASTICITY AND FRACTURE MECHANICS y x Vijay G. Ukadgaonker Theory of Elasticity and Fracture Mechanics VIJAY G. UKADGAONKER Former Professor Indian Institute of Technology Bombay Delhi-110092
More informationEshan V. Dave, Secretary of M&FGM2006 (Hawaii) Research Assistant and Ph.D. Candidate. Glaucio H. Paulino, Chairman of M&FGM2006 (Hawaii)
Asphalt Pavement Aging and Temperature Dependent Properties through a Functionally Graded Viscoelastic Model I: Development, Implementation and Verification Eshan V. Dave, Secretary of M&FGM2006 (Hawaii)
More informationComposite Structural Mechanics using MATLAB
Session 2520 Composite Structural Mechanics using MATLAB Oscar Barton, Jr., Jacob B. Wallace United States Naval Academy Annapolis, Md 21402 Abstract In this paper MATLAB is adopted as the programming
More informationDue Tuesday, September 21 st, 12:00 midnight
Due Tuesday, September 21 st, 12:00 midnight The first problem discusses a plane truss with inclined supports. You will need to modify the MatLab software from homework 1. The next 4 problems consider
More informationSoftware Verification
EXAMPLE 1-026 FRAME MOMENT AND SHEAR HINGES EXAMPLE DESCRIPTION This example uses a horizontal cantilever beam to test the moment and shear hinges in a static nonlinear analysis. The cantilever beam has
More informationDevelopment of discontinuous Galerkin method for linear strain gradient elasticity
Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace
More informationDISPENSA FEM in MSC. Nastran
DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation
More informationCracks Jacques Besson
Jacques Besson Centre des Matériaux UMR 7633 Mines ParisTech PSL Research University Institut Mines Télécom Aγνωστ oς Θεoς Outline 1 Some definitions 2 in a linear elastic material 3 in a plastic material
More informationMulti-scale digital image correlation of strain localization
Multi-scale digital image correlation of strain localization J. Marty a, J. Réthoré a, A. Combescure a a. Laboratoire de Mécanique des Contacts et des Strcutures, INSA Lyon / UMR CNRS 5259 2 Avenue des
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationFRACTURE MECHANICS OF COMPOSITES WITH RESIDUAL STRESSES, TRACTION-LOADED CRACKS, AND IMPERFECT INTERFACES
Proc. 2 nd ESIS TC4 Conference on Polymers and Composites, in press, 1999 Author prepared reprint FRACTURE MECHANICS OF COMPOSITES WITH RESIDUAL STRESSES, TRACTION-LOADED CRACKS, AND IMPERFECT INTERFACES
More informationExample 37 - Analytical Beam
Example 37 - Analytical Beam Summary This example deals with the use of RADIOSS linear and nonlinear solvers. A beam submitted to a concentrated load on one extremity and fixed on the other hand is studied.
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationM.S Comprehensive Examination Analysis
UNIVERSITY OF CALIFORNIA, BERKELEY Spring Semester 2014 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials Name:......................................... M.S Comprehensive
More informationThe Ultimate Load-Carrying Capacity of a Thin-Walled Shuttle Cylinder Structure with Cracks under Eccentric Compressive Force
The Ultimate Load-Carrying Capacity of a Thin-Walled Shuttle Cylinder Structure with Cracks under Eccentric Compressive Force Cai-qin Cao *, Kan Liu, Jun-zhe Dong School of Science, Xi an University of
More informationMicro-meso draping modelling of non-crimp fabrics
Micro-meso draping modelling of non-crimp fabrics Oleksandr Vorobiov 1, Dr. Th. Bischoff 1, Dr. A. Tulke 1 1 FTA Forschungsgesellschaft für Textiltechnik mbh 1 Introduction Non-crimp fabrics (NCFs) are
More informationTRANSVERSE STRESSES IN SHEAR LAG OF BOX-GIRDER BRIDGES. Wang Yuan
TRANSVERSE STRESSES IN SHEAR LAG OF BOX-GIRDER BRIDGES Wang Yuan Bacheloreindwerk Delft University of Technology Faculty of Civil Engineering and Geosciences October 2011 TABLE OF CONTENTS 1 INTRODUCTION...
More informationA 3D Discrete Damage Modeling Methodology for Abaqus for Fatigue Damage Evaluation in Bolted Composite Joints
A 3D Discrete Damage Modeling Methodology for Abaqus for Fatigue Damage Evaluation in Bolted Composite Joints Eugene Fang 1, Ling Liu 1, Michael Stuebner 1, Jim Lua 1 1 Global Engineering and Materials,
More informationStructural Displacements
Structural Displacements Beam Displacement 1 Structural Displacements P Truss Displacements 2 The deflections of civil engineer- ing structures under the action of usual design loads are known to be small
More informationMeshless equilibrium on line method (MELM) for linear elasticity
Structural Engineering and Mechanics, Vol. 35, No. 4 (2010) 000-000 1 Meshless equilibrium on line method (MELM) for linear elasticity A. Sadeghirad*, S. Mohammadi and I. Mahmoudzadeh Kani School of Civil
More informationLecture 8: Assembly of beam elements.
ecture 8: Assembly of beam elements. 4. Example of Assemblage of Beam Stiffness Matrices. Place nodes at the load application points. Assembling the two sets of element equations (note the common elemental
More informationFinite Element Analysis Lecture 1. Dr./ Ahmed Nagib
Finite Element Analysis Lecture 1 Dr./ Ahmed Nagib April 30, 2016 Research and Development Mathematical Model Mathematical Model Mathematical Model Finite Element Analysis The linear equation of motion
More information2 marks Questions and Answers
1. Define the term strain energy. A: Strain Energy of the elastic body is defined as the internal work done by the external load in deforming or straining the body. 2. Define the terms: Resilience and
More informationFLOATING NODE METHOD AND VIRTUAL CRACK CLOSURE TECHNIQUE FOR MODELING MATRIX CRACKING- DELAMINATION MIGRATION
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS FLOATING NODE METHOD AND VIRTUAL CRACK CLOSURE TECHNIQUE FOR MODELING MATRIX CRACKING- DELAMINATION MIGRATION N. V. De Carvalho 1*, B. Y. Chen
More informationSome improvements of Xfem for cracked domains
Some improvements of Xfem for cracked domains E. Chahine 1, P. Laborde 2, J. Pommier 1, Y. Renard 3 and M. Salaün 4 (1) INSA Toulouse, laboratoire MIP, CNRS UMR 5640, Complexe scientifique de Rangueil,
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationPILE SOIL INTERACTION MOMENT AREA METHOD
Pile IGC Soil 2009, Interaction Moment Guntur, INDIA Area Method PILE SOIL INTERACTION MOMENT AREA METHOD D.M. Dewaikar Professor, Department of Civil Engineering, IIT Bombay, Mumbai 400 076, India. E-mail:
More informationWorkshop 8. Lateral Buckling
Workshop 8 Lateral Buckling cross section A transversely loaded member that is bent about its major axis may buckle sideways if its compression flange is not laterally supported. The reason buckling occurs
More informationFracture Mechanics, Damage and Fatigue Linear Elastic Fracture Mechanics - Energetic Approach
University of Liège Aerospace & Mechanical Engineering Fracture Mechanics, Damage and Fatigue Linear Elastic Fracture Mechanics - Energetic Approach Ludovic Noels Computational & Multiscale Mechanics of
More informationHPLP300 - Plate with Young modulus function of the temperature
Titre : HPLP300 - Plaque aec module d Young fonction de l[...] Date : 04/08/2011 Page : 1/6 HPLP300 - Plate with Young modulus function of the temperature Summarized: This thermoelastic test makes it possible
More informationStress Functions. First Semester, Academic Year 2012 Department of Mechanical Engineering Chulalongkorn University
Lecture Note 3 Stress Functions First Semester, Academic Year 01 Department of Mechanical Engineering Chulalongkorn Universit 1 Objectives #1 Describe the equation of equilibrium using D/3D stress elements
More informationCritical applied stresses for a crack initiation from a sharp V-notch
Focussed on: Fracture and Structural Integrity related Issues Critical applied stresses for a crack initiation from a sharp V-notch L. Náhlík, P. Hutař Institute of Physics of Materials, Academy of Sciences
More informationNON-LINEAR FRACTURE BEHAVIOR OF DOUBLE CANTILEVER BEAM
Engineering MECHANICS, Vol., 015, No., p. 95 10 95 NON-LINEAR FRACTURE BEHAVIOR OF DOUBLE CANTILEVER BEAM Viktor Rizov* This article describes a theoretical study of non-linear fracture behavior of the
More informationFEA A Guide to Good Practice. What to expect when you re expecting FEA A guide to good practice
FEA A Guide to Good Practice What to expect when you re expecting FEA A guide to good practice 1. Background Finite Element Analysis (FEA) has transformed design procedures for engineers. Allowing more
More informationACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN XFEM APPROACH
VI International Conference on Adaptive Modeling and Simulation ADMOS 213 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationMaterials and Structures
Journal of Mechanics of Materials and Structures BRITTLE FRACTURE BEYOND THE STRESS INTENSITY FACTOR C. T. Sun and Haiyang Qian Volume 4, Nº 4 April 2009 mathematical sciences publishers JOURNAL OF MECHANICS
More information