2
$\begingroup$

As shown in the following codes,

lamda1 = 1.875; L = 0.2; bL = 0.015; hL = 0.005; rouL = 2550; mL = rouL*bL*hL*L; roL =
mL/L; gama1 = (Cosh[lamda1] + Cos[lamda1])/(Sinh[lamda1] + Sin[lamda1]); 
phi1 = Cosh[lamda1*(x/L)] - Cos[lamda1*(x/L)] - gama1*(Sinh[lamda1*(x/L)] - Sin[lamda1*
(x/L)]); 
fe1 = Cosh[lamda1*(sigmma/L)] - Cos[lamda1*(sigmma/L)] - gama1*(Sinh[lamda1*(sigmma/L)] -
Sin[lamda1*(sigmma/L)]); z1 = D[fe1, sigmma]; z21 = z1^2; 
fec1 = -0.5*Integrate[z21, {sigmma, 0, x}]; ux = fec1*q4[t]^2; uy = phi1*q4[t]; uf = 
{{ux}, {uy}}; AL = {{Cos[q3[t]], -Sin[q3[t]]}, {Sin[q3[t]], Cos[q3[t]]}}; r0 = {{q1[t]},
{q2[t]}}; 
u0 = {{x}, {0}}; RL = r0 + AL . (u0 + uf); REy = Take[RL /. x -> L, -1]; vRL = D[RL, t];
Print["vRL=", MatrixForm[vRL]]; intTL = ExpandAll[Transpose[vRL] . vRL]; Print["intTL=",
intTL]; 
FullSimplify[Re[Integrate[intTL, {x, 0, L}]], Assumptions -> Element[{q1[t], q2[t], 
q3[t], q4[t], Derivative[1][q1][t], Derivative[1][q2][t], Derivative[1][q3][t],
     Derivative[1][q4][t]}, Reals]]

How could I get real numbers instead of imaginary ones?

$\endgroup$
0

1 Answer 1

2
$\begingroup$

You deleted your previous code. Here is a way to deal with the current one:

Chop@FullSimplify[Integrate[intTL, {x, 0, L}], Assumptions -> 
  Element[{q1[t], q2[t], q3[t], q4[t], Derivative[1][q1][t], Derivative[1]q2][t],
                                       Derivative[1][q3][t], Derivative[1][q4][t]}, Reals]]

--- The proposed solution to your previous equation --

Probably

FullSimplify[
 Re@Integrate[-7.79*q4[t]^2*Sin[18.75*x]*Sin[q3[t]]^2*Sinh[18.75*x]* Derivative[1][q4][t]^2 - 
    104*Cos[9.375*x]*Cos[q3[t]]^2*q4[t]^2*Sinh[9.375*x]*Sinh[18.75*x]* Derivative[1][q4][t]^2 - 
    104*Cos[9.375*x]*q4[t]^2*Sin[q3[t]]^2*Sinh[9.375*x]*Sinh[18.75*x]* Derivative[1][q4][t]^2 + 
    13*Cos[q3[t]]^2*q4[t]^2*Sinh[18.75*x]^2*Derivative[1][q4][t]^2 + 
    13*q4[t]^2*Sin[q3[t]]^2*Sinh[18.75*x]^2*Derivative[1][q4][t]^2, {x, 0, 1}], 
 Assumptions -> Element[{q3[t], q4[t], Derivative[1][q4][t]}, Reals]]

(* -> (1.67588*10^15 - 7.90059*10^6 Cos[2 q3[t]]) q4[t]^2 Derivative[1][q4][t]^2 *)
$\endgroup$
4
  • $\begingroup$ I have tried this method to another program, but it failed. and I have sent the code to you by email, could you help me? $\endgroup$
    – SunnySky
    Commented Sep 20, 2012 at 3:11
  • $\begingroup$ @SunnySky Edit your question and add that code. That way, anyone here can help you $\endgroup$ Commented Sep 20, 2012 at 3:13
  • $\begingroup$ Now the post has been modified, could you help me? $\endgroup$
    – SunnySky
    Commented Sep 20, 2012 at 4:13
  • $\begingroup$ @SunnySky You are welcome! $\endgroup$ Commented Sep 20, 2012 at 6:29

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.