资源描述：
《船舶强度与结构设计大作业-剪力弯矩图.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、船舶强度与结构设计大作业（一）船海1301禹宗昕U一．数据与函数准备1.主要数据：计算船长L=200m水密度ρ=1t/m³重力加速度g=9.8m/s²2.原始资料（见附录）（1）全船重量分布（2）全船邦戎曲线数据3.静水平衡参数总重量W=*9.8=.4kN水线面面积A=4800m平均吃水dm=3.9m纵稳心半径R=220m漂心纵向坐标xf=4.3m4.基本函数简介（1）function[Area]=SAREA(DD)%作用：已知各站吃水求各站对应的水下面积,核心方法：三次样条曲线差值法%D：已知水线%br：原始邦戎数据，21*21矩阵，(因数据过大在下面已省略)br(i,:)代表第i行每
2、站在各吃水下的面积，%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=[…..];fori=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：输出，一维矩阵
3、，n站对应吃水c=(df-da)/(n-1);fori=1:1:nD(i)=da+(i-1)*c;end;end（3）function[V]=displacement(area,zhanju)%已知各站面积求排水，核心：辛普森面积法%area：参数，一维矩阵，各站排水面积%zhanju：参数，double，站距%V：输出，double，排水体积S=0;fori=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)%作
4、用：求浮心位置；核心：displacement函数拓展，141辛普森面积法%area：参数，一维矩阵，各站排水面积%zhanju：参数，一维矩阵，站距%xb1：输出，double，浮心到尾垂线的距离S=0;fori=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(a
5、rea)-1);fori=2:2:length(area)-1;mb=mb+area(i)*2*(i-1);end;fori=3:2:length(area)-2;mb=mb+area(i)*1*(i-1);end;xb1=mb*zhanju/V/3*2*zhanju;end一．浮态计算1.总重与重心计算公式：重心距离尾垂线xg1=120Mi120Wi=E23/C24*10=93.2136m2.浮态计算（1）函数说明function[d]=floatingstate(dm,L,A,R,xf,W,xg,zhanju)%作用：通过迭代计算浮态，%dm：正浮时平均吃水，%A：正浮时水线面面积%
6、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/z
7、hanju+1));area=SAREA(D);V=displacement(area,zhanju);xb=flocenter(area,zhanju);%while循环迭代，当误差小于0.5%和0.1%时停止迭代，迭代步长：0.3while(index1>0.005)
8、
9、(index2>0.001)&&(k<20)d.df=d.df+((W-V)/(p*A)+(L/2-xf)*(xg-xb)/R)*0.3;d.da=d.da+(
