Meshfree Inelastic Frame Analysis

Transcription

1 Theory & Results Louie L. Yaw, Sashi Kunnath and N. Sukumar University of California, Davis Department of Civil and Environmental Engineering Minisymposium 47 Recent Advances in Modeling of Engineering Materials/Systems USNCCM9 San Francisco July 24, 2007

2 Acknowledgements Funding provided by Walla Walla College N. Sukumar acknowledges research support of NSF (Grant CMMI ) Helpful discussions with Dr. Michael Puso, Lawrence Livermore National Laboratory Helpful discussions with Professor Boris Jeremic, UC Davis

3 Outline 1 Motivation

4 Motivation Goal to advance collapse simulation technology Current FE technology unsatisfactory for large deformations at collapse limit states To explore feasibility of meshfree approach

5 Outline Motivation 1 Motivation

6 MLS Shape Functions Derivation Start with displacement approximation u h (x) = n φ a (x)d a φ T d a=1 Shape function φ a is of the form (Belytschko et al 1996) φ a (x) = P T (x a )α(x)w(x a ), where P(x) = {1 x y} T is a linear basis in two dimensions, α(x) is a vector of unknowns to be determined and w(x) 0 is a weighting function

7 MLS Shape Functions Derivation The φ s must satisfy reproducing conditions P(x) = n P(x a )φ a (x) a=1 Substitution of φ a into P(x) and solving for α yields α(x) = A 1 (x)p(x) Finally, substituting α into φ a gives φ a (x) = P T (x a )A 1 (x)p(x)w(x a ), where A = n a=1 P(x a)p T (x a )w(x a )

8 Enforcing boundary conditions MLS (meshfree) shape functions do not have the Kronecker-delta property Hence MLS shape functions are blended with quadrilateral FE shape functions at essential B.C. s Then essential boundary conditions are enforced on the finite element nodes in the standard way The blending technique proposed by Huerta and Fernández-Méndez (2004) is adopted

9 1 Weigth function for node 5 1 1D MLS shape functions 0.8 w Shape Function Node 5 Weight function Shape functions Radius of support ρ = (a) (b) 1.2 1D Blending example 1 Shape Function MLS Shape Functions Blended FE Shape Region Functions x

10 Nodal Integration Smoothed strain tensor for node a (per Chen et al (2001)) ε ij (x a ) = 1 (u i,j + u j,i ) dv = 1 (u i n j + u j n i ) ds 2A a V a 2A a S a Strain-displacement relation ε(x a ) = 6 B b (x a )d b Bd b=1

11 Nodal Integration Strain-displacement definitions ε = [ε 11 ε 22 2ε 12 ] T and d a = [d a1 d a2 ] T b b1 (x a ) 0 B b (x a ) = 0 b b2 (x a ) b b2 (x a ) b b1 (x a ) 3 x 2 2 V a a = 1 4 n x 1 S a 6 5 b bi (x a ) = 1 A a S a φ b (x)n i (x) ds Nodally integrated stiffness matrix K bc = n B T b (x a)cb c (x a )A a t a=1

12 Outline Motivation 1 Motivation

13 Stabilization of Nodal Integration Nodal integration w/o stabilization leads to a. hourglass modes b. spurious low energy modes c. and locking Following Puso and Solberg (2006) stabilization is provided to the stiffness matrix as follows: K s = (1 α s )K MLS + α s K FE, where K s is the stabilized matrix and α s = 0.05 is called the stabilization factor

14 Nonlinear analysis Loads applied incrementally At global level a Newton-Raphson scheme is used to iterate the linearized system of equations until equilibrium is achieved K t(ν) (ν) n+1 d n = f ext n+1 f int(ν) n+1 At the constitutive level for J2 plasticity a radial return scheme is used (Simo and Hughes (1998))

15 Meshfree analysis of wide-flange steel sections Based on the meshfree nodal discretization of a beam a Voronoi diagram is generated. A thickness is specified for each Voronoi cell. To get wide-flange behavior a web thickness and a flange thickness is specified.

16 Outline Motivation 1 Motivation

17 Results for normalized tip displacement and maximum bending stress, where δ theor. = in and σ theor. = 25.0 ksi. Grid δ/δ theor. (in) σ xx /σ theor. (ksi) P

18 Outline Motivation 1 Motivation

19 Elasto-plastic cantilever I-beam J2 plasticity with linear hardening σ P (kips) σ m σ y E ε y H ε m ε P 2 Meshfree Analytical Tip Displacement δ (in)

20 Frame corner connection: Load deflection response P 11 ft 3 in Load P (kips) Deflection Gage 11 ft 3 in P Experimental (Beedle 1964) Meshfree - Exponential Hardening Meshfree - Linear Hardening Displacement δ (in)

21 Frame corner connection: Displacement and stress

22 Inelastic frame analysis ft P δ P P (kips) 6 8 ft 4 16 ft 2 Experimental (Baker 1952) MLS Shape Functions Maxent Shape Functions δ (in)

23 Demonstrated feasiblity of wide-flange beam analysis under plane stress A coupled FE and meshfree method shows promise for inelastic frame analysis Demonstrated success of Maxent shape functions Further research is ongoing to extend this for large deformations to enable collapse simulations.