00001 #ifndef LINHEX_H
00002 #define LINHEX_H
00003
00004 #include "alias.h"
00005 struct matrix;
00006 struct vector;
00007 struct ivector;
00008
00009
00010
00011
00012
00013
00014
00015
00016 class linhex
00017 {
00018 public:
00019 linhex ();
00020 ~linhex ();
00021
00022 double approx (double xi,double eta,double zeta,vector &nodval);
00023 void bf_matrix (matrix &n,double xi,double eta,double zeta);
00024 void geom_matrix (matrix &gm,vector &x,vector &y,vector &z,
00025 double xi,double eta,double zeta,double &jac);
00026 void geom_matrix_block (matrix &gm,long ri,vector &x,vector &y,vector &z,
00027 double xi,double eta,double zeta,double &jac);
00028 void bvectors (vector &x,vector &y,vector &z,double xi,double eta,double zeta,double &jac,
00029 vector &b11,vector &b12,vector &b13,
00030 vector &b21,vector &b22,vector &b23,
00031 vector &b31,vector &b32,vector &b33);
00032 void gngeom_matrix (matrix &gm,vector &r,vector &x,vector &y,vector &z,double xi,double eta,double zeta,double &jac);
00033 void gnl_grmatrix (matrix &grm,vector &x,vector &y,vector &z,double xi,double eta,double zeta,double &jac);
00034
00035
00036 void transf_matrix (ivector &nodes,matrix &tmat);
00037
00038 void gl_stiffness_matrix (long eid,long ri,long ci,matrix &sm);
00039 void gnl_stiffness_matrix (long lcid,long eid,long ri,long ci,matrix &sm);
00040 void res_stiffness_matrix (long lcid,long eid,matrix &sm);
00041
00042 void res_mass_matrix (long eid,matrix &mm);
00043 void mass_matrix (long eid,matrix &mm);
00044 void load_matrix (long eid,matrix &lm);
00045 void res_load_matrix (long eid,matrix &lm);
00046
00047 void res_ip_strains (long lcid,long eid);
00048 void gl_ip_strains (long lcid,long eid,long ri,long ci,vector &x,vector &y,vector &z,vector &r);
00049 void gnl_ip_strains (long lcid,long eid,long ri,long ci,vector &x,vector &y,vector &z,vector &r);
00050 void nod_strains_ip (long lcid,long eid,long ri,long ci);
00051 void nod_strains_comp (long lcid,long eid,double **stra);
00052
00053
00054
00055
00056 void appstrain (long lcid,long eid,double xi,double eta,double zeta,long fi,long li,vector &eps);
00057 void strains (long lcid,long eid,long ri,long ci);
00058
00059 void nodecoord (vector &xi,vector &eta,vector &zeta);
00060 void nodipnum (long eid,long ri,long ci,ivector &ipnum);
00061
00062
00063 void res_ip_stresses (long lcid,long eid);
00064 void ip_stresses (long lcid,long eid,long ri,long ci);
00065 void ip_elast_stresses (long lcid,long eid,long ri,long ci);
00066 void nod_stresses_ip (long lcid,long eid,long ri,long ci);
00067 void elem_stresses (double **stra,double **stre,long lcid,long eid,long ri,long ci);
00068 void appstress (long lcid,long eid,double xi,double eta,double zeta,long fi,long li,vector &sig);
00069 void stresses (long lcid,long eid,long ri,long ci);
00070
00071 void nod_others (long lcid,long eid,long ri,long ci);
00072 void nod_eqother_ip (long eid,long ri,long ci);
00073
00074
00075 void gl_internal_forces (long lcid,long eid,long ri,long ci,vector &ifor,vector &x,vector &y,vector &z);
00076 void gnl_internal_forces (long lcid,long eid,long ri,long ci,vector &ifor);
00077 void nonloc_internal_forces (long lcid,long eid,long ri,long ci,vector &ifor,vector &x,vector &y,vector &z);
00078 void incr_internal_forces (long lcid,long eid,long ri,long ci,vector &ifor,vector &x,vector &y,vector &z);
00079 void eigstrain_forces (long lcid,long eid,long ri,long ci,vector &nfor,vector &x,vector &y,vector &z);
00080
00081 void res_internal_forces (long lcid,long eid,vector &ifor);
00082 void res_nonloc_internal_forces (long lcid,long eid,vector &ifor);
00083 void res_incr_internal_forces (long lcid,long eid,vector &ifor);
00084 void res_eigstrain_forces (long lcid,long eid,vector &nfor);
00085
00086 void compute_nlstress (long lcid,long eid,long ri,long ci);
00087 void compute_nlstressincr (long lcid,long eid,long ri,long ci);
00088 void local_values (long lcid,long eid,long ri,long ci);
00089 void compute_nonloc_nlstress (long lcid,long eid,long ri,long ci);
00090 void compute_eigstress (long lcid,long eid,long ri,long ci);
00091 void elem_integration (integratedquant iq,long lcid,long eid,long ri,long ci,vector &nv,vector &x,vector &y,vector &z);
00092
00093 void ipcoord (long eid,long ipp,long ri,long ci,vector &coord);
00094 void inicipval(long eid, long ri, long ci, matrix &nodval, inictype *ictn);
00095 void ipvolume (long eid,long ri,long ci);
00096
00097 void locglob_nodeval (long is,vector &nv,double *tnv,vector &x,vector &y,vector &z);
00098 void node_forces_surf (long lcid,long eid,long *is,double *nv,vector &nf);
00099 void node_forces_surf_old (long lcid,long eid,long *is,double *nv,vector &nf);
00100 void tran_mat(matrix &tran, vector &gx, vector &gy, vector &gz, long is);
00101
00102
00103 void intpointval (long eid,vector &nodval,vector &ipval);
00104
00105 void aver_strains (long lcid,long eid,long ri,long ci,vector &averstra,double &volume);
00106
00107
00108 void surfnodeval (long surf,vector &nodval,double *list);
00109
00110 void find_extreme_strains (vector &min,vector &max,long lcid,long eid);
00111 void find_extreme_stresses (vector &min,vector &max,long lcid,long eid);
00112
00113
00114
00115 long ndofe;
00116
00117 long nne;
00118
00119 long tncomp;
00120
00121 long tnip;
00122
00123 long napfun;
00124
00125 long ned;
00126
00127 long nned;
00128
00129 long nsurf;
00130
00131 long nnsurf;
00132
00133 long **intordsm;
00134
00135 long intordmm;
00136
00137 long intordb;
00138
00139 long **nip;
00140
00141 long nb;
00142
00143 long *ncomp;
00144
00145 long *cncomp;
00146
00147 strastrestate ssst;
00148
00149 };
00150
00151 #endif