How can I obtain timehistory displacement through the SAP2000 OAPI platform?

[maryam-s]

hello.
in this code, i have tried to obtain timehistory displacement out put. according to the Sap OAPI documention, If the user sets the value of 2 in the following function, it will returne timehistory displacement out put.
%get output option
%Value1 = Envelopes; 2 = Step-by-Step; 3 = Last Step)
ret = AnalysisResultsSetup.GetOptionDirectHist(2);
but it seems, there are some problems in my code. because it just results Envelopes displacement
can any one help me?

my code:
%%Record---------------------------------------------------------------------------------------------
%GetRecordNameList
Func= NET.explicitCast(SapModel.Func,'SAP2000v20.cFunction');
NumberNamesF=0;
MyNameF = NET.createArray('System.String',2);
[ret,NumberNamesF, MyNameF] = Func.GetNameList(NumberNamesF, MyNameF,2);

%%LoadCase---------------------------------------------------------------------------------------------
%set data loadcase
LoadCases = NET.explicitCast(SapModel.LoadCases,'SAP2000v20.cLoadCases');
DirHistNonlinear = NET.explicitCast(LoadCases.DirHistNonlinear,'SAP2000v20.cCaseDirectHistoryNonlinear');

%add nonlinear direct history load case
%ret = DirHistNonlinear.SetCase('timehistory');

%%%%%set proportional damping
%get proportional damping
%ret = SapModel.LoadCases.DirHistNonlinear.GetDampProportional("LCASE1", DampType, Dampa, Dampb, Dampf1, Dampf2, Dampd1, Dampd2)
%DampType:
%1 = Mass and stiffness proportional damping by directspecification;
%2 = Mass and stiffness proportional damping by period;
%3 = Mass and stiffness proportional damping by frequency);
%DampType =2;
%ret = DirHistNonlinear.GetDampProportional('timehistory', 0, 0, DampType, Period(1), Period(2), 0.05, 0.05);
%ret = DirHistNonlinear.SetDampProportional('timehistory', 0, 0, DampType, Period(1), Period(2), 0.05, 0.05);
%%%%%set geometric nonlinearity option
%get geometric nonlinearity option
%ret = SapModel.LoadCases.DirHistNonlinear.GetGeometricNonlinearity("LCASE1", NLGeomType)
%NLGeomType0 = None; 1 = P-delta; 2 = P-delta plus large displacements)
NLGeomType = 1;
ret = DirHistNonlinear.SetGeometricNonlinearity('timehistory', NLGeomType);

%%%%%set initial condition
%get initial condition
% ret = SapModel.LoadCases.DirHistNonlinear.GetInitialCase("LCASE1", InitialCase)
%InitialCase:This is blank, None or the name of an existing analysis case. This item specifies if the load case starts from zero initial conditions, that is, an unstressed state, or if it starts from the state at the end of a nonlinear static or nonlinear direct integration time history load case.
ret = DirHistNonlinear.SetInitialCase('timehistory', 'G');

%set load data
%Name:The name of load case.
%NumberLoads:The number of loads assigned to the specified analysis case.
MyNumberLoads=2;
%LoadType:This is an array that includes Load or Accel, indicating the type of each load assigned to the load case.
MyLoadType = {'Accel';'Accel';'Accel'};
%LoadName:This is an array that includes the name of each load assigned to the load case.
%If the LoadType item is Load, this item is the name of a defined load pattern.
%If the LoadType item is Accel, this item is U1, U2, U3, R1, R2 or R3, indicating the direction of the load.
MyLoadName = {'U1';'U2';'U1'};
%Func:This is an array that includes the name of the time history function associated with each load.
MyFunc = {MyNameF(1);MyNameF(1);MyNameF(1)};
%SF:This is an array that includes the scale factor of each load assigned to the load case. [L/s2] for U1 U2 and U3; otherwise unitless
MySF = [500;0;0];
%TF:This is an array that includes the time scale factor of each load assigned to the load case.
MyTF = [1;1;1];
%AT:This is an array that includes the arrival time of each load assigned to the load case.
MyAT = [0;0;0];
%CSys:This is an array that includes the name of the coordinate system associated with each load. If this item is a blank string, the Global coordinate system is assumed.
%This item applies only when the LoadType item is Accel.
MyCSys = {'Global';'Global';'Global'};
%Ang:This is an array that includes the angle between the acceleration local 1 axis and the +X-axis of the coordinate system specified by the CSys item. The rotation is about the Z-axis of the specified coordinate system. [deg]
%This item applies only when the LoadType item is Accel.
MyAng = [0;0;0];

ret = DirHistNonlinear.SetLoads('timehistory', MyNumberLoads, MyLoadType, MyLoadName, MyFunc, MySF, MyTF, MyAT, MyCSys, MyAng);

%%%%%assign a mass source to the static nonlinear load case
ret = DirHistNonlinear.SetMassSource('timehistory','Mass');

%%%%%set solution control parameters
%ret =DirHistNonlinear.SetSolControlParameters('timehistory', 0.01, 0.00001, 15, 50, 0.00005, False, 0.02, 10, 0.2, 1.7);

%%%%%set time integration parameters
%get time integration parameters
%ret = SapModel.LoadCases.DirHistNonlinear.GetTimeIntegration("LCASE1", IntegrationType, Alpha, Beta, Gamma, Theta, m)
%IntegrationType: 1 = Newmark; 2 = Wilson; 3 = Collocation; 4 = Hilber-Hughes-Taylor;5 = Chung and Hulbert
IntegrationType=4;
%Alpha:The alpha factor (-1/3 <= Alpha <= 0)-This item applies only when IntegrationType = 4 or 5
Alpha=0;
%Beta:The beta factor (Beta >= 0)-This item applies only when IntegrationType = 1, 3 or 5. It is returned for informational purposes when IntegrationType = 4
Beta=0.25;
%Gamma:The gamma factor (Gamma >= 0.5)-This item applies only when IntegrationType = 1, 3 or 5. It is returned for informational purposes when IntegrationType = 4.
Gamma=0.5;
%Theta:The theta factor (Theta > 0)-This item applies only when IntegrationType = 2 or 3.
Theta=0;
%m:The alpha-m factor-This item applies only when IntegrationType = 5
m=0;
ret = DirHistNonlinear.SetTimeIntegration('timehistory',IntegrationType , Alpha, Beta, Gamma, Theta, m);

%%%%%get record value
NumberItems=0;
MyTime=NET.createArray('System.Double',2);
Value= NET.createArray('System.Double',2);

[ret, NumberItems, MyTime, Value] = Func.GetValues(MyNameF(1), NumberItems, MyTime, Value);

%%%%%set time step data
%DT:The output time step size.
DT = MyTime(2)-MyTime(1);
%nstep:The number of output time steps.
nstep = (MyTime(NumberItems)/DT)+50;

ret =DirHistNonlinear.SetTimeStep('timehistory', nstep, DT);

%%analysis---------------------------------------------------------------------------------------------
%run analysis
Analyze = NET.explicitCast(SapModel.Analyze,'SAP2000v20.cAnalyze');
ret = Analyze.RunAnalysis;

%%Result------------------------------------------------------------------------------------------------
%SapResult = NET.createArray('System.Double',2);
SapResult = NET.createArray('System.Double',2)
basereactions = NET.createArray('System.Double',2);

%get point
PointName='53';

%% get Sap2000 results for load cases 'timehistory'
AnalysisResults = NET.explicitCast(SapModel.Results,'SAP2000v20.cAnalysisResults');
AnalysisResultsSetup = NET.explicitCast(AnalysisResults.Setup,'SAP2000v20.cAnalysisResultsSetup');

%clear all case and combo output selections
ret = AnalysisResultsSetup.DeselectAllCasesAndCombosForOutput;

%set case and combo output selections
ret = AnalysisResultsSetup.SetCaseSelectedForOutput('timehistory');

%get output option
%Value1 = Envelopes; 2 = Step-by-Step; 3 = Last Step)
Value = 2;

ret = AnalysisResultsSetup.GetOptionDirectHist(2);
% ret = SapModel.Results.Setup.SetOptionModalHist(2)
ret = AnalysisResultsSetup.SetOptionModalHist(2);
NumberResults =0;

Obj = NET.createArray('System.String',2);
Elm = NET.createArray('System.String',2);
ACase = NET.createArray('System.String',2);
StepType = NET.createArray('System.String',2);
StepNum = NET.createArray('System.Double',2);
U1 = NET.createArray('System.Double',2);
U2 = NET.createArray('System.Double',2);
U3 = NET.createArray('System.Double',2);
R1 = NET.createArray('System.Double',2);
R2 = NET.createArray('System.Double',2);
R3 = NET.createArray('System.Double',2);

[ret, NumberResults, Obj, Elm, ACase, StepType, StepNum, U1, U2, U3, R1, R2, R3] = AnalysisResults.JointDispl('53', SAP2000v20.eItemTypeElm.ObjectElm, NumberResults, Obj, Elm, ACase, StepType, StepNum, U1, U2, U3, R1, R2, R3);

SapResult = U1(1)