Constitutive Relations Video Lecture - Civil Engineering (CE)

Okay, we will continue with
laying out the elements of our
linear elliptic differential
equation in one dimension.
Okay?
Let's, let's look at our
differential equation again, so.
We've written it in the form d sigma
dx plus f function of x as we discussed,
is equal to 0 in 0, L.
Right.
Well, one thing you note of course is
that I'm using total derivatives here.
I'm not using partial derivatives, that's
because when we're doing one dimension,
doing things in one dimension,
we don't even have time.
There are,
there are no other variables, right?
So we're free do to this.
Of course when we go into multi,
multi-dimensional problems we will have
to use proper partial derivatives.
Right, so though I've been calling it
a partial differential equation sometimes
and a differential equation,
really they are one and the same thing for
this one dimensional problem
in a single variable.
Okay?
Just a few more things
I need to state here.
One thing you will note is that
I've been careful to say that
we have an open interval over which
the differential equation is specified.
Okay?
So, that is an open interval.
Right, and you recall what this means?
It says that we are interested in
looking at this differential equation
over the range 0 to L, right.
So, we have our domain here.
That is 0, that is L.
This is the domain over which we are
interested in looking at this differential
equation.
But importantly we are seeing that
the differential equation does
not apply at the point x equals 0 and
x equals f, right.
So open interval remember that
this implies the domain between.
X equals 0 and x equals L, but
importantly excluding the end points.
X equals 0 and
x equals L themselves, right?
This is important.
And what I want you to do is think of why
that maybe, why is it that we specify
the differential equation,
the partial differential equation or
just the differential equation,
in this case excluding the point 0 and L.
And to help you think about the answer,
recall that we do need to
know something about 0 and L.
If we don't have the differential
equation specified at those points we
must have something else.
All right?
Well, you probably get the answer.
The an, the, the answer is that we do
already have boundary conditions
specified at x equals 0 and x equals L.
If we were to insist that the differential
equation also held at the boundary points,
we would actually have a,
an over-constrained system.
Okay and this will often manifest
itself in the mathematical or
numerical solution of
problems where you by, by,
by essentially developing instabilities
in your solution mattered
or, or, or situations where you
just can't find a solution, okay?
So let me just state that.
This is because we have
boundary conditions.
Right.
And those boundary conditions
are either on the primal field itself, or
on effectively on u, x, okay, at.
The boundary points.
Okay?
And this is why the partial differential
equation is always specified
over an open interval.
This, by the way, is the case for
any differential equation, right?
And any, any, any,
any partial differential equation,
is not specified at its boundary points.
Likewise, if you had an ordinary
differential equation it would not
be specified at the end points
of the interval of interest.
Okay?
This is an important thing to note and,
and to recognize that it's not purely,
matter of pace.
It really means something.
Okay?
It means something mathematically and
physically as well.
Okay.
That's one thing I wanted to state.
The other thing I also
want to state is I also,
I also want to talk about is
the constitutive relations, right.
Constitutive relation.
And in this problem,
we specify the constitutive
relation as being sigma equals E times u,
x.
All right, in constitutive relation as,
the term suggests,
it tells us something more about
the constitution of the domain
of the problem, okay?
In this particular case,
we started out with a primal field u,
but then we,
introduced another field sigma.
Okay?
And, what we are seeing,
is that we know something more
about the domain of interest,
and in our case the domain
of interest is this one.
Sorry.
That is our domain of interest, right?
We're seeing that, sorry,
it's not x here, it's 0, that is x, okay.
We are saying we know something
more about the domain of this bar,
over the domain 0 and L,
by specifying this particular relation.
Now we are trying to set this
problem up using elasticity as a,
as an economical physical phenomenon,
right?
In that context sigma
is the stress okay and.
And we rewrite it here.
Sigma is the stress.
And it is related,
not to the displacement itself,
but to the gradient of the displacement,
which in this setting of linear or
linearized elasticity is the strain, okay?
The strain is sometime
given its own symbol.
And sometimes, or I guess pretty often
is denoted by epsilon, which is u comma x.
Okay, so in the context of this
particular problem that we're looking at,
the elasticity problem, we are seeing that
there is a relation between the stress and
this which is this field signal and the
strain, which is the displacement reading.
Okay?
And that relation as you can
see here is a linear one.
Okay, so when we see sigma
equals E times u comma x what
we have here is linear,
linearized elasticity.
Right, and
it's obvious why we say it's linear or
linearized because that's a linear
relation between the stress and strain.
In this setting, E is a physical constant,
that often gets called
the Young's modulus,
or you may simply call it the modulus
in a one dimensional setting.
Right, in, in a one dimensional setting,
there is only one modulus,
might as well call it the Young's modulus.
Okay?
So this is the other bit that
I wanted to state about this,
about the physical phenomenon
that we're using for
now as canonical setting for
the problem we're considering.
The probably very last thing I will do
is restate essentially the same problem
but in the context of
a different physical phenomenon.
Okay?
So let me state it here the 1D scalar,
because remember we said that in 1D
everything is essentially scalar.
So the 1D scalar 1D scalar
linear elliptic problem.
Also, models are heat or
diffusive mass transport.
Okay?
Let's suppose we are talking
about heat transport, okay?
In that case we would be,
we would state the problem as following.
We would say find u function of X,
which is a mapping, again,
from this domain 0, L, open interval
to the one dimensional real space.
Okay?
Given, I, I will attempt to use
different symbols here, okay?
And I, I,
I'm going to use different symbols.
But because we're talking
about the heat condition of,
of the heat transport problem we,
I, I will also state here that what
we're talking about when we talk about u,
the primary field, we,
we have the temperature in mind,
okay, for, for heat transport.
Okay, so define the temperature U,
which is mapping from 0,L open
interval to R1, given the following,
okay, given, u0 and ug or.
Now, I want to reserve t for traction.
So for
the other sort of boundary condition
that we could have in the context of heat
conduction, I'm going to use j, j bar.
Okay?
We're given these.
We have our forcing function,
f, and the constitutive relation.
Now in this case the constitutive relation
that I want to use is going to be for
the heat flux.
And want to use a different
symbol instead of sigma here.
Okay.
I'm going to use the, the symbol j for
the heat flux.
J equals minus kappa u comma x.
Okay, all right?
Such that.
We have the following, okay?
Minus d,
dx of j equals f
of x in (0,L),
okay?
All right.
Now let,
let me put down the boundary conditions.
And the boundary conditions are.
'Kay, u at zero
equals u naught and
either u at L equals ug.
Right?
Now both these, u at 0 being u naught, and
u at L being ug would be given
temperature conditions, all right?
Or at X equals L, we may
also have the condition
that j at L, okay,
equals minus j bar.
Okay?
All right.
Right, so this is how we may specify
the heat conduction problem.
Now in doing this observe that
what we have here is the,
is to be be thought of as the divergence
of, of, of the heat flux.
Okay, and in fact G.
In this case,
would be the heat flux.
Okay?
What I've circled here would
be the divergence of the heat flux.
Okay?
In the con, in the context of heat
conduction, what we would be seeing
would be that the, so, sorry this would be
the negative diversions of the heat flux.
Okay.
So the minus diversions will be heat flux.
What we would be seeing is that the
negative of the diversion of the heat flux
in the context of heat conduction is
the amount of heat you're supplying to
a particular little volume element in,
in this domain of interest.
Okay?
I already said that you, you,
that the conditions on you would be,
would be temperature boundary conditions.
So let me specify that also, okay?
So the Dirichlet boundary
conditions would imply in this case,
temperature on great conditions.
Okay?
And this, and finally,
the Neumann boundary condition
would imply a heat flux on recondition.
Okay?
In particular, this condition
that j at L equals minus j bar.
Okay?
This would have the following
significance that j bar,
as I've written it here would
be the heat influx at x equals L.
Okay?
So, if you look at this problem that we've
just described here rather quickly and
you just make substitutions, right?
You just make substitutions
instead of you know, instead of j,
you go back to sigma instead of kappa,
you go back to e.
And instead of j bar, you go back to t.
Okay?
You will see that the problems
are in fact very similar.
Mathematically, it is
exactly the same problem.
But we know that this is a problem
of steady state heat transport or
diffusive mass transport.
So, so, so this is important to,
to note that there
is a multitude of problems
that are described by the same
differential equation or
the same partial differential equation.
It's just a matter of redefining our
fields of interest boundary conditions
of interest and
giving them different ascribing to them
different physical interpretations.
Okay?
And because we're interested
in developing in this class
finite element methods that
apply to a range of problems.
We will tend to focus more upon
the underlying differential
equation of interest.
And, of, we will repeatedly make
connections with the physical problems but
we will also bear in mind that rather
than thinking of the method as
a problem to be developed
as a problem applicable.
Only to structure mechanics and
only to heat conduction,
it truly is a lot more general than that.
All right?
Okay.
So this is where we were and this is where
we are going to stop for this segment.