SDTools code for bending beam error

Siwar Ben Hadj Ali

2022 4 22

0 个回答

2 次查看（30 天）

0 个投票

Please can you help me , I am trying to model bending beam in SDToolbox and an unrecognized function or variable error appear

% Mesh -------------------------------------------------------------
mdl0=struct('Node',[],'Elt',[],'name','primer','unit','SI'); % init model
% initial node declaration
mdl0.Node=[1 0 0 0 0 0 0 ;
2 0 0 0 0 1 0 ;
3 0 0 0 1 1 0 ;
4 0 0 0 1 0 0 ]; % define nodes
% define 4 beam1 elements
prop=[1 1 0 1 0]; % MatId ProId nx ny nz
mdl0.Elt=feutil('ObjectBeamLine 1 2 0 2 3 0 3 4 1',prop);
feutilb('write',mdl0) % display how to write this model
mdl0=feutil('RepeatSel 10 1 0 0',mdl0);
model=mdl0;
cf=feplot(model); % plot model
fecom('ColorDataGroup');
% Material and element properties -------------------------------------
model.pl=[1 fe_mat('m_elastic','SI',1) 7.2e9 .3 2700 2.68e10]; % material
model.il=[1 fe_mat('p_beam','SI',1) 1e-6 1e-7 1e-7 1e-4 0.8];  % beam1 properties
% force MatId and ProId of both groups to be 1
model.Elt=feutil('SetGroup 1 mat1 pro1',model);
model.Elt=feutil('SetGroup 2 mat1 pro1',model);
% Loads and boundary conditions ----------------------------------------
load=struct('DOF',[1.02, 2.02, 3.02, 4.02, 5.02, 6.02, 7.02, 8.02, 9.02, 10.02,],'def',[5]); % Load
model=fe_case(model,'DofLoad','endload',load) ; % apply load 
model=fe_case(model,'FixDof','LeftEdge', 'x==0 -DOF 1 2 3 4 5 6'); % BC
% Solving
% Normal modes ---------------------------------------------------------
cf=feplot(model); model=cf.mdl;
cf.Stack{'info','EigOpt'}=[5 10 1e3];
% deformation calculated for model
% options: 5 - Lanczos solver
% 10 - number of modes
% 1e3 - shift (needed with rigid body modes)
d1=fe_simul('Mode',model); cf.def=d1;
% Time simulation and signal processing --------------------------------
cf.Stack{'curve','q0'}=[];
% initial conditions for time response set to zero
% define options
cf.Stack{'info','TimeOpt'}=...
struct('Method','newmark', ...
'Opt',[.25 .5 0 1e-3 10000], ... % beta,gamma,t0,dt,N
'NeedUVA',[1 1 0]); % compute disp, vel, not acc
% Generate a step over the first ten time steps with the load
cf.CStack{'endload'}.curve={fe_curve('TestStep t1=1e-2')};
% Launch the simulation
def=fe_simul('Time',model);
cf.def=def;
fecom(';colordataEvalZ;view3;viewh+180');
fecom('AnimTime2'); % animate displacement, using 1 every 2 time steps
def.name='Time';ci=iiplot(def); % Do some signal processing
ii_mmif('fft',ci,'Time');
iicom(ci,'iix:fft(Time)Only')