1 : * @LAKAPPA PROCEDUR KK2000 99/11/30 21:15:11 3720 2 : * 3 : ************************************************************************ 4 : * PROCEDURA @LAKAPPA 5 : ************************************************************************ 6 : * 7 : * Descrizione : trova i coefficenti correttivi del taglio K13 e K23, 8 : * lamina per lamina, e fornisce in uscita una tabella 9 : * identica a quella di INPUT, ma con i moduli G13 e G23 10 : * di ogni lamina opportunamente corretti. 11 : * I valori sono moltiplicati per 6/5 in modo da annullare 12 : * la correzione nel caso di materiali isotropi, gi… 13 : * apportata internamente a procedure di calcolo 14 : * 15 : * Sintassi : TG = @LAKAPPA TAB_COMP 16 : * 17 : * Autore/Data : Miliozzi A. 10/04/1997 18 : * 19 : * Revisione : 20 : * 21 : *********************************************************************** 22 : * 23 : *--- INPUT 24 : *========= 25 : * 26 : * TAB_COMP tavola caratteristica dei laminati compositi 27 : * 28 : *--- OUTPUT 29 : *========== 30 : * 31 : * TG tavola caratteristica dei laminati compositi con i valori 32 : * G13 e G23 corretti 33 : * 34 : * 35 : ********************************************************************** 36 : * 37 : DEBPROC @LAKAPPA TCOMP*TABLE ; 38 : * 39 : NZON = (DIME TCOMP) - 2 ; 40 : N1 = 1 ; 41 : REPETER LOOP1 NZON ; 42 : * 43 : TG = TCOMP ; 44 : N2 = 1 ; 45 : PD1 = PROG ; PD2 = PROG ; 46 : PZS = PROG ; PZI = PROG ; 47 : NZN1 = 0; DZN1 = 0 ; 48 : NZN2 = 0; DZN2 = 0 ; 49 : NLAM = (DIME TCOMP.N1) - 6 ; 50 : * 51 : REPETER LOOP2 NLAM ; 52 : * 53 : EX1 = TCOMP.N1.N2.EXC ; 54 : EP1 = TCOMP.N1.N2.EPA ; 55 : AN1 = TCOMP.N1.N2.ANG ; 56 : ID1 = TCOMP.N1.N2.IDM ; 57 : * 58 : TCM = TCOMP.MCOL ; 59 : Y1 = TCM.ID1.YG1 ; Y2 = TCM.ID1.YG2 ; 60 : N12 = TCM.ID1.NU12 ; G1 = TCM.ID1.G12 ; 61 : * 62 : BETA = AN1 ; 63 : * 64 : ******************************************************* 65 : * Ricavo le coordinate estreme di ogni lamina 66 : ******************************************************* 67 : * 68 : ZS = EX1+(EP1 / 2) ; 69 : ZI = EX1-(EP1 / 2) ; 70 : POLD=PZS; 71 : PZS = INSERE POLD N2 ZS ; DETR POLD ; 72 : POLD=PZI; 73 : PZI = INSERE POLD N2 ZI ; DETR POLD ; 74 : * 75 : ******************************************************** 76 : * Calcolo i valori delle costanti ingegneristiche 77 : * nel sistema di riferimento globale 78 : ******************************************************** 79 : * 80 : CB = COS BETA; SB = SIN BETA; 81 : YX = 1/(((CB**4)/Y1)+((1/G1)-((2*N12/Y1))*(SB**2)*(CB**2)) 82 : +((SB**4)/Y2)); 83 : YY = 1/(((SB**4)/Y1)+((1/G1)-((2*N12/Y1))*(SB**2)*(CB**2)) 84 : +((CB**4)/Y2)); 85 : NXY = (((N12/Y1)*((SB**4)+(CB**4)))- 86 : (((1/Y1)+(1/Y2)-(1/G1))*(SB**2)*(CB**2)))*YX ; 87 : GXY = 1/((2*((2/Y1)+(2/Y2)+(4*N12/Y1)-(1/G1))*(SB**2)*(CB**2)) 88 : +(((SB**4)+(CB**4))/G1)) ; 89 : NYX = (YY/YX)*NXY ; 90 : * 91 : ******************************************************* 92 : * Ricavo i valori D1 e D2 per ogni lamina 93 : ******************************************************* 94 : * 95 : D1 = YX/(1-(NXY*NYX)) ; 96 : D2 = YY/(1-(NXY*NYX)) ; 97 : POLD=PD1; 98 : PD1 = INSERE POLD N2 D1 ; DETR POLD ; 99 : POLD=PD2; 100 : PD2 = INSERE POLD N2 D2 ; DETR POLD ; 101 : * 102 : * Ricavo num e denom per il calcolo di z neutra 103 : * 104 : NZN1 = (D1*((ZS**2)-(ZI**2)))+NZN1 ; 105 : DZN1 = (D1*EP1)+DZN1 ; 106 : NZN2 = (D2*((ZS**2)-(ZI**2)))+NZN2 ; 107 : DZN2 = (D2*EP1)+DZN2 ; 108 : * 109 : N2 = N2 + 1; 110 : FIN LOOP2 ; 111 : * 112 : ****************************************************** 113 : * Calcolo la quota delle fibre medie 114 : ****************************************************** 115 : * 116 : ZN1 = 0.5*NZN1/DZN1 ; 117 : ZN2 = 0.5*NZN2/DZN2 ; 118 : * 119 : ****************************************************** 120 : * Calcolo SUM1 e SUM2 per ogni lamina 121 : ****************************************************** 122 : * 123 : PSUM1 = PROG ; PSUM2 = PROG ; 124 : POLD = PSUM1; 125 : PSUM1 = INSERE POLD 1 0. ; DETR POLD ; 126 : POLD=PSUM2; 127 : PSUM2 = INSERE POLD 1 0. ; DETR POLD ; 128 : N2 = NLAM ; 129 : SUM1 = 0. ; SUM2 = 0. ; 130 : SI (NLAM > 1 ) ; 131 : REPETER LOOP3 (NLAM-1) ; 132 : ZS = EXTR PZS N2 ; 133 : ZI = EXTR PZI N2 ; 134 : D1 = EXTR PD1 N2 ; 135 : D2 = EXTR PD2 N2 ; 136 : * 137 : SUM1 = (D1*((((ZS**2)-(ZI**2))/2)-(EP1*ZN1))) 138 : + SUM1 ; 139 : SUM2 = (D2*((((ZS**2)-(ZI**2))/2)-(EP1*ZN2))) 140 : + SUM2 ; 141 : * 142 : N2 = N2-1 ; 143 : POLD=PSUM1; 144 : PSUM1 = INSERE POLD (NLAM+1-N2) SUM1 ; DETR POLD ; 145 : POLD=PSUM2; 146 : PSUM2 = INSERE POLD (NLAM+1-N2) SUM2 ; DETR POLD ; 147 : FIN LOOP3 ; 148 : FINSI; 149 : N2 = 1 ; 150 : R1 = 0 ; R2 = 0 ; 151 : I1TOT = 0 ; I2TOT = 0 ; 152 : PI1 = PROG ; PI2 = PROG ; 153 : REPETER LOOP4 NLAM ; 154 : D1 = EXTR PD1 N2 ; 155 : D2 = EXTR PD2 N2 ; 156 : ZS = EXTR PZS N2 ; 157 : ZI = EXTR PZI N2 ; 158 : SUM1 = EXTR PSUM1 (NLAM+1-N2) ; 159 : SUM2 = EXTR PSUM2 (NLAM+1-N2) ; 160 : * 161 : ******************************************************** 162 : * Calcolo R1 e R2 163 : ******************************************************** 164 : * 165 : R1 = (D1*((((ZS**3)-(ZI**3))/3)-(((ZS**2)-(ZI**2))*ZN1)+ 166 : (EP1*(ZN1**2))))+R1 ; 167 : 168 : R2 = (D2*((((ZS**3)-(ZI**3))/3)-(((ZS**2)-(ZI**2))*ZN2)+ 169 : (EP1*(ZN2**2))))+R2 ; 170 : * 171 : * 172 : ******************************************************** 173 : * Calcolo Gxz eGyz rispetto al riferimento 174 : * globale 175 : ******************************************************** 176 : * 177 : EP1 = TCOMP.N1.N2.EPA ; 178 : AN1 = TCOMP.N1.N2.ANG ; 179 : ID1 = TCOMP.N1.N2.IDM ; 180 : * 181 : TCM = TCOMP.MCOL ; 182 : G2 = TCM.ID1.G23 ; G3 = TCM.ID1.G13 ; 183 : * 184 : BETA = AN1 ; 185 : * 186 : CB = COS BETA; SB = SIN BETA; 187 : GXZ = (G2*G3)/(((SB**2)*G3)+((CB**2)*G2)) ; 188 : GYZ = (G2*G3)/(((CB**2)*G3)+((SB**2)*G2)) ; 189 : * 190 : ******************************************************** 191 : * Calcolo I1 per ogni lamina 192 : ******************************************************** 193 : * 194 : A11 = (-1)*(D1/2) ; 195 : A21 = D1*ZN1 ; 196 : A31 = ((-1)*A11*(ZI**2))-(D1*ZI*ZN1)-SUM1 ; 197 : I1 = (((((ZS**5)-(ZI**5))*(A11**2))/5)+((((ZS**4)-(ZI**4)) 198 : *A11*A21)/2)+((((A21**2)+(2*A11*A31))/3)*((ZS**3)-(ZI**3))) 199 : +(A21*A31*((ZS**2)-(ZI**2)))+(EP1*(A31**2))) 200 : /GXZ ; 201 : * 202 : POLD=PI1; 203 : PI1 = INSERE POLD N2 I1 ; DETR POLD ; 204 : * 205 : ******************************************************** 206 : * Calcolo I2 per ogni lamina 207 : ******************************************************** 208 : * 209 : A12 = (-1)*(D2/2) ; 210 : A22 = D2*ZN2 ; 211 : A32 = ((-1)*A12*(ZI**2))-(D2*ZI*ZN2)-SUM2 ; 212 : I2 = (((((ZS**5)-(ZI**5))*(A12**2))/5)+((((ZS**4)-(ZI**4)) 213 : *A12*A22)/2)+((((A22**2)+(2*A12*A32))/3)*((ZS**3)-(ZI**3))) 214 : +(A22*A32*((ZS**2)-(ZI**2)))+(EP1*(A32**2))) 215 : /GYZ ; 216 : * 217 : POLD=PI2; 218 : PI2 = INSERE POLD N2 I2 ; DETR POLD ; 219 : * 220 : * Calcolo I1TOT e I2TOT 221 : * 222 : I1TOT = I1 + I1TOT ; 223 : I2TOT = I2 + I2TOT ; 224 : * 225 : N2 = N2+1 ; 226 : FIN LOOP4 ; 227 : * 228 : N2 = 1 ; 229 : REPETER LOOP5 NLAM; 230 : I1 = EXTR PI1 N2 ; 231 : I2 = EXTR PI2 N2 ; 232 : * 233 : EP1 = TCOMP.N1.N2.EPA ; 234 : AN1 = TCOMP.N1.N2.ANG ; 235 : * 236 : BETA = AN1 ; 237 : CB = COS BETA; SB = SIN BETA; 238 : * 239 : ******************************************************** 240 : * Calcolo i valori corretti per GXY e GYZ 241 : ******************************************************** 242 : * 243 : GXZ = ((R1**2)*I1)/(EP1*(I1TOT**2)) ; 244 : GYZ = ((R2**2)*I2)/(EP1*(I2TOT**2)) ; 245 : * 246 : * 247 : * e li ruoto per ciascuna lamina nel sistema 248 : * di riferimento ortotropo 249 : * 250 : SI (ABS((CB**4)-(SB**4)) <EG 0.01) ; 251 : G1 = (GXZ + GYZ)*3/5 ; 252 : G2 = G1 ; 253 : SINON ; 254 : G1 = ((((CB**2)*GXZ)-((SB**2)*GYZ))/((CB**4)-(SB**4)))*6/5 ; 255 : G2 = ((((SB**2)*GXZ)-((CB**2)*GYZ))/((SB**4)-(CB**4)))*6/5 ; 256 : FINSI ; 257 : * 258 : TG.N1.N2.G13 = G1 ; 259 : TG.N1.N2.G23 = G2 ; 260 : N2 = N2 + 1 ; 261 : FIN LOOP5 ; 262 : * 263 : N1 = N1 + 1 ; 264 : FIN LOOP1 ; 265 : * 266 : ************************************************************************ 267 : * FINE PROCEDURA @LAKAPPA 268 : ************************************************************************ 269 : * 270 : FINPROC ; 271 : 272 :
© Cast3M 2003 - All rights reserved.
Disclaimer