You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 113 Next »

Error formatting macro: include: java.lang.IllegalArgumentException: No link could be created for 'SIMULATION:Tensile Bar - Panel'.

Pre-Analysis and Start-Up

UNDER CONSTRUCTION

In the Pre-Analysis step, we'll review the following:

  • Mathematical Model: We'll look at the governing equations + boundary conditions and the assumptions contained within the mathematical model.
  • Hand-calculations of expected results: We'll use an analytical solution of the mathematical model to predict the expected stress field from ANSYS. We'll pay close attention to additional assumptions that have to be made in order to obtain an analytical solution.

Mathematical Model

We'll first list the assumptions in the mathematical model. Then, we'll review the governing equations and boundary conditions that form the mathematical model. Note that this type of a mathematical model where you have a set of differential equations together with a set of additional restraints at the boundaries is called a Boundary Value Problem (BVP). A lot of practical problems that are solved using ANSYS and other FEA software are BVP's. You should have encountered simple BVP's in your math courses, problems of the kind that involve solving a differential equation with a set of boundary conditions (I was never good at these math problems and it showed in my math grades to the displeasure of my parents .... fortunately that is now a distant memory!). You can think of the BVP considered in this tutorial as a souped-up version of simpler BVP's you have encountered in math courses (and either liked or hated!).

Assumptions

We'll assume that:

  1. Plane stress conditions apply since the bar is thin, thus we don't expect significant variation of stresses in the z direction:
    Unknown macro: {latex}

    [
    \sigma_

    Unknown macro: {zz}

    = \tau_

    Unknown macro: {xz}

    = \tau_

    Unknown macro: {yz}

    = 0
    ]

  2. Gravity effects can be neglected i.e. no body forces.
    Unknown macro: {latex}

    [
    F_x = F_y =0
    ]

Governing Equations

Since we are assuming plane stress conditions, we can use the 2D version of the equilibrium equations. When the deformed structure reaches equilibrium, the 2D stress components should satisfy the 2D equilibrium equations with zero body forces:

Unknown macro: {latex}

\begin

Unknown macro: {eqnarray}
Unknown macro: {partial sigma_x over partial x}

+ {\partial \tau_

Unknown macro: {yx}

\over \partial y} = 0 \nonumber
{\partial \tau_

Unknown macro: {xy}

\over \partial x} +

Unknown macro: {partial sigma_y over partial y}

= 0 \nonumber
\end

Boundary Conditions

We solve these equations in a rectangular domain and impose the appropriate boundary conditions. At every point on the boundary, either the displacement or the traction must be prescribed.

((Figure))

The bottom and top edges are free. Since traction on free edges are zero, we get

Unknown macro: {latex}

[
\sigma_y = \tau_

Unknown macro: {xy}

= 0 \:\: at \: y = 0 \: and \: y = w
]

The left end is fixed. So both components of displacement are zero at this end:

Unknown macro: {latex}

[
u = v = 0 \:\: at \: x = 0
]

At the right end, the traction is specified at the mid-point where the point load is applied. The applied traction at all other points on the right boundary is zero. For brevity, we won't write out the corresponding equations at the right boundary. We'll simplify this boundary condition in our hand calculations below (to make the problem tractable) but the ANSYS solution provided uses the full set of boundary conditions. Another complication is that since we have a point load, the specified traction at the mid-point of the right end is infinite. We'll later discuss the effect of this in the ANSYS solution. Do keep in mind that there are no point loads in practice, it's just an idealization that can lead to weird behavior that we need to be aware of.

Hand Calculations

Now that we have reviewed the mathematical model for our problem, let's hold off diving into ANSYS just yet and first make some hand calculations of expected results. We'll use these hand calculations to check ANSYS results (like an expert engineer would!). In order to make the problem solvable by hand, we need to make additional assumptions. The ANSYS solution does not make these additional assumptions.

Additional Assumptions in Hand Calculations

#We'll simplify the right boundary condition. Instead of a point load, we'll assume that the load is distributed over the entire right boundary. So the traction condition at the right boundary becomes

Unknown macro: {latex}

[
\sigma_x = \frac

Unknown macro: {P}
Unknown macro: {w , t}

\: \: \tau_

Unknown macro: {xy}

= 0 \: \: \: at \: x = l
]

#Away from the lest and right ends, we expect a uni-axial state of stress with zero shear (OK, this is a bit of a leap of the imagination but it works). So we'll assume that everywhere

Unknown macro: {latex}

[
\tau_

Unknown macro: {xy}

= 0
]

Analytical Solution

With these additional assumptions in hand, we can solve for

Unknown macro: {latex}

$\sigma_y$

first and then

Unknown macro: {latex}

$\sigma_x$

.

Since the length is much larger than the width, we ignore end effects and neglect stress variations in the y direction.

Because they are free surfaces, stress in the y direction equals zero at both the top (y = H) and bottom (y = 0) surfaces of the bar.
If a surface is not constricted and can move freely, it can expand and contract without incurring stress.

All regions other than the left and right ends, experience pure uniaxial stress, thus we assume stress in the y direction and shear in the xy direction equals zero.

Unknown macro: {latex}

\begin

Unknown macro: {eqnarray}

\sigma_y = \tau_

Unknown macro: {xy}

= 0 \: \: at \: x \not = 0\nonumber
\end

Governing equations

Since we are assuming plane stress conditions, we can use the 2D version of the equilibrium equations. When the bar is in equilibrium, the stresses will satisfy the 2D equilibrium equations with no body forces.

Unknown macro: {latex}

\begin

Unknown macro: {eqnarray}
Unknown macro: {partial sigma_x over partial x}

+ {\partial \tau_

Unknown macro: {yx}

\over \partial y} + 0 = 0 \nonumber
{\partial \tau_

Unknown macro: {xy}

\over \partial x} +

Unknown macro: {partial sigma_y over partial y}

+ 0 = 0 \nonumber
\end

Boundary Conditions

The domain over which we'll solve the governing equations is a rectangle.

((Equation and Figure))

The boundary conditions are:

Unknown macro: {latex}

\begin

Unknown macro: {eqnarray}

\sigma_y = \tau_

Unknown macro: {xy}

= 0 \:\: at \: x \not = 0 \:\: (top, bottom, right \: sides)
U = V = 0 \:\: at \: x = 0 \:\: (left side)
\end

Solving the BVP

1. After applying the boundary condition (1) at the top and bottom edges:

Unknown macro: {latex}

\begin

Unknown macro: {eqnarray}
Unknown macro: {partial sigma_x over partial x}

= 0\nonumber
\nonumber
\nonumber
\sigma_y = 0\nonumber
\end


Therefore,

Unknown macro: {latex}

\begin

Unknown macro: {eqnarray}

\sigma_x = constant\nonumber
\end

 

Note that because of the additional assumptions, we are unable to satisfy the boundary condition (2). The additional assumptions were used to simplify the BVP and hand calcuation. Later on, boundary condition (2) will be used to explain edge effects in the Verification and Validation section of the tutorial.

Hand Calculations

Stress is found using the equation:

Unknown macro: {latex}

\begin

Unknown macro: {eqnarray}

\sigma_x = Force/Area \nonumber
\end

Because we chose to neglect end effects, we make a vertical cut in the bar to determine the stress. Therefore,

Unknown macro: {latex}

\begin

Unknown macro: {eqnarray}

\sigma_x = P/A = P/(H*Thickness) \nonumber
\end

 This is of course a well-known result. For this problem, we have

Unknown macro: {latex}

\begin

Unknown macro: {eqnarray}

\sigma_x = 2000/(10*1) = 200 \ N/mm^2 = 200 \ MPa
\nonumber
\end

Numerical Solution using FEA

In the numerical solution using FEA, we solve the 2D BVP directly by dividing the structure into small elements and approximating the solution for these small elements. Unlike the analytical approach, we do not assume that there is no variation in the y direction. Also, end effects are not neglected. The FEA solution is an approximate solution to the 2D BVP. The approximation gets better as the elements become smaller. In contrast, the analytical solution presented above is the exact solution to the 1D BVP obtained by making approximations to the 2D BVP. In other words, in the analytic solution, we have swapped the actual 2D BVP problem for a 1D BVP problem that we can solve in closed form. Both approaches have value in engineering and complement each other. We have checked that the FEA solution presented to you is reasonably accurate.

The following figure summarizes the contrasts between the analytical and numerical approaches.


Load FEA Solution obtained using ANSYS

As mentioned before, we are providing the FEA solution obtained using ANSYS so that you can focus on comparing the analytical and numerical solutions (which is the goal of this exercise). Without further ado, let's download the ANSYS solution and load it into ANSYS.

1. Download "Tensile Bar Demo.zip" by [clicking here|Tensile Bar - Pre-analysis and Start-up^Tensile Bar demo.zip]
The zip should contain a Tensile Bar Demo folder with the following contents:
- Tensile Bar Demo_files folder
- Tensile Bar Demo.wbpj
Please make sure both these are in the folder, otherwise the solution will not load into ANSYS properly. (Note: The solution provided was created using ANSYS workbench 13.0 release, there may be compatibility issues when attempting to open with other versions). Be sure to extract before use.

2. Double click "Tensile Bar Demo.wbpj" - This should automatically open ANSYS Workbench (you have to twiddle your thumbs a bit before it opens up). You will then be presented with the ANSYS solution in the project page.

A tick mark against each step indicates that that step has been completed.

3. To look at the results, double click on Results - This should bring up a new window (again you have to twiddle your thumbs a bit before it opens up).

4. On the left-hand side there should be an Outline toolbar. Look for Solution (A6).


We'll investigate the items listed under Solution in the next step of this tutorial.

[Go to Step 2: Numerical Results]

Go to all ANSYS Learning Modules

  • No labels