船舶强度与结构设计大作业剪力弯矩图.doc

船舶强度与结构设计大作业(一) 船海1301禹宗昕U201312263 :据与函数准备 主要数据： 计算船K: L=200m 水密度P =lt/m3 重力加速度g=9.8m/s2 原始资料(见附录) 全船重暈分布 全船邦戎曲线数据 静水平衡参数 总重量 W=122248 *9.8=120030.4 kN 水线面面积 A=4800m 平均吃水 dm=3.9m 纵稳心半径R=220m 漂心纵句华标xf=4.3m 基本函数简介 function [ Area ] = SAREA( DD ) %作用：已知各站吃水求各站对应的水下面积，核心方法：厂次样条曲线差值法 %D：己知水线 %br：原始邦戎数据，21*21矩阵,(因数据过大在下面已省略)br(i,0代表第i 行每站在各吃水下的面积， %DD：参数，一维矩阵，各站实际吃水 %Area：输出，-维矩阵，各站在DD吃水情况下对应的横截面积 D=[0,0.5,1,1.5,2,2.5,3,3.5,4,4.5,5,5.5,6,6.5,7,7.5,8,8.5,9,9.5,10,10.5,11,11.5,12;] br=[…]； for i=l:length(DD); Area(i)=spline(D,br(i,:),DD(i)); end (2 ) function [ D ] = draft( dfdan } %作用：已知首尾吃水求各站吃水 % df：参数，double,首吃水 % da：参数，double, M吃水 % n:站数 % D：输出，一维矩阵，n站对应吃水 c=(df-da)/(n-l); for i=l:l:n D(i)=da+(i-l)*c; end; end function [ V] = displacement! area,zhanju ) %已知各站面积求排水,核心：辛普森面积法 %area：参数，一维矩阵，各站排水面积 %zhanju：参数，double,站距 %V：输出，double,排水体积 S=0; for i=l:2:length(area)-2; S=S+(area(i)+4*area(i+l)+area(i+2))*2/6; V=S*zhanju; end; function [ xbl] = flocenter( areazhanju } %作用：求浮心位置;核心：displacement函数拓展，141辛普森面积法 %area：参数，一维矩阵，各站排水面积 %zhanju：参数，一维矩阵，站距 %xbl：输出，double,浮心到尾■线的距离 S=0; for i=l:2:length(area)-2; S=S+(area(i)+4*area(i+l)+area(i+2))*2/6; V=S*zhanju; end; %mb：各站体积对尾垂线的矩 %卜‘述两个for循环，是奇数站和偶数站对应的辛普淼系数，第，站和最后-?站 的系数为0.5,其余为1，2交叉 mb=area(l)*0.5*0+area(length(area))*0.5*(length(area)-l); for i=2:2:length(area)-l; mb=mb+area(i)*2*(i-l); end; for i=3:2:length(area)-2; mb=mb+area(i)*l*(i-l); end; xbl= mb*zhanju/V/3*2*zhanju; end 二.浮态计算 1.总重与重心计算 A B C D E 1 垂线间长200m 2 站号 重量(t) 力臂 对尾垂线的矩 3 0-1 791.9 0.5 395. 9735] 4 1-2 810.6 1.5 1215. 9] 5 2-3 855.8 2.5 2139. 39) 6 3-4 900.0 3.5 3150 4-5 399.6 4.5 1798. 3035 8 5-6 488. 1 5.5 2684. 385 9 6-7 525.6 6.5 3416. 7185] 10 7-8 410. 1 7.5 3076. 0425 11 8-9 558.9 8.5 4750. 837] 12 9-10 750. 1 9.5 7125.6175 13 10-11 710.5 10.5 7460. 251 14 11-12 670.8 11.5 7714. 2 15 12-13 540.9 12.5 6761. 25 16 13-14 650.0 13.5 8775 17 14-15 558.6 14.5 8099. 7 18 15-16 610.2 15.5 9458. 1 19 16-17 496.5 16.5 8191.821] 20 17-18 580.0 17.5 10150 21 18-19 520.4