Derive the generalized constitutive matrix \hat{D}{}' for an arbitrary position of the beam axis.
Derive the generalized constitutive matrix \hat{D}{}' for an arbitrary position of the beam axis.
Let us assume that the {x}' axis does not coincide with the neutral axis and {y}' , {z}' are not principal inertia axes. Substituting the expression for the axial stress \varepsilon _{{x}'} (Eq.(4.26)) into the integral for the resultant stress vector (Eq.(4.32)gives
\hat{\sigma {}'}=\int \int_{A}^{}\begin{Bmatrix}E\left ( \frac{\partial u_{0{}'}}{\partial {x}'} +z{}'\frac{\partial\theta _{x{}'}}{\partial {x}'}-{y}'\frac{\partial \theta _{x{}'}}{\partial {x}'}\right )\\ G_{{y}'}\left ( \frac{\partial {v}'_{c}}{\partial x}-\theta _{{x}'} \right )\\ G_{{z}'}\left ( \frac{\partial {w}'_{c}}{\partial x}-\theta _{{y}'} \right )\\ {z}'E\left ( \frac{\partial u_{0{}'}}{\partial {x}'} +z{}'\frac{\partial \theta _{x{}'}}{\partial {x}'}-{y}'\frac{\partial \theta _{x{}'}}{\partial {x}'}\right )\\ \end{Bmatrix}dA
where D_{t} is defined in Eq.(4.33). Simple algebra gives
\hat{{\sigma}' }=\hat{{D}'}\hat{{\varepsilon }'} with \hat{{D}'}=\int \int_{A}^{}\begin{bmatrix}E &0 &0 &{z}'E &{-y}'E &0 \\ 0 &G_{{y}'} &0 &0 &0 &0 \\ 0 &0 &G_{{z}'} &0 &0 &0 \\ {z}'E &0 &0 &{z}'^{2}E &-{y}'{z}'E &0 \\ -{y}'E &0 &0 &-{y}'{z}'E &{y}'^{2}E &0 \\ 0 &0 &0 &0 &0 &D_{t} \end{bmatrix}
with {\hat{\varepsilon} }' given by Eq.(4.27b).
We see that \hat{{D}'} is now a full (still symmetric) matrix. This implies that axial and flexural effects are coupled, i.e. an axial force induces flexural (bending and shear) effects and viceversa. The off-diagonal terms in \hat{{D}'} vanish if {x}' is the neutral axis and {y}' , {z}' are the principal inertial axes .
In practice, either the full form of \hat{{D}'} given above or the diagonal form of Eq.(4.35a) can be used and both yield identical results. The simpler diagonal form requires the “a priori” computation of the position of the neutral axis and the principal inertia axes.
Above considerations do not affect the value of the torsional stiffness \hat{D_{t}}.
\\(Eq.(4.26)): {\varepsilon }'=\begin{Bmatrix}\varepsilon _{{x}'}\\ \gamma _{{x}'{y}'}\\ \gamma _{{x}'{z}'}\end{Bmatrix}=
\overset{\begin{Bmatrix}\frac{\partial {u}'_{0}}{\partial x}\\ 0\\ 0\end{Bmatrix}+\begin{Bmatrix}{z}'\frac{\partial \theta _{{y}'}}{\partial {x}'}\\ 0\\ \frac{\partial {w}'_{c}}{\partial {x}'+\theta _{{y}'}}\end{Bmatrix}+\begin{Bmatrix}-{y}'\frac{\partial \theta _{{z}'}}{\partial {x}'}\\ \frac{\partial {v}'_{c}}{\partial {x}'}-\theta _{{z}'}\\ 0\end{Bmatrix}+\begin{Bmatrix}0\\ \left ( \frac{\partial w}{\partial {y}'}-\left ( {z}'-{z}'_{c} \right ) \right )\frac{\partial \theta _{\hat{x}}{}'}{\partial x}\\ \left ( \frac{\partial w}{\partial {z}'}-\left ( {y}'-{y}'_{c} \right ) \right )\frac{\partial \theta _{\hat{x}}{}'}{\partial x}\end{Bmatrix}}{\begin{matrix}axial & bending in &bending in && & & &free torsion \\ & plane{x}'{z}' &plane{x}'{y}' & & & & & \end{matrix}}
(Eq.(4.32)): \hat{\sigma}{}'=\iint_{A}^{}\begin{Bmatrix}E\varepsilon _{{x}'}\\ G_{{y}'}\left ( \frac{\partial {v}'_{c}}{\partial {x}'}-\theta _{{x}'} \right )\\ G_{{z}'}\left ( \frac{\partial {w}'_{c}}{\partial {z}'}-\theta _{{y}'} \right )\\ {z}'E\varepsilon _{{x}'}\\ -{y}'E\varepsilon _{{x}'}\\ D_{t}\end{Bmatrix}dA
(Eq.(4.33)): D_{t}=\begin{bmatrix}G_{{z}'}\left ( \frac{\partial w}{\partial {z}'}+{y}'-{y}'_{c} \right )\left ( {y}'-{y}'_{c}\right )-G_{{y}'}\left ( \frac{\partial w}{\partial {y}'}+{z}'-{z}'_{c} \right ) \left ( {z}'-{z}'_{c} \right )\end{bmatrix}\frac{\partial \theta _{\hat{x}{}'}}{\partial {x}'}
(Eq.(4.27b)): \hat{\varepsilon} {}'=\begin{bmatrix}\frac{\partial {u}'_{c}}{\partial {x}'},\left ( \frac{\partial {v}'_{c}}{\partial {x}'}-\theta _{{z}'} \right ),\left ( \frac{\partial {w}'_{c}}{\partial {x}'}-\theta _{{y}'} \right),\frac{\partial \theta _{{y}'}}{\partial {x}'}, \frac{\partial \theta _{{z}'}}{\partial {x}'},\frac{\partial\theta _{{\hat{x}}'}}{\partial {x}'}\end{bmatrix}^{T}
(Eq.(4.35a)): \hat{D}{}'=\begin{bmatrix}\hat{D}_{a} &\vdots & 0 &0 &0 &0 &\vdots &0 \\ \cdots &\cdots &\cdots &\cdots &\cdots &\cdots &\cdots &\cdots \\ 0& \vdots &\hat{D}_{s_{{y}'}} &0 &0 &0 &\vdots &0 \\ 0 &\vdots &0 &\hat{D}_{s_{{z}'}} &0 &0 &\vdots &0 \\ 0 &\vdots &0 & 0 &\hat{D}_{b_{{y}'}} & 0 &\vdots &0 \\ 0 &\vdots &0 &0 &0 &\hat{D}_{b_{{z}'}} &\vdots &0 \\ \cdots &\cdots &\cdots &\cdots &\cdots &\cdots &\cdots &\cdots \\ 0 &\vdots &0 &0 &0 &0 &\vdots &\hat{D_{t}} \end{bmatrix}=\begin{bmatrix}\hat{D}_{a} &0 &0 \\ 0 &\hat{D{}'}_{f} & 0\\ 0&0 &\hat{D}_{t} \end{bmatrix}