9/ 29/ 12

St r ess anal si of a pl t e w h a hol y s a t i e

The OpenFOAM Foundation

Home
Le gal

Download

Documentation

Licensing

Code Development

OpenFOAM.com

User Guide
Contents 1 Introduction 2 Tutorials
2.1 Lid-driven cavity flow 2.2 Stress analysis of a plate w ith a hole 2.3 Breaking of a dam

2.2 Stress analysis of a plate with a hole

This tutorial describes how to pre-process, run and post-process a case involving linear-elastic, steady-state stress analysis on a square plate with a circular hole at its centre. The plate dimensions are: side length 4 and radius 0.5 . It is loaded with a uniform traction of 10 over its left and right faces as shown in Figure 2.15. Two symmetry planes can be identified for this geometry and therefore the solution domain need only cover a quarter of the geometry, shown by the shaded area in Figure 2.15.

3 Applications and libraries 4 OpenFOAM cases 5 Mesh generation and conversion 6 Post-processing 7 Models and physical properties Index

Figure 2.15: Geometry of the plate with a hole. The problem can be approximated as 2-dimensional since the load is applied in the plane of the plate. In a Cartesian coordinate system there are two possible assumptions to take in regard to the behaviour of the structure in the third dimension: (1) the plane stress condition, in which the stress components acting out of the 2D plane are assumed to be negligible; (2) the plane strain condition, in which the strain components out of the 2D plane are assumed negligible. The plane stress condition is appropriate for solids whose third dimension is thin as in this case; the plane strain condition is applicable for solids where the third dimension is thick. An analytical solution exists for loading of an infinitely large, thin plate with a circular hole. The solution for the stress normal to the vertical plane of symmetry is

(2.14)

Results from the simulation will be compared with this solution. At the end of the tutorial, the user can: investigate the sensitivity of the solution to mesh resolution and mesh grading; and, increase the size of the plate in comparison to the hole to try to estimate the error in comparing the analytical solution for an infinite plate to the solution of this problem of a finite plate.

2.2.1 Mesh generation

The domain consists of four blocks, some of which have arc-shaped edges. The block structure for the part of the mesh in the plane is shown in Figure 2.16. As already mentioned in section 2.1.1.1, all geometries are generated in 3 dimensions in OpenFOAM even if the case is to be as a 2 dimensional problem. Therefore a dimension of the block in the direction has to be chosen; here, 0.5 is selected. It does not affect the solution since the traction boundary condition is specified as a stress rather than a force, thereby making the solution independent of the cross-sectional area.

w w openf oam . or g/ docs/ user / pl t eH e. php w. a ol

1/ 10

9/ 29/ 12

St r ess anal si of a pl t e w h a hol y s a t i e

Figure 2.16: Block structure of the mesh for the plate with a hole. The user should change into the plateHole case in the $FOAM_RUN/tutorials/stressAnalysis/solidDisplacementFoam directory and open the constant/polyMesh/b lockMeshDict file in an editor, as listed below
17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50

cnetoees1 ovrTMtr ;

vrie etcs ( (. 00 05 ) ( 00 1 ) ( 00 2 ) ( 0770 0 2 .017 ) (.0170770 0 0770 .017 ) (.5530335 0 0335 .553 ) ( 20 2 ) (.01720 0770 ) ( 20 0 ) ( 10 0 ) ( 050 0 . ) (. 005 05 .) ( 005 1 .) ( 005 2 .) ( 0770 05 2 .017 .) (.0170770 05 0770 .017 .) (.5530335 05 0335 .553 .) ( 205 2 .) (.017205 0770 .) ( 205 0 .) ( 105 0 .) ( 0505 0 . .) ) ;

bok lcs ( hx( 491 1 1 2 2)(01 1 smlGaig( 11 e 5 0 6 5 0 1 1 0 ) iperdn 1 ) hx( 1451 1 1 1)(01 1 smlGaig( 11 e 0 1 2 5 6 1 0 ) iperdn 1 ) hx( 2341 1 1 1)(01 1 smlGaig( 11 e 1 2 3 4 5 2 0 ) iperdn 1 ) hx( 3671 1 1 1)(02 1 smlGaig( 11 e 4 5 4 7 8 2 0 ) iperdn 1 )

w w openf oam . or g/ docs/ user / pl t eH e. php w. a ol

2/ 10

9/ 29/ 12

St r ess anal si of a pl t e w h a hol y s a t i e

51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113

hx( 4782 1 1 1)(02 1 smlGaig( 11 e 9 0 5 8 9 1 0 ) iperdn 1 ) ) ;

egs de ( ac05(.68601110 r 0494 .70 ) ac51 (.70 0494 0 r 0 0111 .686 ) ac14(.36303220 r 0999 .40 ) ac49(.40 0999 0 r 0322 .363 ) ac1 1 (.686011105 r 1 6 0494 .70 .) ac1 2 (.70 0494 05 r 6 1 0111 .686 .) ac1 1 (.363032205 r 2 5 0999 .40 .) ac1 2 (.40 0999 05 r 5 0 0322 .363 .) ) ;

budr onay ( lf et { tp smerPae ye ymtyln; fcs ae ( ( 92 1) 8 0 9 ( 1 2 2) 9 0 1 0 ) ; } rgt ih { tp pth ye ac; fcs ae ( ( 31 1) 2 4 3 ( 61 1) 3 7 4 ) ; } dw on { tp smerPae ye ymtyln; fcs ae ( ( 11 1) 0 2 1 ( 21 1) 1 3 2 ) ; } u p { tp pth ye ac; fcs ae ( ( 81 1) 7 9 8 ( 71 1) 6 8 7 ) ; } hl oe { tp pth ye ac; fcs ae ( (051 2) 1 6 1 ( 01 1) 5 1 6 ) ; } fotnBc rnAdak

w w openf oam . or g/ docs/ user / pl t eH e. php w. a ol

3/ 10

9/ 29/ 12

St r ess anal si of a pl t e w h a hol y s a t i e

114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136

{ tp epy ye mt; fcs ae ( (0945 1 ) ( 410 5 ) ( 432 1 ) ( 763 4 ) ( 987 4 ) (11 1 2) 2 6 5 0 (61 1 1) 1 1 2 5 (21 1 1) 1 3 4 5 (51 1 1) 1 4 7 8 (51 1 2) 1 8 9 0 ) ; } ) ;

mreacPis egPthar ( ) ;

/ *************************************/ / ************************************ /

Until now, we have only specified straight edges in the geometries of previous tutorials but here we need to specify curved edges. These are specified under the e g s keyword entry which is a list of non-straight edges. The syntax of each list entry de begins with the type of curve, including a c s m l S l n , p l L n etc., described further in section 5.3.1. In this example, r, ipepie oyie all the edges are circular and so can be specified by the a c keyword entry. The following entries are the labels of the start r and end vertices of the arc and a point vector through which the circular arc passes. The blocks in this b lockMeshDict do not all have the same orientation. As can be seen in Figure 2.16 the direction of block 0 is equivalent to the direction for block 4. This means care must be taken when defining the number and distribution of cells in each block so that the cells match up at the block faces. 6 patches are defined: one for each side of the plate, one for the hole and one for the front and back planes. The l f and et d w patches are both a symmetry plane. Since this is a geometric constraint, it is included in the definition of the mesh, on rather than being purely a specification on the boundary condition of the fields. Therefore they are defined as such using a special symmetryPlane type as shown in the b lockMeshDict. The f o t n B c patch represents the plane which is ignored in a 2D case. Again this is a geometric constraint so is rnAdak defined within the mesh, using the empty type as shown in the b lockMeshDict. For further details of boundary types and geometric constraints, the user should refer to section 5.2.1. The remaining patches are of the regular patch type. The mesh should be generated using b lockMesh and can be viewed in paraFoam as described in section 2.1.2. It should appear as in Figure 2.17.

Figure 2.17: Mesh of the hole in a plate problem.

2.2.1.1 Boundary and initial conditions

Once the mesh generation is complete, the initial field with boundary conditions must be set. For a stress analysis case without thermal stresses, only displacement D needs to be set. The 0/D is as follows:
17

dmnin iesos

[ 100000; 0 ]

w w openf oam . or g/ docs/ user / pl t eH e. php w. a ol

4/ 10

9/ 29/ 12

St r ess anal si of a pl t e w h a hol y s a t i e

18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58

itraFed nenlil

uiom( 00; nfr 0 )

budrFed onayil { lf et { tp ye } rgt ih { tp ye tato rcin pesr rsue vle au } dw on { tp ye } u p { tp ye tato rcin pesr rsue vle au } hl oe { tp ye tato rcin pesr rsue vle au } fotnBc rnAdak { tp ye } } epy mt; tatoDslcmn; rcinipaeet uiom(000) nfr ; uiom0 nfr ; uiom( 00; nfr 0 ) tatoDslcmn; rcinipaeet uiom(000) nfr ; uiom0 nfr ; uiom( 00; nfr 0 ) smerPae ymtyln; tatoDslcmn; rcinipaeet uiom(10000) nfr 00 ; uiom0 nfr ; uiom( 00; nfr 0 ) smerPae ymtyln;

/ *************************************/ / ************************************ /

Firstly, it can be seen that the displacement initial conditions are set to

. The l f and d w patches must be et on

both of symmetryPlane type since they are specified as such in the mesh description in the constant/polyMesh/b oundary file. Similarly the f o t n B c patch is declared empty. rnAdak The other patches are traction boundary conditions, set by a specialist traction boundary type. The traction boundary conditions are specified by a linear combination of: (1) a boundary traction vector under keyword t a t o ; (2) a pressure that rcin produces a traction normal to the boundary surface that is defined as negative when pointing out of the surface, under keyword p e s r . The u and h l patches are zero traction so the boundary traction and pressure are set to zero. For the r g t rsue p oe ih patch the traction should be and the pressure should be 0 .

2.2.1.2 Mechanical properties

The physical properties for the case are set in the mechanicalProperties dictionary in the constant directory. For this problem, we need to specify the mechanical properties of steel given in Table 2.1. In the mechanical properties dictionary, the user must also set p a e t e s to y s lnSrs e. Property Density Youngs modulus Poissons ratio Units Keyword ro h E n u Value 7854 0.3 Table 2.1: Mechanical properties for steel

w w openf oam . or g/ docs/ user / pl t eH e. php w. a ol

5/ 10

9/ 29/ 12

St r ess anal si of a pl t e w h a hol y s a t i e

2.2.1.3 Thermal properties

The temperature field variable T is present in the solidDisplacementFoam solver since the user may opt to solve a thermal equation that is coupled with the momentum equation through the thermal stresses that are generated. The user specifies at run time whether OpenFOAM should solve the thermal equation by the t e m l t e s switch in the thermalProperties hraSrs dictionary. This dictionary also sets the thermal properties for the case, e.g. for steel as listed in Table 2.2. Property Specific heat capacity Thermal conductivity Thermal expansion coeff. Units Keyword C k apa lh Value 434 60.5

Table 2.2: Thermal properties for steel In this case we do not want to solve for the thermal equation. Therefore we must set the t e m l t e s keyword entry to hraSrs n in the thermalProperties dictionary. o

2.2.1.4 Control
As before, the information relating to the control of the solution procedure are read in from the controlDict dictionary. For this case, the s a t i e is 0 . The time step is not important since this is a steady state case; in this situation it is best to set trTm the time step d l a to 1 so it simply acts as an iteration counter for the steady-state case. The e d i e set to 100, then etT nTm, acts as a limit on the number of iterations. The w i e n e v l can be set to rtItra . The controlDict entries are as follows:
17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51

apiain plcto

sldipaeetom oiDslcmnFa;

satrm trFo

satie trTm;

satie trTm

0 ;

soA tpt

edie nTm;

edie nTm

10 0;

dla etT

1 ;

wieoto rtCnrl

tmSe; ietp

wienevl rtItra

2; 0

prert ugWie

0 ;

wieomt rtFra

aci si;

wiercso 6 rtPeiin ;

wieopeso of rtCmrsin f;

tmFra ieomt

gnrl eea;

tmPeiin iercso

6 ;

gahomt rpFra

rw a;

rnieoiibetu; uTmMdfal re

/ *************************************/ / ************************************ /

2.2.1.5 Discretisation schemes and linear-solver control

Let us turn our attention to the fvSchemes dictionary. Firstly, the problem we are analysing is steady-state so the user should

w w openf oam . or g/ docs/ user / pl t eH e. php w. a ol

6/ 10

9/ 29/ 12

St r ess anal si of a pl t e w h a hol y s a t i e

select S e d S a e for the time derivatives in t m S h m . This essentially switches off the time derivative terms. Not all taytt iecee solvers, especially in fluid dynamics, work for both steady-state and transient problems but solidDisplacementFoam does work, since the base algorithm is the same for both types of simulation. The momentum equation in linear-elastic stress analysis includes several explicit terms containing the gradient of displacement. The calculations benefit from accurate and smooth evaluation of the gradient. Normally, in the finite volume method the discretisation is based on Gausss theorem The Gauss method is sufficiently accurate for most purposes but, in this case, the least squares method will be used. The user should therefore open the f S h m sdictionary in the system vcee directory and ensure the l a t q a e method is selected for the g a ( ) gradient discretisation scheme in the esSurs rdU g a S h m ssub-dictionary: rdcee
17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66

dd2cee 2tShms { dfut eal } sedSae taytt;

dtcee dShms { dfut eal } Elr ue;

gaShms rdcee { dfut eal ga() rdD ga() rdT } latqae; esSurs latqae; esSurs latqae; esSurs

dvcee iShms { dfut eal dvsga) i(imD } nn; oe Guslna; as ier

lpaincee alcaShms { dfut eal nn; oe

lpainD,)Guslna cretd alca(DD as ier orce; lpainD,)Guslna cretd alca(TT as ier orce; }

itroaincee nepltoShms { dfut eal } lna; ier

sGaShms nrdcee { dfut eal } nn; oe

fuRqie lxeurd { dfut eal D T } n; o ys e; n; o

/ *************************************/ / ************************************ /

The fvSolution dictionary in the system directory controls the linear equation solvers and algorithms used in the solution. The user should first look at the solvers sub-dictionary and notice that the choice of s l e for D is G M . The solver t l r n e ovr AG oeac should be set to for this problem. The solver relative tolerance, denoted by r l o , sets the required reduction in the eTl

w w openf oam . or g/ docs/ user / pl t eH e. php w. a ol

7/ 10

9/ 29/ 12

St r ess anal si of a pl t e w h a hol y s a t i e

residuals within each iteration. It is uneconomical to set a tight (low) relative tolerance within each iteration since a lot of terms in each equation are explicit and are updated as part of the segregated iterative procedure. Therefore a reasonable value for the relative tolerance is , or possibly even higher, say , or in some cases even (as in this case).
17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41

sles ovr { "DT" (|) { sle ovr tlrne oeac rlo eTl sote mohr GM; AG 1-6 e0; 09 .; Gusedl asSie;

cceglmrto tu; ahAgoeain re nelICassLvl2; Clsnoretee 0 agoeao glmrtr mreees egLvl } } fcAePi; aeraar 1 ;

srsAayi tesnlss { cmatomltesys opcNraSrs e; norcos Cretr D } 1 ; 1-6 e0;

/ *************************************/ / ************************************ /

The fvSolution dictionary contains a sub-dictionary, stressAnalysis that contains some control parameters specific to the application solver. Firstly there is n o r c o swhich specifies the number of outer loops around the complete system of Cretr equations, including traction boundary conditions within each time step. Since this problem is steady-state, we are performing a set of iterations towards a converged solution with the time step acting as an iteration counter. We can therefore set n o r c o s to 1. Cretr The D keyword specifies a convergence tolerance for the outer iteration loop, i.e. sets a level of initial residual below which solving will cease. It should be set to the desired solver tolerance specified earlier, for this problem.

2.2.2 Running the code

The user should run the code here in the background from the command line as specified below, so he/she can look at convergence information in the log file afterwards. c $F A _ U / u o i l / t e s n l s s s l d i p a e e t o m p a e o e d OMRNttrassrsAayi/oiDslcmnFa/ltHl sldipaeetom>lg& oiDslcmnFa o The user should check the convergence information by viewing the generated log file which shows the number of iterations and the initial and final residuals of the displacement in each direction being solved. The final residual should always be less than 0.9 times the initial residual as this iteration tolerance set. Once both initial residuals have dropped below the convergence tolerance of the run has converged and can be stopped by killing the batch job.

2.2.3 Post-processing
Post processing can be performed as in section 2.1.4. The solidDisplacementFoam solver outputs the stress field as a symmetric tensor field s g a This is consistent with the way variables are usually represented in OpenFOAM solvers by the im. mathematical symbol by which they are represented; in the case of Greek symbols, the variable is named phonetically. For post-processing individual scalar field components, , etc., can be generated by running the foamCalc utility as before in section 2.1.5.7, this time on s g a im: faCl cmoet sga omac opnns im Components named s g a x s g a y etc. are written to time directories of the case. The imx, imx paraFoam as shown in Figure 2.18. stresses can be viewed in

w w openf oam . or g/ docs/ user / pl t eH e. php w. a ol

8/ 10

9/ 29/ 12

St r ess anal si of a pl t e w h a hol y s a t i e

Figure 2.18: stress field in the plate with hole. We would like to compare the analytical solution of Equation 2.14 to our solution. We therefore must output a set of data of along the left edge symmetry plane of our domain. The user may generate the required graph data using the sample utility. The utility uses a sampleDict dictionary located in the system directory, whose entries are summarised in Table 6.3. The sample line specified in s t is set between es and , and the fields are specified in the f e d list: ils
17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37

itroainceeclPit nepltoShm elon;

stomt eFra

rw a;

st es ( lfPth etac { tp ye ai xs sat tr ed n uiom nfr; y ; (00502 ) . .5 ; (0202 ) .5 ;

nons10 Pit 0; } ) ;

fed ils

(sgax) imx ;

/ *************************************/ / ************************************ /

The user should execute sample as normal. The w i e o m tis r w 2 column format. The data is written into files within rtFra a time subdirectories of a sets directory, e.g. the data at s is found within the file sets/100/leftPatch_sigmaxx.xy. In an application such as GnuPlot, one could type the following at the command prompt would be sufficient to plot both the numerical data and analytical solution: p o [ . : ] [ : 's t / 0 / e t a c _ i m x . y l t 0 5 2 0 ] e s 1 0 l f P t h s g a x x ', 14(+015(*2)(.97/x*)) e*1(.2/x*)+0035(*4) An example plot is shown in Figure 2.19.

w w openf oam . or g/ docs/ user / pl t eH e. php w. a ol

9/ 10

9/ 29/ 12

St r ess anal si of a pl t e w h a hol y s a t i e

Figure 2.19: Normal stress along the vertical symmetry

2.2.4 Exercises
The user may wish to experiment with solidDisplacementFoam by trying the following exercises:

2.2.4.1 Increasing mesh resolution

Increase the mesh resolution in each of the and directions. Use mapFields to map the final coarse mesh results from section 2.2.3 to the initial conditions for the fine mesh.

2.2.4.2 Introducing mesh grading

Grade the mesh so that the cells near the hole are finer than those away from the hole. Design the mesh so that the ratio of sizes between adjacent cells is no more than 1.1 and so that the ratio of cell sizes between blocks is similar to the ratios within blocks. Mesh grading is described in section 2.1.6. Again use mapFields to map the final coarse mesh results from section 2.2.3 to the initial conditions for the graded mesh. Compare the results with those from the analytical solution and previous calculations. Can this solution be improved upon using the same number of cells with a different solution?

2.2.4.3 Changing the plate size

The analytical solution is for an infinitely large plate with a finite sized hole in it. Therefore this solution is not completely accurate for a finite sized plate. To estimate the error, increase the plate size while maintaining the hole size at the same value.

This User Guide is licensed under a Creative Commons Attribution-NonCommercial-NoDerivs 3.0 Unported License .

Copyright 2011-2012 OpenFOAM Foundation | OPENFOAM and OpenCFD are registered tradem arks of OpenCFD Ltd. Design by 1234.info | Modified by OpenFOAM Foundation | Content generated using tex4ht

w w openf oam . or g/ docs/ user / pl t eH e. php w. a ol

10/ 10