%----------------------------------------------------------------- %%Composite Riser results in coloured bars (needs to be in hatched patterns) close all; clear all; %----------------------------------------------------------------- %%Composite Riser Local Design optimization, PhD Work, by Amaechi Chiemela Victor % Supervised by Prof Jianqiao Ye, Lancaster University, UK, Engineering % Department, July 2018 - local design was carried out in ANSYS ACP 19.0 %-------------------------------------Paper Details---------------------- % Composite risers for deep waters using a numerical modelling approach. % Authors: Chiemela Victor Amaechi, Nathaniel Gillett, Agbomerie Charles % Odijie, Xiaonan Hou, Jianqiao Ye. % Correspondence authors: j.ye2@lancaster.ac.uk; c.amaechi@lancaster.ac.uk. % For supplementary data, https://doi.org/10.1016/j.compstruct.2018.11.057. % Updated MATLAB code and published results in Composite Structures Journal % Published online on 22nd Nov 2018. Authors own rights to the results. %----------------------------------------------------------------------- % The authors acknowledge the support of Lancaster University, UK, Engineering % Department and Niger Delta Development Commission (NDDC) Nigeria. %------------------------------------------------------------------- %%Symbols used include H for Half in FibreHn89H, % n for and, as used in Fibre1n89 for hoop and axial layer angles %fibre454 means ratio of layers in axial:angled:hoop in fibre direction %--------------------------------------------------------------- % Create data for Stresses in Composite Riser under burst load fibre = [1074.36 1074.34 1074.31 1074.3 718.925 695.112 672.356 651.09 630.7 611.687 593.39 576.361 559.941 544.62 1394.57 1377.21 1362.75 1350.84]; transverse = [44.661 42.3581 40.1353 37.998 48.8363 46.6364 44.5233 42.5079 40.5641 38.7112 36.9172 35.2058 33.5478 31.9526 69.0091 61.0296 53.465 45.8869]; inplane = [0.0052427 0.00528007 0.00533829 0.0054269 98.6358 97.1526 95.7591 94.4302 93.18 91.9878 90.8643 89.7921 88.7832 87.8122 0.00583477 0.00594351 0.00602117 0.00610905]; % Data for effect of number of layers in Composite Riser under burst load fibre343 = [1410.9 1410.87 1410.86 857.935 831.33 805.786 782.16 759.366 738.379 718.012 699.38 1813.34 1793.26 1776.93]; fibre344 = [1285.1 1285.08 1285.07 730.4 705.493 681.646 659.429 638.071 618.225 599.09 581.269 1564.93 1545.05 1528.64 1515.28]; fibre443 = [1163.72 1163.7 1163.68 1163.66 896.358 870.396 845.554 822.512 800.383 779.926 760.193 742.066 1744.3 1724.19 1707.61]; fibre444 = [1066.42 1066.4 1066.38 1066.36 770.189 745.851 722.552 700.845 680.039 660.651 642.02 624.75 1512.45 1492.62 1475.93 1462.33]; fibre353 = [1405.5 1405.48 1405.47 788.86 763.065 738.257 715.294 693.113 672.65 652.79 634.531 616.722 600.464 1659.38 1642.1 1628.06]; fibre453 = [1170.23 1170.21 1170.18 1170.17 822.196 797.043 772.961 750.581 729.075 709.149 689.92 672.166 654.952 639.151 1590.68 1573.37 1559.15]; fibre354 = [1282.65 1282.62 1282.61 684.313 659.994 636.587 614.847 593.87 574.393 555.556 538.088 521.121 505.491 1447.79 1430.48 1416.08 1404.54]; fibre454 = [1074.36 1074.34 1074.31 1074.3 718.925 695.112 672.356 651.09 630.7 611.687 593.39 576.361 559.941 544.62 1394.57 1377.21 1362.75 1350.84]; transverse343 = [63.6588 61.143 58.7186 50.9676 48.2845 45.6933 43.2249 40.8285 38.5492 36.3238 34.2168 82.4745 71.9853 61.6072]; transverse344 = [46.6534 44.3282 42.0875 55.48 53.204 51.0178 48.9373 46.9292 45.0196 43.1734 41.4025 84.3036 75.0328 66.2517 57.4617]; transverse443 = [67.0873 64.5329 62.0709 59.6965 39.8621 37.1138 34.4452 31.9024 29.4302 27.0671 24.7679 22.5556 65.7358 55.6032 45.9908]; transverse444 = [50.5076 48.1392 45.8561 43.6549 101.015 99.4949 98.0603 96.6992 95.4158 94.1955 93.0428 91.9531 68.9849 60.147 51.7754 43.4016]; transverse353 = [55.1804 53.1842 50.8468 55.3026 52.8364 50.4592 48.1982 46.0078 43.9275 41.9015 39.98 38.0984 36.3204 82.5234 73.1033 63.8245]; transverse453 = [58.5376 56.0792 53.7094 51.4246 44.9339 42.4073 39.9717 37.6447 35.3914 33.2405 31.1478 29.1519 27.2001 25.3443 66.5862 57.6741 48.899]; transverse354 = [41.3069 39.0417 36.8588 58.4045 56.2745 54.2296 52.29 50.4159 48.6417 46.9228 45.2934 43.707 42.2152 83.4983 75.2508 67.0762 59.347]; transverse454 = [44.661 42.3581 40.1353 37.998 48.8363 46.6364 44.5233 42.5079 40.5641 38.7112 36.9172 35.2058 33.5478 31.9526 69.0091 61.0296 53.465 45.8869]; inplane343 = [0.00985932 0.0099647 0.0100844 124.803 123.125 121.556 120.058 118.657 117.319 116.064 114.87 0.01167 0.0117915 0.0119436]; inplane344 = [0.00673952 0.0068318 0.00691031 110.866 109.318 107.858 106.471 105.164 103.919 102.745 101.633 0.00804193 0.00812101 0.00819936 0.00830183]; inplane443 = [0.00785907 0.00792653 0.00802302 0.00816083 113.305 111.667 110.126 108.664 107.29 105.982 104.756 103.58 0.0085835 0.00874373 0.00882924]; inplane444 = [0.00514985 0.00519359 0.00529589 0.00533614 101.015 99.4949 98.0603 96.6992 95.4158 94.1955 93.0428 91.9531 0.00628575 0.0063602 0.00639891 0.00649462]; inplane353 = [0.0095341 0.00963707 0.00977844 120.708 119.093 117.573 116.128 114.771 113.475 112.259 111.096 110.007 108.956 0.0106585 0.010683 0.0109515]; inplane453 = [0.00743401 0.007530296 0.007629 0.00769204 109.588 108.01 106.523 105.113 103.785 102.52 101.332 100.197 99.133 98.1101 0.0084256 0.00853081 0.00865247]; inplane354 = [0.00704819 0.00712606 0.00723816 108.2 106.69 102.635 103.912 102.635 101.418 100.272 99.1755 98.1427 97.1596 0.00854438 0.00864319 0.00871584 0.00883448]; inplane454 = [0.0052427 0.00528007 0.00533829 0.0054269 98.6358 97.1526 95.7591 94.4302 93.18 91.9878 90.8643 89.7921 88.7832 87.8122 0.00583477 0.00594351 0.00602117 0.00610905]; % Data for effect of layer thickness Stresses in Composite Riser under burst load FibreControl = [1044.82 1044.8 1044.79 1044.76 697.269 674.181 652.074 631.444 611.659 593.191 575.455 558.908 542.965 528.177 1353.91 1337.12 1323.01 1311.54]; FibreAxial = [970.182 970.153 970.129 970.098 734.976 711.558 689.146 668.23 648.177 629.462 611.495 594.736 578.597 563.628 1366.93 1349.67 1335.1 1323.18]; FibreAngled = [526.755 526.735 526.715 526.697 785.244 785.852 733.34 710.34 687.953 667.889 648.214 630.705 613.383 598.123 1177.02 1160.7 1146.69 1134.94]; FibreHoop = [547.595 547.574 547.553 547.535 882.39 858.272 834.7 813.356 792.368 773.474 754.768 738.041 721.364 706.541 1356.95 1340.73 1326.67 1314.46]; TransverseControl = [43.3331 41.0982 38.944 36.867 47.4894 45.3549 43.306 41.3505 39.466 37.6682 35.929 34.2685 32.6614 31.1133 67.1221 59.3737 52.031 44.6715]; TransverseAxial = [46.2125 43.5325 40.9643 38.5021 48.9682 41.7419 39.6024 37.5575 35.5847 33.6995 31.8742 30.1284 28.437 26.8058 61.7055 53.95 46.5971 39.2337]; TransverseAngled = [42.6453 40.3644 38.166 36.046 27.6549 24.9202 22.2691 19.7786 17.3528 15.0726 12.8408 10.7424 8.67731 6.73955 31.0638 24.5989 18.2113 12.1881]; TransverseHoop = [53.9056 51.4984 49.1783 46.9409 22.2355 19.584 16.9894 14.5343 12.123 9.84038 7.58956 5.45946 3.34858 1.35656 29.0501 22.5708 16.0509 10.0255]; InPlaneControl = [0.004517 0.00453611 0.00462467 0.00466271 95.8207 94.3816 93.0289 91.7396 90.5259 89.3693 88.2786 87.2384 86.2589 85.3167 0.00501321 0.00509045 0.00517763 0.00523568]; InPlaneAxial = [0.00449467 0.00454919 0.00461388 0.00466596 93.7788 92.3161 90.9409 89.6304 88.3966 87.221 86.1124 85.0554 84.06 83.1031 0.00500971 0.00510474 0.00517309 0.00524832]; InPlaneAngled = [0.0159416 0.0161821 0.0163862 0.0165912 72.1387 70.4566 68.9 67.4223 66.0565 64.7556 63.555 62.4079 61.3507 60.3393 0.0193553 0.019543 0.0198069 0.0200066]; InPlaneHoop = [0.0162798 0.0164554 0.0166936 0.0168624 78.6339 77.0848 75.6413 74.2593 72.9733 71.7385 70.5914 69.4873 68.4626 67.4765 0.0193405 0.019521 0.0197367 0.0199266]; % Data for effect of axial layer angle Stresses in Composite Riser under burst load %Fibre0toH means Fibre0toHalf or Fibre0to0.5 Fibre0toH = [1074.31 1074.3 1074.28 1074.27 718.948 695.053 672.379 651.028 630.723 611.624 593.413 576.296 559.963 544.553 1394.64 1377.28 1362.81 1350.91]; Fibre0to1 = [1074.19 1074.19 1074.18 1074.17 718.918 694.976 672.346 650.936 630.702 611.517 593.397 576.188 559.921 544.52 1394.81 1377.51 1363.01 1351.1]; Fibre0to2 = [1073.68 1073.73 1073.77 1073.81 718.778 694.669 672.198 650.616 630.544 611.183 593.23 575.842 559.745 544.161 1395.61 1378.3 1363.76 1351.94]; Fibre0to3 = [1072.83 1072.96 1073.08 1073.19 718.489 694.211 671.892 650.137 630.221 610.684 592.891 575.322 559.391 543.621 1396.93 1379.62 1365.08 1353.27]; Fibre0to5 = [1070.06 1070.44 1070.81 1071.16 717.453 692.835 670.798 648.695 629.07 609.179 591.686 573.756 558.133 541.994 1401.17 1383.85 1369.31 1357.51]; Transverse0toH = [44.6515 42.3485 40.1284 37.9879 48.8462 46.6424 44.5333 42.514 40.5743 38.7173 36.9275 35.212 33.5583 31.9589 69.0215 61.0418 53.4767 45.8983]; Transverse0to1 = [44.6232 42.3194 40.0989 37.9573 48.8693 46.6678 44.5573 42.5376 40.6004 38.7401 36.9544 35.2363 33.5804 31.9966 69.0055 61.172 53.3857 46.0671]; Transverse0to2 = [44.5092 42.2029 39.9795 37.8355 48.9725 46.7579 44.6611 42.6306 40.7037 38.8353 37.0586 35.3314 33.6909 32.0797 69.2104 61.2252 53.6541 46.0701]; Transverse0to3 = [44.3186 42.008 39.7803 37.6317 49.1373 46.916 44.8279 42.7906 40.8725 38.9971 37.2295 35.4952 33.8639 32.2456 69.4626 61.4699 53.891 46.2995]; Transverse0to5 = [43.7041 41.38 39.1381 36.9752 49.6614 47.4283 45.3581 43.3088 41.4091 39.5215 37.7727 36.026 34.4139 32.7831 70.272 62.2546 54.651 47.0353]; InPlane0toH = [1.91715 1.88574 1.85547 1.82635 98.6423 97.1649 95.7656 94.4413 93.1877 91.9978 90.8726 89.8024 88.7889 87.83 0.00302205 0.00302364 0.00305626 0.00308837]; InPlane0to1 = [3.83991 3.77727 3.71684 3.65865 98.6691 97.1931 95.7922 94.4694 93.214 92.0259 90.8986 89.8306 88.8146 87.8582 0.000500784 0.000543843 0.000565496 0.000577173]; InPlane0to2 = [7.68856 7.56338 7.44273 7.32637 98.7789 97.3028 95.9013 94.5797 93.3213 92.1365 91.0048 89.9403 88.9228 87.9598 0.00823238 0.00832055 0.00839785 0.00850483]; InPlane0to3 = [11.5439 11.3563 11.1754 11.001 98.9607 97.4877 96.0821 94.7637 93.5011 92.3197 91.1836 90.1226 89.1007 88.1413 0.0151519 0.0153264 0.0154898 0.0156829]; InPlane0to5 = [19.288 18.9761 18.6755 18.3856 99.5442 98.08 96.6624 95.3519 94.0795 92.904 91.7595 90.7045 89.6711 88.7279 0.0278918 0.0282526 0.0286242 0.028966]; % Data for effect of hoop layer angle Stresses in Composite Riser under burst load %Fibre90to89H means Fibre90to89.5 Fibre90to89H = [1074.5 1074.48 1074.45 1074.44 718.796 695.36 672.22 651.34 630.557 611.94 593.242 576.617 559.788 544.877 1394.53 1377.17 1362.71 1350.8]; Fibre90to89 = [1074.89 1074.89 1074.88 1074.85 718.767 695.743 672.182 651.71 630.523 612.299 593.205 576.982 559.713 545.328 1394.4 1377.1 1362.55 1350.73]; Fibre90to88 = [1076.58 1076.56 1076.54 1076.52 719.097 696.812 672.479 652.767 630.792 613.347 593.447 578.026 559.929 546.372 1393.96 1376.66 1362.1 1350.27]; Fibre90to87 = [1079.35 1079.33 1079.31 1079.29 719.9 698.341 673.235 654.269 631.507 614.828 594.124 579.492 560.571 547.83 1393.23 1375.92 1361.35 1349.49]; Transverse90to89H = [44.6705 42.3675 40.1476 38.0072 48.8279 46.6461 44.5143 42.5174 40.5545 38.7206 36.9072 35.2149 33.5375 31.9612 69.0047 61.0248 53.4619 45.884]; Transverse90to89 = [44.6989 42.3956 40.1754 38.0347 48.8201 46.656 44.5052 42.5264 40.5437 38.7288 36.8952 35.2241 33.5185 31.9838 68.9333 61.1072 53.3241 46.0112]; Transverse90to88 = [44.8124 42.5078 40.2864 38.1445 48.8066 46.6765 44.4866 42.543 40.5208 38.7414 36.8678 35.2309 33.4933 31.9721 68.9299 60.9591 53.411 45.8433]; Transverse90to87 = [45.0018 42.6951 40.4716 38.3277 48.7956 46.6977 44.4675 42.558 40.496 38.7489 36.8333 35.2348 33.4443 31.986 68.7661 60.9715 53.2116 45.926]; InPlane90to89H = [0.0192227 0.0194428 0.019714 0.01993 98.6495 97.1612 95.7728 94.4375 93.1948 91.9939 90.8796 89.7986 88.7957 87.8263 1.60711 1.59163 1.5765 1.56546]; InPlane90to89 = [0.0331301 0.0335234 0.0339678 0.034357 98.6864 97.1893 95.8094 94.4651 93.2312 92.0211 90.9158 89.8255 88.8316 87.853 3.20825 3.17727 3.14692 3.12478]; InPlane90to88 = [0.0608411 0.0615796 0.0623691 0.0631035 98.8252 97.3101 95.9469 94.5843 93.3676 92.1388 91.0511 89.9421 88.9659 87.9688 6.41481 6.35288 6.29201 6.2477]; InPlane90to87 = [0.0883202 0.0894008 0.0905326 0.0916092 99.0508 97.5174 96.1702 94.7889 93.5888 92.3411 91.2705 90.1425 89.1835 88.1678 9.63118 9.53835 9.44689 9.38037]; % Data for effect of axial n hoop layer angles Stresses in Composite Riser under burst load %FibreHn89H means Fibre0.5n89.5 FibreHn89H = [1074.46 1074.44 1074.42 1074.41 718.818 695.301 672.242 651.279 630.58 611.877 593.264 576.552 559.811 544.81 1394.6 1377.24 1362.77 1350.86]; Fibre1n89 = [1074.75 1074.75 1074.75 1074.73 718.776 695.587 672.189 651.548 630.529 612.132 593.208 576.809 559.716 545.149 1394.67 1377.36 1362.82 1350.99]; Fibre2n88 = [1075.93 1075.97 1076.02 1076.05 718.97 696.352 672.341 652.287 630.643 612.849 593.287 577.509 559.76 545.837 1395.03 1377.72 1363.16 1351.33]; Fibre3n87 = [1077.88 1078 1078.13 1078.24 719.494 697.43 672.8 653.318 631.043 613.837 593.633 578.463 560.055 546.762 1395.64 1378.26 1363.79 1351.83]; TransverseHn89H = [44.6612 42.358 40.138 37.997 48.8349 46.6551 44.522 42.5245 40.5641 38.7266 36.9171 35.2226 33.542 31.9829 68.9618 61.1316 53.3464 46.0303]; Transverse1n89 = [44.6616 42.3574 40.1364 37.9947 48.8557 46.684 44.5412 42.5547 40.5806 38.7573 36.9326 35.251 33.5623 31.9967 69.0398 61.06 53.4987 45.9212]; Transverse2n88 = [44.6628 42.3548 40.13 37.9841 48.9392 46.7997 44.6216 42.6655 40.6592 38.8643 37.0078 35.3577 33.6282 32.1158 69.0712 61.2508 53.4701 46.1616]; Transverse3n87 = [44.6643 42.3498 40.1183 37.9663 49.0927 46.9782 44.7689 42.8387 40.8009 39.0328 37.144 35.522 33.7591 32.2765 69.2183 61.409 53.6357 46.3363]; InPlaneHn89H = [1.90338 1.87186 1.84142 1.81212 98.6576 97.1711 95.7808 94.4474 95.7808 92.0038 90.8876 89.8085 88.8036 87.8361 1.60372 1.58822 1.57302 1.56196]; InPlane1n89 = [3.81372 3.75072 3.68994 3.63138 98.7208 97.2272 95.8436 94.5029 93.2652 92.0588 90.9495 89.8631 88.8651 87.8904 3.20247 3.17139 3.14096 3.11873]; InPlane2n88 = [7.64609 7.52008 7.39865 7.28151 98.9672 97.459 96.088 94.7324 93.5078 92.2863 91.1905 90.0889 89.1045 88.1149 6.41223 6.34452 6.29209 6.23709]; InPlane3n87 = [11.5052 11.3162 11.134 10.9583 99.3745 97.8511 96.4919 95.1209 93.9087 92.6716 91.5886 90.4715 89.4999 88.4952 9.64524 9.55204 9.4603 9.39345]; % Data for effect of layer thickness Stresses in Composite Riser under burst load fibre454 = [1074.36 1074.34 1074.31 1074.3 718.925 695.112 672.356 651.09 630.7 611.687 593.39 576.361 559.941 544.62 1394.57 1377.21 1362.75 1350.84]; transverse454 = [44.661 42.3581 40.1353 37.998 48.8363 46.6364 44.5233 42.5079 40.5641 38.7112 36.9172 35.2058 33.5478 31.9526 69.0091 61.0296 53.465 45.8869]; inplane454 = [0.0052427 0.00528007 0.00533829 0.0054269 98.6358 97.1526 95.7591 94.4302 93.18 91.9878 90.8643 89.7921 88.7832 87.8122 0.00583477 0.00594351 0.00602117 0.00610905]; %%---------------------plots---------------------------------------------------- %%------------------------------------------------------------------------- %%------------------------------------------------------------------------- % Effect of number of plies in optimization of composite riser % Create a vertical bar chart using the bar function in fibre direction figure h=bar(1:18, [FibreControl' Fibre90to89H' Fibre90to89'], 1); % Add title and axis labels title('(a)') xlabel('Composite Riser Layers') ylabel('Stress in Fibre Direction (MPa)') % Add a legend legend('FibreHoop90', 'FibreHoop90to89.5', 'FibreHoop90to89') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) figure bar(1:18, [FibreControl' Fibre90to88' Fibre90to87'], 1) % Add title and axis labels title('(b)') xlabel('Composite Riser Layers') ylabel('Stress in Fibre Direction (MPa)') % Add a legend legend('FibreHoop90', 'FibreHoop90to88', 'FibreHoop90to87') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) % Create a vertical bar chart using the bar function in transverse direction figure bar(1:18, [TransverseControl' Transverse90to89H' Transverse90to89'], 1) % Add title and axis labels title('(c)') xlabel('Composite Riser Layers') ylabel('Stress in Transverse Direction (MPa)') % Add a legend legend('TransverseHoop90', 'TransverseHoop90to89.5', 'TransverseHoop90to89') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) figure bar(1:18, [TransverseControl' Transverse90to88' Transverse90to87'], 1) % Add title and axis labels title('(d)') xlabel('Composite Riser Layers') ylabel('Stress in Transverse Direction (MPa)') % Add a legend legend('TransverseHoop90', 'TransverseHoop90to88', 'TransverseHoop90to87') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) % Create a vertical bar chart using the bar function in inplane shear direction figure bar(1:18, [InPlaneControl' InPlane90to89H' InPlane90to89'], 1) % Add title and axis labels title('(e)') xlabel('Composite Riser Layers') ylabel('Stress in In-Plane Shear Direction (MPa)') % Add a legend legend('InPlaneHoop90', 'InPlaneHoop90to89.5', 'InPlaneHoop90to89') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) figure bar(1:18, [InPlaneControl' InPlane90to88' InPlane90to87'], 1) % Add title and axis labels title('(f)') xlabel('Composite Riser Layers') ylabel('Stress in In-Plane Shear Direction (MPa)') % Add a legend legend('InPlaneHoop90', 'InPlaneHoop90to88', 'InPlaneHoop90to87') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) %%------------------------------------------------------------------------- % Effect of number of plies in optimization of composite riser % Create a vertical bar chart using the bar function in fibre direction figure bar(1:18, [FibreControl' Fibre0toH' Fibre0to1'], 1) % Add title and axis labels title('(a)') xlabel('Composite Riser Layers') ylabel('Stress in Fibre Direction (MPa)') % Add a legend legend('FibreAxial0', 'FibreAxial0to0.5', 'FibreAxial0to1') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) figure bar(1:18, [Fibre0to2' Fibre0to3' Fibre0to5'], 1) % Add title and axis labels title('(b)') xlabel('Composite Riser Layers') ylabel('Stress in Fibre Direction (MPa)') % Add a legend legend('FibreAxial0to2', 'FibreAxial0to3', 'FibreAxial0to5') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) % Create a vertical bar chart using the bar function in transverse direction figure bar(1:18, [TransverseControl' Transverse0toH' Transverse0to1'], 1) % Add title and axis labels title('(c)') xlabel('Composite Riser Layers') ylabel('Stress in Transverse Direction (MPa)') % Add a legend legend('TransverseAxial0', 'TransverseAxial0to0.5', 'TransverseAxial0to1') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) figure bar(1:18, [Transverse0to2' Transverse0to3' Transverse0to5'], 1) % Add title and axis labels title('(d)') xlabel('Composite Riser Layers') ylabel('Stress in Transverse Direction (MPa)') % Add a legend legend('TransverseAxial0to2', 'TransverseAxial0to3', 'TransverseAxial0to5') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) % Create a vertical bar chart using the bar function in inplane shear direction figure bar(1:18, [InPlaneControl' InPlane0toH' InPlane0to1'], 1) % Add title and axis labels title('(e)') xlabel('Composite Riser Layers') ylabel('Stress in In-Plane Shear Direction (MPa)') % Add a legend legend('InPlaneAxial0', 'InPlaneAxial0to0.5', 'InPlaneAxial0to1') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) figure bar(1:18, [InPlane0to2' InPlane0to3' InPlane0to5'], 1) % Add title and axis labels title('(f)') xlabel('Composite Riser Layers') ylabel('Stress in In-Plane Shear Direction (MPa)') % Add a legend legend('InPlaneAxial0to2', 'InPlaneAxial0to3', 'InPlaneAxial0to5') % Set the axis limits set(gca, 'XTick', 1:18) set(gca, 'XLim', [0,19]) set(gca, 'Fontsize', 8) %-------------------------------------------- edit %----------------------------------------------------------------------- %%------------------------------------------------------------------------- %------------------------------------------------------------------------- % Effect of number of plies in optimization of composite riser % Create a vertical bar chart using the bar function in fibre direction figure bar(1:15, [fibre344' fibre443'], 1) % Add title and axis labels title('(a)') xlabel('Composite Riser Layers') ylabel('Stress in Fibre Direction (MPa)') % Add a legend legend('[(0)_3(+/-53.5)_4(90)_4]', '[(0)_4(+/-53.5)_3(90)_3]') % Set the axis limits set(gca, 'XTick', 1:15) set(gca, 'XLim', [0,16]) set(gca, 'Fontsize', 8) figure bar(1:16, [fibre444' fibre353'], 1) % Add title and axis labels title('(b)') xlabel('Composite Riser Layers') ylabel('Stress in Fibre Direction (MPa)') % Add a legend legend('[(0)_4(+/-53.5)_4(90)_4]', '[(0)_3(+/-53.5)_5(90)_3]') % Set the axis limits set(gca, 'XTick', 1:16) set(gca, 'XLim', [0,17]) set(gca, 'Fontsize', 8) figure bar(1:17, [fibre453' fibre354'], 1) % Add title and axis labels title('(c)') xlabel('Composite Riser Layers') ylabel('Stress in Fibre Direction (MPa)') % Add a legend legend('[(0)_4(+/-53.5)_5(90)_3]', '[(0)_3(+/-53.5)_5(90)_4]') % Set the axis limits set(gca, 'XTick', 1:17) set(gca, 'XLim', [0,18]) set(gca, 'Fontsize', 8) % Create a vertical bar chart using the bar function in transverse direction figure bar(1:15, [transverse344' transverse443'], 1) % Add title and axis labels title('(d)') xlabel('Composite Riser Layers') ylabel('Stress in Transverse Direction (MPa)') % Add a legend legend('[(0)_3(+/-53.5)_4(90)_4]', '[(0)_4(+/-53.5)_4(90)_3]') % Set the axis limits set(gca, 'XTick', 1:15) set(gca, 'XLim', [0,16]) set(gca, 'Fontsize', 8) figure bar(1:16, [transverse444' transverse353'], 1) % Add title and axis labels title('(e)') xlabel('Composite Riser Layers') ylabel('Stress in Transverse Direction (MPa)') % Add a legend legend('[(0)_4(+/-53.5)_4(90)_4]', '[(0)_3(+/-53.5)_5(90)_3]') % Set the axis limits set(gca, 'XTick', 1:16) set(gca, 'XLim', [0,17]) set(gca, 'Fontsize', 8) figure bar(1:17, [transverse453' transverse354'], 1) % Add title and axis labels title('(f)') xlabel('Composite Riser Layers') ylabel('Stress in Transverse Direction (MPa)') % Add a legend legend('[(0)_4(+/-53.5)_5(90)_3]', '[(0)_3(+/-53.5)_5(90)_4]') % Set the axis limits set(gca, 'XTick', 1:17) set(gca, 'XLim', [0,18]) set(gca, 'Fontsize', 8) % Create a vertical bar chart using the bar function in inplane shear direction figure bar(1:15, [inplane344' inplane443'], 1) % Add title and axis labels title('(g)') xlabel('Composite Riser Layers') ylabel('Stress in In-Plane Shear Direction (MPa)') % Add a legend legend('[(0)_3(+/-53.5)_4(90)_4]', '[(0)_4(+/-53.5)_4(90)_3]') % Set the axis limits set(gca, 'XTick', 1:15) set(gca, 'XLim', [0,16]) set(gca, 'Fontsize', 8) figure bar(1:16, [inplane444' inplane353'], 1) % Add title and axis labels title('(h)') xlabel('Composite Riser Layers') ylabel('Stress in In-Plane Shear Direction (MPa)') % Add a legend legend('[(0)_4(+/-53.5)_4(90)_4]', '[(0)_3(+/-53.5)_5(90)_3]') % Set the axis limits set(gca, 'XTick', 1:16) set(gca, 'XLim', [0,17]) set(gca, 'Fontsize', 8) figure bar(1:17, [inplane453' inplane354'], 1) % Add title and axis labels title('(i)') xlabel('Composite Riser Layers') ylabel('Stress in In-Plane Shear Direction (MPa)') % Add a legend legend('[(0)_4(+/-53.5)_5(90)_3]', '[(0)_3(+/-53.5)_5(90)_4]') % Set the axis limits set(gca, 'XTick', 1:17) set(gca, 'XLim', [0,18]) set(gca, 'Fontsize', 8) %------------------------------------------------------------------------- % Effect of angle orientation in optimization of composite riser % Create a vertical bar chart using the bar function in inplane shear direction %%%%unresolved part subplot(2,2,13) % add third plot in 2 x 2 grid bar(1:17, [fibre354'], 1) % Add title and axis labels title('Factor of Safety in Fibre Direction for case1') xlabel('Composite Riser Layers') ylabel('Factor of Safety in Fibre Direction') % Add a legend legend('Fibre Direction', 'Transverse Direction', 'In-Plane Shear') % Set the axis limits set(gca, 'XTick', 1:18) subplot(2,2,13) % add third plot in 2 x 2 grid bar(1:17, [fibre354'], 1) % Add title and axis labels title('Factor of Safety in Fibre Direction for case1') xlabel('Composite Riser Layers') ylabel('Factor of Safety in Fibre Direction') % Add a legend legend('Fibre Direction', 'Transverse Direction', 'In-Plane Shear') % Set the axis limits set(gca, 'XTick', 1:18)