Hi everyone,
Iam just wondering how to export the mathematica nb. file to MS word? Preferably word 2003 version.
I have tried to use cut and paste method into word, however, the presentation in word does not look nice.
Any advice is very much appreciated.
Hi all,
I need help to edit the code so that the arrangment of the outputs match the one that Iam looking for.
For instance:
term[n_Integer, m_Integer] :=
Apply[Dot, Table[matg[[i]], {i, m + 1, n}]].hl[[m]] /; (n >= m);
Subscript[lay, s_] := Sum[term[s, i], {i, s}];
In[1]...
Hi there,
I need help about how to solve the flexural rigidity (D) of a composite. For instance, if I have two plies of a plate each with different material properties, i.e. ply 1 has different Young modulus from ply 2. Therefore, how can I determine the value of the overall D of the plate...
Hi there,
I need help about how to solve the flexural rigidity (D) of a composite. For instance, if I have two plies of a plate each with different material properties, i.e. ply 1 has different Young modulus from ply 2. Therefore, how can I determine the value of the overall D of the plate...
Dear all,
I have difficulty to solve this code which involves looping,
In [1] := m = 1; n = 1 ; A = -1 ; B = 5 ;
In [2] := m = 1; n = 3 ; A = 8 ; B = -2 ;
In [3] := m = 1; n = 5 ; A = 9 ; B = 7 ;
In [4] := m = 3; n = 1 ; A = -2 ; B = 8 ;
In [5] := m = 3; n = 3 ; A = 4 ; B = 3 ;
In...
Dear all,
Iam just wondering whether the term 1/a can be expanded using Fourier expansion. If it does, please let me know how to to do this.
Thank for any kind help.
Dale, first of all, I wanted to say, you are a very nice person and very helpful. Iam very pleased with all your comments and advice from the beginning you helped me until now. You have given me the courage to learn to use mathematica in a such a way that makes me feel enthuastic towards...
Thank you again for your help. It works now but I still can't figure out why I still got the warnings. I have tried many many ways to overcome the warnings but still failed including putting back the hl( hL) to h1 (h one).
Iam still wondering though why did you change from h1(h one) to h1 (hL).
Thank you for the comment, I found that the result of the given code shows the product rather than dot product. I mean after running the code, it gives me for instance
In[1]:= Subscript[mata, 4]
Out[1]:= matg[[4]] matg[[3]] matg[[2]] hl[[1]] + matg[[4]] matg[[3]] hl[[2]] + matg[[4]] hl[[3]] +...
Hm, the concept is correct and works, however, there is minor warnings that I need to fix, warnings says:
' Part::partd:Part specification h1[[1]] is longer than depth of object.>>'
Thank you for your kind help
Dear all,
Iam looking for a generic code that could express the following functions:
Subscript[mata, 1] = hl[[1]];
Subscript[mata, 2] = matg[[2]].hl[[1]] + hl[[2]];
Subscript[mata, 3] = matg[[3]].matg[[2]].hl[[1]]
+ matg[[3]].hl[[2]] + hl[[3]];
Subscript[mata, 4] =...
Finally, I have solved my problems, and I got the results. For now, I just need to refine my codes so that it can works faster and efficiently.
Thank you so much for your kind help.
I have checked that those imaginary parts are not correct. When I solve the eigenvalues of a 6 x 6 matrix, I got 6 eigenvalues which contains 2 reals, 4 complex roots. The imaginary parts are not always true, they are just the result of numerical precision loss as you said.
For the matter of discussion, I run the code for iterators of {m,1,11,2},{n,1,11,2}, (meaning there are 36 cycles in total) since code had failed to give results at the stage m = 11 and n = 5, i.e (11,5).
I encountered two warnings after running the code, first being inverse of badly...
The errors starts from this function (simplified version):
dispmat[variables_]:=Module[{local variables}, psimat = psimatrix[variables];.....; disp = (Inverse[mat A]).(mat B - mat C); disp]
where mat A is 3 x 3 matrix, mat B is 3 x 1 (column vector) and mat C is also a column vector (3 x1 )...
I have successfully created the code separately using Module by breaking the codes into smaller parts and the results are as expected (* through hardwork and difficulties*). However, I have another new problem that it can only run up to about 25 cycles or iterations. When I try to run further, I...
Thank you for your comments.
Iam trying to break my codes now using Module eventhough I have no idea on how to do this. How can I separate the codes into modular functions in which they are inter-related to each other. In another words, I still cannot visualise the separations between them...
Frankly speaking, Iam just a human being who is just about to embark into the world of programming. Being a novice user of Mathematica let alone to the world of programming language, this is the best code I could write so far. Iam pretty sure that the code can be further improved and simplifed...
Dear all,
Please help me to get the right codes of sumu17[m,n] and sumv17[m,n] that should be kept inside the Do-Loop.
You can refer to the attached file to view the code.
I appreciate for any comments.
Thank you.
I've got it finally, the code is working now with Module.
Again, thank you so much for your relentless support and guidance. You are a very nice person...you deserve more medals...
all the expressions of u1, u2..., v1, v2,... consists of non linear equation. In fact, some of them have a high very degree of polynomial expressions. Actually, I was trying to get a single expression of u or v in terms of other variables including m, n and other constants from Mathematica, but...
Apologise for my mistakes. You are absolutely right. Your revised code is working perfectly for the above simplified exercise.
However, due to the real application, working your code with my case, its not practical. One of the reasons, my actual expression for u1, u1, v1 and v2 consisting of...
Thanks for the responses. Actually, the code is quite lengthy and I have no idea how to break it into smaller forms. I am in difficulty to break them since they are related to each other. I am also trying to work with Module, and hope it works. I prefer not to display the code here, as the way I...
Dear all,
I am just wondering why I did not get the right outputs from the Do-Loop iteration. When I run the Do-loop, I got the wrong results as compared to the results I obtained for an individual iteration.
For instance:
Case 1 with Do-Loop run from m = 1,3,5 and n = 1,3,5:
In[1]:=...
Sorry...i have deleted my previous message simply because finally I got the right code already using Module. it's finally working.
Thanks for your kind help.
Thank you, yes, its working now using Module for the above case. However, it did not work for another case as shown below:
f[m_, n_] :=
Module[{u1, v1, u2, v2, totalu1, totalv1, totalu2, totalv2,
totalmat}, u1[a_, b_] := a + b;
v1[a_, b_] := 3 a + 5 b;
u2[a_, b_] := 3 a + 4 b...
Thank you for the advise.
I got it. I have solved my problem based on your first advise. However, I wanted to get the code as per your second advise, i.e. to use Module.
Could you please help me how to rectify the code from Do-loop to Module form:
For example:
(*Do Loop*)...
Dear all,
I have a problem with regard to the summation which is done inside the Do-Loop iteration. I have attached the file which includes further information about the defintions, some examples and etc.
Please help me to get the right codes in order to express totalu1, totalv1, totalu2 and...
Dear all,
Iam trying to figure out how to reduce or simplified the flexural rigidity (D) of composite orthotropic laminates consists of many plies into a single layer isotropic model.
I know that for an isotropic case, D = Et^3/12(1-v^2).
However, for the case of, say, 3 plies...
Thank you for your reply.
The main different between code that works (Code 1) and doesn't work (Code 2) is clearly that in the Code 1, I have expressed the expressions of sub[1], sub[2] and sub[3] with eDzstar.topmat, eDzstar.eDz1.topmat and eDzstar.eDz1.eDz1.topmat, respectively. This code1...
Dear all,
I have a problem with my code using the substitution. The code works perfectly without substitution however it does not work with substitution. I mean, I got the results of sub[1], sub[2] and sub[3].
But using the code replacement, I got sub[1] = sub[2] = sub[3], which are not...
Dear all,
After solving the inverse of 3 x 3 matrix, I got a warning which says 'badly inverse conditioned matrix, the result may contain significant numerical errors'.
My question is that how to treat or improve such matrix so that this warning can be eliminated.
For example: the...
Thank you for the feedback. Frankly speaking, being a new programmer, I can't find other way to simplify the code. But I will try my best to find one, no matter what ever it takes.
All the expressions in the code are expressed interms of m and n right from the top to the end, i.e. from matd...
Dear all,
Iam running a code which has only 2 main variables, i.e. m and n using Module. I have tried to run the code as attached but it is not working.
Please help me to improve the code so that I would be able to get qmatrix[m,n].
Perhaps, is there any other better method to use or...
Dear all,
This forum is really great.
Iam just wondering if there is any dicussion with regard to Abaqus matters. If there is one, please let me know which section that I should post the message to.
For anyone who is not familiar about Abaqus, it is a Finite Element Analysis software that...
The only available data given to me are as follows:
-\frac{dU^{o}}{dz} (1-\frac{x}{a}) - \frac{dU^{a}}{dz}(\frac{x}{a})........(1)
U^{o} = U^{o}_{n} *Sin[\frac{n*π*x}{a}].......(2)
U^{a} = U^{a}_{n} *Sin[\frac{n*π*x}{a}].......(3)
Case 1: m\neq0 and n\neq0:
eqn(1) =...
The good news is that if I use Do[...uvmatrix={Simplify[disp0[[1, 1]]], Simplify[disp0[[2, 1]]]} // Chop;Print["m=", m, ",", "n=", n, "-->", uvmatrix], {m, 1, 7, 2}, {n, 1, 7,
2}], I could get the outputs of uvmatrix={u0, v0} at every m and n.
The error started if I put umatrix and vmatrix...
Apologies for the unclear statements.
f1(x) = (1-x/a) can be transformed using Fourier expansion into
bo +\Sigma bm*Cos[m*Pi*x/a] and
f2(x) = x/a can be transformed using Fourier expansion into
ao + \Sigmaam*Cos[m*Pi*x/a]
Iam aware that the original equation of f(x) = ao...
Thank you for the feedback.
Actually, Iam fully aware and understand about the use of Do-loop, in which at every iteration of different m and n, it would return the latest results of usum and vsum. This process would involve in such a way that uvmatrix, umatrix, vmatrix, usum and vsum results...
The code is working perfectly using uvmatrix = Table[...,{m,1,29,2},{n,1,99,2}];umatrix =...;vmatrix =...;usum[m_,n_]:=...vsum[m_,n_]:=...
Thank you.
However, since, the code is partially displayed, I wish to use Do-Loop and put everything inside the Do. Hence, the code shall look...
Dear all,
Please help me to prove the Fourier expansion for three different cases as follows. I need some help to show that L.H.S = R.H.S
Case1: when m not equal to 0 and n not equal to 0.
dU^0/dz (1- x/a) - dU^a/dz (x/a) = -(2/m^2 Pi^2)(1- Cos[m*Pi]) [dUn^0/dz - dUn^a/dz]
Case2...
I just wanted to say thank you for your relentless support and fantastic guidance. I really begin to enjoy Mathematica and this forum is really great.
Best wishes.
Brilliant...u got them perfectly correct, both usum and vsum are now working great. I hope sum[m_,n_]:= and vsum[m_,n_]:= would pick up both the reals and imaginary parts collectively.
You are genius...
Can I just apply all these codes by pasting them in the Do-Loop?