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

船舶强度与结构设计大作业（一）船海1301 禹宗昕 U201312263数据与函数准备主要数据： 计算船长 L=200m 水密度 ρ=1t/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,:)代表第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=1:length(DD); Area(i)=spline(D,br(i,:),DD(i));end(2 ) function [ D ] = draft( df,da,n )% 作用：已知首尾吃水求各站吃水% df：参数，double，首吃水% da：参数，double，尾吃水% n:站数% D：输出，一维矩阵，n站对应吃水c=(df-da)/(n-1);for i=1:1:n D(i)=da+(i-1)*c;end;end（3）function [ V] = displacement( area,zhanju )% 已知各站面积求排水，核心：辛普森面积法% area：参数，一维矩阵，各站排水面积% zhanju：参数，double，站距% V：输出，double，排水体积S=0;for i=1:2:length(area)-2; S=S+(area(i)+4*area(i+1)+area(i+2))*2/6;V=S*zhanju;end;（4）function [ xb1] = flocenter( area,zhanju )% 作用：求浮心位置；核心：displacement函数拓展，141辛普森面积法% area：参数，一维矩阵，各站排水面积% zhanju：参数，一维矩阵，站距% xb1：输出，double，浮心到尾垂线的距离S=0;for i=1:2:length(area)-2; S=S+(area(i)+4*area(i+1)+area(i+2))*2/6;V=S*zhanju;end;%mb：各站体积对尾垂线的矩%下述两个for循环，是奇数站和偶数站对应的辛普森系数，第一站和最后一站的系数为0.5，其余为1，2交叉mb=area(1)*0.5*0+area(length(area))*0.5*(length(area)-1);for i=2:2:length(area)-1; mb=mb+area(i)*2*(i-1);end;for i=3:2:length(area)-2; mb=mb+area(i)*1*(i-1);end;xb1= mb*zhanju/V/3*2*zhanju;end浮态计算1.总重与重心计算 公式： 重心距离尾垂线xg1 == E23/C24*10 = 93.2136 m2.浮态计算（1）函数说明function [ d] = floatingstate(dm,L,A,R,xf,W,xg,zhanju)% 作用：通过迭代计算浮态，% dm：正浮时平均吃水，% A：正浮时水线面面积% R：纵稳心半径% xf：漂心% W：总重% xg：浮心距尾垂线距离% zhanju：站距% d：数据类型：类，% d.xg重心距舯m% d.xb浮心距舯m% d.weight总重t% d.displacement排水t% d.index1=index1=abs((W-V)/W);判断指标1，当误差小于0.5%时停止迭代% d.index2=index2=abs((xg-xb)/L);判断指标2，当误差小于0.1%时停止迭代% d.k 迭代次数k=0;index1=1;index2=1;%赋初始值，第一次取平均吃水，进行第一次排水与浮心计算d.df=dm; d.da=dm;D=draft(d.df,d.da,(L/zhanju+1));area=SAREA(D);V=displacement(area,zhanju);xb=flocenter(area,zhanju);%while循环迭代，当误差小于0.5%和0.1%时停止迭代，迭代步长