利用PCL語言查看船體剖面特性

唐旭東

（上海交通大學 上海 200030）

利用PCL語言查看船體剖面特性

唐旭東

（上海交通大學 上海 200030）

MSC.Patran；PCL；剖面特性；慣性矩

PCL（Patran Command Language）語言為MSC.Patran的二次開發提供了強大的功能。利用PCL語言編寫程序來計算船體剖面慣性矩、剖面模數及中和軸位置等剖面特性。

0 引言

船體梁剖面特性包括剖面面積、剖面慣性矩、中和軸及剖面模數等，這些參數在船舶結構設計中有著重要的作用，如船體總縱強度計算中需要校核剖面模數，艙段有限元分析中創建端面約束或查看剖面內彎矩時需用到中和軸位置等。

目前，一些船級社的專用軟件可以計算剖面特性，如DNV船級社的Nauticus軟件和BV船級社的Mars 2000軟件等。然而隨著有限元分析方法的應用日益廣泛，我們在進行結構計算，特別是利用MSC.Patran做艙段有限元分析、優化時，可能會頻繁地對列板的寬度、厚度或扶強材的尺寸進行調整，此時，若能直接用MSC.Patran本身來查看船體剖面特性，將會帶來諸多方便。本文正是基于這點，利用PCL語言編寫程序，開發了MSC.Patran的擴展界面，通過選擇組或單元，計算出船體剖面慣性矩、甲板處和船底處的剖面模數、中和軸位置、有效截面積以及簡單的開口薄壁截面在單位剪力作用下的剪應力分布等剖面特性。

1 原理概述

MSC.Patran本身不具備查看模型剖面特性的功能，但可以查詢選定單元的質量、質心位置和轉動慣量等物理特性，而這些物理特性在一定條件下可以推導出剖面特性，本文主要闡述由轉動慣量推出剖面慣性矩的原理，其余如中和軸位置、剖面模數就迎刃而解了。

1.1 由轉動慣量推出剖面慣性矩

轉動慣量與剖面慣性矩是兩個不同的概念。轉動慣量是剛體轉動時慣性的度量，它等于剛體內各質點的質量與質點到軸的垂直距離二次方乘積之和，對于質量連續分布的剛體，轉動慣量可表示為J=m∫r2d m。剖面慣性矩是剖面的幾何特性，指剖面各微元面積與各微元至剖面上某一指定軸線距離二次方乘積的積分，可表示為I=A∫r2d A。從數學上來分析，轉動慣量與剖面慣性矩都是二次矩，在一定條件下可以通過轉動慣量推出剖面慣性矩。

如圖1所示，CG為剛體的質心，xyz為過CG的坐標系，剛體在x方向是等截面的，長度為l，剛體繞y軸的轉動慣量可以表示為：

圖1 剛體的轉動慣量

式中：

d m為微質量，kg；可表示為d A·d x·ρ。d A為yz剖面上的微面積，m2；ρ為剛體密度，kg/m3。

r為微質量到質心的距離，m；表示為r2=x2+z2。

所以，式（1）可以表示為：

對式（2）中的二重積分進一步轉化，可得：

式（3）中：A∫z2d A即為yz剖面內繞y軸的慣性矩，用Iy表示，則：

式（4）的成立必須符合如下條件：均質剛體，沿x方向等厚、等截面。在利用MSC.Patran建立船體模形時，材料一般為鋼材，密度相同，沿x方向一般按肋位劃分網格并且相鄰兩個肋位的截面面積可以認為相等，因此可以用式（4）來計算剖面慣性矩。

1.2 組合剖面的慣性矩

在計算船體剖面慣性矩時，并非所有的構件都是按100%剖面積計入的，船級社對此有特別規定。例如中國船級社規定：在計算耙吸挖泥船剖面慣性矩時，對兩端有良好過渡的箱形龍骨可計入85%的剖面積。若把計入百分比記為f，在計算時，可以按不同的f值將船體構件分為若干組，再用組合剖面的慣性矩公式求出組合慣性矩。根據定積分的性質及平行移軸公式，組合剖面慣性矩可表示為：

式中：

Iyi為第i組剖面相對于組合后中和軸的慣性矩，m4；

Iy0i為第i組剖面相對于自身中和軸的慣性矩，m4；

ai為第i組剖面自身中和軸與組合后的中和軸之間的距離，m；

Ai為第i組剖面面積，m2；

fi為第i組剖面計入組合慣性矩的百分比。

組合后的中和軸：

式中：Sy為組合剖面的靜矩，m3；

A為組合剖面的有效面積，m2；

Ai為第i組剖面面積，m2；

zi為第i組剖面中和軸，m。

1.3 開口薄壁截面的剪應力分布

根據材料力學的理論，開口薄壁構件受彎時的剪應力τ可假定沿壁厚方向均勻分布，其大小為：

式中：Qz為計算截面z軸主平面內的剪力，N；

Sy為計算剪應力處以上（或以下）截面對中和軸y軸的面積矩，m3；

Iy為繞y軸的截面慣性矩，m4；

t為計算點處的板厚，m。

對于閉口薄壁構件受彎時的剪應力分布計算相對復雜，利用PCL語言來實現有一定困難，還有待進一步研究。

2 程序的實現

MSC.Patran的PCL（patran command language）語言為Patran的二次開發提供了強大的功能。用戶可以通過PCL語言調用Patran的數據庫、調用Patran的內部函數、創建用戶面板以及自定義函數等。

本程序主要調用了以下函數：

ui_form_create創建用戶面板“剖面特性”；

ui_switch_create創建選擇開關

ui_text_create創建文本框

ui_button_create創建執行按鈕

mass_prop_create提取質量特性

res_utl_create_elem_result生成結果文件

res_display_fringe_create生成應力云圖

section_characters自定義函數，生成剖面特性

section_shear自定義函數，生成剪力流云圖

section_characters和section_shear兩個自定義函數是把式（4）~（7）程序化，根據mass_prop_create返回所選單元的質量、質心、轉動慣量等參數計算出剖面的慣性矩、中和軸位置等剖面特性及生成剪力流云圖。在p3toolbar.def中創建工具欄按鈕，并指向自定義函數生成程序界面，如圖2所示。

圖2 程序界面

3 程序的應用

圖3 某耙吸挖泥船泥艙典型橫剖面

圖3為某耙吸挖泥船艙段有限元分析中截取的泥艙典型剖面，其中，底部箱形龍骨參與剖面慣性矩的有效比例為85%。在Patran定義兩個組Moment_Mid*0.85和Moment_Mid，并將箱形龍骨處的單元加入Moment_Mid*0.85組中，其余單元加入Moment_Mid中，程序會根據組名稱判斷該組內單元參與剖面慣性矩的有效百分比。點擊工具欄圖標調出程序界面（如圖2），選擇Moment_Mid*0.85和Moment_Mid組，點擊“確定”按鈕查詢剖面特性。對于開口薄壁截面，可輸入截面剪力，點擊“剪力流”按鈕生成剪應力分布圖（如圖4）。

圖4 開口薄壁截面在單位剪力作用下的剪應力分布圖

值得注意的是，在MSC.Patran中使用封閉單位系統，若采用T-mm-s的單位制，則剖面特性中各參數的單位分別為：慣性矩Iy（mm4），中和軸（mm），剖面模數（mm3），有效面積（mm2）。

表1 結果比較

表1為利用MSC.Patran和Mars 2000中查詢的剖面特性對比，其結果非常接近。

4 結語

通過PCL語言編寫的程序，在MSC.Patran中查看剖面特性的操作非常簡單，有一定應用價值。

［1］MSC.Patran User’s Manual［Z］.

［2］劉烈全，梁樞平編.材料力學［M］.華中理工大學出版社.

［3］章偉星，徐立等.利用PCL語言開發船用型材庫［M］.中國船級社.

［4］任慧龍等.基于PCL的船舶自動加載方法（C）.MSC.Software中國用戶論文集.2002.

Examination on hull section characteristics with PCL language

TANG Xu-dong
（Shanghai Jiao tong University，Shanghai 200030，China）

MSC.Patran；PCL；section characteristics；section moment of inertia

PCL（Patr an Command Language）provides a powerful function for secondary development of MSC.Patran.This paper describes the programs based on PCL for calculating the hull section characteristics，such as the hull section moment of inertia，the section modulus，the neutralization axis and so on.

U663

A

1001-9855（2011）03-0027-03

2011-01-19；

2011-04-12

唐旭東（1980-），男，漢族，工程師，上海交通大學在讀工程碩士，主要從事船舶結構研究設計工作。