yusf
Structural
- May 9, 2006
- 58
hi i have a simple question about usage of penalty method.
according to the constraints we are building P matrix and reconstructing stiffness matrix with following structure
and simply putting it into equation KU=F and then find displacments.
that is ok but how about
during calculations for buckling and dynamic problems
do we use same approach
formulation for buckling
and for dynamic solutions
my quesiton is that do we reconstruct K_g stifening term and M mass matrix such as
for buckling
and for dynamic case
C is constraint matrix....
am i wrong ..if yes for buckling and dynamic case what would be the solution way..
thanks in advance
according to the constraints we are building P matrix and reconstructing stiffness matrix with following structure
Code:
K_new=K+C'*P*C
that is ok but how about
during calculations for buckling and dynamic problems
do we use same approach
formulation for buckling
Code:
det(K+vK_g)= 0 eigenvalue problem
Code:
det(K-Mw^2)=0 eigenvalue problem
where M mass matrix
for buckling
Code:
det(K+vK_g) K=K+C'*P*C and K_g=K_g+C'P*C
Code:
det(K-Mw^2) K=K+C'*P*C and M=M+C'*P*C
am i wrong ..if yes for buckling and dynamic case what would be the solution way..
thanks in advance