前項でスカラー関数の定義とその使用方法について説明しました。では、定義したスカラー関数がちゃんと意図したものになっているかどうかを確認するにはどうしたら良いでしょうか。ここでは、その一つの方法を紹介します。ここで紹介する“abstract node”と“genel:clamp”を用いる方法は、スカラー関数の可視化だけでなく、モデル中の任意の変数を可視化する場合にも適用することができます。
“abstract node”は、１つの自由度を持つ汎用のnodeです。nodesブロックの中で次のように定義します。
abstract: <node label>, algebraic, value, <initial value>;
“genel” elementに属する “clamp”は、任意のnode自由度にdriveで指定する値を強制的に与えます。elementsブロックの中で次のように定義します。
genel: <label>,
clamp,
<clamped node>,
<imposed value>;
ここで、<imposed value>は強制的に与える値で、drive callerによって定めます。
定義したスカラー関数 f (x) を -5 ≤ x ≤ 5 の範囲で可視化したい場合、時間0から10までのMBDyn解析（initial value problem）を行い、その間に x = Time - 5 と y = f (x) を各時間ステップで評価して出力すれば求めるデータが得られます。変数 x と y を出力するために、それぞれの変数の値を格納するabstract nodeを定義し、それらのabstract nodeにgenel:crampで値を与えます。次のコード１は、スカラー関数“myfunc2”を可視化するための入力ファイルの例です。
# scalarfunc.mbd
#-----------------------------------------------------------------------------
# Scalar Function
scalar function: "myfunc2",
cubicspline, do not extrapolate,
-2.0, 0.0,
-1.0, 0.0,
0.0, 0.0,
1.0, 1.0,
2.0, 1.0,
3.0, 1.0;
# X Range
set: real XL = -5.; # Lower Limit
set: real XU = 5.; # Upper Limit
set: real dX = 0.01; # Step
#-----------------------------------------------------------------------------
# [Data Block]
begin: data;
problem: initial value;
end: data;
#-----------------------------------------------------------------------------
# [<Problem> Block]
begin: initial value;
initial time: 0.;
final time: XU-XL;
time step: dX;
max iterations: 10;
tolerance: 1.e-6;
end: initial value;
#-----------------------------------------------------------------------------
# [Control Data Block]
begin: control data;
abstract nodes: 2;
genels: 2;
end: control data;
#-----------------------------------------------------------------------------
# Abstract Node Labels
set: integer NoAbs_X = 1;
set: integer NoAbs_Y = 2;
# Genel Labels
set: integer GeClamp_NoAbs_X = 1;
set: integer GeClamp_NoAbs_Y = 2;
#-----------------------------------------------------------------------------
# [Nodes Block]
begin: nodes;
abstract: NoAbs_X, algebraic,
value, XL;
abstract: NoAbs_Y, algebraic,
value, model::sf::myfunc2(XL);
end: nodes;
#-----------------------------------------------------------------------------
# Plugin Variable
set: [dof, X, NoAbs_X, abstract, algebraic];
#-----------------------------------------------------------------------------
# [Elements Block]
begin: elements;
genel: GeClamp_NoAbs_X,
clamp,
NoAbs_X, abstract,
string, "Time+XL";
genel: GeClamp_NoAbs_Y,
clamp,
NoAbs_Y, abstract,
string, "model::sf::myfunc2(X)";
end: elements;
コード１を入力ファイルとしてMBDynを実行すると、出力ファイルとして「scalarfunc.abs」が生成されます。このファイルにabstract nodeのデータが出力されます。absファイルのデータは２列から成り、１列目がnodeラベル、２列目がnodeの値となっています。コード２は、「scalarfunc.abs」のデータを読み込んで、スカラー関数“myfunc2”をプロットするMATLABスクリプトです（コード中の「MBDynLoad」関数については第11項参照）。このスクリプトによって図１のようなグラフが得られます。
% plot_scalarfunc.m
clear; close all;
[LABEL,DATA] = MBDynLoad('scalarfunc.abs');
X = DATA(:,2,LABEL==1);
Y = DATA(:,2,LABEL==2);
Xd = [-2, -1, 0, 1, 2, 3];
Yd = [0, 0, 0, 1, 1, 1];
figure(1)
plot(X,Y,'-b',Xd,Yd,'ob','LineWidth',1)
grid on;
set(gca,'FontSize',16);
xlabel('x','FontSize',16);
ylabel('f(x)','FontSize',16);
xlim([-3,4])
ylim([-0.2,1.4])
