Hi, everybody
I want to create a shell element considering material nonlinear in OS.
Two ways are used to construct $secTag of the shell element as following:
1/define SectionForceDeformation;
2/define some PlateFiberSection, elastic or otherwise;
can you show an example or a process to me for the first way?
for the second way, using the command(section PlateFiber) includes material tag for a previously-defined "plate fiber material" including threeDTag.
my question is what does the threeDtag mean or how to define the three-dimensional material. Can it be replaced by uniaxial material ? i couldn't find the guideline about it in manual.
thank u. looking forward to your reply.
How to model a nonlinear shell element in OS
Moderators: silvia, selimgunay, Moderators
Re: How to model a nonlinear shell element in OS
Hello,anybody answers my questions ?
Re: How to model a nonlinear shell element in OS
to model a nonlinear shell element you can use PlateFiberSection. However, this section has some bugs that are still unresolved and that is why it is not described in the command manual. I would suggest you not to use it until it is fixed.
Re: How to model a nonlinear shell element in OS
vesna wrote:
> to model a nonlinear shell element you can use PlateFiberSection. However,
> this section has some bugs that are still unresolved and that is why it is
> not described in the command manual. I would suggest you not to use it
> until it is fixed.
thank u very much
> to model a nonlinear shell element you can use PlateFiberSection. However,
> this section has some bugs that are still unresolved and that is why it is
> not described in the command manual. I would suggest you not to use it
> until it is fixed.
thank u very much
Re: How to model a nonlinear shell element in OS
vesna wrote:
> to model a nonlinear shell element you can use PlateFiberSection. However,
> this section has some bugs that are still unresolved and that is why it is
> not described in the command manual. I would suggest you not to use it
> until it is fixed.
hi Vesna,
I'd like to know if there is any progress on this matter. I'd like to use the ShellNL element in a geometrically nonlinear analysis, later perhaps with material nonlinearity. I have checked the element on a 1000 mm x 1000 mm recangular plate, single perpendicular concentrated force in the middle. Mesh is 20x20 elements, material is linear elastic. The arrangement is doubly symmetric. ShellMITC4 element was used for the linear (small displacements) and ShellNL for the nonlinear (large displacements) analysis.
Comparison was done by ANSYS Shell181 element (most commonly used 4-node shell element of ANSYS), wich can handle both small and large deformations.
Four cases were investigated: simple support (fix X 1 1 1 0 0 0) or fixed (fix X 1 1 1 1 1 1) along all sides, both with small and large deformations.
In all cases the deflection of the node with the force applied and the deflected shape was checked.
The results are:
OpenSees, single support, ShellMITC4: converges, max. defl: 9.42 mm, shape OK.
OpenSees, single support, ShellNL: matrix reported to be singular, error.
OpenSees, fixed support, ShellMITC4: converges, max. defl: 4.53 mm, shape OK.
OpenSees, fixed support, ShellNL: converges, max. defl: 0.37 mm, non-symmetric deflected shape
ANSYS, single support, Shell181, small deformations: max. defl: 9.43 mm, shape OK.
ANSYS, single support, Shell181, large deformations: max. defl: 4.50 mm, shape OK.
ANSYS, fixed support, Shell181, small deformations: max. defl: 4.55 mm, shape OK.
ANSYS, fixed support, Shell181, large deformations: max. defl: 3.48 mm, shape OK.
I copied the relevant parts of the code here, and included the full source for the model below that.
I must note that due to an error in the source of OS the code does not accept all arguments during object creation: int numData = 6; in line 65 should be int numData = 11; so the .exe must be re-compiled for this to work at all.
I also checked the settings for the solver - they should be OK. Reducing the load and/or step size does not affect the results.
For me the results mean that this element does not work well. I wonder if you could help me with this. I could not find any information on the implemented element, but if you could provide me with these, I'd probably give a shot to find the errors.
The relevant parts of the code for the ShellNL element, fixed supports:
wipe;
set subx 20
set suby 20
model basic -ndm 3 -ndf 6
#defining nodes. distances in [mm]
node 1 0.0000 0.0000 0;
node 2 0.0000 50.0000 0;
node 3 0.0000 100.0000 0;
node 4 0.0000 150.0000 0;
node 5 0.0000 200.0000 0;
node 6 0.0000 250.0000 0;
node 7 0.0000 300.0000 0;
node 8 0.0000 350.0000 0;
node 9 0.0000 400.0000 0;
node 10 0.0000 450.0000 0;
node 11 0.0000 500.0000 0;
node 12 0.0000 550.0000 0;
node 13 0.0000 600.0000 0;
node 14 0.0000 650.0000 0;
node 15 0.0000 700.0000 0;
node 16 0.0000 750.0000 0;
node 17 0.0000 800.0000 0;
node 18 0.0000 850.0000 0;
node 19 0.0000 900.0000 0;
node 20 0.0000 950.0000 0;
node 21 0.0000 1000.0000 0;
node 101 50.0000 0.0000 0;
node 102 50.0000 50.0000 0;
node 103 50.0000 100.0000 0;
.
.
.
node 2018 1000.0000 850.0000 0;
node 2019 1000.0000 900.0000 0;
node 2020 1000.0000 950.0000 0;
node 2021 1000.0000 1000.0000 0;
#slightly modified code from the OS wiki (http://opensees.berkeley.edu/wiki/index.php/ShellNL)
nDMaterial ElasticIsotropic 1 210000. 0.3;#steel
nDMaterial PlateFiber 2 1;
section PlateFiber 3 2 4.;#4mm thickness
#defining elements
element ShellNL 1 1 3 103 201 2 103 202 101 102 3
element ShellNL 3 3 5 105 203 4 105 204 103 104 3
element ShellNL 5 5 7 107 205 6 107 206 105 106 3
.
.
.
element ShellNL 1813 1813 1815 1915 2013 1814 1915 2014 1913 1914 3
element ShellNL 1815 1815 1817 1917 2015 1816 1917 2016 1915 1916 3
element ShellNL 1817 1817 1819 1919 2017 1818 1919 2018 1917 1918 3
element ShellNL 1819 1819 1821 1921 2019 1820 1921 2020 1919 1920 3
#defining constraints
fix 1 1 1 1 1 1 1;
fix 101 1 1 1 1 1 1;
fix 201 1 1 1 1 1 1;
fix 301 1 1 1 1 1 1;
fix 401 1 1 1 1 1 1;
.
.
.
fix 1901 1 1 1 1 1 1;
fix 2001 1 1 1 1 1 1;
fix 21 1 1 1 1 1 1;
fix 121 1 1 1 1 1 1;
fix 221 1 1 1 1 1 1;
.
.
.
fix 1921 1 1 1 1 1 1;
fix 2021 1 1 1 1 1 1;
fix 2 1 1 1 1 1 1;
fix 3 1 1 1 1 1 1;
.
.
.
fix 19 1 1 1 1 1 1;
fix 20 1 1 1 1 1 1;
fix 2002 1 1 1 1 1 1;
fix 2003 1 1 1 1 1 1;
.
.
.
fix 2019 1 1 1 1 1 1;
fix 2020 1 1 1 1 1 1;
#defining loads
pattern Plain 1 Linear {
load 1011 0 0 -1000 0 0 0;#N
}
recorder Node -file {E:\OpenSees\Out\shell_1.out} -time -nodeRange 1 2021 -dof 1 disp
.
.
.
recorder Node -file {E:\OpenSees\Out\shell_6.out} -time -nodeRange 1 2021 -dof 6 disp
#analysis settings
constraints Plain
numberer RCM
system UmfPack -lvalueFact 50
test NormDispIncr 1.0e-6 10
algorithm Newton
integrator Newmark 0.5 0.25
analysis Static
analyze 1
quit
The full code:
#vertical shell element
wipe;
set subx 20
set suby 20
#first example program
model basic -ndm 3 -ndf 6
#defining nodes. distances in [mm]
node 1 0.0000 0.0000 0;
node 2 0.0000 50.0000 0;
node 3 0.0000 100.0000 0;
node 4 0.0000 150.0000 0;
node 5 0.0000 200.0000 0;
node 6 0.0000 250.0000 0;
node 7 0.0000 300.0000 0;
node 8 0.0000 350.0000 0;
node 9 0.0000 400.0000 0;
node 10 0.0000 450.0000 0;
node 11 0.0000 500.0000 0;
node 12 0.0000 550.0000 0;
node 13 0.0000 600.0000 0;
node 14 0.0000 650.0000 0;
node 15 0.0000 700.0000 0;
node 16 0.0000 750.0000 0;
node 17 0.0000 800.0000 0;
node 18 0.0000 850.0000 0;
node 19 0.0000 900.0000 0;
node 20 0.0000 950.0000 0;
node 21 0.0000 1000.0000 0;
node 101 50.0000 0.0000 0;
node 102 50.0000 50.0000 0;
node 103 50.0000 100.0000 0;
node 104 50.0000 150.0000 0;
node 105 50.0000 200.0000 0;
node 106 50.0000 250.0000 0;
node 107 50.0000 300.0000 0;
node 108 50.0000 350.0000 0;
node 109 50.0000 400.0000 0;
node 110 50.0000 450.0000 0;
node 111 50.0000 500.0000 0;
node 112 50.0000 550.0000 0;
node 113 50.0000 600.0000 0;
node 114 50.0000 650.0000 0;
node 115 50.0000 700.0000 0;
node 116 50.0000 750.0000 0;
node 117 50.0000 800.0000 0;
node 118 50.0000 850.0000 0;
node 119 50.0000 900.0000 0;
node 120 50.0000 950.0000 0;
node 121 50.0000 1000.0000 0;
node 201 100.0000 0.0000 0;
node 202 100.0000 50.0000 0;
node 203 100.0000 100.0000 0;
node 204 100.0000 150.0000 0;
node 205 100.0000 200.0000 0;
node 206 100.0000 250.0000 0;
node 207 100.0000 300.0000 0;
node 208 100.0000 350.0000 0;
node 209 100.0000 400.0000 0;
node 210 100.0000 450.0000 0;
node 211 100.0000 500.0000 0;
node 212 100.0000 550.0000 0;
node 213 100.0000 600.0000 0;
node 214 100.0000 650.0000 0;
node 215 100.0000 700.0000 0;
node 216 100.0000 750.0000 0;
node 217 100.0000 800.0000 0;
node 218 100.0000 850.0000 0;
node 219 100.0000 900.0000 0;
node 220 100.0000 950.0000 0;
node 221 100.0000 1000.0000 0;
node 301 150.0000 0.0000 0;
node 302 150.0000 50.0000 0;
node 303 150.0000 100.0000 0;
node 304 150.0000 150.0000 0;
node 305 150.0000 200.0000 0;
node 306 150.0000 250.0000 0;
node 307 150.0000 300.0000 0;
node 308 150.0000 350.0000 0;
node 309 150.0000 400.0000 0;
node 310 150.0000 450.0000 0;
node 311 150.0000 500.0000 0;
node 312 150.0000 550.0000 0;
node 313 150.0000 600.0000 0;
node 314 150.0000 650.0000 0;
node 315 150.0000 700.0000 0;
node 316 150.0000 750.0000 0;
node 317 150.0000 800.0000 0;
node 318 150.0000 850.0000 0;
node 319 150.0000 900.0000 0;
node 320 150.0000 950.0000 0;
node 321 150.0000 1000.0000 0;
node 401 200.0000 0.0000 0;
node 402 200.0000 50.0000 0;
node 403 200.0000 100.0000 0;
node 404 200.0000 150.0000 0;
node 405 200.0000 200.0000 0;
node 406 200.0000 250.0000 0;
node 407 200.0000 300.0000 0;
node 408 200.0000 350.0000 0;
node 409 200.0000 400.0000 0;
node 410 200.0000 450.0000 0;
node 411 200.0000 500.0000 0;
node 412 200.0000 550.0000 0;
node 413 200.0000 600.0000 0;
node 414 200.0000 650.0000 0;
node 415 200.0000 700.0000 0;
node 416 200.0000 750.0000 0;
node 417 200.0000 800.0000 0;
node 418 200.0000 850.0000 0;
node 419 200.0000 900.0000 0;
node 420 200.0000 950.0000 0;
node 421 200.0000 1000.0000 0;
node 501 250.0000 0.0000 0;
node 502 250.0000 50.0000 0;
node 503 250.0000 100.0000 0;
node 504 250.0000 150.0000 0;
node 505 250.0000 200.0000 0;
node 506 250.0000 250.0000 0;
node 507 250.0000 300.0000 0;
node 508 250.0000 350.0000 0;
node 509 250.0000 400.0000 0;
node 510 250.0000 450.0000 0;
node 511 250.0000 500.0000 0;
node 512 250.0000 550.0000 0;
node 513 250.0000 600.0000 0;
node 514 250.0000 650.0000 0;
node 515 250.0000 700.0000 0;
node 516 250.0000 750.0000 0;
node 517 250.0000 800.0000 0;
node 518 250.0000 850.0000 0;
node 519 250.0000 900.0000 0;
node 520 250.0000 950.0000 0;
node 521 250.0000 1000.0000 0;
node 601 300.0000 0.0000 0;
node 602 300.0000 50.0000 0;
node 603 300.0000 100.0000 0;
node 604 300.0000 150.0000 0;
node 605 300.0000 200.0000 0;
node 606 300.0000 250.0000 0;
node 607 300.0000 300.0000 0;
node 608 300.0000 350.0000 0;
node 609 300.0000 400.0000 0;
node 610 300.0000 450.0000 0;
node 611 300.0000 500.0000 0;
node 612 300.0000 550.0000 0;
node 613 300.0000 600.0000 0;
node 614 300.0000 650.0000 0;
node 615 300.0000 700.0000 0;
node 616 300.0000 750.0000 0;
node 617 300.0000 800.0000 0;
node 618 300.0000 850.0000 0;
node 619 300.0000 900.0000 0;
node 620 300.0000 950.0000 0;
node 621 300.0000 1000.0000 0;
node 701 350.0000 0.0000 0;
node 702 350.0000 50.0000 0;
node 703 350.0000 100.0000 0;
node 704 350.0000 150.0000 0;
node 705 350.0000 200.0000 0;
node 706 350.0000 250.0000 0;
node 707 350.0000 300.0000 0;
node 708 350.0000 350.0000 0;
node 709 350.0000 400.0000 0;
node 710 350.0000 450.0000 0;
node 711 350.0000 500.0000 0;
node 712 350.0000 550.0000 0;
node 713 350.0000 600.0000 0;
node 714 350.0000 650.0000 0;
node 715 350.0000 700.0000 0;
node 716 350.0000 750.0000 0;
node 717 350.0000 800.0000 0;
node 718 350.0000 850.0000 0;
node 719 350.0000 900.0000 0;
node 720 350.0000 950.0000 0;
node 721 350.0000 1000.0000 0;
node 801 400.0000 0.0000 0;
node 802 400.0000 50.0000 0;
node 803 400.0000 100.0000 0;
node 804 400.0000 150.0000 0;
node 805 400.0000 200.0000 0;
node 806 400.0000 250.0000 0;
node 807 400.0000 300.0000 0;
node 808 400.0000 350.0000 0;
node 809 400.0000 400.0000 0;
node 810 400.0000 450.0000 0;
node 811 400.0000 500.0000 0;
node 812 400.0000 550.0000 0;
node 813 400.0000 600.0000 0;
node 814 400.0000 650.0000 0;
node 815 400.0000 700.0000 0;
node 816 400.0000 750.0000 0;
node 817 400.0000 800.0000 0;
node 818 400.0000 850.0000 0;
node 819 400.0000 900.0000 0;
node 820 400.0000 950.0000 0;
node 821 400.0000 1000.0000 0;
node 901 450.0000 0.0000 0;
node 902 450.0000 50.0000 0;
node 903 450.0000 100.0000 0;
node 904 450.0000 150.0000 0;
node 905 450.0000 200.0000 0;
node 906 450.0000 250.0000 0;
node 907 450.0000 300.0000 0;
node 908 450.0000 350.0000 0;
node 909 450.0000 400.0000 0;
node 910 450.0000 450.0000 0;
node 911 450.0000 500.0000 0;
node 912 450.0000 550.0000 0;
node 913 450.0000 600.0000 0;
node 914 450.0000 650.0000 0;
node 915 450.0000 700.0000 0;
node 916 450.0000 750.0000 0;
node 917 450.0000 800.0000 0;
node 918 450.0000 850.0000 0;
node 919 450.0000 900.0000 0;
node 920 450.0000 950.0000 0;
node 921 450.0000 1000.0000 0;
node 1001 500.0000 0.0000 0;
node 1002 500.0000 50.0000 0;
node 1003 500.0000 100.0000 0;
node 1004 500.0000 150.0000 0;
node 1005 500.0000 200.0000 0;
node 1006 500.0000 250.0000 0;
node 1007 500.0000 300.0000 0;
node 1008 500.0000 350.0000 0;
node 1009 500.0000 400.0000 0;
node 1010 500.0000 450.0000 0;
node 1011 500.0000 500.0000 0;
node 1012 500.0000 550.0000 0;
node 1013 500.0000 600.0000 0;
node 1014 500.0000 650.0000 0;
node 1015 500.0000 700.0000 0;
node 1016 500.0000 750.0000 0;
node 1017 500.0000 800.0000 0;
node 1018 500.0000 850.0000 0;
node 1019 500.0000 900.0000 0;
node 1020 500.0000 950.0000 0;
node 1021 500.0000 1000.0000 0;
node 1101 550.0000 0.0000 0;
node 1102 550.0000 50.0000 0;
node 1103 550.0000 100.0000 0;
node 1104 550.0000 150.0000 0;
node 1105 550.0000 200.0000 0;
node 1106 550.0000 250.0000 0;
node 1107 550.0000 300.0000 0;
node 1108 550.0000 350.0000 0;
node 1109 550.0000 400.0000 0;
node 1110 550.0000 450.0000 0;
node 1111 550.0000 500.0000 0;
node 1112 550.0000 550.0000 0;
node 1113 550.0000 600.0000 0;
node 1114 550.0000 650.0000 0;
node 1115 550.0000 700.0000 0;
node 1116 550.0000 750.0000 0;
node 1117 550.0000 800.0000 0;
node 1118 550.0000 850.0000 0;
node 1119 550.0000 900.0000 0;
node 1120 550.0000 950.0000 0;
node 1121 550.0000 1000.0000 0;
node 1201 600.0000 0.0000 0;
node 1202 600.0000 50.0000 0;
node 1203 600.0000 100.0000 0;
node 1204 600.0000 150.0000 0;
node 1205 600.0000 200.0000 0;
node 1206 600.0000 250.0000 0;
node 1207 600.0000 300.0000 0;
node 1208 600.0000 350.0000 0;
node 1209 600.0000 400.0000 0;
node 1210 600.0000 450.0000 0;
node 1211 600.0000 500.0000 0;
node 1212 600.0000 550.0000 0;
node 1213 600.0000 600.0000 0;
node 1214 600.0000 650.0000 0;
node 1215 600.0000 700.0000 0;
node 1216 600.0000 750.0000 0;
node 1217 600.0000 800.0000 0;
node 1218 600.0000 850.0000 0;
node 1219 600.0000 900.0000 0;
node 1220 600.0000 950.0000 0;
node 1221 600.0000 1000.0000 0;
node 1301 650.0000 0.0000 0;
node 1302 650.0000 50.0000 0;
node 1303 650.0000 100.0000 0;
node 1304 650.0000 150.0000 0;
node 1305 650.0000 200.0000 0;
node 1306 650.0000 250.0000 0;
node 1307 650.0000 300.0000 0;
node 1308 650.0000 350.0000 0;
node 1309 650.0000 400.0000 0;
node 1310 650.0000 450.0000 0;
node 1311 650.0000 500.0000 0;
node 1312 650.0000 550.0000 0;
node 1313 650.0000 600.0000 0;
node 1314 650.0000 650.0000 0;
node 1315 650.0000 700.0000 0;
node 1316 650.0000 750.0000 0;
node 1317 650.0000 800.0000 0;
node 1318 650.0000 850.0000 0;
node 1319 650.0000 900.0000 0;
node 1320 650.0000 950.0000 0;
node 1321 650.0000 1000.0000 0;
node 1401 700.0000 0.0000 0;
node 1402 700.0000 50.0000 0;
node 1403 700.0000 100.0000 0;
node 1404 700.0000 150.0000 0;
node 1405 700.0000 200.0000 0;
node 1406 700.0000 250.0000 0;
node 1407 700.0000 300.0000 0;
node 1408 700.0000 350.0000 0;
node 1409 700.0000 400.0000 0;
node 1410 700.0000 450.0000 0;
node 1411 700.0000 500.0000 0;
node 1412 700.0000 550.0000 0;
node 1413 700.0000 600.0000 0;
node 1414 700.0000 650.0000 0;
node 1415 700.0000 700.0000 0;
node 1416 700.0000 750.0000 0;
node 1417 700.0000 800.0000 0;
node 1418 700.0000 850.0000 0;
node 1419 700.0000 900.0000 0;
node 1420 700.0000 950.0000 0;
node 1421 700.0000 1000.0000 0;
node 1501 750.0000 0.0000 0;
node 1502 750.0000 50.0000 0;
node 1503 750.0000 100.0000 0;
node 1504 750.0000 150.0000 0;
node 1505 750.0000 200.0000 0;
node 1506 750.0000 250.0000 0;
node 1507 750.0000 300.0000 0;
node 1508 750.0000 350.0000 0;
node 1509 750.0000 400.0000 0;
node 1510 750.0000 450.0000 0;
node 1511 750.0000 500.0000 0;
node 1512 750.0000 550.0000 0;
node 1513 750.0000 600.0000 0;
node 1514 750.0000 650.0000 0;
node 1515 750.0000 700.0000 0;
node 1516 750.0000 750.0000 0;
node 1517 750.0000 800.0000 0;
node 1518 750.0000 850.0000 0;
node 1519 750.0000 900.0000 0;
node 1520 750.0000 950.0000 0;
node 1521 750.0000 1000.0000 0;
node 1601 800.0000 0.0000 0;
node 1602 800.0000 50.0000 0;
node 1603 800.0000 100.0000 0;
node 1604 800.0000 150.0000 0;
node 1605 800.0000 200.0000 0;
node 1606 800.0000 250.0000 0;
node 1607 800.0000 300.0000 0;
node 1608 800.0000 350.0000 0;
node 1609 800.0000 400.0000 0;
node 1610 800.0000 450.0000 0;
node 1611 800.0000 500.0000 0;
node 1612 800.0000 550.0000 0;
node 1613 800.0000 600.0000 0;
node 1614 800.0000 650.0000 0;
node 1615 800.0000 700.0000 0;
node 1616 800.0000 750.0000 0;
node 1617 800.0000 800.0000 0;
node 1618 800.0000 850.0000 0;
node 1619 800.0000 900.0000 0;
node 1620 800.0000 950.0000 0;
node 1621 800.0000 1000.0000 0;
node 1701 850.0000 0.0000 0;
node 1702 850.0000 50.0000 0;
node 1703 850.0000 100.0000 0;
node 1704 850.0000 150.0000 0;
node 1705 850.0000 200.0000 0;
node 1706 850.0000 250.0000 0;
node 1707 850.0000 300.0000 0;
node 1708 850.0000 350.0000 0;
node 1709 850.0000 400.0000 0;
node 1710 850.0000 450.0000 0;
node 1711 850.0000 500.0000 0;
node 1712 850.0000 550.0000 0;
node 1713 850.0000 600.0000 0;
node 1714 850.0000 650.0000 0;
node 1715 850.0000 700.0000 0;
node 1716 850.0000 750.0000 0;
node 1717 850.0000 800.0000 0;
node 1718 850.0000 850.0000 0;
node 1719 850.0000 900.0000 0;
node 1720 850.0000 950.0000 0;
node 1721 850.0000 1000.0000 0;
node 1801 900.0000 0.0000 0;
node 1802 900.0000 50.0000 0;
node 1803 900.0000 100.0000 0;
node 1804 900.0000 150.0000 0;
node 1805 900.0000 200.0000 0;
node 1806 900.0000 250.0000 0;
node 1807 900.0000 300.0000 0;
node 1808 900.0000 350.0000 0;
node 1809 900.0000 400.0000 0;
node 1810 900.0000 450.0000 0;
node 1811 900.0000 500.0000 0;
node 1812 900.0000 550.0000 0;
node 1813 900.0000 600.0000 0;
node 1814 900.0000 650.0000 0;
node 1815 900.0000 700.0000 0;
node 1816 900.0000 750.0000 0;
node 1817 900.0000 800.0000 0;
node 1818 900.0000 850.0000 0;
node 1819 900.0000 900.0000 0;
node 1820 900.0000 950.0000 0;
node 1821 900.0000 1000.0000 0;
node 1901 950.0000 0.0000 0;
node 1902 950.0000 50.0000 0;
node 1903 950.0000 100.0000 0;
node 1904 950.0000 150.0000 0;
node 1905 950.0000 200.0000 0;
node 1906 950.0000 250.0000 0;
node 1907 950.0000 300.0000 0;
node 1908 950.0000 350.0000 0;
node 1909 950.0000 400.0000 0;
node 1910 950.0000 450.0000 0;
node 1911 950.0000 500.0000 0;
node 1912 950.0000 550.0000 0;
node 1913 950.0000 600.0000 0;
node 1914 950.0000 650.0000 0;
node 1915 950.0000 700.0000 0;
node 1916 950.0000 750.0000 0;
node 1917 950.0000 800.0000 0;
node 1918 950.0000 850.0000 0;
node 1919 950.0000 900.0000 0;
node 1920 950.0000 950.0000 0;
node 1921 950.0000 1000.0000 0;
node 2001 1000.0000 0.0000 0;
node 2002 1000.0000 50.0000 0;
node 2003 1000.0000 100.0000 0;
node 2004 1000.0000 150.0000 0;
node 2005 1000.0000 200.0000 0;
node 2006 1000.0000 250.0000 0;
node 2007 1000.0000 300.0000 0;
node 2008 1000.0000 350.0000 0;
node 2009 1000.0000 400.0000 0;
node 2010 1000.0000 450.0000 0;
node 2011 1000.0000 500.0000 0;
node 2012 1000.0000 550.0000 0;
node 2013 1000.0000 600.0000 0;
node 2014 1000.0000 650.0000 0;
node 2015 1000.0000 700.0000 0;
node 2016 1000.0000 750.0000 0;
node 2017 1000.0000 800.0000 0;
node 2018 1000.0000 850.0000 0;
node 2019 1000.0000 900.0000 0;
node 2020 1000.0000 950.0000 0;
node 2021 1000.0000 1000.0000 0;
#code from the OS wiki (http://opensees.berkeley.edu/wiki/index.php/ShellNL)
nDMaterial ElasticIsotropic 1 210000. 0.3;#steel
nDMaterial PlateFiber 2 1
section PlateFiber 3 2 4.;#4 mm thickness
#defining elements
element ShellNL 1 1 3 103 201 2 103 202 101 102 3
element ShellNL 3 3 5 105 203 4 105 204 103 104 3
element ShellNL 5 5 7 107 205 6 107 206 105 106 3
element ShellNL 7 7 9 109 207 8 109 208 107 108 3
element ShellNL 9 9 11 111 209 10 111 210 109 110 3
element ShellNL 11 11 13 113 211 12 113 212 111 112 3
element ShellNL 13 13 15 115 213 14 115 214 113 114 3
element ShellNL 15 15 17 117 215 16 117 216 115 116 3
element ShellNL 17 17 19 119 217 18 119 218 117 118 3
element ShellNL 19 19 21 121 219 20 121 220 119 120 3
element ShellNL 201 201 203 303 401 202 303 402 301 302 3
element ShellNL 203 203 205 305 403 204 305 404 303 304 3
element ShellNL 205 205 207 307 405 206 307 406 305 306 3
element ShellNL 207 207 209 309 407 208 309 408 307 308 3
element ShellNL 209 209 211 311 409 210 311 410 309 310 3
element ShellNL 211 211 213 313 411 212 313 412 311 312 3
element ShellNL 213 213 215 315 413 214 315 414 313 314 3
element ShellNL 215 215 217 317 415 216 317 416 315 316 3
element ShellNL 217 217 219 319 417 218 319 418 317 318 3
element ShellNL 219 219 221 321 419 220 321 420 319 320 3
element ShellNL 401 401 403 503 601 402 503 602 501 502 3
element ShellNL 403 403 405 505 603 404 505 604 503 504 3
element ShellNL 405 405 407 507 605 406 507 606 505 506 3
element ShellNL 407 407 409 509 607 408 509 608 507 508 3
element ShellNL 409 409 411 511 609 410 511 610 509 510 3
element ShellNL 411 411 413 513 611 412 513 612 511 512 3
element ShellNL 413 413 415 515 613 414 515 614 513 514 3
element ShellNL 415 415 417 517 615 416 517 616 515 516 3
element ShellNL 417 417 419 519 617 418 519 618 517 518 3
element ShellNL 419 419 421 521 619 420 521 620 519 520 3
element ShellNL 601 601 603 703 801 602 703 802 701 702 3
element ShellNL 603 603 605 705 803 604 705 804 703 704 3
element ShellNL 605 605 607 707 805 606 707 806 705 706 3
element ShellNL 607 607 609 709 807 608 709 808 707 708 3
element ShellNL 609 609 611 711 809 610 711 810 709 710 3
element ShellNL 611 611 613 713 811 612 713 812 711 712 3
element ShellNL 613 613 615 715 813 614 715 814 713 714 3
element ShellNL 615 615 617 717 815 616 717 816 715 716 3
element ShellNL 617 617 619 719 817 618 719 818 717 718 3
element ShellNL 619 619 621 721 819 620 721 820 719 720 3
element ShellNL 801 801 803 903 1001 802 903 1002 901 902 3
element ShellNL 803 803 805 905 1003 804 905 1004 903 904 3
element ShellNL 805 805 807 907 1005 806 907 1006 905 906 3
element ShellNL 807 807 809 909 1007 808 909 1008 907 908 3
element ShellNL 809 809 811 911 1009 810 911 1010 909 910 3
element ShellNL 811 811 813 913 1011 812 913 1012 911 912 3
element ShellNL 813 813 815 915 1013 814 915 1014 913 914 3
element ShellNL 815 815 817 917 1015 816 917 1016 915 916 3
element ShellNL 817 817 819 919 1017 818 919 1018 917 918 3
element ShellNL 819 819 821 921 1019 820 921 1020 919 920 3
element ShellNL 1001 1001 1003 1103 1201 1002 1103 1202 1101 1102 3
element ShellNL 1003 1003 1005 1105 1203 1004 1105 1204 1103 1104 3
element ShellNL 1005 1005 1007 1107 1205 1006 1107 1206 1105 1106 3
element ShellNL 1007 1007 1009 1109 1207 1008 1109 1208 1107 1108 3
element ShellNL 1009 1009 1011 1111 1209 1010 1111 1210 1109 1110 3
element ShellNL 1011 1011 1013 1113 1211 1012 1113 1212 1111 1112 3
element ShellNL 1013 1013 1015 1115 1213 1014 1115 1214 1113 1114 3
element ShellNL 1015 1015 1017 1117 1215 1016 1117 1216 1115 1116 3
element ShellNL 1017 1017 1019 1119 1217 1018 1119 1218 1117 1118 3
element ShellNL 1019 1019 1021 1121 1219 1020 1121 1220 1119 1120 3
element ShellNL 1201 1201 1203 1303 1401 1202 1303 1402 1301 1302 3
element ShellNL 1203 1203 1205 1305 1403 1204 1305 1404 1303 1304 3
element ShellNL 1205 1205 1207 1307 1405 1206 1307 1406 1305 1306 3
element ShellNL 1207 1207 1209 1309 1407 1208 1309 1408 1307 1308 3
element ShellNL 1209 1209 1211 1311 1409 1210 1311 1410 1309 1310 3
element ShellNL 1211 1211 1213 1313 1411 1212 1313 1412 1311 1312 3
element ShellNL 1213 1213 1215 1315 1413 1214 1315 1414 1313 1314 3
element ShellNL 1215 1215 1217 1317 1415 1216 1317 1416 1315 1316 3
element ShellNL 1217 1217 1219 1319 1417 1218 1319 1418 1317 1318 3
element ShellNL 1219 1219 1221 1321 1419 1220 1321 1420 1319 1320 3
element ShellNL 1401 1401 1403 1503 1601 1402 1503 1602 1501 1502 3
element ShellNL 1403 1403 1405 1505 1603 1404 1505 1604 1503 1504 3
element ShellNL 1405 1405 1407 1507 1605 1406 1507 1606 1505 1506 3
element ShellNL 1407 1407 1409 1509 1607 1408 1509 1608 1507 1508 3
element ShellNL 1409 1409 1411 1511 1609 1410 1511 1610 1509 1510 3
element ShellNL 1411 1411 1413 1513 1611 1412 1513 1612 1511 1512 3
element ShellNL 1413 1413 1415 1515 1613 1414 1515 1614 1513 1514 3
element ShellNL 1415 1415 1417 1517 1615 1416 1517 1616 1515 1516 3
element ShellNL 1417 1417 1419 1519 1617 1418 1519 1618 1517 1518 3
element ShellNL 1419 1419 1421 1521 1619 1420 1521 1620 1519 1520 3
element ShellNL 1601 1601 1603 1703 1801 1602 1703 1802 1701 1702 3
element ShellNL 1603 1603 1605 1705 1803 1604 1705 1804 1703 1704 3
element ShellNL 1605 1605 1607 1707 1805 1606 1707 1806 1705 1706 3
element ShellNL 1607 1607 1609 1709 1807 1608 1709 1808 1707 1708 3
element ShellNL 1609 1609 1611 1711 1809 1610 1711 1810 1709 1710 3
element ShellNL 1611 1611 1613 1713 1811 1612 1713 1812 1711 1712 3
element ShellNL 1613 1613 1615 1715 1813 1614 1715 1814 1713 1714 3
element ShellNL 1615 1615 1617 1717 1815 1616 1717 1816 1715 1716 3
element ShellNL 1617 1617 1619 1719 1817 1618 1719 1818 1717 1718 3
element ShellNL 1619 1619 1621 1721 1819 1620 1721 1820 1719 1720 3
element ShellNL 1801 1801 1803 1903 2001 1802 1903 2002 1901 1902 3
element ShellNL 1803 1803 1805 1905 2003 1804 1905 2004 1903 1904 3
element ShellNL 1805 1805 1807 1907 2005 1806 1907 2006 1905 1906 3
element ShellNL 1807 1807 1809 1909 2007 1808 1909 2008 1907 1908 3
element ShellNL 1809 1809 1811 1911 2009 1810 1911 2010 1909 1910 3
element ShellNL 1811 1811 1813 1913 2011 1812 1913 2012 1911 1912 3
element ShellNL 1813 1813 1815 1915 2013 1814 1915 2014 1913 1914 3
element ShellNL 1815 1815 1817 1917 2015 1816 1917 2016 1915 1916 3
element ShellNL 1817 1817 1819 1919 2017 1818 1919 2018 1917 1918 3
element ShellNL 1819 1819 1821 1921 2019 1820 1921 2020 1919 1920 3
#defining constraints
fix 1 1 1 1 0 0 0;
fix 101 1 1 1 0 0 0;
fix 201 1 1 1 0 0 0;
fix 301 1 1 1 0 0 0;
fix 401 1 1 1 0 0 0;
fix 501 1 1 1 0 0 0;
fix 601 1 1 1 0 0 0;
fix 701 1 1 1 0 0 0;
fix 801 1 1 1 0 0 0;
fix 901 1 1 1 0 0 0;
fix 1001 1 1 1 0 0 0;
fix 1101 1 1 1 0 0 0;
fix 1201 1 1 1 0 0 0;
fix 1301 1 1 1 0 0 0;
fix 1401 1 1 1 0 0 0;
fix 1501 1 1 1 0 0 0;
fix 1601 1 1 1 0 0 0;
fix 1701 1 1 1 0 0 0;
fix 1801 1 1 1 0 0 0;
fix 1901 1 1 1 0 0 0;
fix 2001 1 1 1 0 0 0;
fix 21 1 1 1 0 0 0;
fix 121 1 1 1 0 0 0;
fix 221 1 1 1 0 0 0;
fix 321 1 1 1 0 0 0;
fix 421 1 1 1 0 0 0;
fix 521 1 1 1 0 0 0;
fix 621 1 1 1 0 0 0;
fix 721 1 1 1 0 0 0;
fix 821 1 1 1 0 0 0;
fix 921 1 1 1 0 0 0;
fix 1021 1 1 1 0 0 0;
fix 1121 1 1 1 0 0 0;
fix 1221 1 1 1 0 0 0;
fix 1321 1 1 1 0 0 0;
fix 1421 1 1 1 0 0 0;
fix 1521 1 1 1 0 0 0;
fix 1621 1 1 1 0 0 0;
fix 1721 1 1 1 0 0 0;
fix 1821 1 1 1 0 0 0;
fix 1921 1 1 1 0 0 0;
fix 2021 1 1 1 0 0 0;
fix 2 1 1 1 0 0 0;
fix 3 1 1 1 0 0 0;
fix 4 1 1 1 0 0 0;
fix 5 1 1 1 0 0 0;
fix 6 1 1 1 0 0 0;
fix 7 1 1 1 0 0 0;
fix 8 1 1 1 0 0 0;
fix 9 1 1 1 0 0 0;
fix 10 1 1 1 0 0 0;
fix 11 1 1 1 0 0 0;
fix 12 1 1 1 0 0 0;
fix 13 1 1 1 0 0 0;
fix 14 1 1 1 0 0 0;
fix 15 1 1 1 0 0 0;
fix 16 1 1 1 0 0 0;
fix 17 1 1 1 0 0 0;
fix 18 1 1 1 0 0 0;
fix 19 1 1 1 0 0 0;
fix 20 1 1 1 0 0 0;
fix 2002 1 1 1 0 0 0;
fix 2003 1 1 1 0 0 0;
fix 2004 1 1 1 0 0 0;
fix 2005 1 1 1 0 0 0;
fix 2006 1 1 1 0 0 0;
fix 2007 1 1 1 0 0 0;
fix 2008 1 1 1 0 0 0;
fix 2009 1 1 1 0 0 0;
fix 2010 1 1 1 0 0 0;
fix 2011 1 1 1 0 0 0;
fix 2012 1 1 1 0 0 0;
fix 2013 1 1 1 0 0 0;
fix 2014 1 1 1 0 0 0;
fix 2015 1 1 1 0 0 0;
fix 2016 1 1 1 0 0 0;
fix 2017 1 1 1 0 0 0;
fix 2018 1 1 1 0 0 0;
fix 2019 1 1 1 0 0 0;
fix 2020 1 1 1 0 0 0;
#defining loads
pattern Plain 1 Linear {
load 1011 0 0 -1000 0 0 0;#N
}
recorder Node -file {E:\OpenSees\Out\shell_1.out} -time -nodeRange 1 2021 -dof 1 disp
recorder Node -file {E:\OpenSees\Out\shell_2.out} -time -nodeRange 1 2021 -dof 2 disp
recorder Node -file {E:\OpenSees\Out\shell_3.out} -time -nodeRange 1 2021 -dof 3 disp
recorder Node -file {E:\OpenSees\Out\shell_4.out} -time -nodeRange 1 2021 -dof 4 disp
recorder Node -file {E:\OpenSees\Out\shell_5.out} -time -nodeRange 1 2021 -dof 5 disp
recorder Node -file {E:\OpenSees\Out\shell_6.out} -time -nodeRange 1 2021 -dof 6 disp
#analysis settings
constraints Plain
numberer RCM
system UmfPack -lvalueFact 50
test NormDispIncr 1.0e-6 100
algorithm Newton
integrator Newmark 0.5 0.25
analysis Static
analyze 1
> to model a nonlinear shell element you can use PlateFiberSection. However,
> this section has some bugs that are still unresolved and that is why it is
> not described in the command manual. I would suggest you not to use it
> until it is fixed.
hi Vesna,
I'd like to know if there is any progress on this matter. I'd like to use the ShellNL element in a geometrically nonlinear analysis, later perhaps with material nonlinearity. I have checked the element on a 1000 mm x 1000 mm recangular plate, single perpendicular concentrated force in the middle. Mesh is 20x20 elements, material is linear elastic. The arrangement is doubly symmetric. ShellMITC4 element was used for the linear (small displacements) and ShellNL for the nonlinear (large displacements) analysis.
Comparison was done by ANSYS Shell181 element (most commonly used 4-node shell element of ANSYS), wich can handle both small and large deformations.
Four cases were investigated: simple support (fix X 1 1 1 0 0 0) or fixed (fix X 1 1 1 1 1 1) along all sides, both with small and large deformations.
In all cases the deflection of the node with the force applied and the deflected shape was checked.
The results are:
OpenSees, single support, ShellMITC4: converges, max. defl: 9.42 mm, shape OK.
OpenSees, single support, ShellNL: matrix reported to be singular, error.
OpenSees, fixed support, ShellMITC4: converges, max. defl: 4.53 mm, shape OK.
OpenSees, fixed support, ShellNL: converges, max. defl: 0.37 mm, non-symmetric deflected shape
ANSYS, single support, Shell181, small deformations: max. defl: 9.43 mm, shape OK.
ANSYS, single support, Shell181, large deformations: max. defl: 4.50 mm, shape OK.
ANSYS, fixed support, Shell181, small deformations: max. defl: 4.55 mm, shape OK.
ANSYS, fixed support, Shell181, large deformations: max. defl: 3.48 mm, shape OK.
I copied the relevant parts of the code here, and included the full source for the model below that.
I must note that due to an error in the source of OS the code does not accept all arguments during object creation: int numData = 6; in line 65 should be int numData = 11; so the .exe must be re-compiled for this to work at all.
I also checked the settings for the solver - they should be OK. Reducing the load and/or step size does not affect the results.
For me the results mean that this element does not work well. I wonder if you could help me with this. I could not find any information on the implemented element, but if you could provide me with these, I'd probably give a shot to find the errors.
The relevant parts of the code for the ShellNL element, fixed supports:
wipe;
set subx 20
set suby 20
model basic -ndm 3 -ndf 6
#defining nodes. distances in [mm]
node 1 0.0000 0.0000 0;
node 2 0.0000 50.0000 0;
node 3 0.0000 100.0000 0;
node 4 0.0000 150.0000 0;
node 5 0.0000 200.0000 0;
node 6 0.0000 250.0000 0;
node 7 0.0000 300.0000 0;
node 8 0.0000 350.0000 0;
node 9 0.0000 400.0000 0;
node 10 0.0000 450.0000 0;
node 11 0.0000 500.0000 0;
node 12 0.0000 550.0000 0;
node 13 0.0000 600.0000 0;
node 14 0.0000 650.0000 0;
node 15 0.0000 700.0000 0;
node 16 0.0000 750.0000 0;
node 17 0.0000 800.0000 0;
node 18 0.0000 850.0000 0;
node 19 0.0000 900.0000 0;
node 20 0.0000 950.0000 0;
node 21 0.0000 1000.0000 0;
node 101 50.0000 0.0000 0;
node 102 50.0000 50.0000 0;
node 103 50.0000 100.0000 0;
.
.
.
node 2018 1000.0000 850.0000 0;
node 2019 1000.0000 900.0000 0;
node 2020 1000.0000 950.0000 0;
node 2021 1000.0000 1000.0000 0;
#slightly modified code from the OS wiki (http://opensees.berkeley.edu/wiki/index.php/ShellNL)
nDMaterial ElasticIsotropic 1 210000. 0.3;#steel
nDMaterial PlateFiber 2 1;
section PlateFiber 3 2 4.;#4mm thickness
#defining elements
element ShellNL 1 1 3 103 201 2 103 202 101 102 3
element ShellNL 3 3 5 105 203 4 105 204 103 104 3
element ShellNL 5 5 7 107 205 6 107 206 105 106 3
.
.
.
element ShellNL 1813 1813 1815 1915 2013 1814 1915 2014 1913 1914 3
element ShellNL 1815 1815 1817 1917 2015 1816 1917 2016 1915 1916 3
element ShellNL 1817 1817 1819 1919 2017 1818 1919 2018 1917 1918 3
element ShellNL 1819 1819 1821 1921 2019 1820 1921 2020 1919 1920 3
#defining constraints
fix 1 1 1 1 1 1 1;
fix 101 1 1 1 1 1 1;
fix 201 1 1 1 1 1 1;
fix 301 1 1 1 1 1 1;
fix 401 1 1 1 1 1 1;
.
.
.
fix 1901 1 1 1 1 1 1;
fix 2001 1 1 1 1 1 1;
fix 21 1 1 1 1 1 1;
fix 121 1 1 1 1 1 1;
fix 221 1 1 1 1 1 1;
.
.
.
fix 1921 1 1 1 1 1 1;
fix 2021 1 1 1 1 1 1;
fix 2 1 1 1 1 1 1;
fix 3 1 1 1 1 1 1;
.
.
.
fix 19 1 1 1 1 1 1;
fix 20 1 1 1 1 1 1;
fix 2002 1 1 1 1 1 1;
fix 2003 1 1 1 1 1 1;
.
.
.
fix 2019 1 1 1 1 1 1;
fix 2020 1 1 1 1 1 1;
#defining loads
pattern Plain 1 Linear {
load 1011 0 0 -1000 0 0 0;#N
}
recorder Node -file {E:\OpenSees\Out\shell_1.out} -time -nodeRange 1 2021 -dof 1 disp
.
.
.
recorder Node -file {E:\OpenSees\Out\shell_6.out} -time -nodeRange 1 2021 -dof 6 disp
#analysis settings
constraints Plain
numberer RCM
system UmfPack -lvalueFact 50
test NormDispIncr 1.0e-6 10
algorithm Newton
integrator Newmark 0.5 0.25
analysis Static
analyze 1
quit
The full code:
#vertical shell element
wipe;
set subx 20
set suby 20
#first example program
model basic -ndm 3 -ndf 6
#defining nodes. distances in [mm]
node 1 0.0000 0.0000 0;
node 2 0.0000 50.0000 0;
node 3 0.0000 100.0000 0;
node 4 0.0000 150.0000 0;
node 5 0.0000 200.0000 0;
node 6 0.0000 250.0000 0;
node 7 0.0000 300.0000 0;
node 8 0.0000 350.0000 0;
node 9 0.0000 400.0000 0;
node 10 0.0000 450.0000 0;
node 11 0.0000 500.0000 0;
node 12 0.0000 550.0000 0;
node 13 0.0000 600.0000 0;
node 14 0.0000 650.0000 0;
node 15 0.0000 700.0000 0;
node 16 0.0000 750.0000 0;
node 17 0.0000 800.0000 0;
node 18 0.0000 850.0000 0;
node 19 0.0000 900.0000 0;
node 20 0.0000 950.0000 0;
node 21 0.0000 1000.0000 0;
node 101 50.0000 0.0000 0;
node 102 50.0000 50.0000 0;
node 103 50.0000 100.0000 0;
node 104 50.0000 150.0000 0;
node 105 50.0000 200.0000 0;
node 106 50.0000 250.0000 0;
node 107 50.0000 300.0000 0;
node 108 50.0000 350.0000 0;
node 109 50.0000 400.0000 0;
node 110 50.0000 450.0000 0;
node 111 50.0000 500.0000 0;
node 112 50.0000 550.0000 0;
node 113 50.0000 600.0000 0;
node 114 50.0000 650.0000 0;
node 115 50.0000 700.0000 0;
node 116 50.0000 750.0000 0;
node 117 50.0000 800.0000 0;
node 118 50.0000 850.0000 0;
node 119 50.0000 900.0000 0;
node 120 50.0000 950.0000 0;
node 121 50.0000 1000.0000 0;
node 201 100.0000 0.0000 0;
node 202 100.0000 50.0000 0;
node 203 100.0000 100.0000 0;
node 204 100.0000 150.0000 0;
node 205 100.0000 200.0000 0;
node 206 100.0000 250.0000 0;
node 207 100.0000 300.0000 0;
node 208 100.0000 350.0000 0;
node 209 100.0000 400.0000 0;
node 210 100.0000 450.0000 0;
node 211 100.0000 500.0000 0;
node 212 100.0000 550.0000 0;
node 213 100.0000 600.0000 0;
node 214 100.0000 650.0000 0;
node 215 100.0000 700.0000 0;
node 216 100.0000 750.0000 0;
node 217 100.0000 800.0000 0;
node 218 100.0000 850.0000 0;
node 219 100.0000 900.0000 0;
node 220 100.0000 950.0000 0;
node 221 100.0000 1000.0000 0;
node 301 150.0000 0.0000 0;
node 302 150.0000 50.0000 0;
node 303 150.0000 100.0000 0;
node 304 150.0000 150.0000 0;
node 305 150.0000 200.0000 0;
node 306 150.0000 250.0000 0;
node 307 150.0000 300.0000 0;
node 308 150.0000 350.0000 0;
node 309 150.0000 400.0000 0;
node 310 150.0000 450.0000 0;
node 311 150.0000 500.0000 0;
node 312 150.0000 550.0000 0;
node 313 150.0000 600.0000 0;
node 314 150.0000 650.0000 0;
node 315 150.0000 700.0000 0;
node 316 150.0000 750.0000 0;
node 317 150.0000 800.0000 0;
node 318 150.0000 850.0000 0;
node 319 150.0000 900.0000 0;
node 320 150.0000 950.0000 0;
node 321 150.0000 1000.0000 0;
node 401 200.0000 0.0000 0;
node 402 200.0000 50.0000 0;
node 403 200.0000 100.0000 0;
node 404 200.0000 150.0000 0;
node 405 200.0000 200.0000 0;
node 406 200.0000 250.0000 0;
node 407 200.0000 300.0000 0;
node 408 200.0000 350.0000 0;
node 409 200.0000 400.0000 0;
node 410 200.0000 450.0000 0;
node 411 200.0000 500.0000 0;
node 412 200.0000 550.0000 0;
node 413 200.0000 600.0000 0;
node 414 200.0000 650.0000 0;
node 415 200.0000 700.0000 0;
node 416 200.0000 750.0000 0;
node 417 200.0000 800.0000 0;
node 418 200.0000 850.0000 0;
node 419 200.0000 900.0000 0;
node 420 200.0000 950.0000 0;
node 421 200.0000 1000.0000 0;
node 501 250.0000 0.0000 0;
node 502 250.0000 50.0000 0;
node 503 250.0000 100.0000 0;
node 504 250.0000 150.0000 0;
node 505 250.0000 200.0000 0;
node 506 250.0000 250.0000 0;
node 507 250.0000 300.0000 0;
node 508 250.0000 350.0000 0;
node 509 250.0000 400.0000 0;
node 510 250.0000 450.0000 0;
node 511 250.0000 500.0000 0;
node 512 250.0000 550.0000 0;
node 513 250.0000 600.0000 0;
node 514 250.0000 650.0000 0;
node 515 250.0000 700.0000 0;
node 516 250.0000 750.0000 0;
node 517 250.0000 800.0000 0;
node 518 250.0000 850.0000 0;
node 519 250.0000 900.0000 0;
node 520 250.0000 950.0000 0;
node 521 250.0000 1000.0000 0;
node 601 300.0000 0.0000 0;
node 602 300.0000 50.0000 0;
node 603 300.0000 100.0000 0;
node 604 300.0000 150.0000 0;
node 605 300.0000 200.0000 0;
node 606 300.0000 250.0000 0;
node 607 300.0000 300.0000 0;
node 608 300.0000 350.0000 0;
node 609 300.0000 400.0000 0;
node 610 300.0000 450.0000 0;
node 611 300.0000 500.0000 0;
node 612 300.0000 550.0000 0;
node 613 300.0000 600.0000 0;
node 614 300.0000 650.0000 0;
node 615 300.0000 700.0000 0;
node 616 300.0000 750.0000 0;
node 617 300.0000 800.0000 0;
node 618 300.0000 850.0000 0;
node 619 300.0000 900.0000 0;
node 620 300.0000 950.0000 0;
node 621 300.0000 1000.0000 0;
node 701 350.0000 0.0000 0;
node 702 350.0000 50.0000 0;
node 703 350.0000 100.0000 0;
node 704 350.0000 150.0000 0;
node 705 350.0000 200.0000 0;
node 706 350.0000 250.0000 0;
node 707 350.0000 300.0000 0;
node 708 350.0000 350.0000 0;
node 709 350.0000 400.0000 0;
node 710 350.0000 450.0000 0;
node 711 350.0000 500.0000 0;
node 712 350.0000 550.0000 0;
node 713 350.0000 600.0000 0;
node 714 350.0000 650.0000 0;
node 715 350.0000 700.0000 0;
node 716 350.0000 750.0000 0;
node 717 350.0000 800.0000 0;
node 718 350.0000 850.0000 0;
node 719 350.0000 900.0000 0;
node 720 350.0000 950.0000 0;
node 721 350.0000 1000.0000 0;
node 801 400.0000 0.0000 0;
node 802 400.0000 50.0000 0;
node 803 400.0000 100.0000 0;
node 804 400.0000 150.0000 0;
node 805 400.0000 200.0000 0;
node 806 400.0000 250.0000 0;
node 807 400.0000 300.0000 0;
node 808 400.0000 350.0000 0;
node 809 400.0000 400.0000 0;
node 810 400.0000 450.0000 0;
node 811 400.0000 500.0000 0;
node 812 400.0000 550.0000 0;
node 813 400.0000 600.0000 0;
node 814 400.0000 650.0000 0;
node 815 400.0000 700.0000 0;
node 816 400.0000 750.0000 0;
node 817 400.0000 800.0000 0;
node 818 400.0000 850.0000 0;
node 819 400.0000 900.0000 0;
node 820 400.0000 950.0000 0;
node 821 400.0000 1000.0000 0;
node 901 450.0000 0.0000 0;
node 902 450.0000 50.0000 0;
node 903 450.0000 100.0000 0;
node 904 450.0000 150.0000 0;
node 905 450.0000 200.0000 0;
node 906 450.0000 250.0000 0;
node 907 450.0000 300.0000 0;
node 908 450.0000 350.0000 0;
node 909 450.0000 400.0000 0;
node 910 450.0000 450.0000 0;
node 911 450.0000 500.0000 0;
node 912 450.0000 550.0000 0;
node 913 450.0000 600.0000 0;
node 914 450.0000 650.0000 0;
node 915 450.0000 700.0000 0;
node 916 450.0000 750.0000 0;
node 917 450.0000 800.0000 0;
node 918 450.0000 850.0000 0;
node 919 450.0000 900.0000 0;
node 920 450.0000 950.0000 0;
node 921 450.0000 1000.0000 0;
node 1001 500.0000 0.0000 0;
node 1002 500.0000 50.0000 0;
node 1003 500.0000 100.0000 0;
node 1004 500.0000 150.0000 0;
node 1005 500.0000 200.0000 0;
node 1006 500.0000 250.0000 0;
node 1007 500.0000 300.0000 0;
node 1008 500.0000 350.0000 0;
node 1009 500.0000 400.0000 0;
node 1010 500.0000 450.0000 0;
node 1011 500.0000 500.0000 0;
node 1012 500.0000 550.0000 0;
node 1013 500.0000 600.0000 0;
node 1014 500.0000 650.0000 0;
node 1015 500.0000 700.0000 0;
node 1016 500.0000 750.0000 0;
node 1017 500.0000 800.0000 0;
node 1018 500.0000 850.0000 0;
node 1019 500.0000 900.0000 0;
node 1020 500.0000 950.0000 0;
node 1021 500.0000 1000.0000 0;
node 1101 550.0000 0.0000 0;
node 1102 550.0000 50.0000 0;
node 1103 550.0000 100.0000 0;
node 1104 550.0000 150.0000 0;
node 1105 550.0000 200.0000 0;
node 1106 550.0000 250.0000 0;
node 1107 550.0000 300.0000 0;
node 1108 550.0000 350.0000 0;
node 1109 550.0000 400.0000 0;
node 1110 550.0000 450.0000 0;
node 1111 550.0000 500.0000 0;
node 1112 550.0000 550.0000 0;
node 1113 550.0000 600.0000 0;
node 1114 550.0000 650.0000 0;
node 1115 550.0000 700.0000 0;
node 1116 550.0000 750.0000 0;
node 1117 550.0000 800.0000 0;
node 1118 550.0000 850.0000 0;
node 1119 550.0000 900.0000 0;
node 1120 550.0000 950.0000 0;
node 1121 550.0000 1000.0000 0;
node 1201 600.0000 0.0000 0;
node 1202 600.0000 50.0000 0;
node 1203 600.0000 100.0000 0;
node 1204 600.0000 150.0000 0;
node 1205 600.0000 200.0000 0;
node 1206 600.0000 250.0000 0;
node 1207 600.0000 300.0000 0;
node 1208 600.0000 350.0000 0;
node 1209 600.0000 400.0000 0;
node 1210 600.0000 450.0000 0;
node 1211 600.0000 500.0000 0;
node 1212 600.0000 550.0000 0;
node 1213 600.0000 600.0000 0;
node 1214 600.0000 650.0000 0;
node 1215 600.0000 700.0000 0;
node 1216 600.0000 750.0000 0;
node 1217 600.0000 800.0000 0;
node 1218 600.0000 850.0000 0;
node 1219 600.0000 900.0000 0;
node 1220 600.0000 950.0000 0;
node 1221 600.0000 1000.0000 0;
node 1301 650.0000 0.0000 0;
node 1302 650.0000 50.0000 0;
node 1303 650.0000 100.0000 0;
node 1304 650.0000 150.0000 0;
node 1305 650.0000 200.0000 0;
node 1306 650.0000 250.0000 0;
node 1307 650.0000 300.0000 0;
node 1308 650.0000 350.0000 0;
node 1309 650.0000 400.0000 0;
node 1310 650.0000 450.0000 0;
node 1311 650.0000 500.0000 0;
node 1312 650.0000 550.0000 0;
node 1313 650.0000 600.0000 0;
node 1314 650.0000 650.0000 0;
node 1315 650.0000 700.0000 0;
node 1316 650.0000 750.0000 0;
node 1317 650.0000 800.0000 0;
node 1318 650.0000 850.0000 0;
node 1319 650.0000 900.0000 0;
node 1320 650.0000 950.0000 0;
node 1321 650.0000 1000.0000 0;
node 1401 700.0000 0.0000 0;
node 1402 700.0000 50.0000 0;
node 1403 700.0000 100.0000 0;
node 1404 700.0000 150.0000 0;
node 1405 700.0000 200.0000 0;
node 1406 700.0000 250.0000 0;
node 1407 700.0000 300.0000 0;
node 1408 700.0000 350.0000 0;
node 1409 700.0000 400.0000 0;
node 1410 700.0000 450.0000 0;
node 1411 700.0000 500.0000 0;
node 1412 700.0000 550.0000 0;
node 1413 700.0000 600.0000 0;
node 1414 700.0000 650.0000 0;
node 1415 700.0000 700.0000 0;
node 1416 700.0000 750.0000 0;
node 1417 700.0000 800.0000 0;
node 1418 700.0000 850.0000 0;
node 1419 700.0000 900.0000 0;
node 1420 700.0000 950.0000 0;
node 1421 700.0000 1000.0000 0;
node 1501 750.0000 0.0000 0;
node 1502 750.0000 50.0000 0;
node 1503 750.0000 100.0000 0;
node 1504 750.0000 150.0000 0;
node 1505 750.0000 200.0000 0;
node 1506 750.0000 250.0000 0;
node 1507 750.0000 300.0000 0;
node 1508 750.0000 350.0000 0;
node 1509 750.0000 400.0000 0;
node 1510 750.0000 450.0000 0;
node 1511 750.0000 500.0000 0;
node 1512 750.0000 550.0000 0;
node 1513 750.0000 600.0000 0;
node 1514 750.0000 650.0000 0;
node 1515 750.0000 700.0000 0;
node 1516 750.0000 750.0000 0;
node 1517 750.0000 800.0000 0;
node 1518 750.0000 850.0000 0;
node 1519 750.0000 900.0000 0;
node 1520 750.0000 950.0000 0;
node 1521 750.0000 1000.0000 0;
node 1601 800.0000 0.0000 0;
node 1602 800.0000 50.0000 0;
node 1603 800.0000 100.0000 0;
node 1604 800.0000 150.0000 0;
node 1605 800.0000 200.0000 0;
node 1606 800.0000 250.0000 0;
node 1607 800.0000 300.0000 0;
node 1608 800.0000 350.0000 0;
node 1609 800.0000 400.0000 0;
node 1610 800.0000 450.0000 0;
node 1611 800.0000 500.0000 0;
node 1612 800.0000 550.0000 0;
node 1613 800.0000 600.0000 0;
node 1614 800.0000 650.0000 0;
node 1615 800.0000 700.0000 0;
node 1616 800.0000 750.0000 0;
node 1617 800.0000 800.0000 0;
node 1618 800.0000 850.0000 0;
node 1619 800.0000 900.0000 0;
node 1620 800.0000 950.0000 0;
node 1621 800.0000 1000.0000 0;
node 1701 850.0000 0.0000 0;
node 1702 850.0000 50.0000 0;
node 1703 850.0000 100.0000 0;
node 1704 850.0000 150.0000 0;
node 1705 850.0000 200.0000 0;
node 1706 850.0000 250.0000 0;
node 1707 850.0000 300.0000 0;
node 1708 850.0000 350.0000 0;
node 1709 850.0000 400.0000 0;
node 1710 850.0000 450.0000 0;
node 1711 850.0000 500.0000 0;
node 1712 850.0000 550.0000 0;
node 1713 850.0000 600.0000 0;
node 1714 850.0000 650.0000 0;
node 1715 850.0000 700.0000 0;
node 1716 850.0000 750.0000 0;
node 1717 850.0000 800.0000 0;
node 1718 850.0000 850.0000 0;
node 1719 850.0000 900.0000 0;
node 1720 850.0000 950.0000 0;
node 1721 850.0000 1000.0000 0;
node 1801 900.0000 0.0000 0;
node 1802 900.0000 50.0000 0;
node 1803 900.0000 100.0000 0;
node 1804 900.0000 150.0000 0;
node 1805 900.0000 200.0000 0;
node 1806 900.0000 250.0000 0;
node 1807 900.0000 300.0000 0;
node 1808 900.0000 350.0000 0;
node 1809 900.0000 400.0000 0;
node 1810 900.0000 450.0000 0;
node 1811 900.0000 500.0000 0;
node 1812 900.0000 550.0000 0;
node 1813 900.0000 600.0000 0;
node 1814 900.0000 650.0000 0;
node 1815 900.0000 700.0000 0;
node 1816 900.0000 750.0000 0;
node 1817 900.0000 800.0000 0;
node 1818 900.0000 850.0000 0;
node 1819 900.0000 900.0000 0;
node 1820 900.0000 950.0000 0;
node 1821 900.0000 1000.0000 0;
node 1901 950.0000 0.0000 0;
node 1902 950.0000 50.0000 0;
node 1903 950.0000 100.0000 0;
node 1904 950.0000 150.0000 0;
node 1905 950.0000 200.0000 0;
node 1906 950.0000 250.0000 0;
node 1907 950.0000 300.0000 0;
node 1908 950.0000 350.0000 0;
node 1909 950.0000 400.0000 0;
node 1910 950.0000 450.0000 0;
node 1911 950.0000 500.0000 0;
node 1912 950.0000 550.0000 0;
node 1913 950.0000 600.0000 0;
node 1914 950.0000 650.0000 0;
node 1915 950.0000 700.0000 0;
node 1916 950.0000 750.0000 0;
node 1917 950.0000 800.0000 0;
node 1918 950.0000 850.0000 0;
node 1919 950.0000 900.0000 0;
node 1920 950.0000 950.0000 0;
node 1921 950.0000 1000.0000 0;
node 2001 1000.0000 0.0000 0;
node 2002 1000.0000 50.0000 0;
node 2003 1000.0000 100.0000 0;
node 2004 1000.0000 150.0000 0;
node 2005 1000.0000 200.0000 0;
node 2006 1000.0000 250.0000 0;
node 2007 1000.0000 300.0000 0;
node 2008 1000.0000 350.0000 0;
node 2009 1000.0000 400.0000 0;
node 2010 1000.0000 450.0000 0;
node 2011 1000.0000 500.0000 0;
node 2012 1000.0000 550.0000 0;
node 2013 1000.0000 600.0000 0;
node 2014 1000.0000 650.0000 0;
node 2015 1000.0000 700.0000 0;
node 2016 1000.0000 750.0000 0;
node 2017 1000.0000 800.0000 0;
node 2018 1000.0000 850.0000 0;
node 2019 1000.0000 900.0000 0;
node 2020 1000.0000 950.0000 0;
node 2021 1000.0000 1000.0000 0;
#code from the OS wiki (http://opensees.berkeley.edu/wiki/index.php/ShellNL)
nDMaterial ElasticIsotropic 1 210000. 0.3;#steel
nDMaterial PlateFiber 2 1
section PlateFiber 3 2 4.;#4 mm thickness
#defining elements
element ShellNL 1 1 3 103 201 2 103 202 101 102 3
element ShellNL 3 3 5 105 203 4 105 204 103 104 3
element ShellNL 5 5 7 107 205 6 107 206 105 106 3
element ShellNL 7 7 9 109 207 8 109 208 107 108 3
element ShellNL 9 9 11 111 209 10 111 210 109 110 3
element ShellNL 11 11 13 113 211 12 113 212 111 112 3
element ShellNL 13 13 15 115 213 14 115 214 113 114 3
element ShellNL 15 15 17 117 215 16 117 216 115 116 3
element ShellNL 17 17 19 119 217 18 119 218 117 118 3
element ShellNL 19 19 21 121 219 20 121 220 119 120 3
element ShellNL 201 201 203 303 401 202 303 402 301 302 3
element ShellNL 203 203 205 305 403 204 305 404 303 304 3
element ShellNL 205 205 207 307 405 206 307 406 305 306 3
element ShellNL 207 207 209 309 407 208 309 408 307 308 3
element ShellNL 209 209 211 311 409 210 311 410 309 310 3
element ShellNL 211 211 213 313 411 212 313 412 311 312 3
element ShellNL 213 213 215 315 413 214 315 414 313 314 3
element ShellNL 215 215 217 317 415 216 317 416 315 316 3
element ShellNL 217 217 219 319 417 218 319 418 317 318 3
element ShellNL 219 219 221 321 419 220 321 420 319 320 3
element ShellNL 401 401 403 503 601 402 503 602 501 502 3
element ShellNL 403 403 405 505 603 404 505 604 503 504 3
element ShellNL 405 405 407 507 605 406 507 606 505 506 3
element ShellNL 407 407 409 509 607 408 509 608 507 508 3
element ShellNL 409 409 411 511 609 410 511 610 509 510 3
element ShellNL 411 411 413 513 611 412 513 612 511 512 3
element ShellNL 413 413 415 515 613 414 515 614 513 514 3
element ShellNL 415 415 417 517 615 416 517 616 515 516 3
element ShellNL 417 417 419 519 617 418 519 618 517 518 3
element ShellNL 419 419 421 521 619 420 521 620 519 520 3
element ShellNL 601 601 603 703 801 602 703 802 701 702 3
element ShellNL 603 603 605 705 803 604 705 804 703 704 3
element ShellNL 605 605 607 707 805 606 707 806 705 706 3
element ShellNL 607 607 609 709 807 608 709 808 707 708 3
element ShellNL 609 609 611 711 809 610 711 810 709 710 3
element ShellNL 611 611 613 713 811 612 713 812 711 712 3
element ShellNL 613 613 615 715 813 614 715 814 713 714 3
element ShellNL 615 615 617 717 815 616 717 816 715 716 3
element ShellNL 617 617 619 719 817 618 719 818 717 718 3
element ShellNL 619 619 621 721 819 620 721 820 719 720 3
element ShellNL 801 801 803 903 1001 802 903 1002 901 902 3
element ShellNL 803 803 805 905 1003 804 905 1004 903 904 3
element ShellNL 805 805 807 907 1005 806 907 1006 905 906 3
element ShellNL 807 807 809 909 1007 808 909 1008 907 908 3
element ShellNL 809 809 811 911 1009 810 911 1010 909 910 3
element ShellNL 811 811 813 913 1011 812 913 1012 911 912 3
element ShellNL 813 813 815 915 1013 814 915 1014 913 914 3
element ShellNL 815 815 817 917 1015 816 917 1016 915 916 3
element ShellNL 817 817 819 919 1017 818 919 1018 917 918 3
element ShellNL 819 819 821 921 1019 820 921 1020 919 920 3
element ShellNL 1001 1001 1003 1103 1201 1002 1103 1202 1101 1102 3
element ShellNL 1003 1003 1005 1105 1203 1004 1105 1204 1103 1104 3
element ShellNL 1005 1005 1007 1107 1205 1006 1107 1206 1105 1106 3
element ShellNL 1007 1007 1009 1109 1207 1008 1109 1208 1107 1108 3
element ShellNL 1009 1009 1011 1111 1209 1010 1111 1210 1109 1110 3
element ShellNL 1011 1011 1013 1113 1211 1012 1113 1212 1111 1112 3
element ShellNL 1013 1013 1015 1115 1213 1014 1115 1214 1113 1114 3
element ShellNL 1015 1015 1017 1117 1215 1016 1117 1216 1115 1116 3
element ShellNL 1017 1017 1019 1119 1217 1018 1119 1218 1117 1118 3
element ShellNL 1019 1019 1021 1121 1219 1020 1121 1220 1119 1120 3
element ShellNL 1201 1201 1203 1303 1401 1202 1303 1402 1301 1302 3
element ShellNL 1203 1203 1205 1305 1403 1204 1305 1404 1303 1304 3
element ShellNL 1205 1205 1207 1307 1405 1206 1307 1406 1305 1306 3
element ShellNL 1207 1207 1209 1309 1407 1208 1309 1408 1307 1308 3
element ShellNL 1209 1209 1211 1311 1409 1210 1311 1410 1309 1310 3
element ShellNL 1211 1211 1213 1313 1411 1212 1313 1412 1311 1312 3
element ShellNL 1213 1213 1215 1315 1413 1214 1315 1414 1313 1314 3
element ShellNL 1215 1215 1217 1317 1415 1216 1317 1416 1315 1316 3
element ShellNL 1217 1217 1219 1319 1417 1218 1319 1418 1317 1318 3
element ShellNL 1219 1219 1221 1321 1419 1220 1321 1420 1319 1320 3
element ShellNL 1401 1401 1403 1503 1601 1402 1503 1602 1501 1502 3
element ShellNL 1403 1403 1405 1505 1603 1404 1505 1604 1503 1504 3
element ShellNL 1405 1405 1407 1507 1605 1406 1507 1606 1505 1506 3
element ShellNL 1407 1407 1409 1509 1607 1408 1509 1608 1507 1508 3
element ShellNL 1409 1409 1411 1511 1609 1410 1511 1610 1509 1510 3
element ShellNL 1411 1411 1413 1513 1611 1412 1513 1612 1511 1512 3
element ShellNL 1413 1413 1415 1515 1613 1414 1515 1614 1513 1514 3
element ShellNL 1415 1415 1417 1517 1615 1416 1517 1616 1515 1516 3
element ShellNL 1417 1417 1419 1519 1617 1418 1519 1618 1517 1518 3
element ShellNL 1419 1419 1421 1521 1619 1420 1521 1620 1519 1520 3
element ShellNL 1601 1601 1603 1703 1801 1602 1703 1802 1701 1702 3
element ShellNL 1603 1603 1605 1705 1803 1604 1705 1804 1703 1704 3
element ShellNL 1605 1605 1607 1707 1805 1606 1707 1806 1705 1706 3
element ShellNL 1607 1607 1609 1709 1807 1608 1709 1808 1707 1708 3
element ShellNL 1609 1609 1611 1711 1809 1610 1711 1810 1709 1710 3
element ShellNL 1611 1611 1613 1713 1811 1612 1713 1812 1711 1712 3
element ShellNL 1613 1613 1615 1715 1813 1614 1715 1814 1713 1714 3
element ShellNL 1615 1615 1617 1717 1815 1616 1717 1816 1715 1716 3
element ShellNL 1617 1617 1619 1719 1817 1618 1719 1818 1717 1718 3
element ShellNL 1619 1619 1621 1721 1819 1620 1721 1820 1719 1720 3
element ShellNL 1801 1801 1803 1903 2001 1802 1903 2002 1901 1902 3
element ShellNL 1803 1803 1805 1905 2003 1804 1905 2004 1903 1904 3
element ShellNL 1805 1805 1807 1907 2005 1806 1907 2006 1905 1906 3
element ShellNL 1807 1807 1809 1909 2007 1808 1909 2008 1907 1908 3
element ShellNL 1809 1809 1811 1911 2009 1810 1911 2010 1909 1910 3
element ShellNL 1811 1811 1813 1913 2011 1812 1913 2012 1911 1912 3
element ShellNL 1813 1813 1815 1915 2013 1814 1915 2014 1913 1914 3
element ShellNL 1815 1815 1817 1917 2015 1816 1917 2016 1915 1916 3
element ShellNL 1817 1817 1819 1919 2017 1818 1919 2018 1917 1918 3
element ShellNL 1819 1819 1821 1921 2019 1820 1921 2020 1919 1920 3
#defining constraints
fix 1 1 1 1 0 0 0;
fix 101 1 1 1 0 0 0;
fix 201 1 1 1 0 0 0;
fix 301 1 1 1 0 0 0;
fix 401 1 1 1 0 0 0;
fix 501 1 1 1 0 0 0;
fix 601 1 1 1 0 0 0;
fix 701 1 1 1 0 0 0;
fix 801 1 1 1 0 0 0;
fix 901 1 1 1 0 0 0;
fix 1001 1 1 1 0 0 0;
fix 1101 1 1 1 0 0 0;
fix 1201 1 1 1 0 0 0;
fix 1301 1 1 1 0 0 0;
fix 1401 1 1 1 0 0 0;
fix 1501 1 1 1 0 0 0;
fix 1601 1 1 1 0 0 0;
fix 1701 1 1 1 0 0 0;
fix 1801 1 1 1 0 0 0;
fix 1901 1 1 1 0 0 0;
fix 2001 1 1 1 0 0 0;
fix 21 1 1 1 0 0 0;
fix 121 1 1 1 0 0 0;
fix 221 1 1 1 0 0 0;
fix 321 1 1 1 0 0 0;
fix 421 1 1 1 0 0 0;
fix 521 1 1 1 0 0 0;
fix 621 1 1 1 0 0 0;
fix 721 1 1 1 0 0 0;
fix 821 1 1 1 0 0 0;
fix 921 1 1 1 0 0 0;
fix 1021 1 1 1 0 0 0;
fix 1121 1 1 1 0 0 0;
fix 1221 1 1 1 0 0 0;
fix 1321 1 1 1 0 0 0;
fix 1421 1 1 1 0 0 0;
fix 1521 1 1 1 0 0 0;
fix 1621 1 1 1 0 0 0;
fix 1721 1 1 1 0 0 0;
fix 1821 1 1 1 0 0 0;
fix 1921 1 1 1 0 0 0;
fix 2021 1 1 1 0 0 0;
fix 2 1 1 1 0 0 0;
fix 3 1 1 1 0 0 0;
fix 4 1 1 1 0 0 0;
fix 5 1 1 1 0 0 0;
fix 6 1 1 1 0 0 0;
fix 7 1 1 1 0 0 0;
fix 8 1 1 1 0 0 0;
fix 9 1 1 1 0 0 0;
fix 10 1 1 1 0 0 0;
fix 11 1 1 1 0 0 0;
fix 12 1 1 1 0 0 0;
fix 13 1 1 1 0 0 0;
fix 14 1 1 1 0 0 0;
fix 15 1 1 1 0 0 0;
fix 16 1 1 1 0 0 0;
fix 17 1 1 1 0 0 0;
fix 18 1 1 1 0 0 0;
fix 19 1 1 1 0 0 0;
fix 20 1 1 1 0 0 0;
fix 2002 1 1 1 0 0 0;
fix 2003 1 1 1 0 0 0;
fix 2004 1 1 1 0 0 0;
fix 2005 1 1 1 0 0 0;
fix 2006 1 1 1 0 0 0;
fix 2007 1 1 1 0 0 0;
fix 2008 1 1 1 0 0 0;
fix 2009 1 1 1 0 0 0;
fix 2010 1 1 1 0 0 0;
fix 2011 1 1 1 0 0 0;
fix 2012 1 1 1 0 0 0;
fix 2013 1 1 1 0 0 0;
fix 2014 1 1 1 0 0 0;
fix 2015 1 1 1 0 0 0;
fix 2016 1 1 1 0 0 0;
fix 2017 1 1 1 0 0 0;
fix 2018 1 1 1 0 0 0;
fix 2019 1 1 1 0 0 0;
fix 2020 1 1 1 0 0 0;
#defining loads
pattern Plain 1 Linear {
load 1011 0 0 -1000 0 0 0;#N
}
recorder Node -file {E:\OpenSees\Out\shell_1.out} -time -nodeRange 1 2021 -dof 1 disp
recorder Node -file {E:\OpenSees\Out\shell_2.out} -time -nodeRange 1 2021 -dof 2 disp
recorder Node -file {E:\OpenSees\Out\shell_3.out} -time -nodeRange 1 2021 -dof 3 disp
recorder Node -file {E:\OpenSees\Out\shell_4.out} -time -nodeRange 1 2021 -dof 4 disp
recorder Node -file {E:\OpenSees\Out\shell_5.out} -time -nodeRange 1 2021 -dof 5 disp
recorder Node -file {E:\OpenSees\Out\shell_6.out} -time -nodeRange 1 2021 -dof 6 disp
#analysis settings
constraints Plain
numberer RCM
system UmfPack -lvalueFact 50
test NormDispIncr 1.0e-6 100
algorithm Newton
integrator Newmark 0.5 0.25
analysis Static
analyze 1
-
leopoldo.tesser
- Posts: 4
- Joined: Fri Mar 05, 2010 3:06 pm
- Location: University of Padua
- Contact:
Re: How to model a nonlinear shell element in OS
The ShellNL is a Nine-node Lagrangian element. It does not include geometrical nonlinearity. It can be used with nonlinear material.
I checked few lines of the code posted by jkbgbr and I found problems with node numbering that are sufficient to produce errors and need to be addressed.
As an example, from the definition of the nodes, the geometry of a little square is:
^ y
|
3----103----203
| | |
2----102----202
| | |
1 ---101----201----> x
but the first element is defined as:
element ShellNL 1 1 3 103 201 2 103 202 101 102 3
Two problmes:
1) a node is condensed (203->103) -> maybe you wanted to write 203 instead of the first 103
2) the node numbering should be couterclockwise: first corner nodes, then mid-side nodes and finally central node -> element ShellNL 1 1 201 203 3 101 202 103 2 102 3
Try to address all these problems and let us now.
Leopoldo Tesser
I checked few lines of the code posted by jkbgbr and I found problems with node numbering that are sufficient to produce errors and need to be addressed.
As an example, from the definition of the nodes, the geometry of a little square is:
^ y
|
3----103----203
| | |
2----102----202
| | |
1 ---101----201----> x
but the first element is defined as:
element ShellNL 1 1 3 103 201 2 103 202 101 102 3
Two problmes:
1) a node is condensed (203->103) -> maybe you wanted to write 203 instead of the first 103
2) the node numbering should be couterclockwise: first corner nodes, then mid-side nodes and finally central node -> element ShellNL 1 1 201 203 3 101 202 103 2 102 3
Try to address all these problems and let us now.
Leopoldo Tesser
Leopoldo
Re: How to model a nonlinear shell element in OS
Hello Leopoldo,
thank you for the fast answer!
You are right, node numbering was incorrect and this caused the problem. With the third node number corrected, my example yields results quite close to those of ansys (9.21 vs 9.43 mm, 4.29 vs 4.55 mm). The difference is likely due to discretization, i.e. there are 75% less ShellNL elements.
I'm sorry for this error!
I also checked the numbering rule. Although investigating this using square elements is not necessarily relevant for arbitrary cases, on my problem if a consistent direction is used for the corner and midside nodes (and n5 is the node next to n1 in this direction), the results are the same - as excepted. These node sequences resulted the same deflections:
element ShellNL 1 1 3 203 201 2 103 202 101 102 3 #clockwise
element ShellNL 1 1 201 203 3 101 202 103 2 102 3 #counterclockwise
element ShellNL 1 3 203 201 1 103 202 101 2 102 3 #clockwise, shifted first node
Thanks for your help again, and I'm sorry for my mistake.
If you don't mind, I'll send a PM on other questions regarding this element.
Gábor
thank you for the fast answer!
You are right, node numbering was incorrect and this caused the problem. With the third node number corrected, my example yields results quite close to those of ansys (9.21 vs 9.43 mm, 4.29 vs 4.55 mm). The difference is likely due to discretization, i.e. there are 75% less ShellNL elements.
I'm sorry for this error!
I also checked the numbering rule. Although investigating this using square elements is not necessarily relevant for arbitrary cases, on my problem if a consistent direction is used for the corner and midside nodes (and n5 is the node next to n1 in this direction), the results are the same - as excepted. These node sequences resulted the same deflections:
element ShellNL 1 1 3 203 201 2 103 202 101 102 3 #clockwise
element ShellNL 1 1 201 203 3 101 202 103 2 102 3 #counterclockwise
element ShellNL 1 3 203 201 1 103 202 101 2 102 3 #clockwise, shifted first node
Thanks for your help again, and I'm sorry for my mistake.
If you don't mind, I'll send a PM on other questions regarding this element.
Gábor
-
alessandrofascetti
- Posts: 2
- Joined: Mon Feb 20, 2012 3:13 pm
- Location: Sapienza
Re: How to model a nonlinear shell element in OS
Hello everybody.
Which version of OS is required to run this kind of elements? I am currently using 2.3.2.
The problem is OS cannot even create the elements. The error is
"ERROR: element ShellNL 1section not found
TclElementCommand -- unable to create element of type : ShellNL".
The example from OS Wiki will give me the same error. I also tried to run the script posted by jkbgbr, same error again.
I am attaching the script I'm using, if anybody knows what's the problem I would really appreciate.
# Units: mm/N/sec/MPa/tonne
#
# Previous Models Wipe-ing
#
wipe all
#
# 3D Problem => 6 DOF per node
#
model basicBuilder -ndm 3 -ndf 6
#
# Nodes Definition
#
node 11 0 0 0
node 12 100 0 0
node 13 200 0 0
node 14 300 0 0
node 15 400 0 0
node 21 0 100 0
node 22 100 100 0
node 23 200 100 0
node 24 300 100 0
node 25 400 100 0
node 31 0 200 0
node 32 100 200 0
node 33 200 200 0
node 34 300 200 0
node 35 400 200 0
node 41 0 300 0
node 42 100 300 0
node 43 200 300 0
node 44 300 300 0
node 45 400 300 0
node 51 0 400 0
node 52 100 400 0
node 53 200 400 0
node 54 300 400 0
node 55 400 400 0
#
# Boundary Conditions
#
fix 11 1 1 1 1 1 1
fix 12 1 1 1 1 1 1
fix 13 1 1 1 1 1 1
fix 14 1 1 1 1 1 1
fix 15 1 1 1 1 1 1
fix 21 1 1 1 1 1 1
fix 25 1 1 1 1 1 1
fix 31 1 1 1 1 1 1
fix 35 1 1 1 1 1 1
fix 41 1 1 1 1 1 1
fix 45 1 1 1 1 1 1
fix 51 1 1 1 1 1 1
fix 52 1 1 1 1 1 1
fix 53 1 1 1 1 1 1
fix 54 1 1 1 1 1 1
fix 55 1 1 1 1 1 1
#
# PlateFiber Material Definition
#
nDMaterial ElasticIsotropic 1 200000 0.3
nDMaterial PlateFiber 2 1
section PlateFiber 3 2 20
#
# element ShellNL $eleTag $node1 $node2 ... $node9 $secTag
#
# Non Linear Shell Elements Definition
#
element ShellNL 1 11 13 33 31 12 23 32 21 22 3
element ShellNL 2 13 15 35 33 14 25 34 23 24 3
element ShellNL 3 31 33 53 51 32 43 52 41 42 3
element ShellNL 4 33 35 55 53 34 45 54 43 44 3
#
# Point Load on the Middle-Node
#
pattern Plain 1 Constant {
load 33 0 0 -10000 0 0 0
}
#
# Analysis Settings
#
system BandGeneral
constraints Plain
numberer RCM
test NormDispIncr 1.0e-12 1000 3
algorithm Newton
integrator LoadControl 0.1
analysis Static
#
# Problem Solving
#
analyze 10
Which version of OS is required to run this kind of elements? I am currently using 2.3.2.
The problem is OS cannot even create the elements. The error is
"ERROR: element ShellNL 1section not found
TclElementCommand -- unable to create element of type : ShellNL".
The example from OS Wiki will give me the same error. I also tried to run the script posted by jkbgbr, same error again.
I am attaching the script I'm using, if anybody knows what's the problem I would really appreciate.
# Units: mm/N/sec/MPa/tonne
#
# Previous Models Wipe-ing
#
wipe all
#
# 3D Problem => 6 DOF per node
#
model basicBuilder -ndm 3 -ndf 6
#
# Nodes Definition
#
node 11 0 0 0
node 12 100 0 0
node 13 200 0 0
node 14 300 0 0
node 15 400 0 0
node 21 0 100 0
node 22 100 100 0
node 23 200 100 0
node 24 300 100 0
node 25 400 100 0
node 31 0 200 0
node 32 100 200 0
node 33 200 200 0
node 34 300 200 0
node 35 400 200 0
node 41 0 300 0
node 42 100 300 0
node 43 200 300 0
node 44 300 300 0
node 45 400 300 0
node 51 0 400 0
node 52 100 400 0
node 53 200 400 0
node 54 300 400 0
node 55 400 400 0
#
# Boundary Conditions
#
fix 11 1 1 1 1 1 1
fix 12 1 1 1 1 1 1
fix 13 1 1 1 1 1 1
fix 14 1 1 1 1 1 1
fix 15 1 1 1 1 1 1
fix 21 1 1 1 1 1 1
fix 25 1 1 1 1 1 1
fix 31 1 1 1 1 1 1
fix 35 1 1 1 1 1 1
fix 41 1 1 1 1 1 1
fix 45 1 1 1 1 1 1
fix 51 1 1 1 1 1 1
fix 52 1 1 1 1 1 1
fix 53 1 1 1 1 1 1
fix 54 1 1 1 1 1 1
fix 55 1 1 1 1 1 1
#
# PlateFiber Material Definition
#
nDMaterial ElasticIsotropic 1 200000 0.3
nDMaterial PlateFiber 2 1
section PlateFiber 3 2 20
#
# element ShellNL $eleTag $node1 $node2 ... $node9 $secTag
#
# Non Linear Shell Elements Definition
#
element ShellNL 1 11 13 33 31 12 23 32 21 22 3
element ShellNL 2 13 15 35 33 14 25 34 23 24 3
element ShellNL 3 31 33 53 51 32 43 52 41 42 3
element ShellNL 4 33 35 55 53 34 45 54 43 44 3
#
# Point Load on the Middle-Node
#
pattern Plain 1 Constant {
load 33 0 0 -10000 0 0 0
}
#
# Analysis Settings
#
system BandGeneral
constraints Plain
numberer RCM
test NormDispIncr 1.0e-12 1000 3
algorithm Newton
integrator LoadControl 0.1
analysis Static
#
# Problem Solving
#
analyze 10
Re: How to model a nonlinear shell element in OS
Hi,
You'll have to recompile OpenSees. Check my note in the opening topic about the error in the source code.
Cheers
Gábor
You'll have to recompile OpenSees. Check my note in the opening topic about the error in the source code.
Cheers
Gábor
-
alessandrofascetti
- Posts: 2
- Joined: Mon Feb 20, 2012 3:13 pm
- Location: Sapienza
Re: How to model a nonlinear shell element in OS
Thanks for the quick reply Gàbor.
Unfortunately, I don't seem to able to find the post you're talking about.
Can I contact you via PM to discuss the matter?
Thanks.
Alessandro
Unfortunately, I don't seem to able to find the post you're talking about.
Can I contact you via PM to discuss the matter?
Thanks.
Alessandro
Re: How to model a nonlinear shell element in OS
Hi,
sure you can.
"I must note that due to an error in the source of OS the code does not accept all arguments during object creation: int numData = 6; in line 65 should be int numData = 11; so the .exe must be re-compiled for this to work at all."
This is an error I found in the source code, the cpp file of the ShellNL element if I remember correctly. Due to the error only six parameters are used by the program, even if you provided all eleven. You should fix this and recompile OpenSees to make the element work.
Best,
Gábor
sure you can.
"I must note that due to an error in the source of OS the code does not accept all arguments during object creation: int numData = 6; in line 65 should be int numData = 11; so the .exe must be re-compiled for this to work at all."
This is an error I found in the source code, the cpp file of the ShellNL element if I remember correctly. Due to the error only six parameters are used by the program, even if you provided all eleven. You should fix this and recompile OpenSees to make the element work.
Best,
Gábor