% CDS Web r18 example: closed-end laminated pressure cylinder.
web_request = struct();
web_request.controls.laminate_progressive_failure = [ ...
    1;2;80;1;0.02;0.10;0.10;0.05;12;1;1;1;0.75;0;0;3;0];
web_request.controls.load_definition = [ ...
    1.5e9;1;25000;0;0;0;0;0;0;120;2.5e6;0.1e6;0.5;0.5;0;1];
web_request.geometry.controls = [ ...
    2;0.075;0.80;1;2;1;0.5;1;1;1;1;1;1;0.05;0.01; ...
    0.005;0.002;0.05;0.04;0.04;0.008;1;1];
web_request.controls.output = [0;0;0;0;0;0;0;0;0;0;1;0;10;8;0;1];

web_response = CDS_USolve_Web_v2026_09_05_r18(web_request);
cylinder_shell_summary_1D = web_response.outputs.O61(1:31);

web_request.geometry.controls(4) = 2;
web_response = CDS_USolve_Web_v2026_09_05_r18(web_request);
cylinder_hyer_summary_1D = web_response.outputs.O61(1:31);

fprintf('Shell hoop resultant = %.6g N/m\n',cylinder_shell_summary_1D(18));
fprintf('Hyer outer radial displacement = %.6g m\n',cylinder_hyer_summary_1D(29));
