Why voltage output is equal to zero? Answered

This is my model. I input force at the bottom of the steel block. I can get the displacement at the point in figure. But the  output voltage is equal to zero. I don't know why.

c mem 800 200 /* Allocate 800 megawords of memory - 3 GigaBytes (Not necesarry for Windows Operating Systems)
c NOTE: MEM Command must be first command in file, if used. (Line 1)
rest no 

c ************************************************************************************************************* 
c 
c                 Generated Flex Input File
c 
c ************************************************************************************************************* 
c 
c    DESIGNER            :OnScale - Designer Generated
c    MODEL DESCRIPTION    :
c    DATE CREATED        :26 9ÔÂ 2022
c    VERSION                :1.0
c ************************************************************************************************************* 

mp
    omp * *           /* Number of CPUs to be used in the execution.
    end
 
titl calibration2 

c ************************************************************************************************************* 
c 
c             Define User Variables 
c 
c ************************************************************************************************************* 
c 
c These variables have been set by the user through the interface.
c 
c ************************************************************************************************************* 

symb coordFactor = 1.0         /* Coordinate conversion factor
symb timeFactor = 1.0         /* Time conversion factor
symb dMassFactor = 1.0         /* Mass conversion factor
symb freqint = 1e+06         /* Determined Drive Frequency (Hz)

c ************************************************************************************************************* 
c 
c             Define Meshing 
c 
c ************************************************************************************************************* 
c 
c             Set the variable for the approximate element size for the model. Must be 
c            sufficient to represent the wavelengths of interest. Recommended that at least 
c            15 elements per wavelength are used.

c 
c ************************************************************************************************************* 

symb freqdamp = $freqint
symb box = 0.00025

c ************************************************************************************************************* 
c 
c             Geometry Locations (XYZ)
c 
c ************************************************************************************************************* 

c Scale Parameters
symb xmin = 0.0 * $coordFactor
symb xmax = 0.25 * $coordFactor
symb ymin = 0.0 * $coordFactor
symb ymax = 0.12502 * $coordFactor

c Determine lengths of the model
symb xlen = ( $xmax - $xmin )
symb ylen = ( $ymax - $ymin )

c *************************************************** 
c 
c             Keypoints in the X-Direction
c 
c *************************************************** 

symb #keycord x 1 $xmin $xlen 
symb #get { idx } rootmax x 

c *************************************************** 
c 
c             Keypoints in the Y-Direction
c 
c *************************************************** 

symb #keycord y 1 $ymin $ylen 
symb #get { jdx } rootmax y 

c ************************************************************************************************************* 
c 
c             Indices Locations (IJK)
c 
c ************************************************************************************************************* 

c Grid in I direction, using approximately element size of 'box' and at least 1 element
symb #keyindx i 1 $idx 1 $box 1 
symb indgrd = $i$idx

c Grid in J direction, using approximately element size of 'box' and at least 1 element
symb #keyindx j 1 $jdx 1 $box 1 
symb jndgrd = $j$jdx

c ************************************************************************************************************* 
c 
c             GCON Grid & Geometry Allocation
c 
c ************************************************************************************************************* 

grid $i$idx $j$jdx axiy 
geom
    keypnt $idx $jdx 
    end
 

c ************************************************************************************************************* 
c 
c             Driving Conditions 
c 
c ************************************************************************************************************* 
c 
c             A number of predefined waveform functions can be accessed in OnScale. The DATA 
c            HIST option is used below, other examples include wavelets, gaussians and step 
c            functions. The manual details the function entries

c 
c ************************************************************************************************************* 

symb freqtimefunc_1 = 1.0
data
    hist ud_timefunc_1 * 'E:/softfile/onscale/calibration/time_force.dat' 
    end
 
func
    name timefunc_1 
    hist ud_timefunc_1 * 
    end
 
c -------------------------------------------------------------- 
c                   Project Material List
c -------------------------------------------------------------- 

c -------------------------------------------------------------------------
c Global variables used in all the material definitions
c -------------------------------------------------------------------------
c
symb epvacm = 8.854e-12               /* dielectric constant for vacumn
symb freqdamp = 1.e6 if noexist       /* specified frequency for damping model
symb rmu0 = 1.2566e-6

symb #msg 5
********************************************************
Damping matched at $freqdamp Hz
Redefine variable 'freqdamp' if device centre frequency
varies significantly from this value
********************************************************

c
c -------------------------------------------------------------------------
c    Now define the axis transformation - only posx used in this file
c -------------------------------------------------------------------------

axis
     form vctr
     defn posx car2 0. 0. 0.    1.  0.  0.   0. 1. 0.
     defn negx car2 0. 0. 0.   -1.  0.  0.   0. 1. 0.
     defn posy car2 0. 0. 0.    0.  1.  0.   0. 0. 1.
     defn negy car2 0. 0. 0.    0. -1.  0.   0. 0. 1.
     defn posz car2 0. 0. 0.    0.  0.  1.   1. 0. 0.
     defn negz car2 0. 0. 0.    0.  0. -1.   1. 0. 0.
     end

c -------------------------------------------------------------------------
c    Input material properties to program
c -------------------------------------------------------------------------

matr

c -------------------------------------------------------------- 
c type : METAL : 
c name : steel : 
c desc : Mild Steel , generic : 
c -------------------------------------------------------------- 

    wvsp on 
    type elas 
    prop steel 7900.00 5900.00 3200.00 0.010000
    vdmp steel $freqdamp db 0.3 1.2 1e+06 1 0.01 1
    thrm steel 472 49.8 49.8 49.8 1.0 0 0 27

c -------------------------------------------------------------- 
c type : METAL : 
c name : matching : 
c desc : metal : 
c -------------------------------------------------------------- 

    wvsp on 
    type elas 
    prop matching 3814.00 6818.00 4091.00 0.010000

c -------------------------------------------------------------- 
c type : PIEZO : 
c name : pzt5h : 
c desc : PZT5H Generic : 
c -------------------------------------------------------------- 

    symb rho  = 7500.00                    /* density
    wvsp off 
    type lean

c    define baseline stiffness coefficients (Constant electric field)

    symb c11 = 1.272e+11
    symb c12 = 8.02e+10
    symb c13 = 8.47e+10
    symb c14 = 0
    symb c15 = 0
    symb c16 = 0
    symb c22 = 1.272e+11
    symb c23 = 8.47e+10
    symb c24 = 0
    symb c25 = 0
    symb c26 = 0
    symb c33 = 1.174e+11
    symb c34 = 0
    symb c35 = 0
    symb c36 = 0
    symb c44 = 2.3e+10
    symb c45 = 0
    symb c46 = 0
    symb c55 = 2.3e+10
    symb c56 = 0
    symb c66 = 2.35e+10

    prop pzt5h $rho 
    $c11     $c12     $c13     $c14    $c15    $c16    $c22
    $c23     $c24     $c25     $c26    $c33    $c34    $c35
    $c36     $c44     $c45     $c46    $c55    $c56    $c66

c    define baseline dielectric coefficients

    symb epxx = 1700
    symb epyy = 1700
    symb epzz = 1470

c    scale material properties as specified above
    symb aepxx = $epvacm * $epxx
    symb aepyy = $epvacm * $epyy
    symb aepzz = $epvacm * $epzz

    elec pzt5h $aepxx $aepyy $aepzz

c    define baseline piezoelectric coupling coefficients

    symb ex1 = 0            /* coupling constant
    symb ex2 = 0            /* coupling constant
    symb ex3 = 0            /* coupling constant
    symb ex4 = 0            /* coupling constant
    symb ex5 = 17            /* coupling constant
    symb ex6 = 0            /* coupling constant
    symb ey1 = 0            /* coupling constant
    symb ey2 = 0            /* coupling constant
    symb ey3 = 0            /* coupling constant
    symb ey4 = 17            /* coupling constant
    symb ey5 = 0            /* coupling constant
    symb ey6 = 0            /* coupling constant
    symb ez1 = -6.5            /* coupling constant
    symb ez2 = -6.5            /* coupling constant
    symb ez3 = 23.3            /* coupling constant
    symb ez4 = 0            /* coupling constant
    symb ez5 = 0            /* coupling constant
    symb ez6 = 0            /* coupling constant

    piez pzt5h 1 1 $ex1 1 2 $ex2 1 3 $ex3 1 4 $ex4 1 5 $ex5 1 6 $ex6 &
            2 1 $ey1 2 2 $ey2 2 3 $ey3 2 4 $ey4 2 5 $ey5 2 6 $ey6 &
            3 1 $ez1 3 2 $ez2 3 3 $ez3 3 4 $ez4 3 5 $ez5 3 6 $ez6

    rdmp pzt5h $freqdamp q 65 65 1e+06 * * *
    axis pzt5h posz   /* relate materials local system to global system

c -------------------------------------------------------------- 
c type : METAL : 
c name : alum : 
c desc : Aluminium (damping v. low) : 
c -------------------------------------------------------------- 

    wvsp on 
    type elas 
    prop alum 2690.00 6306.00 3114.00 0.010000
c    define baseline dielectric coefficients

    symb epxx = 3.1
    symb epyy = 3.1
    symb epzz = 3.1

c    scale material properties as specified above
    symb aepxx = $epvacm * $epxx
    symb aepyy = $epvacm * $epyy
    symb aepzz = $epvacm * $epzz

    elec alum $aepxx $aepyy $aepzz

    elec void $epvacm
    end

c ************************************************************************************************************* 
c 
c             Primitive Definitions
c 
c ************************************************************************************************************* 

site
    regn void 
    blok steel 0.0 0.25 0.0 0.1 0.0 0.0 
    blok matching 0.0 0.01005 0.1 0.10116 0.0 0.0 
    blok pzt5h 0.0 0.00674 0.10116 0.11616 0.0 0.0 
    multipoly alum stnd * * * * * * 0.0 0.0 1 
    0.009 0.10116 0.009 0.11502 0.008 0.11602 0.008 0.12402 0.0 0.12402 0.0 0.12502 0.011 0.12502 0.011 0.10056 0.01005 0.10056 0.01005 0.10116 
    end
 

c ************************************************************************************************************* 
c 
c             Boundary Definitions
c 
c ************************************************************************************************************* 

boun
    side xmin symm 
    side xmax free 
    side ymin free 
    side ymax free 
    end
 

c ************************************************************************************************************* 
c 
c             Calculated Properties 
c 
c ************************************************************************************************************* 
c 
c             By default, Flex only calculates the minimum required data set, typically this 
c            means only velocities. This is done for memory efficiency. Should other 
c            properites be required (e.g. displacements, stresses, strains, pressure), then 
c            these must be requested by the CALC command. The manual lists all these options

c 
c ************************************************************************************************************* 

calc
    disp y          /* calculate displacements
    end
 

c ************************************************************************************************************* 
c 
c             Piezoelectric Load Definitions
c 
c ************************************************************************************************************* 

piez
    wndo auto piez 
    defn load_1 
    symb bxmin = 0.0 * $coordFactor
    symb bymin = 0.11616 * $coordFactor
    symb #get { is js * } clsnode $bxmin $bymin *
    symb bxmax = 0.00674 * $coordFactor
    symb bymax = 0.11616 * $coordFactor
    symb #get { ie je * } clsnode $bxmax $bymax *
    nod2 pzt5h void $is $ie $js $je * * 
    bc load_1 off 

    defn load_2 
    nod2 pzt5h matching 
    bc load_2 grnd 
    slvr drct 
    end
 

c ************************************************************************************************************* 
c 
c             Force/Velocity Load Definitions
c 
c ************************************************************************************************************* 

rigd
    defn load_3 
    symb bxmin = 0.0 * $coordFactor - $box
    symb bymin = 0.0 * $coordFactor - $box
    symb #get { is js * } clsnode $bxmin $bymin *
    symb bxmax = 0.0005 * $coordFactor + $box
    symb bymax = 0.0 * $coordFactor + $box
    symb #get { ie je * } clsnode $bxmax $bymax *
    node $is $ie $js $je * * steel 
    bc forc functimefunc_1 0.0 1.0 0.0 * * * 
    end
 

time * * 0.8 

c ************************************************************************************************************* 
c 
c             Process Model 
c 
c ************************************************************************************************************* 
c 
c             Issue process (PRCS) command. Checks model integrity, and calculates stable 
c            time step. NOTE: Process command must always be issued

c 
c ************************************************************************************************************* 

symb #msg 1
Checking Model Integrity......
prcs
grph
    colr user size 4 
    colr user 1 1.0 1.0 0.0 
    colr user 2 0.0 1.0 1.0 
    colr user 3 0.0 1.0 0.0 
    colr user 4 0.0 0.0 1.0 
    colr tabl matr 7 
    map steel 1 
    map matching 2 
    map pzt5h 3 
    map alum 4 
    end
 
grph
    line off 
    nvew 2 1 
    arrow pole 
    plot piez 
    plot matr piez 
    arrow off 
    end
 
term

c ************************************************************************************************************* 
c 
c             Choose Time Histories To Store 
c 
c ************************************************************************************************************* 
c 
c             Save field (such as displacement or pressure) from a node or element for all 
c            time steps with POUT command. Histories are referenced by order of 
c            specification. Histories will be saved in the Flex History file (flxhst).

c 
c ************************************************************************************************************* 

pout
    rate 1 
    symb xloc = 0.0 * $coordFactor
    symb yloc = 0.1 * $coordFactor
    hist xyz ydsp $xloc * * $yloc * * 
    hist functimefunc_1 
    histname electrode vq all 
    end
 

c ************************************************************************************************************* 
c 
c             Run the Model 
c 
c ************************************************************************************************************* 
c 
c             Specify the number of time steps to be run. Can be set to auto by using 
c            'Ringdown'.

c 
c ************************************************************************************************************* 

c User defined runtime for the model
symb #get { step } timestep 
symb ttime = 0.0003
symb nexec = $ttime / $step
exec $nexec 

symb n_extra = nint ( $step / 10 )
exec $n_extra 
term

data
    file out 'calibration2.flxdato' 
    out modl 
    end
 

c ************************************************************************************************************* 
c 
c             Save symbol variables to file for later use
c 
c ************************************************************************************************************* 

symb #get { labl } jobname 
symb #save '$labl.symb' 
stop           /* return to command prompt