Commit e0f0daf2 authored by François Trigaux's avatar François Trigaux
Browse files

Twisted Beam test case is added

parent 260f1400
[Beam]
analysis_flag = 0
nkp = 51
nkp = 31
niter = 1000
nstep = 1
ncond_pt = 2
......
1
0.0
3.999999999999999758e-11 0.000000000000000000e+00 0.000000000000000000e+00 0.000000000000000000e+00 0.000000000000000000e+00 0.000000000000000000e+00
0.000000000000000000e+00 1.210592686002522039e-10 0.000000000000000000e+00 0.000000000000000000e+00 0.000000000000000000e+00 0.000000000000000000e+00
0.000000000000000000e+00 0.000000000000000000e+00 1.210592686002522039e-10 0.000000000000000000e+00 0.000000000000000000e+00 0.000000000000000000e+00
......
2
00.0 00.0
10.0 90.0
\ No newline at end of file
......@@ -24,13 +24,14 @@ int main(int argc, char* argv[])
initBeam(bm);
printf("...Done!\n");
double *twist = malloc(sizeof(double)*bm->nelem);
for(i=0;i<bm->nmemb;i++)
/*double *twist = malloc(sizeof(double)*bm->nkp);
for(i=0;i<bm->nkp;i++)
{
twist[i] = 90.0 * M_PI/180.0 * i / ((double) (bm->nmemb-1.0));
twist[i] = 90.0 * M_PI/180.0 * i / ((double) (bm->nkp - 1.0));
}
beam_setTwist(bm,twist,M_PI/20.0);
free(twist);
beam_setTwist(bm,twist);
free(twist);*/
beam_setTwistFromFile(bm,"Twist.dat");
// Creating an array with the uniform load for each member
printf("Creating the loads...\n");
......
......@@ -38,7 +38,9 @@ k = 5.0/6.0
stiff = np.diag(np.array([E*A, k*G*A, k*G*A, G*J, E*Iy, E*Iz]));
#stiff = 1e8 * np.diag(np.array([250, 92.449, 83.497, 1.498, 5.208, 1.302]));
comp = np.linalg.inv(stiff);
np.savetxt("Compliance.dat",comp);
with open("Compliance.dat","w") as fid:
fid.write("1\n0.0\n");
np.savetxt(fid,comp);
sp.run(["./gebt", "-%e"%(F)]);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment