Cyclic Analysis Failing early shell element

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
stokljos
Posts: 4
Joined: Fri Oct 09, 2020 7:16 pm

Cyclic Analysis Failing early shell element

Post by stokljos » Mon Nov 29, 2021 6:06 pm

Could someone take a look at my code to see why my analysis won't run past the second cycle? I changed from ShellMITC4 to ShellDKGQ. ShellMITC4 ran, but there was no strength loss like many people have also seen.
units are in kip/in.

(Had to shorten the node and shell definitions so it would fit)

wipe

#####################
puts "Initialization"
#####################

model BasicBuilder -ndm 3 -ndf 6;



#####################
puts "Define Nodes"
#####################


node 1 0.0 0 0
node 2 3.15 0 0
node 3 6.3 0 0
node 4 9.45 0 0
node 5 12.6 0 0
node 6 17.95 0 0
node 7 23.46 0 0
node 8 28.98 0 0
node 9 34.49 0 0
node 10 40.0 0 0
node 11 45.51 0 0
node 12 51.02 0 0
node 13 56.54 0 0
node 14 62.05 0 0
node 15 67.56 0 0
node 16 72.91 0 0
node 17 76.06 0 0
node 18 79.21 0 0
node 19 82.36 0 0
node 20 85.51 0 0

node 21 0.0 3.1 0
node 22 3.15 3.1 0
node 23 6.3 3.1 0
node 24 9.45 3.1 0
node 25 12.6 3.1 0
node 26 17.95 3.1 0
node 27 23.46 3.1 0
node 28 28.98 3.1 0
node 29 34.49 3.1 0
node 30 40.0 3.1 0
node 31 45.51 3.1 0
node 32 51.02 3.1 0
node 33 56.54 3.1 0
node 34 62.05 3.1 0
node 35 67.56 3.1 0
node 36 72.91 3.1 0
node 37 76.06 3.1 0
node 38 79.21 3.1 0
node 39 82.36 3.1 0
node 40 85.51 3.1 0

node 41 0.0 6.2 0
node 42 3.15 6.2 0
node 43 6.3 6.2 0
node 44 9.45 6.2 0
node 45 12.6 6.2 0
node 46 17.95 6.2 0
node 47 23.46 6.2 0
node 48 28.98 6.2 0
node 49 34.49 6.2 0
node 50 40.0 6.2 0
node 51 45.51 6.2 0
node 52 51.02 6.2 0
node 53 56.54 6.2 0
node 54 62.05 6.2 0
node 55 67.56 6.2 0
node 56 72.91 6.2 0
node 57 76.06 6.2 0
node 58 79.21 6.2 0
node 59 82.36 6.2 0
node 60 85.51 6.2 0

node 61 0.0 9.3 0
node 62 3.15 9.3 0
node 63 6.3 9.3 0
node 64 9.45 9.3 0
node 65 12.6 9.3 0
node 66 17.95 9.3 0
node 67 23.46 9.3 0
node 68 28.98 9.3 0
node 69 34.49 9.3 0
node 70 40.0 9.3 0
node 71 45.51 9.3 0
node 72 51.02 9.3 0
node 73 56.54 9.3 0
node 74 62.05 9.3 0
node 75 67.56 9.3 0
node 76 72.91 9.3 0
node 77 76.06 9.3 0
node 78 79.21 9.3 0
node 79 82.36 9.3 0
node 80 85.51 9.3 0

node 81 0.0 12.5 0
node 82 3.15 12.5 0
node 83 6.3 12.5 0
node 84 9.45 12.5 0
node 85 12.6 12.5 0
node 86 17.95 12.5 0
node 87 23.46 12.5 0
node 88 28.98 12.5 0
node 89 34.49 12.5 0
node 90 40.0 12.5 0
node 91 45.51 12.5 0
node 92 51.02 12.5 0
node 93 56.54 12.5 0
node 94 62.05 12.5 0
node 95 67.56 12.5 0
node 96 72.91 12.5 0
node 97 76.06 12.5 0
node 98 79.21 12.5 0
node 99 82.36 12.5 0
node 100 85.51 12.5 0
.
.
.
.

node 821 0.0 129.6 0
node 822 3.15 129.6 0
node 823 6.3 129.6 0
node 824 9.45 129.6 0
node 825 12.6 129.6 0
node 826 17.95 129.6 0
node 827 23.46 129.6 0
node 828 28.98 129.6 0
node 829 34.49 129.6 0
node 830 40.0 129.6 0
node 831 45.51 129.6 0
node 832 51.02 129.6 0
node 833 56.54 129.6 0
node 834 62.05 129.6 0
node 835 67.56 129.6 0
node 836 72.91 129.6 0
node 837 76.06 129.6 0
node 838 79.21 129.6 0
node 839 82.36 129.6 0
node 840 85.51 129.6 0

node 841 0.0 132.8 0
node 842 3.15 132.8 0
node 843 6.3 132.8 0
node 844 9.45 132.8 0
node 845 12.6 132.8 0
node 846 17.95 132.8 0
node 847 23.46 132.8 0
node 848 28.98 132.8 0
node 849 34.49 132.8 0
node 850 40.0 132.8 0
node 851 45.51 132.8 0
node 852 51.02 132.8 0
node 853 56.54 132.8 0
node 854 62.05 132.8 0
node 855 67.56 132.8 0
node 856 72.91 132.8 0
node 857 76.06 132.8 0
node 858 79.21 132.8 0
node 859 82.36 132.8 0
node 860 85.51 132.8 0

node 861 0.0 136.0 0
node 862 3.15 136.0 0
node 863 6.3 136.0 0
node 864 9.45 136.0 0
node 865 12.6 136.0 0
node 866 17.95 136.0 0
node 867 23.46 136.0 0
node 868 28.98 136.0 0
node 869 34.49 136.0 0
node 870 40.0 136.0 0
node 871 45.51 136.0 0
node 872 51.02 136.0 0
node 873 56.54 136.0 0
node 874 62.05 136.0 0
node 875 67.56 136.0 0
node 876 72.91 136.0 0
node 877 76.06 136.0 0
node 878 79.21 136.0 0
node 879 82.36 136.0 0
node 880 85.51 136.0 0




#####################
puts "Writing Materials"
#####################


#Confined Concrete
nDMaterial PlaneStressUserMaterial 1 40 7 6.165551810360332 0.2754699885451231 -1.2331103620720665 -0.0052 -0.046600043951363454 0.0010614863899717896 0.3

nDMaterial PlateFromPlaneStress 2 1 1962.7236683840022

#Unconfined Concrete
nDMaterial PlaneStressUserMaterial 3 40 7 4.74273216181564 0.2754699885451231 -0.0474273216181564 -0.0026 -0.03257800452931007 0.0008951670970672521 0.3

nDMaterial PlateFromPlaneStress 4 3 1962.7236683840022

#Horizontal Steel
uniaxialMaterial Steel02 5 73.53410416026082 29000 0.00826426712328766 18.5 0.925 0.15

#Hoop Steel
uniaxialMaterial Steel02 6 56.99980855026136 29000 0.008557326241134752 18.5 0.925 0.15

#Vertical web Steel
uniaxialMaterial Steel02 7 78.75546066868172 29000 0.00791620592592592 18.5 0.925 0.15

#Vertical web Steel
uniaxialMaterial Steel02 8 73.53410416026082 29000 0.0065169649315068405 18.5 0.925 0.15

nDMaterial ElasticIsotropic 9 3925.4473367680043 0.2

uniaxialMaterial MinMax 10 5 -min -0.03257800452931007 -max 0.18

uniaxialMaterial MinMax 11 6 -min -0.046600043951363454 -max 0.18

uniaxialMaterial MinMax 12 8 -min -0.03257800452931007 -max 0.143

uniaxialMaterial MinMax 13 7 -min -0.046600043951363454 -max 0.16

nDMaterial PlateRebar 14 10 0

nDMaterial PlateRebar 15 11 0



#####################
puts "Define Shells"
#####################


#Boundary region
section LayeredShell 1 8 4 0.984251968503937 15 0.4345 2 1.455321643424153 2 1.455321643424153 2 1.455321643424153 2 1.455321643424153 15 0.4345 4 0.984251968503937
#Web region
section LayeredShell 2 8 4 0.984251968503937 14 0.344 4 1.4640434570687344 4 1.4640434570687344 4 1.4640434570687344 4 1.4640434570687344 14 0.344 4 0.984251968503937
#Boundary region Elastic
section LayeredShell 3 8 9 0.984251968503937 15 0.017443627289163023 9 1.455321643424153 9 1.455321643424153 9 1.455321643424153 9 1.455321643424153 15 0.869 9 0.984251968503937
#Web region Elastic
section LayeredShell 4 8 9 0.984251968503937 14 0.344 9 1.4640434570687344 9 1.4640434570687344 9 1.4640434570687344 9 1.4640434570687344 14 0.344 9 0.984251968503937



#####################
puts "Define Elements"
#####################


element ShellDKGQ 1 1 2 22 21 1
element ShellDKGQ 2 2 3 23 22 1
element ShellDKGQ 3 3 4 24 23 1
element ShellDKGQ 4 4 5 25 24 1
element ShellDKGQ 5 5 6 26 25 2
element ShellDKGQ 6 6 7 27 26 2
element ShellDKGQ 7 7 8 28 27 2
element ShellDKGQ 8 8 9 29 28 2
element ShellDKGQ 9 9 10 30 29 2
element ShellDKGQ 10 10 11 31 30 2
element ShellDKGQ 11 11 12 32 31 2
element ShellDKGQ 12 12 13 33 32 2
element ShellDKGQ 13 13 14 34 33 2
element ShellDKGQ 14 14 15 35 34 2
element ShellDKGQ 15 15 16 36 35 2
element ShellDKGQ 16 16 17 37 36 1
element ShellDKGQ 17 17 18 38 37 1
element ShellDKGQ 18 18 19 39 38 1
element ShellDKGQ 19 19 20 40 39 1

element ShellDKGQ 20 21 22 42 41 1
element ShellDKGQ 21 22 23 43 42 1
element ShellDKGQ 22 23 24 44 43 1
element ShellDKGQ 23 24 25 45 44 1
element ShellDKGQ 24 25 26 46 45 2
element ShellDKGQ 25 26 27 47 46 2
element ShellDKGQ 26 27 28 48 47 2
element ShellDKGQ 27 28 29 49 48 2
element ShellDKGQ 28 29 30 50 49 2
element ShellDKGQ 29 30 31 51 50 2
element ShellDKGQ 30 31 32 52 51 2
element ShellDKGQ 31 32 33 53 52 2
element ShellDKGQ 32 33 34 54 53 2
element ShellDKGQ 33 34 35 55 54 2
element ShellDKGQ 34 35 36 56 55 2
element ShellDKGQ 35 36 37 57 56 1
element ShellDKGQ 36 37 38 58 57 1
element ShellDKGQ 37 38 39 59 58 1
element ShellDKGQ 38 39 40 60 59 1

element ShellDKGQ 39 41 42 62 61 1
element ShellDKGQ 40 42 43 63 62 1
element ShellDKGQ 41 43 44 64 63 1
element ShellDKGQ 42 44 45 65 64 1
element ShellDKGQ 43 45 46 66 65 2
element ShellDKGQ 44 46 47 67 66 2
element ShellDKGQ 45 47 48 68 67 2
element ShellDKGQ 46 48 49 69 68 2
element ShellDKGQ 47 49 50 70 69 2
element ShellDKGQ 48 50 51 71 70 2
element ShellDKGQ 49 51 52 72 71 2
element ShellDKGQ 50 52 53 73 72 2
element ShellDKGQ 51 53 54 74 73 2
element ShellDKGQ 52 54 55 75 74 2
element ShellDKGQ 53 55 56 76 75 2
element ShellDKGQ 54 56 57 77 76 1
element ShellDKGQ 55 57 58 78 77 1
element ShellDKGQ 56 58 59 79 78 1
element ShellDKGQ 57 59 60 80 79 1

element ShellDKGQ 58 61 62 82 81 1
element ShellDKGQ 59 62 63 83 82 1
element ShellDKGQ 60 63 64 84 83 1
element ShellDKGQ 61 64 65 85 84 1
element ShellDKGQ 62 65 66 86 85 2
element ShellDKGQ 63 66 67 87 86 2
element ShellDKGQ 64 67 68 88 87 2
element ShellDKGQ 65 68 69 89 88 2
element ShellDKGQ 66 69 70 90 89 2
element ShellDKGQ 67 70 71 91 90 2
element ShellDKGQ 68 71 72 92 91 2
element ShellDKGQ 69 72 73 93 92 2
element ShellDKGQ 70 73 74 94 93 2
element ShellDKGQ 71 74 75 95 94 2
element ShellDKGQ 72 75 76 96 95 2
element ShellDKGQ 73 76 77 97 96 1
element ShellDKGQ 74 77 78 98 97 1
element ShellDKGQ 75 78 79 99 98 1
element ShellDKGQ 76 79 80 100 99 1

element ShellDKGQ 77 81 82 102 101 1
element ShellDKGQ 78 82 83 103 102 1
element ShellDKGQ 79 83 84 104 103 1
element ShellDKGQ 80 84 85 105 104 1
element ShellDKGQ 81 85 86 106 105 2
element ShellDKGQ 82 86 87 107 106 2
element ShellDKGQ 83 87 88 108 107 2
element ShellDKGQ 84 88 89 109 108 2
element ShellDKGQ 85 89 90 110 109 2
element ShellDKGQ 86 90 91 111 110 2
element ShellDKGQ 87 91 92 112 111 2
element ShellDKGQ 88 92 93 113 112 2
element ShellDKGQ 89 93 94 114 113 2
element ShellDKGQ 90 94 95 115 114 2
element ShellDKGQ 91 95 96 116 115 2
element ShellDKGQ 92 96 97 117 116 1
element ShellDKGQ 93 97 98 118 117 1
element ShellDKGQ 94 98 99 119 118 1
element ShellDKGQ 95 99 100 120 119 1
.
.
.
.

element ShellDKGQ 704 741 742 762 761 1
element ShellDKGQ 705 742 743 763 762 1
element ShellDKGQ 706 743 744 764 763 1
element ShellDKGQ 707 744 745 765 764 1
element ShellDKGQ 708 745 746 766 765 2
element ShellDKGQ 709 746 747 767 766 2
element ShellDKGQ 710 747 748 768 767 2
element ShellDKGQ 711 748 749 769 768 2
element ShellDKGQ 712 749 750 770 769 2
element ShellDKGQ 713 750 751 771 770 2
element ShellDKGQ 714 751 752 772 771 2
element ShellDKGQ 715 752 753 773 772 2
element ShellDKGQ 716 753 754 774 773 2
element ShellDKGQ 717 754 755 775 774 2
element ShellDKGQ 718 755 756 776 775 2
element ShellDKGQ 719 756 757 777 776 1
element ShellDKGQ 720 757 758 778 777 1
element ShellDKGQ 721 758 759 779 778 1
element ShellDKGQ 722 759 760 780 779 1

element ShellDKGQ 723 761 762 782 781 1
element ShellDKGQ 724 762 763 783 782 1
element ShellDKGQ 725 763 764 784 783 1
element ShellDKGQ 726 764 765 785 784 1
element ShellDKGQ 727 765 766 786 785 2
element ShellDKGQ 728 766 767 787 786 2
element ShellDKGQ 729 767 768 788 787 2
element ShellDKGQ 730 768 769 789 788 2
element ShellDKGQ 731 769 770 790 789 2
element ShellDKGQ 732 770 771 791 790 2
element ShellDKGQ 733 771 772 792 791 2
element ShellDKGQ 734 772 773 793 792 2
element ShellDKGQ 735 773 774 794 793 2
element ShellDKGQ 736 774 775 795 794 2
element ShellDKGQ 737 775 776 796 795 2
element ShellDKGQ 738 776 777 797 796 1
element ShellDKGQ 739 777 778 798 797 1
element ShellDKGQ 740 778 779 799 798 1
element ShellDKGQ 741 779 780 800 799 1

element ShellDKGQ 742 781 782 802 801 3
element ShellDKGQ 743 782 783 803 802 3
element ShellDKGQ 744 783 784 804 803 3
element ShellDKGQ 745 784 785 805 804 3
element ShellDKGQ 746 785 786 806 805 4
element ShellDKGQ 747 786 787 807 806 4
element ShellDKGQ 748 787 788 808 807 4
element ShellDKGQ 749 788 789 809 808 4
element ShellDKGQ 750 789 790 810 809 4
element ShellDKGQ 751 790 791 811 810 4
element ShellDKGQ 752 791 792 812 811 4
element ShellDKGQ 753 792 793 813 812 4
element ShellDKGQ 754 793 794 814 813 4
element ShellDKGQ 755 794 795 815 814 4
element ShellDKGQ 756 795 796 816 815 4
element ShellDKGQ 757 796 797 817 816 3
element ShellDKGQ 758 797 798 818 817 3
element ShellDKGQ 759 798 799 819 818 3
element ShellDKGQ 760 799 800 820 819 3

element ShellDKGQ 761 801 802 822 821 3
element ShellDKGQ 762 802 803 823 822 3
element ShellDKGQ 763 803 804 824 823 3
element ShellDKGQ 764 804 805 825 824 3
element ShellDKGQ 765 805 806 826 825 4
element ShellDKGQ 766 806 807 827 826 4
element ShellDKGQ 767 807 808 828 827 4
element ShellDKGQ 768 808 809 829 828 4
element ShellDKGQ 769 809 810 830 829 4
element ShellDKGQ 770 810 811 831 830 4
element ShellDKGQ 771 811 812 832 831 4
element ShellDKGQ 772 812 813 833 832 4
element ShellDKGQ 773 813 814 834 833 4
element ShellDKGQ 774 814 815 835 834 4
element ShellDKGQ 775 815 816 836 835 4
element ShellDKGQ 776 816 817 837 836 3
element ShellDKGQ 777 817 818 838 837 3
element ShellDKGQ 778 818 819 839 838 3
element ShellDKGQ 779 819 820 840 839 3

element ShellDKGQ 780 821 822 842 841 3
element ShellDKGQ 781 822 823 843 842 3
element ShellDKGQ 782 823 824 844 843 3
element ShellDKGQ 783 824 825 845 844 3
element ShellDKGQ 784 825 826 846 845 4
element ShellDKGQ 785 826 827 847 846 4
element ShellDKGQ 786 827 828 848 847 4
element ShellDKGQ 787 828 829 849 848 4
element ShellDKGQ 788 829 830 850 849 4
element ShellDKGQ 789 830 831 851 850 4
element ShellDKGQ 790 831 832 852 851 4
element ShellDKGQ 791 832 833 853 852 4
element ShellDKGQ 792 833 834 854 853 4
element ShellDKGQ 793 834 835 855 854 4
element ShellDKGQ 794 835 836 856 855 4
element ShellDKGQ 795 836 837 857 856 3
element ShellDKGQ 796 837 838 858 857 3
element ShellDKGQ 797 838 839 859 858 3
element ShellDKGQ 798 839 840 860 859 3

element ShellDKGQ 799 841 842 862 861 3
element ShellDKGQ 800 842 843 863 862 3
element ShellDKGQ 801 843 844 864 863 3
element ShellDKGQ 802 844 845 865 864 3
element ShellDKGQ 803 845 846 866 865 4
element ShellDKGQ 804 846 847 867 866 4
element ShellDKGQ 805 847 848 868 867 4
element ShellDKGQ 806 848 849 869 868 4
element ShellDKGQ 807 849 850 870 869 4
element ShellDKGQ 808 850 851 871 870 4
element ShellDKGQ 809 851 852 872 871 4
element ShellDKGQ 810 852 853 873 872 4
element ShellDKGQ 811 853 854 874 873 4
element ShellDKGQ 812 854 855 875 874 4
element ShellDKGQ 813 855 856 876 875 4
element ShellDKGQ 814 856 857 877 876 3
element ShellDKGQ 815 857 858 878 877 3
element ShellDKGQ 816 858 859 879 878 3
element ShellDKGQ 817 859 860 880 879 3


element truss 818 1 21 0.6232932290586731 13
element truss 819 21 41 0.6232932290586731 13
element truss 820 41 61 0.6232932290586731 13
element truss 821 61 81 0.6232932290586731 13
element truss 822 81 101 0.6232932290586731 13
element truss 823 101 121 0.6232932290586731 13
element truss 824 121 141 0.6232932290586731 13
element truss 825 141 161 0.6232932290586731 13
element truss 826 161 181 0.6232932290586731 13
element truss 827 181 201 0.6232932290586731 13
element truss 828 201 221 0.6232932290586731 13
element truss 829 221 241 0.6232932290586731 13
element truss 830 241 261 0.6232932290586731 13
element truss 831 261 281 0.6232932290586731 13
element truss 832 281 301 0.6232932290586731 13
element truss 833 301 321 0.6232932290586731 13
element truss 834 321 341 0.6232932290586731 13
element truss 835 341 361 0.6232932290586731 13
element truss 836 361 381 0.6232932290586731 13
element truss 837 381 401 0.6232932290586731 13
element truss 838 401 421 0.6232932290586731 13
element truss 839 421 441 0.6232932290586731 13
element truss 840 441 461 0.6232932290586731 13
element truss 841 461 481 0.6232932290586731 13
element truss 842 481 501 0.6232932290586731 13
element truss 843 501 521 0.6232932290586731 13
element truss 844 521 541 0.6232932290586731 13
element truss 845 541 561 0.6232932290586731 13
element truss 846 561 581 0.6232932290586731 13
element truss 847 581 601 0.6232932290586731 13
element truss 848 601 621 0.6232932290586731 13
element truss 849 621 641 0.6232932290586731 13
element truss 850 641 661 0.6232932290586731 13
element truss 851 661 681 0.6232932290586731 13
element truss 852 681 701 0.6232932290586731 13
element truss 853 701 721 0.6232932290586731 13
element truss 854 721 741 0.6232932290586731 13
element truss 855 741 761 0.6232932290586731 13
element truss 856 761 781 0.6232932290586731 13
element truss 857 781 801 0.6232932290586731 13
element truss 858 801 821 0.6232932290586731 13
element truss 859 821 841 0.6232932290586731 13
element truss 860 841 861 0.6232932290586731 13

.
.
.

element truss 1033 6 26 0.24347391760104417 12
element truss 1034 26 46 0.24347391760104417 12
element truss 1035 46 66 0.24347391760104417 12
element truss 1036 66 86 0.24347391760104417 12
element truss 1037 86 106 0.24347391760104417 12
element truss 1038 106 126 0.24347391760104417 12
element truss 1039 126 146 0.24347391760104417 12
element truss 1040 146 166 0.24347391760104417 12
element truss 1041 166 186 0.24347391760104417 12
element truss 1042 186 206 0.24347391760104417 12
element truss 1043 206 226 0.24347391760104417 12
element truss 1044 226 246 0.24347391760104417 12
element truss 1045 246 266 0.24347391760104417 12
element truss 1046 266 286 0.24347391760104417 12
element truss 1047 286 306 0.24347391760104417 12
element truss 1048 306 326 0.24347391760104417 12
element truss 1049 326 346 0.24347391760104417 12
element truss 1050 346 366 0.24347391760104417 12
element truss 1051 366 386 0.24347391760104417 12
element truss 1052 386 406 0.24347391760104417 12
element truss 1053 406 426 0.24347391760104417 12
element truss 1054 426 446 0.24347391760104417 12
element truss 1055 446 466 0.24347391760104417 12
element truss 1056 466 486 0.24347391760104417 12
element truss 1057 486 506 0.24347391760104417 12
element truss 1058 506 526 0.24347391760104417 12
element truss 1059 526 546 0.24347391760104417 12
element truss 1060 546 566 0.24347391760104417 12
element truss 1061 566 586 0.24347391760104417 12
element truss 1062 586 606 0.24347391760104417 12
element truss 1063 606 626 0.24347391760104417 12
element truss 1064 626 646 0.24347391760104417 12
element truss 1065 646 666 0.24347391760104417 12
element truss 1066 666 686 0.24347391760104417 12
element truss 1067 686 706 0.24347391760104417 12
element truss 1068 706 726 0.24347391760104417 12
element truss 1069 726 746 0.24347391760104417 12
element truss 1070 746 766 0.24347391760104417 12
element truss 1071 766 786 0.24347391760104417 12
element truss 1072 786 806 0.24347391760104417 12
element truss 1073 806 826 0.24347391760104417 12
element truss 1074 826 846 0.24347391760104417 12
element truss 1075 846 866 0.24347391760104417 12

.
.
.

element truss 1635 20 40 0.6232932290586731 13
element truss 1636 40 60 0.6232932290586731 13
element truss 1637 60 80 0.6232932290586731 13
element truss 1638 80 100 0.6232932290586731 13
element truss 1639 100 120 0.6232932290586731 13
element truss 1640 120 140 0.6232932290586731 13
element truss 1641 140 160 0.6232932290586731 13
element truss 1642 160 180 0.6232932290586731 13
element truss 1643 180 200 0.6232932290586731 13
element truss 1644 200 220 0.6232932290586731 13
element truss 1645 220 240 0.6232932290586731 13
element truss 1646 240 260 0.6232932290586731 13
element truss 1647 260 280 0.6232932290586731 13
element truss 1648 280 300 0.6232932290586731 13
element truss 1649 300 320 0.6232932290586731 13
element truss 1650 320 340 0.6232932290586731 13
element truss 1651 340 360 0.6232932290586731 13
element truss 1652 360 380 0.6232932290586731 13
element truss 1653 380 400 0.6232932290586731 13
element truss 1654 400 420 0.6232932290586731 13
element truss 1655 420 440 0.6232932290586731 13
element truss 1656 440 460 0.6232932290586731 13
element truss 1657 460 480 0.6232932290586731 13
element truss 1658 480 500 0.6232932290586731 13
element truss 1659 500 520 0.6232932290586731 13
element truss 1660 520 540 0.6232932290586731 13
element truss 1661 540 560 0.6232932290586731 13
element truss 1662 560 580 0.6232932290586731 13
element truss 1663 580 600 0.6232932290586731 13
element truss 1664 600 620 0.6232932290586731 13
element truss 1665 620 640 0.6232932290586731 13
element truss 1666 640 660 0.6232932290586731 13
element truss 1667 660 680 0.6232932290586731 13
element truss 1668 680 700 0.6232932290586731 13
element truss 1669 700 720 0.6232932290586731 13
element truss 1670 720 740 0.6232932290586731 13
element truss 1671 740 760 0.6232932290586731 13
element truss 1672 760 780 0.6232932290586731 13
element truss 1673 780 800 0.6232932290586731 13
element truss 1674 800 820 0.6232932290586731 13
element truss 1675 820 840 0.6232932290586731 13
element truss 1676 840 860 0.6232932290586731 13
element truss 1677 860 880 0.6232932290586731 13




#####################
puts "Define Constraints"
#####################


fixY 0.0 1 1 1 1 1 1


#####################
puts "Define Gravity"
#####################


pattern Plain 1 Linear {
load 861 0 -15.174603711289922 0 0 0 0
load 862 0 -15.174603711289922 0 0 0 0
load 863 0 -15.174603711289922 0 0 0 0
load 864 0 -15.174603711289922 0 0 0 0
load 865 0 -15.174603711289922 0 0 0 0
load 866 0 -15.174603711289922 0 0 0 0
load 867 0 -15.174603711289922 0 0 0 0
load 868 0 -15.174603711289922 0 0 0 0
load 869 0 -15.174603711289922 0 0 0 0
load 870 0 -15.174603711289922 0 0 0 0
load 871 0 -15.174603711289922 0 0 0 0
load 872 0 -15.174603711289922 0 0 0 0
load 873 0 -15.174603711289922 0 0 0 0
load 874 0 -15.174603711289922 0 0 0 0
load 875 0 -15.174603711289922 0 0 0 0
load 876 0 -15.174603711289922 0 0 0 0
load 877 0 -15.174603711289922 0 0 0 0
load 878 0 -15.174603711289922 0 0 0 0
load 879 0 -15.174603711289922 0 0 0 0
load 880 0 -15.174603711289922 0 0 0 0
}


constraints Plain
numberer RCM
system BandGeneral
test NormDispIncr 1.0e-4 1000 ;
algorithm BFGS -count 500
integrator LoadControl 0.1;
analysis Static
analyze 10;

puts "gravity analysis ok..."
loadConst -time 0.0;





#####################
puts "Define Analysis"
#####################


pattern Plain 2 Linear {

load 861 1 0 0 0 0 0 -13
load 862 1 0 0 0 0 0 -13
load 863 1 0 0 0 0 0 -13
load 864 1 0 0 0 0 0 -13
load 865 1 0 0 0 0 0 -13
load 866 1 0 0 0 0 0 -13
load 867 1 0 0 0 0 0 -13
load 868 1 0 0 0 0 0 -13
load 869 1 0 0 0 0 0 -13
load 870 1 0 0 0 0 0 -13
load 871 1 0 0 0 0 0 -13
load 872 1 0 0 0 0 0 -13
load 873 1 0 0 0 0 0 -13
load 874 1 0 0 0 0 0 -13
load 875 1 0 0 0 0 0 -13
load 876 1 0 0 0 0 0 -13
load 877 1 0 0 0 0 0 -13
load 878 1 0 0 0 0 0 -13
load 879 1 0 0 0 0 0 -13
load 880 1 0 0 0 0 0 -13

}

set Peak [list 0.12 -0.12 0.12 -0.12 0.29 -0.29 0.29 -0.29 0.68 -0.64 0.66 -0.64 1.0 -0.95 1.0 -0.95 1.34 -1.29 1.3 -1.29 1.99 -1.94 1.97 -1.94 2.65 -2.59 2.67 -2.6 4.05 -3.75]
set nodeTag 870

set dofTag 1
set Dtot 0

constraints Plain;
numberer RCM
system BandGeneral
set globalSolutionTol 1.e-2;set numIter 50;set printFlagStatic 2;
test NormDispIncr $globalSolutionTol $numIter $printFlagStatic;
algorithm Newton;
analysis Static;
foreach ii $Peak {
set targetDisp [expr $ii];
if {abs($targetDisp) < 0.9} {set Dpos 0.01};if {abs($targetDisp) > 0.900} {set Dpos 0.01};
if {abs($targetDisp) < 0.9} {set Dneg -0.01};if {abs($targetDisp) > 0.90} {set Dneg -0.01};
if {abs($targetDisp) < 0.09} {set Dneg -0.0002};if {abs($targetDisp) < 0.09} {set Dpos 0.0002};
#set Dpos 0.01;set Dneg -0.01;
set ok 0;
if {$targetDisp > 0} {set Dincr $Dpos} else {set Dincr $Dneg}
if {$targetDisp == 0 && $Dtot > 0} {set Dincr $Dneg};
if {$targetDisp == 0 && $Dtot < 0} {set Dincr $Dpos};
set failTog 0;

while {$ok == 0} {
set tryTog 0;

integrator DisplacementControl $nodeTag $dofTag $Dincr
set ok [analyze 1]
if {$ok == 0} {
set stepSize $Dtot
puts "finished step $stepSize"
} else {puts "analysis failed, lets try smaller step size"
}

if {$ok != 0 && $failTog == 2} {
puts "analysis failed"
set failTog 1
break
if {$ok == 0 && $failTog == 2} {
set stepSize $Dtot
set failTog 0
puts "finished step $stepSize"
}
}

if {$ok != 0} {
set Dtry [expr $Dincr/10.]
set tryTog 1;
integrator DisplacementControl $nodeTag $dofTag $Dtry
set ok [analyze 1]

if {$ok == 0} {
set stepSize $Dtot
puts "finished step $stepSize"
} else {puts "analysis failed, lets try modified Newton"
}
}

if {$ok != 0} {
set tryTog 0;
integrator DisplacementControl $nodeTag $dofTag $Dincr
algorithm Newton -initial
set ok [analyze 1]

if {$ok == 0} {
set stepSize $Dtot
puts "finished step $stepSize"
} else {puts "analysis failed, lets try smaller step"
}
}

if {$ok != 0} {
set Dtry [expr $Dincr/10.]
set tryTog 1;
integrator DisplacementControl $nodeTag $dofTag $Dtry
set ok [analyze 1]

if {$ok == 0} {
set stepSize $Dtot
puts "finished step $stepSize"
} else {puts "analysis failed, try Krylov"
}
}

if {$ok != 0} {
set tryTog 0;
integrator DisplacementControl $nodeTag $dofTag $Dincr
algorithm KrylovNewton
set ok [analyze 1]

if {$ok == 0} {
set stepSize $Dtot
puts "finished step $stepSize"
} else {puts "analysis failed, lets try Krylov-Newton smaller"
}
}

if {$ok != 0} {
set Dtry [expr $Dincr/10.]
set tryTog 1;
integrator DisplacementControl $nodeTag $dofTag $Dtry
set ok [analyze 1]

if {$ok == 0} {
set stepSize $Dtot
puts "finished step $stepSize"
} else {puts "analysis failed"
}
}


if {$ok != 0 && $failTog < 2} {
set Dtry [expr $Dincr*100.]
set ok 0
set failTog 2
puts "analysis failed, try linear step"
}

#if {$ok != 0} {
#set failTog 1
#break
#}

algorithm Newton
incr stepCount
if {$tryTog > 0} {set delta $Dtry} else {set delta $Dincr}

set Dtot [expr $Dtot + $delta]

if {$Dincr > 0 && $Dtot > $ii} {break}
if {$Dincr < 0 && $Dtot < $ii} {break}
if {[expr abs($Dtot - $ii)] < 1e-4} {break}
}

if {$failTog > 0} {
puts "Wall Analysis Failed"
break
}

}
puts "Wall Analysis complete"

exit

Post Reply