Contents     Functions              PDF Index

d_piezo

• Script
• MeshPlate
• Mesh

Purpose

Support function for piezoelectric demos

Syntax

 sdtweb('_taglist','d_piezo')  % display contents

Accepted commands are

Script

These commands group sample scripts. Use d_piezo to display tag list and see available contents.

MeshPlate

Meshing utilities for the placement of piezoelectric patches on a supporting structure (flat plate for now).

The options are specified in a structure with fields

• .list : defines a list of features to be introduced with columns giving name,LamSpec,Geo, name laminate specification and shape options.
• .unit : gives the model unit (needed since patch dimensions are always given in mm).

The laminate specification string is composed of the following

• BaseIdi gives the ProId of the base laminate which is then used to figure out the position of patches.
• +Patch or -Patch to place a patch above or below the base laminate.
• Patch itself is a specification of a patch material and geometry. The list of implemented patch can be obtained using m_piezo Patch
• .In to specify that the patch has an enforced voltage.

The geometry/position specification string Geo can be

• a specification of the patch corner and orientation such as xc=.03 yc=.05 ang=30 if the patch geometry is specified using the laminate specification.
• rect shapes [xc lx nx yc ly ny alpha MatId ProId] where MatId and ProId are filled automatically if not provided.
• circ shapes [xc yc rc lc (MatId ProId)]
• global lx=.4 ly=.3 lc=.02 -Sens. The option -Sens generates a sensor entry corresponding to normal displacement of the initial mesh. Alternatively you can add a sensor configuration SensDOF entry see sdtweb('sensor#scell') or sdtweb('sensor#sstruct').

% Start by defining properties of the underlying laminate
mdl=struct('Node',[],'Elt',[], ... % empty model
    'pl', ... % composite layer property
      [1 fe_mat('m_elastic','SI',1) 42.5e9 .042 1490 3.35e9 .01], ... 
     'il', ... % laminate definition (6 layers at 0,90,0,90,0,90)
     p_shell(['dbval 1 laminate    1 2.167e-4 0    1 2.167e-4 90 ' ...
         '1 2.167e-4 0 1 2.167e-4 90 1 2.167e-4 0 1 2.167e-4 90']), ...
     'unit','SI'); 
RG=struct;
RG.list={'Name','Lam','shape'
   'Main_plate', mdl,'global lx=.4 ly=.3 lc=.02'
   'Act1','BaseId1 +SmartM.MFC-P1.2814 -SmartM.MFC-P1.2814.in','xc=.35 yc=.25 ang=30'
   'Sen2','BaseId1 +SmartM.MFC-P1.2814','xc=.03 yc=.05 ang=30'
   'Sen3','BaseId1 +Noliac.NCE51.OD25TH1','xc=.05 yc=.25'
   };
cf=feplot;d_piezo('MeshPlate',RG);cf.mdl.name='Plate with piezo';
p_piezo('electrodeinfo',cf.mdl.GetData)
matgui('jil',cf);matgui('jpl',cf); % Display properties

The following illustrates transient simulation to a load on a specific piezo

model=cf.mdl.GetData; model=p_piezo('electrode2case',model);
opt=fe_time('TimeOpt Newmark .25 .5 0 .3e-6 200'); 
opt.AssembleCall='assemble -fetime Load';
%opt.FinalCleanupFcn='out.DOF=model.DOF;out.def=[Case.T*out.def+Case.TIn*[0 ft(1:end-1)'']];';
model=stack_set(model,'info','Rayleigh',[0 2*.0025/200e3]);
def=fe_time(opt,cf.mdl);def.name=model.name;
cf.def=def; fecom('colordataevalRadZ-edgealpha0');fecom('scc1e-10');

Mesh

Patch Simple volume patch.

Plate Generic script for arbitrary placement of patches on a flat plate. A list of shapes can be given as a cell array. This is considered as a demo since it currently only supports a rectangular base plate.

GammaS build a weighting for surface control.