船舶强度与结构设计大作业-剪力弯矩图
船舶结构与强度设计 第2章
中和轴处腹板剪应力最大。 中和轴以上剖面对中和轴静矩：
S=[33.4×1.6×33.4/2+36×(33.4+1)]=892.4+1238.4=2130.8cm3
中和轴处腹板剪应力
Q S 380000 2130 .8 6612 N / cm 2 66 N / mm 2 I t 76532 .7 1.6
1. 规范波浪弯矩和剪力公式
中国船级社《钢质海船入级与建造规范》的波浪弯矩 和切力(剪力)计算公式，波浪弯矩和剪力与船长、船宽 方型系数有关。 （1）波浪弯矩 中拱波浪弯矩 中垂波浪弯矩
MW 190MCL2 BCb 103
kN.m
MW 110MCL2 B(Cb 0.7) 10 kN.m
M σ W
船体总强度理论很简单，关键要确定出作用在 船体梁上的弯矩和船体梁的剖面模数。
作用在船体梁上的弯矩：
M = Ms + Mw 式中 M——总弯矩；
Ms——静水弯矩，与重力及其船长分布有关
Mw——波浪弯矩，与波形范围内的船外形和波 浪要素有关。 船体剖面上的剪力： N = Ns + Nw Ns ——静水剪力
W
dw fs 2 A a a 1 10 6 2 A fs
如果不考虑面板和带板自身惯性矩65.3 cm3，则 C1=369463.2 cm4 I1=C1－A*E2=76467.4 cm4
对面板上表面剖面模数： W1=I1/E=2292 cm3 中国规范公式：
重量曲线举例 ①梯形法近似表示船体和舾装重量 已知总重量W, 重心据船中xg，假定梯形分布， 梯形面积等于W, 梯形面积形心纵坐标与重心一致， 并且b=1.195，解出a和c。
船舶结构强度设计第一章
Ns ( x ) Nw( x )
M x N ( x)dx N s ( x)dx N w ( x)dx M s M w
x x 0 0 0
Company Logo
x
Logo
1、重量曲线（Weight curve） W( x ) 2、浮力曲线（Buoyancy curve）
假定:1)纵倾调整过程中,水线面面积A不变；2）漂心位置 不变；3）不计船体变形的影响。
Company Logo
Logo
静水弯矩计算
Company L计算排水体积、纵向坐 标 xb ,漂心 x f ,、水线面面积A、纵稳心M半径、浮心高度 和浮心纵向位置。 检查浮心的纵向位置和重心的纵向位置，精度
(x) W( x ) ( x )
q( x ) 3、负荷曲线（Load curve）
4、剪力曲线（Shear curve） N( x)
x x x
q
0
x
(x)
dx
5、弯矩曲线（Bending moment curve）
示例
M(x)
N ( x )dx
0
q ( x )dx
分布在站外的重量处理 各站间重量叠加，得到各个站间的重量，如下图所示：
Company Logo
Logo
静水弯矩计算
浮力曲线 纵倾调整方法：解析法和逐步近似法。调整船舶在静水 中的纵倾时，用逐步近似法比较方便。 进行纵倾调整时，应具有邦戎曲线、静水力曲线以及船 舶重量重心等资料。
Company Logo
Logo
xg xb V0 V1 d f dm ( L / 2 xf ) R A
船舶强度与结构设计_多媒体课件_第一章
FBi − FAi
ε
ξi
FCi = FAi +
L
FBi − FAi
ε
(ξ 0 + xiψ )
⎫ ∫0 ⎪ ⎬ L ∫0 FC ( x) ⋅ xdx = V ⋅ xb ⎪ ⎭ FC ( x)dx = V
⎫ (ξ 0 + ψx)dx = V ∫0 FA ( x)dx + ∫0 ⎪ ⎪ ε ⎬ L L F ( x) − F ( x) x xb ⎪ x B A (ξ 0 + ψx)dx = V ∫0 FA ( x) L dx + ∫0 ε L L⎪ ⎭
△
ζ 平衡后波浪轴线位置 ζ ψ ζ
静水位置
△
△
Δb( x) = bω ( x) − bs ( x) = ρgΔF ( x)⎫ ⎬ ΔF ( x) = Fω ( x) − Fs ( x) ⎭
⎫ ⎪ o ⎬ x x x M ω ( x) = ∫ N ω ( x)dx = − ∫ ∫ Δb( x)dxdx⎪ 0 0 0 ⎭ N ω ( x) = − ∫ Δb( x)dx
1.2.1 浮态第1次近似计算
纵稳心半径R
'
ψ ψ
'
x g − xb ⎫ L d f1 = dm + ( − x f ) ⎪ R ⎪ 2 ⎬ x g − xb ⎪ L d a1 = d m − ( + x f ) R ⎪ 2 ⎭
⎫ W − Bi ≤ (0.1 ~ 0.5)% ⎪ W ⎪ ⎬ x g − xbi ≤ (0.05 ~ 0.1)%⎪ ⎪ L ⎭
对称的，零点在靠近船舯的某处，而在离 艏、艉端约船长的 1/4 处具有最大正值或负 值。
由于两端的剪力为零，即弯矩曲线在两端 的斜率为零，所以弯矩曲线在两端与纵坐 标轴相切。在计算过程中，常常利用这些 性质来检查计算结果是否正确。 修正：
货运05 船舶强度
三、船舶总体布置对总纵弯曲变形的影响
2.尾机船 尾机船的机舱位于尾部。空航压载时, 有较大的中拱变形状 态。而满载时的拱垂变形因船舶规模的不同而异,大型船 满载时呈中垂变形,而普通船可能处于中垂或中拱变形状 态。尾机船最严重的中拱变形发生在空航压载状态。 措施： (1)合理安排压载水:尾机船空船压载航行时压载水的数量与 位置应能同时达到减少尾倾、增加吃水和减缓中拱变形三 个目的。为此,应尽量使用漂心以前靠近中部的压载水舱, 而不能单独使用首部压载水舱,以免增大中拱变形。杂货 船中部若有深舱,可作为压载水舱使用。油轮属于尾机船, 作为空船压载使用的清洁压载舱应布置在中部略靠前的位 置。 (2)油水分配及使用:油水分配与使用的原则与中机船满载时 相同,即分配时自中部向首尾依次装载,使用时顺序相反。
(kPa)
Hc—甲板设计堆高，重结构取1.5m, 轻结构取1.2m。 μ — 设计舱容系数。
三、用经验方法确定的允许负荷
2．中间甲板和底舱：
允许负荷：
Pa
gHd
(kPa)
H d — 舱高。 无设计值时，取rc=0.72 t/m3, 重结构取rc=1.2 t/m3。 rc =1/μ
§11-2 保证满足船舶的局部强度条件
§11-1 船舶纵向强度
3.船中静水弯矩估算法 中小型船可只算船中弯矩
（1）实际静水弯矩
1 M S ' g ( WLi X i Pi X i ) Bi X i 2
(kN m)
特定船
M S ' f ( P (kN m) i X i , dm )
7．船舶总纵弯曲变形的判断 拱垂值：δ=|dж- dm| (m)
《船舶结构力学》0预备知识-弯矩图、剪力图复习
x
FS x＝qx
0 x l
FS x
ql
M x＝qx2 / 2 0 x l
依方程画出剪力图和弯矩图
x
ql2 / 2 由剪力图、弯矩图可见。最
大剪力和弯矩分别为
ql 2 / 8
FS max＝ql M max＝ql 2 / 2
x
目录
(Internal Forces in Beams)
例题5-3
F
a
b
A
C
x1 x2
FAY
l
FS Fb / l
图示简支梁C点受集中力作用。
B 试写出剪力和弯矩方程，并画 出剪力图和弯矩图。
FBY 解：1．确定约束力
M A＝0， MB＝0
FAy＝Fb/l FBy＝Fa/l
2．写出剪力和弯矩方程
Fa / l
Fab/ l
M
x AC
FS x1＝Fb / l 0 x1 a
M x1＝Fbx1 / l 0 x1 a
CB FS x2 ＝ Fa / l a x2 l
M x2 ＝Fal x2 / l a x2 l
x 3. 依方程画出剪力图和弯矩图。
目录
(Internal Forces in Beams)
例题5-4
a
b
图示简支梁C点受集中力偶作用。
M
A
C
x1
试写出剪力和弯矩方程，并画
M x2 ＝ Mx2 / l 0 x2 b
3. 依方程画出剪力图和弯矩图。
目录
(Internal Forces in Beams)
二、q(x)、Fs(x)图、 M（x）图三者间的关系 (relationships between load,shear force,and bending moment diagrams)
第2章 船舶静置在波浪中的剪力和弯矩计算
根据表格7.2静水平衡计算
1
理论站 号
2
力臂系 数k
3
横剖面 浸水面 积ω/㎡
4
5
6
第一次近似计算
第二次近似计算
力矩函 横剖面 力矩函 数 浸水面 数 k×ω= 积ω/㎡ k×ω= （2）× （2）× （3） （5 ）
Bij
2012年9月
i j
2
L g
船舶工程系 孟巧
第二章 船舶静置在波浪中的剪力和弯矩计算
2.1.2 载荷、剪力、弯矩的基本公式 和计算的一般步骤
基本公式
bw( x) bs ( x) b( x)
q( x) w( x) - bw( x)
N ( x) q( x)dx N s ( x) N w ( x)
0
x
x
M ( x) N ( x)dx M s ( x) M w ( x)
A1 A2 A3 W
A y A
i i
i
L - xg 2
化简后得到：
a c 4b 6
108 x g a-c 7 L
2012年9月
第二章 船舶静置在波浪中的剪力和弯矩计算
船舶工程系 孟巧
根据统计资料 对于瘦形船：b=1.195，于是：
54 xg a 0.61 7 L 54 xg c 0.61- 7 L
0
作用在船体梁上的所有外力是平衡的
2012年9月 第二章 船舶静置在波浪中的剪力和弯矩计算 船舶工程系 孟巧
2.2.4 静水剪力、弯矩曲线
定义：船体梁在静水中所受到的剪力和弯 矩沿船长分布状况的曲线分别称为静水剪 力曲线和静水弯矩曲线。 表达式为
第8章船舶总纵强度ppt课件
3）扭转强度：
– 船舶抵抗扭转变形和破坏的能力。
4）局部强度：
– 船舶抵抗局部变形和破坏的能力。
3）扭转强度（Torsional strength）
定义
– 船舶结构抵抗船体沿船长方向发生扭转变形的能力。
产生原因
– 沿船长方向单位长度重力和浮力横向、纵向不共垂线 造成的。
2676
268 2944 2408
3636
364 4000 3272
2522
252 2774 2270
13408
b、隔舱装载： 应用：合理分配载重提高总纵强度
合理分配载重，提高船舶总纵强度
方法：
货物配置：
➢ 按舱容比分配货物，在舱容允许的条件下，中区货舱应按装货 重量的上限值装，首尾货舱按下限值装；中途港货物不应过分 集中于中区货舱。
位于距首尾L/4的剖面处； 位于船中剖面处。
剪应力最大值τmax： 弯曲应力最大值σmax ：
位于距首尾L/4剖面的中和轴处； 位于船中剖面的上甲板处。
0.5t/m
O
2.5t
12.5t.m
B A
X
C
D
-2.5t
二、 影响船舶总纵强度的因素
1、重力分布对拱垂变形的影响：
– （1）船舶布置 – （2）积载方案
X
危险。
b、船舶中拱，
W
波 峰位于船中，
X
危险
三、标准计算状态
1. 标准波
a. 波长等于船长； b. 波峰（谷）位于船中； c. 坦谷波；
d. h/ 1/ 20
2. 标准装载状态
a. 满载出港 b. 满载到港（剩10％燃料） c. 压载航行
总纵向强度校核校核原理1许用剪力和许用弯矩
三、船体总布置对总弯曲的影响
中机船
满载---中拱 空载---中垂
尾机船
满载---中垂 空载---中拱
中后机船
满载---中垂（中拱） 空载---中拱
四、保证总纵强度满足要求措施
按舱容比分配各舱货物重量；
Pi
Vi Vch
Q
（上下浮动调整值一般取分配量的10%） 根据机舱不同位置适当调整中区货舱货物分
40 36
有利范围
+ |MA| (不包括空船重量) 32 允许范围
图中各曲线的意义：
28 24
MS′=0得点划线
20 中拱 16
中垂
MS′=±MSO得上下两虚线 12
8
MS′=± MS得最外侧两实线 4
3.0 4.0 5.0 6.0 7.0 8.0 9.0 dm (m)
3.站面强度校核法 ⑴许用剪力和许用弯矩 ⑵实际剪力和实际弯矩
P
104.97
下层垫木应沿横向跨甲板下纵骨的方法设置，其跨度 为4倍的纵骨间距。
例题2：某船No.3二层舱高3.5m，舱 容1500m3，均布载荷时的许用负荷量 Pd=24.96kPa，装载如图所示，装货后货 物表面呈水平，试校核二层甲板的局部强度， 若不满足，应如何调整，调整后的局部强度 如何？
相应站面处的剪力为自船尾始到该站面处的重力 和浮力之差。
实际剪力：SFi=Wi-Bi Wi---站面处重力 Bi ---站面处浮力
实际弯矩：BMi=Mi-Wi·Li-Mbi Mi ---站面处重力矩 Mbi ---站面处浮力矩 Li --- 计算点到船中距离 当采用船尾坐标系计算重力矩时，相应站面处
3.船体剪切变形与弯曲变形 ★中拱：
船体受正弯矩作用，中部上拱，中部浮力大 于重力，首、尾部浮力小于重力，上甲板受拉，船 底受压。 ★中垂：
华科船舶结构强度第二次大作业
船体强度与结构设计------第二次大作业班级：姓名：学号：题目：图示为某船舶横剖面结构示意图。
请计算当船舶船舯为波谷，且弯矩值为9.0×107N ·m ，考虑折减系数计算总纵弯矩应力。
解答： 一、计算依据 1、计算载荷计算弯矩 79.010m M N =⨯⋅2、船体材料计算剖面的所有构件均采用低碳钢，屈服极限=235aY MP σ3、许用应力（1）总纵弯曲许用应力[]0.5Yσσ=(2)总纵弯曲与板架局部玩去合成应力的许用应力：在板架跨中 12[+]0.65Y σσσ= 在横仓壁处 12[+]Y σσσ=二、总纵弯曲正应力计算1、总纵弯曲正应力第一次近似计算肋骨剖面计算简图如题图所示。
将图中个强力构件编号并将其尺寸填入表中。
船体剖面要素及第一次近似总纵弯曲应力的计算在下表中完成。
在计算中，参考轴取在基线处。
利用上表中的数据可得第一次近似中和轴距参考轴的距离为：=2748.361702.81=1.614m ∆÷所以，第一次近似中和轴距基线的距离为1.614m 船体剖面对水平中和轴的惯性矩为：222=2(9951.42138.512748.361702.81)11308.1cm m I ⨯+-÷=⋅剖面上各构件的应力为：'i i =/100M Z Iσ 式中'i i Z Z =-∆ 2、临界压力计算由于该计算中船舶船舯处于波谷中，即船舶处于中垂状态，所以下面只列出中和轴以上部分受压板的临界应力。
纵骨架式板格（四边自由支持）按下式计算：210076()cr t bσ=3、船体总纵弯曲应力第二次近似计算（1）剖面折减系数计算已知本船体结构为纵骨架势，因此对于只参加抵抗总纵弯曲的构件criσϕβσ=式中 cr σ——板格的临界应力i σ——与所计算的板在同一水平线上的刚性构件总纵弯曲压应力的绝对值β——系数，275btβ=-，若1β>，则取1β= b ——纵骨间距 计算结果见下表：（2）总纵弯曲应力第二次近似计算由上表知本船在波谷位置时，甲板板发生折减。
船舶结构与强度设计 第3章PPT课件
第7页/共27页
纵向强力构件上有开口，计算时应扣除吗？
16b球扁钢剖面积21.16cm2，自身惯性矩527 cm4， 重心轴9.75 cm。包括带板的剖面惯性矩2489 cm4。
E
2 2.06105 2489
240221.16 801.4
659N
/
m m2
欧拉应力超过屈服极限，修正后得临界（屈曲）应力。
cr
s 1
s 4 E
214N
/
m m2
当压缩应力达到σcr=214N/mm2 时，甲板纵骨连同带 板将发生像压杆一样的失稳。
船体剖面特性计算应包括哪些构件？ 纵向强力构件——纵向连续并能够有效地传递 抗总纵弯曲应力的构件，即船中0.4L范围内纵向连续 构件。如，甲板、外板、内底板、纵骨和纵桁等。
——只能包括纵向强力构件！
第4页/共27页
ቤተ መጻሕፍቲ ባይዱ
船体梁剖面特性计算原理同简单型材（如工字型 材）计算，可列表计算。
第5页/共27页
构件
2.5×14 1.5×16 160×14/40×14 1.0×16 7.2×14 4.0×12 0.5×25 0.4×25 0.5×25 0.5×25
0.035 0.024 0.0084 0.0160 0.1008 0.0480 0.0125 0.0100 0.0125 0.0100
9.0 5.654 9.0 5.654 8.9 5.554 8.5 5.154 4.4 1.054 5.5 2.154 8.75 5.404 8.5 5.154 5.25 1.904 5.0 1.654
《船舶强度与结构设计》课程设计-1500m3耙吸式挖泥船总强度计算书
《船舶强度与结构设计》课程设计题目：1500m3耙吸式挖泥船总强度计算书学院专业年级姓名学号目录第一章计算说明 (1)1.1 计算内容 (1)1.2 主要技术参数 (1)第二章剪力和弯矩计算 (1)2.1 重力分布和浮力分布计算 (1)2.2 静水剪力和弯矩计算 (3)2.3 波浪附加剪力弯矩计算及剪力弯矩合成 (7)第三章总纵弯曲应力计算 (11)3.1 剖面参数 (11)3.2 总纵弯曲应力计算 (13)第四章临界应力计算和构件稳定性校核 (14)4.1 纵骨架式甲板临界应力计算及校核 (14)4.2 甲板纵桁临界应力计算及校核 (15)4.3 纵骨临界应力计算及校核 (15)第五章极限弯矩计算 (16)5.1 极限弯矩下各构件应力计算 (16)5.2 极限弯矩计算 (17)第六章计算结果分析及结论 (18)第一章 计算说明本计算书是1500m 3耙吸式挖泥船总强度计算书，计算出了中拱状态下的船体的静水剪力、弯矩，波浪附加剪力、附加弯矩，合成剪力、合成弯矩，并计算了总弯曲应力，以及校核了是否满足设计要求。
1.1计算内容(1) 静水弯矩、剪力 (2) 波浪附加弯矩、附加剪力 (3) 剪力、弯矩合成(4) 计算总弯曲应力、受压构件的稳定性校核 (5) 计算结果分析及结论 (6) 计算状态：中拱1.2主要技术参数船长：78米；满载排水量：5020吨；平均吃水：5.4米；站距：9.3=∆L 米，波高：4米；重心在舯前：813.0=g x 米;艏吃水:77.5=f T 米;尾吃水:23.5=a T 米。
主尺度：船长：78米，船宽：14.5米，型深：6.3米，设计吃水5.1米，肋距：0.6米，强框架间距：1.8米，纵骨跨度1.8米。
海水密度ρ=1025 kg/m 3，重力加速度g=9.81m/s 2，船体钢材屈服强度σy=235N/mm 2第二章 剪力和弯矩计算2.1重力分布和浮力分布计算静水剪力弯矩计算资料：表2：静水平衡状态各站横剖面浸水面积（m2）根据表1的数据将重量单位t转换为KN，绘制站间重量分布曲线p(x)：图1：重量分布曲线依据表2计算静水浮力，取相邻两站号横剖面浸水面积的平均值作为此站距间的浸水面积A，根据理论站间浮力公式F=ρgΔLA计算出各站间的静水浮力，计算过程如表3所示：表3：静水浮力计算图2：静水浮力分布曲线2.2静水剪力和弯矩计算载荷分布曲线q(x)=p(x)−b(x)，各站间重量减去浮力得载荷值，由此绘制载荷分布曲线：图3：载荷分布曲线静水剪力是载荷分布的一次积分，即N(x)=∫q(x)dxx静水弯矩是剪力的一次积分，载荷的二次积分，即M(x)=∫N(x)dxx0=∫∫q(x)dxdxxx利用表格进行剪力和弯矩的计算，第20站剪力不为0且N(20)N max=2.92%，需要进行剪力修正：ΔN i=i20N20 N i’=N i+ΔN i使船艏艉的剪力为0，同样第20站弯矩也不为0，M(20)M max=3.07%，再次进行弯矩修正：ΔM i=i20M20 M i’=M i+ΔM i计算过程如表4所示：表4：静水剪力和弯矩计算56由上表修正后的数据，绘制站间静水剪力和弯矩图：图4：静水剪力分布曲线图5：静水弯矩分布曲线2.3波浪附加剪力弯矩计算及剪力弯矩合成波浪附加剪力弯矩计算资料：2波浪附加浮力按如下公式计算：Δb(x)=b w (x)−b s (x)其中b w (x)为波浪中平衡时的浮力曲线，b s (x)为静水中平衡时的浮力曲线 波浪附加剪力是附加浮力负数的一次积分：N w (x)=∫[b s (x)−b w (x)]dx x=∫[−Δb(x)]dx x浮力与浸水体积成正比，因此站号上静水浮力与波浪中浮力之差可以用浸水面积之差来表示：N w (x)=γ∫[F s (x)−F c (x)]dx x=γ∫ΔF i (x)dx x其中F si (x)是静水中浸水面积分布曲线，F ci (x)是波面下浸水面积分布曲线 波浪附加弯矩是附加剪力的一次积分，附加浮力的二次积分：M w (x)=∫N w (x)dx x 0=∫∫[−Δb(x)]dxdx x 0x 0同样也可以用浸水面积曲线来表示M w (x)=γ∫∫ΔF i (x)dxdx x 0x 0下面列表计算，在计算过程中，发现波浪附加剪力和附加弯矩同样在船艏处不闭合，N w(20)N wmax=1.59%，M w(20)M wmax=6.59%，采用相同的方法进行线性修正，以确保船艏艉剪力和弯矩都为0。
《船体结构与强度》PPT课件
面上时，从外表上看，重力与浮力 大小相等并处于平衡状态，但事实 上组成船体各分段的重力与浮力的 最终平衡值通常是不相等的。这种
重力与浮力沿船长方向的不均匀分 布，在产生剪切应力的同时，也产 生了总纵弯曲力矩，使船体发生总 纵弯曲。弯矩的最大值在船中附近，
3〕纵向构件：在船体构造中，承担总纵弯曲强度的 构件称为纵向构件，有：甲板纵桁、甲板纵骨、舷 侧纵桁、舷侧纵骨、船底纵桁、船底纵骨、中内龙 骨、旁内龙骨、甲板、内底板、纵向舱壁、船体外 板等。
4〕横向构件：在船体构造中承担横向强度的构件称 为横向构件，有：横梁、强横梁、肋骨、肋板、梁 肘板、舭肘板、横舱壁等。
《船体结构与强度》PPT 课件
本课件PPT仅供大家学习使用 学习完请自行删除，谢谢！ 本课件PPT仅供大家学习使用 学习完请自行删除，谢谢！ 本课件PPT仅供大家学习使用 学习完请自行删除，谢谢！ 本课件PPT仅供大家学习使用 学习完请自行删除，谢谢！
一、船体强度的根本概念
船体强度是指船体构造抵抗各种外力 作用的能力。
内力分为两种，一种是由内部上下移动趋 势构成的内力，称为剪力；另一种是由弯 曲变形趋势引起的内力，称为弯曲力矩。
2〕总纵弯曲力矩和剪力
〔2)船体总纵弯曲变形的形式 船体总纵弯曲变形的形式有中拱和中垂两种。
2〕总纵弯曲力矩和剪力
〔3) 最不利的浮力和重力的大小及分布
① 浮力的大小和分布:当船舶在海上遇到标准波：波的形状 为坦谷波，即波峰较陡而波谷平坦，且波长λ等于船长L，波 高 H 等 于 波 长 λ 的 1 ／ 20( L≥120m 时 ) 或 波 高 等 于 λ ／ 30+2m(L＜120m时)，船与波的相对位置是波峰位于船中 或波谷位于船中时，船舶所受到的浮力分布对船体总纵弯曲 强度是最为不利的。
船舶强度课程设计作业-波浪附加弯矩与总弯矩计算等
1. 波浪附加弯矩与总弯矩计算： 1.1 波浪附加弯矩计算：根据《钢质海船入级建造规范》中给出了船体梁各个横剖面的中拱波浪弯矩)(+w M 与中垂波浪弯矩)-(w M 的计算公式：可根据公式1-1，1.2求得波浪附加弯矩。
m kN BC CL M b w .10M 190)(32-⨯+=+ （中拱） （公式1-1）m kN C B CL M b w .10)7.0(M 110)-(32-⨯+-=（中垂） （公式1-1）其中，M ——弯矩分布系数，对于船中M=1；图1-1 M-弯矩分布系数 b C ——方型系数，取值不得小于0.6，对于本船为0.67；C ——系数，按公式1-2计算，得8.4548。
(按L=126m) 2/310030075.10⎪⎭⎫ ⎝⎛--=L C （m L m 30090≤≤） （公式1-2）1.1.1 中拱情况下波浪附加弯矩计算：1.0M1.00.65 0.4尾1.1.2 中拱情况下总弯矩计算：根据条件，可以近似将船体的静水弯矩当作是波浪附加弯矩的1/3，所以，可以得出总弯矩是波浪附加弯矩的4/3倍。
1.1.3中垂情况下波浪弯矩的计算1.1.4中垂情况下总弯矩计算2. 总弯曲应力计算与稳定性校核根据中横剖面图，计算剖面抗弯几何特性，见下页表1.5。
参考轴取船底板上表面，可求得：A=9390.3cm 2、B=4536030.6cm 3、C=3837423286cm 4。
参考轴距中和轴距离⊿=B/A=483.0549226cm,I=2(C-⊿2 A)= 3292542712cm 4任意构件距中和轴距离Z i ’=Zi-⊿由此求得中横剖面各构件剖面模数W i =I/ Z i ’表1.52.1 计算总弯曲应力：总弯曲应力按下式计算： cws W M M +=σ由表6.1得最小剖面模数Wc 为5469839.086cm 3,因此在中拱情况下最大总弯曲应力为91.63N/mm 2 ，在中垂情况下最大总弯曲应力为108.48N/mm 2，远小于许用应力。
船体强度与结构设计--船体剪力弯矩计算
《船体强度与结构设计》课程作业一船体剪力弯矩计算一、计算过程1、程序代码function [W,Dis,Xg,Xf,Df,Da,Dav,DDDD,Fl,Dvalue,Jdvalue,Mdvalue]= Strength()% 用于计算在已知载荷和邦戎曲线时的船舶首尾吃水；% 绘出载荷、浮力、剪力、弯矩分布曲线；% 输入：% B:邦戎曲线分布矩阵；% L:载荷分布矩阵；% A:水线面面积；% Ps:水的重量密度密度，默认1.0t/m?；% Dx:站距；% Dz:吃水间距；% Eps:计算精度，默认0.001%；% 输出：% W:总重量；% Dis:排水量；% Xg:重心坐标；% Xf:浮心坐标；% Df:首吃水；% Da:尾吃水；% Dav:平均吃水；% DDDD;迭代数据储存；% Fl:每站浮力数据；% Dvalue:每站载荷与浮力差值数据；% Jdvalue:每站中间剪力数据；% Mdvalue:每站中间弯矩数据；% 距离单位：m%**********为方便计算，直接将所需数据写入函数中***********Eps=0.00001;B=[0 0 0.36 0.81 1.28 1.76 2.26 2.79 3.33 3.89 4.46 5.06 5.67 6.3 6.957.62 8.31 9.01 9.74 10.48 11.24 12.02 15.24 21.65 29.990 0 0 0.39 1.28 2.37 3.55 4.73 5.81 6.67 7.27 7.69 8.03 8.33 8.66 9.049.56 10.35 11.83 14.97 20.62 28.73 38.55 49.61 61.60.01 1.3 3.24 5.52 8.02 10.65 13.35 16.06 18.75 21.4 24.0426.73 29.59 32.8 36.59 41.3 47.37 55.01 64.26 74.8686.51 98.97 112.06 125.67 139.680.03 2.17 5.24 8.87 12.93 17.36 22.08 27.07 32.33 37.9443.98 50.54 57.82 66.15 75.6 86.14 97.63 109.89 122.78136.19 150.06 164.28 178.78 193.48 208.310.07 3.77 8.69 14.44 20.88 27.94 35.58 43.79 52.6262.14 72.4 83.41 95.12 107.45 120.32 133.68 147.45 161.6 176.04 190.71 205.54 220.43 235.35 250.27 265.180.13 6.04 13.61 22.26 31.8 42.09 53.06 64.67 76.8589.55 102.68 116.21 130.08 144.27 158.72 173.36 188.15 203.01 217.91 232.83 247.74 262.66 277.57 292.49 307.410.21 8.83 19.22 30.65 42.85 55.65 68.94 82.61 96.59 110.83 125.28 139.89 154.61 169.4 184.23 199.07 213.94 228.82 243.72 258.63 273.55 288.47 303.38 318.3 333.210.29 11.22 23.59 36.74 50.42 64.47 78.78 93.28 107.92 122.63 137.39 152.16 166.95 181.76 196.59 211.44 226.3 241.18 256.08 271 285.91 300.83 315.74 330.66 345.58 0.34 12.58 25.93 39.86 54.14 68.64 83.26 97.96 112.67 127.41 142.17 156.94 171.73 186.54 201.37 216.22 231.08 245.96 260.86 275.78 290.69 305.61 320.52 335.44 350.360.35 12.73 26.19 40.2 54.53 69.07 83.71 98.41 113.13 127.86 142.62 157.39 172.18 186.99 201.82 216.67 231.53 246.41 261.31 276.23 291.14 306.06 320.98 335.89 350.810.35 12.73 26.18 40.19 54.52 69.05 83.69 98.39 113.11 127.85 142.6 157.38 172.17 186.98 201.81 216.65 231.52 246.4 261.3 276.21 291.13 306.04 320.96 335.88 350.790.32 12.03 24.98 38.57 52.58 66.89 81.4 96.05 110.76 125.49 140.25 155.02 169.82 184.62 199.45 214.3 229.16 244.04 258.94 273.86 288.77 303.69 318.61 333.52 348.440.27 10.5 22.24 34.82 47.98 61.57 75.5 89.69 104.09 118.64 133.31 148.07 162.86 177.66 192.49 207.34 222.2 237.08 251.98 266.9 281.81 296.73 311.65 326.56 341.480.2 8.4 18.34 29.32 41.05 53.37 66.17 79.37 92.91 106.73 120.78 135.04 149.46 164.02 178.69 193.44 208.26 223.13 238.03 252.94 267.86 282.78 297.69 312.61 327.520.12 5.95 13.6 22.4 32.05 42.41 53.35 64.8 76.71 89.02 101.69 114.68 127.96 141.48 155.24 169.21 183.35 197.66 212.11 226.68 241.36 256.12 270.95 285.83 300.730.07 3.83 9.18 15.56 22.77 30.67 39.18 48.25 57.8367.88 78.38 89.29 100.59 112.25 124.26 136.59 149.21 162.12 175.28 188.69 202.31 216.14 230.16 244.35 258.680.04 2.37 5.87 10.2 15.2 20.79 26.91 33.5 40.56 48.06 55.96 64.28 72.98 82.07 91.54 101.39 111.6 122.17 133.1 144.37 155.97 167.9 180.14 192.68 205.510.02 1.42 3.58 6.3 9.5 13.11 17.11 21.45 26.13 31.1336.45 42.08 48.03 54.3 60.89 67.81 75.05 82.64 90.5698.83 107.45 116.44 125.8 135.53 145.660.01 0.85 2.14 3.77 5.7 7.91 10.35 13.02 15.91 19.01 22.3125.81 29.49 33.35 37.39 41.61 46.02 50.62 55.4460.49 65.8 71.4 77.3 83.55 90.160 0.48 1.27 2.32 3.63 5.18 6.68 6.68 6.68 6.82 9.22 11.7 14.23 16.7619.24 21.65 23.97 26.19 28.36 30.53 32.78 35.1437.68 40.44 43.460 0 0 0.15 0.71 1.62 2.83 4.3 6 7.9 9.94 12.08 14.28 16.4818.63 20.64 22.4 23.79 24.64 24.95 25.04 25.05 25.0525.05 25.1];L=[791.9 810.6 855.8 900.0 399.6 488.1 525.6 410.1 558.9 750.1 710.5 670.8 540.9 650.0 558.6 610.2 496.5 580.0 520.4 419.3 ]T;A=4800; Ps=1.0; Dx=10; Dz=0.5;Bh=size(B,1); %读取B的行数Ls=size(L,1); %读取L的行数W=0; Wm=0;DDDD=zeros(2,40); %用于储存每次Df，Da迭代数据，便于分析for ii=1:20W=W+L(ii);Wm=Wm+(ii-10.5)*L(ii);endXg=Wm/W*Dx; %重心坐标Df=1.5*W/A;Da=Df;DDDD(1,1)=Df; DDDD(2,1)=Da;jj=1;while(1) %迭代计算首尾吃水准确值Dis=0;Dm=linspace(Da,Df,Bh);Dml=fix(Dm/Dz);Ddd=Dm-Dml/2;Space=zeros(1,Bh);Moment=0;for ii=1:BhSpace(ii)=B(ii,Dml(ii)+1)+(B(ii,Dml(ii)+2)-B(ii,Dml(ii)+1))/Dz*Ddd(ii);%计算各站吃水面积Dis=Dis+Space(ii);Moment=Moment+Space(ii)*(ii-11);endDis=(Dis-0.5*Space(1)-0.5*Space(ii))*Dx; %计算排水量Moment=(Moment-0.5*Space(1)*(-10)-0.5*Space(ii)*(ii-11))*Dx*Dx;%计算对中站面转矩Xf=Moment/Dis; %计算浮心位置if((abs(Xg-Xf)/Dx/(Ls))<Eps&&abs((W-Dis)/W)<Eps)%判断是否满足精度要求breakelseDf=Df+(W-Dis)/A-18*(Xf-Xg)*W/A/Ps/Dx/Ls; %修正首吃水Da=Da+(W-Dis)/A+18*(Xf-Xg)*W/A/Ps/Dx/Ls; %修正尾吃水DDDD(1,jj+1)=Df;DDDD(2,jj+1)=Da;jj=jj+1;endendDav=(Df+Da)/2; %计算平均吃水、% ************画图部分*****************%画图数据处理Lt=-L'; %载荷数据Fl=zeros(1,Ls); %每站浮力数据Dvalue=zeros(1,Ls+2); %每站载荷与浮力差值for ii=1:LsFl(ii)=(Space(ii)+Space(ii+1))/2*Dx;Dvalue(ii+1)=Fl(ii)+Lt(ii);endJdvalue=zeros(1,Ls+2);Mdvalue=zeros(1,Ls+2);%*********************计算剪力********************************Jdvalue(2)=(Dvalue(1)+2*Dvalue(2))/4;for ii=3:Ls+1for jj=2:ii-1Jdvalue(ii)=Jdvalue(ii)+(Dvalue(jj)+Dvalue(jj+1))/2;endJdvalue(ii)=Jdvalue(ii)+Jdvalue(2);endJdvalue(ii+1)=Jdvalue(ii)+(Dvalue(ii+1)+2*Dvalue(ii))/4;%*******************修正剪力********************************JJJJ=Jdvalue;Dvaline=linspace(Jdvalue(1),Jdvalue(22),Ls+2);Jdvalue=Jdvalue-Dvaline;%*******************计算弯矩********************************Mdvalue(2)=(Jdvalue(1)+Jdvalue(2))/4*Dx;for ii=3:Ls+1for jj=2:ii-1Mdvalue(ii)=Mdvalue(ii)+(Jdvalue(jj)+Jdvalue(jj+1))/2*Dx;endMdvalue(ii)=Mdvalue(ii)+Mdvalue(2);endMdvalue(ii+1)=Mdvalue(ii)+(Jdvalue(ii+1)+2*Jdvalue(ii))/4*Dx;%*******************修正弯矩********************************MMMM=Mdvalue;Mvaline=linspace(Mdvalue(1),Mdvalue(22),Ls+2);Mdvalue=Mdvalue-Mvaline;%*******************绘图辅助数据*************************x0=-95:10:95;x1=-105:10:105;x1(1)=-100;x1(22)=100;y0=zeros(1,20);y1=zeros(1,22);%*************************画图部分*************************%**********重量分布与浮力分布图*********axis([-Dx*Bh*0.55,Dx*Bh*0.55,-5,5]);hold onquiver(x0,y0,y0,Lt,0.1,'LineWidth',2),quiver(x0,y0,y0,Fl,0.1,'LineWidth',2); figure%**********重量与浮力之差分布图*********axis([-Dx*Bh*0.55,Dx*Bh*0.55,-5,5]);hold onquiver(x1,y1,y1,Dvalue,0.1,'LineWidth',2); figure%**************剪力分布图**************axis([-Dx*Bh*0.55,Dx*Bh*0.55,-5,5]);hold onquiver(x1,y1,y1,Jdvalue,0.1,'LineWidth',2); figure%**************弯矩分布图**************axis([-Dx*Bh*0.55,Dx*Bh*0.55,-5,5]);hold onquiver(x1,y1,y1,Mdvalue,0.1,'LineWidth',2); figurexxl=-100:1:100;yyl=interp1(x1,Jdvalue,xxl,'pchip');yyyl=interp1(x1,Mdvalue,xxl,'pchip');plot(xxl,yyl); %画连续剪力分布图figureplot(xxl,yyyl); %画连续弯矩分布图end2、输出结果整理（1）、迭代数据：（2）、详细计算结果53、所绘图像（1）、每站载荷和浮力分布示意图（纵坐标均仅表示相对值大小，不表示真实值）（2）、每站载荷与浮力差值分布图（4）、弯矩分布图（6）、弯矩分布连续图像二、计算过程感想1、此次作业的核心是通过迭代求得首尾吃水的值，是之满足精度要求。
船舶剪力弯矩计算表格
船舶剪力弯矩计算表格摘要：一、船舶剪力弯矩计算表格的概念与意义二、船舶剪力弯矩计算表格的内容与结构三、船舶剪力弯矩计算表格的应用与实例四、船舶剪力弯矩计算表格的优缺点分析五、船舶剪力弯矩计算表格的发展趋势正文：一、船舶剪力弯矩计算表格的概念与意义船舶剪力弯矩计算表格是船舶结构设计中常用的一种计算工具，用于计算船舶在各种载荷条件下，船体结构各部位的剪力弯矩。
剪力弯矩是船舶结构在受力时产生的一种内力，它是船舶结构强度计算的重要依据。
通过船舶剪力弯矩计算表格，可以方便、快捷地获取船舶结构各部位的剪力弯矩值，为船舶结构设计提供重要参考。
二、船舶剪力弯矩计算表格的内容与结构船舶剪力弯矩计算表格主要包括以下几个部分：1.船舶的基本参数：包括船舶的类型、尺寸、排水量等，这些参数是计算船舶剪力弯矩的基础。
2.载荷条件：包括船舶在航行、停泊等不同状态下的载荷情况，如货物、燃油、水等。
3.船舶结构部件：包括船体的各个部位，如船首、船尾、船舷、船底等。
4.剪力弯矩计算公式：根据船舶结构部件和载荷条件，选择相应的剪力弯矩计算公式，计算出各部位的剪力弯矩值。
5.剪力弯矩计算结果：将计算出的剪力弯矩值填入表格，便于查看和分析。
三、船舶剪力弯矩计算表格的应用与实例船舶剪力弯矩计算表格在船舶结构设计中有广泛的应用，以下是一个简单的实例：假设有一艘货船，在航行过程中，船上装载了1000 吨货物。
我们需要计算船体在航行过程中，船首部位的剪力弯矩值。
首先，根据船舶的基本参数，确定船体的尺寸、排水量等。
然后，根据载荷条件，将货物的重量纳入考虑。
接着，选择合适的剪力弯矩计算公式，根据船首部位的结构特点进行计算。
最后，将计算出的剪力弯矩值填入表格。
四、船舶剪力弯矩计算表格的优缺点分析船舶剪力弯矩计算表格的优点在于方便、快捷地获取船舶结构各部位的剪力弯矩值，为船舶结构设计提供重要参考。
同时，计算表格具有一定的通用性，可以适用于不同类型的船舶。
船体静波浪剪力和弯矩的计算讲解
（二）坦谷波的绘制方法： 坦谷波为：车轮滚动时，轮盘内任一点的运动轨迹。 1. 按坦谷波面方程原理
其公式如下：【推导】
x r sin 2 y r cos
V
2. 坦谷波曲线的计算表 1-3 【p24】 按波长/波高 比的不同； 求各站的y/λ值，制成表格。
（三）静波浪剪力及弯矩计算
1、船舶在波浪中浮态的轴线 【假设以静水线作为坦谷波的轴线，发现不能平衡】 根据坦谷波的特点，坦谷波在波轴线以上的剖面积 比在轴线以下的剖面积小。 1）船中位于波谷时： 中部较两端丰满，船舶在此位置时的浮力要比在静 水中减少许多，因而不能处于平衡，船舶将下沉ζ值 2）船中位于波峰时：相反，一般船舶要上浮一些。 3）同时，船体首尾线型不对称使船舶发生纵倾变化。
（四）船舶斜置在波浪上的静波浪弯矩计算
对于船长大于波长的内河船，需要用将船舶斜置于波浪上的 方法进行静波浪弯矩计算，斜置的目的在于使船体受力最不利。 斜置的影响： 在各个非船中剖面，浮力沿船宽的分布不是均匀的了，而 是按坦谷曲线。因此船舶除受到总纵弯曲力矩的作用外，还 将受到扭转力距的作用。
（五）波浪浮力修正(或称史密斯修正) 1. 考虑了波浪的惯性力； 2. 修正之后浮力曲线将会变得更平坦（不论波峰或波谷）， 因而静波浪弯矩与剪力也将变小。 3. 修正后反而偏危险！！ 4. 结果表明： 一般船舶在满载吃水时，静波浪弯矩可减少２０％～ ３０％左右，而总纵弯矩大约减少１０％～１５％
20 L V xb 3 0 6 b 7 2 20 L
利用表格计算出上述五个积分系数后，可由上式解出和值，于 是就得到了船舶静置于波浪上的实际平衡位置。 特别提示：在进行上述表格计算时注意各符号代表的意义， 如有可能请自己推导出来。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
船舶强度与结构设计大作业（一）船海1301 禹宗昕U201312263一．数据与函数准备1.主要数据：计算船长L=200m水密度ρ=1t/m³重力加速度g=9.8m/s²2.原始资料（见附录）（1）全船重量分布（2）全船邦戎曲线数据3.静水平衡参数总重量W=122248 *9.8=120030.4 kN水线面面积A=4800m平均吃水dm=3.9m纵稳心半径R=220m漂心纵向坐标xf=4.3m4.基本函数简介（1）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:nD(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 = ∑M i 201∑W i 201= E23/C24*10 = 93.2136 m 2.浮态计算（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%时停止迭代，迭代步长：0.3 while (index1>0.005 )||(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+((W-V)/(p*A)-(L/2+xf)*(xg-xb)/R)*0.3;D=draft(d.df,d.da,(L/zhanju+1));area=SAREA(D);V=displacement(area,zhanju);xb=flocenter(area,zhanju);d.xg=xg-100;d.xb=xb-100;d.weight=W;d.displacement=V;index1=abs((W-V)/W);index2=abs((xg-xb)/L);d.index1=index1;d.index2=index2;k=k+1;d.k=k;end;end（2）运行结果d = floatingstate (3.9,200,4800,220,-4.3,12248,93.214,10)d =df: 3.0666da: 4.6647xg: -6.7860xb: -6.6036weight: 12248displacement: 1.2240e+004index1: 6.5099e-004index2: 9.1211e-004k: 3浮态分析：da>df 尾倾精度分析：index1=abs((W-V)/W)= 0.065% < 0.5%Index2=abs((xg-xb)/L)=0.091% < 0.1%迭代次数：k=3三．剪力与弯矩计算1.初次计算（1）.载荷分布，剪力，弯矩的计算与储存% D：一维矩阵，已修正浮态下各站吃水，% area：一维矩阵，已修正浮态下各站面积% earea：一维矩阵，利用area和三次样条曲线差值，求0.5,1.5..站面积% eachv：各站的体积，*10为站距% qx：载荷分布，重力-浮力，qx1为集中力到分布力的换算% N11：剪力函数，M11：弯矩函数% A，B11，C11：分别为qx1，N11，M11对应的多项式系数矩阵D=draft(d.df,d.da,21)；area=SAREA(D)；earea=spline([0:1:20],area,[0.5:1:19.5])；eachv=10*earea；qx=eweight-eachv；qx1=qx/10*9.8;x=[5:10:195];A11=polyfit(x,qx1,4);syms x q1;q1=poly2sym(A,x);syms N11 M11;N11=int(q1);M11=int(N11);B11=sym2poly(N11); C11=sym2poly(M11);x=[0:10:200];ezplot(M11,x) ;hold on;ezplot(100*N11,x);（2）图形校核原理：1.多项式A的根为N的极值点，B的根为M的极值点，a=roots(A)2.将极值点带入函数求极值，如polyval(B,a(2))3.求第20站的剪力，弯矩值，如polyval(B,20)4.计算k1，k2指标，k1=Nmax/N（20）；k2=Mmax/M（20）；5.k1&&k2<0.05时可进行第二步修正a=roots(A11)a =20.875014.2251-7.78004.0257>> polyval(B11,a(2))ans =-1.6803e+004>> polyval(B11,a(4))ans =1.9640e+004>> polyval(B11,20)ans =6.8987e+001>> k1=ans/ 1.9640e+004k1 =0.0035>> b=roots(B)b =-11.665421.661619.97159.2145>> k2=polyval(C,20)/polyval(C,b(5))k2 =0.0470（3）结果分析k1=Nmax/N（20）=0.0035 <0.05k2=Mmax/M（20）=0.0470 <0.05由图可见，20站处N和M均不为0，即此积分方法存在误差，但均在误差允许范围内，因此可进行二次修正。
2.修正原理：（1.）N1=N-x*N（20）/20 ；M1= M-x*M（20）/20（2）N1对应B11为五次多项式，因而一次项为第五项，M1对应C11为六次多项式，因而一次项为第六项，B22=B11;C22=C11;B22(5)= B22(5)- polyval(B11,200)/200;C22(6)=C22(6)- polyval(C11,200)/200; syms x N1 M1;N1=poly2sym(B22,x);M1=poly2sym(C22,x);x=[0:10:200];ezplot(M1,x) ;hold on;ezplot(100*N1,x);x=[0:1:20];ezplot(M1,x) ;hold on;ezplot(100*N1,x);分析：利用matlab自带的图像分析功能，易得修正后，M1max=1.164*106 KN*mN1max=2*104 KN四．心得体会船舶强度与结构设计已经开课三周，给我最深的感触是，它突破了原先我们结构设计的简单认识，以为设计出来就万事大吉。
课程里，老师无数次强调反复修正的重要性，然而这门课的大作业也是我们上大学以来第一次面临计算的迭代问题。
因为在计算工具上，我选择了matlab。
之前在静力学课设上，同学们普遍使用AutoCAD绘制邦戎曲线，然后描点，量尺寸，求每一站的面积，然而如今面临迭代，手工求法显然不再使用。