Gforge Advanced Server - SDT Support - Answers and discussion http://support.sdtools.com/gf/ Gforge Advanced Server RSS Re: Help to recover my script (function and command)Faraz Sadeghi <faraz.sadeghi86@yahoo.com> Dear Guillaume, Thanks alot for your reply. Yes unfortunately it has been lost, but I have started to write again. Actually I was looking for a way to connect (load) my results (acceleration) to Wavelet apps, and I tried different commands and procedure. Even I do not know how can I solve this issue. To be clear, after running the model I have results such as acc, disp and velo, I am goining to export from sdtool box to apps. I would be greatful if you hint me the way. Regards, Faraz Re: Help to recover my script (function and command)Guillaume Vermot des Roches <vermot@sdtools.com> I am afraid you performed an overwrite operation by directly asking feutilb to replace your script with an m file that generates the model, as it is documented. There is no way of recovering your file like for usual erase/overwrite operations. If you had matlab in auto-save mode and did not do any modification to the overwritten m file, you may check in the same directory for TestA.m~ or TestA.asv that might not have been updated yet. Besides the unfortunate code loss, the question is what did you intend to do ? feutilb('write') is a very specific command that is very different from model saving commands. Kind Regards, Re: Help to recover my script (function and command)Faraz Sadeghi <faraz.sadeghi86@yahoo.com> let me explain again, I wrote a code in Editor window (Matlab SDtool) and several times run the program, and my file is saved as Filename.m (my saved file name is TestA.m) But unfortunately I wrote a function and command name as feutilb(’WriteFileName.m’,model) in command window  and replace my file name (TestA.m) and my model name (modelf) in this function. It changed (actually re-wrote) my code to just number (like matrix) and now I can not see my code (as function and command) in my saved file. My question is how can I restore my script (as function and command)? as I need to see my code (script) as function and command. It was more than 100 line and I do not have backup of my script. Because actually it has been re-wrote and changed. Please find the file in your email. Regards, Faraz Re: Help to recover my script (function and command)Etienne Balmes <balmes@sdtools.com> I am not sur I understand. You are trying to import test results ? or save a FEM model ? feutilb('write') is used to display a clear text version of a FEM model, but in general people will more simple save the model data structure to a .mat file. Help to recover my script (function and command)Faraz Sadeghi <faraz.sadeghi86@yahoo.com> I wrote a code about 100 line in Editor box (script) in Matlab (SDtools), and tried to connect Matlab app to analysis my result (acceleration). I already saved and run several times. Unfortunately when I tried randomly different commands to understand how to import my data to the app, based on the Manual I try this command [[ feutilb('writefilename.m',model) ]] and all my script changed to number like matrix and now I can not go back/undo previous condition. Please let me know how con I recover my script like function and command as I wrote. Thanks Re: Transmissibility for SHMEtienne Balmes <balmes@sdtools.com> Question : What I mean by transmissibility is described in "Damage localization using transmissibility functions: A critical review" Simon Chesné, Arnaud Deraemaeker For example, we may consider a bridge under heavy traffic load and want to do a damage detection on it. As it is not known the accurate quantity of the load, most of the research used TF that is a output response based analysis. Answer : then you should select a loading model (for example a few point loads on the bride) Compute the transfer functions from those loads to the piezo sensors. Evaluate the transmissibility as the ratio of two transfers (you have to remember that for transmissibility computations, the modes found are those where the sensor is fixed, but in the case of a bridge you can probably indeed ignore the change in modal frequencies).   Transmissibility for SHMEtienne Balmes <balmes@sdtools.com> Question : I am going to use the SDtools for my project in Structural Health Monitoring (SHM) and damage detection of civil-structures. Would you please let me know if the SDtoolbox is capable in SHM - damage detection and in particular, Frequency Response Function and Transmissibility Function. If so, could you please give some example (e.g. Damage detection using transmissibility function) or address in the toolbox tutorial? Also, as I could not find any publication related to the aforementioned area in the SDtools web page, please let me know if you have any.   Answer : the documentation for SHM applications is www.sdtools.com/pdf/piezo.pdf you will find there a number of examples of how to compute transfers. When you mention tranmissibility, what you mean exactly is unclear to me. Could you be more specific (from what input to what output).     Re: VoltageMarwa Mallouli <marwa_mel@yahoo.fr> Dear Mr Arnaud Please, can i send you all my work to verify it ? Re: VoltageMarwa Mallouli <marwa_mel@yahoo.fr> Dear Mr Arnaud I want to develop an analytical model of a bimorph cantilever beam with tip mass. Please can you help me to change the equations ?   Re: VoltageArnaud Deraemaeker <aderaema@ulb.ac.be> Dear Marwa, You are using an anaytical formula based on a 1D approximation of the patch, which is not representative of the true behavior of your piezoelectric patch, so that your FRF is probably not representative of the behavior of your system. I advise you to use the piezoelectric models available in SDT in order to do a proper finite element modeling of this system. See the tutorial and examples in http://www.sdtools.com/pdf/piezo.pdf. Arnaud Re: VoltageMarwa Mallouli <marwa_mel@yahoo.fr> Hello The shape of the beam is  figure 1 (the bimorph is composed of  two piezoelectric MFC patches below and above the beam) The input is attached. The output is the voltage formula is mentioned in the figure 1.   Re: VoltageArnaud Deraemaeker <aderaema@ulb.ac.be> Dear Marwa, Can you provide the source code ? It is not possible to answer your question without knowing the shape of your beam as well as what is the input/output of your FRF. Arnaud VoltageMarwa Mallouli <marwa_mel@yahoo.fr>  Hello The response of a bimorph cantilever beam with piezoelectric patches is attached (Voltage). Please can you indicate me if the shape with anti-resonance is true or not ? How to compute geometric stiffness for volumesEtienne Balmes <balmes@sdtools.com> Question : I would like to account for stiffening effects due to pretension in computing modes of a volume structure. Answer : stiffnening effects are due to the existence of geometrically non-linear effects. You must thus Perform a static computation to your loads (calling fe_simul static for example) and store the response in the stack as mdl=stack_set(mdl,'curve','StaticState',def). Assemble matrices that account for prestiffening effects (note the MatDes=5 and not 1). Typically a call of the form  [SE,Case]=fe_case(mdl,'assemble -matdes 2 5 -SE NoT') Compute modes  def=fe_eig({SE.K,Case.T},eigopt) Re: white noise excitationSarehati Umar <sarehati@gmail.com>   I want to simulate a structure with white noise excitations to get the raw time-series output data. How can I do that using fe_time?   Thank you. Global/local matricesEtienne Balmes <balmes@sdtools.com> Question : I have created a model and I use these functions:          mo4=fe_case(mo4,'FixDof','clamped dofs','z==0');           [mo4,Case]=fe_case('assemble -MatDes 2 1 -NoT -SE',mo4);  The output we obtain are here Stiffness and Mass matrices. Are these Global or Local Stiffness/Mass matrices??   Answer : I don’t know what you call local mass matrices. So the answer is global. If you want matrices in local coordinates, then you need to look up the basis trans comman. When importing local coordinate definitions (DID in column 3 of nodes) you will be able to build the local to global coordinate transformation matrix (cGL) which is used as q_global = cGL q_local     Windows / MATLAB > 2014b bug : zoom corrupts screenEtienne Balmes <balmes@sdtools.com> A recent windows update has generated a major bug where calling rbbox for zooming corrupts the screen and shows black squares (or previous screens). This is being investigated by The MathWorks but you can use the following work-around. add the attached rbbox2.m file to the $matlabroot/toolbox/sdt/private directory run rehash toolboxreset to refresh your toolbox path edit iimouse.m too call rbbox2 : in iimouse replace the string rbbox by rbbox2. You can open at the proper line using sdtweb('iimouse','rbbox')   Re: white noise excitationEtienne Balmes <balmes@sdtools.com> Could you be more precise. How do you want the solution to be computed ? time/frequency iterations see script in : sdtweb d_signal shaker time integration in state-space/simulink time integration with SDT fe_time with that info, I will generate an example. white noise excitationSarehati Umar <sarehati_umar09@yahoo.com>  Do you have a demo for a structure which are under white noise excitation? Re: Base excitationEtienne Balmes <balmes@sdtools.com> For harmonic enforced displacement you should use a DofSet entry in fe_simul dfrf. Obviously the node with enforced motion must be part of the model.