---

Hi,

ich arbeite im Moment an einer harmonischen Analyse mit piezoelektrischem Werkstoff und piezoelektrischem ciru94 Element.
Klappt von der Auslenkung des harmonisch angeregten Körpers her gut.
An dem Widerstand möchte ich nun über

esol,13,e_resistor,,smisc,1,VOLT_Resistor
esol,14,e_resistor,,smisc,2,CURRENT_Resistor
esol,15,e_resistor,,smisc,1,POWER_Resistor
plvar,13,14,15

allerhand ausgeben lassen. Es kommt aber immer die Fehlermeldung "No data points saved this storage step"
die Spannungen an den Knoten des Widerstands kann ich mir über nsol ausgeben lassen.
Schätze ich mache etwas bei den Solverinstellungen falsch, so das nur Knotenergebnisse gespeichert werden.
Habe es schon mit 'outres,elem,all' probiert, aber es kommt immer noch die Fehlermeldung.

 

Im Voraus Danke für alle Antworten.

MfG

Simulinho

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo,

1)  you should select element by using of: ESEL. Thus e_resistor should be the variable with the exactly element number in you task.
2)for the presentation of power: NMISC (not smisc).
3)you can get the results of current by using of:
ESOL,2,elment_number,its_node,AMPS,,variable_name
4) Try graphical picking, it should be not the solver problem.

Bye
Serg.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hi,

1) e_resistor is the element-number of the resistance element

2) my mistake, I saw it after this posting

3) I will try

4) In my opinion, too

Thanks for your answer. But now I have another problem with a simlilar analysis.

Just another geometry, more simply to get faster results or errors  . But I get an error message: "Error in numeric factorization.The system matrix, possibly, is numerically singular."     

What's the problem ?

MfG

Simulinho

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Recently I have been solving the 3-d pulse transformer and also had a lot of problems with the convergence. Only by 1e6 elements and 40 calculation hours were OK. The mesh should be fine!
Bye,
Serg.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Check the coupling (CP) such error is it's diagnose, sure

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hi,

I can do what I want.

"Error in numeric factorization.The system matrix, possibly is numerically singular."

I tried everything. Coupling, no coupling. Open, short, it's always the same.

Sometimes I think there is a little gremlin in my PC, who is laughing about me.

MfG

Simulinho

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Could I see the code?

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hi,
at first I make a modal analysis to get resonance frequency.

The bottom electrode is grounded and the top electrode is coupled.

Now I build resistor using circuit element circu94. One node grounded and the other in the coupled set.

/com,*********************************************
/com,*******  PREPROCESSING 2 - HARMONIC      ***
/com,*********************************************

/prep7

n,1,laenge/3,breite/2,hoehe*6
n,2,laenge/3,breite/2,-hoehe*7
n_bot=2
n_top=1
real,2
type,2
e,n_bot,n_top

/com,********************************************
/com,*********  SOLUTION 2 - harmonic  *********
/com,********************************************

/solu
antype,harmic
hropt,full
hrout,off
outpr,nsol
outres,all,all
harfrq,frequency(harm_nmode)*0.9,frequency(harm_nmode)*1.1
nsubst,10

/com,...........................damping..........................
/com,                         
 
mdamp,1,d ! d= 0.01

/com,...........................Randbedingungen

cmsel,s,n_einsp
d,all,ux,0
d,all,uy,0
d,all,uz,0

nsel,s,loc,z,hoehe,hoehe*10
cp,1,volt,all

nsel,s,loc,z,0,-hoehe*10
d,all,volt,0

nsel,s,loc,x,laenge
f,all,fz,-0.001

alls
/pbc,u,,1
eplo
solve

/com,********************************************
/com,*******  POSTPROCESSING 2 -harmonic  ******
/com,********************************************

/post26
numvar,200
/com,

/com,...........................misc.............................
/com,                          2 Kontroll-Knoten finden

/erase
*dim,nt,array,2
nt(1)=node(0.025,0.003,-0.0021 )
nt(2)=node(0.025,0.003,0.0027)
nsol,11,nt(1),u,z,U_1
nsol,12,nt(2),u,z,U_2
nsel,all
/grid,1
/axlab,x,Frequenz [Hz]
/xrange,280,300
/axlab,y,Verschiebung in z-Richtung [m]
plvar,11,12

/erase
/win,1,top
esol,13,e_resistor,,smisc,1,VOLT_Resistor
esol,14,e_resistor,,smisc,2,CURRENT_Resistor
esol,15,e_resistor,,nmisc,1,POWER_Resistor
store
/grid,1
/axlab,x,Frequenz [Hz]
/xrange,280,300
/axlab,y,allerlei [x]

plvar,13,14,15

/noerase
/win,1,bot
nsol,16,n_top,volt,VOLT_Resistor_top
nsol,17,n_bot,volt,VOLT_Resistor_bot
/grid,1
/axlab,x,Frequenz [Hz]
/xrange,280,300
/axlab,y,Spannung [V]

plvar,16,17

/erase
/win,1,full

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Are you using solid5 for the piesoelement? and generally It's quite difficult to invent your preprocessor code. Could you provide the whole code, possible it's realy a devil inside your comp

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

circuit_harmonic.txt

Hi,
sorry, I had a problem with the upload. So I just posted a part of my file.
The harmonic analysis without the circu-element is running with good results using modal superposition.
For an analysis using circu94 I need a full harmonic analysis.

MfG

Simulinho

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

1) If you are using a resistor you should at least connect it, like:

nsel,s,loc,z,0
d,all,volt,spannung_unten
cp,1,volt,all
*GET,NBOT,NODE,0,NUM,MIN      ! GET MASTER NODE ON BOTTOM ELECTRODE
nsel,s,loc,z,hoehe
cp,2,volt,all
*GET,NTOP,NODE,0,NUM,MIN      ! GET MASTER NODE ON BOTTOM ELECTRODE
real,2
type,2
e,NBOT, NTOP

2) You’ve confused CP (coupled degrees of freedom) and D (constrain) by the same labels what it couldn’t be done.
Delete at the harmonic analysis the constrains and the task will be calculated, or complicate electric circuit with the constrain at the another point.

nsel,s,loc,z,0
ddele,all,volt        ! Delete voltage constraints

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hi,

That's the way I made it.

nsel,s,loc,z,0
d,all,volt,spannung_unten

nsel,s,loc,z,hoehe
cp,1,volt,all

I think, the lines like you posted:

nsel,s,loc,z,0
d,all,volt,spannung_unten
cp,1,volt,all

It would not run, because what you wrote under 2).

But now it' running. I think it was a mistake in meshing, I had 4 elements in thickness.
Now I have one. It's running if I do it like this:

block,0,laenge,0,breite,0,hoehe
!lsel,s,tan1,z,1
!lsel,a,tan1,z,-1
!lesize,all,hoehe/4

!lsel,s,tan1,x,1
!lsel,a,tan1,x,-1
lesize,all,laenge/30

!lsel,s,tan1,y,1
!lsel,a,tan1,y,-1
!lesize,all,breite/10

MfG

Simulinho

[Diese Nachricht wurde von Simulinho am 21. Jun. 2007 editiert.]

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP