loop for stresses and strains

sanehatter (47)
Hi There,

This is part of a program I am writing. This part of the code is responsible for computing commulative stresses and strains. The problem here is that the first part of the code works fine (as far as I can tell) but when it reaches the loop segment below it stops functioning and the compiler (code::blocks on Win XP) gives me this error:

Process terminated with status -1073741510 (0 minutes, 10 seconds)

I tried putting a "cout" statement inside this loop but it wouldn't execute. Nothing inside this loops seems to function. This is the problematic part of the code:



1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
 
int xx1 = 0, qq1 = 0, contactpnt1 = 0, l1, i1;

for (i1 = 0; i1 < isobox1; i1= i1+1){

xx1 = isocounterbox1[i1];

loop1:

for ( l1 = 0; l1 < isobox1; l1 = l1+1)

{
    if (l1 == i1) {goto loop1;}

qq1 = isocounterbox1[l1];

if ( (pow((rightisosceles [xx1][0][0] - rightisosceles [qq1][0][0]), 2)) + (pow ((rightisosceles [xx1][0][1] - rightisosceles [qq1][0][1]), 2)) <= 18)

{

contactpnt1 = contactpnt1 + 1;

if (contactpnt1 > 4)  {goto next1; };

ydisplat [xx1][contactpnt1] = ydisplacement2 + (deluy [xx1][contactpnt1] - deluy [qq1][contactpnt1] + deltheta [xx1][contactpnt1] * ( rightisosceles [qq1][1][0] - rightisosceles [xx1][0][0])  - deltheta [qq1][contactpnt1] * (rightisosceles [qq1][0][0] - rightisosceles [qq1][0][0]));
xdisplat [xx1][contactpnt1] = xdisplacement2 + (delux [xx1][contactpnt1] - delux [qq1][contactpnt1] + deltheta [xx1][contactpnt1] * ( rightisosceles [qq1][1][1] - rightisosceles [xx1][0][1])  - deltheta [qq1][contactpnt1] * (rightisosceles [qq1][1][1] -rightisosceles [qq1][0][1]));

delus [xx1][contactpnt1] = (xdisplat [xx1][contactpnt1] * cos (alpha [xx1]) + ydisplat [xx1][contactpnt1] * sin (alpha [xx1]));
delun [xx1][contactpnt1] = (ydisplat [xx1][contactpnt1] * cos (alpha [xx1]) - xdisplat [xx1][contactpnt1] * sin (alpha [xx1]));

fn [xx1][contactpnt1] =  delun [xx1][contactpnt1] * (-1) *  dstiffness_coefficient;
fs [xx1][contactpnt1] =  delus [xx1][contactpnt1] *  dstiffness_coefficient;

dn [xx1][contactpnt1] = delun [xx1][contactpnt1] * (-1) * ddamping_coefficient;
ds [xx1][contactpnt1] = delus [xx1][contactpnt1] * ddamping_coefficient;

yforce [qq1][contactpnt1] =  (((fs [xx1][contactpnt1] + ds [xx1][contactpnt1]) * sin (alpha [xx1])) - ((fn [xx1][contactpnt1] + dn [xx1][contactpnt1]) * cos (alpha [xx1])));
xforce [qq1][contactpnt1] =  (((fs [xx1][contactpnt1] + ds [xx1][contactpnt1]) * cos (alpha [xx1])) + ((fn [xx1][contactpnt1] + dn [xx1][contactpnt1]) * sin (alpha [xx1])));

yforce [xx1][contactpnt1] = (yforce [qq1][contactpnt1] * -1);
xforce [xx1][contactpnt1] = (xforce [qq1][contactpnt1] * -1);

fxsum [xx1][contactpnt1] = xforce [xx1][contactpnt1];
x [xx1] = fxsum [xx1][contactpnt1];
resultfx [xx1] += x[xx1];

fysum [xx1][contactpnt1] = yforce [xx1][contactpnt1]+ p2;
y [xx1] = fysum [xx1][contactpnt1];
resultfy [xx1] += y[xx1];

msum [xx1][contactpnt1] = (yforce [xx1][contactpnt1]* ( rightisosceles [qq1][1][0] - rightisosceles [xx1][0][0]) - xforce [xx1][contactpnt1] * ( rightisosceles [qq1][1][1] - rightisosceles [xx1][0][1]));
m [xx1] = msum [xx1][contactpnt1];
resultm [xx1] += m[xx1];

udoty [xx1][contactpnt1]=  ((fxsum [xx1][contactpnt1] * dtime_increment)/ dmass);
udotysum [xx1] = udoty [xx1][contactpnt1];
resultudoty [xx1] += udotysum [xx1];

udotx [xx1][contactpnt1] =  ((fxsum [xx1][contactpnt1]* dtime_increment)/ dmass);
udotxsum [xx1] = udotx [xx1][contactpnt1];
resultudotx [xx1] += udotxsum [xx1];

thetadot [xx1][contactpnt1] = ((msum [xx1][contactpnt1] * dtime_increment)/ ii);
thetadotsum [xx1] = thetadot [xx1][contactpnt1];
resultthetadot [xx1] += thetadotsum [xx1];

deluy [xx1][contactpnt1] = udoty [xx1][contactpnt1] * dtime_increment;
deluysum [xx1] = deluy [xx1][contactpnt1];
resultdeluy [xx1] += deluysum [xx1];

delux [xx1][contactpnt1] = udotx [xx1][contactpnt1] * dtime_increment;
deluxsum [xx1] = delux [xx1][contactpnt1];
resultdelux [xx1] += deluxsum [xx1];

deltheta [xx1][contactpnt1] = thetadot [xx1][contactpnt1] * dtime_increment;
delthetasum [xx1] = deltheta [xx1][contactpnt1];
resultdeltheta [xx1] += delthetasum [xx1];

uy [xx1][contactpnt1] = deluy [xx1][contactpnt1];
uysum [xx1] = uy [xx1][contactpnt1];
resultuy [xx1]+= uysum [xx1];

ux [xx1][contactpnt1] = delux [xx1][contactpnt1];
uxsum [xx1] = ux [xx1][contactpnt1];
resultux [xx1] += uxsum [xx1];

theta [xx1][contactpnt1] = deltheta [xx1][contactpnt1];
thetasum [xx1] = theta [xx1][contactpnt1];
resulttheta [xx1] += thetasum [xx1];
alpha [xx1] = alpha1 + resulttheta [xx1];

cout<< i1 << endl;
}}

next1:
;}



Any help is highly appreciated.




closed account (D80DSL3A)
I bet it never gets past line 12. You have an infinite loop.
On the 1st iteration line 3 il=0. 1st iteration in for loop on line 9 makes ll=0.
Line 12 tests if ll == il and it does since 0 == 0 so execution goes to line 7. ll and il remain 0 forever.
sanehatter (47)
I did change that. I put "loop1:" at line 2 right now but it still isn't giving me any results.
Moooce (216)
That'll still do the same thing though. It'll just run a bit more code before it gets to this line
if (l1 == i1) {goto loop1;}
which is causing the infinate loop as fun2code said.

If you use the debugger to step though the code and look at the vaues of il and ll you'll see what is happening. (but I don't know how to step through code using code blocks)
sanehatter (47)
Ok I removed the goto statement altogether and tried printing some values from the equations (results of resultuy [xx1] for example) but now the values I am getting are all gibrish (mixture of numbers and letters)...
Last edited on
sanehatter (47)
I even tried printing the results of other equations in this loop, as well, but I still don't have any legitimate (numerical) results.
sanehatter (47)
This is the updated part of the code:


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
 
int xx1 = 0, qq1 = 0, contactpnt1 = 0, l1, i1;



for (i1 = 0; i1 < isobox1; i1= i1+1)

{

xx1 = isocounterbox1[i1];



for ( l1 = 0; l1 < isobox1; l1 = l1+1)

{


qq1 = isocounterbox1[l1];

if ( (pow((rightisosceles [xx1][0][0] - rightisosceles [qq1][0][0]), 2)) + (pow ((rightisosceles [xx1][0][1] - rightisosceles [qq1][0][1]), 2)) <= 18)

{

contactpnt1 = contactpnt1 + 1;

if (contactpnt1 > 4)  {goto next1; };

ydisplat [xx1][contactpnt1] = ydisplacement2 + (deluy [xx1][contactpnt1] - deluy [qq1][contactpnt1] + deltheta [xx1][contactpnt1] * ( rightisosceles [qq1][1][0] - rightisosceles [xx1][0][0])  - deltheta [qq1][contactpnt1] * (rightisosceles [qq1][0][0] - rightisosceles [qq1][0][0]));
xdisplat [xx1][contactpnt1] = xdisplacement2 + (delux [xx1][contactpnt1] - delux [qq1][contactpnt1] + deltheta [xx1][contactpnt1] * ( rightisosceles [qq1][1][1] - rightisosceles [xx1][0][1])  - deltheta [qq1][contactpnt1] * (rightisosceles [qq1][1][1] -rightisosceles [qq1][0][1]));

delus [xx1][contactpnt1] = (xdisplat [xx1][contactpnt1] * cos (alpha [xx1]) + ydisplat [xx1][contactpnt1] * sin (alpha [xx1]));
delun [xx1][contactpnt1] = (ydisplat [xx1][contactpnt1] * cos (alpha [xx1]) - xdisplat [xx1][contactpnt1] * sin (alpha [xx1]));

fn [xx1][contactpnt1] =  delun [xx1][contactpnt1] * (-1) *  dstiffness_coefficient;
fs [xx1][contactpnt1] =  delus [xx1][contactpnt1] *  dstiffness_coefficient;

dn [xx1][contactpnt1] = delun [xx1][contactpnt1] * (-1) * ddamping_coefficient;
ds [xx1][contactpnt1] = delus [xx1][contactpnt1] * ddamping_coefficient;

yforce [qq1][contactpnt1] =  (((fs [xx1][contactpnt1] + ds [xx1][contactpnt1]) * sin (alpha [xx1])) - ((fn [xx1][contactpnt1] + dn [xx1][contactpnt1]) * cos (alpha [xx1])));
xforce [qq1][contactpnt1] =  (((fs [xx1][contactpnt1] + ds [xx1][contactpnt1]) * cos (alpha [xx1])) + ((fn [xx1][contactpnt1] + dn [xx1][contactpnt1]) * sin (alpha [xx1])));

yforce [xx1][contactpnt1] = (yforce [qq1][contactpnt1] * -1);
xforce [xx1][contactpnt1] = (xforce [qq1][contactpnt1] * -1);

fxsum [xx1][contactpnt1] = xforce [xx1][contactpnt1];
x [xx1] = fxsum [xx1][contactpnt1];
resultfx [xx1] += x[xx1];

fysum [xx1][contactpnt1] = yforce [xx1][contactpnt1]+ p2;
y [xx1] = fysum [xx1][contactpnt1];
resultfy [xx1] += y[xx1];

msum [xx1][contactpnt1] = (yforce [xx1][contactpnt1]* ( rightisosceles [qq1][1][0] - rightisosceles [xx1][0][0]) - xforce [xx1][contactpnt1] * ( rightisosceles [qq1][1][1] - rightisosceles [xx1][0][1]));
m [xx1] = msum [xx1][contactpnt1];
resultm [xx1] += m[xx1];

udoty [xx1][contactpnt1]=  ((fxsum [xx1][contactpnt1] * dtime_increment)/ dmass);
udotysum [xx1] = udoty [xx1][contactpnt1];
resultudoty [xx1] += udotysum [xx1];

udotx [xx1][contactpnt1] =  ((fxsum [xx1][contactpnt1]* dtime_increment)/ dmass);
udotxsum [xx1] = udotx [xx1][contactpnt1];
resultudotx [xx1] += udotxsum [xx1];

thetadot [xx1][contactpnt1] = ((msum [xx1][contactpnt1] * dtime_increment)/ ii);
thetadotsum [xx1] = thetadot [xx1][contactpnt1];
resultthetadot [xx1] += thetadotsum [xx1];

deluy [xx1][contactpnt1] = udoty [xx1][contactpnt1] * dtime_increment;
deluysum [xx1] = deluy [xx1][contactpnt1];
resultdeluy [xx1] += deluysum [xx1];

delux [xx1][contactpnt1] = udotx [xx1][contactpnt1] * dtime_increment;
deluxsum [xx1] = delux [xx1][contactpnt1];
resultdelux [xx1] += deluxsum [xx1];

deltheta [xx1][contactpnt1] = thetadot [xx1][contactpnt1] * dtime_increment;
delthetasum [xx1] = deltheta [xx1][contactpnt1];
resultdeltheta [xx1] += delthetasum [xx1];

uy [xx1][contactpnt1] = deluy [xx1][contactpnt1];
uysum [xx1] = uy [xx1][contactpnt1];
resultuy [xx1]+= uysum [xx1];

ux [xx1][contactpnt1] = delux [xx1][contactpnt1];
uxsum [xx1] = ux [xx1][contactpnt1];
resultux [xx1] += uxsum [xx1];

theta [xx1][contactpnt1] = deltheta [xx1][contactpnt1];
thetasum [xx1] = theta [xx1][contactpnt1];
resulttheta [xx1] += thetasum [xx1];
alpha [xx1] = alpha1 + resulttheta [xx1];

}}

next1:
;}


Can anyone tell me what's wrong with it and why am I not getting any results but an error message instead. This is the error message: Process terminated with status -1073741819 (0 minutes, 11 seconds).


Moooce (216)
Do you need to add code to reset these variables when you go to next1 ?
xx1 = 0, qq1 = 0, contactpnt1 = 0

Again, I'd say use the debugger to watch the values of the variables in your code when it's running, sorry I don't know what key to press for code:blocks but look for an option called "step" or "Set breakpoint"

Alternativley, at around line 94, write some code that prints out the values in xx1,contactpnt1,l1,i1,qq1 and then waits for a keypress. You can then look at the values and see if they are what you expect as the program runs.
Last edited on
Topic archived. No new replies allowed.