This is that what I want! please share it
And what is the simulator program's name?
Here is a zip of all the files needed (I think) to run the sim.
The piezo actuator subcircuit is called quartz_crystal because piezo crystals have the same circuit topology, just different component values.
The BU406bd is a subcircuit consisting of a BU406 transistor with an ideal 6V zener in parallel with the base, to represent base-emitter breakdown. This zener is modeled as an ideal diode with a forward voltage of 6V, and 10Ω series resistance. I didn't use a "real" zener model because this would add an actual diode in parallel with the base-emitter junction, with its attendant capacitance, etc. Common transistor models do not include base-emitter reverse breakdown.
You can descend into the schematics of these two parts by right-clicking on the symbol. This will bring up a dialog box. One of the options is
Open Schematic.
Note that, without the base-emitter breakdown being modeled, the simulation results in a 340kHz sinusoidal oscillation at about 8V p-p.
Note also that breakdown of the base-emitter junction causes beta degradation that is proportional to breakdown current, and is cumulative with time.
Here is a thread that contains results of measurements made by one Fred Bartoli that confirms this.
Arvinfx, if you haven't used
LTspice, be prepared for a steep learning curve.