## Calculation Module » History » Version 14

*Andreas Hauffe, 09/15/2021 02:43 PM *

1 | 1 | Andreas Hauffe | h1. Calculation Module |
---|---|---|---|

2 | |||

3 | h2. Generall |
||

4 | |||

5 | 10 | Andreas Hauffe | The calculation module is used to calculate loads or strains, as well as stresses within the specified laminate according to the classical laminated plate theory (!cite{Redd2003},!cite{Schu2004}). In order to be able to directly output the effects of the laminate stacking sequence on the stiffness matrix, the loads or strains and the stresses of the layers, the calculations are updated automatically. Among other things, hygrothermal effects can also be taken into account !cite{Barth2009}. |

6 | 1 | Andreas Hauffe | |

7 | 10 | Andreas Hauffe | During a calculation, the system of linear equations listed below is solved with the help of the Gaussian elimination method. The calculated unknowns are entered in the fields of the loads and strains. In addition, the stresses and strains in each laminate layer are calculated. These can be viewed in the internal forces subwindow in the local fiber coordinate system. |

8 | 1 | Andreas Hauffe | |

9 | 14 | Andreas Hauffe | h2. Window Layout |

10 | 1 | Andreas Hauffe | |

11 | 5 | Andreas Hauffe | p=. {{thumbnail(berechnung.png,size=500,title=calculation window)}} |

12 | 1 | Andreas Hauffe | |

13 | 11 | Andreas Hauffe | h3. 1 - ABD matrix |

14 | 1 | Andreas Hauffe | |

15 | 10 | Andreas Hauffe | Based on the stiffness matrices of the individual layers, the membrane stiffness matrix $\mathbf{A}$, coupling stiffness matrix $\mathbf{B}$ and bending stiffness matrix $\mathbf{D}$ of the laminate are automatically calculated with repect to the classical laminated plate theory and the three matrices are shown separately and color-coded within the ABD matrix. |

16 | 1 | Andreas Hauffe | |

17 | 10 | Andreas Hauffe | h3. 2 - Internal Loads, Strain and hygrothermal loads |

18 | 1 | Andreas Hauffe | |

19 | 11 | Andreas Hauffe | Both membrane internal loads and moments as well as strain and curvatures can be specified as loads on the defined laminate. Either the load or the strain can be specified in the input fields under 2a and 2b for each coordinate direction. In addition, a temperature difference and the percentage change in the relative humidity can be specified in the fields in window section 2c. According to the classical laminated plate theory applies to mechanical and hygrothermal loads and strains |

20 | 1 | Andreas Hauffe | |

21 | 6 | Andreas Hauffe | $$\begin{pmatrix} \underline{n} \\ \underline{m} \end{pmatrix}_{mech} = \begin{bmatrix} \mathbf{A} & \mathbf{B} \\ \mathbf{B} & \mathbf{D} \end{bmatrix} \begin{pmatrix} \underline {\varepsilon} \\ \underline {\kappa} \end{pmatrix} - \begin{pmatrix} \underline{n} \\ \underline{m} \end{pmatrix}_{hygrotherm}$$ |

22 | 1 | Andreas Hauffe | |

23 | If no effects of the hygrothermal loads are to be considered, the fields for the temperature difference and the percentage humidity difference must be filled with the values zero. This view is the standard setting in the calculation module of eLamX². |
||

24 | |||

25 | 11 | Andreas Hauffe | h3. 3 - hygrothermal internal loads |

26 | 1 | Andreas Hauffe | |

27 | 11 | Andreas Hauffe | At this point, the resulting hygrothermal internal loads are output based on the specified temperature difference and the percentage difference in humidity on the entire laminate. They are not input data. The selection of the check boxes for mechanical loads and the resulting distortions has no influence on the hygrothermal cutting loads. They result from the direction-dependent thermal conductivity and swelling coefficients of the individual layers as well as their stiffness matrices and the specified temperature and humidity difference. |

28 | 1 | Andreas Hauffe | |

29 | h3. 4 - Stress and strain distribution button |
||

30 | |||

31 | 12 | Andreas Hauffe | This button can be used to display the stress and strain distribution within the laminate in local and global coordinates for the coordinate axes of the laminate plane. The output is only qualitative. The specified coordinate directions relate to the fiber angle, the specified x-value reflecting the size in the 0° direction and the y-value reflecting the size in the 90° direction of the laminate. The fiber orientations of the individual layers are shown by the hatching. Horizontal lines correspond to zero degrees and vertical lines to 90 degrees fiber angles. If the fiber angles of adjacent layers are very close together, the hatching is output in different colors. |

32 | 1 | Andreas Hauffe | |

33 | 8 | Andreas Hauffe | p=. {{thumbnail(spannungsverteilung.png,size=500,title=call up 3D failure body of the single layer)}} |

34 | 1 | Andreas Hauffe | |

35 | h3. 5 - Delete button |
||

36 | |||

37 | 12 | Andreas Hauffe | Pressing the button deletes the specified loads and the calculated unknown internal loads and distortions, as well as the calculated stress distribution. |

38 | 1 | Andreas Hauffe | |

39 | 12 | Andreas Hauffe | h3. 6 - Representation of the strains and curvatures on the basis of a square plate |

40 | 1 | Andreas Hauffe | |

41 | This button can be used to open a window in which the calculated strains and curvatures are visualized on a square plate. This is used for a better understanding of the coupling effects within the ABD matrix. |
||

42 | |||

43 | 8 | Andreas Hauffe | p=. {{thumbnail(darstellung_eps_kappa.png, size=500, title=call up 3D failure body of the single layer)}} |

44 | 1 | Andreas Hauffe | |

45 | 12 | Andreas Hauffe | h3. 7 - Layer data in the local fiber coordinate system |

46 | 1 | Andreas Hauffe | |

47 | 12 | Andreas Hauffe | After a calculation, the stresses or strains of each layer of the composite are calculated in this part of the window, depending on the selection. The evaluation is carried out on the top and bottom of each individual layer. The stresses in each layer are compared with their strengths and a reserve factor is calculated on the top and bottom of the layer using the selected failure criterion. In addition, the type of failure to be expected is specified with regard to the failure criterion selected for each layer. |

48 | 3 | Andreas Hauffe | |

49 | 8 | Andreas Hauffe | p=. {{thumbnail(aufruf_3d-versagenskoerper_einzelschicht.png, size=500, title=calling up 3D failure body of the single layer)}} |

50 | 3 | Andreas Hauffe | |

51 | In addition to displaying the failure criteria using the reserve factors, right-click on a layer to call up the Failure Criteria module and display the failure body of the individual layer together with the stress state in the stress area (see figure). The current stress on the composite layer is entered as a red point. If the calculated stress vector lies within the selected failure body, only its line is visible. This makes it clear which tension combination prevails in each layer. |
||

52 | |||

53 | h2. Laminate information / engineering constants |
||

54 | |||

55 | 8 | Andreas Hauffe | p=. {{thumbnail(laminatinformationen.png, size=500, title=call up 3D failure body of the single layer)}} |

56 | 3 | Andreas Hauffe | |

57 | This window can be opened by right-clicking on a laminate and then selecting engineering constants. |
||

58 | |||

59 | 11 | Andreas Hauffe | h3. 1 - ABD matrix |

60 | 3 | Andreas Hauffe | |

61 | 13 | Andreas Hauffe | Based on the stiffness matrices of the individual layers, the membrane stiffness matrix $\mathbf{A}$, coupling stiffness matrix $\mathbf{B} $ and bending stiffness matrix $\mathbf{D} $ of the laminate are automatically calculated and color-coded when the calculation module is called up on the basis of classical laminated plate theory. This matrix is shown in the form of a table so that all values can be copied out. |

62 | 3 | Andreas Hauffe | |

63 | 13 | Andreas Hauffe | h3. 2 - Compliance matrix of the laminate |

64 | 3 | Andreas Hauffe | |

65 | The inverse of the stiffness matrix of the laminate is displayed here. It may be that for symmetrical laminates, terms of the inverse coupling stiffness matrix (b) are different from zero, but very small. This results from numerical inaccuracies in the double precision calculation within the program. This matrix is shown in the form of a table so that all values can be copied out. |
||

66 | |||

67 | 13 | Andreas Hauffe | h3. 3 - engineering constants of the laminate |

68 | 3 | Andreas Hauffe | |

69 | 13 | Andreas Hauffe | The engineering constants of the laminate are also automatically calculated according to !cite{Schu2004} for membrane and bending loads. This is done both with and without consideration of Poisson effect. In the case of asymmetrical laminates, where a coupling between membrane and bending deformation is present, the engineering constants have little meaning and should be used with caution. |

70 | 3 | Andreas Hauffe | |

71 | 13 | Andreas Hauffe | h3. 4 - Coefficient of expansion of the laminate for heat and moisture |

72 | 3 | Andreas Hauffe | |

73 | 13 | Andreas Hauffe | In this part of the window the thermal expansion coefficients $\alpha^T_i$ and the swelling coefficient $\beta_i$ of the total laminate are output according to the system of units that was selected for specifying the values of the individual layer. Here, too, the definition according to !cite{Schu2004} applies and thus |

74 | 8 | Andreas Hauffe | \begin{equation} |

75 | \left( |
||

76 | \begin{array}{c} |
||

77 | \epsilon_x \\ |
||

78 | \epsilon_y \\ |
||

79 | \gamma_{xy} |
||

80 | \end{array} |
||

81 | \right)_{T+H} = |
||

82 | \Delta T |
||

83 | \left( |
||

84 | \begin{array}{c} |
||

85 | \alpha^T_x \\ |
||

86 | \alpha^T_y \\ |
||

87 | \alpha^T_{xy} |
||

88 | \end{array} |
||

89 | \right) + |
||

90 | \Delta c |
||

91 | \left( |
||

92 | \begin{array}{c} |
||

93 | \beta_x \\ |
||

94 | \beta_y \\ |
||

95 | 1 | Andreas Hauffe | \beta_{xy} |

96 | \end{array} |
||

97 | 8 | Andreas Hauffe | \right) |

98 | 3 | Andreas Hauffe | \end{equation} |

99 | |||

100 | h2. Notes on the engineering constants |
||

101 | |||

102 | 13 | Andreas Hauffe | The engineering constants can be determined for two types of stress. On the one hand for a membrane stress and on the other hand for a bending stress. In both cases, this is possible with and without Poisson effect. |

103 | 3 | Andreas Hauffe | |

104 | 13 | Andreas Hauffe | The approach is presented on the basis of the explanations in !cite{Schu2004} for symmetrical laminates. For the membrane state, the ABD matrix becomes |

105 | 8 | Andreas Hauffe | \begin{equation} |

106 | \left( |
||

107 | \begin{array}{c} |
||

108 | 3 | Andreas Hauffe | n \\ |

109 | 8 | Andreas Hauffe | m |

110 | \end{array} |
||

111 | \right) = |
||

112 | \left[ |
||

113 | \begin{array}{cc} |
||

114 | A & B \\ |
||

115 | B & D |
||

116 | 1 | Andreas Hauffe | \end{array} |

117 | 8 | Andreas Hauffe | \right] |

118 | \left( |
||

119 | \begin{array}{c} |
||

120 | \epsilon \\ |
||

121 | \kappa |
||

122 | \end{array} |
||

123 | \right) |
||

124 | \end{equation} |
||

125 | 13 | Andreas Hauffe | only the $A$-matrix is used, which is feasible for symmetrical laminates ($B = 0$). This results in |

126 | 8 | Andreas Hauffe | \begin{equation} |

127 | \left( |
||

128 | \begin{array}{c} |
||

129 | 3 | Andreas Hauffe | n_x \\ |

130 | n_y \\ |
||

131 | 8 | Andreas Hauffe | n_{xy} |

132 | \end{array} |
||

133 | \right) = |
||

134 | \left[ |
||

135 | \begin{array}{ccc} |
||

136 | A_{11} & A_{12} & A_{13} \\ |
||

137 | A_{21} & A_{22} & A_{23} \\ |
||

138 | A_{31} & A_{32} & A_{33} |
||

139 | 1 | Andreas Hauffe | \end{array} |

140 | 8 | Andreas Hauffe | \right] |

141 | \left( |
||

142 | \begin{array}{c} |
||

143 | \epsilon_{x} \\ |
||

144 | \epsilon_{y} \\ |
||

145 | \gamma_{xy} |
||

146 | \end{array} |
||

147 | \right) |
||

148 | \end{equation} |
||

149 | 13 | Andreas Hauffe | The aim is now to come up with an equation for the stress $\sigma_{x}$, which corresponds to the Hooke's law $\sigma=E \epsilon$. To do this, the equation is divided by the total thickness of the laminate. This corresponds to a homogenization of the material and one obtains: |

150 | 8 | Andreas Hauffe | \begin{equation} |

151 | \left( |
||

152 | \begin{array}{c} |
||

153 | \sigma_x \\ |
||

154 | \sigma_y \\ |
||

155 | \tau_{xy} |
||

156 | \end{array} |
||

157 | \right) = \frac{1}{t_{ges}} |
||

158 | \left[ |
||

159 | \begin{array}{ccc} |
||

160 | A_{11} & A_{12} & A_{13} \\ |
||

161 | A_{21} & A_{22} & A_{23} \\ |
||

162 | 1 | Andreas Hauffe | A_{31} & A_{32} & A_{33} |

163 | 8 | Andreas Hauffe | \end{array} |

164 | \right] |
||

165 | \left( |
||

166 | \begin{array}{c} |
||

167 | \epsilon_{x} \\ |
||

168 | \epsilon_{y} \\ |
||

169 | \gamma_{xy} |
||

170 | \end{array} |
||

171 | \right) \ . |
||

172 | \end{equation} |
||

173 | 13 | Andreas Hauffe | This equation corresponds to a single anisotropic layer across the entire thickness. There are two ways to come up with an equation according to the Hooke's law. On the one hand, $\epsilon_{y}$ and $\gamma_{xy}$ are set to zero without Poisson effect, and the following equations are obtained: |

174 | 8 | Andreas Hauffe | \begin{equation} |

175 | \left( |
||

176 | \begin{array}{c} |
||

177 | \sigma_x \\ |
||

178 | \sigma_y \\ |
||

179 | \tau_{xy} |
||

180 | \end{array} |
||

181 | \right) = \frac{1}{t_{ges}} |
||

182 | \left[ |
||

183 | \begin{array}{ccc} |
||

184 | A_{11} & A_{12} & A_{13} \\ |
||

185 | A_{21} & A_{22} & A_{23} \\ |
||

186 | A_{31} & A_{32} & A_{33} |
||

187 | \end{array} |
||

188 | 1 | Andreas Hauffe | \right] |

189 | 8 | Andreas Hauffe | \left( |

190 | \begin{array}{c} |
||

191 | \epsilon_{x} \\ |
||

192 | 3 | Andreas Hauffe | 0 \\ |

193 | 8 | Andreas Hauffe | 0 |

194 | \end{array} |
||

195 | \right) \ . |
||

196 | \end{equation} |
||

197 | \begin{equation} |
||

198 | \sigma_x = \frac{A_{11}}{t_{ges}}\epsilon_{x} |
||

199 | \end{equation} |
||

200 | 13 | Andreas Hauffe | In order to obtain a modulus of elasticity with Poisson effect, the equation has to be transformed and $\sigma_{y}$ and $\tau_{xy} $ are set to zero. |

201 | 8 | Andreas Hauffe | \begin{equation} |

202 | \left( |
||

203 | \begin{array}{c} |
||

204 | \epsilon_{x} \\ |
||

205 | \epsilon_{y} \\ |
||

206 | \gamma_{xy} |
||

207 | \end{array} |
||

208 | \right) |
||

209 | = t_{ges} |
||

210 | \left[ |
||

211 | \begin{array}{ccc} |
||

212 | A_{11} & A_{12} & A_{13} \\ |
||

213 | A_{21} & A_{22} & A_{23} \\ |
||

214 | A_{31} & A_{32} & A_{33} |
||

215 | \end{array} |
||

216 | 1 | Andreas Hauffe | \right]^{-1} |

217 | 8 | Andreas Hauffe | \left( |

218 | \begin{array}{c} |
||

219 | \sigma_x \\ |
||

220 | 3 | Andreas Hauffe | 0 \\ |

221 | 8 | Andreas Hauffe | 0 |

222 | \end{array} |
||

223 | 1 | Andreas Hauffe | \right) |

224 | 8 | Andreas Hauffe | \ . |

225 | \end{equation} |
||

226 | \begin{equation} |
||

227 | \sigma_x = \frac{1}{(A^{-1})_{11} \cdot t_{ges}}\epsilon_{x} |
||

228 | \end{equation} |
||

229 | 13 | Andreas Hauffe | In this way, all engineering constants can now be determined. The constants with Poisson effect result from: |

230 | 1 | Andreas Hauffe | \begin{eqnarray} |

231 | 8 | Andreas Hauffe | E_x & = & \frac{1}{(A^{-1})_{11} \cdot t_{ges}} \\ |

232 | E_y & = & \frac{1}{(A^{-1})_{22} \cdot t_{ges}} \\ |
||

233 | G_{xy} & = & \frac{1}{(A^{-1})_{66} \cdot t_{ges}} \\ |
||

234 | \nu_{xy} & = & -\frac{(A^{-1})_{12}}{(A^{-1})_{11}} \\ |
||

235 | \nu_{yx} & = & -\frac{(A^{-1})_{12}}{(A^{-1})_{22}} |
||

236 | 1 | Andreas Hauffe | \end{eqnarray} |

237 | 13 | Andreas Hauffe | The constants without Poisson effect result from: |

238 | 8 | Andreas Hauffe | \begin{eqnarray} |

239 | E_x & = & \frac{A_{11}}{t_{ges}} \\ |
||

240 | E_y & = & \frac{A_{22}}{t_{ges}} \\ |
||

241 | G_{xy} & = & \frac{A_{66}}{t_{ges}} \\ |
||

242 | \end{eqnarray} |
||

243 | 3 | Andreas Hauffe | |

244 | 13 | Andreas Hauffe | If the Poisson effect is not taken into account, Poisson's ratios are not useful and are therefore not included. |

245 | 3 | Andreas Hauffe | |

246 | 13 | Andreas Hauffe | These homogenized material parameters are only feasible for membrane loads. An equivalent procedure using the D-matrix is necessary for bending, since the position of the layers has a great influence in the case of bending stresses in contrast to membrane stresses. The derivation of the formulas is taken from !cite{Schu2004}. |

247 | 3 | Andreas Hauffe | |

248 | The analogy to the beam is used for bending. Since the laminate board has a rectangular cross-section, it is compared with the bending elasticity law of a rectangular beam and the following applies: |
||

249 | 8 | Andreas Hauffe | \begin{equation} |

250 | M = -EI \cdot w'' = - E \frac{b t^3}{12} \cdot w'' \rightarrow m = -E \frac{t^3}{12} \cdot w'' \ . |
||

251 | \end{equation} |
||

252 | 3 | Andreas Hauffe | |

253 | 13 | Andreas Hauffe | In addition, $w'' = -\kappa$ applies. |

254 | 3 | Andreas Hauffe | |

255 | 13 | Andreas Hauffe | Bending loads with the $D$-matrix are described for the laminate. |

256 | 8 | Andreas Hauffe | \begin{equation} |

257 | \left( |
||

258 | \begin{array}{c} |
||

259 | 3 | Andreas Hauffe | m_x \\ |

260 | 1 | Andreas Hauffe | m_y \\ |

261 | 8 | Andreas Hauffe | m_{xy} |

262 | \end{array} |
||

263 | \right) = |
||

264 | \left[ |
||

265 | \begin{array}{ccc} |
||

266 | D_{11} & D_{12} & D_{13} \\ |
||

267 | D_{21} & D_{22} & D_{23} \\ |
||

268 | D_{31} & D_{32} & D_{33} |
||

269 | \end{array} |
||

270 | \right] |
||

271 | \left( |
||

272 | \begin{array}{c} |
||

273 | \kappa_{x} \\ |
||

274 | \kappa_{y} \\ |
||

275 | \kappa_{xy} |
||

276 | \end{array} |
||

277 | \right) |
||

278 | \end{equation} |
||

279 | 13 | Andreas Hauffe | Now the approach is equivalent to stressing the membrane. Assuming no Poisson effect results in: |

280 | 8 | Andreas Hauffe | \begin{equation} |

281 | \left( |
||

282 | \begin{array}{c} |
||

283 | 3 | Andreas Hauffe | m_x \\ |

284 | m_y \\ |
||

285 | 8 | Andreas Hauffe | m_{xy} |

286 | \end{array} |
||

287 | \right) = |
||

288 | \left[ |
||

289 | \begin{array}{ccc} |
||

290 | 1 | Andreas Hauffe | D_{11} & D_{12} & D_{13} \\ |

291 | 8 | Andreas Hauffe | D_{21} & D_{22} & D_{23} \\ |

292 | D_{31} & D_{32} & D_{33} |
||

293 | \end{array} |
||

294 | \right] |
||

295 | 1 | Andreas Hauffe | \left( |

296 | 8 | Andreas Hauffe | \begin{array}{c} |

297 | \kappa_{x} \\ |
||

298 | 3 | Andreas Hauffe | 0 \\ |

299 | 8 | Andreas Hauffe | 0 |

300 | \end{array} |
||

301 | 1 | Andreas Hauffe | \right) |

302 | \end{equation} |
||

303 | 8 | Andreas Hauffe | \begin{equation} |

304 | m_x = D_{11} \kappa_x |
||

305 | \end{equation} |
||

306 | \begin{equation} |
||

307 | E_{x,b} = \frac{12}{t_{ges}^3}D_{11} |
||

308 | \end{equation} |
||

309 | 3 | Andreas Hauffe | |

310 | 13 | Andreas Hauffe | Thus, all engineering constants without Poisson effect result in: |

311 | 8 | Andreas Hauffe | \begin{eqnarray} |

312 | E_{x,b} & = & \frac{12}{t_{ges}^3}D_{11} \\ |
||

313 | E_{y,b} & = & \frac{12}{t_{ges}^3}D_{22} \\ |
||

314 | G_{x,b} & = & \frac{12}{t_{ges}^3}D_{66} \\ |
||

315 | \end{eqnarray} |
||

316 | 13 | Andreas Hauffe | and with Poisson effect: |

317 | 8 | Andreas Hauffe | \begin{eqnarray} |

318 | E_{x,b} & = & \frac{12}{(D^{-1})_{11} \cdot t_{ges}^3} \\ |
||

319 | E_{y,b} & = & \frac{12}{(D^{-1})_{22} \cdot t_{ges}^3} \\ |
||

320 | G_{x,b} & = & \frac{12}{(D^{-1})_{66} \cdot t_{ges}^3} \\ |
||

321 | \end{eqnarray} |
||

322 | 3 | Andreas Hauffe | |

323 | 13 | Andreas Hauffe | At this point, too, Poisson's ratios are not very useful. |

324 | 3 | Andreas Hauffe | |

325 | 13 | Andreas Hauffe | This approach generally only applies to symmetrical laminates. In the case of asymmetrical laminates, the inverse of the entire $ABD$ matrix must be formed and then the corresponding terms must be used. |

326 | 3 | Andreas Hauffe | |

327 | !putbib{} |