question about 3D-shell element

[zhuying0727]

Dear all:

I have three questions about 3D-shell element in file 'shell3d.f'.

1. In subroutine 'jacq3d', the routine compute geometric constants for DKQ element. I want to know what constants are computed?

2. I also want to know the meaning of some arrays b(4), c(4), aa(4), bb(4), cc(4), dd(4) and ee(4) in common block 'shle16.h' and other arrays bm(3,6,9), ii1 and ii2 in common block 'shld16.h'.

3. What is the difference between subroutine 'shl3di' and subroutine' shl3ds'? And what is the meaning of parameter of d(102)?

Thank you!

[Prof. R.L. Taylor]

shl3di integrates through the thickness to compute the force resultants.  It can handle more materials.  shl3ds uses elastic integrated resultants.

The DKQ parameters are used to compute the shape functions.  Coded too long ago to recall all the parameters meaning.  Look where used in routine to see how they affect solution.