Changeset 8303 for issm/trunk/src/c/objects/Elements/TriaRef.cpp
- Timestamp:
- 05/16/11 15:42:41 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Elements/TriaRef.cpp
r6987 r8303 59 59 void TriaRef::GetBMacAyeal(double* B, double* xyz_list, GaussTria* gauss){ 60 60 /*Compute B matrix. B=[B1 B2 B3] where Bi is of size 3*NDOF2. 61 * For gridi, Bi can be expressed in the actual coordinate system61 * For node i, Bi can be expressed in the actual coordinate system 62 62 * by: 63 63 * Bi=[ dh/dx 0 ] 64 64 * [ 0 dh/dy ] 65 65 * [ 1/2*dh/dy 1/2*dh/dx ] 66 * where h is the interpolation function for gridi.66 * where h is the interpolation function for node i. 67 67 * 68 68 * We assume B has been allocated already, of size: 3x(NDOF2*NUMNODES) … … 90 90 91 91 /*Compute B matrix. B=[B1 B2 B3] where Bi is of size 3*NDOF2. 92 * For gridi, Bi can be expressed in the actual coordinate system92 * For node i, Bi can be expressed in the actual coordinate system 93 93 * by: 94 94 * Bi=[ dh/dx 0 ] 95 95 * [ 0 dh/dy ] 96 96 * [ 1/2*dh/dy 1/2*dh/dx ] 97 * where h is the interpolation function for gridi.97 * where h is the interpolation function for node i. 98 98 * 99 99 * We assume B has been allocated already, of size: 3x(NDOF2*NUMNODES) … … 158 158 void TriaRef::GetBPrognostic(double* B_prog, double* xyz_list, GaussTria* gauss){ 159 159 /*Compute B matrix. B=[B1 B2 B3] where Bi is of size 3*NDOF2. 160 * For gridi, Bi can be expressed in the actual coordinate system160 * For node i, Bi can be expressed in the actual coordinate system 161 161 * by: 162 162 * Bi=[ h ] 163 163 * [ h ] 164 * where h is the interpolation function for gridi.164 * where h is the interpolation function for node i. 165 165 * 166 166 * We assume B_prog has been allocated already, of size: 2x(NDOF1*NUMNODES) … … 183 183 184 184 /*Compute B' matrix. B'=[B1' B2' B3'] where Bi' is of size 3*NDOF2. 185 * For gridi, Bi' can be expressed in the actual coordinate system185 * For node i, Bi' can be expressed in the actual coordinate system 186 186 * by: 187 187 * Bi_prime=[ 2*dh/dx dh/dy ] 188 188 * [ dh/dx 2*dh/dy ] 189 189 * [ dh/dy dh/dx ] 190 * where h is the interpolation function for gridi.190 * where h is the interpolation function for node i. 191 191 * 192 192 * We assume B' has been allocated already, of size: 3x(NDOF2*NUMNODES) … … 214 214 215 215 /*Compute Bprime matrix. Bprime=[Bprime1 Bprime2 Bprime3] where Bprimei is of size 3*NDOF2. 216 * For gridi, Bprimei can be expressed in the actual coordinate system216 * For node i, Bprimei can be expressed in the actual coordinate system 217 217 * by: 218 218 * Bprimei=[ dh/dx 0 ] … … 220 220 * [ dh/dy dh/dx ] 221 221 * [ dh/dx dh/dy ] 222 * where h is the interpolation function for gridi.222 * where h is the interpolation function for node i. 223 223 * 224 224 * We assume Bprime has been allocated already, of size: 3x(NDOF2*NUMNODES) … … 247 247 void TriaRef::GetBprimePrognostic(double* Bprime_prog, double* xyz_list, GaussTria* gauss){ 248 248 /*Compute B' matrix. B'=[B1' B2' B3'] where Bi' is of size 3*NDOF2. 249 * For gridi, Bi' can be expressed in the actual coordinate system249 * For node i, Bi' can be expressed in the actual coordinate system 250 250 * by: 251 251 * Bi_prime=[ dh/dx ] 252 252 * [ dh/dy ] 253 * where h is the interpolation function for gridi.253 * where h is the interpolation function for node i. 254 254 * 255 255 * We assume B' has been allocated already, of size: 3x(NDOF2*NUMNODES) … … 272 272 void TriaRef::GetL(double* L, double* xyz_list,GaussTria* gauss,int numdof){ 273 273 /*Compute L matrix. L=[L1 L2 L3] where Li is square and of size numdof. 274 * For gridi, Li can be expressed in the actual coordinate system274 * For node i, Li can be expressed in the actual coordinate system 275 275 * by: 276 276 * numdof=1: … … 279 279 * Li=[ h 0 ] 280 280 * [ 0 h ] 281 * where h is the interpolation function for gridi.281 * where h is the interpolation function for node i. 282 282 * 283 283 * We assume L has been allocated already, of size: NUMNODES (numdof=1), or numdofx(numdof*NUMNODES) (numdof=2)
Note:
See TracChangeset
for help on using the changeset viewer.