00001 double stress[6]={0,0,0,0,0,0};
00002
00003
00004
00005
00006
00007
00008
00009 double dstran[6]={0,0,0,0,0,0};
00010
00011
00012
00013
00014
00015
00016
00017 int ndi=3;
00018 int nshr=3;
00019 int ntens=ndi+nshr;
00020
00021 int nprops=11;
00022 const int maxnprops=50;
00023 double props[maxnprops];
00024 memset(props, 0, sizeof(*props)*50);
00025
00026 for(int i=1; i<(nprops+1); i++) props[i-1]=0.2;
00027
00028 double dtime=Mp->timecon.actualforwtimeincr ();
00029
00030 int nstatv=6;
00031 const int maxnstatv=50;
00032 double statev[maxnstatv];
00033 memset(statev, 0, sizeof(*statev)*maxnstatv);
00034 for(int i=0; i<nstatv; i++) statev[i]=0.3;
00035
00036 double ddum=0;
00037 int idum=0;
00038 double ddum6[6]={0,0,0,0,0,0};
00039 double ddum3[3]={0,0,0};
00040 double ddum36[6*6];
00041 memset(ddum36, 0, sizeof(*ddum36)*36);
00042 char cmname []="tmp";
00043
00044 double dfgrd[3*3]={1,0,0,0,1,0,0,0,1};
00045
00046 static double time=0;
00047
00048 double unsatvar[4]={0,0,0,0};
00049 double stresssuction=0.0;
00050
00051 unsatvar[0]=-stresssuction;
00052 statev[1]=unsatvar[0];
00053
00054 umatunsat_(*stress, *statev, *ddum36,ddum, ddum, ddum, ddum, *ddum36,*ddum36, *ddum36, *ddum6, *dstran,
00055 time, dtime, ddum, ddum, *ddum6, *ddum6, *cmname, ndi,nshr, ntens, nstatv, *props,nprops, *ddum3, *dfgrd, ddum,
00056 ddum, *dfgrd, *dfgrd, idum, idum, idum, idum, idum, idum, *unsatvar);