Hello All,

I am trying to extract the damping matrix for my model, so i start by
building a simple cantilever, and I tried two solutions:
1- Modal analysis
2- Transient Analysis.

I found that people usually use hbmat commands for that, but it so not
working for me in both analysis.

here is the code I am using in case someone can catch the mistake:

fini
/clear
/prep7
ET,1,SOLID185

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,260e6
MPDATA,PRXY,1,,0.25
MPDATA,DAMP,1,,0.1
MPDATA,DENS,1,,1300
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,270e6
MPDATA,PRXY,2,,0.20
MPDATA,DAMP,2,,0.08
MPDATA,DENS,2,,1500
!*
!---- Modeling and meshing
block,0,0.25,0,0.5,0,0.25
TYPE,1
MAT,1
LESIZE,1, , ,2, ,1, , ,0
LESIZE,2, , ,1, ,1, , ,0
LESIZE,9, , ,1, ,1, , ,0
vmesh,all
Mat,2
extopt,esize,2,0 ,0
asel,s,area,,4
vext,all,,,,0.5
!-----BCs
allsel
nummrg,node
numcmp,node
numcmp,elem
Allsel
IC,all,all,0,0
nsel,s,loc,Y,0
D,all,UX,0
D,all,UY,0
D,all,UZ,0
Allsel
!---- Load
F,13,FZ,-1000
Allsel
fini !
!----- Solution ----------------------!
!----------------TRANSIENT -------------!
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,0
NSUBST,1,1,1
TIME,0.1
/solu
SOLVE
fini
!------------  MODAL --------------------!
!/solu
!allsel
!outres,all,all
!ANTYPE,MODAL
!MODOPT,QRDAMP,10
!MXPAND,10
!SOLVE
!fini
!----------------Printing files ---------!
/AUX2
file,file,full
hbmat,HBMstiff,txt,,ascii,stiff,yes,yes
hbmat,HBMmass,txt,,ascii,mass,yes,yes
hbmat,HBMdamp,txt,,ascii,damp,yes,yes

FINISH

Thank you

Ahmed El-Ghandour

PhD Candidate at UIC.

--

Ahmed I. El-Ghandour

Research/Teaching Assistant & Ph.D. Candidate
Dept. of Civil & Materials Engineering
College of Engineering
University of Illinois at Chicago (UIC)
Mobile  : +1 312-912-2617
E-mail  : aelgha2@uic.edu
Skype    : a.i.ghando

Hi Ahmed,

The code to access the various matrices works as written and has produced all three matrices for the modal analysis.  Your code below needs some comments and a fini removed but otherwise gives a result.  Are you looking in the correct location for the file?  By the way, the result of this command is different for the two kinds of analysis.  See the manual for the various caveats.

Jim

James J. Grudzinski, Ph. D., P.E.
Manager, Engineering Operations
Argonne National Laboratory

-----Original Message-----
From: Xansys-temp [mailto:xansys-temp-bounces@xansystest.info] On Behalf Of Ahmed I. El-Ghandour
Sent: Tuesday, November 29, 2016 5:31 PM
To: xansys-temp@xansystest.info
Subject: [Xansys] Extracting Damping matrix

Hello All,

I am trying to extract the damping matrix for my model, so i start by building a simple cantilever, and I tried two solutions:
1- Modal analysis
2- Transient Analysis.

I found that people usually use hbmat commands for that, but it so not working for me in both analysis.

here is the code I am using in case someone can catch the mistake:

fini
/clear
/prep7
ET,1,SOLID185

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,260e6
MPDATA,PRXY,1,,0.25
MPDATA,DAMP,1,,0.1
MPDATA,DENS,1,,1300
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,270e6
MPDATA,PRXY,2,,0.20
MPDATA,DAMP,2,,0.08
MPDATA,DENS,2,,1500
!*
!---- Modeling and meshing
block,0,0.25,0,0.5,0,0.25
TYPE,1
MAT,1
LESIZE,1, , ,2, ,1, , ,0
LESIZE,2, , ,1, ,1, , ,0
LESIZE,9, , ,1, ,1, , ,0
vmesh,all
Mat,2
extopt,esize,2,0 ,0
asel,s,area,,4
vext,all,,,,0.5
!-----BCs
allsel
nummrg,node
numcmp,node
numcmp,elem
Allsel
IC,all,all,0,0
nsel,s,loc,Y,0
D,all,UX,0
D,all,UY,0
D,all,UZ,0
Allsel
!---- Load
F,13,FZ,-1000
Allsel
fini !
!----- Solution ----------------------!
!----------------TRANSIENT -------------!
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,0
NSUBST,1,1,1
TIME,0.1
/solu
SOLVE
fini
!------------  MODAL --------------------!
!/solu
!allsel
!outres,all,all
!ANTYPE,MODAL
!MODOPT,QRDAMP,10
!MXPAND,10
!SOLVE
!fini
!----------------Printing files ---------!
/AUX2
file,file,full
hbmat,HBMstiff,txt,,ascii,stiff,yes,yes
hbmat,HBMmass,txt,,ascii,mass,yes,yes
hbmat,HBMdamp,txt,,ascii,damp,yes,yes

FINISH

Thank you

Ahmed El-Ghandour

PhD Candidate at UIC.

--

Ahmed I. El-Ghandour

Research/Teaching Assistant & Ph.D. Candidate Dept. of Civil & Materials Engineering College of Engineering University of Illinois at Chicago (UIC)
Mobile  : +1 312-912-2617
E-mail  : aelgha2@uic.edu
Skype    : a.i.ghando

Xansys-temp mailing list
Xansys-temp@xansystest.info
http://xansystest.info/mailman/listinfo/xansys-temp_xansystest.info
If you are receiving too many emails from XANSYS please consider changing account settings to Digest mode which will send a single email per day.

Please send administrative requests such as deletion from XANSYS to xansys-mod@tynecomp.co.uk and not to the list

Hi Ahmed,

Stiff MATRIX can be generated for all analysis type. Mass matrix can be obtained only for buckling, substructure and modal analysis. Damp matrix only generated for damped modal analysis.

So, in your code, there is not any error in reading .full file for getting stiff, mass and damp matrices in modal analysis. Because in modal analysis, all can be generated from .full file.

In transient analysis, it gives only stiff matrix. While it gives error while generating mass and damp matrices.  Because, as mentioned above, mass and damp matrices cannot be generated in transient analysis.

So, for transient analysis, below codes will work.

[/AUX2
file,file,full
hbmat,HBMstiff,txt,,ascii,stiff,yes,yes]

but, while giving these commands

[hbmat,HBMmass,txt,,ascii,mass,yes,yes
hbmat,HBMdamp,txt,,ascii,damp,yes,yes]

it will generate error.

I hope this information helps you.

Thanks and Regards,
Varia Bhavesh
Electronic Packaging Mechanical Engineer,
Moog Inc.

-----Original Message-----
From: Xansys-temp [mailto:xansys-temp-bounces@xansystest.info] On Behalf Of Grudzinski, James J.
Sent: Wednesday, November 30, 2016 8:54 AM
To: XANSYS Mailing List Temporary Home xansys-temp@xansystest.info
Subject: Re: [Xansys] Extracting Damping matrix

Hi Ahmed,

The code to access the various matrices works as written and has produced all three matrices for the modal analysis.  Your code below needs some comments and a fini removed but otherwise gives a result.  Are you looking in the correct location for the file?  By the way, the result of this command is different for the two kinds of analysis.  See the manual for the various caveats.

Jim

James J. Grudzinski, Ph. D., P.E.
Manager, Engineering Operations
Argonne National Laboratory

-----Original Message-----
From: Xansys-temp [mailto:xansys-temp-bounces@xansystest.info] On Behalf Of Ahmed I. El-Ghandour
Sent: Tuesday, November 29, 2016 5:31 PM
To: xansys-temp@xansystest.info
Subject: [Xansys] Extracting Damping matrix

Hello All,

I am trying to extract the damping matrix for my model, so i start by building a simple cantilever, and I tried two solutions:
1- Modal analysis
2- Transient Analysis.

I found that people usually use hbmat commands for that, but it so not working for me in both analysis.

here is the code I am using in case someone can catch the mistake:

fini
/clear
/prep7
ET,1,SOLID185

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,260e6
MPDATA,PRXY,1,,0.25
MPDATA,DAMP,1,,0.1
MPDATA,DENS,1,,1300
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,270e6
MPDATA,PRXY,2,,0.20
MPDATA,DAMP,2,,0.08
MPDATA,DENS,2,,1500
!*
!---- Modeling and meshing
block,0,0.25,0,0.5,0,0.25
TYPE,1
MAT,1
LESIZE,1, , ,2, ,1, , ,0
LESIZE,2, , ,1, ,1, , ,0
LESIZE,9, , ,1, ,1, , ,0
vmesh,all
Mat,2
extopt,esize,2,0 ,0
asel,s,area,,4
vext,all,,,,0.5
!-----BCs
allsel
nummrg,node
numcmp,node
numcmp,elem
Allsel
IC,all,all,0,0
nsel,s,loc,Y,0
D,all,UX,0
D,all,UY,0
D,all,UZ,0
Allsel
!---- Load
F,13,FZ,-1000
Allsel
fini !
!----- Solution ----------------------!
!----------------TRANSIENT -------------!
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,0
NSUBST,1,1,1
TIME,0.1
/solu
SOLVE
fini
!------------  MODAL --------------------!
!/solu
!allsel
!outres,all,all
!ANTYPE,MODAL
!MODOPT,QRDAMP,10
!MXPAND,10
!SOLVE
!fini
!----------------Printing files ---------!
/AUX2
file,file,full
hbmat,HBMstiff,txt,,ascii,stiff,yes,yes
hbmat,HBMmass,txt,,ascii,mass,yes,yes
hbmat,HBMdamp,txt,,ascii,damp,yes,yes

FINISH

Thank you

Ahmed El-Ghandour

PhD Candidate at UIC.

--

Ahmed I. El-Ghandour

Research/Teaching Assistant & Ph.D. Candidate Dept. of Civil & Materials Engineering College of Engineering University of Illinois at Chicago (UIC)
Mobile  : +1 312-912-2617
E-mail  : aelgha2@uic.edu
Skype    : a.i.ghando

Xansys-temp mailing list
Xansys-temp@xansystest.info
https://urldefense.proofpoint.com/v2/url?u=http-3A__xansystest.info_mailman_listinfo_xansys-2Dtemp-5Fxansystest.info&d=DgICAg&c=ByDrzdYw8tO08sJlHDO_Vg&r=7iPn9zZRhk1nqZxRpPSLF1-6fRrqioMrbyvgo5pp538&m=mq6r3xVY29XVOmxHi5kb7ZU0fmrnTom61Qtc6oIC4uI&s=xT8QfEhMxhg6f-otOUNur-0k1Lb0ua9_Grtc3DPT7_I&e=
If you are receiving too many emails from XANSYS please consider changing account settings to Digest mode which will send a single email per day.

Please send administrative requests such as deletion from XANSYS to xansys-mod@tynecomp.co.uk and not to the list

Xansys-temp mailing list
Xansys-temp@xansystest.info
https://urldefense.proofpoint.com/v2/url?u=http-3A__xansystest.info_mailman_listinfo_xansys-2Dtemp-5Fxansystest.info&d=DgICAg&c=ByDrzdYw8tO08sJlHDO_Vg&r=7iPn9zZRhk1nqZxRpPSLF1-6fRrqioMrbyvgo5pp538&m=mq6r3xVY29XVOmxHi5kb7ZU0fmrnTom61Qtc6oIC4uI&s=xT8QfEhMxhg6f-otOUNur-0k1Lb0ua9_Grtc3DPT7_I&e=
If you are receiving too many emails from XANSYS please consider changing account settings to Digest mode which will send a single email per day.

Please send administrative requests such as deletion from XANSYS to xansys-mod@tynecomp.co.uk and not to the list

Confidentiality Notice: This electronic mail transmission is intended for the use of the individual or entity to which it is addressed and may contain confidential and/or proprietary information belonging to the sender. If you are not the intended recipient, you are hereby notified that any disclosure, use, copying, distribution, or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this transmission in error, please notify the sender immediately by e-mail and delete the original message. Thank you for your cooperation.