# top1-30x10.dat QQR2-MN-v-v # Original AMPL coding by Sven Leyffer, University of Dundee, April 2000. # # Topology optimization for compliance minimization of statically loaded # structures. # # Reference: O. Sigmund, "A 99 line topology optimization code written # in Matlab", Techincal University of Denmark, October 1999. # # Example 1 (30x10 grid) # ... number of elements in x-y direction param nelx = 30; param nely = 10; # ... other parameters param E := 1.0; param nu := 0.3; param p := 3.0; param volfrac := 0.5; # ... elements are defined by 8 nodes (square, two per vertex) set elements := ( 1 , 2 , 23 , 24 , 25 , 26 , 3 , 4 ) ( 23 , 24 , 45 , 46 , 47 , 48 , 25 , 26 ) ( 45 , 46 , 67 , 68 , 69 , 70 , 47 , 48 ) ( 67 , 68 , 89 , 90 , 91 , 92 , 69 , 70 ) ( 89 , 90 , 111 , 112 , 113 , 114 , 91 , 92 ) ( 111 , 112 , 133 , 134 , 135 , 136 , 113 , 114 ) ( 133 , 134 , 155 , 156 , 157 , 158 , 135 , 136 ) ( 155 , 156 , 177 , 178 , 179 , 180 , 157 , 158 ) ( 177 , 178 , 199 , 200 , 201 , 202 , 179 , 180 ) ( 199 , 200 , 221 , 222 , 223 , 224 , 201 , 202 ) ( 221 , 222 , 243 , 244 , 245 , 246 , 223 , 224 ) ( 243 , 244 , 265 , 266 , 267 , 268 , 245 , 246 ) ( 265 , 266 , 287 , 288 , 289 , 290 , 267 , 268 ) ( 287 , 288 , 309 , 310 , 311 , 312 , 289 , 290 ) ( 309 , 310 , 331 , 332 , 333 , 334 , 311 , 312 ) ( 331 , 332 , 353 , 354 , 355 , 356 , 333 , 334 ) ( 353 , 354 , 375 , 376 , 377 , 378 , 355 , 356 ) ( 375 , 376 , 397 , 398 , 399 , 400 , 377 , 378 ) ( 397 , 398 , 419 , 420 , 421 , 422 , 399 , 400 ) ( 419 , 420 , 441 , 442 , 443 , 444 , 421 , 422 ) ( 441 , 442 , 463 , 464 , 465 , 466 , 443 , 444 ) ( 463 , 464 , 485 , 486 , 487 , 488 , 465 , 466 ) ( 485 , 486 , 507 , 508 , 509 , 510 , 487 , 488 ) ( 507 , 508 , 529 , 530 , 531 , 532 , 509 , 510 ) ( 529 , 530 , 551 , 552 , 553 , 554 , 531 , 532 ) ( 551 , 552 , 573 , 574 , 575 , 576 , 553 , 554 ) ( 573 , 574 , 595 , 596 , 597 , 598 , 575 , 576 ) ( 595 , 596 , 617 , 618 , 619 , 620 , 597 , 598 ) ( 617 , 618 , 639 , 640 , 641 , 642 , 619 , 620 ) ( 639 , 640 , 661 , 662 , 663 , 664 , 641 , 642 ) ( 3 , 4 , 25 , 26 , 27 , 28 , 5 , 6 ) ( 25 , 26 , 47 , 48 , 49 , 50 , 27 , 28 ) ( 47 , 48 , 69 , 70 , 71 , 72 , 49 , 50 ) ( 69 , 70 , 91 , 92 , 93 , 94 , 71 , 72 ) ( 91 , 92 , 113 , 114 , 115 , 116 , 93 , 94 ) ( 113 , 114 , 135 , 136 , 137 , 138 , 115 , 116 ) ( 135 , 136 , 157 , 158 , 159 , 160 , 137 , 138 ) ( 157 , 158 , 179 , 180 , 181 , 182 , 159 , 160 ) ( 179 , 180 , 201 , 202 , 203 , 204 , 181 , 182 ) ( 201 , 202 , 223 , 224 , 225 , 226 , 203 , 204 ) ( 223 , 224 , 245 , 246 , 247 , 248 , 225 , 226 ) ( 245 , 246 , 267 , 268 , 269 , 270 , 247 , 248 ) ( 267 , 268 , 289 , 290 , 291 , 292 , 269 , 270 ) ( 289 , 290 , 311 , 312 , 313 , 314 , 291 , 292 ) ( 311 , 312 , 333 , 334 , 335 , 336 , 313 , 314 ) ( 333 , 334 , 355 , 356 , 357 , 358 , 335 , 336 ) ( 355 , 356 , 377 , 378 , 379 , 380 , 357 , 358 ) ( 377 , 378 , 399 , 400 , 401 , 402 , 379 , 380 ) ( 399 , 400 , 421 , 422 , 423 , 424 , 401 , 402 ) ( 421 , 422 , 443 , 444 , 445 , 446 , 423 , 424 ) ( 443 , 444 , 465 , 466 , 467 , 468 , 445 , 446 ) ( 465 , 466 , 487 , 488 , 489 , 490 , 467 , 468 ) ( 487 , 488 , 509 , 510 , 511 , 512 , 489 , 490 ) ( 509 , 510 , 531 , 532 , 533 , 534 , 511 , 512 ) ( 531 , 532 , 553 , 554 , 555 , 556 , 533 , 534 ) ( 553 , 554 , 575 , 576 , 577 , 578 , 555 , 556 ) ( 575 , 576 , 597 , 598 , 599 , 600 , 577 , 578 ) ( 597 , 598 , 619 , 620 , 621 , 622 , 599 , 600 ) ( 619 , 620 , 641 , 642 , 643 , 644 , 621 , 622 ) ( 641 , 642 , 663 , 664 , 665 , 666 , 643 , 644 ) ( 5 , 6 , 27 , 28 , 29 , 30 , 7 , 8 ) ( 27 , 28 , 49 , 50 , 51 , 52 , 29 , 30 ) ( 49 , 50 , 71 , 72 , 73 , 74 , 51 , 52 ) ( 71 , 72 , 93 , 94 , 95 , 96 , 73 , 74 ) ( 93 , 94 , 115 , 116 , 117 , 118 , 95 , 96 ) ( 115 , 116 , 137 , 138 , 139 , 140 , 117 , 118 ) ( 137 , 138 , 159 , 160 , 161 , 162 , 139 , 140 ) ( 159 , 160 , 181 , 182 , 183 , 184 , 161 , 162 ) ( 181 , 182 , 203 , 204 , 205 , 206 , 183 , 184 ) ( 203 , 204 , 225 , 226 , 227 , 228 , 205 , 206 ) ( 225 , 226 , 247 , 248 , 249 , 250 , 227 , 228 ) ( 247 , 248 , 269 , 270 , 271 , 272 , 249 , 250 ) ( 269 , 270 , 291 , 292 , 293 , 294 , 271 , 272 ) ( 291 , 292 , 313 , 314 , 315 , 316 , 293 , 294 ) ( 313 , 314 , 335 , 336 , 337 , 338 , 315 , 316 ) ( 335 , 336 , 357 , 358 , 359 , 360 , 337 , 338 ) ( 357 , 358 , 379 , 380 , 381 , 382 , 359 , 360 ) ( 379 , 380 , 401 , 402 , 403 , 404 , 381 , 382 ) ( 401 , 402 , 423 , 424 , 425 , 426 , 403 , 404 ) ( 423 , 424 , 445 , 446 , 447 , 448 , 425 , 426 ) ( 445 , 446 , 467 , 468 , 469 , 470 , 447 , 448 ) ( 467 , 468 , 489 , 490 , 491 , 492 , 469 , 470 ) ( 489 , 490 , 511 , 512 , 513 , 514 , 491 , 492 ) ( 511 , 512 , 533 , 534 , 535 , 536 , 513 , 514 ) ( 533 , 534 , 555 , 556 , 557 , 558 , 535 , 536 ) ( 555 , 556 , 577 , 578 , 579 , 580 , 557 , 558 ) ( 577 , 578 , 599 , 600 , 601 , 602 , 579 , 580 ) ( 599 , 600 , 621 , 622 , 623 , 624 , 601 , 602 ) ( 621 , 622 , 643 , 644 , 645 , 646 , 623 , 624 ) ( 643 , 644 , 665 , 666 , 667 , 668 , 645 , 646 ) ( 7 , 8 , 29 , 30 , 31 , 32 , 9 , 10 ) ( 29 , 30 , 51 , 52 , 53 , 54 , 31 , 32 ) ( 51 , 52 , 73 , 74 , 75 , 76 , 53 , 54 ) ( 73 , 74 , 95 , 96 , 97 , 98 , 75 , 76 ) ( 95 , 96 , 117 , 118 , 119 , 120 , 97 , 98 ) ( 117 , 118 , 139 , 140 , 141 , 142 , 119 , 120 ) ( 139 , 140 , 161 , 162 , 163 , 164 , 141 , 142 ) ( 161 , 162 , 183 , 184 , 185 , 186 , 163 , 164 ) ( 183 , 184 , 205 , 206 , 207 , 208 , 185 , 186 ) ( 205 , 206 , 227 , 228 , 229 , 230 , 207 , 208 ) ( 227 , 228 , 249 , 250 , 251 , 252 , 229 , 230 ) ( 249 , 250 , 271 , 272 , 273 , 274 , 251 , 252 ) ( 271 , 272 , 293 , 294 , 295 , 296 , 273 , 274 ) ( 293 , 294 , 315 , 316 , 317 , 318 , 295 , 296 ) ( 315 , 316 , 337 , 338 , 339 , 340 , 317 , 318 ) ( 337 , 338 , 359 , 360 , 361 , 362 , 339 , 340 ) ( 359 , 360 , 381 , 382 , 383 , 384 , 361 , 362 ) ( 381 , 382 , 403 , 404 , 405 , 406 , 383 , 384 ) ( 403 , 404 , 425 , 426 , 427 , 428 , 405 , 406 ) ( 425 , 426 , 447 , 448 , 449 , 450 , 427 , 428 ) ( 447 , 448 , 469 , 470 , 471 , 472 , 449 , 450 ) ( 469 , 470 , 491 , 492 , 493 , 494 , 471 , 472 ) ( 491 , 492 , 513 , 514 , 515 , 516 , 493 , 494 ) ( 513 , 514 , 535 , 536 , 537 , 538 , 515 , 516 ) ( 535 , 536 , 557 , 558 , 559 , 560 , 537 , 538 ) ( 557 , 558 , 579 , 580 , 581 , 582 , 559 , 560 ) ( 579 , 580 , 601 , 602 , 603 , 604 , 581 , 582 ) ( 601 , 602 , 623 , 624 , 625 , 626 , 603 , 604 ) ( 623 , 624 , 645 , 646 , 647 , 648 , 625 , 626 ) ( 645 , 646 , 667 , 668 , 669 , 670 , 647 , 648 ) ( 9 , 10 , 31 , 32 , 33 , 34 , 11 , 12 ) ( 31 , 32 , 53 , 54 , 55 , 56 , 33 , 34 ) ( 53 , 54 , 75 , 76 , 77 , 78 , 55 , 56 ) ( 75 , 76 , 97 , 98 , 99 , 100 , 77 , 78 ) ( 97 , 98 , 119 , 120 , 121 , 122 , 99 , 100 ) ( 119 , 120 , 141 , 142 , 143 , 144 , 121 , 122 ) ( 141 , 142 , 163 , 164 , 165 , 166 , 143 , 144 ) ( 163 , 164 , 185 , 186 , 187 , 188 , 165 , 166 ) ( 185 , 186 , 207 , 208 , 209 , 210 , 187 , 188 ) ( 207 , 208 , 229 , 230 , 231 , 232 , 209 , 210 ) ( 229 , 230 , 251 , 252 , 253 , 254 , 231 , 232 ) ( 251 , 252 , 273 , 274 , 275 , 276 , 253 , 254 ) ( 273 , 274 , 295 , 296 , 297 , 298 , 275 , 276 ) ( 295 , 296 , 317 , 318 , 319 , 320 , 297 , 298 ) ( 317 , 318 , 339 , 340 , 341 , 342 , 319 , 320 ) ( 339 , 340 , 361 , 362 , 363 , 364 , 341 , 342 ) ( 361 , 362 , 383 , 384 , 385 , 386 , 363 , 364 ) ( 383 , 384 , 405 , 406 , 407 , 408 , 385 , 386 ) ( 405 , 406 , 427 , 428 , 429 , 430 , 407 , 408 ) ( 427 , 428 , 449 , 450 , 451 , 452 , 429 , 430 ) ( 449 , 450 , 471 , 472 , 473 , 474 , 451 , 452 ) ( 471 , 472 , 493 , 494 , 495 , 496 , 473 , 474 ) ( 493 , 494 , 515 , 516 , 517 , 518 , 495 , 496 ) ( 515 , 516 , 537 , 538 , 539 , 540 , 517 , 518 ) ( 537 , 538 , 559 , 560 , 561 , 562 , 539 , 540 ) ( 559 , 560 , 581 , 582 , 583 , 584 , 561 , 562 ) ( 581 , 582 , 603 , 604 , 605 , 606 , 583 , 584 ) ( 603 , 604 , 625 , 626 , 627 , 628 , 605 , 606 ) ( 625 , 626 , 647 , 648 , 649 , 650 , 627 , 628 ) ( 647 , 648 , 669 , 670 , 671 , 672 , 649 , 650 ) ( 11 , 12 , 33 , 34 , 35 , 36 , 13 , 14 ) ( 33 , 34 , 55 , 56 , 57 , 58 , 35 , 36 ) ( 55 , 56 , 77 , 78 , 79 , 80 , 57 , 58 ) ( 77 , 78 , 99 , 100 , 101 , 102 , 79 , 80 ) ( 99 , 100 , 121 , 122 , 123 , 124 , 101 , 102 ) ( 121 , 122 , 143 , 144 , 145 , 146 , 123 , 124 ) ( 143 , 144 , 165 , 166 , 167 , 168 , 145 , 146 ) ( 165 , 166 , 187 , 188 , 189 , 190 , 167 , 168 ) ( 187 , 188 , 209 , 210 , 211 , 212 , 189 , 190 ) ( 209 , 210 , 231 , 232 , 233 , 234 , 211 , 212 ) ( 231 , 232 , 253 , 254 , 255 , 256 , 233 , 234 ) ( 253 , 254 , 275 , 276 , 277 , 278 , 255 , 256 ) ( 275 , 276 , 297 , 298 , 299 , 300 , 277 , 278 ) ( 297 , 298 , 319 , 320 , 321 , 322 , 299 , 300 ) ( 319 , 320 , 341 , 342 , 343 , 344 , 321 , 322 ) ( 341 , 342 , 363 , 364 , 365 , 366 , 343 , 344 ) ( 363 , 364 , 385 , 386 , 387 , 388 , 365 , 366 ) ( 385 , 386 , 407 , 408 , 409 , 410 , 387 , 388 ) ( 407 , 408 , 429 , 430 , 431 , 432 , 409 , 410 ) ( 429 , 430 , 451 , 452 , 453 , 454 , 431 , 432 ) ( 451 , 452 , 473 , 474 , 475 , 476 , 453 , 454 ) ( 473 , 474 , 495 , 496 , 497 , 498 , 475 , 476 ) ( 495 , 496 , 517 , 518 , 519 , 520 , 497 , 498 ) ( 517 , 518 , 539 , 540 , 541 , 542 , 519 , 520 ) ( 539 , 540 , 561 , 562 , 563 , 564 , 541 , 542 ) ( 561 , 562 , 583 , 584 , 585 , 586 , 563 , 564 ) ( 583 , 584 , 605 , 606 , 607 , 608 , 585 , 586 ) ( 605 , 606 , 627 , 628 , 629 , 630 , 607 , 608 ) ( 627 , 628 , 649 , 650 , 651 , 652 , 629 , 630 ) ( 649 , 650 , 671 , 672 , 673 , 674 , 651 , 652 ) ( 13 , 14 , 35 , 36 , 37 , 38 , 15 , 16 ) ( 35 , 36 , 57 , 58 , 59 , 60 , 37 , 38 ) ( 57 , 58 , 79 , 80 , 81 , 82 , 59 , 60 ) ( 79 , 80 , 101 , 102 , 103 , 104 , 81 , 82 ) ( 101 , 102 , 123 , 124 , 125 , 126 , 103 , 104 ) ( 123 , 124 , 145 , 146 , 147 , 148 , 125 , 126 ) ( 145 , 146 , 167 , 168 , 169 , 170 , 147 , 148 ) ( 167 , 168 , 189 , 190 , 191 , 192 , 169 , 170 ) ( 189 , 190 , 211 , 212 , 213 , 214 , 191 , 192 ) ( 211 , 212 , 233 , 234 , 235 , 236 , 213 , 214 ) ( 233 , 234 , 255 , 256 , 257 , 258 , 235 , 236 ) ( 255 , 256 , 277 , 278 , 279 , 280 , 257 , 258 ) ( 277 , 278 , 299 , 300 , 301 , 302 , 279 , 280 ) ( 299 , 300 , 321 , 322 , 323 , 324 , 301 , 302 ) ( 321 , 322 , 343 , 344 , 345 , 346 , 323 , 324 ) ( 343 , 344 , 365 , 366 , 367 , 368 , 345 , 346 ) ( 365 , 366 , 387 , 388 , 389 , 390 , 367 , 368 ) ( 387 , 388 , 409 , 410 , 411 , 412 , 389 , 390 ) ( 409 , 410 , 431 , 432 , 433 , 434 , 411 , 412 ) ( 431 , 432 , 453 , 454 , 455 , 456 , 433 , 434 ) ( 453 , 454 , 475 , 476 , 477 , 478 , 455 , 456 ) ( 475 , 476 , 497 , 498 , 499 , 500 , 477 , 478 ) ( 497 , 498 , 519 , 520 , 521 , 522 , 499 , 500 ) ( 519 , 520 , 541 , 542 , 543 , 544 , 521 , 522 ) ( 541 , 542 , 563 , 564 , 565 , 566 , 543 , 544 ) ( 563 , 564 , 585 , 586 , 587 , 588 , 565 , 566 ) ( 585 , 586 , 607 , 608 , 609 , 610 , 587 , 588 ) ( 607 , 608 , 629 , 630 , 631 , 632 , 609 , 610 ) ( 629 , 630 , 651 , 652 , 653 , 654 , 631 , 632 ) ( 651 , 652 , 673 , 674 , 675 , 676 , 653 , 654 ) ( 15 , 16 , 37 , 38 , 39 , 40 , 17 , 18 ) ( 37 , 38 , 59 , 60 , 61 , 62 , 39 , 40 ) ( 59 , 60 , 81 , 82 , 83 , 84 , 61 , 62 ) ( 81 , 82 , 103 , 104 , 105 , 106 , 83 , 84 ) ( 103 , 104 , 125 , 126 , 127 , 128 , 105 , 106 ) ( 125 , 126 , 147 , 148 , 149 , 150 , 127 , 128 ) ( 147 , 148 , 169 , 170 , 171 , 172 , 149 , 150 ) ( 169 , 170 , 191 , 192 , 193 , 194 , 171 , 172 ) ( 191 , 192 , 213 , 214 , 215 , 216 , 193 , 194 ) ( 213 , 214 , 235 , 236 , 237 , 238 , 215 , 216 ) ( 235 , 236 , 257 , 258 , 259 , 260 , 237 , 238 ) ( 257 , 258 , 279 , 280 , 281 , 282 , 259 , 260 ) ( 279 , 280 , 301 , 302 , 303 , 304 , 281 , 282 ) ( 301 , 302 , 323 , 324 , 325 , 326 , 303 , 304 ) ( 323 , 324 , 345 , 346 , 347 , 348 , 325 , 326 ) ( 345 , 346 , 367 , 368 , 369 , 370 , 347 , 348 ) ( 367 , 368 , 389 , 390 , 391 , 392 , 369 , 370 ) ( 389 , 390 , 411 , 412 , 413 , 414 , 391 , 392 ) ( 411 , 412 , 433 , 434 , 435 , 436 , 413 , 414 ) ( 433 , 434 , 455 , 456 , 457 , 458 , 435 , 436 ) ( 455 , 456 , 477 , 478 , 479 , 480 , 457 , 458 ) ( 477 , 478 , 499 , 500 , 501 , 502 , 479 , 480 ) ( 499 , 500 , 521 , 522 , 523 , 524 , 501 , 502 ) ( 521 , 522 , 543 , 544 , 545 , 546 , 523 , 524 ) ( 543 , 544 , 565 , 566 , 567 , 568 , 545 , 546 ) ( 565 , 566 , 587 , 588 , 589 , 590 , 567 , 568 ) ( 587 , 588 , 609 , 610 , 611 , 612 , 589 , 590 ) ( 609 , 610 , 631 , 632 , 633 , 634 , 611 , 612 ) ( 631 , 632 , 653 , 654 , 655 , 656 , 633 , 634 ) ( 653 , 654 , 675 , 676 , 677 , 678 , 655 , 656 ) ( 17 , 18 , 39 , 40 , 41 , 42 , 19 , 20 ) ( 39 , 40 , 61 , 62 , 63 , 64 , 41 , 42 ) ( 61 , 62 , 83 , 84 , 85 , 86 , 63 , 64 ) ( 83 , 84 , 105 , 106 , 107 , 108 , 85 , 86 ) ( 105 , 106 , 127 , 128 , 129 , 130 , 107 , 108 ) ( 127 , 128 , 149 , 150 , 151 , 152 , 129 , 130 ) ( 149 , 150 , 171 , 172 , 173 , 174 , 151 , 152 ) ( 171 , 172 , 193 , 194 , 195 , 196 , 173 , 174 ) ( 193 , 194 , 215 , 216 , 217 , 218 , 195 , 196 ) ( 215 , 216 , 237 , 238 , 239 , 240 , 217 , 218 ) ( 237 , 238 , 259 , 260 , 261 , 262 , 239 , 240 ) ( 259 , 260 , 281 , 282 , 283 , 284 , 261 , 262 ) ( 281 , 282 , 303 , 304 , 305 , 306 , 283 , 284 ) ( 303 , 304 , 325 , 326 , 327 , 328 , 305 , 306 ) ( 325 , 326 , 347 , 348 , 349 , 350 , 327 , 328 ) ( 347 , 348 , 369 , 370 , 371 , 372 , 349 , 350 ) ( 369 , 370 , 391 , 392 , 393 , 394 , 371 , 372 ) ( 391 , 392 , 413 , 414 , 415 , 416 , 393 , 394 ) ( 413 , 414 , 435 , 436 , 437 , 438 , 415 , 416 ) ( 435 , 436 , 457 , 458 , 459 , 460 , 437 , 438 ) ( 457 , 458 , 479 , 480 , 481 , 482 , 459 , 460 ) ( 479 , 480 , 501 , 502 , 503 , 504 , 481 , 482 ) ( 501 , 502 , 523 , 524 , 525 , 526 , 503 , 504 ) ( 523 , 524 , 545 , 546 , 547 , 548 , 525 , 526 ) ( 545 , 546 , 567 , 568 , 569 , 570 , 547 , 548 ) ( 567 , 568 , 589 , 590 , 591 , 592 , 569 , 570 ) ( 589 , 590 , 611 , 612 , 613 , 614 , 591 , 592 ) ( 611 , 612 , 633 , 634 , 635 , 636 , 613 , 614 ) ( 633 , 634 , 655 , 656 , 657 , 658 , 635 , 636 ) ( 655 , 656 , 677 , 678 , 679 , 680 , 657 , 658 ) ( 19 , 20 , 41 , 42 , 43 , 44 , 21 , 22 ) ( 41 , 42 , 63 , 64 , 65 , 66 , 43 , 44 ) ( 63 , 64 , 85 , 86 , 87 , 88 , 65 , 66 ) ( 85 , 86 , 107 , 108 , 109 , 110 , 87 , 88 ) ( 107 , 108 , 129 , 130 , 131 , 132 , 109 , 110 ) ( 129 , 130 , 151 , 152 , 153 , 154 , 131 , 132 ) ( 151 , 152 , 173 , 174 , 175 , 176 , 153 , 154 ) ( 173 , 174 , 195 , 196 , 197 , 198 , 175 , 176 ) ( 195 , 196 , 217 , 218 , 219 , 220 , 197 , 198 ) ( 217 , 218 , 239 , 240 , 241 , 242 , 219 , 220 ) ( 239 , 240 , 261 , 262 , 263 , 264 , 241 , 242 ) ( 261 , 262 , 283 , 284 , 285 , 286 , 263 , 264 ) ( 283 , 284 , 305 , 306 , 307 , 308 , 285 , 286 ) ( 305 , 306 , 327 , 328 , 329 , 330 , 307 , 308 ) ( 327 , 328 , 349 , 350 , 351 , 352 , 329 , 330 ) ( 349 , 350 , 371 , 372 , 373 , 374 , 351 , 352 ) ( 371 , 372 , 393 , 394 , 395 , 396 , 373 , 374 ) ( 393 , 394 , 415 , 416 , 417 , 418 , 395 , 396 ) ( 415 , 416 , 437 , 438 , 439 , 440 , 417 , 418 ) ( 437 , 438 , 459 , 460 , 461 , 462 , 439 , 440 ) ( 459 , 460 , 481 , 482 , 483 , 484 , 461 , 462 ) ( 481 , 482 , 503 , 504 , 505 , 506 , 483 , 484 ) ( 503 , 504 , 525 , 526 , 527 , 528 , 505 , 506 ) ( 525 , 526 , 547 , 548 , 549 , 550 , 527 , 528 ) ( 547 , 548 , 569 , 570 , 571 , 572 , 549 , 550 ) ( 569 , 570 , 591 , 592 , 593 , 594 , 571 , 572 ) ( 591 , 592 , 613 , 614 , 615 , 616 , 593 , 594 ) ( 613 , 614 , 635 , 636 , 637 , 638 , 615 , 616 ) ( 635 , 636 , 657 , 658 , 659 , 660 , 637 , 638 ) ( 657 , 658 , 679 , 680 , 681 , 682 , 659 , 660 ); # ... define loads and support (half MBB-beam) set fixnodes := 1 , 3 , 5 , 7 , 9 , 11 , 13 , 15 , 17 , 19 , 21 , 682 ; let F[2] := -1; # ... entries of element stiffness matrix let k[1] := E / ( 1 - nu^2 ) * ( 1/2 - nu/6 ); let k[2] := E / ( 1 - nu^2 ) * ( 1/8 + nu/8 ); let k[3] := E / ( 1 - nu^2 ) * ( -1/4 - nu/12 ); let k[4] := E / ( 1 - nu^2 ) * ( -1/8 + 3*nu/8 ); let k[5] := E / ( 1 - nu^2 ) * ( -1/4 + nu/12 ); let k[6] := E / ( 1 - nu^2 ) * ( -1/8 - nu/8 ); let k[7] := E / ( 1 - nu^2 ) * ( nu/6 ); let k[8] := E / ( 1 - nu^2 ) * ( 1/8 - 3*nu/8 ); ### set scaling here suffix scale IN; let {j in freenodes} U[j].scale := 1000;